Splicing enhances recruitment of methyltransferase HYPB/Setd2 and methylation of histone H3 Lys36

A computational method for studying the relation between alternative splicing and DNA methylation

Alternative splicing is an important mechanism in eukaryotes that expands the transcriptome and proteome significantly. It plays an important role in a number of biological processes. Understanding its regulation is hence an important challenge. Recently, increasing evidence has been collected that supports an involvement of intragenic DNA methylation in the regulation of alternative splicing. The exact mechanisms of regulation, however, are largely unknown, and speculated to be complex: different methylation profiles might exist, each of which could be associated with a different regulation mechanism. We present a computational technique that is able to determine such stable methylation patterns and allows to correlate these patterns with inclusion propensity of exons. Pattern detection is based on dynamic time warping (DTW) of methylation profiles, a sophisticated similarity measure for signals that can be non-trivially transformed. We design a flexible self-organizing map approach to pattern grouping. Exemplary application on available data sets indicates that stable patterns which correlate non-trivially with exon inclusion do indeed exist. To improve the reliability of these predictions, further studies on larger data sets will be required. We have thus taken great care that our software runs efficiently on modern hardware, so that it can support future studies on large-scale data sets.

Absence of canonical marks of active chromatin in developmentally regulated genes

The interplay of active and repressive histone modifications is assumed to have a key role in the regulation of gene expression. In contrast to this generally accepted view, we show that the transcription of genes temporally regulated during fly and worm development occurs in the absence of canonically active histone modifications. Conversely, strong chromatin marking is related to transcriptional and post-transcriptional stability, an association that we also observe in mammals. Our results support a model in which chromatin marking is associated with the stable production of RNA, whereas unmarked chromatin would permit rapid gene activation and deactivation during development. In the latter case, regulation by transcription factors would have a comparatively more important regulatory role than chromatin marks.

Connecting the dots: chromatin and alternative splicing in EMT

Nature has devised sophisticated cellular machinery to process mRNA transcripts produced by RNA Polymerase II, removing intronic regions and connecting exons together, to produce mature RNAs. This process, known as splicing, is very closely linked to transcription. Alternative splicing, or the ability to produce different combinations of exons that are spliced together from the same genomic template, is a fundamental means of regulating protein complexity. Similar to transcription, both constitutive and alternative splicing can be regulated by chromatin and its associated factors in response to various signal transduction pathways activated by external stimuli. This regulation can vary between different cell types, and interference with these pathways can lead to changes in splicing, often resulting in aberrant cellular states and disease. The epithelial to mesenchymal transition (EMT), which leads to cancer metastasis, is influenced by alternative splicing events of chromatin remodelers and epigenetic factors such as DNA methylation and non-coding RNAs. In this review, we will discuss the role of epigenetic factors including chromatin, chromatin remodelers, DNA methyltransferases, and microRNAs in the context of alternative splicing, and discuss their potential involvement in alternative splicing during the EMT process.

Lessons from modENCODE

The modENCODE (Model Organism Encyclopedia of DNA Elements) Consortium aimed to map functional elements-including transcripts, chromatin marks, regulatory factor binding sites, and origins of DNA replication-in the model organisms Drosophila melanogaster and Caenorhabditis elegans. During its five-year span, the consortium conducted more than 2,000 genome-wide assays in developmentally staged animals, dissected tissues, and homogeneous cell lines. Analysis of these data sets provided foundational insights into genome, epigenome, and transcriptome structure and the evolutionary turnover of regulatory pathways. These studies facilitated a comparative analysis with similar data types produced by the ENCODE Consortium for human cells. Genome organization differs drastically in these distant species, and yet quantitative relationships among chromatin state, transcription, and cotranscriptional RNA processing are deeply conserved. Of the many biological discoveries of the modENCODE Consortium, we highlight insights that emerged from integrative studies. We focus on operational and scientific lessons that may aid future projects of similar scale or aims in other, emerging model systems.

Deciphering the Epigenetic Code of Cardiac Myocyte Transcription

RATIONALE

Epigenetic mechanisms are crucial for cell identity and transcriptional control. The heart consists of different cell types, including cardiac myocytes, endothelial cells, fibroblasts, and others. Therefore, cell type-specific analysis is needed to gain mechanistic insight into the regulation of gene expression in cardiac myocytes. Although cytosolic mRNA represents steady-state levels, nuclear mRNA more closely reflects transcriptional activity. To unravel epigenetic mechanisms of transcriptional control, cell type-specific analysis of nuclear mRNA and epigenetic modifications is crucial.

OBJECTIVE

The aim was to purify cardiac myocyte nuclei from hearts of different species by magnetic- or fluorescent-assisted sorting and to determine the nuclear and cellular RNA expression profiles and epigenetic marks in a cardiac myocyte-specific manner.

METHODS AND RESULTS

