When chromatin meets splicing

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

PARP1 Regulates Circular RNA Biogenesis though Control of Transcriptional Dynamics

Circular RNAs (circRNAs) are a recently discovered class of RNAs derived from protein-coding genes that have important biological and pathological roles. They are formed through backsplicing during co-transcriptional alternative splicing; however, the unified mechanism that accounts for backsplicing decisions remains unclear. Factors that regulate the transcriptional timing and spatial organization of pre-mRNA, including RNAPII kinetics, the availability of splicing factors, and features of gene architecture, have been shown to influence backsplicing decisions. Poly (ADP-ribose) polymerase I (PARP1) regulates alternative splicing through both its presence on chromatin as well as its PARylation activity. However, no studies have investigated PARP1's possible role in regulating circRNA biogenesis. Here, we hypothesized that PARP1's role in splicing extends to circRNA biogenesis. Our results identify many unique circRNAs in PARP1 depletion and PARylation-inhibited conditions compared to the wild type. We found that while all genes producing circRNAs share gene architecture features common to circRNA host genes, genes producing circRNAs in PARP1 knockdown conditions had longer upstream introns than downstream introns, whereas flanking introns in wild type host genes were symmetrical. Interestingly, we found that the behavior of PARP1 in regulating RNAPII pausing is distinct between these two classes of host genes. We conclude that the PARP1 pausing of RNAPII works within the context of gene architecture to regulate transcriptional kinetics, and therefore circRNA biogenesis. Furthermore, this regulation of PARP1 within host genes acts to fine tune their transcriptional output with implications in gene function.

CHD8 suppression impacts on histone H3 lysine 36 trimethylation and alters RNA alternative splicing

Disruptive mutations in the chromodomain helicase DNA-binding protein 8 gene (CHD8) have been recurrently associated with autism spectrum disorders (ASDs). Here we investigated how chromatin reacts to CHD8 suppression by analyzing a panel of histone modifications in induced pluripotent stem cell-derived neural progenitors. CHD8 suppression led to significant reduction (47.82%) in histone H3K36me3 peaks at gene bodies, particularly impacting on transcriptional elongation chromatin states. H3K36me3 reduction specifically affects highly expressed, CHD8-bound genes and correlates with altered alternative splicing patterns of 462 genes implicated in 'regulation of RNA splicing' and 'mRNA catabolic process'. Mass spectrometry analysis uncovered a novel interaction between CHD8 and the splicing regulator heterogeneous nuclear ribonucleoprotein L (hnRNPL), providing the first mechanistic insights to explain the CHD8 suppression-derived splicing phenotype, partly implicating SETD2, a H3K36me3 methyltransferase. In summary, our results point toward broad molecular consequences of CHD8 suppression, entailing altered histone deposition/maintenance and RNA processing regulation as important regulatory processes in ASD.

The SWI/SNF subunit BRG1 affects alternative splicing by changing RNA binding factor interactions with nascent RNA

BRG1 and BRM are ATPase core subunits of the human SWI/SNF chromatin remodelling complexes mainly associated with transcriptional initiation. They also have a role in alternative splicing, which has been shown for BRM-containing SWI/SNF complexes at a few genes. Here, we have identified a subset of genes which harbour alternative exons that are affected by SWI/SNF ATPases by expressing the ATPases BRG1 and BRM in C33A cells, a BRG1- and BRM-deficient cell line, and analysed the effect on splicing by RNA sequencing. BRG1- and BRM-affected sub-sets of genes favouring both exon inclusion and exon skipping, with only a minor overlap between the ATPase. Some of the changes in alternative splicing induced by BRG1 and BRM expression did not require the ATPase activity. The BRG1-ATPase independent included exons displayed an exon signature of a high GC content. By investigating three genes with exons affected by the BRG-ATPase-deficient variant, we show that these exons accumulated phosphorylated RNA pol II CTD, both serine 2 and serine 5 phosphorylation, without an enrichment of the RNA polymerase II. The ATPases were recruited to the alternative exons, together with both core and signature subunits of SWI/SNF complexes, and promoted the binding of RNA binding factors to chromatin and RNA at the alternative exons. The interaction with the nascent RNP, however, did not reflect the association to chromatin. The hnRNPL, hnRNPU and SAM68 proteins associated with chromatin in cells expressing BRG1 and BRM wild type, but the binding of hnRNPU to the nascent RNP was excluded. This suggests that SWI/SNF can regulate alternative splicing by interacting with splicing-RNA binding factor and influence their binding to the nascent pre-mRNA particle.

KDM3A regulates alternative splicing of cell-cycle genes following DNA damage

Changes in the cellular environment result in chromatin structure alteration, which in turn regulates gene expression. To learn about the effect of the cellular environment on the transcriptome, we studied the H3K9 demethylase KDM3A. Using RNA-seq, we found that KDM3A regulates the transcription and alternative splicing of genes associated with cell cycle and DNA damage. We showed that KDM3A undergoes phosphorylation by PKA at serine 265 following DNA damage, and that the phosphorylation is important for proper cell-cycle regulation. We demonstrated that SAT1 alternative splicing, regulated by KDM3A, plays a role in cell-cycle regulation. Furthermore we found that KDM3A's demethylase activity is not needed for SAT1 alternative splicing regulation. In addition, we identified KDM3A's protein partner ARID1A, the SWI/SNF subunit, and SRSF3 as regulators of SAT1 alternative splicing and showed that KDM3A is essential for SRSF3 binding to SAT1 pre-mRNA. These results suggest that KDM3A serves as a sensor of the environment and an adaptor for splicing factor binding. Our work reveals chromatin sensing of the environment in the regulation of alternative splicing.

Dynamic Alternative Splicing During Mouse Preimplantation Embryo Development

