Tissue- and strain-specific effects of a genotoxic carcinogen 1,3-butadiene on chromatin and transcription

Visualization of the landscape of the read alignment shape of ATAC-seq data using Hellinger distance metric

Assay for Transposase-Accessible Chromatin using high-throughput sequencing (ATAC-seq) is the popular technique using next-generation sequencing to measure chromatin accessibility and identify open chromatin regions. While read alignment shape information of next-generation sequencing data with intensity information has been used in various bioinformatics methods, few studies have focused on pure shape information alone. In this study, we investigated what types of ATAC-seq read alignment shapes are observed for the promoter region and whether the pure shape information was related or unrelated to other gene features. We introduced a novel concept and pipeline for handling the pure shape information of NGS data as probability distributions and quantifying their dissimilarities by information theory. Based on this concept, we demonstrate that the pure shape information of ATAC-seq data is correlated with chromatin openness and some gene characteristics. On the other hand, it is suggested that the pure information of ATAC-seq read alignment shape is unlikely to contain additional information to explain differences in RNA expression. Our study suggests that viewing the read alignment shape of NGS data as probability distributions enables us to capture the characteristics of the genome-wide landscape of such data in a non-parametric manner.

Short-term in vivo testing to discriminate genotoxic carcinogens from non-genotoxic carcinogens and non-carcinogens using next-generation RNA sequencing, DNA microarray, and qPCR

Next-generation RNA sequencing (RNA-Seq) has identified more differentially expressed protein-coding genes (DEGs) and provided a wider quantitative range of expression level changes than conventional DNA microarrays. JEMS·MMS·Toxicogenomics group studied DEGs with targeted RNA-Seq on freshly frozen rat liver tissues and on formalin-fixed paraffin-embedded (FFPE) rat liver tissues after 28 days of treatment with chemicals and quantitative real-time PCR (qPCR) on rat and mouse liver tissues after 4 to 48 h treatment with chemicals and analyzed by principal component analysis (PCA) as statics. Analysis of rat public DNA microarray data (Open TG-GATEs) was also performed. In total, 35 chemicals were analyzed [15 genotoxic hepatocarcinogens (GTHCs), 9 non-genotoxic hepatocarcinogens (NGTHCs), and 11 non-genotoxic non-hepatocarcinogens (NGTNHCs)]. As a result, 12 marker genes (Aen, Bax, Btg2, Ccnf, Ccng1, Cdkn1a, Gdf15, Lrp1, Mbd1, Phlda3, Plk2, and Tubb4b) were proposed to discriminate GTHCs from NGTHCs and NGTNHCs. U.S. Environmental Protection Agency studied DEGs induced by 4 known GTHCs in rat liver using DNA microarray and proposed 7 biomarker genes, Bax, Bcmp1, Btg2, Ccng1, Cdkn1a, Cgr19, and Mgmt for GTHCs. Studies involving the use of whole-transcriptome RNA-Seq upon exposure to chemical carcinogens in vivo have also been performed in rodent liver, kidney, lung, colon, and other organs, although discrimination of GTHCs from NGTHCs was not examined. Candidate genes published using RNA-Seq, qPCR, and DNA microarray will be useful for the future development of short-term in vivo studies of environmental carcinogens using RNA-Seq.

Epigenetic alterations induced by genotoxic occupational and environmental human chemical carcinogens: An update of a systematic literature review

Epigenetic alterations, such as changes in DNA methylation, histones/chromatin structure, nucleosome positioning, and expression of non-coding RNAs, are recognized among key characteristics of carcinogens; they may occur independently or concomitantly with genotoxic effects. While data on genotoxicity are collected through standardized guideline tests, data collected on epigenetic effects is far less uniform. In 2016, we conducted a systematic review of published studies of genotoxic carcinogens that reported epigenetic endpoints to better understand the evidence for epigenetic alterations of human carcinogens, and the potential association with genotoxic endpoints. Since then, the number of studies of epigenetic effects of chemicals has nearly doubled. This review stands as an update on epigenetic alterations induced by occupational and environmental human carcinogens that were previously and recently classified as Group 1 by the International Agency for Research on Cancer. We found that the evidence of epigenetic effects remains uneven across agents. Studies of DNA methylation are most abundant, while reports concerning effects on non-coding RNA have increased over the past 5 years. By contrast, mechanistic toxicology studies of histone modifications and chromatin state alterations remain few. We found that most publications of epigenetic effects of carcinogens were studies in exposed humans or human cells. Studies in rodents represent the second most common species used for epigenetic studies in toxicology, in vivo exposures being the most predominant. Future studies should incorporate dose- and time-dependent study designs and also investigate the persistence of effects following cessation of exposure, considering the dynamic nature of most epigenetic alterations.

