Genomic mapping of RNA polymerase II reveals sites of co-transcriptional regulation in human cells

The Functional Relationship Between RNA Splicing and the Chromatin Landscape

Chromatin is a highly regulated and dynamic structure that has been shown to play an essential role in transcriptional and co-transcriptional regulation. In the context of RNA splicing, early evidence suggested a loose connection between the chromatin landscape and splicing. More recently, it has been shown that splicing occurs in a co-transcriptional manner, meaning that the splicing process occurs in the context of chromatin. Experimental and computational evidence have also shown that chromatin dynamics can influence the splicing process and vice versa. However, much of this evidence provides mainly correlative relationships between chromatin and splicing with just a few concrete examples providing defined molecular mechanisms by which these two processes are functionally related. Nevertheless, it is clear that chromatin and RNA splicing are tightly interconnected to one another. In this review, we highlight the current state of knowledge of the relationship between chromatin and splicing.

Histone Marks-Dependent Effect on Alternative Splicing: New Perspectives for Targeted Splicing Modulation in Cancer?

Alternative splicing (AS) is a tightly regulated mechanism that generates the complex human proteome from a small number of genes. Cis-regulatory RNA motifs in exons and introns control AS, recruiting positive and negative trans-acting splicing regulators. At a higher level, chromatin affects splicing events. Growing evidence indicates that the popular histone code hypothesis can be extended to RNA-level processes, such as AS. In addition to nucleosome positioning, which can generate transcriptional barriers to shape the final splicing outcome, histone post-translational modifications can contribute to the detailed regulation of single exon inclusion/exclusion. A histone-based system can identify alternatively spliced chromatin stretches, affecting RNAPII elongation locally or recruiting splicing components via adaptor complexes. In tumor cells, several mechanisms trigger misregulated AS events and produce cancer-associated transcripts. On a genome-wide level, aberrant AS can be the consequence of dysfunctional epigenetic splicing code, including altered enrichment in histone post-translational modifications. This review describes the main findings related to the effect of histone modifications and variants on splicing outcome and how a dysfunctional epigenetic splicing code triggers aberrant AS in cancer. In addition, it highlights recent advances in programmable DNA-targeting technologies and their possible application for AS targeted epigenetic modulation.

Impacts and mechanisms of alternative mRNA splicing in cancer metabolism, immune response, and therapeutics

Alternative pre-mRNA splicing (AS) provides the potential to produce diversity at RNA and protein levels. Disruptions in the regulation of pre-mRNA splicing can lead to diseases. With the development of transcriptome and genome sequencing technology, increasing diseases have been identified to be associated with abnormal splicing of mRNAs. In tumors, abnormal alternative splicing frequently plays critical roles in cancer pathogenesis and may be considered as new biomarkers and therapeutic targets for cancer intervention. Metabolic abnormalities and immune disorders are important hallmarks of cancer. AS produces multiple different isoforms and diversifies protein expression, which is utilized by the immune and metabolic reprogramming systems to expand gene functions. The abnormal splicing events contributed to tumor progression, partially due to effects on immune response and metabolic reprogramming. Herein, we reviewed the vital role of alternative splicing in regulating cancer metabolism and immune response. We discussed how alternative splicing regulates metabolic reprogramming of cancer cells and antitumor immune response, and the possible strategies to targeting alternative splicing pathways or splicing-regulated metabolic pathway in the context of anticancer immunotherapy. Further, we highlighted the challenges and discuss the perspectives for RNA-based strategies for the treatment of cancer with abnormally alternative splicing isoforms.

The long noncoding RNA glycoLINC assembles a lower glycolytic metabolon to promote glycolysis

Non-covalent complexes of glycolytic enzymes, called metabolons, were postulated in the 1970s, but the concept has been controversial. Here we show that a c-Myc-responsive long noncoding RNA (lncRNA) that we call glycoLINC (gLINC) acts as a backbone for metabolon formation between all four glycolytic payoff phase enzymes (PGK1, PGAM1, ENO1, and PKM2) along with lactate dehydrogenase A (LDHA). The gLINC metabolon enhances glycolytic flux, increases ATP production, and enables cell survival under serine deprivation. Furthermore, gLINC overexpression in cancer cells promotes xenograft growth in mice fed a diet deprived of serine, suggesting that cancer cells employ gLINC during metabolic reprogramming. We propose that gLINC makes a functional contribution to cancer cell adaptation and provide the first example of a lncRNA-facilitated metabolon.

RNA polymerase II speed: a key player in controlling and adapting transcriptome composition

RNA polymerase II (RNA Pol II) speed or elongation rate, i.e., the number of nucleotides synthesized per unit of time, is a major determinant of transcriptome composition. It controls co‐transcriptional processes such as splicing, polyadenylation, and transcription termination, thus regulating the production of alternative splice variants, circular RNAs, alternatively polyadenylated transcripts, or read‐through transcripts. RNA Pol II speed itself is regulated in response to intra‐ and extra‐cellular stimuli and can in turn affect the transcriptome composition in response to these stimuli. Evidence points to a potentially important role of transcriptome composition modification through RNA Pol II speed regulation for adaptation of cells to a changing environment, thus pointing to a function of RNA Pol II speed regulation in cellular physiology. Analyzing RNA Pol II speed dynamics may therefore be central to fully understand the regulation of physiological processes, such as the development of multicellular organisms. Recent findings also raise the possibility that RNA Pol II speed deregulation can be detrimental and participate in disease progression. Here, we review initial and current approaches to measure RNA Pol II speed, as well as providing an overview of the factors controlling speed and the co‐transcriptional processes which are affected. Finally, we discuss the role of RNA Pol II speed regulation in cell physiology.

