Variation in DNA-Damage Responses to an Inhalational Carcinogen (1,3-Butadiene) in Relation to Strain-Specific Differences in Chromatin Accessibility and Gene Transcription Profiles in C57BL/6J and CAST/EiJ Mice

Stem cell transcriptional profiles from mouse subspecies reveal cis-regulatory evolution at translation genes

A key goal of evolutionary genomics is to harness molecular data to draw inferences about selective forces that have acted on genomes. The field progresses in large part through the development of advanced molecular-evolution analysis methods. Here we explored the intersection between classical sequence-based tests for selection and an empirical expression-based approach, using stem cells from Mus musculus subspecies as a model. Using a test of directional, cis-regulatory evolution across genes in pathways, we discovered a unique program of induction of translation genes in stem cells of the Southeast Asian mouse M. m. castaneus relative to its sister taxa. We then mined population-genomic sequences to pursue underlying regulatory mechanisms for this expression divergence, finding robust evidence for alleles unique to M. m. castaneus at the upstream regions of the translation genes. We interpret our data under a model of changes in lineage-specific pressures across Mus musculus in stem cells with high translational capacity. Our findings underscore the rigor of integrating expression and sequence-based methods to generate hypotheses about evolutionary events from long ago.

Stem cell transcriptional profiles from mouse subspecies reveal cis -regulatory evolution at translation genes

A key goal of evolutionary genomics is to harness molecular data to draw inferences about selective forces that have acted on genomes. The field progresses in large part through the development of advanced molecular-evolution analysis methods. Here we explored the intersection between classical sequence-based tests for selection and an empirical expression-based approach, using stem cells from Mus musculus subspecies as a model. Using a test of directional, cis -regulatory evolution across genes in pathways, we discovered a unique program of induction of translation genes in stem cells of the Southeast Asian mouse M. m. castaneus relative to its sister taxa. We then mined population-genomic sequences to pursue underlying regulatory mechanisms for this expression divergence, finding robust evidence for alleles unique to M. m. castaneus at the upstream regions of the translation genes. We interpret our data under a model of changes in lineage-specific pressures across Mus musculus in stem cells with high translational capacity. Our findings underscore the rigor of integrating expression and sequence-based methods to generate hypotheses about evolutionary events from long ago.

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.

Integrative analysis reveals mouse strain-dependent responses to acute ozone exposure associated with airway macrophage transcriptional activity

Acute ozone (O3) exposure is associated with multiple adverse cardiorespiratory outcomes, the severity of which varies across individuals in human populations and inbred mouse strains. However, molecular determinants of response, including susceptibility biomarkers that distinguish who will develop severe injury and inflammation, are not well characterized. We and others have demonstrated that airway macrophages (AMs) are an important resident immune cell type that are functionally and transcriptionally responsive to O3 inhalation. Here, we sought to explore influences of strain, exposure, and strain-by-O3 exposure interactions on AM gene expression and identify transcriptional correlates of O3-induced inflammation and injury across six mouse strains, including five Collaborative Cross (CC) strains. We exposed adult mice of both sexes to filtered air (FA) or 2 ppm O3 for 3 h and measured inflammatory and injury parameters 21 h later. Mice exposed to O3 developed airway neutrophilia and lung injury with strain-dependent severity. In AMs, we identified a common core O3 transcriptional response signature across all strains, as well as a set of genes exhibiting strain-by-O3 exposure interactions. In particular, a prominent gene expression contrast emerged between a low- (CC017/Unc) and high-responding (CC003/Unc) strain, as reflected by cellular inflammation and injury. Further inspection indicated that differences in their baseline gene expression and chromatin accessibility profiles likely contribute to their divergent post-O3 exposure transcriptional responses. Together, these results suggest that aspects of O3-induced respiratory responses are mediated through altered AM transcriptional signatures and further confirm the importance of gene-environment interactions in mediating differential responsiveness to environmental agents.

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.

Intra- and Inter-Species Variability in Urinary N7-(1-Hydroxy-3-buten-2-yl)guanine Adducts Following Inhalation Exposure to 1,3-Butadiene

