Promoter directionality is controlled by U1 snRNP and polyadenylation signals

Characterizing nascent transcription patterns of PROMPTs, eRNAs, and readthrough transcripts in the ENCODE4 deeply profiled cell lines

Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.

U1 snRNP Biogenesis Defects in Neurodegenerative Diseases

The U1 small ribonucleoprotein (U1 snRNP) plays a pivotal role in the intricate process of gene expression, specifically within nuclear RNA processing. By initiating the splicing reaction and modulating 3'-end processing, U1 snRNP exerts precise control over RNA metabolism and gene expression. This ribonucleoparticle is abundantly present, and its complex biogenesis necessitates shuttling between the nuclear and cytoplasmic compartments. Over the past three decades, extensive research has illuminated the crucial connection between disrupted U snRNP biogenesis and several prominent human diseases, notably various neurodegenerative conditions. The perturbation of U1 snRNP homeostasis has been firmly established in diseases such as Spinal Muscular Atrophy, Pontocerebellar hypoplasia, and FUS-mediated Amyotrophic Lateral Sclerosis. Intriguingly, compelling evidence suggests a potential correlation in Fronto-temporal dementia and Alzheimer's disease as well. Although the U snRNP biogenesis pathway is conserved across all eukaryotic cells, neurons, in particular, appear to be highly susceptible to alterations in spliceosome homeostasis. In contrast, other cell types exhibit a greater resilience to such disturbances. This vulnerability underscores the intricate relationship between U1 snRNP dynamics and the health of neuronal cells, shedding light on potential avenues for understanding and addressing neurodegenerative disorders.

Sequence basis of transcription initiation in the human genome

Transcription initiation is a process that is essential to ensuring the proper function of any gene, yet we still lack a unified understanding of sequence patterns and rules that explain most transcription start sites in the human genome. By predicting transcription initiation at base-pair resolution from sequences with a deep learning-inspired explainable model called Puffin, we show that a small set of simple rules can explain transcription initiation at most human promoters. We identify key sequence patterns that contribute to human promoter activity, each activating transcription with distinct position-specific effects. Furthermore, we explain the sequence basis of bidirectional transcription at promoters, identify the links between promoter sequence and gene expression variation across cell types, and explore the conservation of sequence determinants of transcription initiation across mammalian species.

Transcription directionality is licensed by Integrator at active human promoters

A universal characteristic of eukaryotic transcription is that the promoter recruits RNA polymerase II (RNAPII) to produce both precursor mRNAs (pre-mRNAs) and short unstable promoter upstream transcripts (PROMPTs) toward the opposite direction. However, how the transcription machinery selects the correct direction to produce pre-mRNAs is largely unknown. Here, through multiple acute auxin-inducible degradation systems, we show that rapid depletion of an RNAPII-binding protein complex, Integrator, results in robust PROMPT accumulation throughout the genome. Interestingly, the accumulation of PROMPTs is compensated by the reduction of pre-mRNA transcripts in actively transcribed genes. Consistently, Integrator depletion alters the distribution of polymerase between the sense and antisense directions, which is marked by increased RNAPII-carboxy-terminal domain Tyr1 phosphorylation at PROMPT regions and a reduced Ser2 phosphorylation level at transcription start sites. Mechanistically, the endonuclease activity of Integrator is critical to suppress PROMPT production. Furthermore, our data indicate that the presence of U1 binding sites on nascent transcripts could counteract the cleavage activity of Integrator. In this process, the absence of robust U1 signal at most PROMPTs allows Integrator to suppress the antisense transcription and shift the transcriptional balance in favor of the sense direction.

Widespread impact of nucleosome remodelers on transcription at cis-regulatory elements

Nucleosome remodeling complexes and other regulatory factors work in concert to build a chromatin environment that directs the expression of a distinct set of genes in each cell using cis-regulatory elements (CREs), such as promoters and enhancers, that drive transcription of both mRNAs and CRE-associated non-coding RNAs (ncRNAs). Two classes of CRE-associated ncRNAs include upstream antisense RNAs (uaRNAs), which are transcribed divergently from a shared mRNA promoter, and enhancer RNAs (eRNAs), which are transcribed bidirectionally from active enhancers. The complicated network of CRE regulation by nucleosome remodelers remains only partially explored, with a focus on a select, limited number of remodelers. We endeavored to elucidate a remodeler-based regulatory network governing CRE-associated transcription (mRNA, eRNA, and uaRNA) in murine embryonic stem (ES) cells to test the hypothesis that many SNF2-family nucleosome remodelers collaborate to regulate the coding and non-coding transcriptome via alteration of underlying nucleosome architecture. Using depletion followed by transient transcriptome sequencing (TT-seq), we identified thousands of misregulated mRNAs and CRE-associated ncRNAs across the remodelers examined, identifying novel contributions by understudied remodelers in the regulation of coding and noncoding transcription. Our findings suggest that mRNA and eRNA transcription are coordinately co-regulated, while mRNA and uaRNAs sharing a common promoter are independently regulated. Subsequent mechanistic studies suggest that while remodelers SRCAP and CHD8 modulate transcription through classical mechanisms such as transcription factors and histone variants, a broad set of remodelers including SMARCAL1 indirectly contribute to transcriptional regulation through maintenance of genomic stability and proper Integrator complex localization. This study systematically examines the contribution of SNF2-remodelers to the CRE-associated transcriptome, identifying at least two classes for remodeler action.

mRNA nuclear export: how mRNA identity features distinguish functional RNAs from junk transcripts

The division of the cellular space into nucleoplasm and cytoplasm promotes quality control mechanisms that prevent misprocessed mRNAs and junk RNAs from gaining access to the translational machinery. Here, we explore how properly processed mRNAs are distinguished from both misprocessed mRNAs and junk RNAs by the presence or absence of various 'identity features'.

Restrictor synergizes with Symplekin and PNUTS to terminate extragenic transcription

Transcription termination pathways mitigate the detrimental consequences of unscheduled promiscuous initiation occurring at hundreds of thousands of genomic cis-regulatory elements. The Restrictor complex, composed of the Pol II-interacting protein WDR82 and the RNA-binding protein ZC3H4, suppresses processive transcription at thousands of extragenic sites in mammalian genomes. Restrictor-driven termination does not involve nascent RNA cleavage, and its interplay with other termination machineries is unclear. Here we show that efficient termination at Restrictor-controlled extragenic transcription units involves the recruitment of the protein phosphatase 1 (PP1) regulatory subunit PNUTS, a negative regulator of the SPT5 elongation factor, and Symplekin, a protein associated with RNA cleavage complexes but also involved in cleavage-independent and phosphatase-dependent termination of noncoding RNAs in yeast. PNUTS and Symplekin act synergistically with, but independently from, Restrictor to dampen processive extragenic transcription. Moreover, the presence of limiting nuclear levels of Symplekin imposes a competition for its recruitment among multiple transcription termination machineries, resulting in mutual regulatory interactions. Hence, by synergizing with Restrictor, Symplekin and PNUTS enable efficient termination of processive, long-range extragenic transcription.

Selected humanization of yeast U1 snRNP leads to global suppression of pre-mRNA splicing and mitochondrial dysfunction in the budding yeast

The recognition of 5' splice site (5' ss) is one of the earliest steps of pre-mRNA splicing. To better understand the mechanism and regulation of 5' ss recognition, we selectively humanized components of the yeast U1 snRNP to reveal the function of these components in 5' ss recognition and splicing. We targeted U1C and Luc7, two proteins that interact with and stabilize the yeast U1 (yU1) snRNA and the 5' ss RNA duplex. We replaced the Zinc-Finger (ZnF) domain of yU1C with its human counterpart, which resulted in cold-sensitive growth phenotype and moderate splicing defects. Next, we added an auxin-inducible degron to yLuc7 protein and found that Luc7-depleted yU1 snRNP resulted in the concomitant loss of PRP40 and Snu71 (two other essential yeast U1 snRNP proteins), and further biochemical analyses suggest a model of how these three proteins interact with each other in the U1 snRNP. The loss of these proteins resulted in a significant growth retardation accompanied by a global suppression of pre-mRNA splicing. The splicing suppression led to mitochondrial dysfunction as revealed by a release of Fe 2+ into the growth medium and an induction of mitochondrial reactive oxygen species. Together, these observations indicate that the human U1C ZnF can substitute that of yeast, Luc7 is essential for the incorporation of the Luc7-Prp40-Snu71 trimer into yeast U1 snRNP, and splicing plays a major role in the regulation of mitochondria function in yeast.

Premature transcription termination complex proteins PCF11 and WDR82 silence HIV-1 expression in latently infected cells

Postintegration transcriptional silencing of HIV-1 leads to the establishment of a pool of latently infected cells. In these cells, mechanisms controlling RNA Polymerase II (RNAPII) pausing and premature transcription termination (PTT) remain to be explored. Here, we found that the cleavage and polyadenylation (CPA) factor PCF11 represses HIV-1 expression independently of the other subunits of the CPA complex or the polyadenylation signal located at the 5' LTR. We show that PCF11 interacts with the RNAPII-binding protein WDR82. Knock-down of PCF11 or WDR82 reactivated HIV-1 expression in latently infected cells. To silence HIV-1 transcription, PCF11 and WDR82 are specifically recruited at the promoter-proximal region of the provirus in an interdependent manner. Codepletion of PCF11 and WDR82 indicated that they act on the same pathway to repress HIV expression. These findings reveal PCF11/WDR82 as a PTT complex silencing HIV-1 expression in latently infected cells.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Dual agonistic and antagonistic roles of ZC3H18 provide for co-activation of distinct nuclear RNA decay pathways

The RNA exosome is a versatile ribonuclease. In the nucleoplasm of mammalian cells, it is assisted by its adaptors the nuclear exosome targeting (NEXT) complex and the poly(A) exosome targeting (PAXT) connection. Via its association with the ARS2 and ZC3H18 proteins, NEXT/exosome is recruited to capped and short unadenylated transcripts. Conversely, PAXT/exosome is considered to target longer and adenylated substrates via their poly(A) tails. Here, mutational analysis of the core PAXT component ZFC3H1 uncovers a separate branch of the PAXT pathway, which targets short adenylated RNAs and relies on a direct ARS2-ZFC3H1 interaction. We further demonstrate that similar acidic-rich short linear motifs of ZFC3H1 and ZC3H18 compete for a common ARS2 epitope. Consequently, while promoting NEXT function, ZC3H18 antagonizes PAXT activity. We suggest that this organization of RNA decay complexes provides co-activation of NEXT and PAXT at loci with abundant production of short exosome substrates.

Long non-coding RNA generated from CDKN1A gene by alternative polyadenylation regulates p21 expression during DNA damage response

Alternative Polyadenylation (APA) is an emerging mechanism for dynamic changes in gene expression. Previously, we described widespread APA occurrence in introns during the DNA damage response (DDR). Here, we show that a DDR-activated APA event occurs in the first intron of CDKN1A, inducing an alternate last exon-containing lncRNA. We named this lncRNA SPUD (Selective Polyadenylation Upon DNA Damage). SPUD localizes to polysomes in the cytoplasm and is detectable as multiple isoforms in available high-throughput studies. SPUD has low abundance compared to the CDKN1A full-length isoform under non-stress conditions, and SPUD is induced in cancer and normal cells under a variety of DNA damaging conditions in part through p53. The RNA binding protein HuR binds to and promotes the stability of SPUD precursor RNA. SPUD induction increases p21 protein, but not mRNA levels, affecting p21 functions in cell-cycle, CDK2 expression and cell growth. Like CDKN1A full-length isoform, SPUD can bind two competitive p21 translational regulators, the inhibitor calreticulin and the activator CUGBP1; SPUD alters their association with CDKN1A full-length in a DDR-dependent manner, promoting CDKN1A translation. Together, these results show a new regulatory mechanism by which a lncRNA controls p21 expression post-transcriptionally, highlighting lncRNA relevance in DDR progression and cell-cycle.

PAPγ associates with PAXT nuclear exosome to control the abundance of PROMPT ncRNAs

Pervasive transcription of the human genome generates an abundance of RNAs that must be processed and degraded. The nuclear RNA exosome is the main RNA degradation machinery in the nucleus. However, nuclear exosome must be recruited to its substrates by targeting complexes, such as NEXT or PAXT. By proteomic analysis, we identify additional subunits of PAXT, including many orthologs of MTREC found in S. pombe. In particular, we show that polyA polymerase gamma (PAPγ) associates with PAXT. Genome-wide mapping of the binding sites of ZFC3H1, RBM27 and PAPγ shows that PAXT is recruited to the TSS of hundreds of genes. Loss of ZFC3H1 abolishes recruitment of PAXT subunits including PAPγ to TSSs and concomitantly increases the abundance of PROMPTs at the same sites. Moreover, PAPγ, as well as MTR4 and ZFC3H1, is implicated in the polyadenylation of PROMPTs. Our results thus provide key insights into the direct targeting of PROMPT ncRNAs by PAXT at their genomic sites.

Genome-wide measurement of RNA dissociation from chromatin classifies transcripts by their dynamics and reveals rapid dissociation of enhancer lncRNAs

Long non-coding RNAs (lncRNAs) are involved in gene expression regulation in cis. Although enriched in the cell chromatin fraction, to what degree this defines their regulatory potential remains unclear. Furthermore, the factors underlying lncRNA chromatin tethering, as well as the molecular basis of efficient lncRNA chromatin dissociation and its impact on enhancer activity and target gene expression, remain to be resolved. Here, we developed chrTT-seq, which combines the pulse-chase metabolic labeling of nascent RNA with chromatin fractionation and transient transcriptome sequencing to follow nascent RNA transcripts from their transcription on chromatin to release and allows the quantification of dissociation dynamics. By incorporating genomic, transcriptomic, and epigenetic metrics, as well as RNA-binding protein propensities, in machine learning models, we identify features that define transcript groups of different chromatin dissociation dynamics. Notably, lncRNAs transcribed from enhancers display reduced chromatin retention, suggesting that, in addition to splicing, their chromatin dissociation may shape enhancer activity.

Enhancers associated with unstable RNAs are rare in plants

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating ("nascent") transcripts across diverse crops and other plants using capped small (cs)RNA-seq. We discovered that unstable transcripts are rare, unlike in vertebrates, and often originate from promoters. Additionally, many "distal" elements in plants initiate tissue-specific stable transcripts and are likely bone fide promoters of yet-unannotated genes or non-coding RNAs, cautioning against using genome annotations to infer "enhancers" or transcript stability. To investigate enhancer function, we integrated STARR-seq data. We found that annotated promoters, and other regions that initiate stable transcripts rather than unstable transcripts, function as stronger enhancers in plants. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements encompass diverse structures and mechanisms in eukaryotes.

Promoters and introns as key drivers for enhanced gene expression in Saccharomyces cerevisiae

The transcription of genes in the yeast Saccharomyces cerevisiae is governed by multiple layers of regulatory elements and proteins, cooperating to ensure optimum expression of the final protein product based on the cellular requirements. Promoters have always been regarded as the most important determinant of gene transcription, but introns also play a key role in the expression of intron-encoding genes. Some introns can enhance transcription when introduced either promoter-proximal or embedded in the open reading frame of genes. However, the outcome is seldom predictable, with some introns increasing or decreasing transcription depending on the promoter and reporter gene employed. This chapter provides an overview of the general structure and function of promoters and introns and how they may cooperate during transcription to allow intron-mediated enhancement of gene expression. Since S. cerevisiae is a suitable host for recombinant protein production on a commercial level, stronger and more controllable promoters are in high demand. Enhanced gene expression can be achieved via promoter engineering, which may include introns that increase the efficacy of recombinant expression cassettes. Different models for the role of introns in transcription are briefly discussed to show how these intervening sequences can actively interact with the transcription machinery. Furthermore, recent examples of improved protein production via the introduction of promoter-proximal introns are highlighted to showcase the potential value of intron-mediated enhancement of gene expression.

Sequence basis of transcription initiation in human genome

Transcription initiation is an essential process for ensuring proper function of any gene, however, a unified understanding of sequence patterns and rules that determine transcription initiation sites in human genome remains elusive. By explaining transcription initiation at basepair resolution from sequence with a deep learning-inspired explainable modeling approach, here we show that simple rules can explain the vast majority of human promoters. We identified key sequence patterns that contribute to human promoter function, each activating transcription with a distinct position-specific effect curve that likely reflects its mechanism of promoting transcription initiation. Most of these position-specific effects have not been previously characterized, and we verified them using experimental perturbations of transcription factors and sequences. We revealed the sequence basis of bidirectional transcription at promoters and links between promoter selectivity and gene expression variation across cell types. Additionally, by analyzing 241 mammalian genomes and mouse transcription initiation site data, we showed that the sequence determinants are conserved across mammalian species. Taken together, we provide a unified model of the sequence basis of transcription initiation at the basepair level that is broadly applicable across mammalian species, and shed new light on basic questions related to promoter sequence and function.

Emerging concepts involving inhibitory and activating RNA functionalization towards the understanding of microcephaly phenotypes and brain diseases in humans

Microcephaly is characterized as a small head circumference, and is often accompanied by developmental disorders. Several candidate risk genes for this disease have been described, and mutations in non-coding regions are occasionally found in patients with microcephaly. Various non-coding RNAs (ncRNAs), such as microRNAs (miRNAs), SINEUPs, telomerase RNA component (TERC), and promoter-associated lncRNAs (pancRNAs) are now being characterized. These ncRNAs regulate gene expression, enzyme activity, telomere length, and chromatin structure through RNA binding proteins (RBPs)-RNA interaction. Elucidating the potential roles of ncRNA-protein coordination in microcephaly pathogenesis might contribute to its prevention or recovery. Here, we introduce several syndromes whose clinical features include microcephaly. In particular, we focus on syndromes for which ncRNAs or genes that interact with ncRNAs may play roles. We discuss the possibility that the huge ncRNA field will provide possible new therapeutic approaches for microcephaly and also reveal clues about the factors enabling the evolutionary acquisition of the human-specific "large brain."

A restrictor complex of ZC3H4, WDR82, and ARS2 integrates with PNUTS to control unproductive transcription

The transcriptional termination of unstable non-coding RNAs (ncRNAs) is poorly understood compared to coding transcripts. We recently identified ZC3H4-WDR82 ("restrictor") as restricting human ncRNA transcription, but how it does this is unknown. Here, we show that ZC3H4 additionally associates with ARS2 and the nuclear exosome targeting complex. The domains of ZC3H4 that contact ARS2 and WDR82 are required for ncRNA restriction, suggesting their presence in a functional complex. Consistently, ZC3H4, WDR82, and ARS2 co-transcriptionally control an overlapping population of ncRNAs. ZC3H4 is proximal to the negative elongation factor, PNUTS, which we show enables restrictor function and is required to terminate the transcription of all major RNA polymerase II transcript classes. In contrast to short ncRNAs, longer protein-coding transcription is supported by U1 snRNA, which shields transcripts from restrictor and PNUTS at hundreds of genes. These data provide important insights into the mechanism and control of transcription by restrictor and PNUTS.

U1 snRNP increases RNA Pol II elongation rate to enable synthesis of long genes

The expansion of introns within mammalian genomes poses a challenge for the production of full-length messenger RNAs (mRNAs), with increasing evidence that these long AT-rich sequences present obstacles to transcription. Here, we investigate RNA polymerase II (RNAPII) elongation at high resolution in mammalian cells and demonstrate that RNAPII transcribes faster across introns. Moreover, we find that this acceleration requires the association of U1 snRNP (U1) with the elongation complex at 5' splice sites. The role of U1 to stimulate elongation rate through introns reduces the frequency of both premature termination and transcriptional arrest, thereby dramatically increasing RNA production. We further show that changes in RNAPII elongation rate due to AT content and U1 binding explain previous reports of pausing or termination at splice junctions and the edge of CpG islands. We propose that U1-mediated acceleration of elongation has evolved to mitigate the risks that long AT-rich introns pose to transcript completion.

Structural basis of a transcription pre-initiation complex on a divergent promoter

Most eukaryotic promoter regions are divergently transcribed. As the RNA polymerase II pre-initiation complex (PIC) is intrinsically asymmetric and responsible for transcription in a single direction, it is unknown how divergent transcription arises. Here, the Saccharomyces cerevisiae Mediator complexed with a PIC (Med-PIC) was assembled on a divergent promoter and analyzed by cryoelectron microscopy. The structure reveals two distinct Med-PICs forming a dimer through the Mediator tail module, induced by a homodimeric activator protein localized near the dimerization interface. The tail dimer is associated with ∼80-bp upstream DNA, such that two flanking core promoter regions are positioned and oriented in a suitable form for PIC assembly in opposite directions. Also, cryoelectron tomography visualized the progress of the PIC assembly on the two core promoter regions, providing direct evidence for the role of the Med-PIC dimer in divergent transcription.

A CpG island-encoded mechanism protects genes from premature transcription termination

Transcription must be tightly controlled to regulate gene expression and development. However, our understanding of the molecular mechanisms that influence transcription and how these are coordinated in cells to ensure normal gene expression remains rudimentary. Here, by dissecting the function of the SET1 chromatin-modifying complexes that bind to CpG island-associated gene promoters, we discover that they play a specific and essential role in enabling the expression of low to moderately transcribed genes. Counterintuitively, this effect can occur independently of SET1 complex histone-modifying activity and instead relies on an interaction with the RNA Polymerase II-binding protein WDR82. Unexpectedly, we discover that SET1 complexes enable gene expression by antagonising premature transcription termination by the ZC3H4/WDR82 complex at CpG island-associated genes. In contrast, at extragenic sites of transcription, which typically lack CpG islands and SET1 complex occupancy, we show that the activity of ZC3H4/WDR82 is unopposed. Therefore, we reveal a gene regulatory mechanism whereby CpG islands are bound by a protein complex that specifically protects genic transcripts from premature termination, effectively distinguishing genic from extragenic transcription and enabling normal gene expression.

The emerging roles of long non-coding RNA in host immune response and intracellular bacterial infections

The long non-coding RNAs (lncRNAs) are evolutionarily conserved classes of non-coding regulatory transcripts of > 200 nucleotides in length. They modulate several transcriptional and post-transcriptional events in the organism. Depending on their cellular localization and interactions, they regulate chromatin function and assembly; and alter the stability and translation of cytoplasmic mRNAs. Although their proposed range of functionality remains controversial, there is increasing research evidence that lncRNAs play a regulatory role in the activation, differentiation and development of immune signaling cascades; microbiome development; and in diseases such as neuronal and cardiovascular disorders; cancer; and pathogenic infections. This review discusses the functional roles of different lncRNAs in regulation of host immune responses, signaling pathways during host-microbe interaction and infection caused by obligate intracellular bacterial pathogens. The study of lncRNAs is assuming significance as it could be exploited for development of alternative therapeutic strategies for the treatment of severe and chronic pathogenic infections caused by Mycobacterium, Chlamydia and Rickettsia infections, as well as commensal colonization. Finally, this review summarizes the translational potential of lncRNA research in development of diagnostic and prognostic tools for human diseases.

De novo genes with an lncRNA origin encode unique human brain developmental functionality

Human de novo genes can originate from neutral long non-coding RNA (lncRNA) loci and are evolutionarily significant in general, yet how and why this all-or-nothing transition to functionality happens remains unclear. Here, in 74 human/hominoid-specific de novo genes, we identified distinctive U1 elements and RNA splice-related sequences accounting for RNA nuclear export, differentiating mRNAs from lncRNAs, and driving the origin of de novo genes from lncRNA loci. The polymorphic sites facilitating the lncRNA-mRNA conversion through regulating nuclear export are selectively constrained, maintaining a boundary that differentiates mRNAs from lncRNAs. The functional new genes actively passing through it thus showed a mode of pre-adaptive origin, in that they acquire functions along with the achievement of their coding potential. As a proof of concept, we verified the regulations of splicing and U1 recognition on the nuclear export efficiency of one of these genes, the ENSG00000205704, in human neural progenitor cells. Notably, knock-out or over-expression of this gene in human embryonic stem cells accelerates or delays the neuronal maturation of cortical organoids, respectively. The transgenic mice with ectopically expressed ENSG00000205704 showed enlarged brains with cortical expansion. We thus demonstrate the key roles of nuclear export in de novo gene origin. These newly originated genes should reflect the novel uniqueness of human brain development.

MPP6 stimulates both RRP6 and DIS3 to degrade a specified subset of MTR4-sensitive substrates in the human nucleus

Recent in vitro reconstitution analyses have proven that the physical interaction between the exosome core and MTR4 helicase, which promotes the exosome activity, is maintained by either MPP6 or RRP6. However, knowledge regarding the function of MPP6 with respect to in vivo exosome activity remains scarce. Here, we demonstrate a facilitative function of MPP6 that composes a specific part of MTR4-dependent substrate decay by the human exosome. Using RNA polymerase II-transcribed poly(A)+ substrate accumulation as an indicator of a perturbed exosome, we found functional redundancy between RRP6 and MPP6 in the decay of these poly(A)+ transcripts. MTR4 binding to the exosome core via MPP6 was essential for MPP6 to exert its redundancy with RRP6. However, at least for the decay of our identified exosome substrates, MTR4 recruitment by MPP6 was not functionally equivalent to recruitment by RRP6. Genome-wide classification of substrates based on their sensitivity to each exosome component revealed that MPP6 deals with a specific range of substrates and highlights the importance of MTR4 for their decay. Considering recent findings of competitive binding to the exosome between auxiliary complexes, our results suggest that the MPP6-incorporated MTR4-exosome complex is one of the multiple alternative complexes rather than the prevailing one.

Computational identification of signals predictive for nuclear RNA exosome degradation pathway targeting

The RNA exosome degrades transcripts in the nucleoplasm of mammalian cells. Its substrate specificity is mediated by two adaptors: the ‘nuclear exosome targeting (NEXT)’ complex and the ‘poly(A) exosome targeting (PAXT)’ connection. Previous studies have revealed some DNA/RNA elements that differ between the two pathways, but how informative these features are for distinguishing pathway targeting, or whether additional genomic features that are informative for such classifications exist, is unknown. Here, we leverage the wealth of available genomic data and develop machine learning models that predict exosome targets and subsequently rank the features the models use by their predictive power. As expected, features around transcript end sites were most predictive; specifically, the lack of canonical 3′ end processing was highly predictive of NEXT targets. Other associated features, such as promoter-proximal G/C content and 5′ splice sites, were informative, but only for distinguishing NEXT and not PAXT targets. Finally, we discovered predictive features not previously associated with exosome targeting, in particular RNA helicase DDX3X binding sites. Overall, our results demonstrate that nucleoplasmic exosome targeting is to a large degree predictable, and our approach can assess the predictive power of previously known and new features in an unbiased way.

Pervasive transcription: a controlled risk

Transcriptome-wide interrogation of eukaryotic genomes has unveiled the pervasive nature of RNA polymerase II transcription. Virtually, any DNA region with an accessible chromatin structure can be transcribed, resulting in a mass production of noncoding RNAs (ncRNAs) with the potential of interfering with gene expression programs. Budding yeast has proved to be a powerful model organism to understand the mechanisms at play to control pervasive transcription and overcome the risks of hazardous disruption of cellular functions. In this review, we focus on the actors and strategies yeasts employ to govern ncRNA production, and we discuss recent findings highlighting the dangers of losing control over pervasive transcription.

ZFC3H1 and U1-70K promote the nuclear retention of mRNAs with 5' splice site motifs within nuclear speckles

Quality control of mRNA represents an important regulatory mechanism for gene expression in eukaryotes. One component of this quality control is the nuclear retention and decay of misprocessed RNAs. Previously, we demonstrated that mature mRNAs containing a 5' splice site (5'SS) motif, which is typically found in misprocessed RNAs such as intronic polyadenylated (IPA) transcripts, are nuclear retained and degraded. Using high-throughput sequencing of cellular fractions, we now demonstrate that IPA transcripts require the zinc finger protein ZFC3H1 for their nuclear retention and degradation. Using reporter mRNAs, we demonstrate that ZFC3H1 promotes the nuclear retention of mRNAs with intact 5'SS motifs by sequestering them into nuclear speckles. Furthermore, we find that U1-70K, a component of the spliceosomal U1 snRNP, is also required for the nuclear retention of these reporter mRNAs and likely functions in the same pathway as ZFC3H1. Finally, we show that the disassembly of nuclear speckles impairs the nuclear retention of reporter mRNAs with 5'SS motifs. Our results highlight a splicing independent role of U1 snRNP and indicate that it works in conjunction with ZFC3H1 in preventing the nuclear export of misprocessed mRNAs by sequestering them into nuclear speckles.

Screening thousands of transcribed coding and non-coding regions reveals sequence determinants of RNA polymerase II elongation potential

Precise regulation of transcription by RNA polymerase II (RNAPII) is critical for organismal growth and development. However, what determines whether an engaged RNAPII will synthesize a full-length transcript or terminate prematurely is poorly understood. Notably, RNAPII is far more susceptible to termination when transcribing non-coding RNAs than when synthesizing protein-coding mRNAs, but the mechanisms underlying this are unclear. To investigate the impact of transcribed sequence on elongation potential, we developed a method to screen the effects of thousands of INtegrated Sequences on Expression of RNA and Translation using high-throughput sequencing (INSERT-seq). We found that higher AT content in non-coding RNAs, rather than specific sequence motifs, drives RNAPII termination. Further, we demonstrate that 5' splice sites autonomously stimulate processive transcription, even in the absence of polyadenylation signals. Our results reveal a potent role for the transcribed sequence in dictating gene output and demonstrate the power of INSERT-seq toward illuminating these contributions.

Inflammatory Immune-Associated eRNA: Mechanisms, Functions and Therapeutic Prospects

The rapid development of multiple high-throughput sequencing technologies has made it possible to explore the critical roles and mechanisms of functional enhancers and enhancer RNAs (eRNAs). The inflammatory immune response, as a fundamental pathological process in infectious diseases, cancers and immune disorders, coordinates the balance between the internal and external environment of the organism. It has been shown that both active enhancers and intranuclear eRNAs are preferentially expressed over inflammation-related genes in response to inflammatory stimuli, suggesting that enhancer transcription events and their products influence the expression and function of inflammatory genes. Therefore, in this review, we summarize and discuss the relevant inflammatory roles and regulatory mechanisms of eRNAs in inflammatory immune cells, non-inflammatory immune cells, inflammatory immune diseases and tumors, and explore the potential therapeutic effects of enhancer inhibitors affecting eRNA production for diseases with inflammatory immune responses.

Promoter antisense RNAs: beyond transcription by-products of active promoters

Omics-based technologies, which have developed rapidly over the last few decades, have generated increasing evidence demonstrating pervasive divergent transcription from RNA polymerase II (Pol II) promoters of eukaryotic genome, and indeed have raised considerable discussion as to their potential physiopathological function. Unlike many other long non-coding RNAs (lncRNAs), promoter antisense RNAs (PAS RNAs) were initially considered to be merely passive transcription by-products of active promoters. However, recent studies have begun to reveal their critical importance in a broad spectrum of biological processes. In this Review, I summarize recent technological advances that enable accurate detection of PAS RNA and discuss the mechanisms of PAS RNA biogenesis emphasizing the functional importance of its structure enabling the diverse functions of PAS RNA in transcription and chromatin regulation.

Control of non-productive RNA polymerase II transcription via its early termination in metazoans

