Integration of curated and high-throughput screening data to elucidate environmental influences on disease pathways

Data-Driven Characterization of Genetic Variability in Disease Pathways and Pesticide-Induced Nervous System Disease in the United States Population

BACKGROUND

Genetic susceptibility to chemicals is incompletely characterized. However, nervous system disease development following pesticide exposure can vary in a population, implying some individuals may have higher genetic susceptibility to pesticide-induced nervous system disease.

OBJECTIVES

We aimed to build a computational approach to characterize single-nucleotide polymorphisms (SNPs) implicated in chemically induced adverse outcomes and used this framework to assess the link between differential population susceptibility to pesticides and human nervous system disease.

METHODS

We integrated publicly available datasets of Chemical-Gene, Gene-Pathway, and SNP-Disease associations to build Chemical-Pathway-Gene-SNP-Disease linkages for humans. As a case study, we integrated these linkages with spatialized pesticide application data for the US from 1992 to 2018 and spatialized nervous system disease rates for 2018. Through this, we characterized SNPs that may be important in states with high disease occurrence based on the pesticides used there.

RESULTS

We found that the number of SNP hits per pesticide in US states positively correlated with disease incidence and prevalence for Alzheimer's disease, Parkinson disease, and multiple sclerosis. We performed frequent itemset mining to differentiate pesticides used over time in states with high and low disease occurrence and found that only 19% of pesticide sets overlapped between 10 states with high disease occurrence and 10 states with low disease occurrence rates, and more SNPs were implicated in pathways in high disease occurrence states. Through a cross-validation of subsets of five high and low disease occurrence states, we characterized SNPs, genes, pathways, and pesticides more frequently implicated in high disease occurrence states.

DISCUSSION

Our findings support that pesticides contribute to nervous system disease, and we developed priority lists of SNPs, pesticides, and pathways for further study. This data-driven approach can be adapted to other chemicals, diseases, and locations to characterize differential population susceptibility to chemical exposures. https://doi.org/10.1289/EHP14108.

CTD tetramers: a new online tool that computationally links curated chemicals, genes, phenotypes, and diseases to inform molecular mechanisms for environmental health

The molecular mechanisms connecting environmental exposures to adverse endpoints are often unknown, reflecting knowledge gaps. At the Comparative Toxicogenomics Database (CTD), we developed a bioinformatics approach that integrates manually curated, literature-based interactions from CTD to generate a "CGPD-tetramer": a 4-unit block of information organized as a step-wise molecular mechanism linking an initiating Chemical, an interacting Gene, a Phenotype, and a Disease outcome. Here, we describe a novel, user-friendly tool called CTD Tetramers that generates these evidence-based CGPD-tetramers for any curated chemical, gene, phenotype, or disease of interest. Tetramers offer potential solutions for the unknown underlying mechanisms and intermediary phenotypes connecting a chemical exposure to a disease. Additionally, multiple tetramers can be assembled to construct detailed modes-of-action for chemical-induced disease pathways. As well, tetramers can help inform environmental influences on adverse outcome pathways (AOPs). We demonstrate the tool's utility with relevant use cases for a variety of environmental chemicals (eg, perfluoroalkyl substances, bisphenol A), phenotypes (eg, apoptosis, spermatogenesis, inflammatory response), and diseases (eg, asthma, obesity, male infertility). Finally, we map AOP adverse outcome terms to corresponding CTD terms, allowing users to query for tetramers that can help augment AOP pathways with additional stressors, genes, and phenotypes, as well as formulate potential AOP disease networks (eg, liver cirrhosis and prostate cancer). This novel tool, as part of the complete suite of tools offered at CTD, provides users with computational datasets and their supporting evidence to potentially fill exposure knowledge gaps and develop testable hypotheses about environmental health.

Comparative Toxicogenomics Database (CTD): update 2023

