Insights into Repeated Renal Injury Using RNA-Seq with Two New RPTEC Cell Lines

Renal proximal tubule epithelial cells (RPTECs) are a primary site for kidney injury. We created two RPTEC lines from CD-1 mice immortalized with hTERT (human telomerase reverse transcriptase) or SV40 LgT antigen (Simian Virus 40 Large T antigen). Our hypothesis was that low-level, repeated exposure to subcytotoxic levels of 0.25-2.5 μM cisplatin (CisPt) or 12.5-100 μM aflatoxin B1 (AFB1) would activate distinctive genes and pathways in these two differently immortalized cell lines. RNA-seq showed only LgT cells responded to AFB1 with 1139 differentially expressed genes (DEGs) at 72 h. The data suggested that AFB1 had direct nephrotoxic properties on the LgT cells. However, both the cell lines responded to 2.5 μM CisPt from 3 to 96 h expressing 2000-5000 total DEGs. For CisPt, the findings indicated a coordinated transcriptional program of injury signals and repair from the expression of immune receptors with cytokine and chemokine secretion for leukocyte recruitment; robust expression of synaptic and substrate adhesion molecules (SAMs) facilitating the expression of neural and hormonal receptors, ion channels/transporters, and trophic factors; and the expression of nephrogenesis transcription factors. Pathway analysis supported the concept of a renal repair transcriptome. In summary, these cell lines provide in vitro models for the improved understanding of repeated renal injury and repair mechanisms. High-throughput screening against toxicant libraries should provide a wider perspective of their capabilities in nephrotoxicity.

Chronic PFOA exposure in vitro causes acquisition of multiple tumor cell characteristics in rat liver cells

Perfluorooctanoic acid (PFOA) is tumorigenic in rats and mice and potentially tumorigenic in humans. Here, we studied long-term PFOA exposure with an in vitro transformation model using the rat liver epithelial cell, TRL 1215. Cells were cultured in 10 μM (T10), 50 μM (T50) and 100 μM (T100) PFOA for 38 weeks and compared to passage-matched control cells. T100 cells showed morphological changes, loss of cell contact inhibition, formation of multinucleated giant and spindle-shaped cells. T10, T50, and T100 cells showed increased LC50 values 20%, 29% to 35% above control with acute PFOA treatment, indicating a resistance to PFOA toxicity. PFOA-treated cells showed increases in Matrix metalloproteinase-9 secretion, cell migration, and developed more and larger colonies in soft agar. Microarray data showed Myc pathway activation at T50 and T100, associating Myc upregulation with PFOA-induced morphological transformation. Western blot confirmed that PFOA produced significant increases in c-MYC protein expression in a time- and concentration-related manner. Tumor invasion indicators MMP-2 and MMP-9, cell cycle regulator cyclin D1, and oxidative stress protein GST were all significantly overexpressed in T100 cells. Taken together, chronic in vitro PFOA exposure produced multiple cell characteristics of malignant progression and differential gene expression changes suggestive of rat liver cell transformation.

Whole exome and transcript profiling of liver following aflatoxin B1 exposure in rats

We recently developed a rat whole exome sequencing (WES) panel and used it to evaluate early somatic mutations in archival liver tissues from F344/N rats exposed to the hepatocarcinogen, Aflatoxin B1 (AFB1), a widely studied, potent mutagen and hepatocarcinogen associated with hepatocellular carcinoma (HCC). Rats were exposed to 1 ppm AFB1 in feed for 14, 90 and 90 days plus a recovery 60-day, non-exposure period (150-day) time point. Isolated liver DNA was exome sequenced. We identified 172 sequence variants across all time points, of which 101 were non-synonymous variants. Well-annotated genes carried a diverse set of 29 non-synonymous mutations at 14-days, increasing to 39 mutations at 90-days and then decreasing to 33 mutations following the 60-day recovery. Gene Set Enrichment Analysis conducted on previously reported, available RNA expression data of the same exome sequenced archival samples identified altered transcripts in pathways associated with malignant transformation. These included HALLMARK gene sets associated with cell proliferation (MYC Targets Version 1 and Version 2, E2F targets), cell cycle (G2M checkpoint, mitotic spindle), cell death (apoptosis), and DNA damage (DNA repair, UV response Up, Reactive oxygen species) pathways. DriverNet Impact analysis integrated exome-seq and expression data to reveal somatic mutations in Mcm8, Bdp1 and Cct6a that may drive cancer formation. Connectivity with transcript expression changes identified these genes as the top-ranked candidate driver genes associated with hepatocellular transformation. In conclusion, exome sequencing revealed early somatic mutations that may play a role in cancer cell transformation that are translatable to aflatoxin-induced HCC.

Single nucleotide polymorphism patterns associated with a cancer resistant phenotype

BACKGROUND AND AIMS

In this study ten mouse strains representing ~90% of genetic diversity in laboratory mice (B6C3F1/J, C57BL/6J, C3H/HeJ, A/J, NOD.B1oSnH2/J, NZO/HILtJ, 129S1/SvImJ, WSB/EiJ, PWK/PhJ, CAST/EiJ) were examined to identify the mouse strain with the lowest incidence of cancer. The unique single polymorphisms (SNPs) associated with this low cancer incidence are reported.

METHODS

Evaluations of cancer incidence in the 10 mouse strains were based on gross and microscopic diagnosis of tumors. Single nucleotide polymorphisms (SNPs) in the coding regions of the genome were derived from the respective mouse strains located in the Sanger mouse sequencing database and the B6C3F1/N genome from the National Toxicology Program (NTP).

RESULTS

The WSB strain had an overall lower incidence of both benign and malignant tumors compared to the other mouse strains. At 2 years, the incidence of total malignant tumors (Poly-3 incidence rate) ranged from 2% (WSB) to 92% (C3H) in males, and 14% (WSB) to 93% (NZO) in females, and the total incidence of benign and malignant tumor incidence ranged from 13% (WSB) to 99% (C3H) in males and 25% (WSB) to 96% (NOD) in females. Single nucleotide polymorphism (SNP) patterns were examined in the following strains: B6C3F1/N, C57BL/6J, C3H/HeJ, 129S1/SvImJ, A/J, NZO/HILtJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ. We identified 7519 SNPs (involving 5751 Ensembl transcripts of 3453 Ensembl Genes) that resulted in a unique amino acid change in the coding region of the WSB strain.

CONCLUSIONS

The inherited genetic patterns in the WSB cancer-resistant mouse strain occurred in genes involved in multiple cell functions including mitochondria, metabolic, immune, and membrane-related cell functions. The unique SNP patterns in a cancer resistant mouse strain provides insights for understanding and developing strategies for cancer prevention.

Deep Learning Image Analysis of High-Throughput Toxicology Assay Images

High-throughput chemical screening approaches often employ microscopy to capture photomicrographs from multi-well cell culture plates, generating thousands of images that require time-consuming human analysis. To automate this subjective and time-consuming manual process, we have developed a method that uses deep learning to automatically classify digital assay images. We have trained a convolutional neural network (CNN) to perform binary and multi-class classification. The binary classifier binned assay images into healthy (comparable to untreated controls) and altered (not comparable to untreated-control) classes with >98% accuracy; the multi-class classifier assigned "Healthy," "Intermediate" and "Altered" labels to assay images with >95% accuracy. Our dataset comprised high-resolution assay images from primary human hepatocytes and undifferentiated (proliferating) and differentiated 2D cultures of HepaRG cells. In this study we have focused on testing and fine-tuning various CNN architectures, including ResNet 34, 50 and 101. To visualize regions in the images that the CNN model used for classification, we employed Class Activation Maps (CAM). This allowed us to better understand the inner workings of the neural network and led to additional optimizations of the algorithm. The results indicate a strong correspondence between dosage and classifier-predicted scores, suggesting that these scores might be useful in further characterizing benchmark dose. Together, these results clearly demonstrate that deep-learning based automated image classification of cell morphology changes upon chemical-induced stress can yield highly accurate and reproducible assessments of cytotoxicity across a variety of cell types.

Whole genome sequencing of low input circulating cell-free DNA obtained from normal human subjects

