Episodic ozone exposure in Long-Evans rats has limited effects on cauda sperm motility and non-coding RNA populations.

Epidemiological evidence suggests the potential for air pollutants to induce male reproductive toxicity. In experimental studies, exposure to ozone during sensitive windows in the sperm lifecycle has been associated with impaired sperm motility. Subsequently, we sought to investigate the effects of episodic exposure to ozone during sperm maturation in the rat. Long-Evans rats were exposed to either filtered air or ozone (0.4 or 0.8 ppm) for five non-consecutive days over two weeks. Ozone exposure did not impact male reproductive organ weights or sperm motility ∼24 hours following the final exposure. Furthermore, circulating sex hormones remained unchanged despite increased T3 and T4 in the 0.8 ppm group. While there was indication of altered adrenergic signaling attributable to ozone exposure in the testis, there were minimal impacts on small non-coding RNAs detected in cauda sperm. Only two piwi-interacting RNAs (piRNAs) were altered in the mature sperm of ozone-exposed rats (piR-rno-346434 and piR-rno-227431). Data across all rats were next analyzed to identify any non-coding RNAs that may be correlated with reduced sperm motility. A total of 7 microRNAs (miRNAs), 8 RNA fragments, and 1682 piRNAs correlated well with sperm motility. Utilizing our exposure paradigm herein, we were unable to substantiate the relationship between ozone exposure during maturation with sperm motility. However, these approaches served to identify a suite of non-coding RNAs that were associated with sperm motility in rats. With additional investigation, these RNAs may prove to have functional roles in the acquisition of motility or be unique biomarkers for male reproductive toxicity.

Persistent gene expression and DNA methylation alterations linked to carcinogenic effects of dichloroacetic acid.

Background: Mechanistic understanding of transient exposures that lead to adverse health outcomes will enhance our ability to recognize biological signatures of disease. Here, we measured the transcriptomic and epigenomic alterations due to exposure to the metabolic reprogramming agent, dichloroacetic acid (DCA). Previously, we showed that exposure to DCA increased liver tumor incidence in B6C3F1 mice after continuous or early life exposures significantly over background level. Methods: Using archived formalin-fixed liver samples, we utilized modern methodologies to measure gene expression and DNA methylation levels to link to previously generated phenotypic measures. Gene expression was measured by targeted RNA sequencing (TempO-seq 1500+ toxicity panel: 2754 total genes) in liver samples collected from 10-, 32-, 57-, and 78-week old mice exposed to deionized water (controls), 3.5 g/L DCA continuously in drinking water ("Direct" group), or DCA for 10-, 32-, or 57-weeks followed by deionized water until sample collection ("Stop" groups). Genome-scaled alterations in DNA methylation were measured by Reduced Representation Bisulfite Sequencing (RRBS) in 78-week liver samples for control, Direct, 10-week Stop DCA exposed mice. Results: Transcriptomic changes were most robust with concurrent or adjacent timepoints after exposure was withdrawn. We observed a similar pattern with DNA methylation alterations where we noted attenuated differentially methylated regions (DMRs) in the 10-week Stop DCA exposure groups compared to the Direct group at 78-weeks. Gene pathway analysis indicated cellular effects linked to increased oxidative metabolism, a primary mechanism of action for DCA, closer to exposure windows especially early in life. Conversely, many gene signatures and pathways reversed patterns later in life and reflected more pro-tumorigenic patterns for both current and prior DCA exposures. DNA methylation patterns correlated to early gene pathway perturbations, such as cellular signaling, regulation and metabolism, suggesting persistence in the epigenome and possible regulatory effects. Conclusion: Liver metabolic reprogramming effects of DCA interacted with normal age mechanisms, increasing tumor burden with both continuous and prior DCA exposure in the male B6C3F1 rodent model.

A 50-gene biomarker identifies estrogen receptor-modulating chemicals in a microarray compendium.