The Comparative Toxicogenomics Database (CTD; http://ctdbase.org/) harmonizes cross-species heterogeneous data for chemical exposures and their biological repercussions by manually curating and interrelating chemical, gene, phenotype, anatomy, disease, taxa, and exposure content from the published literature. This curated information is integrated to generate inferences, providing potential molecular mediators to develop testable hypotheses and fill in knowledge gaps for environmental health. This dual nature, acting as both a knowledgebase and a discoverybase, makes CTD a unique resource for the scientific community. Here, we report a 20% increase in overall CTD content for 17 100 chemicals, 54 300 genes, 6100 phenotypes, 7270 diseases and 202 000 exposure statements. We also present CTD Tetramers, a novel tool that computationally generates four-unit information blocks connecting a chemical, gene, phenotype, and disease to construct potential molecular mechanistic pathways. Finally, we integrate terms for human biological media used in the CTD Exposure module to corresponding CTD Anatomy pages, allowing users to survey the chemical profiles for any tissue-of-interest and see how these environmental biomarkers are related to phenotypes for any anatomical site. These, and other webpage visual enhancements, continue to promote CTD as a practical, user-friendly, and innovative resource for finding information and generating testable hypotheses about environmental health.

Development of the InTelligence And Machine LEarning (TAME) Toolkit for Introductory Data Science, Chemical-Biological Analyses, Predictive Modeling, and Database Mining for Environmental Health Research

Research in environmental health is becoming increasingly reliant upon data science and computational methods that can more efficiently extract information from complex datasets. Data science and computational methods can be leveraged to better identify relationships between exposures to stressors in the environment and human disease outcomes, representing critical information needed to protect and improve global public health. Still, there remains a critical gap surrounding the training of researchers on these in silico methods. We aimed to address this gap by developing the inTelligence And Machine lEarning (TAME) Toolkit, promoting trainee-driven data generation, management, and analysis methods to “TAME” data in environmental health studies. Training modules were developed to provide applications-driven examples of data organization and analysis methods that can be used to address environmental health questions. Target audiences for these modules include students, post-baccalaureate and post-doctorate trainees, and professionals that are interested in expanding their skillset to include recent advances in data analysis methods relevant to environmental health, toxicology, exposure science, epidemiology, and bioinformatics/cheminformatics. Modules were developed by study coauthors using annotated script and were organized into three chapters within a GitHub Bookdown site. The first chapter of modules focuses on introductory data science, which includes the following topics: setting up R/RStudio and coding in the R environment; data organization basics; finding and visualizing data trends; high-dimensional data visualizations; and Findability, Accessibility, Interoperability, and Reusability (FAIR) data management practices. The second chapter of modules incorporates chemical-biological analyses and predictive modeling, spanning the following methods: dose-response modeling; machine learning and predictive modeling; mixtures analyses; -omics analyses; toxicokinetic modeling; and read-across toxicity predictions. The last chapter of modules was organized to provide examples on environmental health database mining and integration, including chemical exposure, health outcome, and environmental justice indicators. Training modules and associated data are publicly available online (https://uncsrp.github.io/Data-Analysis-Training-Modules/). Together, this resource provides unique opportunities to obtain introductory-level training on current data analysis methods applicable to 21st century science and environmental health.

Associating diethylhexyl phthalate to gestational diabetes mellitus via adverse outcome pathways using a network-based approach

Gestational diabetes mellitus (GDM) is a common pregnancy complication that is harmful to both the woman and fetus. Several epidemiological studies have found that exposure to diethylhexyl phthalate (DEHP), an endocrine disruptor ubiquitous in the environment, may be associated with GDM. This study aims to investigate the mechanism between DEHP and GDM using the adverse outcome pathway (AOP) framework, which can integrate information from different sources to elucidate the causal pathways between chemicals and adverse outcomes. We applied a network-based workflow to integrate diverse information to generate computational AOPs and accelerate the AOP development. The interactions among DEHP, genes, phenotypes, and GDM were retrieved from several publicly available databases, including the Comparative Toxicogenomics Database (CTD), Computational Toxicology (CompTox) Chemicals Dashboard, DisGeNET, MalaCards, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG). Based on the above interactions, a DEHP-Gene-Phenotype-GDM network consisting of 52 nodes and 227 edges was formed to support AOP construction. The filtered genes and phenotypes were assembled as molecular initiating events (MIEs) and key events (KEs) according to the upstream and downstream relationships, generating a computational AOP (cAOP) network. Based on the Organization for Economic Co-operation and Development handbook of AOPs, a cAOP was assessed and applied to determine the effects of DEHP on GDM. DEHP could increase TNF-α, downregulate the glucose uptake process, and lead to GDM. Overall, this study revealed the utility of computational methods in integrating a variety of datasets, supporting AOP development, and facilitating a better understanding of the underlying mechanism of exposure to chemicals on human health.

A toxicity pathway-based approach for modeling the mode of action framework of lead-induced neurotoxicity

BACKGROUND

The underlying mechanisms of lead (Pb) toxicity are not fully understood, which makes challenges to the traditional risk assessment. There is growing use of the mode of action (MOA) for risk assessment by integration of experimental data and system biology. The current study aims to develop a new pathway-based MOA for assessing Pb-induced neurotoxicity.

METHODS

The available Comparative Toxicogenomic Database (CTD) was used to search genes associated with Pb-induced neurotoxicity followed by developing toxicity pathways using Ingenuity Pathway Analysis (IPA). The spatiotemporal sequence of disturbing toxicity pathways and key events (KEs) were identified by upstream regulator analysis. The MOA framework was constructed by KEs in biological and chronological order.

RESULTS

There were a total of 71 references showing the relationship between lead exposure and neurotoxicity, which contained 2331 genes. IPA analysis showed that the neuroinflammation signaling pathway was the core toxicity pathway in the enriched pathways relevant to Pb-induced neurotoxicity. The upstream regulator analysis demonstrated that the aryl hydrocarbon receptor (AHR) signaling pathway was the upstream regulator of the neuroinflammation signaling pathway (11.76% overlap with upstream regulators, |Z-score|=1.451). Therefore, AHR activation was recognized as the first key event (KE1) in the MOA framework. The following downstream molecular and cellular key events were also identified. The pathway-based MOA framework of Pb-induced neurotoxicity was built starting with AHR activation, followed by an inflammatory response and neuron apoptosis.

CONCLUSION

Our toxicity pathway-based approach not only advances the development of risk assessment for Pb-induced neurotoxicity but also brings new insights into constructing MOA frameworks of risk assessment for new chemicals.

Distinct genetic regions are associated with differential population susceptibility to chemical exposures

Interactions between environmental factors and genetics underlie the majority of chronic human diseases. Chemical exposures are likely an underestimated contributor, yet gene-environment (GxE) interaction studies rarely assess their modifying effects. Here, we describe a novel method to profile the human genome and identify regions associated with differential population susceptibility to chemical exposures. Single nucleotide polymorphisms (SNPs) implicated in enriched chemical-disease intersections were identified and validated for three chemical classes with expected GxE interaction potential (neuroactive, hepatoactive, and cardioactive compounds). The same approach was then used to characterize consumer product classes with unknown risk for GxE interactions (washing products, cosmetics, and adhesives). Additionally, high-risk variant sets that may confer differential population susceptibility were identified for these consumer product groups through frequent itemset mining and pathway analysis. A dataset of 2454 consumer product chemical-disease linkages, with risk values, SNPs, and pathways for each association was developed, describing the interplay between environmental factors and genetics in human disease progression. We found that genetic hotspots implicated in GxE interactions differ across chemical classes (e.g., washing products had high-risk SNPs implicated in nervous system disease) and illustrate how this approach can discover new associations (e.g., washing product n-butoxyethanol implicated SNPs in the PI3K-Akt signaling pathway for Alzheimer's disease). Hence, our approach can predict high-risk genetic regions for differential population susceptibility to chemical exposures and characterize chemical modifying factors in specific diseases. These methods show promise for describing how chemical exposures can lead to varied health outcomes in a population and for incorporating inter-individual variability into chemical risk assessment.

Investigating Molecular Mechanisms of Immunotoxicity and the Utility of ToxCast for Immunotoxicity Screening of Chemicals Added to Food

The development of high-throughput screening methodologies may decrease the need for laboratory animals for toxicity testing. Here, we investigate the potential of assessing immunotoxicity with high-throughput screening data from the U.S. Environmental Protection Agency ToxCast program. As case studies, we analyzed the most common chemicals added to food as well as per- and polyfluoroalkyl substances (PFAS) shown to migrate to food from packaging materials or processing equipment. The antioxidant preservative tert-butylhydroquinone (TBHQ) showed activity both in ToxCast assays and in classical immunological assays, suggesting that it may affect the immune response in people. From the PFAS group, we identified eight substances that can migrate from food contact materials and have ToxCast data. In epidemiological and toxicological studies, PFAS suppress the immune system and decrease the response to vaccination. However, most PFAS show weak or no activity in immune-related ToxCast assays. This lack of concordance between toxicological and high-throughput data for common PFAS indicates the current limitations of in vitro screening for analyzing immunotoxicity. High-throughput in vitro assays show promise for providing mechanistic data relevant for immune risk assessment. In contrast, the lack of immune-specific activity in the existing high-throughput assays cannot validate the safety of a chemical for the immune system.

CTD Anatomy: analyzing chemical-induced phenotypes and exposures from an anatomical perspective, with implications for environmental health studies

The Comparative Toxicogenomics Database (CTD) is a freely available public resource that curates and interrelates chemical, gene/protein, phenotype, disease, organism, and exposure data. CTD can be used to address toxicological mechanisms for environmental chemicals and facilitate the generation of testable hypotheses about how exposures affect human health. At CTD, manually curated interactions for chemical-induced phenotypes are enhanced with anatomy terms (tissues, fluids, and cell types) to describe the physiological system of the reported event. These same anatomy terms are used to annotate the human media (e.g., urine, hair, nail, blood, etc.) in which an environmental chemical was assayed for exposure. Currently, CTD uses more than 880 unique anatomy terms to contextualize over 255,000 chemical-phenotype interactions and 167,000 exposure statements. These annotations allow chemical-phenotype interactions and exposure data to be explored from a novel, anatomical perspective. Here, we describe CTD's anatomy curation process (including the construction of a controlled, interoperable vocabulary) and new anatomy webpages (that coalesce and organize the curated chemical-phenotype and exposure data sets). We also provide examples that demonstrate how this feature can be used to identify system- and cell-specific chemical-induced toxicities, help inform exposure data, prioritize phenotypes for environmental diseases, survey tissue and pregnancy exposomes, and facilitate data connections with external resources. Anatomy annotations advance understanding of environmental health by providing new ways to explore and survey chemical-induced events and exposure studies in the CTD framework.

Leveraging the Comparative Toxicogenomics Database to Fill in Knowledge Gaps for Environmental Health: A Test Case for Air Pollution-induced Cardiovascular Disease

Environmental health studies relate how exposures (eg, chemicals) affect human health and disease; however, in most cases, the molecular and biological mechanisms connecting an exposure with a disease remain unknown. To help fill in these knowledge gaps, we sought to leverage content from the public Comparative Toxicogenomics Database (CTD) to identify potential intermediary steps. In a proof-of-concept study, we systematically compute the genes, molecular mechanisms, and biological events for the environmental health association linking air pollution toxicants with 2 cardiovascular diseases (myocardial infarction and hypertension) as a test case. Our approach integrates 5 types of curated interactions in CTD to build sets of "CGPD-tetramers," computationally constructed information blocks relating a Chemical- Gene interaction with a Phenotype and Disease. This bioinformatics strategy generates 653 CGPD-tetramers for air pollution-associated myocardial infarction (involving 5 pollutants, 58 genes, and 117 phenotypes) and 701 CGPD-tetramers for air pollution-associated hypertension (involving 3 pollutants, 96 genes, and 142 phenotypes). Collectively, we identify 19 genes and 96 phenotypes shared between these 2 air pollutant-induced outcomes, and suggest important roles for oxidative stress, inflammation, immune responses, cell death, and circulatory system processes. Moreover, CGPD-tetramers can be assembled into extensive chemical-induced disease pathways involving multiple gene products and sequential biological events, and many of these computed intermediary steps are validated in the literature. Our method does not require a priori knowledge of the toxicant, interacting gene, or biological system, and can be used to analyze any environmental chemical-induced disease curated within the public CTD framework. This bioinformatics strategy links and interrelates chemicals, genes, phenotypes, and diseases to fill in knowledge gaps for environmental health studies, as demonstrated for air pollution-associated cardiovascular disease, but can be adapted by researchers for any environmentally influenced disease-of-interest.

New Toxicology Tools and the Emerging Paradigm Shift in Environmental Health Decision-Making

BACKGROUND

Numerous types of rapid toxicity or exposure assays and platforms are providing information relevant to human hazard and exposure identification. They offer the promise of aiding decision-making in a variety of contexts including the regulatory management of chemicals, evaluation of products and environmental media, and emergency response. There is a need to consider both the scientific validity of the new methods and the values applied to a given decision using this new information to ensure that the new methods are employed in ways that enhance public health and environmental protection. In 2018, a National Academies of Sciences, Engineering, and Medicine (NASEM) workshop examined both the toxicological and societal aspects of this challenge.

OBJECTIVES

Our objectives were to explore the challenges of adopting new data streams into regulatory decision-making and highlight the need to align new methods with the information and confidence needs of the decision contexts in which the data may be applied.

METHODS

We go beyond the NASEM workshop to further explore the requirements of different decision contexts. We also call for the new methods to be applied in a manner consistent with the core values of public health and environmental protection. We use the case examples presented in the NASEM workshop to illustrate a range of decision contexts that have applied or could benefit from these new data streams. Organizers of the NASEM workshop came together to further evaluate the main themes from the workshop and develop a joint assessment of the critical needs for improved use of emerging toxicology tools in decision-making. We have drawn from our own experience and individual decision or research contexts as well as from the case studies and panel discussions from the workshop to inform our assessment.

DISCUSSION

Many of the statutes that regulate chemicals in the environment place a high priority on the protection of public health and the environment. Moving away from the sole reliance on traditional approaches and information sources used in hazard, exposure, and risk assessment, toward the more expansive use of rapidly acquired chemical information via in vitro, in silico, and targeted testing strategies will require careful consideration of the information needed and values considerations associated with a particular decision. In this commentary, we explore the ability and feasibility of using emerging data streams, particularly those that allow for the rapid testing of a large number of chemicals across numerous biological targets, to shift the chemical testing paradigm to one in which potentially harmful chemicals are more rapidly identified, prioritized, and addressed. Such a paradigm shift could ultimately save financial and natural resources while ensuring and preserving the protection of public health. https://doi.org/10.1289/EHP4745.

Determination of chemical-disease risk values to prioritize connections between environmental factors, genetic variants, and human diseases

Traditional methods for chemical risk assessment are too time-consuming and resource-intensive to characterize either the diversity of chemicals to which humans are exposed or how that diversity may manifest in population susceptibility differences. The advent of novel toxicological data sources and their integration with bioinformatic databases affords opportunities for modern approaches that consider gene-environment (GxE) interactions in population risk assessment. Here, we present an approach that systematically links multiple data sources to relate chemical risk values to diseases and gene-disease variants. These data sources include high-throughput screening (HTS) results from Tox21/ToxCast, chemical-disease relationships from the Comparative Toxicogenomics Database (CTD), hazard data from resources like the Integrated Risk Information System, exposure data from the ExpoCast initiative, and gene-variant-disease information from the DisGeNET database. We use these integrated data to identify variants implicated in chemical-disease enrichments and develop a new value that estimates the risk of these associations toward differential population responses. Finally, we use this value to prioritize chemical-disease associations by exploring the genomic distribution of variants implicated in high-risk diseases. We offer this modular approach, termed DisQGOS (Disease Quotient Genetic Overview Score), for relating overall chemical-disease risk to potential for population variable responses, as a complement to methods aiming to modernize aspects of risk assessment.