2,3,7,8-Tetrachlorodibenzo-p-dioxin modifies alternative splicing in mouse liver

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.

Integrating gene expression and splicing dynamics across dose-response oxidative modulators

Toxicological risk assessment increasingly utilizes transcriptomics to derive point of departure (POD) and modes of action (MOA) for chemicals. One essential biological process that allows a single gene to generate several different RNA isoforms is called alternative splicing. To comprehensively assess the role of splicing dysregulation in toxicological evaluation and elucidate its potential as a complementary endpoint, we performed RNA-seq on A549 cells treated with five oxidative stress modulators across a wide dose range. Differential gene expression (DGE) showed limited pathway enrichment except at high concentrations. However, alternative splicing analysis revealed variable intron retention events affecting diverse pathways for all chemicals in the absence of significant expression changes. For instance, diazinon elicited negligible gene expression changes but progressive increase in the number of intron retention events, suggesting splicing alterations precede expression responses. Benchmark dose modeling of intron retention data highlighted relevant pathways overlooked by expression analysis. Systematic integration of splicing datasets should be a useful addition to the toxicogenomic toolkit. Combining both modalities paint a more complete picture of transcriptomic dose-responses. Overall, evaluating intron retention dynamics afforded by toxicogenomics may provide biomarkers that can enhance chemical risk assessment and regulatory decision making. This work highlights splicing-aware toxicogenomics as a possible additional tool for examining cellular responses.

Environmentally relevant concentrations of 2,3,7,8-TCDD induced inhibition of multicellular alternative splicing and transcriptional dysregulation

Alternative splicing (AS), found in approximately 95 % of human genes, significantly amplifies protein diversity and is implicated in disease pathogenesis when dysregulated. However, the precise involvement of AS in the toxic mechanisms induced by TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) remains incompletely elucidated. This study conducted a thorough global AS analysis in six human cell lines following TCDD exposure. Our findings revealed that environmentally relevant concentration (0.1 nM) of TCDD significantly suppressed AS events in all cell types, notably inhibiting diverse splicing events and reducing transcript diversity, potentially attributed to modifications in the splicing patterns of the inhibitory factor family, particularly hnRNP. And we identified 151 genes with substantial AS alterations shared among these cell types, particularly enriched in immune and metabolic pathways. Moreover, TCDD induced cell-specific changes in splicing patterns and transcript levels, with increased sensitivity notably in THP-1 monocyte, potentially linked to aberrant expression of pivotal genes within the spliceosome pathway (DDX5, EFTUD2, PUF60, RBM25, SRSF1, and CRNKL1). This study extends our understanding of disrupted alternative splicing and its relation to the multisystem toxicity of TCDD. It sheds light on how environmental toxins affect post-transcriptional regulatory processes, offering a fresh perspective for toxicology and disease etiology investigations.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD)-induced suppression of immunity in THP-1-derived macrophages and the possible mechanisms

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a well-known immunotoxic environmental pollutant. However, most immunotoxicology studies of TCDD were based on the animal models and the inner mechanisms have just focused on a few genes/proteins. In this study, the immune functions of THP-1-derived macrophages was measured with in-vitro bioassays after 24-h exposure of TCDD including environmentally relevant concentrations. RNA-seq and Weighted Gene Co-expression Network Analysis were used to characterize the immunotoxicity molecular mechanisms. Our study is the first report on the TCDD-induced effects of cell adhesion, morphology, and multiple cytokines/chemokines production on THP-1 macrophages. After TCDD treatment, we observed an inhibited cell adherence, probably attributed to the suppressed mRNA levels of adhesion molecules ICAM-1, VCAM-1 and CD11b, and a decrease in cell pseudopodia and expression of F-actin. The inflammatory cytokines TNF-α, IL-10 and other 8 cytokines/chemokines regulating granulocytes/T cells and angiogenesis were disrupted by TCDD. Alternative splicing event was found to be a sensitive target for TCDD. Using WGCNA, we identified 10 hub genes (TNF, SRC, FGF2, PTGS2, CDH2, GNG11, BDNF, WNT5A, CXCR5 and RUNX2) highly relevant to these observed phenotypes, suggesting AhR less important in the effects TCDD have on THP-1 macrophages than in other cells. Our findings broaden the understanding of TCDD immunotoxicity on macrophages and provide new potential targets for clarifying the molecular mechanisms.

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.