Frozen cardiac tissue samples were used to isolate cardiac myocyte nuclei. High sorting purity was confirmed for cardiac myocyte nuclei isolated from mice, rats, and humans. Deep sequencing of nuclear RNA revealed a major fraction of nascent, unspliced RNA in contrast to results obtained from purified cardiac myocytes. Cardiac myocyte nuclear and cellular RNA expression profiles showed differences, especially for metabolic genes. Genome-wide maps of the transcriptional elongation mark H3K36me3 were generated by chromatin-immunoprecipitation. Transcriptome and epigenetic data confirmed the high degree of cardiac myocyte-specificity of our protocol. An integrative analysis of nuclear mRNA and histone mark occurrence indicated a major impact of the chromatin state on transcriptional activity in cardiac myocytes.

CONCLUSIONS

This study establishes cardiac myocyte-specific sorting of nuclei as a universal method to investigate epigenetic and transcriptional processes in cardiac myocytes of different origins. These data sets provide novel insight into cardiac myocyte transcription.

Trimethylation of Lys36 on H3 restricts gene expression change during aging and impacts life span

Pu et al. profiled the genome-wide pattern of trimethylation of Lys36 on histone 3 (H3K36me3) in the somatic cells of young and old C. elegans. Genes with dramatic expression change during aging are marked with low or undetectable levels of H3K36me3 in their gene bodies. A similar negative correlation between H3K36me3 marking and mRNA expression change during aging is also observed in Drosophila melanogaster. Inactivation of the methyltransferase met-1 resulted in a decrease in global H3K36me3 marks, an increase in mRNA expression change with age, and a shortened life span.

Functional data indicate that specific histone modification enzymes can be key to longevity in Caenorhabditis elegans, but the molecular basis of how chromatin structure modulates longevity is not well understood. In this study, we profiled the genome-wide pattern of trimethylation of Lys36 on histone 3 (H3K36me3) in the somatic cells of young and old Caenorhabditis elegans. We revealed a new role of H3K36me3 in maintaining gene expression stability through aging with important consequences on longevity. We found that genes with dramatic expression change during aging are marked with low or even undetectable levels of H3K36me3 in their gene bodies irrespective of their corresponding mRNA abundance. Interestingly, 3′ untranslated region (UTR) length strongly correlates with H3K36me3 levels and age-dependent mRNA expression stability. A similar negative correlation between H3K36me3 marking and mRNA expression change during aging was also observed in Drosophila melanogaster, suggesting a conserved mechanism for H3K36me3 in suppressing age-dependent mRNA expression change. Importantly, inactivation of the methyltransferase met-1 resulted in a decrease in global H3K36me3 marks, an increase in mRNA expression change with age, and a shortened life span, suggesting a causative role of the H3K36me3 marking in modulating age-dependent gene expression stability and longevity.

A chromatin code for alternative splicing involving a putative association between CTCF and HP1α proteins

Background

Alternative splicing is primarily controlled by the activity of splicing factors and by the elongation of the RNA polymerase II (RNAPII). Recent experiments have suggested a new complex network of splicing regulation involving chromatin, transcription and multiple protein factors. In particular, the CCCTC-binding factor (CTCF), the Argonaute protein AGO1, and members of the heterochromatin protein 1 (HP1) family have been implicated in the regulation of splicing associated with chromatin and the elongation of RNAPII. These results raise the question of whether these proteins may associate at the chromatin level to modulate alternative splicing.

Results

Using chromatin immunoprecipitation sequencing (ChIP-Seq) data for CTCF, AGO1, HP1α, H3K27me3, H3K9me2, H3K36me3, RNAPII, total H3 and 5metC and alternative splicing arrays from two cell lines, we have analyzed the combinatorial code of their binding to chromatin in relation to the alternative splicing patterns between two cell lines, MCF7 and MCF10. Using Machine Learning techniques, we identified the changes in chromatin signals that are most significantly associated with splicing regulation between these two cell lines. Moreover, we have built a map of the chromatin signals on the pre-mRNA, that is, a chromatin-based RNA-map, which can explain 606 (68.55%) of the regulated events between MCF7 and MCF10. This chromatin code involves the presence of HP1α, CTCF, AGO1, RNAPII and histone marks around regulated exons and can differentiate patterns of skipping and inclusion. Additionally, we found a significant association of HP1α and CTCF activities around the regulated exons and a putative DNA binding site for HP1α.

Conclusions

Our results show that a considerable number of alternative splicing events could have a chromatin-dependent regulation involving the association of HP1α and CTCF near regulated exons. Additionally, we find further evidence for the involvement of HP1α and AGO1 in chromatin-related splicing regulation.

Electronic supplementary material

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

Quantification of co-transcriptional splicing from RNA-Seq data

During gene expression, protein-coding transcripts are shaped by multiple processing events: 5' end capping, pre-mRNA splicing, RNA editing, and 3' end cleavage and polyadenylation. These events are required to produce mature mRNA, which can be subsequently translated. Nearly all of these RNA processing steps occur during transcription, while the nascent RNA is still attached to the DNA template by RNA polymerase II (i.e. co-transcriptionally). Polyadenylation occurs after 3' end cleavage or post-transcriptionally. Pre-mRNA splicing - the removal of introns and ligation of exons - can be initiated and concluded co-transcriptionally, although this is not strictly required. Recently, a number of studies using global methods have shown that the majority of splicing is co-transcriptional, yet not all published studies agree in their conclusions. Short read sequencing of RNA (RNA-Seq) is the prevailing approach to measuring splicing levels in nascent RNA, mRNA or total RNA. Here, we compare four different strategies for analyzing and quantifying co-transcriptional splicing. To do so, we reanalyze two nascent RNA-Seq datasets of the same species, but different cell type and RNA isolation procedure. Average co-transcriptional splicing values calculated on a per intron basis are similar, independent of the strategy used. We emphasize the technical requirements for identifying co-transcriptional splicing events with high confidence, e.g. how to calculate co-transcriptional splicing from nascent RNA- versus mRNA-Seq data, the number of biological replicates needed, depletion of polyA+RNA, and appropriate normalization. Finally, we present guidelines for planning a nascent RNA-Seq experiment.

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.

Co-transcriptional mRNP formation is coordinated within a molecular mRNP packaging station in S. cerevisiae

In eukaryotes, the messenger RNA (mRNA), the blueprint of a protein‐coding gene, is processed and packaged into a messenger ribonucleoprotein particle (mRNP) by mRNA‐binding proteins in the nucleus. The steps of mRNP formation – transcription, processing, packaging, and the orchestrated release of the export‐competent mRNP from the site of transcription for nuclear mRNA export – are tightly coupled to ensure a highly efficient and regulated process. The importance of highly accurate nuclear mRNP formation is illustrated by the fact that mutations in components of this pathway lead to cellular inviability or to severe diseases in metazoans. We hypothesize that efficient mRNP formation is realized by a molecular mRNP packaging station, which is built by several recruitment platforms and coordinates the individual steps of mRNP formation.

Alternative splicing regulates the expression of G9A and SUV39H2 methyltransferases, and dramatically changes SUV39H2 functions

Alternative splicing is the main source of proteome diversity. Here, we have investigated how alternative splicing affects the function of two human histone methyltransferases (HMTase): G9A and SUV39H2. We show that exon 10 in G9A and exon 3 in SUV39H2 are alternatively included in a variety of tissues and cell lines, as well as in a different species. The production of these variants is likely tightly regulated because both constitutive and alternative splicing factors control their splicing profiles. Based on this evidence, we have assessed the link between the inclusion of these exons and the activity of both enzymes. We document that these HMTase genes yield several protein isoforms, which are likely issued from alternative splicing regulation. We demonstrate that inclusion of SUV39H2 exon 3 is a determinant of the stability, the sub-nuclear localization, and the HMTase activity. Genome-wide expression analysis further revealed that alternative inclusion of SUV39H2 exon 3 differentially modulates the expression of target genes. Our data also suggest that a variant of G9A may display a function that is independent of H3K9 methylation. Our work emphasizes that expression and function of genes are not collinear; therefore alternative splicing must be taken into account in any functional study.

Chromatin modification by the RNA Polymerase II elongation complex

Transcription elongation by RNA polymerase II (RNAP II) involves the coordinated action of numerous regulatory factors. Among these are chromatin-modifying enzymes, which generate a stereotypic and conserved pattern of histone modifications along transcribed genes. This pattern implies a precise coordination between regulators of histone modification and the RNAP II elongation complex. Here I review the pathways and molecular events that regulate co-transcriptional histone modifications. Insight into these events will illuminate the assembly of functional RNAP II elongation complexes and how the chromatin landscape influences their composition and function.

Selective inhibitors of protein methyltransferases

Mounting evidence suggests that protein methyltransferases (PMTs), which catalyze methylation of histone and nonhistone proteins, play a crucial role in diverse biological processes and human diseases. In particular, PMTs have been recognized as major players in regulating gene expression and chromatin state. PMTs are divided into two categories: protein lysine methyltransferases (PKMTs) and protein arginine methyltransferases (PRMTs). There has been a steadily growing interest in these enzymes as potential therapeutic targets and therefore discovery of PMT inhibitors has also been pursued increasingly over the past decade. Here, we present a perspective on selective, small-molecule inhibitors of PMTs with an emphasis on their discovery, characterization, and applicability as chemical tools for deciphering the target PMTs' physiological functions and involvement in human diseases. We highlight the current state of PMT inhibitors and discuss future directions and opportunities for PMT inhibitor discovery.

Histone H3.3 and cancer: A potential reader connection

The building block of chromatin is nucleosome, which consists of 146 base pairs of DNA wrapped around a histone octamer composed of two copies of histone H2A, H2B, H3, and H4. Significantly, the somatic missense mutations of the histone H3 variant, H3.3, are associated with childhood and young-adult tumors, such as pediatric high-grade astrocytomas, as well as chondroblastoma and giant-cell tumors of the bone. The mechanisms by which these histone mutations cause cancer are by and large unclear. Interestingly, two recent studies identified BS69/ZMYND11, which was proposed to be a candidate tumor suppressor, as a specific reader for a modified form of H3.3 (H3.3K36me3). Importantly, some H3.3 cancer mutations are predicted to abrogate the H3.3K36me3/BS69 interaction, suggesting that this interaction may play an important role in tumor suppression. These new findings also raise the question of whether H3.3 cancer mutations may lead to the disruption and/or gain of interactions of additional cellular factors that contribute to tumorigenesis.

Identification by high-throughput imaging of the histone methyltransferase EHMT2 as an epigenetic regulator of VEGFA alternative splicing

