Search

The National Health and Nutrition Examination Survey (NHANES) provides data on the health and environmental exposure of the non-institutionalized US population. Such data have considerable potential to understand how the environment and behaviors impact human health. These data are also currently leveraged to answer public health questions such as prevalence of disease. However, these data need to first be processed before new insights can be derived through large-scale analyses. NHANES data are stored across hundreds of files with multiple inconsistencies. Correcting such inconsistencies takes systematic cross examination and considerable efforts but is required for accurately and reproducibly characterizing the associations between the exposome and diseases. Thus, we developed a set of curated and unified datasets and accompanied code by merging 614 separate files and harmonizing unrestricted data across NHANES III (1988-1994) and Continuous (1999-2018), totaling 134,310 participants and 4,740 variables. The variables convey 1) demographic information, 2) dietary consumption, 3) physical examination results, 4) occupation, 5) questionnaire items (e.g., physical activity, general health status, medical conditions), 6) medications, 7) mortality status linked from the National Death Index, 8) survey weights, 9) environmental exposure biomarker measurements, and 10) chemical comments that indicate which measurements are below or above the lower limit of detection. We also provide a data dictionary listing the variables and their descriptions to help researchers browse the data. We also provide R markdown files to show example codes on calculating summary statistics and running regression models to help accelerate high-throughput analysis and secular trends of the exposome.

Summary Background Mortality risk stratification based on dichotomising a physiological indicator with a cutoff point might not adequately capture increased mortality risk and might not account for non-linear associations. We aimed to characterise the linear and non-linear relationships of 27 physiological indicators with all-cause mortality to evaluate whether the current clinical thresholds are suitable in distinguishing patients at high risk for mortality from those at low risk. Methods For this observational cohort study of the US non-institutionalised population, we used data from adults (≥18 years) included in the 1999–2014 National Health and Nutrition Examination Survey (NHANES) linked with National Death Index mortality data collected from Jan 1, 1999, up until Dec 31, 2015. We used Cox proportional hazards regression models adjusted for age, sex, and race or ethnicity to assess associations of physiological indicators with all-cause mortality. We assessed non-linear associations by discretising the physiological indicator into nine quantiles (termed novemtiles) and by using a weighted sum of cubic polynomials (spline). We used ten-fold cross validation to select the most appropriate model using the concordance index, Nagelkerke R2, and Akaike Information Criterion. We identified the level of each physiological indicator that led to a 10% increase in mortality risk to define our cutoffs used to compare with the current clinical thresholds. Findings We included 47 266 adults of 82 091 assessed for eligibility. 25 (93%) of 27 indicators showed non-linear associations with substantial increases compared with linear models in mortality risk (1·5–2·5-times increase). Height and 60 s pulse were the only physiological indicators to show linear associations. For example, participants with an estimated glomerular filtration rate (GFR) of less than 65 mL/min per 1·73 m2 or between 90–116 mL/min per 1·73 m2 are at moderate (hazard ratio 1–2) mortality risk. Those with a GFR greater than 117 mL/min per 1·73 m2 show substantial (hazard ratio ≥2) mortality risk. Both lower and higher values of cholesterol are associated with increased mortality risk. The current clinical thresholds do not align with our mortality-based cutoffs for fat deposition indices, 60 s pulse, triglycerides, cholesterol-related indicators, alkaline phosphatase, glycohaemoglobin, homoeostatic model assessment of insulin resistance, and GFR. For these indicators, the misalignment suggests the need to consider an additional bound when only one is provided. Interpretation Most clinical indicators were shown to have non-linear associations with all-cause mortality. Furthermore, considering these non-linear associations can help derive reliable cutoffs to complement risk stratification and help inform clinical care delivery. Given the poor alignment with our proposed cutoffs, the current clinical thresholds might not adequately capture mortality risk. Funding Ravitz Family Foundation, Forbes Institute for Cancer Discovery, and National Institutes of Health.