1,3-Butadiene is a known carcinogen primarily targeting lymphoid tissues, lung, and liver. Cytochrome P450 activates butadiene to epoxides which form covalent DNA adducts that are thought to be a key mechanistic event in cancer. Previous studies suggested that inter-species, -tissue, and -individual susceptibility to adverse health effects of butadiene exposure may be due to differences in metabolism and other mechanisms. In this study, we aimed to examine the extent of inter-individual and inter-species variability in the urinary N7-(1-hydroxy-3-buten-2-yl)guanine (EB-GII) DNA adduct, a well-known biomarker of exposure to butadiene. For a population variability study in mice, we used the collaborative cross model. Female and male mice from five strains were exposed to filtered air or butadiene (590 ppm, 6 h/day, 5 days/week for 2 weeks) by inhalation. Urine samples were collected, and the metabolic activation of butadiene by DNA-reactive species was quantified as urinary EB-GII adducts. We quantified the degree of EB-GII variation across mouse strains and sexes; then, we compared this variation with the data from rats (exposed to 62.5 or 200 ppm butadiene) and humans (0.004-2.2 ppm butadiene). We show that sex and strain are significant contributors to the variability in urinary EB-GII levels in mice. In addition, we find that the degree of variability in urinary EB-GII in collaborative cross mice, when expressed as an uncertainty factor for the inter-individual variability (UF_H), is relatively modest (≤threefold) possibly due to metabolic saturation. By contrast, the variability in urinary EB-GII (adjusted for exposure) observed in humans, while larger than the default value of 10-fold, is largely consistent with UF_H estimates for other chemicals based on human data for non-cancer endpoints. Overall, these data demonstrate that urinary EB-GII levels, particularly from human studies, may be useful for quantitative characterization of human variability in cancer risks to butadiene.

Transcriptomic analyses of gastrulation-stage mouse embryos with differential susceptibility to alcohol

Genetics are a known contributor to differences in alcohol sensitivity in humans with fetal alcohol spectrum disorders (FASDs) and in animal models. Our study profiled gene expression in gastrulation-stage embryos from two commonly used, genetically similar mouse substrains, C57BL/6J (6J) and C57BL/6NHsd (6N), that differ in alcohol sensitivity. First, we established normal gene expression patterns at three finely resolved time points during gastrulation and developed a web-based interactive tool. Baseline transcriptional differences across strains were associated with immune signaling. Second, we examined the gene networks impacted by alcohol in each strain. Alcohol caused a more pronounced transcriptional effect in the 6J versus 6N mice, matching the increased susceptibility of the 6J mice. The 6J strain exhibited dysregulation of pathways related to cell death, proliferation, morphogenic signaling and craniofacial defects, while the 6N strain showed enrichment of hypoxia and cellular metabolism pathways. These datasets provide insight into the changing transcriptional landscape across mouse gastrulation, establish a valuable resource that enables the discovery of candidate genes that may modify alcohol susceptibility that can be validated in humans, and identify novel pathogenic mechanisms of alcohol. This article has an associated First Person interview with the first author of the paper.

SARS-CoV-2 infection and smoking: What is the association? A brief review

•

The link between smoking and the expression of SARS-CoV-2 key entry genes is discussed.

•

Smoking-related cardiac and respiratory diseases are risk factors for COVID-19.

•

The impact of smoking on ACE-2 and TMPRSS2 receptors expression is controversial.

Susceptibility to severe illness from COVID-19 is anticipated to be associated with cigarette smoking as it aggravates the risk of cardiovascular and respiratory illness, including infections. This is particularly important with the advent of a new strain of coronaviruses, the severe acute respiratory syndrome coronavirus (SARS-CoV-2) that has led to the present pandemic, coronavirus disease 2019 (COVID-19). Although, the effects of smoking on COVID-19 are less described and controversial, we presume a link between smoking and COVID-19. Smoking has been shown to enhance the expression of the angiotensin-converting enzyme-2 (ACE-2) and transmembrane serine protease 2 (TMPRSS2) key entry genes utilized by SARS-CoV-2 to infect cells and induce a ‘cytokine storm’, which further increases the severity of COVID-19 clinical course. Nevertheless, the impact of smoking on ACE-2 and TMPRSS2 receptors expression remains paradoxical. Thus, further research is necessary to unravel the association between smoking and COVID-19 and to pursue the development of potential novel therapies that are able to constrain the morbidity and mortality provoked by this infectious disease. Herein we present a brief overview of the current knowledge on the correlation between smoking and the expression of SARS-CoV-2 key entry genes, clinical manifestations, and disease progression.

