Guided optimization of ToxPi model weights using a Semi-Automated approach.

The Toxicological Prioritization Index (ToxPi) is a visual analysis and decision support tool for dimension reduction and visualization of high throughput, multi-dimensional feature data. ToxPi was originally developed for assessing the relative toxicity of multiple chemicals or stressors by synthesizing complex toxicological data to provide a single comprehensive view of the potential health effects. It continues to be used for profiling chemicals and has since been applied to other types of "sample" entities, including geospatial (e.g. county-level Covid-19 risk and sites of historical PFAS exposure) and other profiling applications. For any set of features (data collected on a set of sample entities), ToxPi integrates the data into a set of weighted slices that provide a visual profile and a score metric for comparison. This scoring system is highly dependent on user-provided feature weights, yet users often lack knowledge of how to define these feature weights. Common methods for predicting feature weights are generally unusable due to inappropriate statistical assumptions and lack of global distributional expectation. However, users often have an inherent understanding of expected results for a small subset of samples. For example, in chemical toxicity, prior knowledge can often place subsets of chemicals into categories of low, moderate or high toxicity (reference chemicals). Ordinal regression can be used to predict weights based on these response levels that are applicable to the entire feature set, analogous to using positive and negative controls to contextualize an empirical distribution. We propose a semi-supervised method utilizing ordinal regression to predict a set of feature weights that produces the best fit for the known response ("reference") data and subsequently fine-tunes the weights via a customized genetic algorithm. We conduct a simulation study to show when this method can improve the results of ordinal regression, allowing for accurate feature weight prediction and sample ranking in scenarios with minimal response data. To ground-truth the guided weight optimization, we test this method on published data to build a ToxPi model for comparison against expert-knowledge-driven weight assignments.

Polygenic subtype identified in ACCORD trial displays a favorable type 2 diabetes phenotype in the UKBiobank population.

INTRODUCTION: We previously identified a genetic subtype (C4) of type 2 diabetes (T2D), benefitting from intensive glycemia treatment in the Action to Control Cardiovascular Risk in Diabetes (ACCORD) trial. Here, we characterized the population of patients that met the C4 criteria in the UKBiobank cohort. RESEARCH DESIGN AND METHODS: Using our polygenic score (PS), we identified C4 individuals in the UKBiobank and tested C4 status with risk of developing T2D, cardiovascular disease (CVD) outcomes, and differences in T2D medications. RESULTS: C4 individuals were less likely to develop T2D, were slightly older at T2D diagnosis, had lower HbA1c values, and were less likely to be prescribed T2D medications (P < .05). Genetic variants in MAS1 and IGF2R, major components of the C4 PS, were associated with fewer overall T2D prescriptions. CONCLUSION: We have confirmed C4 individuals are a lower risk subpopulation of patients with T2D.

Gene-environment interactions within a precision environmental health framework.

Understanding the complex interplay of genetic and environmental factors in disease etiology and the role of gene-environment interactions (GEIs) across human development stages is important. We review the state of GEI research, including challenges in measuring environmental factors and advantages of GEI analysis in understanding disease mechanisms. We discuss the evolution of GEI studies from candidate gene-environment studies to genome-wide interaction studies (GWISs) and the role of multi-omics in mediating GEI effects. We review advancements in GEI analysis methods and the importance of large-scale datasets. We also address the translation of GEI findings into precision environmental health (PEH), showcasing real-world applications in healthcare and disease prevention. Additionally, we highlight societal considerations in GEI research, including environmental justice, the return of results to participants, and data privacy. Overall, we underscore the significance of GEI for disease prediction and prevention and advocate for integrating the exposome into PEH omics studies.

Adaptive selection at G6PD and disparities in diabetes complications.