Genome-wide profiling reveals functional interplay of DNA sequence composition, transcriptional activity, and nucleosome positioning in driving DNA supercoiling and helix destabilization in C. elegans

DNA topology and alternative DNA structures are implicated in regulating diverse biological processes. Although biomechanical properties of these structures have been studied extensively in vitro, characterization in vivo, particularly in multicellular organisms, is limited. We devised new methods to map DNA supercoiling and single-stranded DNA in Caenorhabditis elegans embryos and diapause larvae. To map supercoiling, we quantified the incorporation of biotinylated psoralen into DNA using high-throughput sequencing. To map single-stranded DNA, we combined permanganate treatment with genome-wide sequencing of induced double-stranded breaks. We found high levels of negative supercoiling at transcription start sites (TSSs) in embryos. GC-rich regions flanked by a sharp GC-to-AT transition delineate boundaries of supercoil propagation. In contrast to TSSs in embryos, TSSs in diapause larvae showed dramatic reductions in negative supercoiling without concomitant attenuation of transcription, suggesting developmental-stage-specific regulation. To assess whether alternative DNA structures control chromosome architecture and gene expression, we examined DNA supercoiling in the context of X-Chromosome dosage compensation. We showed that the condensin dosage compensation complex creates negative supercoils locally at its highest-occupancy binding sites but found no evidence for large-scale supercoiling domains along X Chromosomes. In contrast to transcription-coupled negative supercoiling, single-strandedness, which is most pronounced at transcript end sites, is dependent on high AT content and symmetrically positioned nucleosomes. We propose that sharp transitions in sequence composition at functional genomic elements constitute a common regulatory code and that DNA structure and propagation of torsional stress at regulatory elements are critical parameters in shaping important developmental events.

Splicing-associated chromatin signatures: a combinatorial and position-dependent role for histone marks in splicing definition

Alternative splicing relies on the combinatorial recruitment of splicing regulators to specific RNA binding sites. Chromatin has been shown to impact this recruitment. However, a limited number of histone marks have been studied at a global level. In this work, a machine learning approach, applied to extensive epigenomics datasets in human H1 embryonic stem cells and IMR90 foetal fibroblasts, has identified eleven chromatin modifications that differentially mark alternatively spliced exons depending on the level of exon inclusion. These marks act in a combinatorial and position-dependent way, creating characteristic splicing-associated chromatin signatures (SACS). In support of a functional role for SACS in coordinating splicing regulation, changes in the alternative splicing of SACS-marked exons between ten different cell lines correlate with changes in SACS enrichment levels and recruitment of the splicing regulators predicted by RNA motif search analysis. We propose the dynamic nature of chromatin modifications as a mechanism to rapidly fine-tune alternative splicing when necessary.

Chromatin is known to regulate splicing by modulating recruitment of splicing factors. Using machine learning approaches, the authors have underlined a chromatin code for alternative splicing regulation that is conserved amongst cell lines.

Transcription and splicing: A two-way street

RNA synthesis by RNA polymerase II and RNA processing are closely coupled during the transcription cycle of protein-coding genes. This coupling affords opportunities for quality control and regulation of gene expression and the effects can go in both directions. For example, polymerase speed can affect splice site selection and splicing can increase transcription and affect the chromatin landscape. Here we review the many ways that transcription and splicing influence one another, including how splicing "talks back" to transcription. We will also place the connections between transcription and splicing in the context of other RNA processing events that define the exons that will make up the final mRNA. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Promoter-proximal pausing of RNA polymerase II: a nexus of gene regulation

In this review, Core et al. discuss the recent advances in our understanding of the early steps in Pol II transcription, highlighting the events and factors involved in the establishment and release of paused Pol II. They also discuss a number of unanswered questions about the regulation and function of Pol II pausing.

Precise spatio–temporal control of gene activity is essential for organismal development, growth, and survival in a changing environment. Decisive steps in gene regulation involve the pausing of RNA polymerase II (Pol II) in early elongation, and the controlled release of paused polymerase into productive RNA synthesis. Here we describe the factors that enable pausing and the events that trigger Pol II release into the gene. We also discuss open questions in the field concerning the stability of paused Pol II, nucleosomes as obstacles to elongation, and potential roles of pausing in defining the precision and dynamics of gene expression.

Epigenetic regulation of transcriptional plasticity associated with developmental song learning

Ethologists discovered over 100 years ago that some lifelong behavioural patterns were acquired exclusively during restricted developmental phases called critical periods (CPs). Developmental song learning in zebra finches is one of the most striking examples of a CP for complex learned behaviour. After post-hatch day 65, whether or not a juvenile male can memorize the song of a 'tutor' depends on his experiences in the month prior. If he experienced a tutor, he can no longer learn, but if he has been isolated from hearing a tutor the learning period is extended. We aimed to identify how tutor experience alters the brain and controls the ability to learn. Epigenetic landscapes are modulated by experience and are able to regulate the transcription of sets of genes, thereby affecting cellular function. Thus, we hypothesized that tutor experiences determine the epigenetic landscape in the auditory forebrain, a region required for tutor song memorization. Using ChIPseq, RNAseq and molecular biology, we provide evidence that naturalistic experiences associated with the ability to learn can induce epigenetic changes, and propose transcriptional plasticity as a mediator of CP learning potential.

Coupling of PARP1-mediated chromatin structural changes to transcriptional RNA polymerase II elongation and cotranscriptional splicing

Background

Recently, we showed that PARP1 is involved in cotranscriptional splicing, possibly by bridging chromatin to RNA and recruiting splicing factors. It also can influence alternative splicing decisions through the regulation of RNAPII elongation. In this study, we investigated the effect of PARP1-mediated chromatin changes on RNAPII movement, during transcription and alternative splicing.

Results

