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.

Transcriptomic re-analyses of human hepatocyte spheroids treated with PFAS reveals chain length and dose-dependent modes of action.

To identify pathway perturbations and examine biological modes of action (MOAs) for various perfluoroalkyl substances, we re-analyzed published in vitro gene expression studies from human primary liver spheroids. With treatment times ranging from 10 to 14 days, shorter-chain PFAS (those with 6 or fewer fluorinated carbon atoms in the alkyl chain) showed enrichment for pathways of fatty acid metabolism and fatty acid beta-oxidation with upregulated genes. Longer-chain PFAS compounds, specifically PFOS (perfluorooctane sulfonate), PFDS (perfluorodecane sulfonate), and higher doses of PFOA (perfluorooctanoic acid), had enrichment for pathways involved in steroid metabolism, fatty acid metabolism, and biological oxidation for downregulated genes. Although PFNA (perfluorononanoic acid), PFDA (perfluorodecanoic acid), and PFUnDA (perfluoroundecanoic acid) were more toxic and could only be examined after a 1-day treatment, all three had enrichment patterns similar to those observed with PFOS. With PFOA there were dose-dependent changes in pathway enrichment, shifting from upregulation of fatty acid metabolism and downregulation of steroid metabolism to downregulation of both at higher doses. The response to PFHpS (perfluoroheptanesulfonic acid) was similar to the PFOA pattern at the lower treatment dose. Based on results of transcription factor binding sites analyses, we propose that downregulation of pathways of lipid metabolism by longer chain PFAS may be due to inhibitory interactions of PPARD on genes controlled by PPARA and PPARG. In conclusion, our transcriptomic analysis indicates that the biological MOAs of PFAS compounds differ according to chain length and dose, and that risk assessments for PFAS should consider these differences in biological MOAs when evaluating mixtures of these compounds.

Transcriptomic analysis of AHR wildtype and Knock-out rat livers supports TCDD's role in AHR/ARNT-mediated circadian disruption and hepatotoxicity.

Single, high doses of TCDD in rats are known to cause wasting, a progressive loss of 30 to 50% body weight and death within several weeks. To identify pathway perturbations at or near doses causing wasting, we examined differentially gene expression (DGE) and pathway enrichment in centrilobular (CL) and periportal (PP) regions of female rat livers following 6 dose levels of TCDD - 0, 3, 22, 100, 300, and 1000 ng/kg/day, 5 days/week for 4 weeks. At the higher doses, rats lost weight, had increased liver/body weight ratios and nearly complete cessation of liver cell proliferation, signs consistent with wasting. DGE curves were left shifted for the CL versus the PP regions. Canonical Phase I and Phase II genes were maximally increased at lower doses and remained elevated at all doses. At lower doses, ≤ 22 ng/kg/day in the CL and ≤ 100 ng/kg/day, upregulated genes showed transcription factor (TF) enrichment for AHR and ARNT. At the mid- and high-dose doses, there was a large number of downregulated genes and pathway enrichment for DEGs which showed downregulation of many cellular metabolism processes including those for steroids, fatty acid metabolism, pyruvate metabolism and citric acid cycle. There was significant TF enrichment of the hi-dose downregulated genes for RXR, ESR1, LXR, PPARalpha. At the highest dose, there was also pathway enrichment with upregulated genes for extracellular matrix organization, collagen formation, hemostasis and innate immune system. TCDD demonstrates most of its effects through binding the aryl hydrocarbon receptor (AHR) while the downregulation of metabolism genes at higher TCDD doses is known to be independent of AHR binding to DREs. Based on our results with DEG, we provide a hypothesis for wasting in which high doses of TCDD shift circadian processes away from the resting state, leading to greatly reduced synthesis of steroids and complex lipids needed for cell growth, and producing gene expression signals consistent with an epithelial-to-mesenchymal transition in hepatocytes.

Assessment of TGx-DDI genes for genotoxicity in a comprehensive panel of chemicals.

BACKGROUND: The TGx-DDI biomarker identifies transcripts specifically induced by primary DNA damage. Profiling similarity of TGx-DDI signatures can allow clustering compounds by genotoxic mechanism. This transcriptomics-based approach complements conventional toxicology testing by enhancing mechanistic resolution. METHODS: Unsupervised hierarchical clustering and t-distributed stochastic neighbor embedding (tSNE) were utilized to assess similarity of publicly-available per- and polyfluoroalkyl substances (PFAS) and ToxCast chemicals based on TGx-DDI modulation. TempO-seq transcriptomic data after highest chemical concentrations were analyzed. RESULTS: Clustering discriminated between genotoxic and non-genotoxic compounds while drawing similarity among chemicals with shared mechanisms. PFAS largely clustered distinctly from classical mutagens. However, dynamic range across PFAS types and durations indicated variable potential for DNA damage. tSNE visualization reinforced phenotypic groupings, with genotoxins clustering separately from non-DNA damaging agents. DISCUSSION: Unsupervised learning approaches applied to TGx-DDI profiles effectively categorizes chemical genotoxicity potential, aiding elucidation of biological response pathways. This transcriptomics-based strategy gives further insight into the role and effect of individual TGx-DDI biomarker genes and complements existing assays by enhancing mechanistic resolution. Overall, TGx-DDI biomarker profiling holds promise for predictive safety screening.

Mining toxicogenomic data for dose-responsive pathways: implications in advancing next-generation risk assessment.

Introduction: While targeted investigation of key toxicity pathways has been instrumental for biomarker discovery, unbiased and holistic analysis of transcriptomic data provides a complementary systems-level perspective. However, in a systematic context, this approach has yet to receive comprehensive and methodical implementation. Methods: Here, we took an integrated bioinformatic approach by re-analyzing publicly available MCF7 cell TempO-seq data for 44 ToxCast chemicals using an alternative pipeline to demonstrate the power of this approach. The original study has focused on analyzing the gene signature approach and comparing them to in vitro biological pathway altering concentrations determined from ToxCast HTS assays. Our workflow, in comparison, involves sequential differential expression, gene set enrichment, benchmark dose modeling, and identification of commonly perturbed pathways by network visualization. Results: Using this approach, we identified dose-responsive molecular changes, biological pathways, and points of departure in an untargeted manner. Critically, benchmark dose modeling based on pathways recapitulated points of departure for apical endpoints, while also revealing additional perturbed mechanisms missed by single endpoint analyses. Discussion: This systems-toxicology approach provides multifaceted insights into the complex effects of chemical exposures. Our work highlights the importance of unbiased data-driven techniques, alongside targeted methods, for comprehensively evaluating molecular initiating events, dose-response relationships, and toxicity pathways. Overall, integrating omics assays with robust bioinformatics holds promise for improving chemical risk assessment and advancing new approach methodologies (NAMs).

Increasing confidence in new approach methodologies for inhalation risk assessment with multiple end point assays using 5-day repeated exposure to 1,3-dichloropropene.

New Approach Methodologies (NAMs) are being widely used to reduce, refine, and replace, animal use in studying toxicology. For respiratory toxicology, this includes both in silico and in vitro alternatives to replace traditional in vivo inhalation studies. 1,3-Dichloropropene (1,3-DCP) is a volatile organic compound that is widely used in agriculture as a pre-planting fumigant. Short-term exposure of humans to 1,3-DCP can result in mucous membrane irritation, chest pain, headache, and dizziness. In our previous work, we exposed differentiated cells representing different parts of the respiratory epithelium to 1,3-DCP vapor, measured cytotoxicity, and did In Vitro to In Vivo Extrapolation (IVIVE). We have extended our previous study with 1,3-DCP vapors by conducting transcriptomics on acutely exposed nasal cultures and have implemented a separate 5-day repeated exposure with multiple endpoints to gain further molecular insight into our model. MucilAir™ Nasal cell culture models, representing the nasal epithelium, were exposed to six sub-cytotoxic concentrations of 1,3-DCP vapor at the air-liquid interface, and the nasal cultures were analyzed by different methodologies, including histology, transcriptomics, and glutathione (GSH) -depletion assays. We observed the dose-dependent effect of 1,3-DCP in terms of differential gene expression, change in cellular morphology from pseudostratified columnar epithelium to squamous epithelium, and depletion of GSH in MucilAir™ nasal cultures. The MucilAir™ nasal cultures were also exposed to 3 concentrations of 1,3-DCP using repeated exposure 4 h per day for 5 days and the histological analyses indicated changes in cellular morphology and a decrease in ciliated bodies and an increase in apoptotic bodies, with increasing concentrations of 1,3-DCP. Altogether, our results suggest that sub-cytotoxic exposures to 1,3-DCP lead to several molecular and cellular perturbations, providing significant insight into the mode-of-action (MoA) of 1,3-DCP using an innovative NAM model.

Prolonged FOS activity disrupts a global myogenic transcriptional program by altering 3D chromatin architecture in primary muscle progenitor cells.

BACKGROUND: The AP-1 transcription factor, FBJ osteosarcoma oncogene (FOS), is induced in adult muscle satellite cells (SCs) within hours following muscle damage and is required for effective stem cell activation and muscle repair. However, why FOS is rapidly downregulated before SCs enter cell cycle as progenitor cells (i.e., transiently expressed) remains unclear. Further, whether boosting FOS levels in the proliferating progeny of SCs can enhance their myogenic properties needs further evaluation. METHODS: We established an inducible, FOS expression system to evaluate the impact of persistent FOS activity in muscle progenitor cells ex vivo. We performed various assays to measure cellular proliferation and differentiation, as well as uncover changes in RNA levels and three-dimensional (3D) chromatin interactions. RESULTS: Persistent FOS activity in primary muscle progenitor cells severely antagonizes their ability to differentiate and form myotubes within the first 2 weeks in culture. RNA-seq analysis revealed that ectopic FOS activity in muscle progenitor cells suppressed a global pro-myogenic transcriptional program, while activating a stress-induced, mitogen-activated protein kinase (MAPK) transcriptional signature. Additionally, we observed various FOS-dependent, chromosomal re-organization events in A/B compartments, topologically associated domains (TADs), and genomic loops near FOS-regulated genes. CONCLUSIONS: Our results suggest that elevated FOS activity in recently activated muscle progenitor cells perturbs cellular differentiation by altering the 3D chromosome organization near critical pro-myogenic genes. This work highlights the crucial importance of tightly controlling FOS expression in the muscle lineage and suggests that in states of chronic stress or disease, persistent FOS activity in muscle precursor cells may disrupt the muscle-forming process.

Systematic mapping of nuclear domain-associated transcripts reveals speckles and lamina as hubs of functionally distinct retained introns.

The nucleus is highly compartmentalized through the formation of distinct classes of membraneless domains. However, the composition and function of many of these structures are not well understood. Using APEX2-mediated proximity labeling and RNA sequencing, we surveyed human transcripts associated with nuclear speckles, several additional domains, and the lamina. Remarkably, speckles and lamina are associated with distinct classes of retained introns enriched in genes that function in RNA processing, translation, and the cell cycle, among other processes. In contrast to the lamina-proximal introns, retained introns associated with speckles are relatively short, GC-rich, and enriched for functional sites of RNA-binding proteins that are concentrated in these domains. They are also highly differentially regulated across diverse cellular contexts, including the cell cycle. Thus, our study provides a resource of nuclear domain-associated transcripts and further reveals speckles and lamina as hubs of distinct populations of retained introns linked to gene regulation and cell cycle progression.

Differential contribution of steady-state RNA and active transcription in chromatin organization.

Nuclear RNA and the act of transcription have been implicated in nuclear organization. However, their global contribution to shaping fundamental features of higher-order chromatin organization such as topologically associated domains (TADs) and genomic compartments remains unclear. To investigate these questions, we perform genome-wide chromatin conformation capture (Hi-C) analysis in the presence and absence of RNase before and after crosslinking, or a transcriptional inhibitor. TAD boundaries are largely unaffected by RNase treatment, although a subtle disruption of compartmental interactions is observed. In contrast, transcriptional inhibition leads to weaker TAD boundary scores. Collectively, our findings demonstrate differences in the relative contribution of RNA and transcription to the formation of TAD boundaries detected by the widely used Hi-C methodology.

Interchromosomal interactions: A genomic love story of kissing chromosomes.

Nuclei require a precise three- and four-dimensional organization of DNA to establish cell-specific gene-expression programs. Underscoring the importance of DNA topology, alterations to the nuclear architecture can perturb gene expression and result in disease states. More recently, it has become clear that not only intrachromosomal interactions, but also interchromosomal interactions, a less studied feature of chromosomes, are required for proper physiological gene-expression programs. Here, we review recent studies with emerging insights into where and why cross-chromosomal communication is relevant. Specifically, we discuss how long noncoding RNAs (lncRNAs) and three-dimensional gene positioning are involved in genome organization and how low-throughput (live-cell imaging) and high-throughput (Hi-C and SPRITE) techniques contribute to understand the fundamental properties of interchromosomal interactions.

Enhancers in the Peril lincRNA locus regulate distant but not local genes.

BACKGROUND: Recently, it has become clear that some promoters function as long-range regulators of gene expression. However, direct and quantitative assessment of enhancer activity at long intergenic noncoding RNA (lincRNA) or mRNA gene bodies has not been performed. To unbiasedly assess the enhancer capacity across lincRNA and mRNA loci, we performed a massively parallel reporter assay (MPRA) on six lincRNA loci and their closest protein-coding neighbors. RESULTS: For both gene classes, we find significantly more MPRA activity in promoter regions than in gene bodies. However, three lincRNA loci, Lincp21, LincEnc1, and Peril, and one mRNA locus, Morc2a, display significant enhancer activity within their gene bodies. We hypothesize that such peaks may mark long-range enhancers, and test this in vivo using RNA sequencing from a knockout mouse model and high-throughput chromosome conformation capture (Hi-C). We find that ablation of a high-activity MPRA peak in the Peril gene body leads to consistent dysregulation of Mccc1 and Exosc9 in the neighboring topologically associated domain (TAD). This occurs irrespective of Peril lincRNA expression, demonstrating this regulation is DNA-dependent. Hi-C confirms long-range contacts with the neighboring TAD, and these interactions are altered upon Peril knockout. Surprisingly, we do not observe consistent regulation of genes within the local TAD. Together, these data suggest a long-range enhancer-like function for the Peril gene body. CONCLUSIONS: A multi-faceted approach combining high-throughput enhancer discovery with genetic models can connect enhancers to their gene targets and provides evidence of inter-TAD gene regulation.

Reorganization of inter-chromosomal interactions in the 2q37-deletion syndrome.

Chromosomes occupy distinct interphase territories in the three-dimensional nucleus. However, how these chromosome territories are arranged relative to one another is poorly understood. Here, we investigated the inter-chromosomal interactions between chromosomes 2q, 12, and 17 in human mesenchymal stem cells (MSCs) and MSC-derived cell types by DNA-FISH We compared our findings in normal karyotypes with a three-generation family harboring a 2q37-deletion syndrome, featuring a heterozygous partial deletion of histone deacetylase 4 (HDAC4) on chr2q37. In normal karyotypes, we detected stable, recurring arrangements and interactions between the three chromosomal territories with a tissue-specific interaction bias at certain loci. These inter-chromosomal interactions were confirmed by Hi-C. Interestingly, the disease-related HDAC4 deletion resulted in displaced inter-chromosomal arrangements and altered interactions between the deletion-affected chromosome 2 and chromosome 12 and/or 17 in 2q37-deletion syndrome patients. Our findings provide evidence for a direct link between a structural chromosomal aberration and altered interphase architecture that results in a nuclear configuration, supporting a possible molecular pathogenesis.

A TAD boundary is preserved upon deletion of the CTCF-rich Firre locus.

The binding of the transcriptional regulator CTCF to the genome has been implicated in the formation of topologically associated domains (TADs). However, the general mechanisms of folding the genome into TADs are not fully understood. Here we test the effects of deleting a CTCF-rich locus on TAD boundary formation. Using genome-wide chromosome conformation capture (Hi-C), we focus on one TAD boundary on chromosome X harboring ~ 15 CTCF binding sites and located at the long non-coding RNA (lncRNA) locus Firre. Specifically, this TAD boundary is invariant across evolution, tissues, and temporal dynamics of X-chromosome inactivation. We demonstrate that neither the deletion of this locus nor the ectopic insertion of Firre cDNA or its ectopic expression are sufficient to alter TADs in a sex-specific or allele-specific manner. In contrast, Firre's deletion disrupts the chromatin super-loop formation of the inactive X-chromosome. Collectively, our findings suggest that apart from CTCF binding, additional mechanisms may play roles in establishing TAD boundary formation.

Inter-chromosomal Contact Properties in Live-Cell Imaging and in Hi-C.

Imaging (fluorescence in situ hybridization [FISH]) and genome-wide chromosome conformation capture (Hi-C) are two major approaches to the study of higher-order genome organization in the nucleus. Intra-chromosomal and inter-chromosomal interactions (referred to as non-homologous chromosomal contacts [NHCCs]) have been observed by several FISH-based studies, but locus-specific NHCCs have not been detected by Hi-C. Due to crosslinking, neither of these approaches assesses spatiotemporal properties. Toward resolving the discrepancies between imaging and Hi-C, we sought to understand the spatiotemporal properties of NHCCs in living cells by CRISPR/Cas9 live-cell imaging (CLING). In mammalian cells, we find that NHCCs are stable and occur as frequently as intra-chromosomal interactions, but NHCCs occur at farther spatial distance that could explain their lack of detection in Hi-C. By revealing the spatiotemporal properties in living cells, our study provides fundamental insights into the biology of NHCCs.

Spatiotemporal allele organization by allele-specific CRISPR live-cell imaging (SNP-CLING).

Imaging and chromatin capture techniques have provided important insights into our understanding of nuclear organization. A limitation of these techniques is the inability to resolve allele-specific spatiotemporal properties of genomic loci in living cells. Here, we describe an allele-specific CRISPR live-cell DNA imaging technique (SNP-CLING) to provide the first comprehensive insights into allelic positioning across space and time in mouse embryonic stem cells and fibroblasts. With 3D imaging, we studied alleles on different chromosomes in relation to one another and relative to nuclear substructures such as the nucleolus. We find that alleles maintain similar positions relative to each other and the nucleolus; however, loci occupy unique positions. To monitor spatiotemporal dynamics by SNP-CLING, we performed 4D imaging and determined that alleles are either stably positioned or fluctuating during cell state transitions, such as apoptosis. SNP-CLING is a universally applicable technique that enables the dissection of allele-specific spatiotemporal genome organization in live cells.

The connection between BRG1, CTCF and topoisomerases at TAD boundaries.

The eukaryotic genome is partitioned into topologically associating domains (TADs). Despite recent advances characterizing TADs and TAD boundaries, the organization of these structures is an important dimension of genome architecture and function that is not well understood. Recently, we demonstrated that knockdown of BRG1, an ATPase driving the chromatin remodeling activity of mammalian SWI/SNF enzymes, globally alters long-range genomic interactions and results in a reduction of TAD boundary strength. We provided evidence suggesting that this effect may be due to BRG1 affecting nucleosome occupancy around CTCF sites present at TAD boundaries. In this review, we elaborate on our findings and speculate that BRG1 may contribute to the regulation of the structural and functional properties of chromatin at TAD boundaries by affecting the function or the recruitment of CTCF and DNA topoisomerase complexes.

Identifying Nuclear Matrix-Attached DNA Across the Genome.

Experimental approaches to define the relationship between gene expression and nuclear matrix attachment regions (MARs) have given contrasting and method-specific results. We have developed a next generation sequencing strategy to identify MARs across the human genome (MAR-Seq). The method is based on crosslinking chromatin to its nuclear matrix attachment sites to minimize changes during biochemical processing. We used this method to compare nuclear matrix organization in MCF-10A mammary epithelial-like cells and MDA-MB-231 breast cancer cells and evaluated the results in the context of global gene expression (array analysis) and positional enrichment of gene-regulatory histone modifications (ChIP-Seq). In the normal-like cells, nuclear matrix-attached DNA was enriched in expressed genes, while in the breast cancer cells, it was enriched in non-expressed genes. In both cell lines, the chromatin modifications that mark transcriptional activation or repression were appropriately associated with gene expression. Using this new MAR-Seq approach, we provide the first genome-wide characterization of nuclear matrix attachment in mammalian cells and reveal that the nuclear matrix-associated genome is highly cell-context dependent. J. Cell. Physiol. 232: 1295-1305, 2017. © 2016 Wiley Periodicals, Inc.

In Vivo Characterization of Linc-p21 Reveals Functional cis-Regulatory DNA Elements.

The Linc-p21 locus, encoding a long non-coding RNA, plays an important role in p53 signaling, cell-cycle regulation, and tumor suppression. However, despite extensive study, confusion exists regarding its mechanism of action: is activity driven by the transcript acting in trans, in cis, or by an underlying functional enhancer? Here, using a knockout mouse model and a massively parallel enhancer assay, we delineate the functional elements at this locus. We observe that, even in tissues with no detectable Linc-p21 transcript, deletion of the locus significantly affects local gene expression, including of the cell-cycle regulator Cdkn1a. To characterize this RNA-independent regulatory effect, we systematically interrogated the underlying DNA sequence for enhancer activity at nucleotide resolution and confirmed the existence of multiple enhancer elements. Together, these data suggest that, in vivo, the cis-regulatory effects mediated by Linc-p21, in the presence or absence of transcription, are due to DNA enhancer elements.