Transcription establishes the universal first step of gene expression where RNA is produced by a DNA-dependent RNA polymerase. The most versatile of eukaryotic RNA polymerases, RNA polymerase II (Pol II), transcribes a broad range of DNA including protein-coding and a variety of non-coding transcription units. Although Pol II can be configured as a durable enzyme capable of transcribing hundreds of kilobases, there is reliable evidence of widespread abortive Pol II transcription termination shortly after initiation, which is often followed by rapid degradation of the associated RNA. The molecular details underlying this phenomenon are still vague but likely reflect the action of quality control mechanisms on the early Pol II complex. Here, we summarize current knowledge of how and when such promoter-proximal quality control is asserted on metazoan Pol II.

The activity of human enhancers is modulated by the splicing of their associated lncRNAs

Pervasive enhancer transcription is at the origin of more than half of all long noncoding RNAs in humans. Transcription of enhancer-associated long noncoding RNAs (elncRNA) contribute to their cognate enhancer activity and gene expression regulation in cis. Recently, splicing of elncRNAs was shown to be associated with elevated enhancer activity. However, whether splicing of elncRNA transcripts is a mere consequence of accessibility at highly active enhancers or if elncRNA splicing directly impacts enhancer function, remains unanswered. We analysed genetically driven changes in elncRNA splicing, in humans, to address this outstanding question. We showed that splicing related motifs within multi-exonic elncRNAs evolved under selective constraints during human evolution, suggesting the processing of these transcripts is unlikely to have resulted from transcription across spurious splice sites. Using a genome-wide and unbiased approach, we used nucleotide variants as independent genetic factors to directly assess the causal relationship that underpin elncRNA splicing and their cognate enhancer activity. We found that the splicing of most elncRNAs is associated with changes in chromatin signatures at cognate enhancers and target mRNA expression. We provide evidence that efficient and conserved processing of enhancer-associated elncRNAs contributes to enhancer activity.

Most, if not all, active enhancers are transcribed, giving rise to a plethora of transcripts, including enhancer-associated long noncoding RNAs (elncRNAs). Changes in elncRNA levels impacts cognate enhancer activity. Recently splicing of elncRNA has also been found to associate with enhancer activity. Whether this associations reflects a contribution of elncRNA splicing to increased enhancer activity or else is simply the consequence of increased chromatin accessibility that promotes transcriptional elongation and allows for spurious splicing events remains unknown. We show that natural selection has acted, at the species and population level, to preserve DNA elements required for frequent and efficient elncRNA splicing Importantly, using a genome-wide and unbiased statistical population genomics approach, we demonstrate that elncRNA splicing is associated with cognate enhancer function, contributing to chromatin status and enhancer activity. Our results provides strong evidence that efficient elncRNA splicing contributes to enhancer activity genome-wide.

Nuclear RNA Exosome and Pervasive Transcription: Dual Sculptors of Genome Function

The genome is pervasively transcribed across various species, yielding numerous non-coding RNAs. As a counterbalance for pervasive transcription, various organisms have a nuclear RNA exosome complex, whose structure is well conserved between yeast and mammalian cells. The RNA exosome not only regulates the processing of stable RNA species, such as rRNAs, tRNAs, small nucleolar RNAs, and small nuclear RNAs, but also plays a central role in RNA surveillance by degrading many unstable RNAs and misprocessed pre-mRNAs. In addition, associated cofactors of RNA exosome direct the exosome to distinct classes of RNA substrates, suggesting divergent and/or multi-layer control of RNA quality in the cell. While the RNA exosome is essential for cell viability and influences various cellular processes, mutations and alterations in the RNA exosome components are linked to the collection of rare diseases and various diseases including cancer, respectively. The present review summarizes the relationships between pervasive transcription and RNA exosome, including evolutionary crosstalk, mechanisms of RNA exosome-mediated RNA surveillance, and physiopathological effects of perturbation of RNA exosome.

The Hda1 histone deacetylase limits divergent non-coding transcription and restricts transcription initiation frequency

Nucleosome-depleted regions (NDRs) at gene promoters support initiation of RNA polymerase II transcription. Interestingly, transcription often initiates in both directions, resulting in an mRNA and a divergent non-coding (DNC) transcript of unclear purpose. Here, we characterized the genetic architecture and molecular mechanism of DNC transcription in budding yeast. Using high-throughput reverse genetic screens based on quantitative single-cell fluorescence measurements, we identified the Hda1 histone deacetylase complex (Hda1C) as a repressor of DNC transcription. Nascent transcription profiling showed a genome-wide role of Hda1C in repression of DNC transcription. Live-cell imaging of transcription revealed that mutations in the Hda3 subunit increased the frequency of DNC transcription. Hda1C contributed to decreased acetylation of histone H3 in DNC transcription regions, supporting DNC transcription repression by histone deacetylation. Our data support the interpretation that DNC transcription results as a consequence of the NDR-based architecture of eukaryotic promoters, but that it is governed by locus-specific repression to maintain genome fidelity.

Beyond the canonical role of TFIIB in eukaryotic transcription

The role of general transcription factor TFIIB in transcription extends well beyond its evolutionarily conserved function in initiation. Chromatin localization studies demonstrating binding of TFIIB to both the 5’ and 3’ ends of genes in a diverse set of eukaryotes strongly suggested a rather unexpected role of the factor in termination. TFIIB indeed plays a role in termination of transcription. TFIIB occupancy of the 3’ end is possibly due to its interaction with the termination factors residing there. Interaction of the promoter-bound TFIIB with factors occupying the 3’ end of a gene may be the basis of transcription-dependent gene looping. The proximity of the terminator-bound factors with the promoter in a gene loop has the potential to terminate promoter-initiated upstream anti-sense transcription thereby conferring promoter directionality. TFIIB, therefore, is emerging as a factor with pleiotropic roles in the transcription cycle. This could be the reason for preferential targeting of TFIIB by viruses. Further studies are needed to understand the critical role of TFIIB in viral pathogenesis in the context of its newly identified roles in termination, gene looping and promoter directionality.

Heterogeneity among enhancer RNAs: origins, consequences and perspectives

Enhancer RNAs (eRNAs) are non-coding RNAs transcribed from distal cis-regulatory elements (i.e. enhancers), which are stereotyped as short, rarely spliced and unstable. In fact, a non-negligible fraction of eRNAs seems to be longer, spliced and more stable, and their cognate enhancers are epigenomically and functionally distinguishable from typical enhancers. In this review, we first summarized the genomic and molecular origins underlying the observed heterogeneity among eRNAs. Then, we discussed how their heterogeneous properties (e.g. stability) affect the modes of interaction with their regulatory partners, from promiscuous cis-interactions to specific trans-interactions. Finally, we highlighted the existence of a seemingly continuous spectrum of eRNA properties and its implications in the genomic origins of non-coding RNA genes from an evolutionary perspective.