We show that RNAPII pauses at PARP1–chromatin structures within the gene body. Knockdown of PARP1 abolishes this RNAPII pausing, suggesting that PARP1 may regulate RNAPII elongation. Additionally, PARP1 alters nucleosome deposition and histone post-translational modifications at specific exon–intron boundaries, thereby affecting RNAPII movement. Lastly, genome-wide analyses confirmed that PARP1 influences changes in RNAPII elongation by either reducing or increasing the rate of RNAPII elongation depending on the chromatin context.

Conclusions

These studies suggest a context-specific effect of PARP1–chromatin binding on RNA polymerase movement and provide a platform to delineate PARP1’s role in RNA biogenesis and processing.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0261-1) contains supplementary material, which is available to authorized users.

Emerging Properties and Functional Consequences of Noncoding Transcription

Eukaryotic genomes are rich in transcription units encoding "long noncoding RNAs" (lncRNAs). The purpose of all this transcription is unclear since most lncRNAs are quickly targeted for destruction during synthesis or shortly thereafter. As debates continue over the functional significance of many specific lncRNAs, support grows for the notion that the act of transcription rather than the RNA product itself is functionally important in many cases. Indeed, this alternative mechanism might better explain how low-abundance lncRNAs transcribed from noncoding DNA function in organisms. Here, we highlight some of the recently emerging features that distinguish coding from noncoding transcription and discuss how these differences might have important implications for the functional consequences of noncoding transcription.

Microprocessor Recruitment to Elongating RNA Polymerase II Is Required for Differential Expression of MicroRNAs

The cellular abundance of mature microRNAs (miRNAs) is dictated by the efficiency of nuclear processing of primary miRNA transcripts (pri-miRNAs) into pre-miRNA intermediates. The Microprocessor complex of Drosha and DGCR8 carries this out, but it has been unclear what controls Microprocessor's differential processing of various pri-miRNAs. Here, we show that Drosophila DGCR8 (Pasha) directly associates with the C-terminal domain of the RNA polymerase II elongation complex when it is phosphorylated by the Cdk9 kinase (pTEFb). When association is blocked by loss of Cdk9 activity, a global change in pri-miRNA processing is detected. Processing of pri-miRNAs with a UGU sequence motif in their apical junction domain increases, while processing of pri-miRNAs lacking this motif decreases. Therefore, phosphorylation of RNA polymerase II recruits Microprocessor for co-transcriptional processing of non-UGU pri-miRNAs that would otherwise be poorly processed. In contrast, UGU-positive pri-miRNAs are robustly processed by Microprocessor independent of RNA polymerase association.

Transcription-coupled changes to chromatin underpin gene silencing by transcriptional interference

Long non-coding RNA (lncRNA) transcription into a downstream promoter frequently results in transcriptional interference. However, the mechanism of this repression is not fully understood. We recently showed that drug tolerance in fission yeast Schizosaccharomyces pombe is controlled by lncRNA transcription upstream of the tgp1⁺ permease gene. Here we demonstrate that transcriptional interference of tgp1⁺ involves several transcription-coupled chromatin changes mediated by conserved elongation factors Set2, Clr6CII, Spt6 and FACT. These factors are known to travel with RNAPII and establish repressive chromatin in order to limit aberrant transcription initiation from cryptic promoters present in gene bodies. We therefore conclude that conserved RNAPII-associated mechanisms exist to both suppress intragenic cryptic promoters during genic transcription and to repress gene promoters by transcriptional interference. Our analyses also demonstrate that key mechanistic features of transcriptional interference are shared between S. pombe and the highly divergent budding yeast Saccharomyces cerevisiae. Thus, transcriptional interference is an ancient, conserved mechanism for tightly controlling gene expression. Our mechanistic insights allowed us to predict and validate a second example of transcriptional interference involving the S. pombe pho1⁺ gene. Given that eukaryotic genomes are pervasively transcribed, transcriptional interference likely represents a more general feature of gene regulation than is currently appreciated.

How Are Short Exons Flanked by Long Introns Defined and Committed to Splicing?

The splice sites (SSs) delimiting an intron are brought together in the earliest step of spliceosome assembly yet it remains obscure how SS pairing occurs, especially when introns are thousands of nucleotides long. Splicing occurs in vivo in mammals within minutes regardless of intron length, implying that SS pairing can instantly follow transcription. Also, factors required for SS pairing, such as the U1 small nuclear ribonucleoprotein (snRNP) and U2AF65, associate with RNA polymerase II (RNAPII), while nucleosomes preferentially bind exonic sequences and associate with U2 snRNP. Based on recent publications, we assume that the 5' SS-bound U1 snRNP can remain tethered to RNAPII until complete synthesis of the downstream intron and exon. An additional U1 snRNP then binds the downstream 5' SS, whereas the RNAPII-associated U2AF65 binds the upstream 3' SS to facilitate SS pairing along with exon definition. Next, the nucleosome-associated U2 snRNP binds the branch site to advance splicing complex assembly. This may explain how RNAPII and chromatin are involved in spliceosome assembly and how introns lengthened during evolution with a relatively minimal compromise in splicing.

Global intron retention mediated gene regulation during CD4+ T cell activation

T cell activation is a well-established model for studying cellular responses to exogenous stimulation. Using strand-specific RNA-seq, we observed that intron retention is prevalent in polyadenylated transcripts in resting CD4⁺ T cells and is significantly reduced upon T cell activation. Several lines of evidence suggest that intron-retained transcripts are less stable than fully spliced transcripts. Strikingly, the decrease in intron retention (IR) levels correlate with the increase in steady-state mRNA levels. Further, the majority of the genes upregulated in activated T cells are accompanied by a significant reduction in IR. Of these 1583 genes, 185 genes are predominantly regulated at the IR level, and highly enriched in the proteasome pathway, which is essential for proper T cell proliferation and cytokine release. These observations were corroborated in both human and mouse CD4⁺ T cells. Our study revealed a novel post-transcriptional regulatory mechanism that may potentially contribute to coordinated and/or quick cellular responses to extracellular stimuli such as an acute infection.

Chromatin structure and pre-mRNA processing work together

Chromatin is the natural context for transcription elongation. However, the elongating RNA polymerase II (RNAPII) is forced to pause by the positioned nucleosomes present in gene bodies. Here, we briefly discuss the current results suggesting that those pauses could serve as a mechanism to coordinate transcription elongation with pre-mRNA processing. Further, histone post-translational modifications have been found to regulate the recruitment of factors involved in pre-mRNA processing. This view highlights the important regulatory role of the chromatin context in the whole process of the mature mRNA synthesis.

Localization of RNAPII and 3' end formation factor CstF subunits on C. elegans genes and operons

Transcription termination is mechanistically coupled to pre-mRNA 3' end formation to prevent transcription much beyond the gene 3' end. C. elegans, however, engages in polycistronic transcription of operons in which 3' end formation between genes is not accompanied by termination. We have performed RNA polymerase II (RNAPII) and CstF ChIP-seq experiments to investigate at a genome-wide level how RNAPII can transcribe through multiple poly-A signals without causing termination. Our data shows that transcription proceeds in some ways as if operons were composed of multiple adjacent single genes. Total RNAPII shows a small peak at the promoter of the gene cluster and a much larger peak at 3' ends. These 3' peaks coincide with maximal phosphorylation of Ser2 within the C-terminal domain (CTD) of RNAPII and maximal localization of the 3' end formation factor CstF. This pattern occurs at all 3' ends including those at internal sites in operons where termination does not occur. Thus the normal mechanism of 3' end formation does not always result in transcription termination. Furthermore, reduction of CstF50 by RNAi did not substantially alter the pattern of CstF64, total RNAPII, or Ser2 phosphorylation at either internal or terminal 3' ends. However, CstF50 RNAi did result in a subtle reduction of CstF64 binding upstream of the site of 3' cleavage, suggesting that the CstF50/CTD interaction may facilitate bringing the 3' end machinery to the transcription complex.

Architectural and Functional Commonalities between Enhancers and Promoters

With the explosion of genome-wide studies of regulated transcription, it has become clear that traditional definitions of enhancers and promoters need to be revisited. These control elements can now be characterized in terms of their local and regional architecture, their regulatory components, including histone modifications and associated binding factors, and their functional contribution to transcription. This Review discusses unifying themes between promoters and enhancers in transcriptional regulatory mechanisms.

Defective histone supply causes changes in RNA polymerase II elongation rate and cotranscriptional pre-mRNA splicing

RNA polymerase II (RNAPII) transcription elongation is a highly regulated process that greatly influences mRNA levels as well as pre-mRNA splicing. Despite many studies in vitro, how chromatin modulates RNAPII elongation in vivo is still unclear. Here, we show that a decrease in the level of available canonical histones leads to more accessible chromatin with decreased levels of canonical histones and variants H2A.X and H2A.Z and increased levels of H3.3. With this altered chromatin structure, the RNAPII elongation rate increases, and the kinetics of pre-mRNA splicing is delayed with respect to RNAPII elongation. Consistent with the kinetic model of cotranscriptional splicing, the rapid RNAPII elongation induced by histone depletion promotes the skipping of variable exons in the CD44 gene. Indeed, a slowly elongating mutant of RNAPII was able to rescue this defect, indicating that the defective splicing induced by histone depletion is a direct consequence of the increased elongation rate. In addition, genome-wide analysis evidenced that histone reduction promotes widespread alterations in pre-mRNA processing, including intron retention and changes in alternative splicing. Our data demonstrate that pre-mRNA splicing may be regulated by chromatin structure through the modulation of the RNAPII elongation rate.

The Hox cluster microRNA miR-615: a case study of intronic microRNA evolution

BACKGROUND

Introns represent a potentially rich source of existing transcription for the evolution of novel microRNAs (miRNAs). Within the Hox gene clusters, a miRNA gene, miR-615, is located within the intron of the Hoxc5 gene. This miRNA has a restricted phylogenetic distribution, providing an opportunity to examine the origin and evolution of a new miRNA within the intron of a developmentally-important homeobox gene.

RESULTS

Alignment and structural analyses show that the sequence is highly conserved across eutherian mammals and absent in non-mammalian tetrapods. Marsupials possess a similar sequence which we predict will not be efficiently processed as a miRNA. Our analyses suggest that transcription of HOXC5 in humans is accompanied by expression of miR-615 in all cases, but that the miRNA can also be transcribed independently of its host gene through the use of an intragenic promoter. We present scenarios for the evolution of miR-615 through intronic exaptation, and speculate on the acquisition of independent transcriptional regulation. Target prediction and transcriptomic analyses suggest that the dominant product of miR-615 is involved in the regulation of growth and a range of developmental processes.

CONCLUSIONS

The miR-615 gene evolved within the intron of Hoxc5 in the ancestor of placental mammals. Using miR-615 as a case study, we propose a model by which a functional miRNA can emerge within an intron gradually, by selection on secondary structure followed by evolution of an independent miRNA promoter. The location within a Hox gene intron is of particular interest as the miRNA is specific to placental mammals, is co-expressed with its host gene and may share complementary functions.

An Overview of Genome Organization and How We Got There: from FISH to Hi-C