The mechanism of alternative pre-mRNA splicing (AS) during preimplantation development is largely unknown. In order to capture the dynamic changes of AS occurring during embryogenesis, we carried out bioinformatics analysis based on scRNA-seq data over the time-course preimplantation development in mouse. We detected numerous previously-unreported differentially expressed genes at specific developmental stages and investigated the nature of AS at both minor and major zygotic genome activation (ZGA). The AS and differential AS atlas over preimplantation development were established. The differentially alternatively spliced genes (DASGs) are likely to be key splicing factors (SFs) during preimplantation development. We also demonstrated that there is a regulatory cascade of AS events in which some key SFs are regulated by differentially AS of their own gene transcripts. Moreover, 212 isoform switches (ISs) during preimplantation development were detected, which may be critical for decoding the mechanism of early embryogenesis. Importantly, we uncovered that zygotic AS activation (ZASA) is in conformity with ZGA and revealed that AS is coupled with transcription during preimplantation development. Our results may provide a deeper insight into the regulation of early embryogenesis.

DNA methylation directs microRNA biogenesis in mammalian cells

MicroRNA (miRNA) biogenesis initiates co-transcriptionally, but how the Microprocessor machinery pinpoints the locations of short precursor miRNA sequences within long flanking regions of the transcript is not known. Here we show that miRNA biogenesis depends on DNA methylation. When the regions flanking the miRNA coding sequence are highly methylated, the miRNAs are more highly expressed, have greater sequence conservation, and are more likely to drive cancer-related phenotypes than miRNAs encoded by unmethylated loci. We show that the removal of DNA methylation from miRNA loci leads to their downregulation. Further, we found that MeCP2 binding to methylated miRNA loci halts RNA polymerase II elongation, leading to enhanced processing of the primary miRNA by Drosha. Taken together, our data reveal that DNA methylation directly affects miRNA biogenesis.

Histone H1.5 binds over splice sites in chromatin and regulates alternative splicing

Chromatin organization and epigenetic markers influence splicing, though the magnitudes of these effects and the mechanisms are largely unknown. Here, we demonstrate that linker histone H1.5 influences mRNA splicing. We observed that linker histone H1.5 binds DNA over splice sites of short exons in human lung fibroblasts (IMR90 cells). We found that association of H1.5 with these splice sites correlated with the level of inclusion of alternatively spliced exons. Exons marked by H1.5 had more RNA polymerase II (RNAP II) stalling near the 3' splice site than did exons not associated with H1.5. In cells depleted of H1.5, we showed that the inclusion of five exons evaluated decreased and that RNAP II levels over these exons were also reduced. Our findings indicate that H1.5 is involved in regulation of splice site selection and alternative splicing, a function not previously demonstrated for linker histones.

Circular exonic RNAs: When RNA structure meets topology

Although RNA circularization was first documented in the 1990s, the extent to which it occurs was not known until recent advances in high-throughput sequencing enabled the widespread identification of circular RNAs (circRNAs). Despite this, many aspects of circRNA biogenesis, structure, and function yet remain obscure. This review focuses on circular exonic RNAs, a subclass of circRNAs that are generated through backsplicing. Here, I hypothesize that RNA secondary structure can be the common factor that promotes both exon skipping and spliceosomal RNA circularization, and that backsplicing of double-stranded regions could generate topologically linked circRNA molecules. CircRNAs manifest themselves by the presence of tail-to-head exon junctions, which were previously attributed to post-transcriptional exon permutation and repetition. I revisit these observations and argue that backsplicing does not automatically imply RNA circularization because tail-to-head exon junctions give only local information about transcript architecture and, therefore, they are in principle insufficient to determine globally circular topology. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

The chromatin remodeler ZmCHB101 impacts alternative splicing contexts in response to osmotic stress

Maize SWI3-type chromatin remodeler impacts alternative splicing contexts in response to osmotic stress by altering nucleosome density and affecting transcriptional elongation rate. Alternative splicing (AS) is commonly found in higher eukaryotes and is an important posttranscriptional regulatory mechanism to generate transcript diversity. AS has been widely accepted as playing essential roles in different biological processes including growth, development, signal transduction and responses to biotic and abiotic stresses in plants. However, whether and how chromatin remodeling complex functions in AS in plant under osmotic stress remains unknown. Here, we show that a maize SWI3D protein, ZmCHB101, impacts AS contexts in response to osmotic stress. Genome-wide analysis of mRNA contexts in response to osmotic stress using ZmCHB101-RNAi lines reveals that ZmCHB101 impacts alternative splicing contexts of a subset of osmotic stress-responsive genes. Intriguingly, ZmCHB101-mediated regulation of gene expression and AS is largely uncoupled, pointing to diverse molecular functions of ZmCHB101 in transcriptional and posttranscriptional regulation. We further found ZmCHB101 impacts the alternative splicing contexts by influencing alteration of chromatin and histone modification status as well as transcriptional elongation rates mediated by RNA polymerase II. Taken together, our findings suggest a novel insight of how plant chromatin remodeling complex impacts AS under osmotic stress .

Androgen-regulated transcription of ESRP2 drives alternative splicing patterns in prostate cancer

Prostate is the most frequent cancer in men. Prostate cancer progression is driven by androgen steroid hormones, and delayed by androgen deprivation therapy (ADT). Androgens control transcription by stimulating androgen receptor (AR) activity, yet also control pre-mRNA splicing through less clear mechanisms. Here we find androgens regulate splicing through AR-mediated transcriptional control of the epithelial-specific splicing regulator ESRP2. Both ESRP2 and its close paralog ESRP1 are highly expressed in primary prostate cancer. Androgen stimulation induces splicing switches in many endogenous ESRP2-controlled mRNA isoforms, including splicing switches correlating with disease progression. ESRP2 expression in clinical prostate cancer is repressed by ADT, which may thus inadvertently dampen epithelial splice programmes. Supporting this, treatment with the AR antagonist bicalutamide (Casodex) induced mesenchymal splicing patterns of genes including FLNB and CTNND1. Our data reveals a new mechanism of splicing control in prostate cancer with important implications for disease progression.

Alternative splicing links histone modifications to stem cell fate decision

BACKGROUND

Understanding the embryonic stem cell (ESC) fate decision between self-renewal and proper differentiation is important for developmental biology and regenerative medicine. Attention has focused on mechanisms involving histone modifications, alternative pre-messenger RNA splicing, and cell-cycle progression. However, their intricate interrelations and joint contributions to ESC fate decision remain unclear.

RESULTS

We analyze the transcriptomes and epigenomes of human ESC and five types of differentiated cells. We identify thousands of alternatively spliced exons and reveal their development and lineage-dependent characterizations. Several histone modifications show dynamic changes in alternatively spliced exons and three are strongly associated with 52.8% of alternative splicing events upon hESC differentiation. The histone modification-associated alternatively spliced genes predominantly function in G2/M phases and ATM/ATR-mediated DNA damage response pathway for cell differentiation, whereas other alternatively spliced genes are enriched in the G1 phase and pathways for self-renewal. These results imply a potential epigenetic mechanism by which some histone modifications contribute to ESC fate decision through the regulation of alternative splicing in specific pathways and cell-cycle genes. Supported by experimental validations and extended datasets from Roadmap/ENCODE projects, we exemplify this mechanism by a cell-cycle-related transcription factor, PBX1, which regulates the pluripotency regulatory network by binding to NANOG. We suggest that the isoform switch from PBX1a to PBX1b links H3K36me3 to hESC fate determination through the PSIP1/SRSF1 adaptor, which results in the exon skipping of PBX1.

CONCLUSION

We reveal the mechanism by which alternative splicing links histone modifications to stem cell fate decision.

Identification of H4K20me3- and H3K4me3-associated RNAs using CARIP-Seq expands the transcriptional and epigenetic networks of embryonic stem cells

RNA has been shown to interact with various proteins to regulate chromatin dynamics and gene expression. However, it is unknown whether RNAs associate with epigenetic marks such as post-translational modifications of histones, including histone 4 lysine 20 trimethylation (H4K20me3) or trimethylated histone 3 lysine 4 (H3K4me3), to regulate chromatin and gene expression. Here, we used chromatin-associated RNA immunoprecipitation (CARIP) followed by next-generation sequencing (CARIP-Seq) to survey RNAs associated with H4K20me3- and H3K4me3-marked chromatin on a global scale in embryonic stem (ES) cells. We identified thousands of mRNAs and noncoding RNAs that associate with H4K20me3- and H3K4me3-marked chromatin. H4K20me3- and H3K4me3-interacting RNAs are involved in chromatin organization and modification and RNA processing, whereas H4K20me3-only RNAs are involved in cell motility and differentiation, and H3K4me3-only RNAs are involved in metabolic processes and RNA processing. Expression of H3K4me3-associated RNAs is enriched in ES cells, whereas expression of H4K20me3-associated RNAs is enriched in ES cells and differentiated cells. H4K20me3- and H3K4me3-interacting RNAs originate from genes that co-localize with features of active chromatin, including transcriptional machinery and active promoter regions, and the histone modification H3K36me3 in gene body regions. We also found that H4K20me3 and H3K4me3 are associated with distinct gene features including transcripts of greater length and exon number relative to unoccupied transcripts. H4K20me3- and H3K4me3-marked chromatin is also associated with processed RNAs (exon transcripts) relative to unspliced pre-mRNA and ncRNA transcripts. In summary, our results provide evidence that H4K20me3- and H3K4me3-associated RNAs represent a distinct subnetwork of the ES cell transcriptional repertoire.

Cigarette smoke and chewing tobacco alter expression of different sets of miRNAs in oral keratinocytes

Carcinogenic effect of tobacco in oral cancer is through chewing and/or smoking. Significant differences exist in development of oral cancer between tobacco users and non-users. However, molecular alterations induced by different forms of tobacco are yet to be fully elucidated. We developed cellular models of chronic exposure to chewing tobacco and cigarette smoke using immortalized oral keratinocytes. Chronic exposure to tobacco resulted in increased cell scattering and invasiveness in immortalized oral keratinocytes. miRNA sequencing using Illumina HiSeq 2500 resulted in the identification of 10 significantly dysregulated miRNAs (4 fold; p ≤ 0.05) in chewing tobacco treated cells and 6 in cigarette smoke exposed cells. We integrated this data with global proteomic data and identified 36 protein targets that showed inverse expression pattern in chewing tobacco treated cells and 16 protein targets that showed inverse expression in smoke exposed cells. In addition, we identified 6 novel miRNAs in chewing tobacco treated cells and 18 novel miRNAs in smoke exposed cells. Integrative analysis of dysregulated miRNAs and their targets indicates that signaling mechanisms leading to oncogenic transformation are distinct between both forms of tobacco. Our study demonstrates alterations in miRNA expression in oral cells in response to two frequently used forms of tobacco.

Histone posttranslational modifications predict specific alternative exon subtypes in mammalian brain

A compelling body of literature, based on next generation chromatin immunoprecipitation and RNA sequencing of reward brain regions indicates that the regulation of the epigenetic landscape likely underlies chronic drug abuse and addiction. It is now critical to develop highly innovative computational strategies to reveal the relevant regulatory transcriptional mechanisms that may underlie neuropsychiatric disease. We have analyzed chromatin regulation of alternative splicing, which is implicated in cocaine exposure in mice. Recent literature has described chromatin-regulated alternative splicing, suggesting a novel function for drug-induced neuroepigenetic remodeling. However, the extent of the genome-wide association between particular histone modifications and alternative splicing remains unexplored. To address this, we have developed novel computational approaches to model the association between alternative splicing and histone posttranslational modifications in the nucleus accumbens (NAc), a brain reward region. Using classical statistical methods and machine learning to combine ChIP-Seq and RNA-Seq data, we found that specific histone modifications are strongly associated with various aspects of differential splicing. H3K36me3 and H3K4me1 have the strongest association with splicing indicating they play a significant role in alternative splicing in brain reward tissue.

MiasDB: A Database of Molecular Interactions Associated with Alternative Splicing of Human Pre-mRNAs

Alternative splicing (AS) is pervasive in human multi-exon genes and is a major contributor to expansion of the transcriptome and proteome diversity. The accurate recognition of alternative splice sites is regulated by information contained in networks of protein-protein and protein-RNA interactions. However, the mechanisms leading to splice site selection are not fully understood. Although numerous databases have been built to describe AS, molecular interaction databases associated with AS have only recently emerged. In this study, we present a new database, MiasDB, that provides a description of molecular interactions associated with human AS events. This database covers 938 interactions between human splicing factors, RNA elements, transcription factors, kinases and modified histones for 173 human AS events. Every entry includes the interaction partners, interaction type, experimental methods, AS type, tissue specificity or disease-relevant information, a simple description of the functionally tested interaction in the AS event and references. The database can be queried easily using a web server (http://47.88.84.236/Miasdb). We display some interaction figures for several genes. With this database, users can view the regulation network describing AS events for 12 given genes.

Involvement of PARP1 in the regulation of alternative splicing

Specialized chromatin structures such as nucleosomes with specific histone modifications decorate exons in eukaryotic genomes, suggesting a functional connection between chromatin organization and the regulation of pre-mRNA splicing. Through profiling the functional location of Poly (ADP) ribose polymerase, we observed that it is associated with the nucleosomes at exon/intron boundaries of specific genes, suggestive of a role for this enzyme in alternative splicing. Poly (ADP) ribose polymerase has previously been implicated in the PARylation of splicing factors as well as regulation of the histone modification H3K4me3, a mark critical for co-transcriptional splicing. In light of these studies, we hypothesized that interaction of the chromatin-modifying factor, Poly (ADP) ribose polymerase with nucleosomal structures at exon–intron boundaries, might regulate pre-mRNA splicing. Using genome-wide approaches validated by gene-specific assays, we show that depletion of PARP1 or inhibition of its PARylation activity results in changes in alternative splicing of a specific subset of genes. Furthermore, we observed that PARP1 bound to RNA, splicing factors and chromatin, suggesting that Poly (ADP) ribose polymerase serves as a gene regulatory hub to facilitate co-transcriptional splicing. These studies add another function to the multi-functional protein, Poly (ADP) ribose polymerase, and provide a platform for further investigation of this protein’s function in organizing chromatin during gene regulatory processes.

Widespread RNA binding by chromatin-associated proteins

Background

Recent evidence suggests that RNA interaction can regulate the activity and localization of chromatin-associated proteins. However, it is unknown if these observations are specialized instances for a few key RNAs and chromatin factors in specific contexts, or a general mechanism underlying the establishment of chromatin state and regulation of gene expression.

Results

Here, we perform formaldehyde RNA immunoprecipitation (fRIP-Seq) to survey the RNA associated with a panel of 24 chromatin regulators and traditional RNA binding proteins. For each protein that reproducibly bound measurable quantities of bulk RNA (90 % of the panel), we detect enrichment for hundreds to thousands of both noncoding and mRNA transcripts.

Conclusion

For each protein, we find that the enriched sets of RNAs share distinct biochemical, functional, and chromatin properties. Thus, these data provide evidence for widespread specific and relevant RNA association across diverse classes of chromatin-modifying complexes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-0878-3) contains supplementary material, which is available to authorized users.

Global analysis of physical and functional RNA targets of hnRNP L reveals distinct sequence and epigenetic features of repressed and enhanced exons

HnRNP L is a ubiquitous splicing-regulatory protein that is critical for the development and function of mammalian T cells. Previous work has identified a few targets of hnRNP L-dependent alternative splicing in T cells and has described transcriptome-wide association of hnRNP L with RNA. However, a comprehensive analysis of the impact of hnRNP L on mRNA expression remains lacking. Here we use next-generation sequencing to identify transcriptome changes upon depletion of hnRNP L in a model T-cell line. We demonstrate that hnRNP L primarily regulates cassette-type alternative splicing, with minimal impact of hnRNP L depletion on transcript abundance, intron retention, or other modes of alternative splicing. Strikingly, we find that binding of hnRNP L within or flanking an exon largely correlates with exon repression by hnRNP L. In contrast, exons that are enhanced by hnRNP L generally lack proximal hnRNP L binding. Notably, these hnRNP L-enhanced exons share sequence and context features that correlate with poor nucleosome positioning, suggesting that hnRNP may enhance inclusion of a subset of exons via a cotranscriptional or epigenetic mechanism. Our data demonstrate that hnRNP L controls inclusion of a broad spectrum of alternative cassette exons in T cells and suggest both direct RNA regulation as well as indirect mechanisms sensitive to the epigenetic landscape.

Let there be light: regulation of gene expression in plants

Gene expression regulation relies on a variety of molecular mechanisms affecting different steps of a messenger RNA (mRNA) life: transcription, processing, splicing, alternative splicing, transport, translation, storage and decay. Light induces massive reprogramming of gene expression in plants. Differences in alternative splicing patterns in response to environmental stimuli suggest that alternative splicing plays an important role in plant adaptation to changing life conditions. In a recent publication, our laboratories showed that light regulates alternative splicing of a subset of Arabidopsis genes encoding proteins involved in RNA processing by chloroplast retrograde signals. The light effect on alternative splicing is also observed in roots when the communication with the photosynthetic tissues is not interrupted, suggesting that a signaling molecule travels through the plant. These results point at alternative splicing regulation by retrograde signals as an important mechanism for plant adaptation to their environment.

Epigenome-based personalized medicine in human cancer

Cancer genome sequencing has created an opportunity for precision medicine. Thus far, genetic alterations can only be used to guide treatment for small subsets of certain cancer types with these key alterations. Similar to mutations, epigenetic events are equally suitable for personalized medicine. DNA methylation alterations have been used to identify tumor-specific drug responsive markers. Methylation of MGMT sensitizes gliomas to alkylating agents is an example of epigenetic personalized medicine. Recent studies have revealed that 5-azacytidine and decitabine show activity in myelodysplasia, lung and other cancers. There are currently at least 20 kinds of histone deacetylase inhibitors in clinical testing. Inhibitors targeting other epigenetic regulators are being clinically tested, such as EZH2 inhibitor EPZ-6438.

Imprecise intron losses are less frequent than precise intron losses but are not rare in plants

Abstract

In this study, we identified 19 intron losses, including 11 precise intron losses (PILs), six imprecise intron losses (IILs), one de-exonization, and one exon deletion in tomato and potato, and 17 IILs in Arabidopsis thaliana. Comparative analysis of related genomes confirmed that all of the IILs have been fixed during evolution. Consistent with previous studies, our results indicate that PILs are a major type of intron loss. However, at least in plants, IILs are unlikely to be as rare as previously reported.

Reviewers

This article was reviewed by Jun Yu and Zhang Zhang. For complete reviews, see the Reviewers’ Reports section.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0056-7) contains supplementary material, which is available to authorized users.

Genetic interaction mapping reveals a role for the SWI/SNF nucleosome remodeler in spliceosome activation in fission yeast

Although numerous regulatory connections between pre-mRNA splicing and chromatin have been demonstrated, the precise mechanisms by which chromatin factors influence spliceosome assembly and/or catalysis remain unclear. To probe the genetic network of pre-mRNA splicing in the fission yeast Schizosaccharomyces pombe, we constructed an epistatic mini-array profile (E-MAP) and discovered many new connections between chromatin and splicing. Notably, the nucleosome remodeler SWI/SNF had strong genetic interactions with components of the U2 snRNP SF3 complex. Overexpression of SF3 components in ΔSWI/SNF cells led to inefficient splicing of many fission yeast introns, predominantly those with non-consensus splice sites. Deletion of SWI/SNF decreased recruitment of the splicing ATPase Prp2, suggesting that SWI/SNF promotes co-transcriptional spliceosome assembly prior to first step catalysis. Importantly, defects in SWI/SNF as well as SF3 overexpression each altered nucleosome occupancy along intron-containing genes, illustrating that the chromatin landscape both affects—and is affected by—co-transcriptional splicing.

It has recently become apparent that most introns are removed from pre-mRNA while the transcript is still engaged with RNA polymerase II (RNAPII). To gain insight into possible roles for chromatin in co-transcriptional splicing, we generated a genome-wide genetic interaction map in fission yeast and uncovered numerous connections between splicing and chromatin. The SWI/SNF remodeling complex is typically thought to activate gene expression by relieving barriers to polymerase elongation imposed by nucleosomes. Here we show that this remodeler is important for an early step in splicing in which Prp2, an RNA-dependent ATPase, is recruited to the assembling spliceosome to promote catalytic activation. Interestingly, introns with sub-optimal splice sites are particularly dependent on SWI/SNF, suggesting the impact of nucleosome dynamics on the kinetics of spliceosome assembly and catalysis. By monitoring nucleosome occupancy, we show significant alterations in nucleosome density in particular splicing and chromatin mutants, which generally paralleled the levels of RNAPII. Taken together, our findings challenge the notion that nucleosomes simply act as barriers to elongation; rather, we suggest that polymerase pausing at nucleosomes can activate gene expression by allowing more time for co-transcriptional splicing.

DNA-Encoded Chromatin Structural Intron Boundary Signals Identify Conserved Genes with Common Function

The regulation of metazoan gene expression occurs in part by pre-mRNA splicing into mature RNAs. Signals affecting the efficiency and specificity with which introns are removed have not been completely elucidated. Splicing likely occurs cotranscriptionally, with chromatin structure playing a key regulatory role. We calculated DNA encoded nucleosome occupancy likelihood (NOL) scores at the boundaries between introns and exons across five metazoan species. We found that (i) NOL scores reveal a sequence-based feature at the introns on both sides of the intron-exon boundary; (ii) this feature is not part of any recognizable consensus sequence; (iii) this feature is conserved throughout metazoa; (iv) this feature is enriched in genes sharing similar functions: ATPase activity, ATP binding, helicase activity, and motor activity; (v) genes with these functions exhibit different genomic characteristics; (vi) in vivo nucleosome positioning data confirm ontological enrichment at this feature; and (vii) genes with this feature exhibit unique dinucleotide distributions at the intron-exon boundary. The NOL scores point toward a physical property of DNA that may play a role in the mechanism of pre-mRNA splicing. These results provide a foundation for identification of a new set of regulatory DNA elements involved in splicing regulation.

Intronic retroelements: Not just "speed bumps" for RNA polymerase II

Two well-known retroelements, L1 and Alu, comprise about one third of the human genome and are nearly equally distributed between the intergenic and intragenic regions. They carry different regulatory elements and contribute structurally and functionally to the expression of our genes. Recent data also suggest that hundreds of intronic L1s and Alus interfere with the transcription of human genes by inducing intron retention, forcing exonization and cryptic polyadenylation. These novel features can be explained with the RNA polymerase kinetic model and suggest that intronic L1s and Alus are not just "speed bumps" in regulation of RNA polymerase traffic. Here we discuss the complexity of the regulation of gene transcription imposed by intronic retroelements and predict that in addition to transcriptional activity, transcription factor binding and nucleosomal occupancy play a significant role in the transcriptional interference effects of the host genes.