The evolutionary acquisition and mode of functions of promoter-associated non-coding RNAs (pancRNAs) for mammalian development

Increasing evidence has shown that many long non-coding RNAs (lncRNAs) are involved in gene regulation in a variety of ways such as transcriptional, post-transcriptional and epigenetic regulation. Promoter-associated non-coding RNAs (pancRNAs), which are categorized into the most abundant single-copy lncRNA biotype, play vital regulatory roles in finely tuning cellular specification at the epigenomic level. In short, pancRNAs can directly or indirectly regulate downstream genes to participate in the development of organisms in a cell-specific manner. In this review, we will introduce the evolutionarily acquired characteristics of pancRNAs as determined by comparative epigenomics and elaborate on the research progress on pancRNA-involving processes in mammalian embryonic development, including neural differentiation.

The upstream 5' splice site remains associated to the transcription machinery during intron synthesis

In the earliest step of spliceosome assembly, the two splice sites flanking an intron are brought into proximity by U1 snRNP and U2AF along with other proteins. The mechanism that facilitates this intron looping is poorly understood. Using a CRISPR interference-based approach to halt RNA polymerase II transcription in the middle of introns in human cells, we discovered that the nascent 5′ splice site base pairs with a U1 snRNA that is tethered to RNA polymerase II during intron synthesis. This association functionally corresponds with splicing outcome, involves bona fide 5′ splice sites and cryptic intronic sites, and occurs transcriptome-wide. Overall, our findings reveal that the upstream 5′ splice sites remain attached to the transcriptional machinery during intron synthesis and are thus brought into proximity of the 3′ splice sites; potentially mediating the rapid splicing of long introns.

We know that most splicing reactions take place co-transcriptionally, but how the transcription machinery facilitate splicing of introns is unknown. Here the authors show that the 5′ splice site remains associated with the transcription machinery during intron synthesis through U1 snRNP, providing a basis for the rapid splicing reaction of introns.

An in vitro reconstituted U1 snRNP allows the study of the disordered regions of the particle and the interactions with proteins and ligands

U1 small nuclear ribonucleoparticle (U1 snRNP) plays a central role during RNA processing. Previous structures of U1 snRNP revealed how the ribonucleoparticle is organized and recognizes the pre-mRNA substrate at the exon–intron junction. As with many other ribonucleoparticles involved in RNA metabolism, U1 snRNP contains extensions made of low complexity sequences. Here, we developed a protocol to reconstitute U1 snRNP in vitro using mostly full-length components in order to perform liquid-state NMR spectroscopy. The accuracy of the reconstitution was validated by probing the shape and structure of the particle by SANS and cryo-EM. Using an NMR spectroscopy-based approach, we probed, for the first time, the U1 snRNP tails at atomic detail and our results confirm their high degree of flexibility. We also monitored the labile interaction between the splicing factor PTBP1 and U1 snRNP and validated the U1 snRNA stem loop 4 as a binding site for the splicing regulator on the ribonucleoparticle. Altogether, we developed a method to probe the intrinsically disordered regions of U1 snRNP and map the interactions controlling splicing regulation. This approach could be used to get insights into the molecular mechanisms of alternative splicing and screen for potential RNA therapeutics.

Characterizing Promoter and Enhancer Sequences by a Deep Learning Method

Promoters and enhancers are well-known regulatory elements modulating gene expression. As confirmed by high-throughput sequencing technologies, these regulatory elements are bidirectionally transcribed. That is, promoters produce stable mRNA in the sense direction and unstable RNA in the antisense direction, while enhancers transcribe unstable RNA in both directions. Although it is thought that enhancers and promoters share a similar architecture of transcription start sites (TSSs), how the transcriptional machinery distinctly uses these genomic regions as promoters or enhancers remains unclear. To address this issue, we developed a deep learning (DL) method by utilizing a convolutional neural network (CNN) and the saliency algorithm. In comparison with other classifiers, our CNN presented higher predictive performance, suggesting the overarching importance of the high-order sequence features, captured by the CNN. Moreover, our method revealed that there are substantial sequence differences between the enhancers and promoters. Remarkably, the 20–120 bp downstream regions from the center of bidirectional TSSs seemed to contribute to the RNA stability. These regions in promoters tend to have a larger number of guanines and cytosines compared to those in enhancers, and this feature contributed to the classification of the regulatory elements. Our CNN-based method can capture the complex TSS architectures. We found that the genomic regions around TSSs for promoters and enhancers contribute to RNA stability and show GC-biased characteristics as a critical determinant for promoter TSSs.

High-resolution mapping of function and protein binding in an RNA nuclear enrichment sequence

The functions of long RNAs, including mRNAs and long noncoding RNAs (lncRNAs), critically depend on their subcellular localization. The identity of the sequences that dictate subcellular localization and their high-resolution anatomy remain largely unknown. We used a suite of massively parallel RNA assays and libraries containing thousands of sequence variants to pinpoint the functional features within the SIRLOIN element, which dictates nuclear enrichment through hnRNPK recruitment. In addition, we profiled the endogenous SIRLOIN RNA-nucleoprotein complex and identified the nuclear RNA-binding proteins SLTM and SNRNP70 as novel SIRLOIN binders. Taken together, using massively parallel assays, we identified the features that dictate binding of hnRNPK, SLTM, and SNRNP70 to SIRLOIN and found that these factors are jointly required for SIRLOIN activity. Our study thus provides a roadmap for high-throughput dissection of functional sequence elements in long RNAs.

Discovery of widespread transcription initiation at microsatellites predictable by sequence-based deep neural network

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

Known sequence features can explain half of all human gene ends

Cleavage and polyadenylation (CPA) sites define eukaryotic gene ends. CPA sites are associated with five key sequence recognition elements: the upstream UGUA, the polyadenylation signal (PAS), and U-rich sequences; the CA/UA dinucleotide where cleavage occurs; and GU-rich downstream elements (DSEs). Currently, it is not clear whether these sequences are sufficient to delineate CPA sites. Additionally, numerous other sequences and factors have been described, often in the context of promoting alternative CPA sites and preventing cryptic CPA site usage. Here, we dissect the contributions of individual sequence features to CPA using standard discriminative models. We show that models comprised only of the five primary CPA sequence features give highest probability scores to constitutive CPA sites at the ends of coding genes, relative to the entire pre-mRNA sequence, for 41% of all human genes. U1-hybridizing sequences provide a small boost in performance. The addition of all known RBP RNA binding motifs to the model, however, increases this figure to 49%, and suggests an involvement of both known and suspected CPA regulators as well as potential new factors in delineating constitutive CPA sites. To our knowledge, this high effectiveness of established features to predict human gene ends has not previously been documented.

Comprehensive RNP profiling in cells identifies U1 snRNP complexes with cleavage and polyadenylation factors active in telescripting

Full-length transcription in the majority of protein-coding and other genes transcribed by RNA polymerase II in complex eukaryotes requires U1 snRNP (U1) to co-transcriptionally suppress transcription-terminating premature 3'-end cleavage and polyadenylation (PCPA) from cryptic polyadenylation signals (PASs). This U1 activity, termed telescripting, requires U1 to base-pair with the nascent RNA and inhibit usage of a downstream PAS. Here we describe experimental methods to determine the mechanism of U1 telescripting, involving mapping of U1 and CPA factors (CPAFs) binding locations in relation to PCPA sites, and identify U1 and CPAFs interactomes. The methods which utilizes rapid reversible protein-RNA and protein-protein chemical crosslinking, immunoprecipitations (XLIPs) of components of interest, and RNA-seq and quantitative proteomic mass spectrometry, captured U1-CPAFs complexes in cells, providing important insights into telescripting mechanism. XLIP profiling can be used for comprehensive molecular definition of diverse RNPs.

Mouse totipotent stem cells captured and maintained through spliceosomal repression

Since establishment of the first embryonic stem cells (ESCs), in vitro culture of totipotent cells functionally and molecularly comparable with in vivo blastomeres with embryonic and extraembryonic developmental potential has been a challenge. Here we report that spliceosomal repression in mouse ESCs drives a pluripotent-to-totipotent state transition. Using the splicing inhibitor pladienolide B, we achieve stable in vitro culture of totipotent ESCs comparable at molecular levels with 2- and 4-cell blastomeres, which we call totipotent blastomere-like cells (TBLCs). Mouse chimeric assays combined with single-cell RNA sequencing (scRNA-seq) demonstrate that TBLCs have a robust bidirectional developmental capability to generate multiple embryonic and extraembryonic cell lineages. Mechanically, spliceosomal repression causes widespread splicing inhibition of pluripotent genes, whereas totipotent genes, which contain few short introns, are efficiently spliced and transcriptionally activated. Our study provides a means for capturing and maintaining totipotent stem cells.

Intergenic RNA mainly derives from nascent transcripts of known genes

BACKGROUND

Eukaryotic genomes undergo pervasive transcription, leading to the production of many types of stable and unstable RNAs. Transcription is not restricted to regions with annotated gene features but includes almost any genomic context. Currently, the source and function of most RNAs originating from intergenic regions in the human genome remain unclear.

RESULTS

We hypothesize that many intergenic RNAs can be ascribed to the presence of as-yet unannotated genes or the "fuzzy" transcription of known genes that extends beyond the annotated boundaries. To elucidate the contributions of these two sources, we assemble a dataset of more than 2.5 billion publicly available RNA-seq reads across 5 human cell lines and multiple cellular compartments to annotate transcriptional units in the human genome. About 80% of transcripts from unannotated intergenic regions can be attributed to the fuzzy transcription of existing genes; the remaining transcripts originate mainly from putative long non-coding RNA loci that are rarely spliced. We validate the transcriptional activity of these intergenic RNAs using independent measurements, including transcriptional start sites, chromatin signatures, and genomic occupancies of RNA polymerase II in various phosphorylation states. We also analyze the nuclear localization and sensitivities of intergenic transcripts to nucleases to illustrate that they tend to be rapidly degraded either on-chromatin by XRN2 or off-chromatin by the exosome.

CONCLUSIONS

We provide a curated atlas of intergenic RNAs that distinguishes between alternative processing of well-annotated genes from independent transcriptional units based on the combined analysis of chromatin signatures, nuclear RNA localization, and degradation pathways.

Genome-wide strand asymmetry in massively parallel reporter activity favors genic strands

Massively parallel reporter assays (MPRAs) are useful tools to characterize regulatory elements in human genomes. An aspect of MPRAs that is not typically the focus of analysis is their intrinsic ability to differentiate activity levels for a given sequence element when placed in both of its possible orientations relative to the reporter construct. Here, we describe pervasive strand asymmetry of MPRA signals in data sets from multiple reporter configurations in both published and newly reported data. These effects are reproducible across different cell types and in different treatments within a cell type and are observed both within and outside of annotated regulatory elements. From elements in gene bodies, MPRA strand asymmetry favors the sense strand, suggesting that function related to endogenous transcription is driving the phenomenon. Similarly, we find that within Alu mobile element insertions, strand asymmetry favors the transcribed strand of the ancestral retrotransposon. The effect is consistent across the multiplicity of Alu elements in human genomes and is more pronounced in less diverged Alu elements. We find sequence features driving MPRA strand asymmetry and show its prediction from sequence alone. We see some evidence for RNA stabilization and transcriptional activation mechanisms and hypothesize that the effect is driven by natural selection favoring efficient transcription. Our results indicate that strand asymmetry is a pervasive and reproducible feature in MPRA data. More importantly, the fact that MPRA asymmetry favors naturally transcribed strands suggests that it stems from preserved biological functions that have a substantial, global impact on gene and genome evolution.

A first exon termination checkpoint preferentially suppresses extragenic transcription

Interactions between the splicing machinery and RNA polymerase II increase protein-coding gene transcription. Similarly, exons and splicing signals of enhancer-generated long noncoding RNAs (elncRNAs) augment enhancer activity. However, elncRNAs are inefficiently spliced, suggesting that, compared with protein-coding genes, they contain qualitatively different exons with a limited ability to drive splicing. We show here that the inefficiently spliced first exons of elncRNAs as well as promoter-antisense long noncoding RNAs (pa-lncRNAs) in human and mouse cells trigger a transcription termination checkpoint that requires WDR82, an RNA polymerase II-binding protein, and its RNA-binding partner of previously unknown function, ZC3H4. We propose that the first exons of elncRNAs and pa-lncRNAs are an intrinsic component of a regulatory mechanism that, on the one hand, maximizes the activity of these cis-regulatory elements by recruiting the splicing machinery and, on the other, contains elements that suppress pervasive extragenic transcription.