Diabetes complications occur at higher rates in individuals of African ancestry. Glucose-6-phosphate dehydrogenase deficiency (G6PDdef), common in some African populations, confers malaria resistance, and reduces hemoglobin A1c (HbA1c) levels by shortening erythrocyte lifespan. In a combined-ancestry genome-wide association study of diabetic retinopathy, we identified nine loci including a G6PDdef causal variant, rs1050828 -T (Val98Met), which was also associated with increased risk of other diabetes complications. The effect of rs1050828 -T on retinopathy was fully mediated by glucose levels. In the years preceding diabetes diagnosis and insulin prescription, glucose levels were significantly higher and HbA1c significantly lower in those with versus without G6PDdef. In the Action to Control Cardiovascular Risk in Diabetes (ACCORD) trial, participants with G6PDdef had significantly higher hazards of incident retinopathy and neuropathy. At the same HbA1c levels, G6PDdef participants in both ACCORD and the Million Veteran Program had significantly increased risk of retinopathy. We estimate that 12% and 9% of diabetic retinopathy and neuropathy cases, respectively, in participants of African ancestry are due to this exposure. Across continentally defined ancestral populations, the differences in frequency of rs1050828 -T and other G6PDdef alleles contribute to disparities in diabetes complications. Diabetes management guided by glucose or potentially genotype-adjusted HbA1c levels could lead to more timely diagnoses and appropriate intensification of therapy, decreasing the risk of diabetes complications in patients with G6PDdef alleles.

Questionnaire-based exposome-wide association studies for common diseases in the Personalized Environment and Genes Study.

The exposome collectively refers to all exposures, beginning in utero and continuing throughout life, and comprises not only standard environmental exposures such as point source pollution and ozone levels but also exposures from diet, medication, lifestyle factors, stress, and occupation. The exposome interacts with individual genetic and epigenetic characteristics to affect human health and disease, but large-scale studies that characterize the exposome and its relationships with human disease are limited. To address this gap, we used extensive questionnaire data from the diverse North Carolina-based Personalized Environment and Genes Study (PEGS, n = 9, 429) to evaluate exposure associations in relation to common diseases. We performed an exposome-wide association study (ExWAS) to examine single exposure models and their associations with 11 common complex diseases, namely allergic rhinitis, asthma, bone loss, fibroids, high cholesterol, hypertension, iron-deficient anemia, ovarian cysts, lower GI polyps, migraines, and type 2 diabetes. Across diseases, we found associations with lifestyle factors and socioeconomic status as well as asbestos, various dust types, biohazardous material, and textile-related exposures. We also found disease-specific associations such as fishing with lead weights and migraines. To differentiate between a replicated result and a novel finding, we used an AI-based literature search and database tool that allowed us to examine the current literature. We found both replicated findings, especially for lifestyle factors such as sleep and smoking across diseases, and novel findings, especially for occupational exposures and multiple diseases.

Interactive data sharing for multiple questionnaire-based exposome-wide association studies and exposome correlations in the Personalized Environment and Genes Study.

The correlations among individual exposures in the exposome, which refers to all exposures an individual encounters throughout life, are important for understanding the landscape of how exposures co-occur, and how this impacts health and disease. Exposome-wide association studies (ExWAS), which are analogous to genome-wide association studies (GWAS), are increasingly being used to elucidate links between the exposome and disease. Despite increased interest in the exposome, tools and publications that characterize exposure correlations and their relationships with human disease are limited, and there is a lack of data and results sharing in resources like the GWAS catalog. To address these gaps, we developed the PEGS Explorer web application to explore exposure correlations in data from the diverse North Carolina-based Personalized Environment and Genes Study (PEGS) that were rigorously calculated to account for differing data types and previously published results from ExWAS. Through globe visualizations, PEGS Explorer allows users to explore correlations between exposures found to be associated with complex diseases. The exposome data used for analysis includes not only standard environmental exposures such as point source pollution and ozone levels but also exposures from diet, medication, lifestyle factors, stress, and occupation. The web application addresses the lack of accessible data and results sharing, a major challenge in the field, and enables users to put results in context, generate hypotheses, and, importantly, replicate findings in other cohorts. PEGS Explorer will be updated with additional results as they become available, ensuring it is an up-to-date resource in exposome science.