Recent evidence points to a role of chromatin in regulation of alternative pre-mRNA splicing (AS). In order to identify novel chromatin regulators of AS, we screened an RNAi library of chromatin proteins using a cell-based high-throughput in vivo assay. We identified a set of chromatin proteins that regulate AS. Using simultaneous genome-wide expression and AS analysis, we demonstrate distinct and non-overlapping functions of these chromatin modifiers on transcription and AS. Detailed mechanistic characterization of one dual function chromatin modifier, the H3K9 methyltransferase EHMT2 (G9a), identified VEGFA as a major chromatin-mediated AS target. Silencing of EHMT2, or its heterodimer partner EHMT1, affects AS by promoting exclusion of VEGFA exon 6a, but does not alter total VEGFA mRNA levels. The epigenetic regulatory mechanism of AS by EHMT2 involves an adaptor system consisting of the chromatin modulator HP1γ, which binds methylated H3K9 and recruits splicing regulator SRSF1. The epigenetic regulation of VEGFA is physiologically relevant since EHMT2 is transcriptionally induced in response to hypoxia and triggers concomitant changes in AS of VEGFA. These results characterize a novel epigenetic regulatory mechanism of AS and they demonstrate separate roles of epigenetic modifiers in transcription and alternative splicing.

Introns and gene expression: cellular constraints, transcriptional regulation, and evolutionary consequences

A gene's “expression profile” denotes the number of transcripts present relative to all other transcripts. The overall rate of transcript production is determined by transcription and RNA processing rates. While the speed of elongating RNA polymerase II has been characterized for many different genes and organisms, gene-architectural features – primarily the number and length of exons and introns – have recently emerged as important regulatory players. Several new studies indicate that rapidly cycling cells constrain gene-architecture toward short genes with a few introns, allowing efficient expression during short cell cycles. In contrast, longer genes with long introns exhibit delayed expression, which can serve as timing mechanisms for patterning processes. These findings indicate that cell cycle constraints drive the evolution of gene-architecture and shape the transcriptome of a given cell type. Furthermore, a tendency for short genes to be evolutionarily young hints at links between cellular constraints and the evolution of animal ontogeny.

Argonaute-1 binds transcriptional enhancers and controls constitutive and alternative splicing in human cells

The roles of Argonaute proteins in cytoplasmic microRNA and RNAi pathways are well established. However, their implication in small RNA-mediated transcriptional gene silencing in the mammalian cell nucleus is less understood. We have recently shown that intronic siRNAs cause chromatin modifications that inhibit RNA polymerase II elongation and modulate alternative splicing in an Argonaute-1 (AGO1)-dependent manner. Here we used chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) to investigate the genome-wide distribution of AGO1 nuclear targets. Unexpectedly, we found that about 80% of AGO1 clusters are associated with cell-type-specific transcriptional enhancers, most of them (73%) overlapping active enhancers. This association seems to be mediated by long, rather than short, enhancer RNAs and to be more prominent in intragenic, rather than intergenic, enhancers. Paradoxically, crossing ChIP-seq with RNA-seq data upon AGO1 depletion revealed that enhancer-bound AGO1 is not linked to the global regulation of gene transcription but to the control of constitutive and alternative splicing, which was confirmed by an individual gene analysis explaining how AGO1 controls inclusion levels of the cassette exon 107 in the SYNE2 gene.

Cotranscriptional histone H2B monoubiquitylation is tightly coupled with RNA polymerase II elongation rate

Various histone modifications decorate nucleosomes within transcribed genes. Among these, monoubiquitylation of histone H2B (H2Bub1) and methylation of histone H3 on lysines 36 (H3K36me2/3) and 79 (H3K79me2/3) correlate positively with gene expression. By measuring the progression of the transcriptional machinery along genes within live cells, we now report that H2B monoubiquitylation occurs cotranscriptionally and accurately reflects the advance of RNA polymerase II (Pol II). In contrast, H3K36me3 and H3K79me2 are less dynamic and represent Pol II movement less faithfully. High-resolution ChIP-seq reveals that H2Bub1 levels are selectively reduced at exons and decrease in an exon-dependent stepwise manner toward the 3' end of genes. Exonic depletion of H2Bub1 in gene bodies is highly correlated with Pol II pausing at exons, suggesting elongation rate changes associated with intron-exon structure. In support of this notion, H2Bub1 levels were found to be significantly correlated with transcription elongation rates measured in various cell lines. Overall, our data shed light on the organization of H2Bub1 within transcribed genes and single out H2Bub1 as a reliable marker for ongoing transcription elongation.

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.

Histone modifications are associated with transcript isoform diversity in normal and cancer cells

