Functional association between promoter structure and transcript alternative splicing

Predicting splicing patterns from the transcription factor binding sites in the promoter with deep learning

BACKGROUND

Alternative splicing is a pivotal mechanism of post-transcriptional modification that contributes to the transcriptome plasticity and proteome diversity in metazoan cells. Although many splicing regulations around the exon/intron regions are known, the relationship between promoter-bound transcription factors and the downstream alternative splicing largely remains unexplored.

RESULTS

In this study, we present computational approaches to unravel the regulatory relationship between promoter-bound transcription factor binding sites (TFBSs) and the splicing patterns. We curated a fine dataset that includes DNase I hypersensitive site sequencing and transcriptomes across fifteen human tissues from ENCODE. Specifically, we proposed different representations of TF binding context and splicing patterns to examine the associations between the promoter and downstream splicing events. While machine learning models demonstrated potential in predicting splicing patterns based on TFBS occupancies, the limitations in the generalization of predicting the splicing forms of singleton genes across diverse tissues was observed with carefully examination using different cross-validation methods. We further investigated the association between alterations in individual TFBS at promoters and shifts in exon splicing efficiency. Our results demonstrate that the convolutional neural network (CNN) models, trained on TF binding changes in the promoters, can predict the changes in splicing patterns. Furthermore, a systemic in silico substitutions analysis on the CNN models highlighted several potential splicing regulators. Notably, using empirical validation using K562 CTCFL shRNA knock-down data, we showed the significant role of CTCFL in splicing regulation.

CONCLUSION

In conclusion, our finding highlights the potential role of promoter-bound TFBSs in influencing the regulation of downstream splicing patterns and provides insights for discovering alternative splicing regulations.

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.

A group 3 medulloblastoma stem cell program is maintained by OTX2-mediated alternative splicing

OTX2 is a transcription factor and known driver in medulloblastoma (MB), where it is amplified in a subset of tumours and overexpressed in most cases of group 3 and group 4 MB. Here we demonstrate a noncanonical role for OTX2 in group 3 MB alternative splicing. OTX2 associates with the large assembly of splicing regulators complex through protein-protein interactions and regulates a stem cell splicing program. OTX2 can directly or indirectly bind RNA and this may be partially independent of its DNA regulatory functions. OTX2 controls a pro-tumorigenic splicing program that is mirrored in human cerebellar rhombic lip origins. Among the OTX2-regulated differentially spliced genes, PPHLN1 is expressed in the most primitive rhombic lip stem cells, and targeting PPHLN1 splicing reduces tumour growth and enhances survival in vivo. These findings identify OTX2-mediated alternative splicing as a major determinant of cell fate decisions that drive group 3 MB progression.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

(Alternative) transcription start sites as regulators of RNA processing

Alternative transcription start site usage (ATSS) is a widespread regulatory strategy that enables genes to choose between multiple genomic loci for initiating transcription. This mechanism is tightly controlled during development and is often altered in disease states. In this review, we examine the growing evidence highlighting a role for transcription start sites (TSSs) in the regulation of mRNA isoform selection during and after transcription. We discuss how the choice of transcription initiation sites influences RNA processing and the importance of this crosstalk for cell identity and organism function. We also speculate on possible mechanisms underlying the integration of transcriptional and post-transcriptional processes.

Alternative splicing is coupled to gene expression in a subset of variably expressed genes

Numerous factors regulate alternative splicing of human genes at a co-transcriptional level. However, how alternative splicing depends on the regulation of gene expression is poorly understood. We leveraged data from the Genotype-Tissue Expression (GTEx) project to show a significant association of gene expression and splicing for 6874 (4.9%) of 141,043 exons in 1106 (13.3%) of 8314 genes with substantially variable expression in ten GTEx tissues. About half of these exons demonstrate higher inclusion with higher gene expression, and half demonstrate higher exclusion, with the observed direction of coupling being highly consistent across different tissues and in external datasets. The exons differ with respect to sequence characteristics, enriched sequence motifs, RNA polymerase II binding, and inferred transcription rate of downstream introns. The exons were enriched for hundreds of isoform-specific Gene Ontology annotations, suggesting that the coupling of expression and alternative splicing described here may provide an important gene regulatory mechanism that might be used in a variety of biological contexts. In particular, higher inclusion exons could play an important role during cell division.

Coalescent RNA-localizing and transcriptional activities of SAM68 modulate adhesion and subendothelial basement membrane assembly

Endothelial cell interactions with their extracellular matrix are essential for vascular homeostasis and expansion. Large-scale proteomic analyses aimed at identifying components of integrin adhesion complexes have revealed the presence of several RNA binding proteins (RBPs) of which the functions at these sites remain poorly understood. Here, we explored the role of the RBP SAM68 (Src associated in mitosis, of 68 kDa) in endothelial cells. We found that SAM68 is transiently localized at the edge of spreading cells where it participates in membrane protrusive activity and the conversion of nascent adhesions to mechanically loaded focal adhesions by modulation of integrin signaling and local delivery of β-actin mRNA. Furthermore, SAM68 depletion impacts cell-matrix interactions and motility through induction of key matrix genes involved in vascular matrix assembly. In a 3D environment SAM68-dependent functions in both tip and stalk cells contribute to the process of sprouting angiogenesis. Altogether, our results identify the RBP SAM68 as a novel actor in the dynamic regulation of blood vessel networks.

Isoforms of the TAL1 transcription factor have different roles in hematopoiesis and cell growth

T-cell acute lymphoblastic leukemia (T-ALL) protein 1 (TAL1) is a central transcription factor in hematopoiesis. The timing and level of TAL1 expression orchestrate the differentiation to specialized blood cells and its overexpression is a common cause of T-ALL. Here, we studied the 2 protein isoforms of TAL1, short and long, which are generated by the use of alternative promoters as well as by alternative splicing. We analyzed the expression of each isoform by deleting an enhancer or insulator, or by opening chromatin at the enhancer location. Our results show that each enhancer promotes expression from a specific TAL1 promoter. Expression from a specific promoter gives rise to a unique 5' UTR with differential regulation of translation. Moreover, our study suggests that the enhancers regulate TAL1 exon 3 alternative splicing by inducing changes in the chromatin at the splice site, which we demonstrate is mediated by KMT2B. Furthermore, our results indicate that TAL1-short binds more strongly to TAL1 E-protein partners and functions as a stronger transcription factor than TAL1-long. Specifically TAL1-short has a unique transcription signature promoting apoptosis. Finally, when we expressed both isoforms in mice bone marrow, we found that while overexpression of both isoforms prevents lymphoid differentiation, expression of TAL-short alone leads to hematopoietic stem cell exhaustion. Furthermore, we found that TAL1-short promoted erythropoiesis and reduced cell survival in the CML cell line K562. While TAL1 and its partners are considered promising therapeutic targets in the treatment of T-ALL, our results show that TAL1-short could act as a tumor suppressor and suggest that altering TAL1 isoform's ratio could be a preferred therapeutic approach.

Dnmt3bas coordinates transcriptional induction and alternative exon inclusion to promote catalytically active Dnmt3b expression

Embryonic expression of DNMT3B is critical for establishing de novo DNA methylation. This study uncovers the mechanism through which the promoter-associated long non-coding RNA (lncRNA) Dnmt3bas controls the induction and alternative splicing of Dnmt3b during embryonic stem cell (ESC) differentiation. Dnmt3bas recruits the PRC2 (polycomb repressive complex 2) at cis-regulatory elements of the Dnmt3b gene expressed at a basal level. Correspondingly, Dnmt3bas knockdown enhances Dnmt3b transcriptional induction, whereas overexpression of Dnmt3bas dampens it. Dnmt3b induction coincides with exon inclusion, switching the predominant isoform from the inactive Dnmt3b6 to the active Dnmt3b1. Intriguingly, overexpressing Dnmt3bas further enhances the Dnmt3b1:Dnmt3b6 ratio, attributed to its interaction with hnRNPL (heterogeneous nuclear ribonucleoprotein L), a splicing factor that promotes exon inclusion. Our data suggest that Dnmt3bas coordinates alternative splicing and transcriptional induction of Dnmt3b by facilitating the hnRNPL and RNA polymerase II (RNA Pol II) interaction at the Dnmt3b promoter. This dual mechanism precisely regulates the expression of catalytically active DNMT3B, ensuring fidelity and specificity of de novo DNA methylation.

From words to complete phrases: insight into single-cell isoforms using short and long reads

The profiling of gene expression patterns to glean biological insights from single cells has become commonplace over the last few years. However, this approach overlooks the transcript contents that can differ between individual cells and cell populations. In this review, we describe early work in the field of single-cell short-read sequencing as well as full-length isoforms from single cells. We then describe recent work in single-cell long-read sequencing wherein some transcript elements have been observed to work in tandem. Based on earlier work in bulk tissue, we motivate the study of combination patterns of other RNA variables. Given that we are still blind to some aspects of isoform biology, we suggest possible future avenues such as CRISPR screens which can further illuminate the function of RNA variables in distinct cell populations.

Sites of transcription initiation drive mRNA isoform selection

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

Integrating transcription and splicing into cell fate: Transcription factors on the block

Transcription factors (TFs) are present in all life forms and conserved across great evolutionary distances in eukaryotes. From yeast to complex multicellular organisms, they are pivotal players of cell fate decision by orchestrating gene expression at diverse molecular layers. Notably, TFs fine-tune gene expression by coordinating RNA fate at both the expression and splicing levels. They regulate alternative splicing, an essential mechanism for cell plasticity, allowing the production of many mRNA and protein isoforms in precise cell and tissue contexts. Despite this apparent role in splicing, how TFs integrate transcription and splicing to ultimately orchestrate diverse cell functions and cell fate decisions remains puzzling. We depict substantial studies in various model organisms underlining the key role of TFs in alternative splicing for promoting tissue-specific functions and cell fate. Furthermore, we emphasize recent advances describing the molecular link between the transcriptional and splicing activities of TFs. As TFs can bind both DNA and/or RNA to regulate transcription and splicing, we further discuss their flexibility and compatibility for DNA and RNA substrates. Finally, we propose several models integrating transcription and splicing activities of TFs in the coordination and diversification of cell and tissue identities. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > Splicing Mechanisms.

Single-nuclei isoform RNA sequencing unlocks barcoded exon connectivity in frozen brain tissue

Single-nuclei RNA sequencing characterizes cell types at the gene level. However, compared to single-cell approaches, many single-nuclei cDNAs are purely intronic, lack barcodes and hinder the study of isoforms. Here we present single-nuclei isoform RNA sequencing (SnISOr-Seq). Using microfluidics, PCR-based artifact removal, target enrichment and long-read sequencing, SnISOr-Seq increased barcoded, exon-spanning long reads 7.5-fold compared to naive long-read single-nuclei sequencing. We applied SnISOr-Seq to adult human frontal cortex and found that exons associated with autism exhibit coordinated and highly cell-type-specific inclusion. We found two distinct combination patterns: those distinguishing neural cell types, enriched in TSS-exon, exon-polyadenylation-site and non-adjacent exon pairs, and those with multiple configurations within one cell type, enriched in adjacent exon pairs. Finally, we observed that human-specific exons are almost as tightly coordinated as conserved exons, implying that coordination can be rapidly established during evolution. SnISOr-Seq enables cell-type-specific long-read isoform analysis in human brain and in any frozen or hard-to-dissociate sample.

RBFOX2-regulated TEAD1 alternative splicing plays a pivotal role in Hippo-YAP signaling

Alternative pre-mRNA splicing is key to proteome diversity; however, the biological roles of alternative splicing (AS) in signaling pathways remain elusive. Here, we focus on TEA domain transcription factor 1 (TEAD1), a YAP binding factor in the Hippo signaling pathway. Public database analyses showed that expression of YAP-TEAD target genes negatively correlated with the expression of a TEAD1 isoform lacking exon 6 (TEAD1ΔE6) but did not correlate with overall TEAD1 expression. We confirmed that the transcriptional activity and oncogenic properties of the full-length TEAD1 isoform were greater than those of TEAD1ΔE6, with the difference in transcription related to YAP interaction. Furthermore, we showed that RNA-binding Fox-1 homolog 2 (RBFOX2) promoted the inclusion of TEAD1 exon 6 via binding to the conserved GCAUG element in the downstream intron. These results suggest a regulatory mechanism of RBFOX2-mediated TEAD1 AS and provide insight into AS-specific modulation of signaling pathways.

The Hox transcription factor Ultrabithorax binds RNA and regulates co-transcriptional splicing through an interplay with RNA polymerase II

Transcription factors (TFs) play a pivotal role in cell fate decision by coordinating gene expression programs. Although most TFs act at the DNA layer, few TFs bind RNA and modulate splicing. Yet, the mechanistic cues underlying TFs activity in splicing remain elusive. Focusing on the Drosophila Hox TF Ultrabithorax (Ubx), our work shed light on a novel layer of Ubx function at the RNA level. Transcriptome and genome-wide binding profiles in embryonic mesoderm and Drosophila cells indicate that Ubx regulates mRNA expression and splicing to promote distinct outcomes in defined cellular contexts. Our results demonstrate a new RNA-binding ability of Ubx. We find that the N51 amino acid of the DNA-binding Homeodomain is non-essential for RNA interaction in vitro, but is required for RNA interaction in vivo and Ubx splicing activity. Moreover, mutation of the N51 amino acid weakens the interaction between Ubx and active RNA Polymerase II (Pol II). Our results reveal that Ubx regulates elongation-coupled splicing, which could be coordinated by a dynamic interplay with active Pol II on chromatin. Overall, our work uncovered a novel role of the Hox TFs at the mRNA regulatory layer. This could be an essential function for other classes of TFs to control cell diversity.

Post-transcriptional regulation in spermatogenesis: all RNA pathways lead to healthy sperm

Multiple RNA pathways are required to produce functional sperm. Here, we review RNA post-transcriptional regulation during spermatogenesis with particular emphasis on the role of 3' end modifications. From early studies in the 1970s, it became clear that spermiogenesis transcripts could be stored for days only to be translated at advanced stages of spermatid differentiation. The transition between the translationally repressed and active states was observed to correlate with the shortening of the transcripts' poly(A) tail, establishing a link between RNA 3' end metabolism and male germ cell differentiation. Since then, numerous RNA metabolic pathways have been implicated not only in the progression through spermatogenesis, but also in the maintenance of genomic integrity. Recent studies have characterized the elusive 3' biogenesis of Piwi-interacting RNAs (piRNAs), identified a critical role for messenger RNA (mRNA) 3' uridylation in meiotic progression, established the mechanisms that destabilize transcripts with long 3' untranslated regions (3'UTRs) in post-mitotic cells, and defined the physiological relevance of RNA exonucleases and deadenylases in male germ cells. In this review, we discuss RNA processing in the male germline in the light of the most recent findings. A brief recollection of different RNA-processing events will aid future studies exploring post-transcriptional regulation in spermatogenesis.

VEGFA's distal enhancer regulates its alternative splicing in CML

Enhancer demethylation in leukemia has been shown to lead to overexpression of genes which promote cancer characteristics. The vascular endothelial growth factor A (VEGFA) enhancer, located 157 Kb downstream of its promoter, is demethylated in chronic myeloid leukemia (CML). VEGFA has several alternative splicing isoforms with different roles in cancer progression. Since transcription and splicing are coupled, we wondered whether VEGFA enhancer activity can also regulate the gene's alternative splicing to contribute to the pathology of CML. Our results show that mutating the VEGFA +157 enhancer promotes exclusion of exons 6a and 7 and activating the enhancer by tethering a chromatin activator has the opposite effect. In line with these results, CML patients present with high expression of +157 eRNA and inclusion of VEGFA exons 6a and 7. In addition, our results show that the positive regulator of RNAPII transcription elongation, CCNT2, binds VEGFA's promoter and enhancer, and its silencing promotes exclusion of exons 6a and 7 as it slows down RNAPII elongation rate. Thus our results suggest that VEGFA's +157 enhancer regulates its alternative splicing by increasing RNAPII elongation rate via CCNT2. Our work demonstrates for the first time a connection between an endogenous enhancer and alternative splicing regulation of its target gene.

Efficient RNA polymerase II pause release requires U2 snRNP function

Transcription by RNA polymerase II (Pol II) is coupled to pre-mRNA splicing, but the underlying mechanisms remain poorly understood. Co-transcriptional splicing requires assembly of a functional spliceosome on nascent pre-mRNA, but whether and how this influences Pol II transcription remains unclear. Here we show that inhibition of pre-mRNA branch site recognition by the spliceosome component U2 snRNP leads to a widespread and strong decrease in new RNA synthesis from human genes. Multiomics analysis reveals that inhibition of U2 snRNP function increases the duration of Pol II pausing in the promoter-proximal region, impairs recruitment of the pause release factor P-TEFb, and reduces Pol II elongation velocity at the beginning of genes. Our results indicate that efficient release of paused Pol II into active transcription elongation requires the formation of functional spliceosomes and that eukaryotic mRNA biogenesis relies on positive feedback from the splicing machinery to the transcription machinery.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

An emerging role of chromatin-interacting RNA-binding proteins in transcription regulation

Transcription factors (TFs) are well-established key factors orchestrating gene transcription, and RNA-binding proteins (RBPs) are mainly thought to participate in post-transcriptional control of gene. In fact, these two steps are functionally coupled, offering a possibility for reciprocal communications between transcription and regulatory RNAs and RBPs. Recently, a series of exploratory studies, utilizing functional genomic strategies, have revealed that RBPs are prevalently involved in transcription control genome-wide through their interactions with chromatin. Here, we present a refined census of RBPs to grope for such an emerging role and discuss the global view of RBP-chromatin interactions and their functional diversities in transcription regulation.

Computational identification of eukaryotic promoters based on cascaded deep capsule neural networks

A promoter is a region in the DNA sequence that defines where the transcription of a gene by RNA polymerase initiates, which is typically located proximal to the transcription start site (TSS). How to correctly identify the gene TSS and the core promoter is essential for our understanding of the transcriptional regulation of genes. As a complement to conventional experimental methods, computational techniques with easy-to-use platforms as essential bioinformatics tools can be effectively applied to annotate the functions and physiological roles of promoters. In this work, we propose a deep learning-based method termed Depicter (Deep learning for predicting promoter), for identifying three specific types of promoters, i.e. promoter sequences with the TATA-box (TATA model), promoter sequences without the TATA-box (non-TATA model), and indistinguishable promoters (TATA and non-TATA model). Depicter is developed based on an up-to-date, species-specific dataset which includes Homo sapiens, Mus musculus, Drosophila melanogaster and Arabidopsis thaliana promoters. A convolutional neural network coupled with capsule layers is proposed to train and optimize the prediction model of Depicter. Extensive benchmarking and independent tests demonstrate that Depicter achieves an improved predictive performance compared with several state-of-the-art methods. The webserver of Depicter is implemented and freely accessible at https://depicter.erc.monash.edu/.

Role of Viral Ribonucleoproteins in Human Papillomavirus Type 16 Gene Expression

Human papillomaviruses (HPVs) depend on the cellular RNA-processing machineries including alternative RNA splicing and polyadenylation to coordinate HPV gene expression. HPV RNA processing is controlled by cis-regulatory RNA elements and trans-regulatory factors since the HPV splice sites are suboptimal. The definition of HPV exons and introns may differ between individual HPV mRNA species and is complicated by the fact that many HPV protein-coding sequences overlap. The formation of HPV ribonucleoproteins consisting of HPV pre-mRNAs and multiple cellular RNA-binding proteins may result in the different outcomes of HPV gene expression, which contributes to the HPV life cycle progression and HPV-associated cancer development. In this review, we summarize the regulation of HPV16 gene expression at the level of RNA processing with focus on the interactions between HPV16 pre-mRNAs and cellular RNA-binding factors.

Revisiting the window of opportunity for cotranscriptional splicing in budding yeast

We reported previously that, in budding yeast, transcription rate affects both the efficiency and fidelity of pre-mRNA splicing, especially of ribosomal protein transcripts. Here, we report that the majority of ribosomal protein transcripts with non-consensus 5' splice sites are spliced less efficiently when transcription is faster, and more efficiently with slower transcription. These results support the "window of opportunity" model, and we suggest a possible mechanism to explain these findings.

Linking transcription, RNA polymerase II elongation and alternative splicing

Gene expression is an intricately regulated process that is at the basis of cell differentiation, the maintenance of cell identity and the cellular responses to environmental changes. Alternative splicing, the process by which multiple functionally distinct transcripts are generated from a single gene, is one of the main mechanisms that contribute to expand the coding capacity of genomes and help explain the level of complexity achieved by higher organisms. Eukaryotic transcription is subject to multiple layers of regulation both intrinsic — such as promoter structure — and dynamic, allowing the cell to respond to internal and external signals. Similarly, alternative splicing choices are affected by all of these aspects, mainly through the regulation of transcription elongation, making it a regulatory knob on a par with the regulation of gene expression levels. This review aims to recapitulate some of the history and stepping-stones that led to the paradigms held today about transcription and splicing regulation, with major focus on transcription elongation and its effect on alternative splicing.

Light in the transcription landscape: chromatin, RNA polymerase II and splicing throughout Arabidopsis thaliana's life cycle

Plants have a high level of developmental plasticity that allows them to respond and adapt to changes in the environment. Among the environmental cues, light controls almost every aspect of A. thaliana's life cycle, including seed maturation, seed germination, seedling de-etiolation and flowering time. Light signals induce massive reprogramming of gene expression, producing changes in RNA polymerase II transcription, alternative splicing, and chromatin state. Since splicing reactions occur mainly while transcription takes place, the regulation of RNAPII transcription has repercussions in the splicing outcomes. This cotranscriptional nature allows a functional coupling between transcription and splicing, in which properties of the splicing reactions are affected by the transcriptional process. Chromatin landscapes influence both transcription and splicing. In this review, we highlight, summarize and discuss recent progress in the field to gain a comprehensive insight on the cross-regulation between chromatin state, RNAPII transcription and splicing decisions in plants, with a special focus on light-triggered responses. We also introduce several examples of transcription and splicing factors that could be acting as coupling factors in plants. Unravelling how these connected regulatory networks operate, can help in the design of better crops with higher productivity and tolerance.

Reciprocal regulation between alternative splicing and the DNA damage response

Splicing, the process that catalyzes intron removal and flanking exon ligation, can occur in different ways (alternative splicing) in immature RNAs transcribed from a single gene. In order to adapt to a particular context, cells modulate not only the quantity but also the quality (alternative isoforms) of their transcriptome. Since 95% of the human coding genome is subjected to alternative splicing regulation, it is expected that many cellular pathways are modulated by alternative splicing, as is the case for the DNA damage response. Moreover, recent evidence demonstrates that upon a genotoxic insult, classical DNA damage response kinases such as ATM, ATR and DNA-PK orchestrate the gene expression response therefore modulating alternative splicing which, in a reciprocal way, shapes the response to a damaging agent.

Mechanistic Dissection of RNA-Binding Proteins in Regulated Gene Expression at Chromatin Levels

Eukaryotic genomes are known to prevalently transcribe diverse classes of RNAs, virtually all of which, including nascent RNAs from protein-coding genes, are now recognized to have regulatory functions in gene expression, suggesting that RNAs are both the products and the regulators of gene expression. Their functions must enlist specific RNA-binding proteins (RBPs) to execute their regulatory activities, and recent evidence suggests that nearly all biochemically defined chromatin regions in the human genome, whether defined for gene activation or silencing, have the involvement of specific RBPs. Interestingly, the boundary between RNA- and DNA-binding proteins is also melting, as many DNA-binding proteins traditionally studied in the context of transcription are able to bind RNAs, some of which may simultaneously bind both DNA and RNA to facilitate network interactions in three-dimensional (3D) genome. In this review, we focus on RBPs that function at chromatin levels, with particular emphasis on their mechanisms of action in regulated gene expression, which is intended to facilitate future functional and mechanistic dissection of chromatin-associated RBPs.

Longitudinal RNA-Seq Analysis of the Repeatability of Gene Expression and Splicing in Human Platelets Identifies a Platelet SELP Splice QTL

RATIONALE

Longitudinal studies are required to distinguish within versus between-individual variation and repeatability of gene expression. They are uniquely positioned to decipher genetic signal from environmental noise, with potential application to gene variant and expression studies. However, longitudinal analyses of gene expression in healthy individuals-especially with regards to alternative splicing-are lacking for most primary cell types, including platelets.

OBJECTIVE

To assess repeatability of gene expression and splicing in platelets and use repeatability to identify novel platelet expression quantitative trait loci (QTLs) and splice QTLs.

METHODS AND RESULTS

We sequenced the transcriptome of platelets isolated repeatedly up to 4 years from healthy individuals. We examined within and between individual variation and repeatability of platelet RNA expression and exon skipping, a readily measured alternative splicing event. We find that platelet gene expression is generally stable between and within-individuals over time-with the exception of a subset of genes enriched for the inflammation gene ontology. We show an enrichment among repeatable genes for associations with heritable traits, including known and novel platelet expression QTLs. Several exon skipping events were also highly repeatable, suggesting heritable patterns of splicing in platelets. One of the most repeatable was exon 14 skipping of SELP. Accordingly, we identify rs6128 as a platelet splice QTL and define an rs6128-dependent association between SELP exon 14 skipping and race. In vitro experiments demonstrate that this single nucleotide variant directly affects exon 14 skipping and changes the ratio of transmembrane versus soluble P-selectin protein production.

CONCLUSIONS

We conclude that the platelet transcriptome is generally stable over 4 years. We demonstrate the use of repeatability of gene expression and splicing to identify novel platelet expression QTLs and splice QTLs. rs6128 is a platelet splice QTL that alters SELP exon 14 skipping and soluble versus transmembrane P-selectin protein production.

ARGLU1 is a transcriptional coactivator and splicing regulator important for stress hormone signaling and development

Stress hormones bind and activate the glucocorticoid receptor (GR) in many tissues including the brain. We identified arginine and glutamate rich 1 (ARGLU1) in a screen for new modulators of glucocorticoid signaling in the CNS. Biochemical studies show that the glutamate rich C-terminus of ARGLU1 coactivates multiple nuclear receptors including the glucocorticoid receptor (GR) and the arginine rich N-terminus interacts with splicing factors and binds to RNA. RNA-seq of neural cells depleted of ARGLU1 revealed significant changes in the expression and alternative splicing of distinct genes involved in neurogenesis. Loss of ARGLU1 is embryonic lethal in mice, and knockdown in zebrafish causes neurodevelopmental and heart defects. Treatment with dexamethasone, a GR activator, also induces changes in the pattern of alternatively spliced genes, many of which were lost when ARGLU1 was absent. Importantly, the genes found to be alternatively spliced in response to glucocorticoid treatment were distinct from those under transcriptional control by GR, suggesting an additional mechanism of glucocorticoid action is present in neural cells. Our results thus show that ARGLU1 is a novel factor for embryonic development that modulates basal transcription and alternative splicing in neural cells with consequences for glucocorticoid signaling.

FLAM-seq: full-length mRNA sequencing reveals principles of poly(A) tail length control

Although messenger RNAs are key molecules for understanding life, until now, no method has existed to determine the full-length sequence of endogenous mRNAs including their poly(A) tails. Moreover, although non-A nucleotides can be incorporated in poly(A) tails, there also exists no method to accurately sequence them. Here, we present full-length poly(A) and mRNA sequencing (FLAM-seq), a rapid and simple method for high-quality sequencing of entire mRNAs. We report a complementary DNA library preparation method coupled to single-molecule sequencing to perform FLAM-seq. Using human cell lines, brain organoids and Caenorhabditis elegans we show that FLAM-seq delivers high-quality full-length mRNA sequences for thousands of different genes per sample. We find that 3' untranslated region length is correlated with poly(A) tail length, that alternative polyadenylation sites and alternative promoters for the same gene are linked to different tail lengths, and that tails contain a substantial number of cytosines.

Getting the Entire Message: Progress in Isoform Sequencing

The advent of second-generation sequencing and its application to RNA sequencing have revolutionized the field of genomics by allowing quantification of gene expression, as well as the definition of transcription start/end sites, exons, splice sites and RNA editing sites. However, due to the sequencing of fragments of cDNAs, these methods have not given a reliable picture of complete RNA isoforms. Third-generation sequencing has filled this gap and allows end-to-end sequencing of entire RNA/cDNA molecules. This approach to transcriptomics has been a "niche" technology for a couple of years but now is becoming mainstream with many different applications. Here, we review the background and progress made to date in this rapidly growing field. We start by reviewing the progressive realization that alternative splicing is omnipresent. We then focus on long-noncoding RNA isoforms and the distinct combination patterns of exons in noncoding and coding genes. We consider the implications of the recent technologies of direct RNA sequencing and single-cell isoform RNA sequencing. Finally, we discuss the parameters that define the success of long-read RNA sequencing experiments and strategies commonly used to make the most of such data.

FLAM-seq: full-length mRNA sequencing reveals principles of poly(A) tail length control

Although messenger RNAs are key molecules for understanding life, until now, no method has existed to determine the full-length sequence of endogenous mRNAs including their poly(A) tails. Moreover, although non-A nucleotides can be incorporated in poly(A) tails, there also exists no method to accurately sequence them. Here, we present full-length poly(A) and mRNA sequencing (FLAM-seq), a rapid and simple method for high-quality sequencing of entire mRNAs. We report a complementary DNA library preparation method coupled to single-molecule sequencing to perform FLAM-seq. Using human cell lines, brain organoids and Caenorhabditis elegans we show that FLAM-seq delivers high-quality full-length mRNA sequences for thousands of different genes per sample. We find that 3' untranslated region length is correlated with poly(A) tail length, that alternative polyadenylation sites and alternative promoters for the same gene are linked to different tail lengths, and that tails contain a substantial number of cytosines.

RBM47-regulated alternative splicing of TJP1 promotes actin stress fiber assembly during epithelial-to-mesenchymal transition

Morphological and functional changes in cells during the epithelial-mesenchymal transition (EMT) process are known to be regulated by alternative splicing. However, only a few splicing factors involved in EMT have been reported and their underlying mechanisms remain largely unknown. Here, we showed that an isoform of tight junction protein 1 (TJP1) lacking exon 20 (TJP1-α-) is predominantly expressed in tumor tissues and in A549 cells during transforming growth factor-β (TGF-β)-induced EMT. RBM47 promoted the inclusion of exon 20 of TJP1, the alternative exon encoding the α-domain, by which RBM47 recognizes to (U)GCAUG in the downstream intronic region of exon 20. We also found that the first RNA recognition motif (RRM) domain of RBM47 is critical in the regulation of alternative splicing and its recognition to pre-mRNA of TJP1. Furthermore, we demonstrated that the TJP1-α- isoform enhances the assembly of actin stress fibers, thereby promoting cellular migration in a wound healing assay. Our results suggest the regulatory mechanism for the alternative splicing of TJP1 pre-mRNA by RBM47 during EMT, providing a basis for studies related to the modulation of EMT via alternative splicing.

Combinatorial regulation of alternative splicing

The generation of protein coding mRNAs from pre-mRNA is a fundamental biological process that is required for gene expression. Alternative pre-mRNA splicing is responsible for much of the transcriptomic and proteomic diversity observed in higher order eukaryotes. Aberrations that disrupt regular alternative splicing patterns are known to cause human diseases, including various cancers. Alternative splicing is a combinatorial process, meaning many factors affect which two splice sites are ligated together. The features that dictate exon inclusion are comprised of splice site strength, intron-exon architecture, RNA secondary structure, splicing regulatory elements, promoter use and transcription speed by RNA polymerase and the presence of post-transcriptional nucleotide modifications. A comprehensive view of all of the factors that influence alternative splicing decisions is necessary to predict splicing outcomes and to understand the molecular basis of disease. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Regulation of alternative splicing by p300-mediated acetylation of splicing factors

Splicing of precursor mRNA (pre-mRNA) is an important regulatory step in gene expression. Recent evidence points to a regulatory role of chromatin-related proteins in alternative splicing regulation. Using an unbiased approach, we have identified the acetyltransferase p300 as a key chromatin-related regulator of alternative splicing. p300 promotes genome-wide exon inclusion in both a transcription-dependent and -independent manner. Using CD44 as a paradigm, we found that p300 regulates alternative splicing by modulating the binding of splicing factors to pre-mRNA. Using a tethering strategy, we found that binding of p300 to the CD44 promoter region promotes CD44v exon inclusion independently of RNAPII transcriptional elongation rate. Promoter-bound p300 regulates alternative splicing by acetylating splicing factors, leading to exclusion of hnRNP M from CD44 pre-mRNA and activation of Sam68. p300-mediated CD44 alternative splicing reduces cell motility and promotes epithelial features. Our findings reveal a chromatin-related mechanism of alternative splicing regulation and demonstrate its impact on cellular function.

Pervasive Chromatin-RNA Binding Protein Interactions Enable RNA-Based Regulation of Transcription

Increasing evidence suggests that transcriptional control and chromatin activities at large involve regulatory RNAs, which likely enlist specific RNA-binding proteins (RBPs). Although multiple RBPs have been implicated in transcription control, it has remained unclear how extensively RBPs directly act on chromatin. We embarked on a large-scale RBP ChIP-seq analysis, revealing widespread RBP presence in active chromatin regions in the human genome. Like transcription factors (TFs), RBPs also show strong preference for hotspots in the genome, particularly gene promoters, where their association is frequently linked to transcriptional output. Unsupervised clustering reveals extensive co-association between TFs and RBPs, as exemplified by YY1, a known RNA-dependent TF, and RBM25, an RBP involved in splicing regulation. Remarkably, RBM25 depletion attenuates all YY1-dependent activities, including chromatin binding, DNA looping, and transcription. We propose that various RBPs may enhance network interaction through harnessing regulatory RNAs to control transcription.

Evolutionary Divergence of Duplicated Hsf Genes in Populus

Heat shock transcription factors (Hsfs), which function as the activator of heat shock proteins (Hsps), play multiple roles in response to environmental stress and the development of plants. The Hsf family had experienced gene expansion via whole-genome duplication from a single cell algae to higher plants. However, how the Hsf gene family went through evolutionary divergence after genome duplication is unknown. As a model wood species, Populus trichocarpa is widely distributed in North America with various ecological and climatic environments. In this study, we used P. trichocarpa as materials and identified the expression divergence of the PtHsf gene family in developmental processes, such as dormant bud formation and opening, catkins development, and in response to environments. Through the co-expression network, we further discovered the divergent co-expressed genes that related to the functional divergence of PtHsfs. Then, we studied the alternative splicing events, single nucleotide polymorphism distribution and tertiary structures of members of the PtHsf gene family. In addition to expression divergence, we uncovered the evolutionary divergence in the protein level which may be important to new function formations and for survival in changing environments. This study comprehensively analyzed the evolutionary divergence of a member of the PtHsf gene family after genome duplication, paving the way for further gene function analysis and genetic engineering.

Pre-mRNA Splicing in the Nuclear Landscape

Eukaryotic gene expression requires the cumulative activity of multiple molecular machines to synthesize and process newly transcribed pre-messenger RNA. Introns, the noncoding regions in pre-mRNA, must be removed by the spliceosome, which assembles on the pre-mRNA as it is transcribed by RNA polymerase II (Pol II). The assembly and activity of the spliceosome can be modulated by features including the speed of transcription elongation, chromatin, post-translational modifications of Pol II and histone tails, and other RNA processing events like 5'-end capping. Here, we review recent work that has revealed cooperation and coordination among co-transcriptional processing events and speculate on new avenues of research. We anticipate new mechanistic insights capable of unraveling the relative contribution of coupled processing to gene expression.

Alternative Splicing and Transcription Elongation in Plants

Alternative splicing and transcription elongation by RNA polymerase II (RNAPII) are two processes which are tightly connected. Splicing is a co-transcriptional process, and different experimental approaches show that splicing is coupled to transcription in Drosophila, yeast and mammals. However, little is known about coupling of transcription and alternative splicing in plants. The kinetic coupling explains how changes in RNAPII elongation rate influence alternative splicing choices. Recent work in Arabidopsis shows that expression of a dominant negative transcription elongation factor, TFIIS, enhances exon inclusion. Furthermore, the Arabidopsis transcription elongation complex has been recently described, providing new information about elongation factors that interact with elongating RNAPII. Light regulates alternative splicing in plants through a chloroplast retrograde signaling. We have recently shown that light promotes RNAPII elongation in the affected genes, while in darkness elongation is lower. These changes in transcription are consistent with elongation causing the observed changes in alternative splicing. Altogether, these findings provide evidence that coupling between transcription and alternative splicing is an important layer of gene expression regulation in plants.

Epigenetic regulation of alternative splicing

Alternative splicing (AS) serves as an additional regulatory process for gene expression after transcription, and it generates distinct mRNA species, and even noncoding RNAs (ncRNAs), from one primary transcript. Generally, AS can be coupled with transcription and subjected to epigenetic regulation, such as DNA methylation and histone modifications. In addition, ncRNAs, especially long noncoding RNAs (lncRNAs), can be generated from AS and function as splicing factors ("interactors" or "hijackers") in AS. Recently, RNA modifications, such as the RNA N6-methyladenosine (m6A) modification, have been found to regulate AS. In this review, we summarize recent achievements related to the epigenetic regulation of AS.

Beyond Transcription: Roles of Transcription Factors in Pre-mRNA Splicing

Whereas individual steps of protein-coding gene expression in eukaryotes can be studied in isolation in vitro, it has become clear that these steps are intimately connected within cells. Connections not only ensure quality control but also fine-tune the gene expression process, which must adapt to environmental changes while remaining robust. In this review, we systematically present proven and potential mechanisms by which sequence-specific DNA-binding transcription factors can alter gene expression beyond transcription initiation and regulate pre-mRNA splicing, and thereby mRNA isoform production, by (i) influencing transcription elongation rates, (ii) binding to pre-mRNA to recruit splicing factors, and/or (iii) blocking the association of splicing factors with pre-mRNA. We propose various mechanistic models throughout the review, in some cases without explicit supportive evidence, in hopes of providing fertile ground for future studies.

Microfluidic isoform sequencing shows widespread splicing coordination in the human transcriptome

Understanding transcriptome complexity is crucial for understanding human biology and disease. Technologies such as Synthetic long-read RNA sequencing (SLR-RNA-seq) delivered 5 million isoforms and allowed assessing splicing coordination. Pacific Biosciences and Oxford Nanopore increase throughput also but require high input amounts or amplification. Our new droplet-based method, sparse isoform sequencing (spISO-seq), sequences 100k–200k partitions of 10–200 molecules at a time, enabling analysis of 10–100 million RNA molecules. SpISO-seq requires less than 1 ng of input cDNA, limiting or removing the need for prior amplification with its associated biases. Adjusting the number of reads devoted to each molecule reduces sequencing lanes and cost, with little loss in detection power. The increased number of molecules expands our understanding of isoform complexity. In addition to confirming our previously published cases of splicing coordination (e.g., BIN1), the greater depth reveals many new cases, such as MAPT. Coordination of internal exons is found to be extensive among protein coding genes: 23.5%–59.3% (95% confidence interval) of highly expressed genes with distant alternative exons exhibit coordination, showcasing the need for long-read transcriptomics. However, coordination is less frequent for noncoding sequences, suggesting a larger role of splicing coordination in shaping proteins. Groups of genes with coordination are involved in protein–protein interactions with each other, raising the possibility that coordination facilitates complex formation and/or function. We also find new splicing coordination types, involving initial and terminal exons. Our results provide a more comprehensive understanding of the human transcriptome and a general, cost-effective method to analyze it.

Alberto Kornblihtt: Coupling alternative splicing with transcription

Kornblihtt studies how gene expression is controlled in mammals and plants.

A Broad Set of Chromatin Factors Influences Splicing

Several studies propose an influence of chromatin on pre-mRNA splicing, but it is still unclear how widespread and how direct this phenomenon is. We find here that when assembled in vivo, the U2 snRNP co-purifies with a subset of chromatin-proteins, including histones and remodeling complexes like SWI/SNF. Yet, an unbiased RNAi screen revealed that the outcome of splicing is influenced by a much larger variety of chromatin factors not all associating with the spliceosome. The availability of this broad range of chromatin factors impacting splicing further unveiled their very context specific effect, resulting in either inclusion or skipping, depending on the exon under scrutiny. Finally, a direct assessment of the impact of chromatin on splicing using an in vitro co-transcriptional splicing assay with pre-mRNAs transcribed from a nucleosomal template, demonstrated that chromatin impacts nascent pre-mRNP in their competence for splicing. Altogether, our data show that numerous chromatin factors associated or not with the spliceosome can affect the outcome of splicing, possibly as a function of the local chromatin environment that by default interferes with the efficiency of splicing.

Splicing is an RNA editing step allowing to produce multiple transcripts from a single gene. The gene itself is organized in chromatin, associating DNA and multiple proteins. Some proteins regulating the compaction of the chromatin also affect RNA splicing. Yet, it was unclear whether these chromatin proteins were exceptions or whether chromatin very generally affected the outcome of splicing. Here, we show that a subset of chromatin proteins is physically in interaction with the enzyme responsible for RNA splicing. In addition, several chromatin proteins not found directly associated with the splicing machinery were also able to influence RNA splicing, suggesting that chromatin compaction very globally plays a role in splicing. This finding was confirmed using the first in vitro assay combining transcription and splicing in the context of chromatin; this assay showed that assembling DNA with chromatin proteins influences the efficiency of splicing.

Connections Underlying Translation and mRNA Stability

Gene expression and regulation in organisms minimally depends on transcription by RNA polymerase and on the stability of the RNA product (for both coding and non-coding RNAs). For coding RNAs, gene expression is further influenced by the amount of translation by the ribosome and by the stability of the protein product. The stabilities of these two classes of RNA, non-coding and coding, vary considerably: tRNAs and rRNAs tend to be long lived while mRNAs tend to be more short lived. Even among mRNAs, however, there is a considerable range in stability (ranging from seconds to hours in bacteria and up to days in metazoans), suggesting a significant role for stability in the regulation of gene expression. Here, we review recent experiments from bacteria, yeast and metazoans indicating that the stability of most mRNAs is broadly impacted by the actions of ribosomes that translate them. Ribosomal recognition of defective mRNAs triggers "mRNA surveillance" pathways that target the mRNA for degradation [Shoemaker and Green (2012) ]. More generally, even the stability of perfectly functional mRNAs appears to be dictated by overall rates of translation by the ribosome [Herrick et al. (1990), Presnyak et al. (2015) ]. Given that mRNAs are synthesized for the purpose of being translated into proteins, it is reassuring that such intimate connections between mRNA and the ribosome can drive biological regulation. In closing, we consider the likelihood that these connections between protein synthesis and mRNA stability are widespread or whether other modes of regulation dominate the mRNA stability landscape in higher organisms.

Characterization of the Regulation of CD46 RNA Alternative Splicing

Here we present a detailed analysis of the alternative splicing regulation of human CD46, which generates different isoforms with distinct functions. CD46 is a ubiquitous membrane protein that protects host cells from complement and plays other roles in immunity, autophagy, and cell adhesion. CD46 deficiency causes an autoimmune disorder, and this protein is also involved in pathogen infection and cancer. Before this study, the mechanisms of CD46 alternative splicing remained unexplored even though dysregulation of this process has been associated with autoimmune diseases. We proved that the 5' splice sites of CD46 cassette exons 7 and 8 encoding extracellular domains are defined by noncanonical mechanisms of base pairing to U1 small nuclear RNA. Next we characterized the regulation of CD46 cassette exon 13, whose inclusion or skipping generates different cytoplasmic tails with distinct functions. Using splicing minigenes, we identified multiple exonic and intronic splicing enhancers and silencers that regulate exon 13 inclusion via trans-acting splicing factors like PTBP1 and TIAL1. Interestingly, a common splicing activator such as SRSF1 appears to repress CD46 exon 13 inclusion. We also report that expression of CD46 mRNA isoforms is further regulated by non-sense-mediated mRNA decay and transcription speed. Finally, we successfully manipulated CD46 exon 13 inclusion using antisense oligonucleotides, opening up opportunities for functional studies of the isoforms as well as for therapeutics for autoimmune diseases. This study provides insight into CD46 alternative splicing regulation with implications for its function in the immune system and for genetic disease.

Coupling of RNA Polymerase II Transcription Elongation with Pre-mRNA Splicing

Pre-mRNA maturation frequently occurs at the same time and place as transcription by RNA polymerase II. The co-transcriptionality of mRNA processing has permitted the evolution of mechanisms that functionally couple transcription elongation with diverse events that occur on the nascent RNA. This review summarizes the current understanding of the relationship between transcriptional elongation through a chromatin template and co-transcriptional splicing including alternative splicing decisions that affect the expression of most human genes.