Dataset Information

Differential reconstructed gene interaction networks for deriving toxicity threshold in chemical risk assessment.

ABSTRACT:

Background

Pathway alterations reflected as changes in gene expression regulation and gene interaction can result from cellular exposure to toxicants. Such information is often used to elucidate toxicological modes of action. From a risk assessment perspective, alterations in biological pathways are a rich resource for setting toxicant thresholds, which may be more sensitive and mechanism-informed than traditional toxicity endpoints. Here we developed a novel differential networks (DNs) approach to connect pathway perturbation with toxicity threshold setting.

Methods

Our DNs approach consists of 6 steps: time-series gene expression data collection, identification of altered genes, gene interaction network reconstruction, differential edge inference, mapping of genes with differential edges to pathways, and establishment of causal relationships between chemical concentration and perturbed pathways. A one-sample Gaussian process model and a linear regression model were used to identify genes that exhibited significant profile changes across an entire time course and between treatments, respectively. Interaction networks of differentially expressed (DE) genes were reconstructed for different treatments using a state space model and then compared to infer differential edges/interactions. DE genes possessing differential edges were mapped to biological pathways in databases such as KEGG pathways.

Results

Using the DNs approach, we analyzed a time-series Escherichia coli live cell gene expression dataset consisting of 4 treatments (control, 10, 100, 1000 mg/L naphthenic acids, NAs) and 18 time points. Through comparison of reconstructed networks and construction of differential networks, 80 genes were identified as DE genes with a significant number of differential edges, and 22 KEGG pathways were altered in a concentration-dependent manner. Some of these pathways were perturbed to a degree as high as 70% even at the lowest exposure concentration, implying a high sensitivity of our DNs approach.

Conclusions

Findings from this proof-of-concept study suggest that our approach has a great potential in providing a novel and sensitive tool for threshold setting in chemical risk assessment. In future work, we plan to analyze more time-series datasets with a full spectrum of concentrations and sufficient replications per treatment. The pathway alteration-derived thresholds will also be compared with those derived from apical endpoints such as cell growth rate.

SUBMITTER: Yang Y

PROVIDER: S-EPMC3851258 | biostudies-literature | 2013

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Differential reconstructed gene interaction networks for deriving toxicity threshold in chemical risk assessment.

Yang Yi Y Maxwell Andrew A Zhang Xiaowei X Wang Nan N Perkins Edward J EJ Zhang Chaoyang C Gong Ping P

BMC bioinformatics 20131009

<h4>Background</h4>Pathway alterations reflected as changes in gene expression regulation and gene interaction can result from cellular exposure to toxicants. Such information is often used to elucidate toxicological modes of action. From a risk assessment perspective, alterations in biological pathways are a rich resource for setting toxicant thresholds, which may be more sensitive and mechanism-informed than traditional toxicity endpoints. Here we developed a novel differential networks (DNs) ...[more]

PMID: 24268022

Similar Datasets

Project description:Quantifying spatiotemporally explicit interactions within animal populations facilitates the understanding of social structure and its relationship with ecological processes. Data from animal tracking technologies (Global Positioning Systems ["GPS"]) can circumvent longstanding challenges in the estimation of spatiotemporally explicit interactions, but the discrete nature and coarse temporal resolution of data mean that ephemeral interactions that occur between consecutive GPS locations go undetected. Here, we developed a method to quantify individual and spatial patterns of interaction using continuous-time movement models (CTMMs) fit to GPS tracking data. We first applied CTMMs to infer the full movement trajectories at an arbitrarily fine temporal scale before estimating interactions, thus allowing inference of interactions occurring between observed GPS locations. Our framework then infers indirect interactions-individuals occurring at the same location, but at different times-while allowing the identification of indirect interactions to vary with ecological context based on CTMM outputs. We assessed the performance of our new method using simulations and illustrated its implementation by deriving disease-relevant interaction networks for two behaviorally differentiated species, wild pigs (Sus scrofa) that can host African Swine Fever and mule deer (Odocoileus hemionus) that can host chronic wasting disease. Simulations showed that interactions derived from observed GPS data can be substantially underestimated when temporal resolution of movement data exceeds 30-min intervals. Empirical application suggested that underestimation occurred in both interaction rates and their spatial distributions. CTMM-Interaction method, which can introduce uncertainties, recovered majority of true interactions. Our method leverages advances in movement ecology to quantify fine-scale spatiotemporal interactions between individuals from lower temporal resolution GPS data. It can be leveraged to infer dynamic social networks, transmission potential in disease systems, consumer-resource interactions, information sharing, and beyond. The method also sets the stage for future predictive models linking observed spatiotemporal interaction patterns to environmental drivers.