Cell-free DNA circulates in plasma at low levels as a normal by-product of cellular apoptosis. Multiple clinical pathologies, as well as environmental stressors can lead to increased circulating cell-free DNA (ccfDNA) levels. Plasma DNA studies frequently employ targeted amplicon deep sequencing platforms due to limited concentrations (ng/ml) of ccfDNA in the blood. Here, we report whole genome sequencing (WGS) and read distribution across chromosomes of ccfDNA extracted from two human plasma samples from normal, healthy subjects, representative of limited clinical samples at <1 ml. Amplification was sufficiently robust with ~90% of the reference genome (GRCh38.p2) exhibiting 10X coverage. Chromosome read coverage was uniform and directly proportional to the number of reads for each chromosome across both samples. Almost 99% of the identified genomic sequence variants were known annotated dbSNP variants in the hg38 reference genome. A high prevalence of C>T and T>C mutations was present along with a strong concordance of variants shared between the germline genome databases; gnomAD (81.1%) and the 1000 Genome Project (93.6%). This study demonstrates isolation and amplification procedures from low input ccfDNA samples that can detect sequence variants across the whole genome from amplified human plasma ccfDNA that can translate to multiple clinical research disciplines.

Utility of Extrapolating Human S1500+ Genes to the Whole Transcriptome: Tunicamycin Case Study

The TempO-Seq S1500+ platform(s), now available for human, mouse, rat, and zebrafish, measures a discrete number of genes that are representative of biological and pathway co-regulation across the entire genome in a given species. While measurement of these genes alone provides a direct assessment of gene expression activity, extrapolating expression values to the whole transcriptome (~26 000 genes in humans) can estimate measurements of non-measured genes of interest and increases the power of pathway analysis algorithms by using a larger background gene expression space. Here, we use data from primary hepatocytes of 54 donors that were treated with the endoplasmic reticulum (ER) stress inducer tunicamycin and then measured on the human S1500+ platform containing ~3000 representative genes. Measurements for the S1500+ genes were then used to extrapolate expression values for the remaining human transcriptome. As a case study of the improved downstream analysis achieved by extrapolation, the “measured only” and “whole transcriptome” (measured + extrapolated) gene sets were compared. Extrapolation increased the number of significant genes by 49%, bringing to the forefront many that are known to be associated with tunicamycin exposure. The extrapolation procedure also correctly identified established tunicamycin-related functional pathways reflected by coordinated changes in interrelated genes while maintaining the sample variability observed from the “measured only” genes. Extrapolation improved the gene- and pathway-level biological interpretations for a variety of downstream applications, including differential expression analysis, gene set enrichment pathway analysis, DAVID keyword analysis, Ingenuity Pathway Analysis, and NextBio correlated compound analysis. The extrapolated data highlight the role of metabolism/metabolic pathways, the ER, immune response, and the unfolded protein response, each of which are key activities associated with tunicamycin exposure that were unrepresented or underrepresented in one or more of the analyses of the original “measured only” dataset. Furthermore, the inclusion of the extrapolated genes raised “tunicamycin” from third to first upstream regulator in Ingenuity Pathway Analysis and from sixth to second most correlated compound in NextBio analysis. Therefore, our case study suggests an approach to extend and enhance data from the S1500+ platform for improved insight into biological mechanisms and functional outcomes of diseases, drugs, and other perturbations.

Quantitative Proteomic Profiling of Mitochondrial Toxicants in a Human Cardiomyocyte Cell Line

Mitochondria are essential cellular organelles that participate in important cellular processes, including bioenergetics, metabolism, and signaling. Recent functional and proteomic studies have revealed the remarkable complexity of mitochondrial protein organization. Mitochondrial protein machineries with diverse functions such as protein translocation, respiration, metabolite transport, protein quality control and the control of membrane architecture interact with each other in dynamic networks. The goal of this study was to identify protein expression changes in a human cardiomyocyte cell line treated with several mitochondrial toxicants which inhibit mitochondrial membrane potential (MMP) and mitochondrial respiration. AC16 human cardiomyocyte cells were treated with carbonyl cyanide p-(trifluoromethoxy)phenylhydrazone (FCCP), dinoterb, picoxystrobin, pinacyanol, and triclocarban for 18 h around the IC₅₀ values generated from MMP assay. The samples were harvested and labeled with tandem mass tags with different mass isotopes. Peptide assignment was performed in Proteome Discoverer. Each dataset was analyzed in Ingenuity Pathway Analysis (IPA). In the proteomic profile, these compounds showed dysregulation of a group of mitochondrial proteins (e.g., NDUA, NDUB, BCS1, CYB5B, and SDHF2), as well as proteins involved in lipid metabolism (e.g., CPT, MECR, and LPGAT1), cytoskeleton protein changes (e.g., CROCC, LAMC3, FBLN1, and FMN2) and stress response (e.g., IKBKG, IKBB, SYVN1, SOD2, and CPIN1). Proteomic data from the current study provides key insights into chemical induced cellular pathway dysregulation, supporting the use of proteomic profiling as a sensitive method to further explore molecular functions and disease pathogenesis upon exposure to environmental chemicals.

SWIFT-Active Screener: Accelerated document screening through active learning and integrated recall estimation

BACKGROUND

In the screening phase of systematic review, researchers use detailed inclusion/exclusion criteria to decide whether each article in a set of candidate articles is relevant to the research question under consideration. A typical review may require screening thousands or tens of thousands of articles in and can utilize hundreds of person-hours of labor.

METHODS

Here we introduce SWIFT-Active Screener, a web-based, collaborative systematic review software application, designed to reduce the overall screening burden required during this resource-intensive phase of the review process. To prioritize articles for review, SWIFT-Active Screener uses active learning, a type of machine learning that incorporates user feedback during screening. Meanwhile, a negative binomial model is employed to estimate the number of relevant articles remaining in the unscreened document list. Using a simulation involving 26 diverse systematic review datasets that were previously screened by reviewers, we evaluated both the document prioritization and recall estimation methods.

RESULTS

On average, 95% of the relevant articles were identified after screening only 40% of the total reference list. In the 5 document sets with 5,000 or more references, 95% recall was achieved after screening only 34% of the available references, on average. Furthermore, the recall estimator we have proposed provides a useful, conservative estimate of the percentage of relevant documents identified during the screening process.

CONCLUSION

SWIFT-Active Screener can result in significant time savings compared to traditional screening and the savings are increased for larger project sizes. Moreover, the integration of explicit recall estimation during screening solves an important challenge faced by all machine learning systems for document screening: when to stop screening a prioritized reference list. The software is currently available in the form of a multi-user, collaborative, online web application.

KRAS-retroviral fusion transcripts and gene amplification in arsenic-transformed, human prostate CAsE-PE cancer cells

CAsE-PE cells are an arsenic-transformed, human prostate epithelial line containing oncogenic mutations in KRAS compared to immortalized, normal KRAS parent cells, RWPE-1. We previously reported increased copy number of mutated KRAS in CAsE-PE cells, suggesting gene amplification. Here, KRAS flanking genomic and transcriptomic regions were sequenced in CAsE-PE cells for insight into KRAS amplification. Comparison of DNA-Seq and RNA-Seq showed increased reads from background aligning to all KRAS exons in CAsE-PE cells, while a uniform DNA-Seq read distribution occurred in RWPE-1 cells with normal transcript expression. We searched for KRAS fusions in DNA and RNA sequencing data finding a portion of reads aligning to KRAS and viral sequence. After generation of cDNA from total RNA, short and long KRAS probes were generated to hybridize cDNA and KRAS enriched fragments were PacBio sequenced. More KRAS reads were captured from CAsE-PE cDNA versus RWPE-1 by each probe set. Only CAsE-PE cDNA showed KRAS viral fusion transcripts, primarily mapping to LTR and endogenous retrovirus sequences on either 5'- or 3'-ends of KRAS. Most KRAS viral fusion transcripts contained 4 to 6 exons but some PacBio sequences were in unusual orientations, suggesting viral insertions within the gene body. Additionally, conditioned media was extracted for potential retroviral particles. RNA-Seq of culture media isolates identified KRAS retroviral fusion transcripts in CAsE-PE media only. Truncated KRAS transcripts suggested multiple retroviral integration sites occurred within the KRAS gene producing KRAS retroviral fusions of various lengths. Findings suggest activation of endogenous retroviruses in arsenic carcinogenesis should be explored.

Next generation sequencing data for use in risk assessment

Next generation sequencing (NGS) represents several powerful platforms that have revolutionized RNA and DNA analysis. The parallel sequencing of millions of DNA molecules can provide mechanistic insights into toxicology and provide new avenues for biomarker discovery with growing relevance for risk assessment. The evolution of NGS technologies has improved over the last decade with increased sensitivity and accuracy to foster new biomarker assays from tissue, blood and other biofluids. NGS sequencing technologies can identify transcriptional changes and genomic targets with base pair precision in response to chemical exposure. Further, there are several exciting movements within the toxicology community that incorporate NGS platforms into new strategies for more rapid toxicological characterizations. These include the Tox21 in vitro high throughput transcriptomic screening program, development of organotypic spheroids, alternative animal models, mining archival tissues, liquid biopsy and epigenomics. This review will describe NGS-based technologies, demonstrate how they can be used as tools for target discovery in tissue and blood, and suggest how they might be applied for risk assessment.