Plasma protein signatures of adult asthma.

BACKGROUND: Adult asthma is complex and incompletely understood. Plasma proteomics is an evolving technique that can both generate biomarkers and provide insights into disease mechanisms. We aimed to identify plasma proteomic signatures of adult asthma. METHODS: Protein abundance in plasma was measured in individuals from the Agricultural Lung Health Study (ALHS) (761 asthma, 1095 non-case) and the Atherosclerosis Risk in Communities study (470 asthma, 10,669 non-case) using the SOMAScan 5K array. Associations with asthma were estimated using covariate adjusted logistic regression and meta-analyzed using inverse-variance weighting. Additionally, in ALHS, we examined phenotypes based on both asthma and seroatopy (asthma with atopy (n = 207), asthma without atopy (n = 554), atopy without asthma (n = 147), compared to neither (n = 948)). RESULTS: Meta-analysis of 4860 proteins identified 115 significantly (FDR<0.05) associated with asthma. Multiple signaling pathways related to airway inflammation and pulmonary injury were enriched (FDR<0.05) among these proteins. A proteomic score generated using machine learning provided predictive value for asthma (AUC = 0.77, 95% CI = 0.75-0.79 in training set; AUC = 0.72, 95% CI = 0.69-0.75 in validation set). Twenty proteins are targeted by approved or investigational drugs for asthma or other conditions, suggesting potential drug repurposing. The combined asthma-atopy phenotype showed significant associations with 20 proteins, including five not identified in the overall asthma analysis. CONCLUSION: This first large-scale proteomics study identified over 100 plasma proteins associated with current asthma in adults. In addition to validating previous associations, we identified many novel proteins that could inform development of diagnostic biomarkers and therapeutic targets in asthma management.

Comprehensive evaluation of smoking exposures and their interactions on DNA methylation.

BACKGROUND: Smoking impacts DNA methylation, but data are lacking on smoking-related differential methylation by sex or dietary intake, recent smoking cessation (<1 year), persistence of differential methylation from in utero smoking exposure, and effects of environmental tobacco smoke (ETS). METHODS: We meta-analysed data from up to 15,014 adults across 5 cohorts with DNA methylation measured in blood using Illumina's EPIC array for current smoking (2560 exposed), quit < 1 year (500 exposed), in utero (286 exposed), and ETS exposure (676 exposed). We also evaluated the interaction of current smoking with sex or diet (fibre, folate, and vitamin C). FINDINGS: Using false discovery rate (FDR < 0.05), 65,857 CpGs were differentially methylated in relation to current smoking, 4025 with recent quitting, 594 with in utero exposure, and 6 with ETS. Most current smoking CpGs attenuated within a year of quitting. CpGs related to in utero exposure in adults were enriched for those previously observed in newborns. Differential methylation by current smoking at 4-71 CpGs may be modified by sex or dietary intake. Nearly half (35-50%) of differentially methylated CpGs on the 450 K array were associated with blood gene expression. Current smoking and in utero smoking CpGs implicated 3049 and 1067 druggable targets, including chemotherapy drugs. INTERPRETATION: Many smoking-related methylation sites were identified with Illumina's EPIC array. Most signals revert to levels observed in never smokers within a year of cessation. Many in utero smoking CpGs persist into adulthood. Smoking-related druggable targets may provide insights into cancer treatment response and shared mechanisms across smoking-related diseases. FUNDING: Intramural Research Program of the National Institutes of Health, Norwegian Ministry of Health and Care Services and the Ministry of Education and Research, Chief Scientist Office of the Scottish Government Health Directorates and the Scottish Funding Council, Medical Research Council UK and the Wellcome Trust.

Occupational farm work activities influence workers' indoor home microbiome.