High throughput transcriptomics (HTTr) profiling has the potential to rapidly and comprehensively identify molecular targets of environmental chemicals that can be linked to adverse outcomes. We describe here the construction and characterization of a 50-gene expression biomarker designed to identify estrogen receptor (ER) active chemicals in HTTr datasets. Using microarray comparisons, the genes in the biomarker were identified as those that exhibited consistent directional changes when ER was activated (4 ER agonists; 4 ESR1 gene constitutively active mutants) and opposite directional changes when ER was suppressed (4 antagonist treatments; 4 ESR1 knockdown experiments). The biomarker was evaluated as a predictive tool using the Running Fisher algorithm by comparison to annotated gene expression microarray datasets including those evaluating the transcriptional effects of hormones and chemicals in MCF-7 cells. Depending on the reference dataset used, the biomarker had a predictive accuracy for activation of up to 96%. To demonstrate applicability for HTTr data analysis, the biomarker was used to identify ER activators in a set of 15 chemicals that are considered potential bisphenol A (BPA) alternatives examined at up to 10 concentrations in MCF-7 cells and analyzed by full-genome TempO-Seq. Using benchmark dose (BMD) modeling, the biomarker genes stratified the ER potency of BPA alternatives consistent with previous studies. These results demonstrate that the ER biomarker can be used to accurately identify ER activators in transcript profile data derived from MCF-7 cells.

Early microRNA responses in rodent liver mediated by furan exposure establish dose thresholds for later adverse outcomes.

To reduce reliance on long-term in vivo studies, short-term data linking early molecular-based measurements to later adverse health effects is needed. Although transcriptional-based benchmark dose (BMDT) modeling has been used to estimate potencies and stratify chemicals based on potential to induce later-life effects, dose-responsive epigenetic alterations have not been routinely considered. Here, we evaluated the utility of microRNA (miRNA) profiling in mouse liver and blood, as well as in mouse primary hepatocytes in vitro, to indicate mechanisms of liver perturbation due to short-term exposure of the known rodent liver hepatotoxicant and carcinogen, furan. Benchmark dose modeling of miRNA measurements (BMDmiR) were compared to the referent transcriptional (BMDT) and apical (BMDA) estimates. These analyses indicate a robust dose response for 34 miRNAs to furan and involvement of p53-linked pathways in furan-mediated hepatotoxicity, supporting mRNA and apical measurements. Liver-sourced miRNAs were also altered in the blood and primary hepatocytes. Overall, these results indicate mechanistic involvement of miRNA in furan carcinogenicity and provide evidence of their potential utility as accessible biomarkers of exposure and disease.

Vinyl chloride enhances high-fat diet-induced proteome alterations in the mouse pancreas related to metabolic dysfunction.

Alterations in physiological processes in pancreas have been associated with various metabolic dysfunctions and can result from environmental exposures, such as chemicals and diet. It was reported that environmental vinyl chloride (VC) exposure, a common industrial organochlorine and environmental pollutant, significantly exacerbated metabolic-related phenotypes in mice fed concurrently with high-fat diet (HFD) but not low-fat diet (LFD). However, little is known about the role of the pancreas in this interplay, especially at a proteomic level. The present study was undertaken to examine the protein responses to VC exposure in pancreas tissues of C57BL/6J mice fed LFD or HFD, with focus on the investigation of protein expression and/or phosphorylation levels of key protein biomarkers of carbohydrate, lipid, and energy metabolism, oxidative stress and detoxification, insulin secretion and regulation, cell growth, development, and communication, immunological responses and inflammation, and biomarkers of pancreatic diseases and cancers. We found that the protein alterations may indicate diet-mediated susceptibility in mouse pancreas induced by HFD to concurrent exposure of low levels of inhaled VC. These proteome biomarkers may lead to a better understanding of pancreas-mediated adaptive or adverse response and susceptibility to metabolic disease.

Circulating MicroRNAs, Polychlorinated Biphenyls, and Environmental Liver Disease in the Anniston Community Health Survey.