Harnessing natural variation to identify cis regulators of sex-biased gene expression in a multi-strain mouse liver model

Sex differences in gene expression are widespread in the liver, where many autosomal factors act in tandem with growth hormone signaling to regulate individual variability of sex differences in liver metabolism and disease. Here, we compare hepatic transcriptomic and epigenetic profiles of mouse strains C57BL/6J and CAST/EiJ, representing two subspecies separated by 0.5-1 million years of evolution, to elucidate the actions of genetic factors regulating liver sex differences. We identify 144 protein coding genes and 78 lncRNAs showing strain-conserved sex bias; many have gene ontologies relevant to liver function, are more highly liver-specific and show greater sex bias, and are more proximally regulated than genes whose sex bias is strain-dependent. The strain-conserved genes include key growth hormone-dependent transcriptional regulators of liver sex bias; however, three other transcription factors, Trim24, Tox, and Zfp809, lose their sex-biased expression in CAST/EiJ mouse liver. To elucidate the observed strain specificities in expression, we characterized the strain-dependence of sex-biased chromatin opening and enhancer marks at cis regulatory elements (CREs) within expression quantitative trait loci (eQTL) regulating liver sex-biased genes. Strikingly, 208 of 286 eQTLs with strain-specific, sex-differential effects on expression were associated with a complete gain, loss, or reversal of the sex differences in expression between strains. Moreover, 166 of the 286 eQTLs were linked to the strain-dependent gain or loss of localized sex-biased CREs. Remarkably, a subset of these CREs apparently lacked strain-specific genetic variants yet showed coordinated, strain-dependent sex-biased epigenetic regulation. Thus, we directly link hundreds of strain-specific genetic variants to the high variability in CRE activity and expression of sex-biased genes and uncover underlying genetically-determined epigenetic states controlling liver sex bias in genetically diverse mouse populations.

A systematic approach to evaluate plausible modes of actions for mouse lung tumors in mice exposed to 4-methylimidozole

The National Toxicology Program (NTP) reported that chronic dietary exposure to 4-methylimidazole (4-MeI) increased the incidence of lung adenomas/carcinomas beyond the normally high spontaneous rate in B6C3F1 mice. To examine plausible modes of action (MoAs) for mouse lung tumors (MLTs) upon exposure to high levels of 4-MeI, and their relevance in assessing human risk, a systematic approach was used to identify and evaluate mechanistic data (in vitro and in vivo) in the primary and secondary literature, along with high-throughput screening assay data. Study quality, relevance, and activity of mechanistic data identified across the evidence-base were organized according to key characteristics of carcinogens (KCCs) to identify potential key events in known or novel MLT MoAs. Integration of these evidence streams provided confirmation that 4-MeI lacks genotoxic and cytotoxic activity with some evidence to support a lack of mitogenic activity. Further evaluation of contextual and chemical-specific characteristics of 4-MeI was consequently undertaken. Due to lack of genotoxicity, along with transcriptomic and histopathological lung changes up to 28 and 90 days of exposure, the collective evidence suggests MLTs observed following exposure to high levels of 4-MeI develop at a late stage in the mouse chronic bioassay, albeit the exact MoA remains unclear.

Genotoxicity of Nanomaterials: Advanced In Vitro Models and High Throughput Methods for Human Hazard Assessment-A Review

Changes in the genetic material can lead to serious human health defects, as mutations in somatic cells may cause cancer and can contribute to other chronic diseases. Genotoxic events can appear at both the DNA, chromosomal or (during mitosis) whole genome level. The study of mechanisms leading to genotoxicity is crucially important, as well as the detection of potentially genotoxic compounds. We consider the current state of the art and describe here the main endpoints applied in standard human in vitro models as well as new advanced 3D models that are closer to the in vivo situation. We performed a literature review of in vitro studies published from 2000–2020 (August) dedicated to the genotoxicity of nanomaterials (NMs) in new models. Methods suitable for detection of genotoxicity of NMs will be presented with a focus on advances in miniaturization, organ-on-a-chip and high throughput methods.

Comparison of differential accessibility analysis strategies for ATAC-seq data