BACKGROUND: Farm work entails a heterogeneous mixture of exposures that vary considerably across farms and farmers. Farm work is associated with various health outcomes, both adverse and beneficial. One mechanism by which farming exposures can impact health is through the microbiome, including the indoor home environment microbiome. It is unknown how individual occupational exposures shape the microbial composition in workers' homes. OBJECTIVES: We investigated associations between farm work activities, including specific tasks and pesticide use, and the indoor microbiome in the homes of 468 male farmers. METHODS: Participants were licensed pesticide applicators, mostly farmers, enrolled in the Agricultural Lung Health Study from 2008 to 2011. Vacuumed dust from participants' bedrooms underwent whole-genome shotgun sequencing for indoor microbiome assessment. Using questionnaire data, we evaluated 6 farm work tasks (processing of either hay, silage, animal feed, fertilizer, or soy/grains, and cleaning grain bins) and 19 pesticide ingredients currently used in the past year, plus 7 banned persistent pesticide ingredients ever used. RESULTS: All 6 work tasks were associated with increased microbial diversity levels, with a positive dose-response for the total number of tasks performed (P = 0.001). All tasks were associated with altered microbial compositions (weighted UniFrac P = 0.001) and with higher abundance of specific microbes, including soil-based commensal microbes such as Haloterrigena. Among the 19 pesticides, current use of glyphosate and past use of lindane were associated with increased microbial diversity (P = 0.02-0.04). Ten currently used pesticides and all 7 banned pesticides were associated with altered microbial composition (P = 0.001-0.04). Six pesticides were associated with differential abundance of certain microbes. DISCUSSION: Different farm activities and exposures can uniquely impact the dust microbiome inside homes. Our work suggests that changes to the home microbiome could serve as one pathway for how occupational exposures impact the health of workers and their cohabitating family members, offering possible future intervention targets.

Model-based evaluation of spatiotemporal data reduction methods with unknown ground truth through optimal visualization and interpretability metrics.

Optimizing and benchmarking data reduction methods for dynamic or spatial visualization and interpretation (DSVI) face challenges due to many factors, including data complexity, lack of ground truth, time-dependent metrics, dimensionality bias and different visual mappings of the same data. Current studies often focus on independent static visualization or interpretability metrics that require ground truth. To overcome this limitation, we propose the MIBCOVIS framework, a comprehensive and interpretable benchmarking and computational approach. MIBCOVIS enhances the visualization and interpretability of high-dimensional data without relying on ground truth by integrating five robust metrics, including a novel time-ordered Markov-based structural metric, into a semi-supervised hierarchical Bayesian model. The framework assesses method accuracy and considers interaction effects among metric features. We apply MIBCOVIS using linear and nonlinear dimensionality reduction methods to evaluate optimal DSVI for four distinct dynamic and spatial biological processes captured by three single-cell data modalities: CyTOF, scRNA-seq and CODEX. These data vary in complexity based on feature dimensionality, unknown cell types and dynamic or spatial differences. Unlike traditional single-summary score approaches, MIBCOVIS compares accuracy distributions across methods. Our findings underscore the joint evaluation of visualization and interpretability, rather than relying on separate metrics. We reveal that prioritizing average performance can obscure method feature performance. Additionally, we explore the impact of data complexity on visualization and interpretability. Specifically, we provide optimal parameters and features and recommend methods, like the optimized variational contractive autoencoder, for targeted DSVI for various data complexities. MIBCOVIS shows promise for evaluating dynamic single-cell atlases and spatiotemporal data reduction models.

Occupational Farm Work Activities Influence Workers' Indoor Home Microbiome.