Project description:The number of legacy chemicals without toxicity reference values combined with the rate of new chemical development are overwhelming the capacity of the traditional risk assessment paradigm. More efficient approaches are needed to quantitatively estimate chemical risks. In this study, rats were dosed orally with multiple doses of six chemicals for 5 days, 2, 4, and 13 weeks. Target organs were analyzed for traditional histological and organ weight changes and transcriptional changes using microarrays. Histological and organ weight changes in this study and the tumor incidences in the original cancer bioassays were analyzed using benchmark dose (BMD) methods to identify noncancer and cancer points-of-departure. The dose-response changes in gene expression were also analyzed using BMD methods and the responses grouped based on signaling pathways. A comparison of transcriptional BMD values for the most sensitive pathway with BMD values for the noncancer and cancer apical endpoints showed a high degree of correlation at all time points. When the analysis included data from an earlier study with 8 additional chemicals, transcriptional BMD values for the most sensitive pathway were significantly correlated with noncancer (r = 0.827, p = 0.0031) and cancer-related (r = 0.940, p = 0.0002) BMD values at 13 weeks. The average ratio of apical-to-transcriptional BMD values was less than two suggesting that for the current chemicals, transcriptional perturbation did not occur at significantly lower doses than apical responses. Based on our results, we propose a practical framework for application of transcriptomic data to chemical risk assessment. Four to five week-old rats were exposed in dose-response and for multiple durations to 6 chemicals. The chemicals were 1,2,4-tribromobenzene (CAS No. 615-54-3), 2,3,4,6-tetrachlorophenol (CAS No. 58-90-2), bromobenzene (CAS No. 108-86-1), 4,4'-methylenebis(N,N-dimethyl)benzenamine (CAS No. 101-61-1), hydrazobenzene (CAS No. 122-66-7), N-nitrosodiphenylamine (CAS No. 86-30-6). For 1,2,4-tribromobenzene and 2,3,4,6-tetrachlorophenol, male Sprague Dawley (Rat/Crl:CD(SD)) rats were used. For bromobenzene and hydroazobenzene, male F344/DuCrl rats were used. For 4,4'-methylenebis(N,N-dimethyl)benzenamine and N-nitrosodiphenylamine, female F344/DuCrl rats were used. The exposure durations were 5 days, 2 weeks, 4 weeks, and 13 weeks. The dose and route of exposure are provided with the individual sample annotations. With each chemical treatment, a matched vehicle control group was run concurrently with the exposure. Gavage exposures were administered 7 days per week and feed exposures were provided 7 days per week. At the appropriate time point, animals were euthanized with a lethal intraperitoneal injection of sodium pentobarbital. For the liver, slices from each of the lobes were mixed together and used for microarray analysis. For the thyroid, the gland was transversely sectioned and the bottom portion was used for microarray analysis. For the bladder, the organ was longitudinally bisected, a slice was removed for histology and the remainder of the organ was minced and used for microarray analysis. Microarrays were run on five rats per dose per time point.

Project description:Toxicology is traditionally divided between human and eco-toxicology. In the shared pursuit of environmental health, this separation does not account for discoveries made in the comparative studies of animal genomes. Here, we provide evidence on the feasibility of understanding the health impact of chemicals on all animals, including ecological keystone species and humans, based on a significant number of conserved genes and their functional associations to health-related outcomes across much of animal diversity. We test four conditions to understand the value of comparative genomics data to inform mechanism-based human and environmental hazard assessment: (1) genes that are most fundamental for health evolved early during animal evolution; (2) the molecular functions of pathways are better conserved among distantly related species than the individual genes that are members of these pathways; (3) the most conserved pathways among animals are those that cause adverse health outcomes when disrupted; (4) gene sets that serve as molecular signatures of biological processes or disease-states are largely enriched by evolutionarily conserved genes across the animal phylogeny. The concept of homology is applied in a comparative analysis of gene families and pathways among invertebrate and vertebrate species compared with humans. Results show that over 70% of gene families associated with disease are shared among the greatest variety of animal species through evolution. Pathway conservation between invertebrates and humans is based on the degree of conservation within vertebrates and the number of interacting genes within the human network. Human gene sets that already serve as biomarkers are enriched by evolutionarily conserved genes across the animal phylogeny. By implementing a comparative method for chemical hazard assessment, human and eco-toxicology converge towards a more holistic and mechanistic understanding of toxicity disrupting biological processes that are important for health and shared among animals (including humans).

Dataset Information

Differential reconstructed gene interaction networks for deriving toxicity threshold in chemical risk assessment.

Background

Methods

Results

Conclusions

Publications

Differential reconstructed gene interaction networks for deriving toxicity threshold in chemical risk assessment.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets