PCB126 Exposure Revealed Alterations in m6A RNA Modifications in Transcripts Associated With AHR Activation

Extractable organic matter from PM2.5 inhibits cardiomyocyte differentiation via AHR-mediated m6A RNA methylation

An ever-increasing body of research has established a link between maternal PM2.5 exposure and congenital heart diseases in the offspring, but the underlying mechanisms remain elusive. We recently reported that activation of the aryl hydrocarbon receptor (AHR) by PM2.5 causes aberrant m6A RNA methylation, leading to cardiac malformations in zebrafish embryos. We hypothesized that PM2.5 can disrupt heart development by inducing m6A methylation changes through AHR in mammals. In this study, we observed that extractable organic matters (EOM) from PM2.5 significantly impaired cardiomyocyte differentiation in embryonic rat cardiomyoblasts H9c2. Importantly, EOM exposure reduced global m6A methylation levels, which was reversed by AHR inhibition. Moreover, AHR, activated by EOM directly promoted the transcription of the demethylase, FTO, leading to global m6A hypomethylation. Specifically, AHR-induced FTO overexpression decreased the m6A methylation levels of Nox4 mRNA, resulting in NOX4 overexpression and subsequent oxidative stress in EOM samples. We then demonstrated that oxidative stress contributes to the inhibition of cardiomyocyte differentiation by EOM through suppression of Wnt/β-catenin signaling. In summary, our findings indicate that AHR activation by PM2.5 directly enhances the expression of the demethylase, FTO, which increases NOX4 expression by reducing its m6A methylation. The oxidative stress caused by NOX4 overexpression inhibits Wnt/β-catenin signaling, thereby compromising cardiomyocyte differentiation.

Global transcriptome modulation by xenobiotics: the role of alternative splicing in adaptive responses to chemical exposures

BACKGROUND

Xenobiotic exposures can extensively influence the expression and alternative splicing of drug-metabolizing enzymes, including cytochromes P450 (CYPs), though their transcriptome-wide impact on splicing remains underexplored. This study used a well-characterized splicing event in the Cyp2b2 gene to validate a sandwich-cultured primary rat hepatocyte model for studying global splicing in vitro. Using endpoint PCR, RNA sequencing, and bioinformatics tools (rSeqDiff, rMATs, IGV), we analyzed differential gene expression and splicing in CYP and nuclear receptor genes, as well as the entire transcriptome, to understand how xenobiotic exposures shape alternative splicing and activate xenosensors.

METHODS

Primary rat hepatocytes in sandwich culture were exposed to two methylenedioxybenzene (MDB) congeners and carbamazepine, with gene expression and splicing assessed. A 3D-clustergram integrating KEGG pathway analysis with differential gene expression provided distinct splicing landscapes for each xenobiotic, showing that splicing diversity does not always align with gene expression changes.

RESULTS

Endpoint PCR revealed a Cyp2b2v to wild-type Cyp2b2 splicing ratio near 1:1 (100%) under most conditions, while RNA-seq showed a stable baseline closer to 40%. C6-MDB reduced this ratio to ~ 50% by PCR and ~ 39% by RNA-seq, showing slight method-dependent variations yet consistent trends. In contrast, exon 6 skipping in Cyp1a1 occurred only with MDB exposure, implicating AHR activation. Xenobiotic treatments also induced alternative splicing in defensome and stress-responsive genes, including the phase II enzyme Gstm3, Albumin, Orm1, and Fxyd1, highlighting their roles in xenobiotic response modulation. Significant splicing changes in factors such as SRSF1, SRSF7, and METTL3 suggest a coordinated feedback loop involving epitranscriptomic modulation and cross-talk within SR protein networks, refining splice site selection, transcript stability, and cellular fate.

CONCLUSIONS