Background: Farm work entails a heterogeneous mixture of exposures that vary considerably across farms and farmers. Farm work is associated with various health outcomes, both adverse and beneficial. One mechanism by which farming exposures can impact health is through the microbiome, including the indoor built environment microbiome. It is unknown how individual occupational exposures shape the microbial composition in workers' homes. Objectives: We investigated associations between farm work activities, including specific tasks and pesticide use, and the indoor microbiome in the homes of 468 male farmers. Methods: Participants were licensed pesticide applicators, mostly farmers, enrolled in the Agricultural Lung Health Study from 2008-2011. Vacuumed dust from participants' bedrooms underwent whole-genome shotgun sequencing for indoor microbiome assessment. Using questionnaire data, we evaluated 6 farm work tasks (processing of either hay, silage, animal feed, fertilizer, or soy/grains, and cleaning grain bins) and 19 pesticide ingredients currently used in the past year, plus 7 persistent banned pesticide ingredients ever used. Results: All 6 work tasks were associated with increased within-sample microbial diversity, with a positive dose-response for the sum of tasks (p=0.001). All tasks were associated with altered overall microbial compositions (weighted UniFrac p=0.001) and with higher abundance of specific microbes, including soil-based microbes such as Haloterrigena. Among the 19 pesticides, only current use of glyphosate and past use of lindane were associated with increased within-sample diversity (p=0.02-0.04). Ten currently used pesticides and all 7 banned pesticides were associated with altered microbial composition (p=0.001-0.04). Six pesticides were associated with differential abundance of certain microbes. Discussion: Specific farm activities and exposures can impact the dust microbiome inside homes. Our work suggests that occupational farm exposures could impact the health of workers and their families through modifying the indoor environment, specifically the microbial composition of house dust, offering possible future intervention targets.

Race/ethnicity-stratified fine-mapping of the MHC locus reveals genetic variants associated with late-onset asthma.

Introduction: Asthma is a chronic disease of the airways that impairs normal breathing. The etiology of asthma is complex and involves multiple factors, including the environment and genetics, especially the distinct genetic architecture associated with ancestry. Compared to early-onset asthma, little is known about genetic predisposition to late-onset asthma. We investigated the race/ethnicity-specific relationship among genetic variants within the major histocompatibility complex (MHC) region and late-onset asthma in a North Carolina-based multiracial cohort of adults. Methods: We stratified all analyses by self-reported race (i.e., White and Black) and adjusted all regression models for age, sex, and ancestry. We conducted association tests within the MHC region and performed fine-mapping analyses conditioned on the race/ethnicity-specific lead variant using whole-genome sequencing (WGS) data. We applied computational methods to infer human leukocyte antigen (HLA) alleles and residues at amino acid positions. We replicated findings in the UK Biobank. Results: The lead signals, rs9265901 on the 5' end of HLA-B, rs55888430 on HLA-DOB, and rs117953947 on HCG17, were significantly associated with late-onset asthma in all, White, and Black participants, respectively (OR = 1.73, 95%CI: 1.31 to 2.14, p = 3.62 × 10-5; OR = 3.05, 95%CI: 1.86 to 4.98, p = 8.85 × 10-6; OR = 19.5, 95%CI: 4.37 to 87.2, p = 9.97 × 10-5, respectively). For the HLA analysis, HLA-B*40:02 and HLA-DRB1*04:05, HLA-B*40:02, HLA-C*04:01, and HLA-DRB1*04:05, and HLA-DRB1*03:01 and HLA-DQB1 were significantly associated with late-onset asthma in all, White, and Black participants. Conclusion: Multiple genetic variants within the MHC region were significantly associated with late-onset asthma, and the associations were significantly different by race/ethnicity group.

Metagenomics reveals novel microbial signatures of farm exposures in house dust.