In humans, nearly two meters of genomic material must be folded to fit inside each micrometer-scale cell nucleus while remaining accessible for gene transcription, DNA replication, and DNA repair. This fact highlights the need for mechanisms governing genome organization during any activity and to maintain the physical organization of chromosomes at all times. Insight into the functions and three-dimensional structures of genomes comes mostly from the application of visual techniques such as fluorescence in situ hybridization (FISH) and molecular approaches including chromosome conformation capture (3C) technologies. Recent developments in both types of approaches now offer the possibility of exploring the folded state of an entire genome and maybe even the identification of how complex molecular machines govern its shape. In this review, we present key methodologies used to study genome organization and discuss what they reveal about chromosome conformation as it relates to transcription regulation across genomic scales in mammals.

Native elongating transcript sequencing reveals human transcriptional activity at nucleotide resolution

Major features of transcription by human RNA polymerase II (Pol II) remain poorly defined due to a lack of quantitative approaches for visualizing Pol II progress at nucleotide resolution. We developed a simple and powerful approach for performing native elongating transcript sequencing (NET-seq) in human cells that globally maps strand-specific Pol II density at nucleotide resolution. NET-seq exposes a mode of antisense transcription that originates downstream and converges on transcription from the canonical promoter. Convergent transcription is associated with a distinctive chromatin configuration and is characteristic of lower-expressed genes. Integration of NET-seq with genomic footprinting data reveals stereotypic Pol II pausing coincident with transcription factor occupancy. Finally, exons retained in mature transcripts display Pol II pausing signatures that differ markedly from skipped exons, indicating an intrinsic capacity for Pol II to recognize exons with different processing fates. Together, human NET-seq exposes the topography and regulatory complexity of human gene expression.

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.

The many faces of Dicer: the complexity of the mechanisms regulating Dicer gene expression and enzyme activities

There is increasing evidence indicating that the production of small regulatory RNAs is not the only process in which ribonuclease Dicer can participate. For example, it has been demonstrated that this enzyme is also involved in chromatin structure remodelling, inflammation and apoptotic DNA degradation. Moreover, it has become increasingly clear that cellular transcript and protein levels of Dicer must be strictly controlled because even small changes in their accumulation can initiate various pathological processes, including carcinogenesis. Accordingly, in recent years, a number of studies have been performed to identify the factors regulating Dicer gene expression and protein activity. As a result, a large amount of complex and often contradictory data has been generated. None of these data have been subjected to an exhaustive review or critical discussion. This review attempts to fill this gap by summarizing the current knowledge of factors that regulate Dicer gene transcription, primary transcript processing, mRNA translation and enzyme activity. Because of the high complexity of this topic, this review mainly concentrates on human Dicer. This review also focuses on an additional regulatory layer of Dicer activity involving the interactions of protein and RNA factors with Dicer substrates.

Genome-wide mapping of matrix attachment regions in Drosophila melanogaster

Background

Eukaryotic genome acquires functionality upon proper packaging within the nucleus. This process is facilitated by the structural framework of Nuclear Matrix, a nucleo-proteinaceous meshwork. Matrix Attachment Regions (MARs) in the genome serve as anchoring sites to this framework.

Results

Here we report direct sequencing of the MAR preparation from Drosophila melanogaster embryos and identify >7350 MARs. This amounts to ~2.5% of the fly genome and often coincide with AT rich non-coding regions. We find significant association of MARs with the origins of replication, transcription start sites, paused RNA Polymerase II sites and exons, but not introns, of highly expressed genes. We also identified sequence motifs and repeats that constitute MARs.

Conclusion

Our data reveal the contact points of genome to the nuclear architecture and provide a link between nuclear functions and genomic packaging.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1022) contains supplementary material, which is available to authorized users.

Map of open and closed chromatin domains in Drosophila genome

BACKGROUND

Chromatin compactness has been considered a major determinant of gene activity and has been associated with specific chromatin modifications in studies on a few individual genetic loci. At the same time, genome-wide patterns of open and closed chromatin have been understudied, and are at present largely predicted from chromatin modification and gene expression data. However the universal applicability of such predictions is not self-evident, and requires experimental verification.

RESULTS

We developed and implemented a high-throughput analysis for general chromatin sensitivity to DNase I which provides a comprehensive epigenomic assessment in a single assay. Contiguous domains of open and closed chromatin were identified by computational analysis of the data, and correlated to other genome annotations including predicted chromatin "states", individual chromatin modifications, nuclear lamina interactions, and gene expression. While showing that the widely trusted predictions of chromatin structure are correct in the majority of cases, we detected diverse "exceptions" from the conventional rules. We found a profound paucity of chromatin modifications in a major fraction of closed chromatin, and identified a number of loci where chromatin configuration is opposite to that expected from modification and gene expression patterns. Further, we observed that chromatin of large introns tends to be closed even when the genes are expressed, and that a significant proportion of active genes including their promoters are located in closed chromatin.

CONCLUSIONS

These findings reveal limitations of the existing predictive models, indicate novel mechanisms of epigenetic regulation, and provide important insights into genome organization and function.

Genome-wide dynamics of Pol II elongation and its interplay with promoter proximal pausing, chromatin, and exons

Production of mRNA depends critically on the rate of RNA polymerase II (Pol II) elongation. To dissect Pol II dynamics in mouse ES cells, we inhibited Pol II transcription at either initiation or promoter-proximal pause escape with Triptolide or Flavopiridol, and tracked Pol II kinetically using GRO-seq. Both inhibitors block transcription of more than 95% of genes, showing that pause escape, like initiation, is a ubiquitous and crucial step within the transcription cycle. Moreover, paused Pol II is relatively stable, as evidenced from half-life measurements at ∼3200 genes. Finally, tracking the progression of Pol II after drug treatment establishes Pol II elongation rates at over 1000 genes. Notably, Pol II accelerates dramatically while transcribing through genes, but slows at exons. Furthermore, intergenic variance in elongation rates is substantial, and is influenced by a positive effect of H3K79me2 and negative effects of exon density and CG content within genes.DOI: http://dx.doi.org/10.7554/eLife.02407.001.

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.

Argonaute proteins couple chromatin silencing to alternative splicing

Argonaute proteins play a major part in transcriptional gene silencing in many organisms, but their role in the nucleus of somatic mammalian cells remains elusive. Here, we have immunopurified human Argonaute-1 and Argonaute-2 (AGO1 and AGO2) chromatin-embedded proteins and found them associated with chromatin modifiers and, notably, with splicing factors. Using the CD44 gene as a model, we show that AGO1 and AGO2 facilitate spliceosome recruitment and modulate RNA polymerase II elongation rate, thereby affecting alternative splicing. Proper AGO1 and AGO2 recruitment to CD44 transcribed regions required the endonuclease Dicer and the chromobox protein HP1γ, and resulted in increased histone H3 lysine 9 methylation on variant exons. Our data thus uncover a new model for the regulation of alternative splicing, in which Argonaute proteins couple RNA polymerase II elongation to chromatin modification.

MAnorm: a robust model for quantitative comparison of ChIP-Seq data sets

ChIP-Seq is widely used to characterize genome-wide binding patterns of transcription factors and other chromatin-associated proteins. Although comparison of ChIP-Seq data sets is critical for understanding cell type-dependent and cell state-specific binding, and thus the study of cell-specific gene regulation, few quantitative approaches have been developed. Here, we present a simple and effective method, MAnorm, for quantitative comparison of ChIP-Seq data sets describing transcription factor binding sites and epigenetic modifications. The quantitative binding differences inferred by MAnorm showed strong correlation with both the changes in expression of target genes and the binding of cell type-specific regulators.

High resolution methylome map of rat indicates role of intragenic DNA methylation in identification of coding region

DNA methylation is crucial for gene regulation and maintenance of genomic stability. Rat has been a key model system in understanding mammalian systemic physiology, however detailed rat methylome remains uncharacterized till date. Here, we present the first high resolution methylome of rat liver generated using Methylated DNA immunoprecipitation and high throughput sequencing (MeDIP-Seq) approach. We observed that within the DNA/RNA repeat elements, simple repeats harbor the highest degree of methylation. Promoter hypomethylation and exon hypermethylation were common features in both RefSeq genes and expressed genes (as evaluated by proteomic approach). We also found that although CpG islands were generally hypomethylated, about 6% of them were methylated and a large proportion (37%) of methylated islands fell within the exons. Notably, we obeserved significant differences in methylation of terminal exons (UTRs); methylation being more pronounced in coding/partially coding exons compared to the non-coding exons. Further, events like alternate exon splicing (cassette exon) and intron retentions were marked by DNA methylation and these regions are retained in the final transcript. Thus, we suggest that DNA methylation could play a crucial role in marking coding regions thereby regulating alternative splicing. Apart from generating the first high resolution methylome map of rat liver tissue, the present study provides several critical insights into methylome organization and extends our understanding of interplay between epigenome, gene expression and genome stability.

Co-transcriptional regulation of alternative pre-mRNA splicing

While studies of alternative pre-mRNA splicing regulation have typically focused on RNA-binding proteins and their target sequences within nascent message, it is becoming increasingly evident that mRNA splicing, RNA polymerase II (pol II) elongation and chromatin structure are intricately intertwined. The majority of introns in higher eukaryotes are excised prior to transcript release in a manner that is dependent on transcription through pol II. As a result of co-transcriptional splicing, variations in pol II elongation influence alternative splicing patterns, wherein a slower elongation rate is associated with increased inclusion of alternative exons within mature mRNA. Physiological barriers to pol II elongation, such as repressive chromatin structure, can thereby similarly impact splicing decisions. Surprisingly, pre-mRNA splicing can reciprocally influence pol II elongation and chromatin structure. Here, we highlight recent advances in co-transcriptional splicing that reveal an extensive network of coupling between splicing, transcription and chromatin remodeling complexes. This article is part of a Special Issue entitled: Chromatin in time and space.

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.

Suppression of Myc oncogenic activity by nucleostemin haploinsufficiency

Nucleostemin (NS), a nucleolar GTPase, is highly expressed in stem/progenitor cells and in most cancer cells. However, little is known about the regulation of its expression. Here, we identify the NS gene as a novel direct transcriptional target of the c-Myc oncoprotein. We show that Myc overexpression enhances NS transcription in cultured cells and in pre-neoplastic B cells from Eμ-myc transgenic mice. Consistent with NS being downstream of Myc, NS expression parallels that of Myc in a large panel of human cancer cell lines. Using chromatin immunoprecipitation we show that c-Myc binds to a well-conserved E-box in the NS promoter. Critically, we show NS haploinsufficiency profoundly delays Myc-induced cancer formation in vivo. NS+/-Eμ-myc transgenic mice have much slower rates of B-cell lymphoma development, with life spans twice that of their wild-type littermates. Moreover, we demonstrate that NS is essential for the proliferation of Myc-overexpressing cells in cultured cells and in vivo: impaired lymphoma development was associated with a drastic decrease of c-Myc-induced proliferation of pre-tumoural B cells. Finally, we provide evidence that in cell culture NS controls cell proliferation independently of p53 and that NS haploinsufficiency significantly delays lymphomagenesis in p53-deficient mice. Together these data indicate that NS functions downstream of Myc as a rate-limiting regulator of cell proliferation and transformation, independently from its putative role within the p53 pathway. Targeting NS is therefore expected to compromise early tumour development irrespectively of the p53 status.

RNA Polymerase II Elongation at the Crossroads of Transcription and Alternative Splicing

The elongation phase of transcription lies at the core of several simultaneous and coupled events leading to alternative splicing regulation. Although underestimated in the past, it is at this phase of the transcription cycle where complexes affecting the transcription machinery itself, chromatin structure, posttranscriptional gene regulation and pre-mRNA processing converge to regulate each other or simply to consolidate higher-order complexes and functions. This paper focuses on the multiple processes that take place during transcription elongation which ultimately regulate the outcome of alternative splicing decisions.

Occupancy of chromatin organizers in the Epstein-Barr virus genome

The human CCCTC-binding factor, CTCF, regulates transcription of the double-stranded DNA genomes of herpesviruses. The architectural complex cohesin and RNA Polymerase II also contribute to this organization. We profiled the occupancy of CTCF, cohesin, and RNA Polymerase II on the episomal genome of the Epstein-Barr virus in a cell culture model of latent infection. CTCF colocalizes with cohesin but not RNA Polymerase II. CTCF and cohesin bind specific sequences throughout the genome that are found not just proximal to the regulatory elements of latent genes, but also near lytic genes. In addition to tracking with known transcripts, RNA Polymerase II appears at two unannotated positions, one of which lies within the latent origin of replication. The widespread occupancy profile of each protein reveals binding near or at a myriad of regulatory elements and suggests context-dependent functions.

Pause locally, splice globally

Most eukaryotic protein-coding transcripts contain introns, which vary in number and position along the transcript body. Intron removal through pre-mRNA splicing is tightly linked to transcription by RNA polymerase II as it translocates along each gene. Here, we review recent evidence that transcription and splicing are functionally coupled. We focus on how RNA polymerase II elongation rates impact splicing through local regulation and transcriptional pausing within genes. Emerging concepts of how splicing-related changes in elongation might be achieved are highlighted. We place the interplay between transcription and splicing in the context of chromatin where nucleosome positioning influences elongation, and histone modifications participate directly in the recruitment of splicing regulators to nascent transcripts.

More than a splicing code: integrating the role of RNA, chromatin and non-coding RNA in alternative splicing regulation

Large portions of the genome undergo alternative pre-mRNA splicing in often intricate patterns. Alternative splicing regulation requires extensive control mechanisms since errors can have deleterious consequences and may lead to developmental defects and disease. Recent work has identified a complex network of regulatory RNA elements which guide splicing decisions. In addition, the discovery that transcription and splicing are intimately coupled has opened up new directions into alternative splicing regulation. Work at the interface of chromatin and RNA biology has revealed unexpected molecular links between histone modifications, the transcription machinery, and non-coding RNAs (ncRNAs) in the determination of alternative splicing patterns.

Epigenetics in alternative pre-mRNA splicing

Alternative splicing plays critical roles in differentiation, development, and disease and is a major source for protein diversity in higher eukaryotes. Analysis of alternative splicing regulation has traditionally focused on RNA sequence elements and their associated splicing factors, but recent provocative studies point to a key function of chromatin structure and histone modifications in alternative splicing regulation. These insights suggest that epigenetic regulation determines not only what parts of the genome are expressed but also how they are spliced.

Genomic-wide methods to evaluate transcription rates in yeast

Gene transcription is a dynamic process in which the desired amount of an mRNA is obtained by the equilibrium between its transcription (TR) and degradation (DR) rates. The control mechanism at the RNA polymerase level primarily causes changes in TR. Despite their importance, TRs have been rarely measured. In the yeast Saccharomyces cerevisiae, we have implemented two techniques to evaluate TRs: run-on and chromatin immunoprecipitation of RNA polymerase II. These techniques allow the discrimination of the relative importance of TR and DR in gene regulation for the first time in a eukaryote.

Splicing-dependent RNA polymerase pausing in yeast

In eukaryotic cells, there is evidence for functional coupling between transcription and processing of pre-mRNAs. To better understand this coupling, we performed a high-resolution kinetic analysis of transcription and splicing in budding yeast. This revealed that shortly after induction of transcription, RNA polymerase accumulates transiently around the 3′ end of the intron on two reporter genes. This apparent transcriptional pause coincides with splicing factor recruitment and with the first detection of spliced mRNA and is repeated periodically thereafter. Pausing requires productive splicing, as it is lost upon mutation of the intron and restored by suppressing the splicing defect. The carboxy-terminal domain of the paused polymerase large subunit is hyperphosphorylated on serine 5, and phosphorylation of serine 2 is first detected here. Phosphorylated polymerase also accumulates around the 3′ splice sites of constitutively expressed, endogenous yeast genes. We propose that transcriptional pausing is imposed by a checkpoint associated with cotranscriptional splicing.

Global analysis of nascent RNA reveals transcriptional pausing in terminal exons

Pre-mRNA splicing is catalyzed by the spliceosome, which can assemble on pre-mRNA cotranscriptionally. However, whether splicing generally occurs during transcription has not been addressed. Indeed, splicing catalysis is expected to occur posttranscriptionally in yeast, where the shortness of terminal exons should leave insufficient time for splicing. Here, we isolate endogenous S. cerevisiae nascent RNA and determine gene-specific splicing efficiencies and transcription profiles, using high-density tiling microarrays. Surprisingly, we find that splicing occurs cotranscriptionally for the majority of intron-containing genes. Analysis of transcription profiles reveals Pol II pausing within the terminal exons of these genes. Intronless and inefficiently spliced genes lack this pause. In silico simulations of transcription and splicing kinetics confirm that this pausing event provides sufficient time for splicing before termination. The discovery of terminal exon pausing demonstrates functional coupling of transcription and splicing near gene ends.

Noncanonical transcript forms in yeast and their regulation during environmental stress

Surveys of transcription in many organisms have observed widespread expression of RNAs with no known function, encoded within and between canonical coding genes. The search to distinguish functional RNAs from transcriptional noise represents one of the great challenges in genomic biology. Here we report a next-generation sequencing technique designed to facilitate the inference of function of uncharacterized transcript forms by improving their coverage in sequencing libraries, in parallel with the detection of canonical mRNAs. We piloted this protocol, which is based on the capture of 3' ends of polyadenylated RNAs, in budding yeast. Analysis of transcript ends in coding regions uncovered hundreds of alternative-length coding forms, which harbored a unique sequence motif and showed signatures of regulatory function in particular gene categories; independent single-gene measurements confirmed the differential regulation of short coding forms during heat shock. In addition, our 3'-end RNA-seq method applied to wild-type strains detected putative noncoding transcripts previously reported only in RNA surveillance mutants, and many such transcripts showed differential expression in yeast cultures grown under chemical stress. Our results underscore the power of the 3'-end protocol to improve detection of noncanonical transcript forms in a sequencing experiment of standard depth, and our findings strongly suggest that many unannotated, polyadenylated RNAs may have as yet uncharacterized regulatory functions.

Widespread transcription at neuronal activity-regulated enhancers

We used genome-wide sequencing methods to study stimulus-dependent enhancer function in neurons. We identified ∼12,000 neuronal activity-regulated enhancers that are bound by the general transcriptional co-activator CBP in an activity-dependent manner. A function of CBP at enhancers may be to recruit RNA polymerase II (RNAPII), as we also observed activity-regulated RNAPII binding to thousands of enhancers. Remarkably, RNAPII at enhancers transcribes bi-directionally a novel class of enhancer RNAs (eRNAs) within enhancer domains defined by the presence of histone H3 that is mono-methylated at lysine 4 (H3K4me1). The level of eRNA expression at neuronal enhancers positively correlates with the level of mRNA synthesis at nearby genes, suggesting that eRNA synthesis occurs specifically at enhancers that are actively engaged in promoting mRNA synthesis. These findings reveal that a widespread mechanism of enhancer activation involves RNAPII binding and eRNA synthesis.

PPS, a large multidomain protein, functions with sex-lethal to regulate alternative splicing in Drosophila

Alternative splicing controls the expression of many genes, including the Drosophila sex determination gene Sex-lethal (Sxl). Sxl expression is controlled via a negative regulatory mechanism where inclusion of the translation-terminating male exon is blocked in females. Previous studies have shown that the mechanism leading to exon skipping is autoregulatory and requires the SXL protein to antagonize exon inclusion by interacting with core spliceosomal proteins, including the U1 snRNP protein Sans-fille (SNF). In studies begun by screening for proteins that interact with SNF, we identified PPS, a previously uncharacterized protein, as a novel component of the machinery required for Sxl male exon skipping. PPS encodes a large protein with four signature motifs, PHD, BRK, TFS2M, and SPOC, typically found in proteins involved in transcription. We demonstrate that PPS has a direct role in Sxl male exon skipping by showing first that loss of function mutations have phenotypes indicative of Sxl misregulation and second that the PPS protein forms a complex with SXL and the unspliced Sxl RNA. In addition, we mapped the recruitment of PPS, SXL, and SNF along the Sxl gene using chromatin immunoprecipitation (ChIP), which revealed that, like many other splicing factors, these proteins bind their RNA targets while in close proximity to the DNA. Interestingly, while SNF and SXL are specifically recruited to their predicted binding sites, PPS has a distinct pattern of accumulation along the Sxl gene, associating with a region that includes, but is not limited to, the SxlPm promoter. Together, these data indicate that PPS is different from other splicing factors involved in male-exon skipping and suggest, for the first time, a functional link between transcription and SXL-mediated alternative splicing. Loss of zygotic PPS function, however, is lethal to both sexes, indicating that its role may be of broad significance.

The organization of nucleosomes around splice sites

The occupancy of nucleosomes along chromosome is a key factor for gene regulation. However, except promoter regions, genome-wide properties and functions of nucleosome organization remain unclear in mammalian genomes. Using the computational model of Increment of Diversity with Quadratic Discriminant (IDQD) trained from the microarray data, the nucleosome occupancy score (NOScore) was defined and applied to splice junction regions of constitutive, cassette exon, alternative 3′ and 5′ splicing events in the human genome. We found an interesting relation between NOScore and RNA splicing: exon regions have higher NOScores compared with their flanking intron sequences in both constitutive and alternative splicing events, indicating the stronger nucleosome occupation potential of exon regions. In addition, NOScore valleys present at ∼25 bp upstream of the acceptor site in all splicing events. By defining folding diversity-to-energy ratio to describe RNA structural flexibility, we demonstrated that primary RNA transcripts from nucleosome occupancy regions are relatively rigid and those from nucleosome depleted regions are relatively flexible. The negative correlation between nucleosome occupation/depletion of DNA sequence and structural flexibility/rigidity of its primary transcript around splice junctions may provide clues to the deeper understanding of the unexpected role for nucleosome organization in the regulation of RNA splicing.

Chromatin organization marks exon-intron structure

An increasing body of evidence indicates that transcription and splicing are coupled, and it is accepted that chromatin organization regulates transcription. Little is known about the cross-talk between chromatin structure and exon-intron architecture. By analysis of genome-wide nucleosome-positioning data sets from humans, flies and worms, we found that exons show increased nucleosome-occupancy levels with respect to introns, a finding that we link to differential GC content and nucleosome-disfavoring elements between exons and introns. Analysis of genome-wide chromatin immunoprecipitation data in humans and mice revealed four specific post-translational histone modifications enriched in exons. Our findings indicate that previously described enrichment of H3K36me3 modifications in exons reflects a more fundamental phenomenon, namely increased nucleosome occupancy along exons. Our results suggest an RNA polymerase II-mediated cross-talk between chromatin structure and exon-intron architecture, implying that exon selection may be modulated by chromatin structure.