ATAC-seq is widely used to measure chromatin accessibility and identify open chromatin regions (OCRs). OCRs usually indicate active regulatory elements in the genome and are directly associated with the gene regulatory network. The identification of differential accessibility regions (DARs) between different biological conditions is critical in determining the differential activity of regulatory elements. Differential analysis of ATAC-seq shares many similarities with differential expression analysis of RNA-seq data. However, the distribution of ATAC-seq signal intensity is different from that of RNA-seq data, and higher sensitivity is required for DARs identification. Many different tools can be used to perform differential analysis of ATAC-seq data, but a comprehensive comparison and benchmarking of these methods is still lacking. Here, we used simulated datasets to systematically measure the sensitivity and specificity of six different methods. We further discussed the statistical and signal density cut-offs in the differential analysis of ATAC-seq by applying them to real data. Batch effects are very common in high-throughput sequencing experiments. We illustrated that batch-effect correction can dramatically improve sensitivity in the differential analysis of ATAC-seq data. Finally, we developed a user-friendly package, BeCorrect, to perform batch effect correction and visualization of corrected ATAC-seq signals in a genome browser.

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.

Drug-induced chromatin accessibility changes associate with sensitivity to liver tumor promotion

Liver cancer susceptibility varies amongst humans and between experimental animal models because of multiple genetic and epigenetic factors. The molecular characterization of such susceptibilities has the potential to enhance cancer risk assessment of xenobiotic exposures and disease prevention strategies. Here, using DNase I hypersensitivity mapping coupled with transcriptomic profiling, we investigate perturbations in cis-acting gene regulatory elements associated with the early stages of phenobarbital (PB)-mediated liver tumor promotion in susceptible versus resistant mouse strains (B6C3F1 versus C57BL/6J). Integrated computational analyses of strain-selective changes in liver chromatin accessibility underlying PB response reveal differential epigenetic regulation of molecular pathways associated with PB-mediated tumor promotion, including Wnt/β-catenin signaling. Complementary transcription factor motif analyses reveal mouse strain-selective gene regulatory networks and a novel role for Stat, Smad, and Fox transcription factors in the early stages of PB-mediated tumor promotion. Mapping perturbations in cis-acting gene regulatory elements provides novel insights into the molecular basis for susceptibility to xenobiotic-induced rodent liver tumor promotion and has the potential to enhance mechanism-based cancer risk assessments of xenobiotic exposures.

Using Collaborative Cross Mouse Population to Fill Data Gaps in Risk Assessment: A Case Study of Population-Based Analysis of Toxicokinetics and Kidney Toxicodynamics of Tetrachloroethylene

BACKGROUND

Interindividual variability in susceptibility remains poorly characterized for environmental chemicals such as tetrachloroethylene (PERC). Development of population-based experimental models provide a potential approach to fill this critical need in human health risk assessment.

OBJECTIVES

In this study, we aimed to better characterize the contribution of glutathione (GSH) conjugation to kidney toxicity of PERC and the degree of associated interindividual toxicokinetic (TK) and toxicodynamic (TD) variability by using the Collaborative Cross (CC) mouse population.

METHODS

Male mice from 45 strains were intragastrically dosed with PERC ([Formula: see text]) or vehicle (5% Alkamuls EL-620 in saline), and time-course samples were collected for up to 24 h. Population variability in TK of S-(1,2,2-trichlorovinyl)GSH (TCVG), S-(1,2,2-trichlorovinyl)-L-cysteine (TCVC), and N-acetyl-S-(1,2,2-trichlorovinyl)-L-cysteine (NAcTCVC) was quantified in serum, liver, and kidney, and analyzed using a toxicokinetic model. Effects of PERC on kidney weight, fatty acid metabolism-associated genes [ Acot1 (Acyl-CoA thioesterase 1), Fabp1 (fatty acid-binding protein 1), and Ehhadh (enoyl-coenzyme A, hydratase/3-hydroxyacyl coenzyme A dehydrogenase)], and a marker of proximal tubular injury [KIM-1 (kidney injury molecule-1)/Hepatitis A virus cellular receptor 1 ( Havcr1)] were evaluated. Finally, quantitative data on interstrain variability in both formation of GSH conjugation metabolites of PERC and its kidney effects was used to calculate adjustment factors for the interindividual variability in both TK and TD.

RESULTS

Mice treated with PERC had significantly lower kidney weight, higher kidney-to-body weight (BW) ratio, and higher expression of fatty acid metabolism-associated genes ( Acot1, Fabp1, and Ehhadh) and a marker of proximal tubular injury (KIM-1/ Havcr1). Liver levels of TCVG were significantly correlated with KIM-1/ Havcr1 in kidney, consistent with kidney injury being associated with GSH conjugation. We found that the default uncertainty factor for human variability may be marginally adequate to protect 95%, but not more, of the population for kidney toxicity mediated by PERC.

DISCUSSION

Overall, this study demonstrates the utility of the CC mouse population in characterizing metabolism-toxicity interactions and quantifying interindividual variability. Further refinement of the characterization of interindividual variability can be accomplished by incorporating these data into in silico population models both for TK (such as a physiologically based pharmacokinetic model), as well as for toxicodynamic responses. https://doi.org/10.1289/EHP5105.