CRE promoter sites modulate alternative splicing via p300-mediated histone acetylation

Histone acetylation modulates alternative splicing of several hundred genes. Here, we tested the role of the histone acetyltransferase p300 in alternative splicing and showed that knockdown of p300 promotes inclusion of the fibronectin (FN1) alternative EDB exon. p300 associates with CRE sites in the promoter via the CREB transcription factor. We created mini-gene reporters driven by an artificial promoter containing CRE sites. Both deletion and mutation of the CRE site affected EDB alternative splicing in the same manner as p300 knockdown. Next we showed that p300 controls histone H4 acetylation along the FN1 gene. Consistently, p300 depletion and CRE deletion/mutation both reduced histone H4 acetylation on mini-gene reporters. Finally, we provide evidence that the effect of CRE inactivation on H4 acetylation and alternative splicing is counteracted by the inhibition of histone deacetylases. Together, these data suggest that histone acetylation could be one of the mechanisms how promoter and promoter binding proteins influence alternative splicing.

Illuminating the Transcriptome through the Genome

Sequencing the human genome was a huge milestone in genetic research that revealed almost the total DNA sequence required to create a human being. However, in order to function, the DNA genome needs to be expressed as an RNA transcriptome. This article reviews how knowledge of genome sequence information has led to fundamental discoveries in how the transcriptome is processed, with a focus on new system-wide insights into how pre-mRNAs that are encoded by split genes in the genome are rearranged by splicing into functional mRNAs. These advances have been made possible by the development of new post-genome technologies to probe splicing patterns. Transcriptome-wide approaches have characterised a "splicing code" that is embedded within and has a significant role in deciphering the genome, and is deciphered by RNA binding proteins. These analyses have also found that most human genes encode multiple mRNA isoforms, and in some cases proteins, leading in turn to a re-assessment of what exactly a gene is. Analysis of the transcriptome has given insights into how the genome is packaged and transcribed, and is helping to explain important aspects of genome evolution.

Message ends: RNA 3' processing and flowering time control

Plants control the time at which they flower in order to ensure reproductive success. This control is underpinned by precision in gene regulation acting through genetically separable pathways. The genetic dissection of this process in the model plant Arabidopsis thaliana has led to the recurrent identification of plant-specific and highly conserved RNA 3' end processing factors required to control flowering by specifically controlling transcription of mRNA encoding the floral repressor FLOWERING LOCUS C (FLC). Here, we review the features of these RNA-processing and RNA-associated proteins, and the complex architecture of coding and non-coding RNA transcription at the FLC locus. We discuss alternative concepts that might explain how these RNA-processing events regulate FLC transcription and hence control flowering time.

Epigenomics and the concept of degeneracy in biological systems

Researchers in the field of epigenomics are developing more nuanced understandings of biological complexity, and exploring the multiple pathways that lead to phenotypic expression. The concept of degeneracy-referring to the multiple pathways that a system recruits to achieve functional plasticity-is an important conceptual accompaniment to the growing body of knowledge in epigenomics. Distinct from degradation, redundancy and dilapidation; degeneracy refers to the plasticity of traits whose function overlaps in some environments, but diverges in others. While a redundant system is composed of repeated identical elements performing the same function, a degenerate system is composed of different elements performing similar or overlapping functions. Here, we describe the degenerate structure of gene regulatory systems from the basic genetic code to flexible epigenomic modifications, and discuss how these structural features have contributed to organism complexity, robustness, plasticity and evolvability.

Reprogramming cellular events by poly(ADP-ribose)-binding proteins

Poly(ADP-ribosyl)ation is a posttranslational modification catalyzed by the poly(ADP-ribose) polymerases (PARPs). These enzymes covalently modify glutamic, aspartic and lysine amino acid side chains of acceptor proteins by the sequential addition of ADP-ribose (ADPr) units. The poly(ADP-ribose) (pADPr) polymers formed alter the physico-chemical characteristics of the substrate with functional consequences on its biological activities. Recently, non-covalent binding to pADPr has emerged as a key mechanism to modulate and coordinate several intracellular pathways including the DNA damage response, protein stability and cell death. In this review, we describe the basis of non-covalent binding to pADPr that has led to the emerging concept of pADPr-responsive signaling pathways. This review emphasizes the structural elements and the modular strategies developed by pADPr-binding proteins to exert a fine-tuned control of a variety of pathways. Poly(ADP-ribosyl)ation reactions are highly regulated processes, both spatially and temporally, for which at least four specialized pADPr-binding modules accommodate different pADPr structures and reprogram protein functions. In this review, we highlight the role of well-characterized and newly discovered pADPr-binding modules in a diverse set of physiological functions.

Coupling pre-mRNA processing to transcription on the RNA factory assembly line

It has been well-documented that nuclear processing of primary transcripts of RNA polymerase II occurs co-transcriptionally and is functionally coupled to transcription. Moreover, increasing evidence indicates that transcription influences pre-mRNA splicing and even several post-splicing RNA processing events. In this review, we discuss the issues of how RNA polymerase II modulates co-transcriptional RNA processing events via its carboxyl terminal domain, and the protein domains involved in coupling of transcription and RNA processing events. In addition, we describe how transcription influences the expression or stability of mRNAs through the formation of distinct mRNP complexes. Finally, we delineate emerging findings that chromatin modifications function in the regulation of RNA processing steps, especially splicing, in addition to transcription. Overall, we provide a comprehensive view that transcription could integrate different control systems, from epigenetic to post-transcriptional control, for efficient gene expression.

Function of alternative splicing