This study demonstrates how xenobiotic structural features influence gene expression and splicing, revealing splicing patterns that expand our understanding of transcriptome diversity and function. By identifying regulatory mechanisms, including AHR activation, epitranscriptomic modulation, and crosstalk within SR protein networks, that shape adaptive responses to xenobiotic stress, this work offers insights into the splicing and transcriptional networks that maintain cellular homeostasis. These findings provide predictive biomarkers for toxic exposures and highlight the potential of splicing profiles as diagnostic tools for assessing the health impacts of chemical exposure.

Targeting aryl hydrocarbon receptor to prevent cancer in barrier organs

The skin, lung, and gut are important barrier organs that control how the body reacts to environmental stressors such as ultraviolet (UV) radiation, air pollutants, dietary components, and microorganisms. The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that plays an important role in maintaining homeostasis of barrier organs. AhR was initially discovered as a receptor for environmental chemical carcinogens such as polycyclic aromatic hydrocarbons (PAHs). Activation of AhR pathways by PAHs leads to increased DNA damage and mutations which ultimately lead to carcinogenesis. Ongoing evidence reveals an ever-expanding role of AhR. Recently, AhR has been linked to immune systems by the interaction with the development of natural killer (NK) cells, regulatory T (Treg) cells, and T helper 17 (Th17) cells, as well as the production of immunosuppressive cytokines. However, the role of AhR in carcinogenesis is not as straightforward as we initially thought. Although AhR activation has been shown to promote carcinogenesis in some studies, others suggest that it may act as a tumor suppressor. In this review, we aim to explore the role of AhR in the development of cancer that originates from barrier organs. We also examined the preclinical efficacy data of AhR agonists and antagonists on carcinogenesis to determine whether AhR modulation can be a viable option for cancer chemoprevention.

Discovering Consensus Regions for Interpretable Identification of RNA N6-Methyladenosine Modification Sites via Graph Contrastive Clustering

As a pivotal post-transcriptional modification of RNA, N6-methyladenosine (m6A) has a substantial influence on gene expression modulation and cellular fate determination. Although a variety of computational models have been developed to accurately identify potential m6A modification sites, few of them are capable of interpreting the identification process with insights gained from consensus knowledge. To overcome this problem, we propose a deep learning model, namely M6A-DCR, by discovering consensus regions for interpretable identification of m6A modification sites. In particular, M6A-DCR first constructs an instance graph for each RNA sequence by integrating specific positions and types of nucleotides. The discovery of consensus regions is then formulated as a graph clustering problem in light of aggregating all instance graphs. After that, M6A-DCR adopts a motif-aware graph reconstruction optimization process to learn high-quality embeddings of input RNA sequences, thus achieving the identification of m6A modification sites in an end-to-end manner. Experimental results demonstrate the superior performance of M6A-DCR by comparing it with several state-of-the-art identification models. The consideration of consensus regions empowers our model to make interpretable predictions at the motif level. The analysis of cross validation through different species and tissues further verifies the consistency between the identification results of M6A-DCR and the evolutionary relationships among species.

AHR-mediated m6A RNA methylation contributes to PM2.5-induced cardiac malformations in zebrafish larvae

A growing body of evidence indicates that ambient fine particle matter (PM_2.5) exposure inhibits heart development, but the underlying mechanisms remain elusive. We hypothesized that m⁶A RNA methylation plays an important role in the cardiac developmental toxicity of PM_2.5. In this study, we demonstrated that extractable organic matter (EOM) from PM2.5 significantly decreased global m⁶A RNA methylation levels in the heart of zebrafish larvae, which were restored by the methyl donor, betaine. Betaine also attenuated EOM-induced ROS overgeneration, mitochondrial damage, apoptosis and heart defects. Furthermore, we found that the aryl hydrocarbon receptor (AHR), which was activated by EOM_, directly repressed the transcription of methyltransferases mettl14 and mettl3. EOM also induced genome-wide m⁶A RNA methylation changes, which led us to focus more on the aberrant m⁶A methylation changes that were subsequently alleviated by the AHR inhibitor, CH223191. In addition, we found that the expression levels of traf4a and bbc3, two apoptosis related genes, were upregulated by EOM but restored to control levels by the forced expression of mettl14. Moreover, knockdown of either traf4a or bbc3 attenuated EOM-induced ROS overproduction and apoptosis. In conclusion, our results indicate that PM_2.5 induces m⁶A RNA methylation changes via AHR-mediated mettl14 downregulation, which upregulates traf4a and bbc3, leading to apoptosis and cardiac malformations.