BACKGROUND: Polychlorinated biphenyl (PCB) exposures have been associated with liver injury in human cohorts, and steatohepatitis with liver necrosis in model systems. MicroRNAs (miRs) maintain cellular homeostasis and may regulate the response to environmental stress. OBJECTIVES: We tested the hypothesis that specific miRs are associated with liver disease and PCB exposures in a residential cohort. METHODS: Sixty-eight targeted hepatotoxicity miRs were measured in archived serum from 734 PCB-exposed participants in the cross-sectional Anniston Community Health Survey. Necrotic and other liver disease categories were defined by serum keratin 18 (K18) biomarkers. Associations were determined between exposure biomarkers (35 ortho-substituted PCB congeners) and disease biomarkers (highly expressed miRs or previously measured cytokines), and Ingenuity Pathway Analysis was performed. RESULTS: The necrotic liver disease category was associated with four up-regulated miRs (miR-99a-5p, miR-122-5p, miR-192-5p, and miR-320a) and five down-regulated miRs (let-7d-5p, miR-17-5p, miR-24-3p, miR-197-3p, and miR-221-3p). Twenty-two miRs were associated with the other liver disease category or with K18 measurements. Eleven miRs were associated with 24 PCBs, most commonly congeners with anti-estrogenic activities. Most of the exposure-associated miRs were associated with at least one serum hepatocyte death, pro-inflammatory cytokine or insulin resistance bioarker, or with both. Within each biomarker category, associations were strongest for the liver-specific miR-122-5p. Pathways of liver toxicity that were identified included inflammation/hepatitis, hyperplasia/hyperproliferation, cirrhosis, and hepatocellular carcinoma. Tumor protein p53 and tumor necrosis factor α were well integrated within the top identified networks. DISCUSSION: These results support the human hepatotoxicity of environmental PCB exposures while elucidating potential modes of PCB action. The MiR-derived liquid liver biopsy represents a promising new technique for environmental hepatology cohort studies. https://doi.org/10.1289/EHP9467.

Genomic comparisons between hepatocarcinogenic and non-hepatocarcinogenic organophosphate insecticides in the mouse liver.

Short-term biomarkers of toxicity have an increasingly important role in the screening and prioritization of new chemicals. In this study, we examined early indicators of liver toxicity for three reference organophosphate (OP) chemicals, which are among the most widely used insecticides in the world. The OP methidathion was previously shown to increase the incidence of liver toxicity, including hepatocellular tumors, in male mice. To provide insights into the adverse outcome pathway (AOP) that underlies these tumors, effects of methidathion in the male mouse liver were examined after 7 and 28 day exposures and compared to those of two other OPs that either do not increase (fenthion) or possibly suppress liver cancer (parathion) in mice. None of the chemicals caused increases in liver weight/body weight or histopathological changes in the liver. Parathion decreased liver cell proliferation after 7 and 28 days while the other chemicals had no effects. There was no evidence for hepatotoxicity in any of the treatment groups. Full-genome microarray analysis of the livers from the 7 and 28 day treatments demonstrated that methidathion and fenthion regulated a large number of overlapping genes, while parathion regulated a unique set of genes. Examination of cytochrome P450 enzyme activities and use of predictive gene expression biomarkers found no consistent evidence for activation of AhR, CAR, PXR, or PPARα. Parathion suppressed the male-specific gene expression pattern through STAT5b, similar to genetic and dietary conditions that decrease liver tumor incidence in mice. Overall, these findings indicate that methidathion causes liver cancer by a mechanism that does not involve common mechanisms of liver cancer induction.

Integrated Omic Analyses Identify Pathways and Transcriptomic Regulators Associated With Chemical Alterations of In Vitro Neural Network Formation.