Transcriptome and DNA Methylome Analysis in a Mouse Model of Diet-Induced Obesity Predicts Increased Risk of Colorectal Cancer

Colorectal cancer (CRC) tends to occur at older age; however, CRC incidence rates have been rising sharply among young age groups. The increasing prevalence of obesity is recognized as a major risk, yet the mechanistic underpinnings remain poorly understood. Using a diet-induced obesity mouse model, we identified obesity-associated molecular changes in the colonic epithelium of young and aged mice, and we further investigated whether the changes were reversed after weight loss. Transcriptome analysis indicated that obesity-related colonic cellular metabolic switch favoring long-chain fatty acid oxidation happened in young mice, while obesity-associated downregulation of negative feedback regulators of pro-proliferative signaling pathways occurred in older mice. Strikingly, colonic DNA methylome was pre-programmed by obesity at young age, priming for a tumor-prone gene signature after aging. Furthermore, obesity-related changes were substantially preserved after short-term weight loss, but they were largely reversed after long-term weight loss. We provided mechanistic insights into increased CRC risk in obesity.

Identifying Compounds with Genotoxicity Potential Using Tox21 High-Throughput Screening Assays

Genotoxicity is a critical component of a comprehensive toxicological profile. The Tox21 Program used five quantitative high-throughput screening (qHTS) assays measuring some aspect of DNA damage/repair to provide information on the genotoxic potential of over 10 000 compounds. Included were assays detecting activation of p53, increases in the DNA repair protein ATAD5, phosphorylation of H2AX, and enhanced cytotoxicity in DT40 cells deficient in DNA-repair proteins REV3 or KU70/RAD54. Each assay measures a distinct component of the DNA damage response signaling network; >70% of active compounds were detected in only one of the five assays. When qHTS results were compared with results from three standard genotoxicity assays (bacterial mutation, in vitro chromosomal aberration, and in vivo micronucleus), a maximum of 40% of known, direct-acting genotoxicants were active in one or more of the qHTS genotoxicity assays, indicating low sensitivity. This suggests that these qHTS assays cannot in their current form be used to replace traditional genotoxicity assays. However, despite the low sensitivity, ranking chemicals by potency of response in the qHTS assays revealed an enrichment for genotoxicants up to 12-fold compared with random selection, when allowing a 1% false positive rate. This finding indicates these qHTS assays can be used to prioritize chemicals for further investigation, allowing resources to focus on compounds most likely to induce genotoxic effects. To refine this prioritization process, models for predicting the genotoxicity potential of chemicals that were active in Tox21 genotoxicity assays were constructed using all Tox21 assay data, yielding a prediction accuracy up to 0.83. Data from qHTS assays related to stress-response pathway signaling (including genotoxicity) were the most informative for model construction. By using the results from qHTS genotoxicity assays, predictions from models based on qHTS data, and predictions from commercial bacterial mutagenicity QSAR models, we prioritized Tox21 chemicals for genotoxicity characterization.

Development of a Zebrafish S1500+ Sentinel Gene Set for High-Throughput Transcriptomics

Sentinel gene sets have been developed with the purpose of maximizing the information from targeted transcriptomic platforms. We recently described the development of an S1500+ sentinel gene set, which was built for the human transcriptome, utilizing a data- and knowledge-driven hybrid approach to select a small subset of genes that optimally capture transcriptional diversity, correlation with other genes based on large-scale expression profiling, and known pathway annotation within the human genome. While this detailed bioinformatics approach for gene selection can in principle be applied to other species, the reliability of the resulting gene set depends on availability of a large body of transcriptomics data. For the model organism zebrafish, we aimed to create a similar sentinel gene set (Zf S1500+ gene set); however, there is insufficient standardized expression data in the public domain to train the gene correlation model. Therefore, our strategy was to use human-zebrafish ortholog mapping of the human S1500+ genes and nominations from experts in the zebrafish scientific community. In this study, we present the bioinformatics curation and refinement process to produce the final Zf S1500+ gene set, explore whole transcriptome extrapolation using this gene set, and assess pathway-level inference. This gene set will add value to targeted high-throughput transcriptomics in zebrafish for toxicogenomic screening and other research domains.

BMDExpress 2: enhanced transcriptomic dose-response analysis workflow

SUMMARY

A new version (version 2) of the genomic dose-response analysis software, BMDExpress, has been created. The software addresses the increasing use of transcriptomic dose-response data in toxicology, drug design, risk assessment and translational research. In this new version, we have implemented additional statistical filtering options (e.g. Williams' trend test), curve fitting models, Linux and Macintosh compatibility and support for additional transcriptomic platforms with up-to-date gene annotations. Furthermore, we have implemented extensive data visualizations, on-the-fly data filtering, and a batch-wise analysis workflow. We have also significantly re-engineered the code base to reflect contemporary software engineering practices and streamline future development. The first version of BMDExpress was developed in 2007 to meet an unmet demand for easy-to-use transcriptomic dose-response analysis software. Since its original release, however, transcriptomic platforms, technologies, pathway annotations and quantitative methods for data analysis have undergone a large change necessitating a significant re-development of BMDExpress. To that end, as of 2016, the National Toxicology Program assumed stewardship of BMDExpress. The result is a modernized and updated BMDExpress 2 that addresses the needs of the growing toxicogenomics user community.

AVAILABILITY AND IMPLEMENTATION

BMDExpress 2 is available at https://github.com/auerbachs/BMDExpress-2/releases.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Arsenite malignantly transforms human prostate epithelial cells in vitro by gene amplification of mutated KRAS

Inorganic arsenic is an environmental human carcinogen of several organs including the urinary tract. RWPE-1 cells are immortalized, non-tumorigenic, human prostate epithelia that become malignantly transformed into the CAsE-PE line after continuous in vitro exposure to 5μM arsenite over a period of months. For insight into in vitro arsenite transformation, we performed RNA-seq for differential gene expression and targeted sequencing of KRAS. We report >7,000 differentially expressed transcripts in CAsE-PE cells compared to RWPE-1 cells at >2-fold change, q<0.05 by RNA-seq. Notably, KRAS expression was highly elevated in CAsE-PE cells, with pathway analysis supporting increased cell proliferation, cell motility, survival and cancer pathways. Targeted DNA sequencing of KRAS revealed a mutant specific allelic imbalance, ‘MASI’, frequently found in primary clinical tumors. We found high expression of a mutated KRAS transcript carrying oncogenic mutations at codons 12 and 59 and many silent mutations, accompanied by lower expression of a wild-type allele. Parallel cultures of RWPE-1 cells retained a wild-type KRAS genotype. Copy number analysis and sequencing showed amplification of the mutant KRAS allele. KRAS is expressed as two splice variants, KRAS4a and KRAS4b, where variant 4b is more prevalent in normal cells compared to greater levels of variant 4a seen in tumor cells. 454 Roche sequencing measured KRAS variants in each cell type. We found KRAS4a as the predominant transcript variant in CAsE-PE cells compared to KRAS4b, the variant expressed primarily in RWPE-1 cells and in normal prostate, early passage, primary epithelial cells. Overall, gene expression data were consistent with KRAS-driven proliferation pathways found in spontaneous tumors and malignantly transformed cell lines. Arsenite is recognized as an important environmental carcinogen, but it is not a direct mutagen. Further investigations into this in vitro transformation model will focus on genomic events that cause arsenite-mediated mutation and overexpression of KRAS in CAsE-PE cells.

Identification of Compounds That Inhibit Estrogen-Related Receptor Alpha Signaling Using High-Throughput Screening Assays

The nuclear receptor, estrogen-related receptor alpha (ERRα; NR3B1), plays a pivotal role in energy homeostasis. Its expression fluctuates with the demands of energy production in various tissues. When paired with the peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α), the PGC/ERR pathway regulates a host of genes that participate in metabolic signaling networks and in mitochondrial oxidative respiration. Unregulated overexpression of ERRα is found in many cancer cells, implicating a role in cancer progression and other metabolism-related diseases. Using high throughput screening assays, we screened the Tox21 10K compound library in stably transfected HEK293 cells containing either the ERRα-reporter or the reporter plus PGC-1α expression plasmid. We identified two groups of antagonists that were potent inhibitors of ERRα activity and/or the PGC/ERR pathway: nine antineoplastic agents and thirteen pesticides. Results were confirmed using gene expression studies. These findings suggest a novel mechanism of action on bioenergetics for five of the nine antineoplastic drugs. Nine of the thirteen pesticides, which have not been investigated previously for ERRα disrupting activity, were classified as such. In conclusion, we demonstrated that high-throughput screening assays can be used to reveal new biological properties of therapeutic and environmental chemicals, broadening our understanding of their modes of action.