Indoor home dust microbial communities, important contributors to human health, are shaped by environmental factors, including farm-related exposures. Advanced metagenomic whole genome shotgun sequencing (WGS) improves detection and characterization of microbiota in the indoor built-environment dust microbiome, compared to conventional 16S rRNA amplicon sequencing (16S). We hypothesized that the improved characterization of indoor dust microbial communities by WGS will enhance detection of exposure-outcome associations. The objective of this study was to identify novel associations of environmental exposures with the dust microbiome from the homes of 781 farmers and farm spouses enrolled in the Agricultural Lung Health Study. We examined various farm-related exposures, including living on a farm, crop versus animal production, and type of animal production, as well as non-farm exposures, including home cleanliness and indoor pets. We assessed the association of the exposures on within-sample alpha diversity and between-sample beta diversity, and the differential abundance of specific microbes by exposure. Results were compared to previous findings using 16S. We found most farm exposures were significantly positively associated with both alpha and beta diversity. Many microbes exhibited differential abundance related to farm exposures, mainly in the phyla Actinobacteria, Bacteroidetes, Firmicutes, and Proteobacteria. The identification of novel differential taxa associated with farming at the genera level, including Rhodococcus, Bifidobacterium, Corynebacterium, and Pseudomonas, was a benefit of WGS compared to 16S. Our findings indicate that characterization of dust microbiota, an important component of the indoor environment relevant to human health, is heavily influenced by sequencing techniques. WGS is a powerful tool to survey the microbial community that provides novel insights on the impact of environmental exposures on indoor dust microbiota. These findings can inform the design of future studies in environmental health.

Pharmacogenomic Analyses Implicate B Cell Developmental Status and MKL1 as Determinants of Sensitivity toward Anti-CD20 Monoclonal Antibody Therapy.

Monoclonal antibody (mAb) therapy directed against CD20 is an important tool in the treatment of B cell disorders. However, variable patient response and acquired resistance remain important clinical challenges. To identify genetic factors that may influence sensitivity to treatment, the cytotoxic activity of three CD20 mAbs: rituximab; ofatumumab; and obinutuzumab, were screened in high-throughput assays using 680 ethnically diverse lymphoblastoid cell lines (LCLs) followed by a pharmacogenomic assessment. GWAS analysis identified several novel gene candidates. The most significant SNP, rs58600101, in the gene MKL1 displayed ethnic stratification, with the variant being significantly more prevalent in the African cohort and resulting in reduced transcript levels as measured by qPCR. Functional validation of MKL1 by shRNA-mediated knockdown of MKL1 resulted in a more resistant phenotype. Gene expression analysis identified the developmentally associated TGFB1I1 as the most significant gene associated with sensitivity. qPCR among a panel of sensitive and resistant LCLs revealed immunoglobulin class-switching as well as differences in the expression of B cell activation markers. Flow cytometry showed heterogeneity within some cell lines relative to surface Ig isotype with a shift to more IgG+ cells among the resistant lines. Pretreatment with prednisolone could partly reverse the resistant phenotype. Results suggest that the efficacy of anti-CD20 mAb therapy may be influenced by B cell developmental status as well as polymorphism in the MKL1 gene. A clinical benefit may be achieved by pretreatment with corticosteroids such as prednisolone followed by mAb therapy.

RYK Gene Expression Associated with Drug Response Variation of Temozolomide and Clinical Outcomes in Glioma Patients.

Temozolomide (TMZ) chemotherapy is an important tool in the treatment of glioma brain tumors. However, variable patient response and chemo-resistance remain exceptionally challenging. Our previous genome-wide association study (GWAS) identified a suggestively significant association of SNP rs4470517 in the RYK (receptor-like kinase) gene with TMZ drug response. Functional validation of RYK using lymphocytes and glioma cell lines resulted in gene expression analysis indicating differences in expression status between genotypes of the cell lines and TMZ dose response. We conducted univariate and multivariate Cox regression analyses using publicly available TCGA and GEO datasets to investigate the impact of RYK gene expression status on glioma patient overall (OS) and progression-free survival (PFS). Our results indicated that in IDH mutant gliomas, RYK expression and tumor grade were significant predictors of survival. In IDH wildtype glioblastomas (GBM), MGMT status was the only significant predictor. Despite this result, we revealed a potential benefit of RYK expression in IDH wildtype GBM patients. We found that a combination of RYK expression and MGMT status could serve as an additional biomarker for improved survival. Overall, our findings suggest that RYK expression may serve as an important prognostic or predictor of TMZ response and survival for glioma patients.