Development of in vitro new approach methodologies has been driven by the need for developmental neurotoxicity (DNT) hazard data on thousands of chemicals. The network formation assay characterizes DNT hazard based on changes in network formation but provides no mechanistic information. This study investigated nervous system signaling pathways and upstream physiological regulators underlying chemically induced neural network dysfunction. Rat primary cortical neural networks grown on microelectrode arrays were exposed for 12 days in vitro to cytosine arabinoside, 5-fluorouracil, domoic acid, cypermethrin, deltamethrin, or haloperidol as these exposures altered network formation in previous studies. RNA-seq from cells and gas chromatography/mass spectrometry analysis of media extracts collected on days in vitro 12 provided gene expression and metabolomic identification, respectively. The integration of differentially expressed genes and metabolites for each neurotoxicant was analyzed using ingenuity pathway analysis. All 6 compounds altered gene expression that linked to developmental disorders and neurological diseases. Other enriched canonical pathways overlapped among compounds of the same class; eg, genes and metabolites altered by both cytosine arabinoside and 5-fluorouracil exposures are enriched in axonal guidance pathways. Integrated analysis of upstream regulators was heterogeneous across compounds, but identified several transcriptomic regulators including CREB1, SOX2, NOTCH1, and PRODH. These results demonstrate that changes in network formation are accompanied by transcriptomic and metabolomic changes and that different classes of compounds produce differing responses. This approach can enhance information obtained from new approach methodologies and contribute to the identification and development of adverse outcome pathways associated with DNT.

Progress towards an OECD reporting framework for transcriptomics and metabolomics in regulatory toxicology.

Omics methodologies are widely used in toxicological research to understand modes and mechanisms of toxicity. Increasingly, these methodologies are being applied to questions of regulatory interest such as molecular point-of-departure derivation and chemical grouping/read-across. Despite its value, widespread regulatory acceptance of omics data has not yet occurred. Barriers to the routine application of omics data in regulatory decision making have been: 1) lack of transparency for data processing methods used to convert raw data into an interpretable list of observations; and 2) lack of standardization in reporting to ensure that omics data, associated metadata and the methodologies used to generate results are available for review by stakeholders, including regulators. Thus, in 2017, the Organisation for Economic Co-operation and Development (OECD) Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) launched a project to develop guidance for the reporting of omics data aimed at fostering further regulatory use. Here, we report on the ongoing development of the first formal reporting framework describing the processing and analysis of both transcriptomic and metabolomic data for regulatory toxicology. We introduce the modular structure, content, harmonization and strategy for trialling this reporting framework prior to its publication by the OECD.

Methodological considerations for measuring biofluid-based microRNA biomarkers.

MicroRNAs (miRNAs) are small non-coding RNA that regulate the expression of messenger RNA and are implicated in almost all cellular processes. Importantly, miRNAs can be released extracellularly and are stable in these matrices where they may serve as indicators of organ or cell-specific toxicity, disease, and biological status. There has thus been great enthusiasm for developing miRNAs as biomarkers of adverse outcomes for scientific, regulatory, and clinical purposes. Despite advances in measurement capabilities for miRNAs, miRNAs are still not routinely employed as noninvasive biomarkers. This is in part due to the lack of standard approaches for sample preparation and miRNA measurement and uncertainty in their biological interpretation. Members of the microRNA Biomarkers Workgroup within the Health and Environmental Sciences Institute's (HESI) Committee on Emerging Systems Toxicology for the Assessment of Risk (eSTAR) are a consortium of private- and public-sector scientists dedicated to developing miRNAs as applied biomarkers. Here, we explore major impediments to routine acceptance and use of miRNA biomarkers and case examples of successes and deficiencies in development. Finally, we provide insight on miRNA measurement, collection, and analysis tools to provide solid footing for addressing knowledge gaps toward routine biomarker use.

Urinary miRNA Biomarkers of Drug-Induced Kidney Injury and Their Site Specificity Within the Nephron.