Almost all polymerase II transcripts undergo alternative pre-mRNA splicing. Here, we review the functions of alternative splicing events that have been experimentally determined. The overall function of alternative splicing is to increase the diversity of mRNAs expressed from the genome. Alternative splicing changes proteins encoded by mRNAs, which has profound functional effects. Experimental analysis of these protein isoforms showed that alternative splicing regulates binding between proteins, between proteins and nucleic acids as well as between proteins and membranes. Alternative splicing regulates the localization of proteins, their enzymatic properties and their interaction with ligands. In most cases, changes caused by individual splicing isoforms are small. However, cells typically coordinate numerous changes in 'splicing programs', which can have strong effects on cell proliferation, cell survival and properties of the nervous system. Due to its widespread usage and molecular versatility, alternative splicing emerges as a central element in gene regulation that interferes with almost every biological function analyzed.

Connections between chromatin signatures and splicing

Splicing and alternative splicing are involved in the expression of most human genes, playing key roles in differentiation, cell cycle progression, and development. Misregulation of splicing is frequently associated to disease, which imposes a better understanding of the mechanisms underlying splicing regulation. Accumulated evidence suggests that multiple trans-acting factors and cis-regulatory elements act together to determine tissue-specific splicing patterns. Besides, as splicing is often cotranscriptional, a complex picture emerges in which splicing regulation not only depends on the balance of splicing factor binding to their pre-mRNA target sites but also on transcription-associated features such as protein recruitment to the transcribing machinery and elongation kinetics. Adding more complexity to the splicing regulation network, recent evidence shows that chromatin structure is another layer of regulation that may act through various mechanisms. These span from regulation of RNA polymerase II elongation, which ultimately determines splicing decisions, to splicing factor recruitment by specific histone marks. Chromatin may not only be involved in alternative splicing regulation but in constitutive exon recognition as well. Moreover, splicing was found to be necessary for the proper 'writing' of particular chromatin signatures, giving further mechanistic support to functional interconnections between splicing, transcription and chromatin structure. These links between chromatin configuration and splicing raise the intriguing possibility of the existence of a memory for splicing patterns to be inherited through epigenetic modifications.

Differential GC content between exons and introns establishes distinct strategies of splice-site recognition

During evolution segments of homeothermic genomes underwent a GC content increase. Our analyses reveal that two exon-intron architectures have evolved from an ancestral state of low GC content exons flanked by short introns with a lower GC content. One group underwent a GC content elevation that abolished the differential exon-intron GC content, with introns remaining short. The other group retained the overall low GC content as well as the differential exon-intron GC content, and is associated with longer introns. We show that differential exon-intron GC content regulates exon inclusion level in this group, in which disease-associated mutations often lead to exon skipping. This group's exons also display higher nucleosome occupancy compared to flanking introns and exons of the other group, thus "marking" them for spliceosomal recognition. Collectively, our results reveal that differential exon-intron GC content is a previously unidentified determinant of exon selection and argue that the two GC content architectures reflect the two mechanisms by which splicing signals are recognized: exon definition and intron definition.

Epigenetic features are significantly associated with alternative splicing

BACKGROUND

While alternative splicing (AS) contributes greatly to protein diversities, the relationship between various types of AS and epigenetic factors remains largely unknown.

RESULTS

In this study, we discover that a number of epigenetic features, including DNA methylation, nucleosome occupancy, specific histone modifications and protein features, are strongly associated with AS. To further enhance our understanding of the association between these features and AS, we cluster our investigated features based on their association patterns with each AS type into four groups, with H3K36me3, EGR1, GABP, SRF, SIN3A and RNA Pol II grouped together and showing strongest association with AS. In addition, we find that the AS types can be classified into two general classes, namely the exon skipping related process (ESRP), and the alternative splice site selection process (ASSP), based on their association levels with the epigenetic features.

CONCLUSION

Our analysis thus suggests that epigenetic features are likely to play important roles in regulating AS.

Nuclear quality control of RNA polymerase II transcripts