MKX-AS1 Gene Expression Associated with Variation in Drug Response to Oxaliplatin and Clinical Outcomes in Colorectal Cancer Patients.

Oxaliplatin (OXAL) is a commonly used chemotherapy for treating colorectal cancer (CRC). A recent genome wide association study (GWAS) showed that a genetic variant (rs11006706) in the lncRNA gene MKX-AS1 and partnered sense gene MKX could impact the response of genetically varied cell lines to OXAL treatment. This study found that the expression levels of MKX-AS1 and MKX in lymphocytes (LCLs) and CRC cell lines differed between the rs11006706 genotypes, indicating that this gene pair could play a role in OXAL response. Further analysis of patient survival data from the Cancer Genome Atlas (TCGA) and other sources showed that patients with high MKX-AS1 expression status had significantly worse overall survival (HR = 3.2; 95%CI = (1.17-9); p = 0.024) compared to cases with low MKX-AS1 expression status. Alternatively, high MKX expression status had significantly better overall survival (HR = 0.22; 95%CI = (0.07-0.7); p = 0.01) compared to cases with low MKX expression status. These results suggest an association between MKX-AS1 and MKX expression status that could be useful as a prognostic marker of response to OXAL and potential patient outcomes in CRC.

Metagenomics reveals novel microbial signatures of farm exposures in house dust.

Indoor home dust microbial communities, important contributors to human health outcomes, are shaped by environmental factors, including farm-related exposures. Detection and characterization of microbiota are influenced by sequencing methodology; however, it is unknown if advanced metagenomic whole genome shotgun sequencing (WGS) can detect novel associations between environmental exposures and the indoor built-environment dust microbiome, compared to conventional 16S rRNA amplicon sequencing (16S). This study aimed to better depict indoor dust microbial communities using WGS to investigate novel associations with environmental risk factors from the homes of 781 farmers and farm spouses enrolled in the Agricultural Lung Health Study. We examined various farm-related exposures, including living on a farm, crop versus animal production, and type of animal production, as well as non-farm exposures, including home cleanliness and indoor pets. We assessed the association of the exposures on within-sample alpha diversity and between-sample beta diversity, and the differential abundance of specific microbes by exposure. Results were compared to previous findings using 16S. We found most farm exposures were significantly positively associated with both alpha and beta diversity. Many microbes exhibited differential abundance related to farm exposures, mainly in the phyla Actinobacteria, Bacteroidetes, Firmicutes , and Proteobacteria . The identification of novel differential taxa associated with farming at the genera level, including Rhodococcus, Bifidobacterium, Corynebacterium , and Pseudomonas , was a benefit of WGS compared to 16S. Our findings indicate that characterization of dust microbiota, an important component of the indoor environment relevant to human health, is heavily influenced by sequencing techniques. WGS is a powerful tool to survey the microbial community that provides novel insights on the impact of environmental exposures on indoor dust microbiota, and should be an important consideration in designing future studies in environmental health.

ToxicR: A computational platform in R for computational toxicology and dose-response analyses.

The need to analyze the complex relationships observed in high-throughput toxicogenomic and other omic platforms has resulted in an explosion of methodological advances in computational toxicology. However, advancements in the literature often outpace the development of software researchers can implement in their pipelines, and existing software is frequently based on pre-specified workflows built from well-vetted assumptions that may not be optimal for novel research questions. Accordingly, there is a need for a stable platform and open-source codebase attached to a programming language that allows users to program new algorithms. To fill this gap, the Biostatistics and Computational Biology Branch of the National Institute of Environmental Health Sciences, in cooperation with the National Toxicology Program (NTP) and US Environmental Protection Agency (EPA), developed ToxicR, an open-source R programming package. The ToxicR platform implements many of the standard analyses used by the NTP and EPA, including dose-response analyses for continuous and dichotomous data that employ Bayesian, maximum likelihood, and model averaging methods, as well as many standard tests the NTP uses in rodent toxicology and carcinogenicity studies, such as the poly-K and Jonckheere trend tests. ToxicR is built on the same codebase as current versions of the EPA's Benchmark Dose software and NTP's BMDExpress software but has increased flexibility because it directly accesses this software. To demonstrate ToxicR, we developed a custom workflow to illustrate its capabilities for analyzing toxicogenomic data. The unique features of ToxicR will allow researchers in other fields to add modules, increasing its functionality in the future.