Drug-induced kidney injury (DIKI) is a major concern in both drug development and clinical practice. There is an unmet need for biomarkers of glomerular damage and more distal renal injury in the loop of Henle and the collecting duct (CD). A cross-laboratory program to identify and characterize urinary microRNA (miRNA) patterns reflecting tissue- or pathology-specific DIKI was conducted. The overall goal was to propose miRNA biomarker candidates for DIKI that could supplement information provided by protein kidney biomarkers in urine. Rats were treated with nephrotoxicants causing injury to distinct nephron segments: the glomerulus, proximal tubule, thick ascending limb (TAL) of the loop of Henle and CD. Meta-analysis identified miR-192-5p as a potential proximal tubule-specific urinary miRNA candidate. This result was supported by data obtained in laser capture microdissection nephron segments showing that miR-192-5p expression was enriched in the proximal tubule. Discriminative miRNAs including miR-221-3p and -222-3p were increased in urine from rats treated with TAL versus proximal tubule toxicants in accordance with their expression localization in the kidney. Urinary miR-210-3p increased up to 40-fold upon treatment with TAL toxicants and was also enriched in laser capture microdissection samples containing TAL and/or CD versus proximal tubule. miR-23a-3p was enriched in the glomerulus and was increased in urine from rats treated with doxorubicin, a glomerular toxicant, but not with toxicants affecting other nephron segments. Taken together these results suggest that urinary miRNA panels sourced from specific nephron regions may be useful to discriminate the pathology of toxicant-induced lesions in the kidney, thereby contributing to DIKI biomarker development needs for industry, clinical, and regulatory use.

Ozone-induced fetal growth restriction in rats is associated with sexually dimorphic placental and fetal metabolic adaptation.

OBJECTIVE: The importance of the placenta in mediating the pre- and post-natal consequences of fetal growth restriction has been increasingly recognized. However, the influence of placental sexual dimorphism on driving these outcomes has received little attention. The purpose of this study was to characterize how sex contributes to the relationship between placental metabolism and fetal programming utilizing a novel rodent model of growth restriction. METHODS: Fetal growth restriction was induced by maternal inhalation of 0.8 ppm ozone (4 h/day) during implantation receptivity (gestation days [GDs] 5 and 6) in Long-Evans rats. Control rats were exposed to filtered air. At GD 21, placental and fetal tissues were obtained for metabolic and genomic assessments. RESULTS: Growth-restricted male placentae exhibited increased mitochondrial biogenesis, increased oxygen consumption, and reduced nutrient storage. Male growth-restricted fetuses also had evidence of reduced adiposity and downregulation of hepatic metabolic signaling. In contrast, placentae from growth-restricted females had elevated markers of autophagy accompanied by an observed protection against hepatic metabolic perturbations. Despite this, growth restriction in females induced a greater number of hypothalamic gene and pathway alterations compared to growth-restricted males. CONCLUSIONS: Increases in mitochondrial metabolism in growth-restricted male placentae likely initiates a sequela of adaptations that promote poor nutrient availability and adiposity. Divergently, the female placenta expresses protective mechanisms that may serve to increase nutrient availability to support fetal metabolic development. Collectively, this work emphasizes the importance of sex in mediating alterations in placental metabolism and fetal programming.

Fish oil and olive oil-enriched diets alleviate acute ozone-induced cardiovascular effects in rats.

Fish oil (FO) and olive oil (OO) supplementations attenuate the cardiovascular responses to inhaled concentrated ambient particles in human volunteers. This study was designed to examine the cardiovascular effects of ozone (O3) exposure and the efficacy of FO and OO-enriched diets in attenuating the cardiovascular effects from O3 exposure in rats. Rats were fed either a normal diet (ND), a diet enriched with 6% FO or OO starting at 4 weeks of age. Eight weeks following the start of these diet, animals were exposed to filtered air (FA) or 0.8 ppm O3, 4 h/day for 2 consecutive days. Immediately after exposure, cardiac function was measured as the indices of left-ventricular developed pressure (LVDP) and contractility (dP/dtmax and dP/dtmin) before ischemia. In addition, selective microRNAs (miRNAs) of inflammation, endothelial function, and cardiac function were assessed in cardiac tissues to examine the molecular alterations of diets and O3 exposure. Pre-ischemic LVDP and dP/dtmax were lower after O3 exposure in rats fed ND but not FO and OO. Cardiac miRNAs expressions were altered by both diet and O3 exposure. Specifically, O3-induced up-regulation of miR-150-5p and miR-208a-5p were attenuated by FO and/or OO. miR-21 was up-regulated by both FO and OO after O3 exposure. This study demonstrated that O3-induced cardiovascular responses appear to be blunted by FO and OO diets. O3-induced alterations in miRNAs linked to inflammation, cardiac function, and endothelial dysfunction support these pathways are involved, and dietary supplementation with FO or OO may alleviate these adverse cardiovascular effects in rats.