Saliva is a widely used biological sample, especially in pediatric research, containing a heterogenous mixture of immune and epithelial cells. Associations of exposure or disease with saliva DNA methylation can be influenced by cell-type proportions. Here, we developed a saliva cell-type DNA methylation reference panel to estimate interindividual cell-type heterogeneity in whole saliva studies. Saliva was collected from 22 children (7–16 years) and sorted into immune and epithelial cells, using size exclusion filtration and magnetic bead sorting. DNA methylation was measured using the Illumina MethylationEPIC BeadChip. We assessed cell-type differences in DNA methylation profiles and tested for enriched biological pathways. Immune and epithelial cells differed at 181,577 (22.8%) DNA methylation sites (t-test p

BackgroundAccording to the World Health Organization, occupational exposures to hazardous chemicals are estimated to cause over 370,000 premature annual deaths. The risks due to multiple workplace chemical exposures, and those occupations most susceptible to the resulting health effects, remain poorly characterized.ObjectivesThe aim of this study is to identify occupations with elevated toxicant biomarker concentrations and increased health risk associated with toxicant exposures in a working US population from diverse categories of occupation. More specifically, we aim to 1) define differences in chemical exposures based on occupation description, 2) identify occupational groups with similar chemical exposure profiles, and 3) identify occupational groups with chemical biomarker levels exceeding acceptable health-based biomarker levels.MethodsFor this observational study of 51,008 participants, we used data from the 1999-2014 National Health and Nutrition Examination Survey. We characterized differences in chemical exposures by occupational group for 129 chemicals by applying a series of generalized linear models with the outcome as biomarker concentrations and the main predictor as the occupational groups, adjusting for age, sex, race/ethnicity, poverty income ratio, study period, and biomarker of tobacco use. We identified groups of occupations with similar chemical exposure profiles via hierarchical clustering. For each occupational group, we calculated percentages of participants with chemical biomarker levels exceeding acceptable health-based guidelines.ResultsBlue collar workers from “Construction”, “Professional, Scientific, Technical Services”, “Real Estate, Rental, Leasing”, “Manufacturing”, and “Wholesale Trade” have higher biomarker levels of toxic chemicals such as several heavy metals, acrylamide, glycideamide, and several volatile organic compounds compared to their white-collar counterparts. For these toxicants, 1-58% of blue-collar workers from these industries have toxicant concentrations exceeding acceptable levels.DiscussionBlue collar workers have toxicant levels higher relative to their white-collar counterparts, often exceeding acceptable levels associated with noncancer effects. Our findings identify multiple occupations to prioritize for targeted interventions and health policies to monitor and reduce high toxicant exposures.

CYP3A5 and CYP3A4 are the predominant enzymes responsible for tacrolimus metabolism; however only a proportion of the population expresses CYP3A5 secondary to genetic variation. CYP3A5 is expressed in both the intestine and the liver and has been shown to impact both the bioavailability and metabolism of orally administered tacrolimus. Increasing the initial tacrolimus dose by 50% to 100% is recommended in patients who are known CYP3A5 expressers; however, whether this dose adjustment is appropriate for i.v. tacrolimus administration is unclear. The objective of this study was to evaluate the impact of CYP3A5 genotype as well as other pharmacogenes on i.v. tacrolimus exposure to determine whether the current genotype-guided dosing recommendations are appropriate for this formulation. In addition, this study aimed to investigate dose conversion requirements among CYP3A5 genotypes when converting from i.v. to p.o. tacrolimus. This study is a retrospective chart review of all patients who underwent allogeneic stem cell transplantation at Michigan Medicine between June 1, 2014, and March 1, 2018, who received i.v. tacrolimus at the time of their transplantation. Secondary use samples were obtained for genotyping CYP3A5, CYP3A4, and ABCB1. Patient demographic information, tacrolimus dosing and trough levels, and concomitant medications received at the time of tacrolimus trough were collected retrospectively from the patients’ medical records. The i.v. dose-controlled concentration (C/D) and the i.v.:p.o. exposure ratio was calculated for all tacrolimus doses and patients, respectively. The impact of CYP3A5, CYP3A4, and ABCB1 genotypes on the i.v. C/D were evaluated with linear mixed modeling. The impact of CYP3A5 genotype on the i.v.:p.o. ratio was evaluated while controlling for age and concomitant use of an azole inhibitor. CYP3A5 and CYP3A4 genotypes were significantly associated with the i.v. C/D, with CYP3A5 expressers and CYP3A4 rapid metabolizers having 20% lower tacrolimus exposure. Neither genotype remained significant in the multivariable model, although age, hematocrit, and concomitant use of strong azole inhibitors were associated with increased i.v. C/D. When controlling for patient age and sex, CYP3A5 expressers had significantly higher i.v.:p.o. ratios than CYP3A5 nonexpressers (3.42 versus 2.78; P = .04). Post hoc analysis showed that the i.v.:p.o. ratio may differ among different CYP3A5 genotypes and azole inhibitor combinations. This study demonstrates that the current genotype-guided tacrolimus dose adjustment recommendations are inappropriate for CYP3A5 expressers receiving i.v. tacrolimus. Although CYP3A5 genotype is likely a minor contributor to i.v. tacrolimus exposure, genotype, in addition to capturing concomitant CYP3A inhibitors, would likely improve i.v.:p.o. dose conversion selection. © 2021 American Society for Transplantation and Cellular Therapy. Published by Elsevier Inc.