Implementing machine learning methods with complex survey data: Lessons learned on the impacts of accounting sampling weights in gradient boosting.

Despite the prominent use of complex survey data and the growing popularity of machine learning methods in epidemiologic research, few machine learning software implementations offer options for handling complex samples. A major challenge impeding the broader incorporation of machine learning into epidemiologic research is incomplete guidance for analyzing complex survey data, including the importance of sampling weights for valid prediction in target populations. Using data from 15, 820 participants in the 1988-1994 National Health and Nutrition Examination Survey cohort, we determined whether ignoring weights in gradient boosting models of all-cause mortality affected prediction, as measured by the F1 score and corresponding 95% confidence intervals. In simulations, we additionally assessed the impact of sample size, weight variability, predictor strength, and model dimensionality. In the National Health and Nutrition Examination Survey data, unweighted model performance was inflated compared to the weighted model (F1 score 81.9% [95% confidence interval: 81.2%, 82.7%] vs 77.4% [95% confidence interval: 76.1%, 78.6%]). However, the error was mitigated if the F1 score was subsequently recalculated with observed outcomes from the weighted dataset (F1: 77.0%; 95% confidence interval: 75.7%, 78.4%). In simulations, this finding held in the largest sample size (N = 10,000) under all analytic conditions assessed. For sample sizes <5,000, sampling weights had little impact in simulations that more closely resembled a simple random sample (low weight variability) or in models with strong predictors, but findings were inconsistent under other analytic scenarios. Failing to account for sampling weights in gradient boosting models may limit generalizability for data from complex surveys, dependent on sample size and other analytic properties. In the absence of software for configuring weighted algorithms, post-hoc re-calculations of unweighted model performance using weighted observed outcomes may more accurately reflect model prediction in target populations than ignoring weights entirely.

Questionnaire-Based Polyexposure Assessment Outperforms Polygenic Scores for Classification of Type 2 Diabetes in a Multiancestry Cohort.

OBJECTIVE: Environmental exposures may have greater predictive power for type 2 diabetes than polygenic scores (PGS). Studies examining environmental risk factors, however, have included only individuals with European ancestry, limiting the applicability of results. We conducted an exposome-wide association study in the multiancestry Personalized Environment and Genes Study to assess the effects of environmental factors on type 2 diabetes. RESEARCH DESIGN AND METHODS: Using logistic regression for single-exposure analysis, we identified exposures associated with type 2 diabetes, adjusting for age, BMI, household income, and self-reported sex and race. To compare cumulative genetic and environmental effects, we computed an overall clinical score (OCS) as a weighted sum of BMI and prediabetes, hypertension, and high cholesterol status and a polyexposure score (PXS) as a weighted sum of 13 environmental variables. Using UK Biobank data, we developed a multiancestry PGS and calculated it for participants. RESULTS: We found 76 significant associations with type 2 diabetes, including novel associations of asbestos and coal dust exposure. OCS, PXS, and PGS were significantly associated with type 2 diabetes. PXS had moderate power to determine associations, with larger effect size and greater power and reclassification improvement than PGS. For all scores, the results differed by race. CONCLUSIONS: Our findings in a multiancestry cohort elucidate how type 2 diabetes odds can be attributed to clinical, genetic, and environmental factors and emphasize the need for exposome data in disease-risk association studies. Race-based differences in predictive scores highlight the need for genetic and exposome-wide studies in diverse populations.