Mutational analysis of pentabrominated diphenyl-induced hepatocellular tumors in rats and mice, tissue levels of PBDE congeners in rats and mice, and AhR genotyping of Wistar Han rats

This article describes data related to the research article entitled “Carcinogenic activity of pentabrominated diphenyl ether mixture (DE-71) in rats and mice” (Dunnick et al., 2018). PBDE-induced hepatocellular tumors harbored Hras and Ctnnb1 mutations and the methods for these studies are provided. Tissue levels of PBDE congeners in rats and mice after oral exposure to PBDE mixture increased with increasing dose of PBDE. There was no correlation between AhR status and the incidence of hepatocellular tumors in female Wistar Han rats. This manuscript provides additional information on the methods for conducting mutational analysis, PBDE tissue level determinations, and AhR genotyping.

Leucine-rich repeats and calponin homology containing 4 (Lrch4) regulates the innate immune response

Toll-like receptors (TLRs) are pathogen-recognition receptors that trigger the innate immune response. Recent reports have identified accessory proteins that provide essential support to TLR function through ligand delivery and receptor trafficking. Herein, we introduce leucine-rich repeats (LRRs) and calponin homology containing 4 (Lrch4) as a novel TLR accessory protein. Lrch4 is a membrane protein with nine LRRs in its predicted ectodomain. It is widely expressed across murine tissues and has two expression variants that are both regulated by lipopolysaccharide (LPS). Predictive modeling indicates that Lrch4 LRRs conform to the horseshoe-shaped structure typical of LRRs in pathogen-recognition receptors and that the best structural match in the protein database is to the variable lymphocyte receptor of the jawless vertebrate hagfish. Silencing Lrch4 attenuates cytokine induction by LPS and multiple other TLR ligands and dampens the in vivo innate immune response. Lrch4 promotes proper docking of LPS in lipid raft membrane microdomains. We provide evidence that this is through regulation of lipid rafts as Lrch4 silencing reduces cell surface gangliosides, a metric of raft abundance, as well as expression and surface display of CD14, a raft-resident LPS co-receptor. Taken together, we identify Lrch4 as a broad-spanning regulator of the innate immune response and a potential molecular target in inflammatory disease.

Exome Sequencing of Fresh-frozen or Formalin-fixed Paraffin-embedded B6C3F1/N Mouse Hepatocellular Carcinomas Arising Either Spontaneously or due to Chronic Chemical Exposure

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide; however, the mutational properties of HCC-associated carcinogens remain largely uncharacterized. We hypothesized that mechanisms underlying chemical-induced HCC can be characterized by evaluating the mutational spectra of these tumors. To test this hypothesis, we performed exome sequencing of B6C3F1/N HCCs that arose either spontaneously in vehicle controls ( n = 3) or due to chronic exposure to gingko biloba extract (GBE; n = 4) or methyleugenol (MEG; n = 3). Most archived tumor samples are available as formalin-fixed paraffin-embedded (FFPE) blocks, rather than fresh-frozen (FF) samples; hence, exome sequencing from paired FF and FFPE samples was compared. FF and FFPE samples showed 63% to 70% mutation concordance. Multiple known (e.g., Ctnnb1T41A, BrafV637E) and novel (e.g., Erbb4C559S, Card10A700V, and Klf11P358L) mutations in cancer-related genes were identified. The overall mutational burden was greater for MEG than for GBE or spontaneous HCC samples. To characterize the mutagenic mechanisms, we analyzed the mutational spectra in the HCCs according to their trinucleotide motifs. The MEG tumors clustered closest to Catalogue of Somatic Mutations in Cancer signatures 4 and 24, which are, respectively, associated with benzo(a)pyrene- and aflatoxin-induced HCCs in humans. These results establish a novel approach for classifying liver carcinogens and understanding the mechanisms of hepatocellular carcinogenesis.

Carcinogenic activity of pentabrominated diphenyl ether mixture (DE-71) in rats and mice

•

Pentabrominated diphenyl ether (PBDE) mixture was a multispecies carcinogen causing liver tumors in male and female rats and mice.

•

Hras or Ctnnb1 mutations characterized the PBDE-induced liver tumors.

•

PBDE-induced liver tumors increased with increasing PBDE exposure.

Pentabrominated diphenyl ether (PBDE) flame retardants have been phased out in Europe and in the United States, but these lipid soluble chemicals persist in the environment and are found human and animal tissues. PBDEs have limited genotoxic activity. However, in a 2-year cancer study of a PBDE mixture (DE-71) (0, 3, 15, or 50 mg/kg (rats); 0, 3, 30, or 100 mg/kg (mice)) there were treatment-related liver tumors in male and female Wistar Han rats [Crl:WI(Han) after in utero/postnatal/adult exposure, and in male and female B6C3F1 mice, after adult exposure. In addition, there was evidence for a treatment-related carcinogenic effect in the thyroid and pituitary gland tumor in male rats, and in the uterus (stromal polyps/stromal sarcomas) in female rats. The treatment-related liver tumors in female rats were unrelated to the AhR genotype status, and occurred in animals with wild, mutant, or heterozygous Ah receptor. The liver tumors in rats and mice had treatment-related Hras and Ctnnb mutations, respectively. The PBDE carcinogenic activity could be related to oxidative damage, disruption of hormone homeostasis, and molecular and epigenetic changes in target tissue. Further work is needed to compare the PBDE toxic effects in rodents and humans.

Identification of Estrogen-Related Receptor α Agonists in the Tox21 Compound Library

The estrogen-related receptor α (ERRα) is an orphan nuclear receptor (NR) that plays a role in energy homeostasis and controls mitochondrial oxidative respiration. Increased expression of ERRα in certain ovarian, breast, and colon cancers has a negative prognosis, indicating an important role for ERRα in cancer progression. An interaction between ERRα and peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α) has also recently been shown to regulate an enzyme in the β-oxidation of free fatty acids, thereby suggesting that ERRα plays an important role in obesity and type 2 diabetes. Therefore, it would be prudent to identify compounds that can act as activators of ERRα. In this study, we screened ∼10,000 (8311 unique) compounds, known as the Tox21 10K collection, to identify agonists of ERRα. We performed this screen using two stably transfected HEK 293 cell lines, one with the ERRα-reporter alone and the other with both ERRα-reporter and PGC-1α expression vectors. After the primary screening, we identified more than five agonist clusters based on compound structural similarity analysis (e.g., statins). By examining the activities of the confirmed ERRα modulators in other Tox21 NR assays, eliminating those with promiscuous NR activity, and performing follow-up assays (e.g., small interfering RNA knockdown), we identified compounds that might act as endocrine disrupters through effects on ERRα signaling. To our knowledge, this study is the first comprehensive analysis in discovering potential endocrine disrupters that affect the ERRα signaling pathway.

HAfTs are novel lncRNA transcripts from aflatoxin exposure

The transcriptome can reveal insights into precancer biology. We recently conducted RNA-Seq analysis on liver RNA from male rats exposed to the carcinogen, aflatoxin B1 (AFB1), for 90 days prior to liver tumor onset. Among >1,000 differentially expressed transcripts, several novel, unannotated Cufflinks-assembled transcripts, or HAfTs (Hepatic Aflatoxin Transcripts) were found. We hypothesized PCR-cloning and RACE (rapid amplification of cDNA ends) could further HAfT identification. Sanger data was obtained for 6 transcripts by PCR and 16 transcripts by 5’- and 3’-RACE. BLAST alignments showed, with two exceptions, HAfT transcripts were lncRNAs, >200nt without apparent long open reading frames. Six rat HAfT transcripts were classified as ‘novel’ without RefSeq annotation. Sequence alignment and genomic synteny showed each rat lncRNA had a homologous locus in the mouse genome and over half had homologous loci in the human genome, including at least two loci (and possibly three others) that were previously unannotated. While HAfT functions are not yet clear, coregulatory roles may be possible from their adjacent orientation to known coding genes with altered expression that include 8 HAfT-gene pairs. For example, a unique rat HAfT, homologous to Pvt1, was adjacent to known genes controlling cell proliferation. Additionally, PCR and RACE Sanger sequencing showed many alternative splice variants and refinements of exon sequences compared to Cufflinks assembled transcripts and gene prediction algorithms. Presence of multiple splice variants and short tandem repeats found in some HAfTs may be consequential for secondary structure, transcriptional regulation, and function. In summary, we report novel, differentially expressed lncRNAs after exposure to the genotoxicant, AFB1, prior to neoplastic lesions. Complete cloning and sequencing of such transcripts could pave the way for a new set of sensitive and early prediction markers for chemical hepatocarcinogens.