Mechanisms that generate transcript diversity are of fundamental importance in eukaryotes. Although a large fraction of human protein-coding genes and lincRNAs produce more than one mRNA isoform each, the regulation of this phenomenon is still incompletely understood. Much progress has been made in deciphering the role of sequence-specific features as well as DNA-and RNA-binding proteins in alternative splicing. Recently, however, several experimental studies of individual genes have revealed a direct involvement of epigenetic factors in alternative splicing and transcription initiation. While histone modifications are generally correlated with overall gene expression levels, it remains unclear how histone modification enrichment affects relative isoform abundance. Therefore, we sought to investigate the associations between histone modifications and transcript diversity levels measured by the rates of transcription start-site switching and alternative splicing on a genome-wide scale across protein-coding genes and lincRNAs. We found that the relationship between enrichment levels of epigenetic marks and transcription start-site switching is similar for protein-coding genes and lincRNAs. Furthermore, we found associations between splicing rates and enrichment levels of H2az, H3K4me1, H3K4me2, H3K4me3, H3K9ac, H3K9me3, H3K27ac, H3K27me3, H3K36me3, H3K79me2, and H4K20me, marks traditionally associated with enhancers, transcription initiation, transcriptional repression, and others. These patterns were observed in both normal and cancer cell lines. Additionally, we developed a novel computational method that identified 840 epigenetically regulated candidate genes and predicted transcription start-site switching and alternative exon splicing with up to 92% accuracy based on epigenetic patterning alone. Our results suggest that the epigenetic regulation of transcript isoform diversity may be a relatively common genome-wide phenomenon representing an avenue of deregulation in tumor development.

Traditionally, the regulation of gene expression was thought to be largely based on DNA and RNA sequence motifs. However, this dogma has recently been challenged as other factors, such as epigenetic patterning of the genome, have become better understood. Sparse but convincing experimental evidence suggests that the epigenetic background, in the form of histone modifications, acts as an additional layer of regulation determining how transcripts are processed. Here we developed a computational approach to investigate the genome-wide prevalence and the level of association between the enrichment of epigenetic marks and transcript diversity generated via alternative transcription start sites and splicing. We found that the role of epigenetic patterning in alternative transcription start-site switching is likely to be the same for all genes whereas the role of epigenetic patterns in splicing is likely gene-specific. Furthermore, we show that epigenetic data alone can be used to predict the inclusion pattern of an exon. These findings have significant implications for a better understanding of the regulation of transcript diversity in humans as well as the modifications arising during tumor development.

T cell activation regulates CD6 alternative splicing by transcription dynamics and SRSF1

The T cell-surface glycoprotein CD6 is a modulator of cellular responses and has been implicated in several autoimmune diseases such as multiple sclerosis, rheumatoid arthritis, and psoriasis. During Ag presentation, CD6 is targeted to the immunological synapse in a ligand binding-dependent manner, in which CD6 domain 3 directly contacts CD166, expressed on the APC. T cell activation results in the induction of CD6Δd3, an alternatively spliced isoform that lacks the ligand-binding domain and thus no longer localizes at the immunological synapse. In this study, we investigated the molecular mechanisms regulating the expression of CD6Δd3 upon human primary T cell activation. Using chromatin immunoprecipitation, we observed an increase in RNA polymerase II occupancy along the CD6 gene and augmented CD6 transcription. We showed that activation leads to transcription-related chromatin modifications, revealed by higher CD6 acetylation levels. Modulation of chromatin conformation using a histone deacetylase inhibitor that increases transcription rate causes an increase of exon 5 skipping. We further showed that the splicing factor SRSF1 binds to a regulatory element in CD6 intron 4, activating exon 5 splicing and promoting exon 5 inclusion. Concomitant with T cell activation-induced exon 5 skipping, we observed a downregulation of SRSF1. Using RNA immunoprecipitation, we showed that in activated T cells, SRSF1 recruitment to the CD6 transcript is impaired by increased chromatin acetylation levels. We propose that upon T cell activation, SRSF1 becomes limiting, and its function in CD6 exon 5 splicing is countered by an increase in CD6 transcription, dependent on chromatin acetylation.

Repressive histone methylation: a case study in deterministic versus stochastic gene regulation

Transcriptionally repressive histone lysine methylation is used by eukaryotes to tightly control cell fate. Here we explore the importance of this form of regulation in the control of clustered genes in the genome. Two distinctly regulated gene families with important roles in vertebrates are discussed, namely the Hox genes and olfactory receptor genes. Major recent advances in these two fields are compared and contrasted, with an emphasis on the roles of the two different forms of histone trimethylation. We discuss how this repression may impact both the transcriptional output of these loci and the way higher-order chromatin organization is related to their unique control.

Tyrosine phosphorylation of RNA polymerase II CTD is associated with antisense promoter transcription and active enhancers in mammalian cells