Mining a human transcriptome database for chemical modulators of NRF2.

Nuclear factor erythroid-2 related factor 2 (NRF2) encoded by the NFE2L2 gene is a transcription factor critical for protecting cells from chemically-induced oxidative stress. We developed computational procedures to identify chemical modulators of NRF2 in a large database of human microarray data. A gene expression biomarker was built from statistically-filtered gene lists derived from microarray experiments in primary human hepatocytes and cancer cell lines exposed to NRF2-activating chemicals (oltipraz, sulforaphane, CDDO-Im) or in which the NRF2 suppressor Keap1 was knocked down by siRNA. Directionally consistent biomarker genes were further filtered for those dependent on NRF2 using a microarray dataset from cells after NFE2L2 siRNA knockdown. The resulting 143-gene biomarker was evaluated as a predictive tool using the correlation-based Running Fisher algorithm. Using 59 gene expression comparisons from chemically-treated cells with known NRF2 activating potential, the biomarker gave a balanced accuracy of 93%. The biomarker was comprised of many well-known NRF2 target genes (AKR1B10, AKR1C1, NQO1, TXNRD1, SRXN1, GCLC, GCLM), 69% of which were found to be bound directly by NRF2 using ChIP-Seq. NRF2 activity was assessed across ~9840 microarray comparisons from ~1460 studies examining the effects of ~2260 chemicals in human cell lines. A total of 260 and 43 chemicals were found to activate or suppress NRF2, respectively, most of which have not been previously reported to modulate NRF2 activity. Using a NRF2-responsive reporter gene in HepG2 cells, we confirmed the activity of a set of chemicals predicted using the biomarker. The biomarker will be useful for future gene expression screening studies of environmentally-relevant chemicals.

Early microRNA indicators of PPARα pathway activation in the liver.

MicroRNAs (miRNAs) are short non-coding RNA species that play key roles in post-transcriptional regulation of gene expression. MiRNAs also serve as a promising source of early biomarkers for different environmental exposures and health effects, although there is limited information linking miRNA changes to specific target pathways. In this study, we measured liver miRNAs in male B6C3F1 mice exposed to a known chemical activator of the peroxisome proliferator-activated receptor alpha (PPARα) pathway, di(2-ethylhexyl) phthalate (DEHP), for 7 and 28 days at concentrations of 0, 750, 1500, 3000, or 6000 ppm in feed. At the highest dose tested, DEHP altered 61 miRNAs after 7 days and 171 miRNAs after 28 days of exposure, with 48 overlapping miRNAs between timepoints. Analysis of these 48 common miRNAs indicated enrichment in PPARα-related targets and other pathways related to liver injury and cancer. Four of the 10 miRNAs exhibiting a clear dose trend were linked to the PPARα pathway: mmu-miRs-125a-5p, -182-5p, -20a-5p, and -378a-3p. mmu-miRs-182-5p and -378a-3p were subsequently measured using digital drop PCR across a dose range for DEHP and two related phthalates with weaker PPARα activity, di-n-octyl phthalate and n-butyl benzyl phthalate, following 7-day exposures. Analysis of mmu-miRs-182-5p and -378a-3p by transcriptional benchmark dose analysis correctly identified DEHP as having the greatest potency. However, benchmark dose estimates for DEHP based on these miRNAs (average 163; range 126-202 mg/kg-day) were higher on average than values for PPARα target genes (average 74; range 29-183 mg/kg-day). These findings identify putative miRNA biomarkers of PPARα pathway activity and suggest that early miRNA changes may be used to stratify chemical potency.

Transplacental arsenic exposure produced 5-methylcytosine methylation changes and aberrant microRNA expressions in livers of male fetal mice.