Development of Novel Cell Lines for High-Throughput Screening to Detect Estrogen-Related Receptor Alpha Modulators

Estrogen-related receptor alpha (ERRα), the first orphan nuclear receptor discovered, is crucial for the control of cellular energy metabolism. ERRα and its coactivator, peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC-1α), are required for rapid energy production in response to environmental challenges. They have been implicated in the etiology of metabolic disorders such as type 2 diabetes and metabolic syndrome. ERRα also plays a role in the pathogenesis of breast cancer. Identification of compounds that modulate ERRα signaling may elucidate environmental factors associated with these diseases. Therefore, we developed stable cell lines containing an intact ERRα signaling pathway, with and without the coactivator PGC-1α, to use as high-throughput screening tools to detect ERRα modulators. The lentiviral PGC-1α expression constructs and ERRα multiple hormone response element (MHRE) reporters were introduced into HEK293T cells that express endogenous ERRα. A cell line expressing the reporter alone was designated "ERR." A second cell line expressing both reporter and PGC-1α was named "PGC/ERR." Initial screenings of the Library of Pharmacologically Active Compounds (LOPAC) identified 33 ERR and 22 PGC/ERR agonists, and 54 ERR and 15 PGC/ERR antagonists. Several potent ERRα agonists were dietary plant compounds (e.g., genistein). In conclusion, these cell lines are suitable for high-throughput screens to identify environmental chemicals affecting metabolic pathways and breast cancer progression.

Tetrabromobisphenol A activates the hepatic interferon pathway in rats

Tetrabromobisphenol A (TBBPA) is a widely used flame retardant in printed circuit boards, paper, and textiles. In a two-year study, TBBPA showed evidence of uterine tumors in female Wistar-Han rats and liver and colon tumors in B6C3F1 mice. In order to gain further insight into early gene and pathway changes leading to cancer, we exposed female Wistar Han rats to TBBPA at 0, 25, 250, or 1000mg/kg (oral gavage in corn oil, 5×/week) for 13 weeks. Because at the end of the TBBPA exposure period, there were no treatment-related effects on body weights, liver or uterus lesions, and liver and uterine organ weights were within 10% of controls, only the high dose animals were analyzed. Analysis of the hepatic and uterine transcriptomes showed TBBPA-induced changes primarily in the liver (1000mg/kg), with 159 transcripts corresponding to 132 genes differentially expressed compared to controls (FDR=0.05). Pathway analysis showed activation of interferon (IFN) and metabolic networks. TBBPA induced few molecular changes in the uterus. Activation of the interferon pathway in the liver occurred after 13-weeks of TBBPA exposure, and with longer term TBBPA exposure this may lead to immunomodulatory changes that contribute to carcinogenic processes.

Invited Review: Toxicoproteomics: Proteomics Applied to Toxicology and Pathology

Global measurement of proteins and their many attributes in tissues and biofluids defines the field of proteomics. Toxicoproteomics, as part of the larger field of toxicogenomics, seeks to identify critical proteins and pathways in biological systems that are affected by and respond to adverse chemical and environmental exposures using global protein expression technologies. Toxicoproteomics integrates 3 disciplinary areas: traditional toxicology and pathology, differential protein and gene expression analysis, and systems biology. Key topics to be reviewed are the evolution of proteomics, proteomic technology platforms and their capabilities with exemplary studies from biology and medicine, a review of over 50 recent studies applying proteomic analysis to toxicological research, and the recent development of databases designed to integrate -Omics technologies with toxicology and pathology. Proteomics is examined for its potential in discovery of new biomarkers and toxicity signatures, in mapping serum, plasma, and other biofluid proteomes, and in parallel proteomic and transcriptomic studies. The new field of toxicoproteomics is uniquely positioned toward an expanded understanding of protein expression during toxicity and environmental disease for the advancement of public health.

Molecular Changes in the Nasal Cavity after N, N-dimethyl-p-toluidine Exposure

N, N-dimethyl-p-toluidine (DMPT; Cas No. 99-97-8), an accelerant for methyl methacrylate monomers in medical devices, is a nasal cavity carcinogen according to a 2-yr cancer study of male and female F344/N rats, with the nasal tumors arising from the transitional cell epithelium. In this study, we exposed male F344/N rats for 5 days to DMPT (0, 1, 6, 20, 60, or 120 mg/kg [oral gavage]) to explore the early changes in the nasal cavity after short-term exposure. Lesions occurred in the nasal cavity including hyperplasia of transitional cell epithelium (60 and 120 mg/kg). Nasal tissue was rapidly removed and preserved for subsequent laser capture microdissection and isolation of the transitional cell epithelium (0 and 120 mg/kg) for transcriptomic studies. DMPT transitional cell epithelium gene transcript patterns were characteristic of an antioxidative damage response (e.g., Akr7a3, Maff, and Mgst3), cell proliferation, and decrease in signals for apoptosis. The transcripts of amino acid transporters were upregulated (e.g., Slc7a11). The DMPT nasal transcript expression pattern was similar to that found in the rat nasal cavity after formaldehyde exposure, with over 1,000 transcripts in common. Molecular changes in the nasal cavity after DMPT exposure suggest that oxidative damage is a mechanism of the DMPT toxic and/or carcinogenic effects.

Differential DNA methylation profile of key genes in malignant prostate epithelial cells transformed by inorganic arsenic or cadmium

Previous work shows altered methylation patterns in inorganic arsenic (iAs)- or cadmium (Cd)-transformed epithelial cells. Here, the methylation status near the transcriptional start site was assessed in the normal human prostate epithelial cell line (RWPE-1) that was malignantly transformed by 10μM Cd for 11weeks (CTPE) or 5μM iAs for 29weeks (CAsE-PE), at which time cells showed multiple markers of acquired cancer phenotype. Next generation sequencing of the transcriptome of CAsE-PE cells identified multiple dysregulated genes. Of the most highly dysregulated genes, five genes that can be relevant to the carcinogenic process (S100P, HYAL1, NTM, NES, ALDH1A1) were chosen for an in-depth analysis of the DNA methylation profile. DNA was isolated, bisulfite converted, and combined bisulfite restriction analysis was used to identify differentially methylated CpG sites, which was confirmed with bisulfite sequencing. Four of the five genes showed differential methylation in transformants relative to control cells that was inversely related to altered gene expression. Increased expression of HYAL1 (>25-fold) and S100P (>40-fold) in transformants was correlated with hypomethylation near the transcriptional start site. Decreased expression of NES (>15-fold) and NTM (>1000-fold) in transformants was correlated with hypermethylation near the transcriptional start site. ALDH1A1 expression was differentially expressed in transformed cells but was not differentially methylated relative to control. In conclusion, altered gene expression observed in Cd and iAs transformed cells may result from altered DNA methylation status.

Cell-Based High-Throughput Screening for Aromatase Inhibitors in the Tox21 10K Library

Multiple mechanisms exist for endocrine disruption; one nonreceptor-mediated mechanism is via effects on aromatase, an enzyme critical for maintaining the normal in vivo balance of androgens and estrogens. We adapted the AroER tri-screen 96-well assay to 1536-well format to identify potential aromatase inhibitors (AIs) in the U.S. Tox21 10K compound library. In this assay, screening with compound alone identifies estrogen receptor alpha (ERα) agonists, screening in the presence of testosterone (T) identifies AIs and/or ERα antagonists, and screening in the presence of 17β-estradiol (E2) identifies ERα antagonists. Screening the Tox-21 library in the presence of T resulted in finding 302 potential AIs. These compounds, along with 31 known AI actives and inactives, were rescreened using all 3 assay formats. Of the 333 compounds tested, 113 (34%; 63 actives, 50 marginal actives) were considered to be potential AIs independent of cytotoxicity and ER antagonism activity. Structure-activity analysis suggested the presence of both conventional (eg, 1, 2, 4, - triazole class) and novel AI structures. Due to their novel structures, 14 of the 63 potential AI actives, including both drugs and fungicides, were selected for confirmation in the biochemical tritiated water-release aromatase assay. Ten compounds were active in the assay; the remaining 4 were only active in high-throughput screen assay, but with low efficacy. To further characterize these 10 novel AIs, we investigated their binding characteristics. The AroER tri-screen, in high-throughput format, accurately and efficiently identified chemicals in a large and diverse chemical library that selectively interact with aromatase.

Gene expression in obliterative bronchiolitis-like lesions in 2,3-pentanedione-exposed rats