Cigarette Smoke Exposure and Inflammatory Signaling Increase the Expression of the SARS-CoV-2 Receptor ACE2 in the Respiratory Tract

The factors mediating fatal SARS-CoV-2 infections are poorly understood. Here, we show that cigarette smoke causes a dose-dependent upregulation of angiotensin converting enzyme 2 (ACE2), the SARS-CoV-2 receptor, in rodent and human lungs. Using single-cell sequencing data, we demonstrate that ACE2 is expressed in a subset of secretory cells in the respiratory tract. Chronic smoke exposure triggers the expansion of this cell population and a concomitant increase in ACE2 expression. In contrast, quitting smoking decreases the abundance of these secretory cells and reduces ACE2 levels. Finally, we demonstrate that ACE2 expression is responsive to inflammatory signaling and can be upregulated by viral infections or interferon treatment. Taken together, these results may partially explain why smokers are particularly susceptible to severe SARS-CoV-2 infections. Furthermore, our work identifies ACE2 as an interferon-stimulated gene in lung cells, suggesting that SARS-CoV-2 infections could create positive feedback loops that increase ACE2 levels and facilitate viral dissemination.

Assessment of the Mode of Action Underlying the Effects of GenX in Mouse Liver and Implications for Assessing Human Health Risks

GenX is an alternative to environmentally persistent long-chain perfluoroalkyl and polyfluoroalkyl substances. Mice exposed to GenX exhibit liver hypertrophy, elevated peroxisomal enzyme activity, and other apical endpoints consistent with peroxisome proliferators. To investigate the potential role of peroxisome proliferator-activated receptor alpha (PPARα) activation in mice, and other molecular signals potentially related to observed liver changes, RNA sequencing was conducted on paraffin-embedded liver sections from a 90-day subchronic toxicity study of GenX conducted in mice. Differentially expressed genes were identified for each treatment group, and gene set enrichment analysis was conducted using gene sets that represent biological processes and known canonical pathways. Peroxisome signaling and fatty acid metabolism were among the most significantly enriched gene sets in both sexes at 0.5 and 5 mg/kg GenX; no pathways were enriched at 0.1 mg/kg. Gene sets specific to the PPARα subtype were significantly enriched. These findings were phenotypically anchored to histopathological changes in the same tissue blocks: hypertrophy, mitoses, and apoptosis. In vitro PPARα transactivation assays indicated that GenX activates mouse PPARα. These results indicate that the liver changes observed in GenX-treated mice occur via a mode of action (MOA) involving PPARα, an important finding for human health risk assessment as this MOA has limited relevance to humans.

Comparison of Gene Expression Responses in the Small Intestine of Mice Following Exposure to 3 Carcinogens Using the S1500+ Gene Set Informs a Potential Common Adverse Outcome Pathway

Carcinogenesis of the small intestine is rare in humans and rodents. Oral exposure to hexavalent chromium (Cr(VI)) and the fungicides captan and folpet induce intestinal carcinogenesis in mice. Previously (Toxicol Pathol. 330:48-52), we showed that B6C3F1 mice exposed to carcinogenic concentrations of Cr(VI), captan, or folpet for 28 days exhibited similar histopathological responses including villus enterocyte cytotoxicity and regenerative crypt epithelial hyperplasia. Herein, we analyze transcriptomic responses from formalin-fixed, paraffin-embedded duodenal sections from the aforementioned study. TempO-Seq technology and the S1500+ gene set were used to analyze transcription responses. Transcriptional responses were similar between all 3 agents; gene-level comparison identified 126/546 (23%) differentially expressed genes altered in the same direction, with a total of 25 upregulated pathways. These changes were related to cellular metabolism, stress, inflammatory/immune cell response, and cell proliferation, including upregulation in hypoxia inducible factor 1 (HIF-1) and activator protein 1 (AP1) signaling pathways, which have also been shown to be related to intestinal injury and angiogenesis/carcinogenesis. The similar molecular-, cellular-, and tissue-level changes induced by these 3 carcinogens can be informative for the development of an adverse outcome pathway for intestinal cancer.

Chemically induced carcinogenesis in rodent models of aging: assessing organismal resilience to genotoxic stressors in geroscience research

