The sensitivity of transcriptomics BMD modeling to the methods used for microarray data normalization

Metabolomics Simultaneously Derives Benchmark Dose Estimates and Discovers Metabolic Biotransformations in a Rat Bioassay

Benchmark dose (BMD) modeling estimates the dose of a chemical that causes a perturbation from baseline. Transcriptional BMDs have been shown to be relatively consistent with apical end point BMDs, opening the door to using molecular BMDs to derive human health-based guidance values for chemical exposure. Metabolomics measures the responses of small-molecule endogenous metabolites to chemical exposure, complementing transcriptomics by characterizing downstream molecular phenotypes that are more closely associated with apical end points. The aim of this study was to apply BMD modeling to in vivo metabolomics data, to compare metabolic BMDs to both transcriptional and apical end point BMDs. This builds upon our previous application of transcriptomics and BMD modeling to a 5-day rat study of triphenyl phosphate (TPhP), applying metabolomics to the same archived tissues. Specifically, liver from rats exposed to five doses of TPhP was investigated using liquid chromatography-mass spectrometry and 1H nuclear magnetic resonance spectroscopy-based metabolomics. Following the application of BMDExpress2 software, 2903 endogenous metabolic features yielded viable dose-response models, confirming a perturbation to the liver metabolome. Metabolic BMD estimates were similarly sensitive to transcriptional BMDs, and more sensitive than both clinical chemistry and apical end point BMDs. Pathway analysis of the multiomics data sets revealed a major effect of TPhP exposure on cholesterol (and downstream) pathways, consistent with clinical chemistry measurements. Additionally, the transcriptomics data indicated that TPhP activated xenobiotic metabolism pathways, which was confirmed by using the underexploited capability of metabolomics to detect xenobiotic-related compounds. Eleven biotransformation products of TPhP were discovered, and their levels were highly correlated with multiple xenobiotic metabolism genes. This work provides a case study showing how metabolomics and transcriptomics can estimate mechanistically anchored points-of-departure. Furthermore, the study demonstrates how metabolomics can also discover biotransformation products, which could be of value within a regulatory setting, for example, as an enhancement of OECD Test Guideline 417 (toxicokinetics).

Characterizing Variability and Uncertainty Associated with Transcriptomic Dose-Response Modeling

Transcriptomics dose-response analysis (TDRA) has emerged as a promising approach for integrating toxicogenomics data into a risk assessment context; however, variability and uncertainty associated with experimental design are not well understood. Here, we evaluated n = 55 RNA-seq profiles derived from Japanese quail liver tissue following exposure to chlorpyrifos (0, 0.04, 0.1, 0.2, 0.4, 1, 2, 4, 10, 20, and 40 μg/g; n = 5 replicates per group) via egg injection. The full dataset was subsampled 637 times to generate smaller datasets with different dose ranges and spacing (designs A-E) and number of replicates (n = 2-5). TDRA of the 637 datasets revealed substantial variability in the gene and pathway benchmark doses, but relative stability in overall transcriptomic point-of-departure (tPOD) values when tPODs were calculated with the "pathway" and "mode" methods. Further, we found that tPOD values were more dependent on the dose range and spacing than on the number of replicates, suggesting that optimal experimental designs should use fewer replicates (n = 2 or 3) and more dose groups to reduce uncertainty in the results. Finally, tPOD values ranged by over ten times for all surveyed experimental designs and tPOD types, suggesting that tPODs should be interpreted as order-of-magnitude estimates.

Adverse Outcome Pathways and Linkages to Transcriptomic Effects Relevant to Ionizing Radiation Injury

BACKGROUND