Obliterative bronchiolitis (OB) is an irreversible lung disease characterized by progressive fibrosis in the small airways with eventual occlusion of the airway lumens. OB is most commonly associated with lung transplant rejection; however, OB has also been diagnosed in workers exposed to artificial butter flavoring (ABF) vapors. Research has been limited by the lack of an adequate animal model of OB, and as a result the mechanism(s) is unclear and there are no effective treatments for this condition. Exposure of rats to the ABF component, 2,3-pentanedione (PD) results in airway lesions that are histopathologically similar to those in human OB. We used this animal model to evaluate changes in gene expression in the distal bronchi of rats with PD-induced OB. Male Wistar Han rats were exposed to 200 ppm PD or air 6 h/d, 5 d/wk for 2-wks. Bronchial tissues were laser microdissected from serial sections of frozen lung. In exposed lungs, both fibrotic and non-fibrotic airways were collected. Following RNA extraction and microarray analysis, differential gene expression was evaluated. In non-fibrotic bronchi of exposed rats, 4683 genes were significantly altered relative to air-exposed controls with notable down-regulation of many inflammatory cytokines and chemokines. In contrast, in fibrotic bronchi, 3807 genes were significantly altered with a majority of genes being up-regulated in affected pathways. Tgf-β2 and downstream genes implicated in fibrosis were significantly up-regulated in fibrotic lesions. Genes for collagens and extracellular matrix proteins were highly up-regulated. In addition, expression of genes for peptidases and peptidase inhibitors were significantly altered, indicative of the tissue remodeling that occurs during airway fibrosis. Our data provide new insights into the molecular mechanisms of OB. This new information is of potential significance with regard to future therapeutic targets for treatment.

Intersection of toxicogenomics and high throughput screening in the Tox21 program: an NIEHS perspective

Humans are exposed to thousands of chemicals with inadequate toxicological data. Advances in computational toxicology, robotic high throughput screening (HTS), and genome-wide expression have been integrated into the Tox21 program to better predict the toxicological effects of chemicals. Tox21 is a collaboration among US government agencies initiated in 2008 that aims to shift chemical hazard assessment from traditional animal toxicology to target-specific, mechanism-based, biological observations using in vitro assays and lower organism models. HTS uses biocomputational methods for probing thousands of chemicals in in vitro assays for gene-pathway response patterns predictive of adverse human health outcomes. In 1999, NIEHS began exploring the application of toxicogenomics to toxicology and recent advances in NextGen sequencing should greatly enhance the biological content obtained from HTS platforms. We foresee an intersection of new technologies in toxicogenomics and HTS as an innovative development in Tox21. Tox21 goals, priorities, progress, and challenges will be reviewed.

Development of a stable cell line with an intact PGC-1α/ERRα axis for screening environmental chemicals

The estrogen-related receptor α (ERRα) and the peroxisome proliferator-activated receptor γ (PPARγ) coactivator 1α (PGC-1α) play critical roles in the control of several physiological functions, including the regulation of genes involved in energy homeostasis. However, little is known about the ability of environmental chemicals to disrupt or modulate this important bioenergetics pathway in humans. The goal of this study was to develop a cell-based assay system with an intact PGC-1α/ERRα axis that could be used as a screening assay for detecting such chemicals. To this end, we successfully generated several stable cell lines expressing PGC-1α and showed that the reporter driven by the native ERRα hormone response unit (AAB-Luc) is active in these cell lines and that the activation is PGC-1α-dependent. Furthermore, we show that this activation can be blocked by the ERRα selective inverse agonist, XCT790. In addition, we find that genistein and bisphenol A further stimulate the reporter activity, while kaempferol has minimal effect. These cell lines will be useful for identifying environmental chemicals that modulate this important pathway.

In vivo molecular markers for pro-inflammatory cytokine M1 stage and resident microglia in trimethyltin-induced hippocampal injury

Microglia polarization to the classical M1 activation state is characterized by elevated pro-inflammatory cytokines; however, a full profile has not been generated in the early stages of a sterile inflammatory response recruiting only resident microglia. We characterized the initial M1 state in a hippocampal injury model dependent upon tumor necrosis factor (TNF) receptor signaling for dentate granule cell death. Twenty-one-day-old CD1 male mice were injected with trimethyltin (TMT 2.3 mg/kg, i.p.) and the hippocampus was examined at an early stage (24-h post-dosing) of neuronal death. Glia activation was assessed using a custom quantitative nuclease protection assay. We report elevated mRNA levels for glia response such as ionizing calcium-binding adapter molecule-1 and glial fibrillary acidic protein (Gfap); Fas, hypoxia inducible factor alpha, complement component 1qb, TNF-related genes (Tnf, Tnfaip3, Tnfrsfla); interleukin-1 alpha, Cd44, chemokine (C-C motif) ligand (Ccl)2, Cc14, integrin alpha M, lipocalin (Lcn2), and secreted phosphoprotein 1 (Spp1). These changes occurred in the absence of changes in matrix metalloproteinase 9 and 12, neural cell adhesion molecule, metabotropic glutamate receptor (Grm)3, and Ly6/neurotoxin 1 (Lynx1), as well as, a decrease in neurotrophin 3, glutamate receptor subunit epsilon (Grin)-2b, and neurotrophic tyrosine kinase receptor, type 3. The M2 anti-inflammatory marker, transforming growth factor beta-1 (Tgfb1) was elevated. mRNAs associated with early stage of injury-induced neurogenesis including fibroblast growth factor 21 and Mki67 were elevated. In the "non-injured" temporal cortex receiving projections from the hippocampus, Lynx1, Grm3, and Grin2b were decreased and Gfap increased. Formalin fixed-paraffin-embedded tissue did not generate a comparable profile.

p53 Integrates host defense and cell fate during bacterial pneumonia

p53 deletion augments neutrophil-mediated bacterial clearance in the lung at the expense of tissue homeostasis, leading to increased mortality.

Cancer and infection are predominant causes of human mortality and derive, respectively, from inadequate genomic and host defenses against environmental agents. The transcription factor p53 plays a central role in human tumor suppression. Despite its expression in immune cells and broad responsiveness to stressors, it is virtually unknown whether p53 regulates host defense against infection. We report that the lungs of naive p53^−/− mice display genome-wide induction of NF-κB response element–enriched proinflammatory genes, suggestive of type 1 immune priming. p53-null and p53 inhibitor–treated mice clear Gram-negative and -positive bacteria more effectively than controls after intrapulmonary infection. This is caused, at least in part, by cytokines produced by an expanded population of apoptosis-resistant, TLR-hyperresponsive alveolar macrophages that enhance airway neutrophilia. p53^−/− neutrophils, in turn, display heightened phagocytosis, Nox-dependent oxidant generation, degranulation, and bacterial killing. p53 inhibition boosts bacterial killing by mouse neutrophils and oxidant generation by human neutrophils. Despite enhanced bacterial clearance, infected p53^−/− mice suffer increased mortality associated with aggravated lung injury. p53 thus modulates host defense through regulating microbicidal function and fate of phagocytes, revealing a fundamental link between defense of genome and host during environmental insult.

RNA-Seq profiling reveals novel hepatic gene expression pattern in aflatoxin B1 treated rats

Deep sequencing was used to investigate the subchronic effects of 1 ppm aflatoxin B1 (AFB1), a potent hepatocarcinogen, on the male rat liver transcriptome prior to onset of histopathological lesions or tumors. We hypothesized RNA-Seq would reveal more differentially expressed genes (DEG) than microarray analysis, including low copy and novel transcripts related to AFB1’s carcinogenic activity compared to feed controls (CTRL). Paired-end reads were mapped to the rat genome (Rn4) with TopHat and further analyzed by DESeq and Cufflinks-Cuffdiff pipelines to identify differentially expressed transcripts, new exons and unannotated transcripts. PCA and cluster analysis of DEGs showed clear separation between AFB1 and CTRL treatments and concordance among group replicates. qPCR of eight high and medium DEGs and three low DEGs showed good comparability among RNA-Seq and microarray transcripts. DESeq analysis identified 1,026 differentially expressed transcripts at greater than two-fold change (p<0.005) compared to 626 transcripts by microarray due to base pair resolution of transcripts by RNA-Seq, probe placement within transcripts or an absence of probes to detect novel transcripts, splice variants and exons. Pathway analysis among DEGs revealed signaling of Ahr, Nrf2, GSH, xenobiotic, cell cycle, extracellular matrix, and cell differentiation networks consistent with pathways leading to AFB1 carcinogenesis, including almost 200 upregulated transcripts controlled by E2f1-related pathways related to kinetochore structure, mitotic spindle assembly and tissue remodeling. We report 49 novel, differentially-expressed transcripts including confirmation by PCR-cloning of two unique, unannotated, hepatic AFB1-responsive transcripts (HAfT’s) on chromosomes 1.q55 and 15.q11, overexpressed by 10 to 25-fold. Several potentially novel exons were found and exon refinements were made including AFB1 exon-specific induction of homologous family members, Ugt1a6 and Ugt1a7c. We find the rat transcriptome contains many previously unidentified, AFB1-responsive exons and transcripts supporting RNA-Seq’s capabilities to provide new insights into AFB1-mediated gene expression leading to hepatocellular carcinoma.