Molecular mechanisms of environmental exposures and human disease

A substantial proportion of disease risk for common complex disorders is attributable to environmental exposures and pollutants. An appreciation of how environmental pollutants act on our cells to produce deleterious health effects has led to advances in our understanding of the molecular mechanisms underlying the pathogenesis of chronic diseases, including cancer and cardiovascular, neurodegenerative and respiratory diseases. Here, we discuss emerging research on the interplay of environmental pollutants with the human genome and epigenome. We review evidence showing the environmental impact on gene expression through epigenetic modifications, including DNA methylation, histone modification and non-coding RNAs. We also highlight recent studies that evaluate recently discovered molecular processes through which the environment can exert its effects, including extracellular vesicles, the epitranscriptome and the mitochondrial genome. Finally, we discuss current challenges when studying the exposome - the cumulative measure of environmental influences over the lifespan - and its integration into future environmental health research.

Role of RNA modifications in carcinogenesis and carcinogen damage response

The field of epitranscriptomics encompasses the study of post-transcriptional RNA modifications and their regulatory enzymes. Among the numerous RNA modifications, N6 -methyladenosine (m6 A) has been identified as the most common internal modification of messenger RNA (mRNA). Although m6 A modifications were first discovered in the 1970s, advances in technology have revived interest in this field, driving an abundance of research into the role of RNA modifications in various biological processes, including cancer. As analogs to epigenetic modifications, RNA modifications also play an important role in carcinogenesis by regulating gene expression post-transcriptionally. A growing body of evidence suggests that carcinogens can modulate RNA modifications to alter the expression of oncogenes or tumor suppressors during cellular transformation. Additionally, the expression and activity of the enzymes that regulate RNA modifications can be dysregulated and contribute to carcinogenesis, making these enzymes promising targets of drug discovery. Here we summarize the roles of RNA modifications during carcinogenesis induced by exposure to various environmental carcinogens, with a main focus on the roles of the most widely studied m6 A mRNA methylation.

Polychlorinated biphenyls alter hepatic m6A mRNA methylation in a mouse model of environmental liver disease

Exposure to polychlorinated biphenyls (PCBs) has been associated with liver injury in human cohorts and with nonalcoholic steatohepatitis (NASH) in mice fed a high fat diet (HFD). N (6)-methyladenosine (m6A) modification of mRNA regulates transcript fate, but the contribution of m6A modification on the regulation of transcripts in PCB-induced steatosis and fibrosis is unknown. This study tested the hypothesis that PCB and HFD exposure alters the levels of m6A modification in transcripts that play a role in NASH in vivo. Male C57Bl6/J mice were fed a HFD (12 wks) and administered a single oral dose of Aroclor1260, PCB126, or Aroclor1260 + PCB126. Genome-wide identification of m6A peaks was accomplished by m6A mRNA immunoprecipitation sequencing (m6A-RIP) and the mRNA transcriptome identified by RNA-seq. Exposure of HFD-fed mice to Aroclor1260 decreased the number of m6A peaks and m6A-containing genes relative to PCB vehicle control whereas PCB126 or the combination of Aroclor1260 + PCB126 increased m6A modification frequency. ∼41% of genes had one m6A peak and ∼49% had 2-4 m6A peaks. 117 m6A peaks were common in the four experimental groups. The Aroclor1260 + PCB126 exposure group showed the highest number (52) of m6A-peaks. qRT-PCR confirmed enrichment of m6A-containing fragments of the Apob transcript with PCB exposure. A1cf transcript abundance, m6A peak count, and protein abundance was increased with Aroclor1260 + PCB126 co-exposure. Irrespective of the PCB type, all PCB groups exhibited enriched pathways related to lipid/lipoprotein metabolism and inflammation through the m6A modification. Integrated analysis of m6A-RIP-seq and mRNA-seq identified 242 differentially expressed genes (DEGs) with increased or reduced number of m6A peaks. These data show that PCB exposure in HFD-fed mice alters the m6A landscape offering an additional layer of regulation of gene expression affecting a subset of gene responses in NASH.