In mammals, the carboxy-terminal domain (CTD) of RNA polymerase (Pol) II consists of 52 conserved heptapeptide repeats containing the consensus sequence Tyr1-Ser2-Pro3-Thr4-Ser5-Pro6-Ser7. Post-translational modifications of the CTD coordinate the transcription cycle and various steps of mRNA maturation. Here we describe Tyr1 phosphorylation (Tyr1P) as a hallmark of promoter (5' associated) Pol II in mammalian cells, in contrast to what was described in yeast. Tyr1P is predominantly found in antisense orientation at promoters but is also specifically enriched at active enhancers. Mutation of Tyr1 to phenylalanine (Y1F) prevents the formation of the hyper-phosphorylated Pol IIO form, induces degradation of Pol II to the truncated Pol IIB form, and results in a lethal phenotype. Our results suggest that Tyr1P has evolved specialized and essential functions in higher eukaryotes associated with antisense promoter and enhancer transcription, and Pol II stability.DOI: http://dx.doi.org/10.7554/eLife.02105.001.

The epigenetic landscape of T-cell acute lymphoblastic leukemia

The genetic landscape of T-ALL has been very actively explored during the past decades. This leads to an overwhelming body of exciting novel findings providing insight into (1) the genetic heterogeneity of the disease with marked genetic subsets, (2) the mechanisms by which aberrant T-cell development drive leukemogenesis and (3) emerging opportunities for novel therapeutic interventions. Of further interest, recent genome wide sequencing studies identified proteins that actively participate in the regulation of the T-cell epigenome as novel oncogenes and tumor suppressor genes in T-ALL. The identification of these perturbed molecular epigenetic events in the pathogenesis of T-ALL will contribute to the further exploration of novel therapies in this cancer type. As some epigenetic therapies have recently been approved for a number of hematological neoplasms, one could speculate that targeted therapies against epigenetic regulators might offer good prospects for T-ALL treatment in the near future. In this review, we summarize the epigenetic discoveries made in T-ALL hitherto and discuss possible new venues for epigenetic therapeutic intervention in this aggressive subtype of human leukemia. This article is part of a Directed Issue entitled: Rare Cancers.

Reciprocal regulatory links between cotranscriptional splicing and chromatin

Here we review recent findings showing that chromatin organization adds another layer of complexity to the already intricate network of splicing regulatory mechanisms. Chromatin structure can impact splicing by either affecting the elongation rate of RNA polymerase II or by signaling the recruitment of splicing regulatory proteins. The C-terminal domain of the RNA polymerase II largest subunit serves as a coordination platform that binds factors required for adapting chromatin structure to both efficient transcription and processing of the newly synthesized RNA. Reciprocal interconnectivity of steps required for gene activation plays a critical role ensuring efficiency and fidelity of gene expression.

The alternative splicing side of cancer

Alternative splicing emerges as a potent and pervasive mechanism of gene expression regulation that expands the coding capacity of the genome and forms an intermediate layer of regulation between transcriptional and post-translational networks. Indeed, alternative splicing occupies a pivotal position in developmental programs and in the cell response to external and internal stimuli. Not surprisingly, therefore, its deregulation frequently leads to human disease. In this review we provide an updated overview of the impact of alternative splicing on tumorigenesis. Moreover, we discuss the intricacy of the reciprocal interactions between alternative splicing programs and signal transduction pathways, which appear to be crucially linked to cancer progression in response to the tumor microenvironment. Finally, we focus on the recently described interplay between splicing and chromatin organization which is expected to shed new lights into gene expression regulation in normal and cancer cells.

Sudemycin E influences alternative splicing and changes chromatin modifications

Sudemycin E is an analog of the pre-messenger RNA splicing modulator FR901464 and its derivative spliceostatin A. Sudemycin E causes the death of cancer cells through an unknown mechanism. We found that similar to spliceostatin A, sudemycin E binds to the U2 small nuclear ribonucleoprotein (snRNP) component SF3B1. Native chromatin immunoprecipitations showed that U2 snRNPs physically interact with nucleosomes. Sudemycin E induces a dissociation of the U2 snRNPs and decreases their interaction with nucleosomes. To determine the effect on gene expression, we performed genome-wide array analysis. Sudemycin E first causes a rapid change in alternative pre-messenger RNA splicing, which is later followed by changes in overall gene expression and arrest in the G2 phase of the cell cycle. The changes in alternative exon usage correlate with a loss of the H3K36me3 modification in chromatin encoding these exons. We propose that sudemycin E interferes with the ability of U2 snRNP to maintain an H3K36me3 modification in actively transcribed genes. Thus, in addition to the reversible changes in alternative splicing, sudemycin E causes changes in chromatin modifications that result in chromatin condensation, which is a likely contributing factor to cancer cell death.

Coupling mRNA processing with transcription in time and space

Maturation of mRNA precursors often occurs simultaneously with their synthesis by RNA polymerase II (Pol II). The co-transcriptional nature of mRNA processing has permitted the evolution of coupling mechanisms that coordinate transcription with mRNA capping, splicing, editing and 3' end formation. Recent experiments using sophisticated new methods for analysis of nascent RNA have provided important insights into the relative amount of co-transcriptional and post-transcriptional processing, the relationship between mRNA elongation and processing, and the role of the Pol II carboxy-terminal domain (CTD) in regulating these processes.

3' end formation of pre-mRNA and phosphorylation of Ser2 on the RNA polymerase II CTD are reciprocally coupled in human cells

Pre-mRNA 3′ end formation is coupled to transcription via the RNA Pol II C-terminal domain (CTD). However, how transcription and pre-mRNA maturation are coordinated in humans is poorly understood. Here, West and colleagues show that Pol II pausing promotes Ser2p by Cdk12, which recruits CPA factor CstF77 and is required for optimal 3′ end processing. This study delineates a reciprocal relationship between early steps in poly(A) site processing and Pol II Ser2p that ensures efficient pre-mRNA 3′ end formation in human cells.

3′ end formation of pre-mRNAs is coupled to their transcription via the C-terminal domain (CTD) of RNA polymerase II (Pol II). Nearly all protein-coding transcripts are matured by cleavage and polyadenylation (CPA), which is frequently misregulated in disease. Understanding how transcription is coordinated with CPA in human cells is therefore very important. We found that the CTD is heavily phosphorylated on Ser2 (Ser2p) at poly(A) (pA) signals coincident with recruitment of the CstF77 CPA factor. Depletion of the Ser2 kinase Cdk12 impairs Ser2p, CstF77 recruitment, and CPA, strongly suggesting that the processes are linked, as they are in budding yeast. Importantly, we additionally show that the high Ser2p signals at the 3′ end depend on pA signal function. Down-regulation of CPA results in the loss of a 3′ Ser2p peak, whereas a new peak is formed when CPA is induced de novo. Finally, high Ser2p signals are generated by Pol II pausing, which is a well-known feature of pA site recognition. Thus, a reciprocal relationship between early steps in pA site processing and Ser2p ensures efficient 3′ end formation.

Introduction to cotranscriptional RNA splicing

The discovery that many intron-containing genes can be cotranscriptionally spliced has led to an increased understanding of how splicing and transcription are intricately intertwined. Cotranscriptional splicing has been demonstrated in a number of different organisms and has been shown to play roles in coordinating both constitutive and alternative splicing. The nature of cotranscriptional splicing suggests that changes in transcription can dramatically affect splicing, and new evidence suggests that splicing can, in turn, influence transcription. In this chapter, we discuss the mechanisms and consequences of cotranscriptional splicing and introduce some of the tools used to measure this process.

Chromatin and splicing

In the past several years, the relationship between chromatin structure and mRNA processing has been the source of significant investigation across diverse disciplines. Central to these efforts was an unanticipated nonrandom distribution of chromatin marks across transcribed regions of protein-coding genes. In addition to the presence of specific histone modifications at the 5' and 3' ends of genes, exonic DNA was demonstrated to present a distinct chromatin landscape relative to intronic DNA. As splicing in higher eukaryotes predominantly occurs co-transcriptionally, these studies raised the possibility that chromatin modifications may aid the spliceosome in the detection of exons amidst vast stretches of noncoding intronic sequences. Recent investigations have supported a direct role for chromatin in splicing regulation and have suggested an intriguing role for splicing in the establishment of chromatin modifications. Here we will summarize an accumulating body of data that begins to reveal extensive coupling between chromatin structure and pre-mRNA splicing.

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.

Variation in chromatin accessibility in human kidney cancer links H3K36 methyltransferase loss with widespread RNA processing defects

Comprehensive sequencing of human cancers has identified recurrent mutations in genes encoding chromatin regulatory proteins. For clear cell renal cell carcinoma (ccRCC), three of the five commonly mutated genes encode the chromatin regulators PBRM1, SETD2, and BAP1. How these mutations alter the chromatin landscape and transcriptional program in ccRCC or other cancers is not understood. Here, we identified alterations in chromatin organization and transcript profiles associated with mutations in chromatin regulators in a large cohort of primary human kidney tumors. By associating variation in chromatin organization with mutations in SETD2, which encodes the enzyme responsible for H3K36 trimethylation, we found that changes in chromatin accessibility occurred primarily within actively transcribed genes. This increase in chromatin accessibility was linked with widespread alterations in RNA processing, including intron retention and aberrant splicing, affecting ∼25% of all expressed genes. Furthermore, decreased nucleosome occupancy proximal to misspliced exons was observed in tumors lacking H3K36me3. These results directly link mutations in SETD2 to chromatin accessibility changes and RNA processing defects in cancer. Detecting the functional consequences of specific mutations in chromatin regulatory proteins in primary human samples could ultimately inform the therapeutic application of an emerging class of chromatin-targeted compounds.

Epigenetic and transcriptional signatures of stable versus plastic differentiation of proinflammatory γδ T cell subsets

Two distinct subsets of γδ T cells that produce interleukin 17 (IL-17) (CD27(-) γδ T cells) or interferon-γ (IFN-γ) (CD27(+) γδ T cells) develop in the mouse thymus, but the molecular determinants of their functional potential in the periphery remain unknown. Here we conducted a genome-wide characterization of the methylation patterns of histone H3, along with analysis of mRNA encoding transcription factors, to identify the regulatory networks of peripheral IFN-γ-producing or IL-17-producing γδ T cell subsets in vivo. We found that CD27(+) γδ T cells were committed to the expression of Ifng but not Il17, whereas CD27(-) γδ T cells displayed permissive chromatin configurations at loci encoding both cytokines and their regulatory transcription factors and differentiated into cells that produced both IL-17 and IFN-γ in a tumor microenvironment.

Regulation of alternative splicing by local histone modifications: potential roles for RNA-guided mechanisms

The molecular mechanisms through which alternative splicing and histone modifications regulate gene expression are now understood in considerable detail. Here, we discuss recent studies that connect these two previously separate avenues of investigation, beginning with the unexpected discoveries that nucleosomes are preferentially positioned over exons and DNA methylation and certain histone modifications also show exonic enrichment. These findings have profound implications linking chromatin structure, histone modification and splicing regulation. Complementary single gene studies provided insight into the mechanisms through which DNA methylation and histones modifications modulate alternative splicing patterns. Here, we review an emerging theme resulting from these studies: RNA-guided mechanisms integrating chromatin modification and splicing. Several groundbreaking papers reported that small noncoding RNAs affect alternative exon usage by targeting histone methyltransferase complexes to form localized facultative heterochromatin. More recent studies provided evidence that pre-messenger RNA itself can serve as a guide to enable precise alternative splicing regulation via local recruitment of histone-modifying enzymes, and emerging evidence points to a similar role for long noncoding RNAs. An exciting challenge for the future is to understand the impact of local modulation of transcription elongation rates on the dynamic interplay between histone modifications, alternative splicing and other processes occurring on chromatin.

The C-terminal domain of Brd2 is important for chromatin interaction and regulation of transcription and alternative splicing

This study determines genes that are regulated by Brd2 and finds that, in addition to expression control, Brd2 modulates the alternative splicing of several hundred genes. The in vivo interaction of Brd2 with chromatin is analyzed, and the contributions of individual Brd2 domains to the chromatin interaction are determined.

Brd2 is a member of the bromodomain extra terminal (BET) protein family, which consists of four chromatin-interacting proteins that regulate gene expression. Each BET protein contains two N-terminal bromodomains, which recognize acetylated histones, and the C-terminal protein–protein interaction domain. Using a genome-wide screen, we identify 1450 genes whose transcription is regulated by Brd2. In addition, almost 290 genes change their alternative splicing pattern upon Brd2 depletion. Brd2 is specifically localized at promoters of target genes, and our data show that Brd2 interaction with chromatin cannot be explained solely by histone acetylation. Using coimmunoprecipitation and live-cell imaging, we show that the C-terminal part is crucial for Brd2 association with chromatin. Live-cell microscopy also allows us to map the average binding time of Brd2 to chromatin and quantify the contributions of individual Brd2 domains to the interaction with chromatin. Finally, we show that bromodomains and the C-terminal domain are equally important for transcription and splicing regulation, which correlates with the role of these domains in Brd2 binding to chromatin.

lincRNAs: genomics, evolution, and mechanisms

Long intervening noncoding RNAs (lincRNAs) are transcribed from thousands of loci in mammalian genomes and might play widespread roles in gene regulation and other cellular processes. This Review outlines the emerging understanding of lincRNAs in vertebrate animals, with emphases on how they are being identified and current conclusions and questions regarding their genomics, evolution and mechanisms of action.

Chromatin's thread to alternative splicing regulation

Intron removal (pre-mRNA splicing) is a necessary step for expression of most genes in higher eukaryotes. Alternative splice site selection is a prevalent mechanism that diversifies genome outputs and offers ample opportunities for gene regulation in these organisms. Pre-mRNA splicing occurs co-transcriptionally and is influenced by features in chromatin structure, including nucleosome density and epigenetic modifications. We review here the molecular mechanisms by which the reciprocal interplay between chromatin and RNA processing can contribute to alternative splicing regulation.

Ancient cis-regulatory constraints and the evolution of genome architecture

The order of genes along metazoan chromosomes has generally been thought to be largely random, with few implications for organismal function. However, two recent studies, reporting hundreds of pairs of genes that have remained linked in diverse metazoan species over hundreds of millions of years of evolution, suggest widespread functional implications for gene order. These associations appear to largely reflect cis-regulatory constraints, with either (i) multiple genes sharing transcriptional regulatory elements, or (ii) regulatory elements for a developmental gene being found within a neighboring 'bystander' gene (known as a genomic regulatory block). We discuss implications, questions raised, and new research directions arising from these studies, as well as evidence for similar phenomena in other eukaryotic groups.

Histone modifications for human epigenome analysis

Histones function both positively and negatively in the regulation of gene expression, mainly governed by post-translational modifications on specific amino acid residues. Although histone modifications are not necessarily prerequisite codes, they may still serve as good epigenetic indicators of chromatin state associated with gene activation or repression. In particular, six emerging classes of histone H3 modifications are subjected for epigenome profiling by the International Human Epigenome Consortium. In general, transcription start sites of actively transcribed genes are marked by trimethylated H3K4 (H3K4me3) and acetylated H3K27 (H3K27ac), and active enhancers can be identified by enrichments of both monomethylated H3K4 (H3K4me1) and H3K27ac. Gene bodies of actively transcribed genes are associated with trimethylated H3K36 (H3K36me3). Gene repression can be mediated through two distinct mechanisms involving trimethylated H3K9 (H3K9me3) and trimethylated H3K27 (H3K27me3). Enrichments of these histone modifications on specific loci, or in genome wide, in given cells can be analyzed by chromatin immunoprecipitation (ChIP)-based methods using an antibody directed against the site-specific modification. When performing ChIP experiments, one should be careful about the specificity of antibody, as this affects the data interpretation. If cell samples with preserved histone-DNA contacts are available, evaluation of histone modifications, in addition to DNA methylaion, at specific gene loci would be useful for deciphering the epigenome state for human genetics studies.