Arsenic is a known human carcinogen. Early-life exposure to inorganic arsenic induces tumors in humans and in C3H mice. We hypothesized that arsenic exposure in utero may induce epigenetic changes at the level of DNA methylation and miRNA alterations that could lead to greater postnatal susceptibility to cancer. To test this hypothesis, pregnant C3H mice were given sodium arsenite at doses known to cause liver cancer (42.5 and 85 ppm in the drinking water) from gestation day 8-19, and the livers from male fetal mice were collected for analysis. The antibody against 5-methylcytosine was used to perform chromatin-immunoprecipitation coupled with sequencing (ChIP-Seq) to determine genome-wide methylation alterations. In utero arsenic exposure produced global DNA hypomethylation and an array of gene-specific DNA methylation changes, including hypomethylation of Cyclin D1 and hypermethylation of Tp53. Illumina Correlation Engine analysis revealed 260 methylation alterations that would affect 143 microRNAs. MicroRNA array further revealed 140 aberrantly expressed miRNAs out of the 718 miRNAs. The increased expression of miR-205, miR-203, miR-215, miR-34a, and decreased expression of miR-217 were confirmed by qPCR. Comparison of the methylation changes to those of microarray analyses indicates little if any correspondence between gene methylation and gene expression. The increased expression of Xist, Prrc2, Krit1, Nish, and decreased expression of Prss2, Spp1, Col1a2, and Lox were confirmed by qPCR. In summary, in utero arsenic exposure induced global alterations in DNA methylation and aberrant miRNA expression that might contribute to adult adverse outcomes including liver cancer.

Identification of differentially expressed genes and networks related to hepatic lipid dysfunction.

A range of chemical exposures that resulted in the specific pathology of hepatic lipid dysfunction in rats were selected from DrugMatrix, a publicly available toxicogenomic database. Raw microarray data collected from these exposures were further analyzed using bioinformatic tools to generate a differentially expressed genes (DEGs) dataset associated with hepatic lipid dysfunction. Further analysis of the DEGs dataset resulted in 324 upregulated genes, and 275 genes that were down regulated. Meanwhile, 36 genes were either up regulated or down regulated in different chemical treatments. All identified genes were uploaded in the web application for Database for Annotation, Visualization and Integrated Discovery (DAVID) for gene ontology enrichments and to identify Kyoto Encyclopedia of Genes and Genome (KEGG) pathways. Some of the identified pathways included glycolysis/gluconeogenesis, steroid hormone biosynthesis, retinol metabolism, and metabolism of xenobiotics by cytochrome P450. The same DEGs dataset was also analyzed using Ingenuity Pathway Analysis (IPA) software. IPA identified several pathways including PXR/RXR activation, Aryl hydrocarbon receptor signaling, and xenobiotic metabolism signaling. Furthermore, the generated DEGs lists were uploaded into NCATS BioPlanet platform. Some of the identified pathways were related to fatty acid omega oxidation, lipid and lipoprotein metabolism, and adipogenesis. The enrichment and clarification of the pathways and biological networks obtained from the DEGs dataset provide prior knowledge on the underlying biological key events and molecular mechanisms for the computational development of putative adverse outcome pathways (AOPs) for hepatic lipid dysfunction as a precursor to hepatic steatosis.

Enhanced Quality Metrics for Assessing RNA Derived From Archival Formalin-Fixed Paraffin-Embedded Tissue Samples.

Formalin-fixed paraffin-embedded (FFPE) tissues provide an important resource for toxicogenomic research. However, variability in the integrity or quality of RNA obtained from archival FFPE specimens can lead to unreliable data and wasted resources, and standard protocols for measuring RNA integrity do not adequately assess the suitability of FFPE RNA. The main goal of this study was to identify improved methods for evaluating FFPE RNA quality for whole-genome sequencing. We examined RNA quality metrics conducted prior to RNA-sequencing in paired frozen and FFPE samples with varying levels of quality based on age in block and time in formalin. RNA quality was measured by the RNA integrity number (RIN), a modified RIN called the paraffin-embedded RNA metric, the percentage of RNA fragments >100-300 nucleotides in size (DV100-300), and 2 quantitative PCR-based methods. This information was correlated to sequencing read quality, mapping, and gene detection. Among fragmentation-based methods, DV and PCR-based metrics were more informative than RIN or paraffin-embedded RNA metric in determining sequencing success. Across low- and high-quality FFPE samples, a minimum of 80% of RNA fragments >100 nucleotides (DV100 > 80) provided the best indication of gene diversity and read counts upon sequencing. The PCR-based methods further showed quantitative reductions in amplifiable RNA of target genes related to sample age and time in formalin that inform input quantity of FFPE RNA for sequencing. These results should aid in screening and prioritizing archival FFPE samples for retrospective analyses of gene expression.