Population-Based Analysis of DNA Damage and Epigenetic Effects of 1,3-Butadiene in the Mouse

Metabolism of 1,3-butadiene, a known human and rodent carcinogen, results in formation of reactive epoxides, a key event in its carcinogenicity. Although mice exposed to 1,3-butadiene present DNA adducts in all tested tissues, carcinogenicity is limited to liver, lung, and lymphoid tissues. Previous studies demonstrated that strain- and tissue-specific epigenetic effects in response to 1,3-butadiene exposure may influence susceptibly to DNA damage and serve as a potential mechanism of tissue-specific carcinogenicity. This study aimed to investigate interindividual variability in the effects of 1,3-butadiene using a population-based mouse model. Male mice from 20 Collaborative Cross strains were exposed to 0 or 635 ppm 1,3-butadiene by inhalation (6 h/day, 5 days/week) for 2 weeks. We evaluated DNA damage and epigenetic effects in target (lung and liver) and nontarget (kidney) tissues of 1,3-butadiene-induced carcinogenesis. DNA damage was assessed by measuring N-7-(2,3,4-trihydroxybut-1-yl)-guanine (THB-Gua) adducts. To investigate global histone modification alterations, we evaluated the trimethylation and acetylation of histones H3 and H4 across tissues. Changes in global cytosine DNA methylation were evaluated from the levels of methylation of LINE-1 and SINE B1 retrotransposons. We quantified the degree of variation across strains, deriving a chemical-specific human variability factor to address population variability in carcinogenic risk, which is largely ignored in current cancer risk assessment practice. Quantitative trait locus mapping identified four candidate genes related to chromatin remodeling whose variation was associated with interstrain susceptibility. Overall, this study uses 1,3-butadiene to demonstrate how the Collaborative Cross mouse population can be used to identify the mechanisms for and quantify the degree of interindividual variability in tissue-specific effects that are relevant to chemically induced carcinogenesis.

Sex-specific differences in genotoxic and epigenetic effects of 1,3-butadiene among mouse tissues

Exposure to environmental chemicals has been shown to have an impact on the epigenome. One example is a known human carcinogen 1,3-butadiene which acts primarily by a genotoxic mechanism, but also disrupts the chromatin structure by altering patterns of cytosine DNA methylation and histone modifications. Sex-specific differences in 1,3-butadiene-induced genotoxicity and carcinogenicity are well established; however, it remains unknown whether 1,3-butadiene-associated epigenetic alterations are also sex dependent. Therefore, we tested the hypothesis that inhalational exposure to 1,3-butadiene will result in sex-specific epigenetic alterations. DNA damage and epigenetic effects of 1,3-butadiene were evaluated in liver, lung, and kidney tissues of male and female mice of two inbred strains (C57BL/6J and CAST/EiJ). Mice were exposed to 0 or 425 ppm of 1,3-butadiene by inhalation (6 h/day, 5 days/week) for 2 weeks. Strain- and tissue-specific differences in 1,3-butadiene-induced DNA adducts and crosslinks were detected in the liver, lung and kidney; however, significant sex-specific differences in DNA damage were observed in the lung of C57BL/6J mice only. In addition, we assessed expression of the DNA repair genes and observed a marked upregulation of Mgmt in the kidney in female C57BL/6J mice. Sex-specific epigenetic effects of 1,3-butadiene exposure were evident in alterations of cytosine DNA methylation and histone modifications in the liver and lung in both strains. Specifically, we observed a loss of cytosine DNA methylation in the liver and lung of male and female 1,3-butadiene-exposed C57BL/6J mice, whereas hypermethylation was found in the liver and lung in 1,3-butadiene-exposed female CAST/EiJ mice. Our findings suggest that strain- and sex-specific effects of 1,3-butadiene on the epigenome may contribute to the known differences in cancer susceptibility.

DNA damaging agents trigger the expression of the HML silent mating type locus in Saccharomyces cerevisiae

Many DNA damaging agents also react with RNA and protein, and could thus cause epigenetic as well as genotoxic changes. To investigate which DNA damaging agents alter epigenetic states, we studied the chemical-induced changes in expression of the yeast silent mating type locus HMLα, which can be triggered by inhibiting yeast Sir2. We observed that the alkylating agent methyl methane sulfonate (MMS) can result in HMLα expression, using a colony sector assay that results from expression of a HML-positioned cre gene. Using single-cell imaging we also observed that alkylating agents, including MMS and methyl-3-nitro-1-nitrosoguanidine (MNNG), as well as short-wave UV, also decreased HML silencing. We suggest that chemical-induced alterations in heterochromatin structure could confer transient phenotypic changes that affect the cellular responses to DNA damaging agents.