SummaryWe are bioinformatics trainees at the University of Michigan who started a local chapter of Girls Who Code to provide a fun and supportive environment for high school women to learn the power of coding. Our goal was to cover basic coding topics and data science concepts through live coding and hands-on practice. However, we could not find a resource that exactly met our needs. Therefore, over the past three years, we have developed a curriculum and instructional format using Jupyter notebooks to effectively teach introductory Python for data science. This method, inspired by The Carpentries organization, uses bite-sized lessons followed by independent practice time to reinforce coding concepts, and culminates in a data science capstone project using real-world data. We believe our open curriculum is a valuable resource to the wider education community and hope that educators will use and improve our lessons, practice problems, and teaching best practices. Anyone can contribute to our educational materials on GitHub.

Background: According to the World Health Organization, occupational exposures to hazardous chemicals are estimated to cause over 370,000 premature annual deaths. The risks due to multiple workplace chemical exposures, and those occupations most susceptible to the resulting health effects, remain poorly characterized. The aim of this study is to identify occupations with elevated toxicant biomarker concentrations and increased health risk associated with exposure to combinations of toxicants in a working US population from diverse categories of occupation. Methods: For this observational study of 51,008 participants, we used data from the 1999-2014 National Health and Nutrition Examination Survey. We characterized differences in chemical exposures by occupational group for 129 chemicals by applying a series of generalized linear models with the outcome as biomarker concentrations and the main predictor as the occupational groups, adjusting for age, sex, race/ethnicity, poverty income ratio, study period, and biomarker of tobacco use. We identified groups of occupations with similar chemical exposure profiles via hierarchical clustering. For each occupational group, we calculated percentages of participants with chemical biomarker levels exceeding acceptable health-based guidelines. Findings: Blue collar workers from “Construction”, “Other Services”, “Professional, Scientific, Technical Services”, “Real Estate, Rental, Leasing”, “Manufacturing”, and “Wholesale Trade” along with several unemployed groups have higher biomarker levels of toxic chemicals such as lead, cadmium, arsenic, mercury, acrylamide, glycideamide, xylenes, 1,4-dichlorobenzene, tetrachloroethane, methyl tert-butyl ether, and benzene compared to their white-collar counterparts. For these toxicants, 1-58% of blue collar workers from these industries have chemical biomarker concentrations exceeding acceptable levels. Interpretation: Blue collar workers and several unemployed groups have toxicant levels higher relative to their white-collar counterparts, often exceeding acceptable levels associated with noncancer effects. Our findings identify multiple occupations to prioritize for targeted interventions and health policies to monitor and reduce high toxicant exposures in susceptible individuals. Funding Information: Ravitz Family Foundation, University of Michigan Forbes Institute for Cancer Discovery, Harvard Data Science Initiative, University of Michigan Center for Occupational Health and Safety Engineering, and National Institutes of Health (R01 ES028802, P30 ES017885, P30 CA046592, UG3 CA267907). Declaration of Interests: Drs. Nguyen, Colacino, Sartor, and Jolliet has nothing to disclose. At the time of the writing of the manuscript, Dr. Patel was a consultant and shareholder of XY.health. Dr. Patel is also paid by J&J.