Activation of Nrf2 in the liver is associated with stress resistance mediated by suppression of the growth hormone-regulated STAT5b transcription factor.

The transcription factor Nrf2 (encoded by Nfe2l2) induces expression of numerous detoxifying and antioxidant genes in response to oxidative stress. The cytoplasmic protein Keap1 interacts with and represses Nrf2 function. Computational approaches were developed to identify factors that modulate Nrf2 in a mouse liver gene expression compendium. Forty-eight Nrf2 biomarker genes were identified using profiles from the livers of mice in which Nrf2 was activated genetically in Keap1-null mice or chemically by a potent activator of Nrf2 signaling. The rank-based Running Fisher statistical test was used to determine the correlation between the Nrf2 biomarker genes and a test set of 81 profiles with known Nrf2 activation status demonstrating a balanced accuracy of 96%. For a large number of factors examined in the compendium, we found consistent relationships between activation of Nrf2 and feminization of the liver transcriptome through suppression of the male-specific growth hormone (GH)-regulated transcription factor STAT5b. The livers of female mice exhibited higher Nrf2 activation than male mice in untreated or chemical-treated conditions. In male mice, Nrf2 was activated by treatment with ethinyl estradiol, whereas in female mice, Nrf2 was suppressed by treatment with testosterone. Nrf2 was activated in 5 models of disrupted GH signaling containing mutations in Pit1, Prop1, Ghrh, Ghrhr, and Ghr. Out of 59 chemical treatments that activated Nrf2, 36 exhibited STAT5b suppression in the male liver. The Nrf2-STAT5b coupling was absent in in vitro comparisons of chemical treatments. Treatment of male and female mice with 11 chemicals that induce oxidative stress led to activation of Nrf2 to greater extents in females than males. The enhanced basal and inducible levels of Nrf2 activation in females relative to males provides a molecular explanation for the greater resistance often seen in females vs. males to age-dependent diseases and chemical-induced toxicity.

Identification of Androgen Receptor Modulators in a Prostate Cancer Cell Line Microarray Compendium.

High-throughput transcriptomic (HTTr) technologies are increasingly being used to screen environmental chemicals in vitro to identify molecular targets and provide mechanistic context for regulatory testing. Here, we describe the development and validation of a novel gene expression biomarker to identify androgen receptor (AR)-modulating chemicals using a pattern matching method. Androgen receptor biomarker genes were identified by their consistent expression after exposure to 4 AR agonists and 4 AR antagonists and included only those genes that were regulated by AR. The 51 gene biomarker was evaluated as a predictive tool using the fold-change, rank-based Running Fisher algorithm. Using 158 comparisons from cells treated with 95 chemicals, the biomarker gave balanced accuracies for prediction of AR activation or AR suppression of 97% or 98%, respectively. The biomarker correctly classified 16 out of the 17 AR reference antagonists including those that are "weak" and "very weak". Predictions based on microarray profiles from AR-positive LAPC-4 cells treated with 28 chemicals in antagonist mode were compared with those from an AR pathway model which used 11 in vitro HT assays. The balanced accuracy for suppression was 93%. Using our approach, we identified conditions in which AR was modulated in a large collection of microarray profiles from prostate cancer cell lines including (1) constitutively active mutants or knockdown of AR, (2) decreases in availability of androgens by castration or removal from media, and (3) exposure to chemical modulators that work through indirect mechanisms including suppression of AR expression. These results demonstrate that the AR gene expression biomarker could be a useful tool in HTTr to identify AR modulators.