There is significant overlap between the cellular and molecular mechanisms of aging and pathways contributing to carcinogenesis, including the role of genome maintenance pathways. In the field of geroscience analysis of novel genetic mouse models with either a shortened, or an extended, lifespan provides a unique opportunity to evaluate the synergistic roles of longevity assurance pathways in cancer resistance and regulation of lifespan and to develop novel targets for interventions that both delay aging and prevent carcinogenesis. There is a growing need for robust assays to assess the susceptibility of cancer in these models. The present review focuses on a well-characterized method frequently used in cancer research, which can be adapted to study resilience to genotoxic stress and susceptibility to genotoxic stress-induced carcinogenesis in geroscience research namely, chemical carcinogenesis induced by treatment with 7,12-dimethylbenz(a)anthracene (DMBA). Recent progress in understanding how longer-living mice may achieve resistance to chemical carcinogenesis and how these pathways are modulated by anti-aging interventions is reviewed. Strain-specific differences in sensitivity to DMBA-induced carcinogenesis are also explored and contrasted with mouse lifespan. The clinical relevance of inhibition of DMBA-induced carcinogenesis for the pathogenesis of mammary adenocarcinomas in older human subjects is discussed. Finally, the potential role of insulin-like growth factor-1 (IGF-1) in the regulation of pathways responsible for cellular resilience to DMBA-induced mutagenesis is discussed.

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.

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

Epigenetic effects of environmental chemicals are under intense investigation to fill existing knowledge gaps between environmental/occupational exposures and adverse health outcomes. Chromatin accessibility is one prominent mechanism of epigenetic control of transcription, and understanding of the chemical effects on both could inform the causal role of epigenetic alterations in disease mechanisms. In this study, we hypothesized that baseline variability in chromatin organization and transcription profiles among various tissues and mouse strains influence the outcome of exposure to the DNA damaging chemical 1,3-butadiene. To test this hypothesis, we evaluated DNA damage along with comprehensive quantification of RNA transcripts (RNA-seq), identification of accessible chromatin (ATAC-seq), and characterization of regions with histone modifications associated with active transcription (ChIP-seq for acetylation at histone 3 lysine 27, H3K27ac). We collected these data in the lung, liver, and kidney of mice from two genetically divergent strains, C57BL/6J and CAST/EiJ, that were exposed to clean air or to 1,3-butadiene (~600 ppm) for 2 weeks. We found that tissue effects dominate differences in both gene expression and chromatin states, followed by strain effects. At baseline, xenobiotic metabolism was consistently more active in CAST/EiJ, while immune system pathways were more active in C57BL/6J across tissues. Surprisingly, even though all three tissues in both strains harbored butadiene-induced DNA damage, little transcriptional effect of butadiene was observed in liver and kidney. Toxicologically relevant effects of butadiene in the lung were on the pathways of xenobiotic metabolism and inflammation. We also found that variability in chromatin accessibility across individuals (i.e., strains) only partially explains the variability in transcription. This study showed that variation in the basal states of epigenome and transcriptome may be useful indicators for individuals or tissues susceptible to genotoxic environmental chemicals.

Genetic and epigenetic determinants of inter-individual variability in responses to toxicants

It is well established that genetic variability has a major impact on susceptibility to common diseases, responses to drugs and toxicants, and influences disease-related outcomes. The appreciation that epigenetic marks also vary across the population is growing with more data becoming available from studies in humans and model organisms. In addition, the links between genetic variability, toxicity outcomes and epigenetics are being actively explored. Recent studies demonstrate that gene-by-environment interactions involve both chromatin states and transcriptional regulation, and that epigenetics provides important mechanistic clues to connect expression-related quantitative trait loci (QTL) and disease outcomes. However, studies of Gene×Environment×Epigenetics further extend the complexity of the experimental designs and create a challenge for selecting the most informative epigenetic readouts that can be feasibly performed to interrogate multiple individuals, exposures, tissue types and toxicity phenotypes. We propose that among the many possible epigenetic experimental methodologies, assessment of chromatin accessibility coupled with total RNA levels provides a cost-effective and comprehensive option to sufficiently characterize the complexity of epigenetic and regulatory activity in the context of understanding the inter-individual variability in responses to toxicants.