Eukaryotic RNA polymerase II produces an astounding diversity of transcripts. These may need to be 5(') capped, spliced, polyadenylated, and packaged with proteins before their export to the cytoplasm. Unscheduled accumulation of any RNA species can interfere with normal RNA metabolism and poses a serious hazard to cells. Yet, given the amount of primary transcripts and the complexity of the RNA maturation process, production of aberrant RNA species is unavoidable. Cells, therefore, employ nuclear RNA quality control mechanisms to rapidly degrade, actively retain, or transcriptionally silence unwanted RNAs. Pathways that monitor mRNA production are best understood and similar pathways are employed to destroy transcriptional noise. Finally, related mechanisms also contribute to gene regulation during normal growth.

Proteome analysis of protein partners to nucleosomes containing canonical H2A or the variant histones H2A.Z or H2A.X

Although the existence of histone variants has been known for quite some time, only recently are we grasping the breadth and diversity of the cellular processes in which they are involved. Of particular interest are the two variants of histone H2A, H2A.Z and H2A.X because of their roles in regulation of gene expression and in DNA double-strand break repair, respectively. We hypothesize that nucleosomes containing these variants may perform their distinct functions by interacting with different sets of proteins. Here, we present our proteome analysis aimed at identifying protein partners that interact with nucleosomes containing H2A.Z, H2A.X or their canonical H2A counterpart. Our development of a nucleosome-pull down assay and analysis of the recovered nucleosome-interacting proteins by mass spectrometry allowed us to directly compare nuclear partners of these variant-containing nucleosomes to those containing canonical H2A. To our knowledge, our data represent the first systematic analysis of the H2A.Z and H2A.X interactome in the context of nucleosome structure.

Depolarization-mediated regulation of alternative splicing

Alternative splicing in eukaryotes plays an important role in regulating gene expression by selectively including alternative exons. A wealth of information has been accumulated that explains how alternative exons are selected in a developmental stage- or tissue-specific fashion. However, our knowledge of how cells respond to environmental changes to alter alternative splicing is very limited. For example, although a number of alternative exons have been shown to be regulated by calcium level alterations, the underlying mechanisms are not well understood. As calcium signaling in neurons plays a crucial role in essential neuronal functions such as learning and memory formation, it is important to understand how this process is regulated at every level in gene expression. The significance of the dynamic control of alternative splicing in response to changes of calcium levels has been largely unappreciated. In this communication, we will summarize the recent advances in calcium signaling-mediated alternative splicing that have provided some insights into the important regulatory mechanisms. In addition to describing the cis-acting RNA elements on the pre-mRNA molecules that respond to changes of intracellular calcium levels, we will summarize how splicing regulators change and affect alternative splicing in this process. We will also discuss a novel mode of calcium-mediated splicing regulation at the level of chromatin structure and transcription.

Son maintains accurate splicing for a subset of human pre-mRNAs

Serine-arginine-rich (SR) proteins play a key role in alternative pre-mRNA splicing in eukaryotes. We recently showed that a large SR protein called Son has unique repeat motifs that are essential for maintaining the subnuclear organization of pre-mRNA processing factors in nuclear speckles. Motif analysis of Son highlights putative RNA interaction domains that suggest a direct role for Son in pre-mRNA splicing. Here, we used in situ approaches to show that Son localizes to a reporter minigene transcription site, and that RNAi-mediated Son depletion causes exon skipping on reporter transcripts at this transcription site. A genome-wide exon microarray analysis was performed to identify human transcription and splicing targets of Son. Our data show that Son-regulated splicing encompasses all known types of alternative splicing, the most common being alternative splicing of cassette exons. We confirmed that knockdown of Son leads to exon skipping in pre-mRNAs for chromatin-modifying enzymes, including ADA, HDAC6 and SetD8. This study reports a comprehensive view of human transcription and splicing targets for Son in fundamental cellular pathways such as integrin-mediated cell adhesion, cell cycle regulation, cholesterol biosynthesis, apoptosis and epigenetic regulation of gene expression.

Genome-wide "re"-modeling of nucleosome positions

Two recent studies mapped nucleosomes across the yeast and human genomes, teasing apart the relative contributions of DNA sequence and chromatin remodelers to nucleosome organization. These data suggest two emerging models: chromatin remodelers position nucleosomes around transcriptional start sites in yeast, and a few "locked" nucleosomes may serve as barriers from which nucleosome arrays emanate in human genomes.

Histone H3 lysine 9 trimethylation and HP1γ favor inclusion of alternative exons

Pre-messenger RNAs (pre-mRNAs) maturation is initiated cotranscriptionally. It is therefore conceivable that chromatin-borne information participates in alternative splicing. Here we find that elevated levels of trimethylation of histone H3 on Lys9 (H3K9me3) are a characteristic of the alternative exons of several genes including CD44. On this gene the chromodomain protein HP1γ, frequently defined as a transcriptional repressor, facilitates inclusion of the alternative exons via a mechanism involving decreased RNA polymerase II elongation rate. In addition, accumulation of HP1γ on the variant region of the CD44 gene stabilizes association of the pre-mRNA with the chromatin. Altogether, our data provide evidence for localized histone modifications impacting alternative splicing. They further implicate HP1γ as a possible bridging molecule between the chromatin and the maturating mRNA, with a general impact on splicing decisions.

Cross-talk in transcription, splicing and chromatin: who makes the first call?

The complex processes of mRNA transcription and splicing were traditionally studied in isolation. In vitro studies showed that splicing could occur independently of transcription and the perceived wisdom was that, to a large extent, it probably did. However, there is now abundant evidence for functional interactions between transcription and splicing, with important consequences for splicing regulation. In the present paper, we summarize the evidence that transcription affects splicing and vice versa, and the more recent indications of epigenetic effects on splicing, through chromatin modifications. We end by discussing the potential for a systems biology approach to obtain better insight into how these processes affect each other.

Charting histone modifications and the functional organization of mammalian genomes

A succession of technological advances over the past decade have enabled researchers to chart maps of histone modifications and related chromatin structures with increasing accuracy, comprehensiveness and throughput. The resulting data sets highlight the interplay between chromatin and genome function, dynamic variations in chromatin structure across cellular conditions, and emerging roles for large-scale domains and higher-ordered chromatin organization. Here we review a selection of recent studies that have probed histone modifications and successive layers of chromatin structure in mammalian genomes, the patterns that have been identified and future directions for research.

How do RNA sequence, DNA sequence, and chromatin properties regulate splicing?

Recent genome-wide studies have revealed a remarkable correspondence between nucleosome positions and exon-intron boundaries, and several studies have implicated specific histone modifications in regulating alternative splicing. In addition, recent progress in cracking the ‘splicing code’ shows that sequence motifs carried on the nascent RNA molecule itself are sufficient to accurately predict tissue-specific alternative splicing patterns. Together, these studies shed light on the complex interplay between RNA sequence, DNA sequence, and chromatin properties in regulating splicing.

Splicing factor and exon profiling across human tissues

It has been shown that alternative splicing is especially prevalent in brain and testis when compared to other tissues. To test whether there is a specific propensity of these tissues to generate splicing variants, we used a single source of high-density microarray data to perform both splicing factor and exon expression profiling across 11 normal human tissues. Paired comparisons between tissues and an original exon-based statistical group analysis demonstrated after extensive RT-PCR validation that the cerebellum, testis, and spleen had the largest proportion of differentially expressed alternative exons. Variations at the exon level correlated with a larger number of splicing factors being expressed at a high level in the cerebellum, testis and spleen than in other tissues. However, this splicing factor expression profile was similar to a more global gene expression pattern as a larger number of genes had a high expression level in the cerebellum, testis and spleen. In addition to providing a unique resource on expression profiling of alternative splicing variants and splicing factors across human tissues, this study demonstrates that the higher prevalence of alternative splicing in a subset of tissues originates from the larger number of genes, including splicing factors, being expressed than in other tissues.