Phthalates with potential adverse health effects are being replaced by other phthalates or phthalate alternatives. Little is known about temporal trends of phthalate exposure in pregnant women in the United States. We quantified 16 metabolites of eight phthalates and di(isononyl) cyclohexane-1,2-dicarboxylate (DINCH) in 656 urine samples collected from 192 California pregnant women in 2007-2013 during their second and third trimesters of pregnancy who participated in the MARBLES (Markers of Autism Risk in Babies-Learning Early Signs) study. We used multiple regression to estimate least squares geometric means of phthalate biomarker concentrations and annual percent changes over the study period. Biomarker concentrations of diethyl phthalate (DEP) and three phthalates with known toxicity and adverse health effects (i.e., butyl benzyl phthalate [BBzP], dibutyl phthalate [DBP], di(2-ethylhexyl) phthalate [DEHP]) decreased, while those of di-isobutyl phthalate [DiBP], di-isononyl phthalate [DiNP], and di-n-octyl phthalate [DOP] increased in California pregnant women during our study period. To understand broad social forces that may influence temporal trends and geographic variations in phthalate exposure across countries, we compared our phthalate biomarker concentrations with those of other populations. We observed over a factor of 2 differences in exposure across countries for some phthalate biomarkers and between pregnant and nonpregnant women for DEP.

Saliva is a widely used biological sample, especially in pediatric research, containing a heterogenous mixture of immune and epithelial cells. Associations of exposure or disease with saliva DNA methylation can be influenced by cell-type proportions. Here, we developed a saliva cell-type DNA methylation reference panel to estimate interindividual cell-type heterogeneity in whole saliva studies. Saliva was collected from 22 children (7-16 years) and sorted into immune and epithelial cells, using size exclusion filtration and magnetic bead sorting. DNA methylation was measured using the Illumina MethylationEPIC BeadChip. We assessed cell-type differences in DNA methylation profiles and tested for enriched biological pathways. Immune and epithelial cells differed at 164,793 (20.7%) DNA methylation sites (t-test p < 10-8). Immune cell hypomethylated sites mapped to genes enriched for immune pathways (p < 3.2 × 10-5). Epithelial cell hypomethylated sites were enriched for cornification (p = 5.2 × 10-4), a key process for hard palette formation. Saliva immune and epithelial cells have distinct DNA methylation profiles which can drive whole saliva DNA methylation measures. A primary saliva DNA methylation reference panel, easily implemented with an R package, will allow estimates of cell proportions from whole saliva samples and improve epigenetic epidemiology studies by accounting for measurement heterogeneity by cell-type proportions.

Background: Chemical biomarker concentrations are driven by complex interactions between chemical use patterns, exposure pathways, and toxicokinetic parameters such as biological half-lives. Criteria to differentiate legacy from current exposures are helpful for interpreting variation in age-based and time trends of chemical exposure and identifying chemicals to which children are highly exposed. A systematic approach is needed to study temporal trends for a wide range of chemicals in the US population. Objectives: Using National Health and Nutrition Examination Survey (NHANES) data on measured biomarker concentrations for 141 chemicals from 1999 to 2014, we aim to 1) understand the influence of temporal determinants, in particular time trends, biological half-lives, and restriction dates on age-based trends, 2) systematically define an age-based pattern to identify chemicals with ongoing and high exposure in children, and 3) characterize how age-based trends for six Per- and Polyfluoroalkyl Substances (PFASs) are changing over time. Methods: We performed an integrated analysis of biological half-lives and restriction dates, compared distributions of chemical biomarker concentrations by age group, and then applied a series of regression models to evaluate the linear (βage) and nonlinear (βage2) relationships between age and chemical biomarker levels. Results: For restricted chemicals, a minimum persistence of 1 year in the human body is needed to observe substantial differences between the less exposed young population and historically exposed adults. We define a metric (βage2/βage>126.9) that identifies several phthalates, brominated flame retardants, pesticides, and metals such as lead and tungsten as elevated and ongoing exposures in children. While a substantial reduction in children's exposures was reflected in PFOS and PFOA, levels of PFNA and PFHxS in children were higher in 2013–2014 compared to those in 1999–2000. Conclusions: Integrating a series of regression models with systemized stratified analyses by age group enabled us to define an age-based pattern to identify chemicals that are of higher levels in children. Keywords: Age-based exposure, Temporal trends, Biomonitoring, Environmental chemicals