Testing an aflatoxin B1 gene signature in rat archival tissues

Archival tissues from laboratory studies represent a unique opportunity to explore the relationship between genomic changes and agent-induced disease. In this study, we evaluated the applicability of qPCR for detecting genomic changes in formalin-fixed, paraffin-embedded (FFPE) tissues by determining if a subset of 14 genes from a 90-gene signature derived from microarray data and associated with eventual tumor development could be detected in archival liver, kidney, and lung of rats exposed to aflatoxin B1 (AFB1) for 90 days in feed at 1 ppm. These tissues originated from the same rats used in the microarray study. The 14 genes evaluated were Adam8, Cdh13, Ddit4l, Mybl2, Akr7a3, Akr7a2, Fhit, Wwox, Abcb1b, Abcc3, Cxcl1, Gsta5, Grin2c, and the C8orf46 homologue. The qPCR FFPE liver results were compared to the original liver microarray data and to qPCR results using RNA from fresh frozen liver. Archival liver paraffin blocks yielded 30 to 50 μg of degraded RNA that ranged in size from 0.1 to 4 kB. qPCR results from FFPE and fresh frozen liver samples were positively correlated (p ≤ 0.05) by regression analysis and showed good agreement in direction and proportion of change with microarray data for 11 of 14 genes. All 14 transcripts could be amplified from FFPE kidney RNA except the glutamate receptor gene Grin2c; however, only Abcb1b was significantly upregulated from control. Abundant constitutive transcripts, S18 and β-actin, could be amplified from lung FFPE samples, but the narrow RNA size range (25-500 bp length) prevented consistent detection of target transcripts. Overall, a discrete gene signature derived from prior transcript profiling and representing cell cycle progression, DNA damage response, and xenosensor and detoxication pathways was successfully applied to archival liver and kidney by qPCR and indicated that gene expression changes in response to subchronic AFB1 exposure occurred predominantly in the liver, the primary target for AFB1-induced tumors. We conclude that an evaluation of gene signatures in archival tissues can be an important toxicological tool for evaluating critical molecular events associated with chemical exposures.

Proteomic profiling of S-acylated macrophage proteins identifies a role for palmitoylation in mitochondrial targeting of phospholipid scramblase 3

S-Palmitoylation, the reversible post-translational acylation of specific cysteine residues with the fatty acid palmitate, promotes the membrane tethering and subcellular localization of proteins in several biological pathways. Although inhibiting palmitoylation holds promise as a means for manipulating protein targeting, advances in the field have been hampered by limited understanding of palmitoylation enzymology and consensus motifs. In order to define the complement of S-acylated proteins in the macrophage, we treated RAW 264.7 macrophage membranes with hydroxylamine to cleave acyl thioesters, followed by biotinylation of newly exposed sulfhydryls and streptavidin-agarose affinity chromatography. Among proteins identified by LC-MS/MS, S-acylation status was established by spectral counting to assess enrichment under hydroxylamine versus mock treatment conditions. Of 1183 proteins identified in four independent experiments, 80 proteins were significant for S-acylation at false discovery rate = 0.05, and 101 significant at false discovery rate = 0.10. Candidate S-acylproteins were identified from several functional categories, including membrane trafficking, signaling, transporters, and receptors. Among these were 29 proteins previously biochemically confirmed as palmitoylated, 45 previously reported as putative S-acylproteins in proteomic screens, 24 not previously associated with palmitoylation, and three presumed false-positives. Nearly half of the candidates were previously identified by us in macrophage detergent-resistant membranes, suggesting that palmitoylation promotes lipid raft-localization of proteins in the macrophage. Among the candidate novel S-acylproteins was phospholipid scramblase 3 (Plscr3), a protein that regulates apoptosis through remodeling the mitochondrial membrane. Palmitoylation of Plscr3 was confirmed through (3)H-palmitate labeling. Moreover, site-directed mutagenesis of a cluster of five cysteines (Cys159-161-163-164-166) abolished palmitoylation, caused Plscr3 mislocalization from mitochondrion to nucleus, and reduced macrophage apoptosis in response to etoposide, together suggesting a role for palmitoylation at this site for mitochondrial targeting and pro-apoptotic function of Plscr3. Taken together, we propose that manipulation of protein palmitoylation carries great potential for intervention in macrophage biology via reprogramming of protein localization.

Platforms for biomarker analysis using high-throughput approaches in genomics, transcriptomics, proteomics, metabolomics, and bioinformatics

Global biological responses that reflect disease or exposure biology are kinetic and highly dynamic phenomena. While high-throughput DNA sequencing continues to drive genomics, the possibility of more broadly measuring changes in gene expression has been a recent development manifested by a diversity of technical platforms. Such technologies measure transcripts, proteins and small biological molecules, or metabolites, and respectively define the fields of transcriptomics, proteomics and metabolomics that can be performed at a cell-, tissue-, or organism-wide basis. Bioinformatics is the discipline that derives knowledge from the large quantity and diversity of biological, genetic, genomic and gene expression data by integrating computer science, mathematics, statistics and graphic arts. Gene, protein and metabolite expression profiles can be thought of as snapshots of the current, poorly-mapped molecular landscape. The ultimate aim of genomic platforms is to fully map this landscape to more completely describe all of the biological interactions within a living system, during disease and toxicity, and define the behaviour and relationships of all the components of a biological system. The development of databases and knowledge bases will support the integration of data from multiple domains, as well as computational modelling. This chapter will describe the technical platform methods involving DNA sequencing, mass spectrometry, nuclear magnetic resonance combined with separation systems, and bioinformatics to derive genomic and gene expression data and include the relevant bioinformatic tools for analysis. These genomic, or omics platforms should have wide application to epidemiological studies.

Myeloid differentiation primary response protein 88 couples reverse cholesterol transport to inflammation

Crosstalk exists in mammalian cells between cholesterol trafficking and innate immune signaling. Apolipoprotein A-I (apoA-I), a serum apolipoprotein that induces antiatherogenic efflux of macrophage cholesterol, is widely described as anti-inflammatory because it neutralizes bacterial lipopolysaccharide. Conversely, lipopolysaccharide-induced inflammation is proatherogenic. However, whether innate immunity plays an endogenous, physiological role in host cholesterol homeostasis in the absence of infection is undetermined. We report that apoA-I signals in the macrophage through Toll-like receptor (TLR)2, TLR4, and CD14, utilizing myeloid differentiation primary response protein 88 (MyD88)-dependent and -independent pathways, to activate nuclear factor-kappaB and induce cytokines. MyD88 plays a critical role in reverse cholesterol transport in vitro and in vivo, in part through promoting ATP-binding cassette A1 transporter upregulation. Taken together, this work identifies apoA-I as an endogenous stimulus of innate immunity that couples cholesterol trafficking to inflammation through MyD88 and identifies innate immunity as a physiologic signal in cholesterol homeostasis.

The role of toxicoproteomics in assessing organ specific toxicity

Aims of this chapter on the role of toxicoproteomics in assessing organ-specific toxicity are to define the field of toxicoproteomics, describe its development among global technologies, and show potential uses in experimental toxicological research, preclinical testing and mechanistic biological research. Disciplines within proteomics deployed in preclinical research are described as Tier I analysis, involving global protein mapping and protein profiling for differential expression, and Tier II proteomic analysis, including global methods for description of function, structure, interactions and post-translational modification of proteins. Proteomic platforms used in toxicoproteomics research are briefly reviewed. Preclinical toxicoproteomic studies with model liver and kidney toxicants are critically assessed for their contributions toward understanding pathophysiology and in biomarker discovery. Toxicoproteomics research conducted in other organs and tissues are briefly discussed as well. The final section suggests several key developments involving new approaches and research focus areas for the field of toxicoproteomics as a new tool for toxicological pathology.

Quantitative proteomics analysis of macrophage rafts reveals compartmentalized activation of the proteasome and of proteasome-mediated ERK activation in response to lipopolysaccharide