The role of N6-methyladenosine methylation in environmental exposure-induced health damage

The health risks caused by environmental pollution have long been of substantial concern. With the development of epigenetics, a large number of studies have demonstrated that N6-methyladenosine (m6A) modification is involved in the regulation of various important life activities associated with various diseases. Recent studies have revealed that m6A plays a key role in health damage caused by environmental exposure by regulating post-transcriptional gene expression. Therefore, our study outlined the effects of environmental pollutant exposure on m6A methylation and its regulator levels. Moreover, we found that m6A methylation modifications were involved in the development of various health damages by regulating important life activities in vivo, such as reactive oxygen species imbalance, apoptosis, epithelial-mesenchymal transition (EMT), and inflammatory processes. More importantly, we delved into the regulatory mechanisms of m6A methylation dysregulation in environmental pollution-induced diseases. Finally, by examining the published literature, we found that methyltransferase-like protein 3 (METTL3) and fat mass- and obesity-associated protein (FTO) were potentially used as biomarkers of health damage induced by particulate matter exposure and heavy metal exposure, respectively. The current studies on regulators of METTL3 and FTO were more promising to bring new perspectives for the treatment of environmental health-related diseases.

Critical role of m6A modification in T-helper cell disorders

Diseases with T-helper cell subset imbalance involve multiple systems and organs. In addition to this, the pathogenesis of these diseases is always complex, and involves Th1, Th2, Th9, Th17, Th22, and Tfh cells. T-helper cell subset imbalance mediates immune responses to various pathogenic factors, by secreting specific cytokines. Although several studies have revealed the specific mechanisms of the occurrence and development of these diseases from different aspects, there is still a need for more comprehensive and in-depth studies that can compensate for the corresponding gaps in the diagnosis, targeted therapy, and prognosis of these diseases. N6-methyladenosine(m⁶A) modification is the most prevalent and abundant post-transcriptional modification in eukaryotic RNAs. In recent years, the critical role of m⁶A modification has been confirmed in multiple diseases with T-helper cell subset imbalance. m⁶A modification affects the immune cell development, inflammatory processes, biological behaviour of tumours, and immune response in these diseases. In this review, we focussed on how the enzymes involved in m⁶A modification, directly or indirectly, influence the pathogenesis and phenotype of various diseases with T-helper cell subset imbalance, and could therefore, serve as potential diagnostic markers and therapeutic targets for these diseases. In addition, this review also discusses the focus of future research in this area. Finally, we summarise the prospects of m⁶A modification in immunotherapy and chemotherapy.

OUP accepted manuscript

Exposure to environmental toxicants during preconception has been shown to affect offspring health and epigenetic mechanisms such as DNA methylation are hypothesized to be involved in adverse outcomes. However, studies addressing the effects of exposure to environmental toxicants during preconception on epigenetic changes in gametes are limited. The objective of this study is to determine the effect of preconceptional exposure to a dioxin-like polychlorinated biphenyl (3,3',4,4',5-pentachlorobiphenyl [PCB126]) on DNA methylation and gene expression in testis. Adult zebrafish were exposed to 3 and 10 nM PCB126 for 24 h and testis tissue was sampled at 7 days postexposure for histology, DNA methylation, and gene expression profiling. Reduced representation bisulfite sequencing revealed 37 and 92 differentially methylated regions (DMRs) in response to 3 and 10 nM PCB126 exposures, respectively. Among them, 19 DMRs were found to be common between both PCB126 treatment groups. Gene ontology (GO) analysis of DMRs revealed that enrichment of terms such as RNA processing, iron-sulfur cluster assembly, and gluconeogenesis. Gene expression profiling showed differential expression of 40 and 1621 genes in response to 3 and 10 nM PCB126 exposures, respectively. GO analysis of differentially expressed genes revealed enrichment of terms related to xenobiotic metabolism, oxidative stress, and immune function. There is no overlap in the GO terms or individual genes between DNA methylation and RNA sequencing results, but functionally many of the altered pathways have been shown to cause spermatogenic defects.