In the past three decades, a large body of data on the effects of exposure to ionizing radiation and the ensuing changes in gene expression has been generated. These data have allowed for an understanding of molecular-level events and shown a level of consistency in response despite the vast formats and experimental procedures being used across institutions. However, clarity on how this information may inform strategies for health risk assessment needs to be explored. An approach to bridge this gap is the adverse outcome pathway (AOP) framework. AOPs represent an illustrative framework characterizing a stressor associated with a sequential set of causally linked key events (KEs) at different levels of biological organization, beginning with a molecular initiating event (MIE) and culminating in an adverse outcome (AO). Here, we demonstrate the interpretation of transcriptomic datasets in the context of the AOP framework within the field of ionizing radiation by using a lung cancer AOP (AOP 272: https://www.aopwiki.org/aops/272) as a case example.

METHODS

Through the mining of the literature, radiation exposure-related transcriptomic studies in line with AOP 272 related to lung cancer, DNA damage response, and repair were identified. The differentially expressed genes within relevant studies were collated and subjected to the pathway and network analysis using Reactome and GeneMANIA platforms. Identified pathways were filtered (p < 0.001, ≥ 3 genes) and categorized based on relevance to KEs in the AOP. Gene connectivities were identified and further grouped by gene expression-informed associated events (AEs). Relevant quantitative dose-response data were used to inform the directionality in the expression of the genes in the network across AEs.

RESULTS

Reactome analyses identified 7 high-level biological processes with multiple pathways and associated genes that mapped to potential KEs in AOP 272. The gene connectivities were further represented as a network of AEs with associated expression profiles that highlighted patterns of gene expression levels.

CONCLUSIONS

This study demonstrates the application of transcriptomics data in AOP development and provides information on potential data gaps. Although the approach is new and anticipated to evolve, it shows promise for improving the understanding of underlying mechanisms of disease progression with a long-term vision to be predictive of adverse outcomes.

Deep learning identifies erroneous microarray-based, gene-level conclusions in literature

More than 110 000 publications have used microarrays to decipher phenotype-associated genes, clinical biomarkers and gene functions. Microarrays rely on digital assaying the fluorescence signals of arrays. In this study, we retrospectively constructed raw images for 37 724 published microarray data, and developed deep learning algorithms to automatically detect systematic defects. We report that an alarming amount of 26.73% of the microarray-based studies are affected by serious imaging defects. By literature mining, we found that publications associated with these affected microarrays have reported disproportionately more biological discoveries on the genes in the contaminated areas compared to other genes. 28.82% of the gene-level conclusions reported in these publications were based on measurements falling into the contaminated area, indicating severe, systematic problems caused by such contaminations. We provided the identified published, problematic datasets, affected genes and the imputed arrays as well as software tools for scanning such contamination that will become essential to future studies to scrutinize and critically analyze microarray data.

Classifying polycyclic aromatic hydrocarbons by carcinogenic potency using in vitro biosignatures

One of the most difficult challenges for risk assessment is evaluation of chemicals that predominately co-occur in mixtures like polycyclic aromatic hydrocarbons (PAHs). We previously developed a classification model in which systems biology data collected from mice short-term after chemical exposure accurately predict tumor outcome. The present study demonstrates translation of this approach into a human in vitro model in which chemical-specific bioactivity profiles from 3D human bronchial epithelial cells (HBEC) classify PAHs by carcinogenic potency. Gene expression profiles were analyzed from HBEC exposed to carcinogenic and non-carcinogenic PAHs and classification accuracies were identified for individual pathway-based gene sets. Posterior probabilities of best performing gene sets were combined via Bayesian integration resulting in a classifier with four gene sets, including aryl hydrocarbon receptor signaling, regulation of epithelial mesenchymal transition, regulation of angiogenesis, and cell cycle G2-M. In addition, transcriptional benchmark dose modeling of benzo[a]pyrene (BAP) showed that the most sensitive gene sets to BAP regulation were largely dissimilar from those that best classified PAH carcinogenicity challenging current assumptions that BAP carcinogenicity (and subsequent mode of action) is reflective of overall PAH carcinogenicity. These results illustrate utility of using systems toxicology approaches to analyze global gene expression towards carcinogenic hazard assessment.

Comparative toxicity and liver transcriptomics of legacy and emerging brominated flame retardants following 5-day exposure in the rat

The relative toxicity of three legacy and six emerging brominated flame retardants* was studied in the male Harlan Sprague Dawley rat. The hepatocellular and thyroid toxicity of each flame retardant was evaluated following five-day exposure to each of the nine flame retardants (oral gavage in corn oil) at 0.1-1000 μmol/kg body weight per day. Histopathology and transcriptomic analysis were performed on the left liver lobe. Centrilobular hypertrophy of hepatocytes and increases in liver weight were seen following exposure to two legacy (PBDE-47, HBCD) and to one emerging flame retardant (HCDBCO). Total thyroxine (TT4) concentrations were reduced to the greatest extent after PBDE-47 exposure. The PBDE-47, decaBDE, and HBCD liver transcriptomes were characterized by upregulation of liver disease-related and/or metabolic transcripts. Fewer liver disease or metabolic transcript changes were detected for the other flame retardants studied (TBB, TBPH, TBBPA-DBPE, BTBPE, DBDPE, or HCDBCO). PBDE-47 exhibited the most disruption of hepatocellular toxic endpoints, with the Nrf2 antioxidant pathway transcripts upregulated to the greatest extent, although some activation of this pathway also occurred after decaBDE, HBCD, TBB, and HCBCO exposure. These studies provide information that can be used for prioritizing the need for more in-depth brominated flame retardant toxicity studies.