Lipopolysaccharide (LPS), a glycolipid component of the outer membrane of Gram-negative bacteria, is a potent initiator of the innate immune response of the macrophage. LPS triggers downstream signaling by selectively recruiting and activating proteins in cholesterol-rich membrane microdomains called lipid rafts. We applied proteomics analysis to macrophage detergent-resistant membranes (DRMs) during an LPS exposure time course in an effort to identify and validate novel events occurring in macrophage rafts. Following metabolic incorporation in cell culture of heavy isotopes of amino acids arginine and lysine ([(13)C(6)]Arg and [(13)C(6)]Lys) or their light counterparts, a SILAC (stable isotope labeling with amino acids in cell culture)-based quantitative, liquid chromatography-tandem mass spectrometry proteomics approach was used to profile LPS-induced changes in the lipid raft proteome of RAW 264.7 macrophages. Unsupervised network analysis of the proteomics data set revealed a marked representation of the ubiquitin-proteasome system as well as changes in proteasome subunit composition following LPS challenge. Functional analysis of DRMs confirmed that LPS causes selective activation of the proteasome in macrophage rafts and proteasome inactivation outside of rafts. Given previous reports of an essential role for proteasomal degradation of IkappaB kinase-phosphorylated p105 in LPS activation of ERK mitogen-activated protein kinase, we tested for a role of rafts in compartmentalization of these events. Immunoblotting of DRMs revealed proteasome-dependent activation of MEK and ERK specifically occurring in lipid rafts as well as proteasomal activity upon raft-localized p105 that was enhanced by LPS. Cholesterol extraction from the intact macrophage with methyl-beta-cyclodextrin was sufficient to activate ERK, recapitulating the LPS-IkappaB kinase-p105-MEK-ERK cascade, whereas both it and the alternate raft-disrupting agent nystatin blocked subsequent LPS activation of the ERK cascade. Taken together, our findings indicate a critical, selective role for raft compartmentalization and regulation of proteasome activity in activation of the MEK-ERK pathway.

Complementary gene and protein expression studies and integrative approaches in toxicogenomics

Parallel transcript and protein profiling is a strategy to gain further insight into the mechanisms of toxicity and disease. The technologies used to measure expression at the transcript and protein levels each convey different information and have different technical capabilities that can complement each other when combined. In this review, over twenty studies are considered for the use of -Omics platform, the chemical or disease being profiled, tissues, the number of genes and proteins found by each platform and common expression products. A strategy is suggested for toxicant expression profiling that combines the transcriptomics and proteomics of both the target tissue and blood/serum that could provide a more complete characterization of toxicity as well as synergize expression technologies toward biomarker discovery.

Gene expression response in target organ and whole blood varies as a function of target organ injury phenotype

This report details the standardized experimental design and the different data streams that were collected (histopathology, clinical chemistry, hematology and gene expression from the target tissue (liver) and a bio-available tissue (blood)) after treatment with eight known hepatotoxicants (at multiple time points and doses with multiple biological replicates). The results of the study demonstrate the classification of histopathological differences, likely reflecting differences in mechanisms of cell-specific toxicity, using either liver tissue or blood transcriptomic data.

The plasma proteome, adductome and idiosyncratic toxicity in toxicoproteomics research

Toxicoproteomics uses the discovery potential of proteomics in toxicology research by applying global protein measurement technologies to biofluids and tissues after host exposure to injurious agents. Toxicoproteomic studies thus far have focused on protein profiling of major organs and biofluids such as liver and blood in preclinical species exposed to model toxicants. The slow pace of discovery for new biomarkers, toxicity signatures and mechanistic insights is partially due to the limited proteome coverage derived from analysis of native organs, tissues and body fluids by traditional proteomic platforms. Improved toxicoproteomic analysis would result by combining higher data density LC-MS/MS platforms with stable isotope labelled peptides and parallel use of complementary platforms. Study designs that remove abundant proteins from biofluids, enrich subcellular structures and include cell specific isolation from heterogeneous tissues would greatly increase differential expression capabilities. By leveraging resources from immunology, cell biology and nutrition research communities, toxicoproteomics could make particular contributions in three inter-related areas to advance mechanistic insights and biomarker development: the plasma proteome and circulating microparticles, the adductome and idiosyncratic toxicity.

The p53 transcriptional target gene wnt7b contributes to NGF-inducible neurite outgrowth in neuronal PC12 cells

Differentiation of PC12 cells by nerve growth factor (NGF) is characterized by changes in signal transduction pathways leading to growth arrest and neurite extension. The transcription factor p53, involved in regulating cell cycle and apoptosis, is also activated during PC12 differentiation and contributes to each of these processes but the mechanisms are incompletely understood. NGF signaling stabilizes p53 protein expression, which enables its transcriptional regulation of target genes, including the newly identified target, wnt7b, a member of the wnt family of secreted morphogens. We tested the hypothesis that wnt7b expression is a factor in NGF-dependent neurite outgrowth of differentiating PC12 cells. Wnt7b transcript and protein levels are increased following NGF treatment in a p53-dependent manner, as demonstrated by a reduction in wnt7b protein levels following stable shRNA-mediated silencing of p53. In addition, overexpressed human tp53 was capable of inducing marked wnt7b expression in neuronal PC12 cells but tp53 overexpression did not elevate wnt7b levels in several tested human tumor cell lines. Ectopic wnt7b overexpression was sufficient to rescue neurite outgrowth in NGF-treated p53-silenced PC12 cells, which could be blocked by c-Jun N-terminal kinase (JNK) inhibition with SP600125 and did not involve beta-catenin nuclear translocation. Addition of sFRP1 to differentiation medium inhibited wnt7b-dependent phosphorylation of JNK, demonstrating that wnt7b is secreted and signals through a JNK-dependent mechanism in PC12 cells. We further identify an NGF-inducible subset of wnt receptors that likely supports wnt7b-mediated neurite extension in PC12 cells. In conclusion, wnt7b is a novel p53-regulated neuritogenic factor in PC12 cells that in conjunction with NGF-regulated Fzd expression is involved in p53-dependent neurite outgrowth through noncanonical JNK signaling.

CEBS--Chemical Effects in Biological Systems: a public data repository integrating study design and toxicity data with microarray and proteomics data

CEBS (Chemical Effects in Biological Systems) is an integrated public repository for toxicogenomics data, including the study design and timeline, clinical chemistry and histopathology findings and microarray and proteomics data. CEBS contains data derived from studies of chemicals and of genetic alterations, and is compatible with clinical and environmental studies. CEBS is designed to permit the user to query the data using the study conditions, the subject responses and then, having identified an appropriate set of subjects, to move to the microarray module of CEBS to carry out gene signature and pathway analysis. Scope of CEBS: CEBS currently holds 22 studies of rats, four studies of mice and one study of Caenorhabditis elegans. CEBS can also accommodate data from studies of human subjects. Toxicogenomics studies currently in CEBS comprise over 4000 microarray hybridizations, and 75 2D gel images annotated with protein identification performed by MALDI and MS/MS. CEBS contains raw microarray data collected in accordance with MIAME guidelines and provides tools for data selection, pre-processing and analysis resulting in annotated lists of genes of interest. Additionally, clinical chemistry and histopathology findings from over 1500 animals are included in CEBS. CEBS/BID: The BID (Biomedical Investigation Database) is another component of the CEBS system. BID is a relational database used to load and curate study data prior to export to CEBS, in addition to capturing and displaying novel data types such as PCR data, or additional fields of interest, including those defined by the HESI Toxicogenomics Committee (in preparation). BID has been shared with Health Canada and the US Environmental Protection Agency. CEBS is available at http://cebs.niehs.nih.gov. BID can be accessed via the user interface from https://dir-apps.niehs.nih.gov/arc/. Requests for a copy of BID and for depositing data into CEBS or BID are available at http://www.niehs.nih.gov/cebs-df/.

Nerve growth factor potentiates p53 DNA binding but inhibits nitric oxide-induced apoptosis in neuronal PC12 cells

NGF is recognized for its role in neuronal differentiation and maintenance. Differentiation of PC12 cells by NGF involves p53, a transcription factor that controls growth arrest and apoptosis. We investigated NGF influence over p53 activity during NO-induced apoptosis by sodium nitroprusside in differentiated and mitotic PC12 cells. NGF-differentiation produced increased p53 levels, nuclear localization and sequence-specific DNA binding. Apoptosis in mitotic cells also produced these events but the accompanying activation of caspases 1-10 and mitochondrial depolarization were inhibited during NGF differentiation and could be reversed in p53-silenced cells. Transcriptional regulation of PUMA and survivin expression were not inhibited by NGF, although NO-induced mitochondrial depolarization was dependent upon de novo gene transcription and only occurred in mitotic cells. We conclude that NGF mediates prosurvival signaling by increasing factors such as Bcl-2 and p21(Waf1/Cip1) without altering p53 transcriptional activity to inhibit mitochondrial depolarization, caspase activation and apoptosis.