Combined exposure to polychlorinated biphenyls and high-fat diet modifies the global epitranscriptomic landscape in mouse liver

Exposure to a single dose of polychlorinated biphenyls (PCBs) and a 12-week high-fat diet (HFD) results in nonalcoholic steatohepatitis (NASH) in mice by altering intracellular signaling and inhibiting epidermal growth factor receptor signaling. Post-transcriptional chemical modification (PTM) of RNA regulates biological processes, but the contribution of epitranscriptomics to PCB-induced steatosis remains unknown. This study tested the hypothesis that PCB and HFD exposure alters the global RNA epitranscriptome in male mouse liver. C57BL/6J male mice were fed a HFD for 12 weeks and exposed to a single dose of Aroclor 1260 (20 mg/kg), PCB 126 (20 µg/kg), both Aroclor 1260 and PCB 126 or vehicle control after 2 weeks on HFD. Chemical RNA modifications were identified at the nucleoside level by liquid chromatography-mass spectrometry. From 22 PTM global RNA modifications, we identified 10 significant changes in RNA modifications in liver with HFD and PCB 126 exposure. Only two modifications were significantly different from HFD control liver in all three PCB exposure groups: 2'-O-methyladenosine (Am) and N(6)-methyladenosine (m6A). Exposure to HFD + PCB 126 + Aroclor 1260 increased the abundance of N(6), O(2)-dimethyladenosine (m6Am), which is associated with the largest number of transcript changes. Increased m6Am and pseudouridine were associated with increased protein expression of the writers of these modifications: Phosphorylated CTD Interacting Factor 1 (PCIF1) and Pseudouridine Synthase 10 (PUS10), respectively, in HFD + PCB 126- + Aroclor 1260-exposed mouse liver. Increased N1-methyladenosine (m1A) and m6A were associated with increased transcript levels of the readers of these modifications: YTH N6-Methyladenosine RNA Binding Protein 2 (YTHDF2), YTH Domain Containing 2 (YTHDC2), and reader FMRP Translational Regulator 1 (FMR1) transcript and protein abundance. The results demonstrate that PCB exposure alters the global epitranscriptome in a mouse model of NASH; however, the mechanism for these changes requires further investigation.

Identification of molecular toxicity pathways across early life-stages of zebrafish exposed to PCB126 using a whole transcriptomics approach

Although withdrawn from the market in the 1980s, polychlorinated biphenyls (PCBs) are still found ubiquitously in the aquatic environment and pose a serious risk to biota due to their teratogenic potential. In fish, early life-stages are often considered most sensitive with regard to their exposure to PCBs and other dioxin-like compounds. However, little is known about the molecular drivers of the frequently observed teratogenic effects. Therefore, the aims of our study were to: (1) characterize the baseline transcriptome profiles at different embryonic life-stages in zebrafish (Danio rerio); and (2) to identify the molecular response to PCB exposure and life-stage specific-effects of the chemical on associated processes. For both objectives, embryos were sampled at 12, 48, and 96 h post-fertilization (hpf) and subjected to Illumina sequence-by-synthesis and RNAseq analysis. Results revealed that with increasing age more genes and related pathways were upregulated both in terms of number and magnitude. Yet, other transcripts followed an opposite pattern with greater transcript abundance at the earlier time points. Additionally, embryos were exposed to PCB126, a potent agonist of the aryl hydrocarbon receptor (AHR). ClueGO network analysis revealed significant enrichment of genes associated with basic cell metabolism, communication, and homeostasis as well as eye development, muscle formation, and skeletal formation. We selected eight genes involved in the affected pathways for an in-depth characterization of their regulation throughout normal embryogenesis and after exposure to PCB126 by quantification of transcript abundances every 12 h until 118 hpf. Among these, fgf7 and c9 stood out because of their strong upregulation by PCB126 exposure at 48 and 96 hpf, respectively. Cyp2aa12 was upregulated from 84 hpf on. Fabp10ab, myhz1.1, col8a1a, sulf1, and opn1sw1 displayed specific regulation depending on the developmental stage. Overall, we demonstrate that (1) the developmental transcriptome of zebrafish is highly dynamic, and (2) dysregulation of gene expression by exposure to PCB126 was significant and in several cases not directly connected to AHR-signaling. Hence, this study improves the understanding of linkages between molecular events and apical outcomes that are of regulatory relevance.