Background: Stark racial disparities in disease incidence among American women remain a persistent public health challenge. These disparities likely result from complex interactions between genetic, social, lifestyle, and environmental risk factors. The influence of environmental risk factors, such as chemical exposure, however, may be substantial and is poorly understood. Objectives: We quantitatively evaluated chemical-exposure disparities by race/ethnicity, life stage, and time in United States (US) women (n = 38,080) by using biomarker data for 143 chemicals from the National Health and Nutrition Examination Survey (NHANES) 1999–2014. Methods: We applied a series of survey-weighted, generalized linear models using data from the entire NHANES women population along with cycle and age-group stratified subpopulations. The outcome was chemical biomarker concentration, and the main predictor was race/ethnicity with adjustment for age, socioeconomic status, smoking habits, and NHANES cycle. Results: Compared to non-Hispanic White women, the highest disparities were observed for non-Hispanic Black, Mexican American, Other Hispanic, and Other Race/Multi-Racial women with higher levels of pesticides and their metabolites, including 2,5-dichlorophenol, o,p’-DDE, beta-hexachlorocyclohexane, and 2,4-dichlorophenol, along with personal care and consumer product compounds, including parabens and monoethyl phthalate, as well as several metals, such as mercury and arsenic. Moreover, for Mexican American, Other Hispanic, and non-Hispanic black women, there were several exposure disparities that persisted across age groups, such as higher 2,4- and 2,5-dichlorophenol concentrations. Exposure levels for methyl and propyl parabens, however, were the highest in non-Hispanic black compared to non-Hispanic white children with average differences exceeding 4-fold. Exposure disparities for methyl and propyl parabens are increasing over time in Other Race/Multi-Racial women while fluctuating for non-Hispanic Black, Mexican American, and Other Hispanic. Cotinine levels are among the highest in Non-Hispanic White women compared to Mexican American and Other Hispanic women with disparities plateauing and increasing, respectively. Discussion: We systematically evaluated differences in chemical exposures across women of various race/ethnic groups and across age groups and time. Our findings could help inform chemical prioritization in designing epidemiological and toxicological studies. In addition, they could help guide public health interventions to reduce environmental and health disparities across populations. Keywords: Racial disparities, Age-based trends, Temporal trends, Environmental justice, Biomonitoring, Environmental chemicals

Among women, breast cancer is the most prevalent form of cancer worldwide and has the second highest mortality rate of any cancer in the United States. The breast cancer related death rate is 40 % higher in non-Hispanic Black women compared to non-Hispanic White women. The incidence of triple negative breast cancer (TNBC), an aggressive subtype of breast cancer for which there is no targeted therapy, is also approximately three times higher for Black, relative to, White women. The drivers of these differences are poorly understood. Here, we aimed to identify chemical exposures which play a role in breast cancer disparities. Using chemical biomonitoring data from the National Health and Nutrition Examination Survey (NHANES) and biological activity data from the EPA's ToxCast program, we assessed the toxicological profiles of chemicals to which US Black women are disproportionately exposed. We conducted a literature search to identify breast cancer targets in ToxCast to analyze the response of chemicals with exposure disparities in these assays. Forty-three chemical biomarkers are significantly higher in Black women. Investigation of these chemicals in ToxCast resulted in 32,683 assays for analysis, 5172 of which contained nonzero values for the concentration at which the dose-response fitted model reaches the cutoff considered "active". Of these chemicals BPA, PFOS, and thiram are most comprehensively assayed. 2,5-dichlorophenol, 1,4-dichlorobenzene, and methyl and propyl parabens had higher biomarker concentrations in Black women and moderate testing and activity in ToxCast. The distribution of active concentrations for these chemicals in ToxCast assays are comparable to biomarker concentrations in Black women NHANES participants. Through this integrated analysis, we identify that multiple chemicals, including thiram, propylparaben, and p,p' DDE, have disproportionate exposures in Black women and have breast cancer associated biological activity at human exposure relevant doses.

Background: Billions of dollars are spent on environmental dredging (ED) to remediate contaminated sediments, with one goal being reduced human health risks. However, ED may increase health risks i...

Development of sustainable products requires comparing direct impacts on the user and indirect impacts such as environmental pollution in a consistent exposome-based approach. We currently lack too...