PCB126 Exposure Revealed Alterations in m6A RNA Modifications in Transcripts Associated With AHR Activation

Chemical modifications of proteins, DNA, and RNA moieties play critical roles in regulating gene expression. Emerging evidence suggests the RNA modifications (epitranscriptomics) have substantive roles in basic biological processes. One of the most common modifications in mRNA and noncoding RNAs is N6-methyladenosine (m6A). In a subset of mRNAs, m6A sites are preferentially enriched near stop codons, in 3' UTRs, and within exons, suggesting an important role in the regulation of mRNA processing and function including alternative splicing and gene expression. Very little is known about the effect of environmental chemical exposure on m6A modifications. As many of the commonly occurring environmental contaminants alter gene expression profiles and have detrimental effects on physiological processes, it is important to understand the effects of exposure on this important layer of gene regulation. Hence, the objective of this study was to characterize the acute effects of developmental exposure to PCB126, an environmentally relevant dioxin-like PCB, on m6A methylation patterns. We exposed zebrafish embryos to PCB126 for 6 h starting from 72 h post fertilization and profiled m6A RNA using methylated RNA immunoprecipitation followed by sequencing (MeRIP-seq). Our analysis revealed 117 and 217 m6A peaks in the DMSO and PCB126 samples (false discovery rate 5%), respectively. The majority of the peaks were preferentially located around the 3' UTR and stop codons. Statistical analysis revealed 15 m6A marked transcripts to be differentially methylated by PCB126 exposure. These include transcripts that are known to be activated by AHR agonists (eg, ahrra, tiparp, nfe2l2b) as well as others that are important for normal development (vgf, cebpd, sned1). These results suggest that environmental chemicals such as dioxin-like PCBs could affect developmental gene expression patterns by altering m6A levels. Further studies are necessary to understand the functional consequences of exposure-associated alterations in m6A levels.

A Review of the Functional Roles of the Zebrafish Aryl Hydrocarbon Receptors

Over the last 2 decades, the zebrafish (Danio rerio) has emerged as a stellar model for unraveling molecular signaling events mediated by the aryl hydrocarbon receptor (AHR), an important ligand-activated receptor found in all eumetazoan animals. Zebrafish have 3 AHRs-AHR1a, AHR1b, and AHR2, and studies have demonstrated the diversity of both the endogenous and toxicological functions of the zebrafish AHRs. In this contemporary review, we first highlight the evolution of the zebrafish ahr genes, and the characteristics of the receptors including developmental and adult expression, their endogenous and inducible roles, and the predicted ligands from homology modeling studies. We then review the toxicity of a broad spectrum of AHR ligands across multiple life stages (early stage, and adult), discuss their transcriptomic and epigenetic mechanisms of action, and report on any known interactions between the AHRs and other signaling pathways. Through this article, we summarize the promising research that furthers our understanding of the complex AHR pathway through the extensive use of zebrafish as a model, coupled with a large array of molecular techniques. As much of the research has focused on the functions of AHR2 during development and the mechanism of TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) toxicity, we illustrate the need to address the considerable knowledge gap in our understanding of both the mechanistic roles of AHR1a and AHR1b, and the diverse modes of toxicity of the various AHR ligands.