Transcription elongation rate has a tissue-specific impact on alternative cleavage and polyadenylation in Drosophila melanogaster

Tumor suppressor miR-317 and lncRNA Peony are expressed from a polycistronic non-coding RNA locus that regulates germline differentiation and testis morphology

This research focuses on investigating the impact of non-coding RNAs on stem cell biology and differentiation processes. We found that miR-317 plays a role in germline stem cell progeny differentiation. miR-317 and its neighbor, the lncRNA Peony, originate and are co-expressed from a singular polycistronic non-coding RNA locus. Alternative polyadenylation is implicated in regulation of their differential expression. While the increased expression of the lncRNA Peony results in the disruption of the muscle sheath covering the testis, the absence of miR-317 leads to the emergence of germline tumors in young flies. The deficiency of miR-317 increases Notch signaling activity in the somatic cyst cells, which drives germline tumorigenesis. Germline tumors also arise from upregulation of several predicted targets of miR-317, among which are regulators of the Notch pathway. This implicates miR-317 as a novel tumor suppressor that modulates Notch signaling strength.

Light regulates widespread plant alternative polyadenylation through the chloroplast

Transcription of eukaryotic protein-coding genes generates immature mRNAs that are subjected to a series of processing events, including capping, splicing, cleavage, and polyadenylation (CPA), and chemical modifications of bases. Alternative polyadenylation (APA) greatly contributes to mRNA diversity in the cell. By determining the length of the 3' untranslated region, APA generates transcripts with different regulatory elements, such as miRNA and RBP binding sites, which can influence mRNA stability, turnover, and translation. In the model plant Arabidopsis thaliana, APA is involved in the control of seed dormancy and flowering. In view of the physiological importance of APA in plants, we decided to investigate the effects of light/dark conditions and compare the underlying mechanisms to those elucidated for alternative splicing (AS). We found that light controls APA in approximately 30% of Arabidopsis genes. Similar to AS, the effect of light on APA requires functional chloroplasts, is not affected in mutants of the phytochrome and cryptochrome photoreceptor pathways, and is observed in roots only when the communication with the photosynthetic tissues is not interrupted. Furthermore, mitochondrial and TOR kinase activities are necessary for the effect of light. However, unlike AS, coupling with transcriptional elongation does not seem to be involved since light-dependent APA regulation is neither abolished in mutants of the TFIIS transcript elongation factor nor universally affected by chromatin relaxation caused by histone deacetylase inhibition. Instead, regulation seems to correlate with changes in the abundance of constitutive CPA factors, also mediated by the chloroplast.

Neuronal RNA processing: cross-talk between transcriptional regulation and RNA-binding proteins

In the nervous system, alternative RNA processing is particularly prevalent, which results in the expression of thousands of transcript variants found in no other tissue. Neuron-specific RNA-binding proteins co-transcriptionally regulate alternative splicing, alternative polyadenylation, and RNA editing, thereby shaping the RNA identity of nervous system cells. Recent evidence suggests that interactions between RNA-binding proteins and cis-regulatory elements such as promoters and enhancers play a role in the determination of neuron-specific expression profiles. Here, we discuss possible mechanisms through which transcription and RNA processing cross-talk to generate the uniquely complex neuronal transcriptome, with a focus on alternative 3'-end formation.

Chromatin regulates alternative polyadenylation via the RNA polymerase II elongation rate

The RNA polymerase II (Pol II) elongation rate influences poly(A) site selection, with slow and fast Pol II derivatives causing upstream and downstream shifts, respectively, in poly(A) site utilization. In yeast, depletion of either of the histone chaperones FACT or Spt6 causes an upstream shift of poly(A) site use that strongly resembles the poly(A) profiles of slow Pol II mutant strains. Like slow Pol II mutant strains, FACT- and Spt6-depleted cells exhibit Pol II processivity defects, indicating that both Spt6 and FACT stimulate the Pol II elongation rate. Poly(A) profiles of some genes show atypical downstream shifts; this subset of genes overlaps well for FACT- or Spt6-depleted strains but is different from the atypical genes in Pol II speed mutant strains. In contrast, depletion of histone H3 or H4 causes a downstream shift of poly(A) sites for most genes, indicating that nucleosomes inhibit the Pol II elongation rate in vivo. Thus, chromatin-based control of the Pol II elongation rate is a potential mechanism, distinct from direct effects on the cleavage/polyadenylation machinery, to regulate alternative polyadenylation in response to genetic or environmental changes.

How is polyadenylation restricted to 3'-untranslated regions?

Polyadenylation occurs at numerous sites within 3'-untranslated regions (3'-UTRs) but rarely within coding regions. How does Pol II travel through long coding regions without generating poly(A) sites, yet then permits promiscuous polyadenylation once it reaches the 3'-UTR? The cleavage/polyadenylation (CpA) machinery preferentially associates with 3'-UTRs, but it is unknown how its recruitment is restricted to 3'-UTRs during Pol II elongation. Unlike coding regions, 3'-UTRs have long AT-rich stretches of DNA that may be important for restricting polyadenylation to 3'-UTRs. Recognition of the 3'-UTR could occur at the DNA (AT-rich), RNA (AU-rich), or RNA:DNA hybrid (rU:dA- and/or rA:dT-rich) level. Based on the nucleic acid critical for 3'-UTR recognition, there are three classes of models, not mutually exclusive, for how the CpA machinery is selectively recruited to 3'-UTRs, thereby restricting where polyadenylation occurs: (1) RNA-based models suggest that the CpA complex directly (or indirectly through one or more intermediary proteins) binds long AU-rich stretches that are exposed after Pol II passes through these regions. (2) DNA-based models suggest that the AT-rich sequence affects nucleosome depletion or the elongating Pol II machinery, resulting in dissociation of some elongation factors and subsequent recruitment of the CpA machinery. (3) RNA:DNA hybrid models suggest that preferential destabilization of the Pol II elongation complex at rU:dA- and/or rA:dT-rich duplexes bridging the nucleotide addition and RNA exit sites permits preferential association of the CpA machinery with 3'-UTRs. Experiments to provide evidence for one or more of these models are suggested.

(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.

RNA binding proteins in cardiovascular development and disease

Congenital heart disease (CHD) is the most common birth defect affecting>1.35 million newborn babies worldwide. CHD can lead to prenatal, neonatal, postnatal lethality or life-long cardiac complications. RNA binding protein (RBP) mutations or variants are emerging as contributors to CHDs. RBPs are wizards of gene regulation and are major contributors to mRNA and protein landscape. However, not much is known about RBPs in the developing heart and their contributions to CHD. In this chapter, we will discuss our current knowledge about specific RBPs implicated in CHDs. We are in an exciting era to study RBPs using the currently available and highly successful RNA-based therapies and methodologies. Understanding how RBPs shape the developing heart will unveil their contributions to CHD. Identifying their target RNAs in the embryonic heart will ultimately lead to RNA-based treatments for congenital heart disease.

The association of the RSC remodeler complex with chromatin is influenced by the prefoldin-like Bud27 and determines nucleosome positioning and polyadenylation sites usage in Saccharomyces cerevisiae

The tripartite interaction between the chromatin remodeler complex RSC, RNA polymerase subunit Rpb5 and prefoldin-like Bud27 is necessary for proper RNA pol II elongation. Indeed lack of Bud27 alters this association and affects transcription elongation. This work investigates the consequences of lack of Bud27 on the chromatin association of RSC and RNA pol II, and on nucleosome positioning. Our results demonstrate that RSC binds chromatin in gene bodies and lack of Bud27 alters this association, mainly around polyA sites. This alteration impacts chromatin organization and leads to the accumulation of RNA pol II molecules around polyA sites, likely due to pausing or arrest. Our data suggest that RSC is necessary to maintain chromatin organization around those sites, and any alteration of this organization results in the widespread use of alternative polyA sites. Finally, we also find a similar molecular phenotype that occurs upon TOR inhibition with rapamycin, which suggests that alternative polyadenylation observed upon TOR inhibition is likely Bud27-dependent.

Transcriptome Responses to Different Salinity Conditions in Litoditis marina, Revealed by Long-Read Sequencing

The marine nematode Litoditis marina is widely distributed in intertidal zones around the globe, yet the mechanisms underlying its broad adaptation to salinity remain elusive. In this study, we applied ONT long-read sequencing technology to unravel the transcriptome responses to different salinity conditions in L. marina. Through ONT sequencing under 3‱, 30‱ and 60‱ salinity environments, we obtained 131.78 G clean data and 26,647 non-redundant long-read transcripts, including 6464 novel transcripts. The DEGs obtained from the current ONT lrRNA-seq were highly correlated with those identified in our previously reported Illumina short-read RNA sequencing data. When we compared the 30‱ to the 3‱ salinity condition, we found that GO terms such as oxidoreductase activity, cation transmembrane transport and ion transmembrane transport were shared between the ONT lrRNA-seq and Illumina data. Similarly, GO terms including extracellular space, structural constituents of cuticle, substrate-specific channel activity, ion transport and substrate-specific transmembrane transporter activity were shared between the ONT and Illumina data under 60‱ compared to 30‱ salinity. In addition, we found that 79 genes significantly increased, while 119 genes significantly decreased, as the salinity increased. Furthermore, through the GO enrichment analysis of 214 genes containing DAS, in 30‱ compared to 3‱ salinity, we found that GO terms such as cellular component assembly and coenzyme biosynthetic process were enriched. Additionally, we observed that GO terms such as cellular component assembly and coenzyme biosynthetic process were also enriched in 60‱ compared to 30‱ salinity. Moreover, we found that 86, 125, and 81 genes that contained DAS were also DEGs, in comparisons between 30‱ and 3‱, 60‱ and 30‱, and 60‱ and 3‱ salinity, respectively. In addition, we demonstrated the landscape of alternative polyadenylation in marine nematode under different salinity conditions This report provides several novel insights for the further study of the mechanisms by which euryhalinity formed and evolved, and it might also contribute to the investigation of salinity dynamics induced by global climate change.

Elongation rate of RNA polymerase II affects pausing patterns across 3' UTRs

Yeast mRNAs are polyadenylated at multiple sites in their 3' untranslated regions (3' UTRs), and poly(A) site usage is regulated by the rate of transcriptional elongation by RNA polymerase II (Pol II). Slow Pol II derivatives favor upstream poly(A) sites, and fast Pol II derivatives favor downstream poly(A) sites. Transcriptional elongation and polyadenylation are linked at the nucleotide level, presumably reflecting Pol II dwell time at each residue that influences the level of polyadenylation. Here, we investigate the effect of Pol II elongation rate on pausing patterns and the relationship between Pol II pause sites and poly(A) sites within 3' UTRs. Mutations that affect Pol II elongation rate alter sequence preferences at pause sites within 3' UTRs, and pausing preferences differ between 3' UTRs and coding regions. In addition, sequences immediately flanking the pause sites show preferences that are largely independent of Pol II speed. In wild-type cells, poly(A) sites are preferentially located < 50 nucleotides upstream from Pol II pause sites, but this spatial relationship is diminished in cells harboring Pol II speed mutants. Based on a random forest classifier, Pol II pause sites are modestly predicted by the distance to poly(A) sites but are better predicted by the chromatin landscape in Pol II speed derivatives. Transcriptional regulatory proteins can influence the relationship between Pol II pausing and polyadenylation but in a manner distinct from Pol II elongation rate derivatives. These results indicate a complex relationship between Pol II pausing and polyadenylation.

Elongation roadblocks mediated by dCas9 across human genes modulate transcription and nascent RNA processing

Non-cleaving Cas9 (dCas9) is widely employed to manipulate specific gene loci, often with scant regard for unintended transcriptional effects. We demonstrate here that dCas9 mediates precise RNA polymerase II transcriptional pausing followed by transcription termination and potential alternative polyadenylation. By contrast, alternative splicing is unaffected, likely requiring more sustained alteration to elongation speed. The effect on transcription is orientation specific, with pausing only being induced when dCas9-associated guide RNA anneals to the non-template strand. Targeting the template strand induces minimal effects on transcription elongation and thus provides a neutral approach to recruit dCas9-linked effector domains to specific gene regions. In essence, we evaluate molecular effects of targeting dCas9 to mammalian transcription units. In so doing, we also provide new information on elongation by RNA polymerase II and coupled pre-mRNA processing.

Elevated pre-mRNA 3' end processing activity in cancer cells renders vulnerability to inhibition of cleavage and polyadenylation

Cleavage and polyadenylation (CPA) is responsible for 3' end processing of eukaryotic poly(A)+ RNAs and preludes transcriptional termination. JTE-607, which targets CPSF-73, is the first known CPA inhibitor (CPAi) in mammalian cells. Here we show that JTE-607 perturbs gene expression through both transcriptional readthrough and alternative polyadenylation (APA). Sensitive genes are associated with features similar to those previously identified for PCF11 knockdown, underscoring a unified transcriptomic signature of CPAi. The degree of inhibition of an APA site by JTE-607 correlates with its usage level and, consistently, cells with elevated CPA activities, such as those with induced overexpression of FIP1, display greater transcriptomic disturbances when treated with JTE-607. Moreover, JTE-607 causes S phase crisis and is hence synergistic with inhibitors of DNA damage repair pathways. Together, our data reveal CPA activity and proliferation rate as determinants of CPAi-mediated cell death, raising the possibility of using CPAi as an adjunct therapy to suppress certain cancers.

Sequential Polyadenylation to Enable Alternative mRNA 3' End Formation

In eukaryotic cells, a key RNA processing step to generate mature mRNA is the coupled reaction for cleavage and polyadenylation (CPA) at the 3' end of individual transcripts. Many transcripts are alternatively polyadenylated (APA) to produce mRNAs with different 3' ends that may either alter protein coding sequence (CDS-APA) or create different lengths of 3'UTR (tandem-APA). As the CPA reaction is intimately associated with transcriptional termination, it has been widely assumed that APA is regulated cotranscriptionally. Isoforms terminated at different regions may have distinct RNA stability under different conditions, thus altering the ratio of APA isoforms. Such differential impacts on different isoforms have been considered as post-transcriptional APA, but strictly speaking, this can only be considered "apparent" APA, as the choice is not made during the CPA reaction. Interestingly, a recent study reveals sequential APA as a new mechanism for post-transcriptional APA. This minireview will focus on this new mechanism to provide insights into various documented regulatory paradigms.

Regulation and function of alternative polyadenylation in development and differentiation

Alternative processing of nascent mRNAs is widespread in eukaryotic organisms and greatly impacts the output of gene expression. Specifically, alternative cleavage and polyadenylation (APA) is a co-transcriptional molecular process that switches the polyadenylation site (PAS) at which a nascent mRNA is cleaved, resulting in mRNA isoforms with different 3'UTR length and content. APA can potentially affect mRNA translation efficiency, localization, stability, and mRNA seeded protein-protein interactions. APA naturally occurs during development and cellular differentiation, with around 70% of human genes displaying APA in particular tissues and cell types. For example, neurons tend to express mRNAs with long 3'UTRs due to preferential processing at PASs more distal than other PASs used in other cell types. In addition, changes in APA mark a variety of pathological states, including many types of cancer, in which mRNAs are preferentially cleaved at more proximal PASs, causing expression of mRNA isoforms with short 3'UTRs. Although APA has been widely reported, both the function of APA in development and the mechanisms that regulate the choice of 3'end cut sites in normal and pathogenic conditions are still poorly understood. In this review, we summarize current understanding of how APA is regulated during development and cellular differentiation and how the resulting change in 3'UTR content affects multiple aspects of gene expression. With APA being a widespread phenomenon, the advent of cutting-edge scientific techniques and the pressing need for in-vivo studies, there has never been a better time to delve into the intricate mechanisms of alternative cleavage and polyadenylation.

Nucleotide-level linkage of transcriptional elongation and polyadenylation

Alternative polyadenylation yields many mRNA isoforms whose 3' termini occur disproportionately in clusters within 3' untranslated regions. Previously, we showed that profiles of poly(A) site usage are regulated by the rate of transcriptional elongation by RNA polymerase (Pol) II (Geisberg et al., 2020). Pol II derivatives with slow elongation rates confer an upstream-shifted poly(A) profile, whereas fast Pol II strains confer a downstream-shifted poly(A) profile. Within yeast isoform clusters, these shifts occur steadily from one isoform to the next across nucleotide distances. In contrast, the shift between clusters - from the last isoform of one cluster to the first isoform of the next - is much less pronounced, even over large distances. GC content in a region 13-30 nt downstream from isoform clusters correlates with their sensitivity to Pol II elongation rate. In human cells, the upstream shift caused by a slow Pol II mutant also occurs continuously at single nucleotide resolution within clusters but not between them. Pol II occupancy increases just downstream of poly(A) sites, suggesting a linkage between reduced elongation rate and cluster formation. These observations suggest that (1) Pol II elongation speed affects the nucleotide-level dwell time allowing polyadenylation to occur, (2) poly(A) site clusters are linked to the local elongation rate, and hence do not arise simply by intrinsically imprecise cleavage and polyadenylation of the RNA substrate, (3) DNA sequence elements can affect Pol II elongation and poly(A) profiles, and (4) the cleavage/polyadenylation and Pol II elongation complexes are spatially, and perhaps physically, coupled so that polyadenylation occurs rapidly upon emergence of the nascent RNA from the Pol II elongation complex.

Tubular-specific CDK12 knockout causes a defect in urine concentration due to premature cleavage of the slc12a1 gene

Cyclin-dependent kinase 12 (CDK12) plays a critical role in regulating gene transcription. CDK12 inhibition is a potential anticancer therapeutic strategy. However, several clinical trials have shown that CDK inhibitors might cause renal dysfunction and electrolyte disorders. CDK12 is abundant in renal tubular epithelial cells (RTECs), but the exact role of CDK12 in renal physiology remains unclear. Genetic knockout of CDK12 in mouse RTECs causes polydipsia, polyuria, and hydronephrosis. This phenotype is caused by defects in water reabsorption that are the result of reduced Na-K-2Cl cotransporter 2 (NKCC2) levels in the kidney. In addition, CKD12 knockout causes an increase in Slc12a1 (which encodes NKCC2) intronic polyadenylation events, which results in Slc12a1 truncated transcript production and NKCC2 downregulation. These findings provide novel insight into CDK12 being necessary for maintaining renal homeostasis by regulating NKCC2 transcription, which explains the critical water and electrolyte disturbance that occurs during the application of CDK12 inhibitors for cancer treatment. Therefore, there are safety concerns about the clinical use of these new anticancer drugs.

Alternative polyadenylation regulation: insights from sequential polyadenylation

The processing of the proximal and distal poly(A) sites in alternative polyadenylation (APA) has long been thought to independently occur on pre-mRNAs during transcription. However, a recent study by our groups demonstrated that the proximal sites for many genes could be activated sequentially following the distal ones, suggesting a multi-cleavage-same-transcript mode beyond the canonical one-cleavage-per-transcript view. Here, we review the established mechanisms for APA regulation and then discuss the additional insights into APA regulation from the perspective of sequential polyadenylation, resulting in a unified leverage model for understanding the mechanisms of regulated APA.

Precise gene models using long-read sequencing reveal a unique poly(A) signal in Giardia lamblia

During pre-mRNA processing, the poly(A) signal is recognized by a protein complex that ensures precise cleavage and polyadenylation of the nascent transcript. The location of this cleavage event establishes the length and sequence of the 3' UTR of an mRNA, thus determining much of its post-transcriptional fate. Using long-read sequencing, we characterize the polyadenylation signal and related sequences surrounding Giardia lamblia cleavage sites for over 2600 genes. We find that G. lamblia uses an AGURAA poly(A) signal, which differs from the mammalian AAUAAA. We also describe how G. lamblia lacks common auxiliary elements found in other eukaryotes, along with the proteins that recognize them. Further, we identify 133 genes with evidence of alternative polyadenylation. These results suggest that despite pared-down cleavage and polyadenylation machinery, 3' end formation still appears to be an important regulatory step for gene expression in G. lamblia.

Identification of Alternative Polyadenylation in Cyanidioschyzon merolae Through Long-Read Sequencing of mRNA

Alternative polyadenylation (APA) is widespread among metazoans and has been shown to have important impacts on mRNA stability and protein expression. Beyond a handful of well-studied organisms, however, its existence and consequences have not been well investigated. We therefore turned to the deep-branching red alga, Cyanidioschyzon merolae, to study the biology of polyadenylation in an organism highly diverged from humans and yeast. C. merolae is an acidothermophilic alga that lives in volcanic hot springs. It has a highly reduced genome (16.5 Mbp) and has lost all but 27 of its introns and much of its splicing machinery, suggesting that it has been under substantial pressure to simplify its RNA processing pathways. We used long-read sequencing to assess the key features of C. merolae mRNAs, including splicing status and polyadenylation cleavage site (PAS) usage. Splicing appears to be less efficient in C. merolae compared with yeast, flies, and mammalian cells. A high proportion of transcripts (63%) have at least two distinct PAS's, and 34% appear to utilize three or more sites. The apparent polyadenylation signal UAAA is used in more than 90% of cases, in cells grown in both rich media or limiting nitrogen. Our documentation of APA for the first time in this non-model organism highlights its conservation and likely biological importance of this regulatory step in gene expression.

Oncogenic dysregulation of pre-mRNA processing by protein kinases: challenges and therapeutic opportunities

Alternative splicing and polyadenylation represent two major steps in pre-mRNA-processing, which ensure proper gene expression and diversification of human transcriptomes. Deregulation of these processes contributes to oncogenic programmes involved in the onset, progression and evolution of human cancers, which often result in the acquisition of resistance to existing therapies. On the other hand, cancer cells frequently increase their transcriptional rate and develop a transcriptional addiction, which imposes a high stress on the pre-mRNA-processing machinery and establishes a therapeutically exploitable vulnerability. A prominent role in fine-tuning pre-mRNA-processing mechanisms is played by three main families of protein kinases: serine arginine protein kinase (SRPK), CDC-like kinase (CLK) and cyclin-dependent kinase (CDK). These kinases phosphorylate the RNA polymerase, splicing factors and regulatory proteins involved in cleavage and polyadenylation of the nascent transcripts. The activity of SRPKs, CLKs and CDKs can be altered in cancer cells, and their inhibition was shown to exert anticancer effects. In this review, we describe key findings that have been reported on these topics and discuss challenges and opportunities of developing therapeutic approaches targeting splicing factor kinases.

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

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

LABRAT reveals association of alternative polyadenylation with transcript localization, RNA binding protein expression, transcription speed, and cancer survival

BACKGROUND

The sequence content of the 3' UTRs of many mRNA transcripts is regulated through alternative polyadenylation (APA). The study of this process using RNAseq data, though, has been historically challenging.

RESULTS

To combat this problem, we developed LABRAT, an APA isoform quantification method. LABRAT takes advantage of newly developed transcriptome quantification techniques to accurately determine relative APA site usage and how it varies across conditions. Using LABRAT, we found consistent relationships between gene-distal APA and subcellular RNA localization in multiple cell types. We also observed connections between transcription speed and APA site choice as well as tumor-specific transcriptome-wide shifts in APA isoform abundance in hundreds of patient-derived tumor samples that were associated with patient prognosis. We investigated the effects of APA on transcript expression and found a weak overall relationship, although many individual genes showed strong correlations between relative APA isoform abundance and overall gene expression. We interrogated the roles of 191 RNA-binding proteins in the regulation of APA isoforms, finding that dozens promote broad, directional shifts in relative APA isoform abundance both in vitro and in patient-derived samples. Finally, we find that APA site shifts in the two classes of APA, tandem UTRs and alternative last exons, are strongly correlated across many contexts, suggesting that they are coregulated.

CONCLUSIONS

We conclude that LABRAT has the ability to accurately quantify APA isoform ratios from RNAseq data across a variety of sample types. Further, LABRAT is able to derive biologically meaningful insights that connect APA isoform regulation to cellular and molecular phenotypes.

Real-time single-cell characterization of the eukaryotic transcription cycle reveals correlations between RNA initiation, elongation, and cleavage

The eukaryotic transcription cycle consists of three main steps: initiation, elongation, and cleavage of the nascent RNA transcript. Although each of these steps can be regulated as well as coupled with each other, their in vivo dissection has remained challenging because available experimental readouts lack sufficient spatiotemporal resolution to separate the contributions from each of these steps. Here, we describe a novel application of Bayesian inference techniques to simultaneously infer the effective parameters of the transcription cycle in real time and at the single-cell level using a two-color MS2/PP7 reporter gene and the developing fruit fly embryo as a case study. Our method enables detailed investigations into cell-to-cell variability in transcription-cycle parameters as well as single-cell correlations between these parameters. These measurements, combined with theoretical modeling, suggest a substantial variability in the elongation rate of individual RNA polymerase molecules. We further illustrate the power of this technique by uncovering a novel mechanistic connection between RNA polymerase density and nascent RNA cleavage efficiency. Thus, our approach makes it possible to shed light on the regulatory mechanisms in play during each step of the transcription cycle in individual, living cells at high spatiotemporal resolution.

On the function and relevance of alternative 3'-UTRs in gene expression regulation

Messanger RNA (mRNA) isoforms with alternative 3'-untranslated regions (3'-UTRs) are produced by alternative polyadenylation (APA), which occurs during transcription in most eukaryotic genes. APA fine-tunes gene expression in a cell-type- and cellular state-dependent manner. Selection of an APA site entails the binding of core cleavage and polyadenylation factors to a particular polyadenylation site localized in the pre-mRNA and is controlled by multiple regulatory determinants, including transcription, pre-mRNA cis-regulatory sequences, and protein factors. Alternative 3'-UTRs serve as platforms for specific RNA binding proteins and microRNAs, which regulate gene expression in a coordinated manner by controlling mRNA fate and function in the cell. Genome-wide studies illustrated the full extent of APA prevalence and revealed that specific 3'-UTR profiles are associated with particular cellular states and diseases. Generally, short 3'-UTRs are associated with proliferative and cancer cells, and long 3'-UTRs are mostly found in polarized and differentiated cells. Fundamental new insights on the physiological consequences of this widespread event and the molecular mechanisms involved have been revealed through single-cell studies. Publicly available comprehensive databases that cover all APA mRNA isoforms identified in many cellular states and diseases reveal specific APA signatures. Therapies tackling APA mRNA isoforms or APA regulators may be regarded as innovative and attractive tools for diagnostics or treatment of several pathologies. We highlight the function of APA and alternative 3'-UTRs in gene expression regulation, the control of these mechanisms, their physiological consequences, and their potential use as new biomarkers and therapeutic tools. This article is categorized under: RNA Processing > 3' End Processing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA in Disease and Development > RNA in Disease.

FlyBase: updates to the Drosophila melanogaster knowledge base

FlyBase (flybase.org) is an essential online database for researchers using Drosophila melanogaster as a model organism, facilitating access to a diverse array of information that includes genetic, molecular, genomic and reagent resources. Here, we describe the introduction of several new features at FlyBase, including Pathway Reports, paralog information, disease models based on orthology, customizable tables within reports and overview displays ('ribbons') of expression and disease data. We also describe a variety of recent important updates, including incorporation of a developmental proteome, upgrades to the GAL4 search tab, additional Experimental Tool Reports, migration to JBrowse for genome browsing and improvements to batch queries/downloads and the Fast-Track Your Paper tool.

A missense mutation in the CSTF2 gene that impairs the function of the RNA recognition motif and causes defects in 3' end processing is associated with intellectual disability in humans

CSTF2 encodes an RNA-binding protein that is essential for mRNA cleavage and polyadenylation (C/P). No disease-associated mutations have been described for this gene. Here, we report a mutation in the RNA recognition motif (RRM) of CSTF2 that changes an aspartic acid at position 50 to alanine (p.D50A), resulting in intellectual disability in male patients. In mice, this mutation was sufficient to alter polyadenylation sites in over 1300 genes critical for brain development. Using a reporter gene assay, we demonstrated that C/P efficiency of CSTF2^D50A was lower than wild type. To account for this, we determined that p.D50A changed locations of amino acid side chains altering RNA binding sites in the RRM. The changes modified the electrostatic potential of the RRM leading to a greater affinity for RNA. These results highlight the significance of 3′ end mRNA processing in expression of genes important for brain plasticity and neuronal development.

Cold-induced chromatin compaction and nuclear retention of clock mRNAs resets the circadian rhythm

Cooling patients to sub‐physiological temperatures is an integral part of modern medicine. We show that cold exposure induces temperature‐specific changes to the higher‐order chromatin and gene expression profiles of human cells. These changes are particularly dramatic at 18°C, a temperature synonymous with that experienced by patients undergoing controlled deep hypothermia during surgery. Cells exposed to 18°C exhibit largely nuclear‐restricted transcriptome changes. These include the nuclear accumulation of mRNAs encoding components of the negative limbs of the core circadian clock, most notably REV‐ERBα. This response is accompanied by compaction of higher‐order chromatin and hindrance of mRNPs from engaging nuclear pores. Rewarming reverses chromatin compaction and releases the transcripts into the cytoplasm, triggering a pulse of negative limb gene proteins that reset the circadian clock. We show that cold‐induced upregulation of REV‐ERBα is sufficient to trigger this reset. Our findings uncover principles of the cellular cold response that must be considered for current and future applications involving therapeutic deep hypothermia.

So close, no matter how far: multiple paths connecting transcription to mRNA translation in eukaryotes

Transcription of DNA into mRNA and translation of mRNA into proteins are two major processes underlying gene expression. Due to the distinct molecular mechanisms, timings, and locales of action, these processes are mainly considered to be independent. During the last two decades, however, multiple factors and elements were shown to coordinate transcription and translation, suggesting an intricate level of synchronization. This review discusses the molecular mechanisms that impact both processes in eukaryotic cells of different origins. The emerging global picture suggests evolutionarily conserved regulation and coordination between transcription and mRNA translation, indicating the importance of this phenomenon for the fine-tuning of gene expression and the adjustment to constantly changing conditions.

The transcriptional elongation rate regulates alternative polyadenylation in yeast

Yeast cells undergoing the diauxic response show a striking upstream shift in poly(A) site utilization, with increased use of ORF-proximal poly(A) sites resulting in shorter 3' mRNA isoforms for most genes. This altered poly(A) pattern is extremely similar to that observed in cells containing Pol II derivatives with slow elongation rates. Conversely, cells containing derivatives with fast elongation rates show a subtle downstream shift in poly(A) sites. Polyadenylation patterns of many genes are sensitive to both fast and slow elongation rates, and a global shift of poly(A) utilization is strongly linked to increased purine content of sequences flanking poly(A) sites. Pol II processivity is impaired in diauxic cells, but strains with reduced processivity and normal Pol II elongation rates have normal polyadenylation profiles. Thus, Pol II elongation speed is important for poly(A) site selection and for regulating poly(A) patterns in response to environmental conditions.

SETX (senataxin), the helicase mutated in AOA2 and ALS4, functions in autophagy regulation

SETX (senataxin) is an RNA/DNA helicase that has been implicated in transcriptional regulation and the DNA damage response through resolution of R-loop structures. Mutations in SETX result in either of two distinct neurodegenerative disorders. SETX dominant mutations result in a juvenile form of amyotrophic lateral sclerosis (ALS) called ALS4, whereas recessive mutations are responsible for ataxia called ataxia with oculomotor apraxia type 2 (AOA2). How mutations in the same protein can lead to different phenotypes is still unclear. To elucidate AOA2 disease mechanisms, we first examined gene expression changes following SETX depletion. We observed the effects on both transcription and RNA processing, but surprisingly observed decreased R-loop accumulation in SETX-depleted cells. Importantly, we discovered a strong connection between SETX and the macroautophagy/autophagy pathway, reflecting a direct effect on transcription of autophagy genes. We show that SETX depletion inhibits the progression of autophagy, leading to an accumulation of ubiquitinated proteins, decreased ability to clear protein aggregates, as well as mitochondrial defects. Analysis of AOA2 patient fibroblasts also revealed a perturbation of the autophagy pathway. Our work has thus identified a novel function for SETX in the regulation of autophagy, whose modulation may have a therapeutic impact for AOA2.Abbreviations: 3'READS: 3' region extraction and deep sequencing; ACTB: actin beta; ALS4: amyotrophic lateral sclerosis type 4; AOA2: ataxia with oculomotor apraxia type 2; APA: alternative polyadenylation; AS: alternative splicing; ATG7: autophagy-related 7; ATP6V0D2: ATPase H+ transporting V0 subunit D2; BAF: bafilomycin A₁; BECN1: beclin 1; ChIP: chromatin IP; Chloro: chloroquine; CPT: camptothecin; DDR: DNA damage response; DNMT1: DNA methyltransferase 1; DRIP: DNA/RNA IP; DSBs: double strand breaks; EBs: embryoid bodies; FTD: frontotemporal dementia; GABARAP: GABA type A receptor-associated protein; GO: gene ontology; HR: homologous recombination; HTT: huntingtin; IF: immunofluorescence; IP: immunoprecipitation; iPSCs: induced pluripotent stem cells; KD: knockdown; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; MN: motor neuron; MTORC1: mechanistic target of rapamycin kinase complex 1; PASS: PolyA Site Supporting; PFA: paraformaldehyde; RNAPII: RNA polymerase II; SCA: spinocerebellar ataxia; SETX: senataxin; SMA: spinal muscular atrophy; SMN1: survival of motor neuron 1, telomeric; SQSTM1/p62: sequestosome 1; TFEB: transcription factor EB; TSS: transcription start site; TTS: transcription termination site; ULK1: unc-51 like autophagy activating kinase 1; WB: western blot; WIPI2: WD repeat domain, phosphoinositide interacting 2; XRN2: 5'-3' exoribonuclease 2.

Nutrient-dependent control of RNA polymerase II elongation rate regulates specific gene expression programs by alternative polyadenylation

In this study from Yague-Sanz et al., the authors investigated the physiological relevance of variations in RNAPII elongation kinetics, and show in yeast that a RNAPII mutant that reduces the transcription elongation rate causes widespread changes in alternative polyadenylation (APA). Their findings indicate that RNAPII is a sensor of nucleotide availability and that genes important for nucleotide pool maintenance have adopted regulatory mechanisms responsive to reduced rates of transcription elongation.

Transcription by RNA polymerase II (RNAPII) is a dynamic process with frequent variations in the elongation rate. However, the physiological relevance of variations in RNAPII elongation kinetics has remained unclear. Here we show in yeast that a RNAPII mutant that reduces the transcription elongation rate causes widespread changes in alternative polyadenylation (APA). We unveil two mechanisms by which APA affects gene expression in the slow mutant: 3′ UTR shortening and gene derepression by premature transcription termination of upstream interfering noncoding RNAs. Strikingly, the genes affected by these mechanisms are enriched for functions involved in phosphate uptake and purine synthesis, processes essential for maintenance of the intracellular nucleotide pool. As nucleotide concentration regulates transcription elongation, our findings argue that RNAPII is a sensor of nucleotide availability and that genes important for nucleotide pool maintenance have adopted regulatory mechanisms responsive to reduced rates of transcription elongation.

Alternative polyadenylation mediates genetic regulation of gene expression

Little is known about co-transcriptional or post-transcriptional regulatory mechanisms linking noncoding variation to variation in organismal traits. To begin addressing this gap, we used 3’ Seq to study the impact of genetic variation on alternative polyadenylation (APA) in the nuclear and total mRNA fractions of 52 HapMap Yoruba human lymphoblastoid cell lines. We mapped 602 APA quantitative trait loci (apaQTLs) at 10% FDR, of which 152 were nuclear specific. Effect sizes at intronic apaQTLs are negatively correlated with eQTL effect sizes. These observations suggest genetic variants can decrease mRNA expression levels by increasing usage of intronic PAS. We also identified 24 apaQTLs associated with protein levels, but not mRNA expression. Finally, we found that 19% of apaQTLs can be associated with disease. Thus, our work demonstrates that APA links genetic variation to variation in gene expression, protein expression, and disease risk, and reveals uncharted modes of genetic regulation.

Emerging Roles for 3' UTRs in Neurons

The 3′ untranslated regions (3′ UTRs) of mRNAs serve as hubs for post-transcriptional control as the targets of microRNAs (miRNAs) and RNA-binding proteins (RBPs). Sequences in 3′ UTRs confer alterations in mRNA stability, direct mRNA localization to subcellular regions, and impart translational control. Thousands of mRNAs are localized to subcellular compartments in neurons—including axons, dendrites, and synapses—where they are thought to undergo local translation. Despite an established role for 3′ UTR sequences in imparting mRNA localization in neurons, the specific RNA sequences and structural features at play remain poorly understood. The nervous system selectively expresses longer 3′ UTR isoforms via alternative polyadenylation (APA). The regulation of APA in neurons and the neuronal functions of longer 3′ UTR mRNA isoforms are starting to be uncovered. Surprising roles for 3′ UTRs are emerging beyond the regulation of protein synthesis and include roles as RBP delivery scaffolds and regulators of alternative splicing. Evidence is also emerging that 3′ UTRs can be cleaved, leading to stable, isolated 3′ UTR fragments which are of unknown function. Mutations in 3′ UTRs are implicated in several neurological disorders—more studies are needed to uncover how these mutations impact gene regulation and what is their relationship to disease severity.

Pan-tissue analysis of allelic alternative polyadenylation suggests widespread functional regulation

Alternative polyadenylation (APA) is a major layer of gene regulation. However, it has recently been argued that most APA represents molecular noise. To clarify their functional relevance and evolution, we quantified allele-specific APA patterns in multiple tissues from an F1 hybrid mouse. We found a clearly negative correlation between gene expression and APA diversity for the 2,866 genes (24.9%) with a dominant polyadenylation site (PAS) usage above or equal to 90%, suggesting that their other PASs represent molecular errors. Among the remaining genes with multiple PASs, 3,971 genes (34.5%) express two or more isoforms with potentially functional importance. Interestingly, the genes with potentially functional minor PASs specific to neuronal tissues often express two APA isoforms with distinct subcellular localizations. Furthermore, our analysis of cis-APA divergence shows its pattern across tissues is distinct from that of gene expression. Finally, we demonstrate that the relative usage of alternative PASs is not only affected by their cis-regulatory elements, but also by potential coupling between transcriptional and APA regulation as well as competition kinetics between alternative sites.

RNA Polymerase II CTD phosphatase Rtr1 fine-tunes transcription termination

RNA Polymerase II (RNAPII) transcription termination is regulated by the phosphorylation status of the C-terminal domain (CTD). The phosphatase Rtr1 has been shown to regulate serine 5 phosphorylation on the CTD; however, its role in the regulation of RNAPII termination has not been explored. As a consequence of RTR1 deletion, interactions within the termination machinery and between the termination machinery and RNAPII were altered as quantified by Disruption-Compensation (DisCo) network analysis. Of note, interactions between RNAPII and the cleavage factor IA (CF1A) subunit Pcf11 were reduced in rtr1Δ, whereas interactions with the CTD and RNA-binding termination factor Nrd1 were increased. Globally, rtr1Δ leads to decreases in numerous noncoding RNAs that are linked to the Nrd1, Nab3 and Sen1 (NNS) -dependent RNAPII termination pathway. Genome-wide analysis of RNAPII and Nrd1 occupancy suggests that loss of RTR1 leads to increased termination at noncoding genes. Additionally, premature RNAPII termination increases globally at protein-coding genes with a decrease in RNAPII occupancy occurring just after the peak of Nrd1 recruitment during early elongation. The effects of rtr1Δ on RNA expression levels were lost following deletion of the exosome subunit Rrp6, which works with the NNS complex to rapidly degrade a number of noncoding RNAs following termination. Overall, these data suggest that Rtr1 restricts the NNS-dependent termination pathway in WT cells to prevent premature termination of mRNAs and ncRNAs. Rtr1 facilitates low-level elongation of noncoding transcripts that impact RNAPII interference thereby shaping the transcriptome.

Many cellular RNAs including those that encode for proteins are produced by the enzyme RNA Polymerase II. In this work, we have defined a new role for the phosphatase Rtr1 in the regulation of RNA Polymerase II progression from the start of transcription to the 3’ end of the gene where the nascent RNA from protein-coding genes is typically cleaved and polyadenylated. Deletion of the gene that encodes RTR1 leads to changes in the interactions between RNA polymerase II and the termination machinery. Rtr1 loss also causes early termination of RNA Polymerase II at many of its target gene types, including protein coding genes and noncoding RNAs. Evidence suggests that the premature termination observed in RTR1 knockout cells occurs through the termination factor and RNA binding protein Nrd1 and its binding partner Nab3. Deletion of RRP6, a known component of the Nrd1-Nab3 termination coupled RNA degradation pathway, is epistatic to RTR1 suggesting that Rrp6 is required to terminate and/or degrade many of the noncoding RNAs that have increased turnover in RTR1 deletion cells. These findings suggest that Rtr1 normally promotes elongation of RNA Polymerase II transcripts through prevention of Nrd1-directed termination.

Transcriptional regulation in model organisms: recent progress and clinical implications

In this review, we will summarize model organisms used by scientists in the laboratory, including Escherichia coli, yeast, Arabidopsis thaliana, nematodes, Drosophila, zebrafish, mice and other animals. We focus on the progress in research exploring different types of E. coli in the human body, and the specific molecular mechanisms by which they play a role in humans. First, we discuss the specific transcriptional regulation mechanism of E. coli in cell development, maturation, ageing and longevity, as well as tumorigenesis and development. Then, we discuss how the synthesis of some important substances in cells is regulated and how this affects biological behaviour. Understanding and applying these mechanisms, presumably, can greatly improve the quality of people's lives as well as increase their lifespan. For example, some E. coli can activate certain cells by secreting insulin-like growth factor-1, thus activating the inflammatory response of the body, while other E. coli can inactivate the immune response of the body by secreting toxic factors.

Targeting mRNA processing as an anticancer strategy

Discoveries in the past decade have highlighted the potential of mRNA as a therapeutic target for cancer. Specifically, RNA sequencing revealed that, in addition to gene mutations, alterations in mRNA can contribute to the initiation and progression of cancer. Indeed, precursor mRNA processing, which includes the removal of introns by splicing and the formation of 3' ends by cleavage and polyadenylation, is frequently altered in tumours. These alterations result in numerous cancer-specific mRNAs that generate altered levels of normal proteins or proteins with new functions, leading to the activation of oncogenes or the inactivation of tumour-suppressor genes. Abnormally spliced and polyadenylated mRNAs are also associated with resistance to cancer treatment and, unexpectedly, certain cancers are highly sensitive to the pharmacological inhibition of splicing. This Review summarizes recent progress in our understanding of how splicing and polyadenylation are altered in cancer and highlights how this knowledge has been translated for drug discovery, resulting in the production of small molecules and oligonucleotides that modulate the spliceosome and are in clinical trials for the treatment of cancer.

New means to an end: mRNA export activity impacts alternative polyadenylation

Gene expression involves multiple co- and post-transcriptional processes that have been increasingly found intertwined. A recent work by our groups (Chen et al. Mol Cell, 2019) indicates that expression of alternative polyadenylation isoforms in mammalian cells can be controlled by nuclear export activities. This regulation has distinct impacts on genes having different sizes and nucleotide contents, and involves RNA polymerase II distribution toward the 3' end of genes. This work raises a number of intriguing questions concerning how 3' end processing and nuclear export are integrated and how their regulation feeds back to transcription.

Transcriptional Control by Premature Termination: A Forgotten Mechanism

The concept of early termination as an important means of transcriptional control has long been established. Even so, its role in metazoan gene expression is underappreciated. Recent technological advances provide novel insights into premature transcription termination (PTT). This process is frequent, widespread, and can occur close to the transcription start site (TSS), or within the gene body. Stable prematurely terminated transcripts contribute to the transcriptome as instances of alternative polyadenylation (APA). Independently of transcript stability and function, premature termination opposes the formation of full-length transcripts, thereby negatively regulating gene expression, especially of transcriptional regulators. Premature termination can be beneficial or harmful, depending on its context. As a result, multiple factors have evolved to control this process.

Paf1C regulates RNA polymerase II progression by modulating elongation rate

Elongation factor Paf1C regulates several stages of the RNA polymerase II (Pol II) transcription cycle, although it is unclear how it modulates Pol II distribution and progression in mammalian cells. We found that conditional ablation of Paf1 resulted in the accumulation of unphosphorylated and Ser5 phosphorylated Pol II around promoter-proximal regions and within the first 20 to 30 kb of gene bodies, respectively. Paf1 ablation did not impact the recruitment of other key elongation factors, namely, Spt5, Spt6, and the FACT complex, suggesting that Paf1 function may be mechanistically distinguishable from each of these factors. Moreover, loss of Paf1 triggered an increase in TSS-proximal nucleosome occupancy, which could impose a considerable barrier to Pol II elongation past TSS-proximal regions. Remarkably, accumulation of Ser5P in the first 20 to 30 kb coincided with reductions in histone H2B ubiquitylation within this region. Furthermore, we show that nascent RNA species accumulate within this window, suggesting a mechanism whereby Paf1 loss leads to aberrant, prematurely terminated transcripts and diminution of full-length transcripts. Importantly, we found that loss of Paf1 results in Pol II elongation rate defects with significant rate compression. Our findings suggest that Paf1C is critical for modulating Pol II elongation rates by functioning beyond the pause-release step as an "accelerator" over specific early gene body regions.

Regulation of Intronic Polyadenylation by PCF11 Impacts mRNA Expression of Long Genes

Regulation of cleavage and polyadenylation (CPA) affects gene expression and polyadenylation site (PAS) choice. Here, we report that the CPA and termination factor PCF11 modulates gene expression on the basis of gene size. Although downregulation of PCF11 leads to inhibition of short gene expression, long genes are upregulated because of suppressed intronic polyadenylation (IPA) enriched in large introns. We show that this regulatory scheme, named PCF11-mediated expression regulation through IPA (PEIPA), takes place in cell differentiation, during which downregulation of PCF11 is coupled with upregulation of long genes with functions in cell morphology, adhesion, and migration. PEIPA targets distinct gene sets in different cell contexts with similar rules. Furthermore, PCF11 is autoregulated through a conserved IPA site, the removal of which leads to global activation of PASs close to gene promotors. Therefore, PCF11 uses distinct mechanisms to regulate genes of different sizes, and its autoregulation maintains homeostasis of PAS usage in the cell.

Tissue-specific mechanisms of alternative polyadenylation: Testis, brain, and beyond (2018 update)

Alternative polyadenylation (APA) is how genes choose different sites for 3′ end formation for mRNAs during transcription. APA often occurs in a tissue‐ or developmental stage‐specific manner that can significantly affect gene activity by changing the protein product generated, the stability of the transcript, its localization within the cell, or its translatability. Despite the important regulatory effects that APA has on tissue‐specific gene expression, only a few examples have been characterized mechanistically. In this 2018 update to our 2010 review, we examine mechanisms for the control of APA and update our understanding of the older mechanisms since 2010. We once postulated the existence of tissue‐specific factors in APA. However, while a few tissue‐specific polyadenylation factors are known, the emerging conclusion is that the majority of APA is accomplished by altering levels of core polyadenylation proteins. Examples of those core proteins include CSTF2, CPSF1, and subunits of mammalian cleavage factor I. But despite support for these mechanisms, no one has yet documented any of these proteins changing in either a tissue‐specific or developmental manner. Given the profound effect that APA can have on gene expression and human health, improved understanding of tissue‐specific APA could lead to numerous advances in gene activity control.

This article is categorized under:

RNA Processing > 3′ End Processing

RNA in Disease and Development > RNA in Development

The mRNA Export Receptor NXF1 Coordinates Transcriptional Dynamics, Alternative Polyadenylation, and mRNA Export

Alternative polyadenylation (APA) produces mRNA isoforms with different 3' UTR lengths. Previous studies indicated that 3' end processing and mRNA export are intertwined in gene regulation. Here, we show that mRNA export factors generally facilitate usage of distal cleavage and polyadenylation sites (PASs), leading to long 3' UTR isoform expression. By focusing on the export receptor NXF1, which exhibits the most potent effect on APA in this study, we reveal several gene features that impact NXF1-dependent APA, including 3' UTR size, gene size, and AT content. Surprisingly, NXF1 downregulation results in RNA polymerase II (Pol II) accumulation at the 3' end of genes, correlating with its role in APA regulation. Moreover, NXF1 cooperates with CFI-68 to facilitate nuclear export of long 3' UTR isoform with UGUA motifs. Together, our work reveals important roles of NXF1 in coordinating transcriptional dynamics, 3' end processing, and nuclear export of long 3' UTR transcripts, implicating NXF1 as a nexus of gene regulation.

CDK12 regulates DNA repair genes by suppressing intronic polyadenylation

Mutations that attenuate homologous recombination (HR)-mediated repair promote tumorigenesis and sensitize cells to chemotherapeutics that cause replication fork collapse, a phenotype known as 'BRCAness'¹. BRCAness tumours arise from loss-of-function mutations in 22 genes¹. Of these genes, all but one (CDK12) function directly in the HR repair pathway¹. CDK12 phosphorylates serine 2 of the RNA polymerase II C-terminal domain heptapeptide repeat^2-7, a modification that regulates transcription elongation, splicing, and cleavage and polyadenylation^8,9. Genome-wide expression studies suggest that depletion of CDK12 abrogates the expression of several HR genes relatively specifically, thereby blunting HR repair^3-7,10,11. This observation suggests that the mutational status of CDK12 may predict sensitivity to targeted treatments against BRCAness, such as PARP1 inhibitors, and that CDK12 inhibitors may induce sensitization of HR-competent tumours to these treatments^6,7,10,11. Despite growing clinical interest, the mechanism by which CDK12 regulates HR genes remains unknown. Here we show that CDK12 globally suppresses intronic polyadenylation events in mouse embryonic stem cells, enabling the production of full-length gene products. Many HR genes harbour more intronic polyadenylation sites than other expressed genes, and these sites are particularly sensitive to loss of CDK12. The cumulative effect of these sites accounts for the enhanced sensitivity of HR gene expression to CDK12 loss, and we find that this mechanism is conserved in human tumours that contain loss-of-function CDK12 mutations. This work clarifies the function of CDK12 and underscores its potential both as a chemotherapeutic target and as a tumour biomarker.

Decoding Transcriptome Dynamics of Genome-Encoded Polyadenylation and Autoregulation with Small-Molecule Modulators of Alternative Polyadenylation

Alternative polyadenylation (APA) plays a critical role in regulating gene expression. However, the balance between genome-encoded APA processing and autoregulation by APA modulating RNA binding protein (RBP) factors is not well understood. We discovered two potent small-molecule modulators of APA (T4 and T5) that promote distal-to-proximal (DtoP) APA usage in multiple transcripts. Monotonically responsive APA events, induced by short exposure to T4 or T5, were defined in the transcriptome, allowing clear isolation of the genomic sequence features and RBP motifs associated with DtoP regulation. We found that longer vulnerable introns, enriched with distinctive A-rich motifs, were preferentially affected by DtoP APA, thus defining a core set of genes with genomically encoded DtoP regulation. Through APA response pattern and compound-small interfering RNA epistasis analysis of APA-associated RBP factors, we further demonstrated that DtoP APA usage is partly modulated by altered autoregulation of polyadenylate binding nuclear protein-1 signaling.

Codon usage biases co-evolve with transcription termination machinery to suppress premature cleavage and polyadenylation

Codon usage biases are found in all genomes and influence protein expression levels. The codon usage effect on protein expression was thought to be mainly due to its impact on translation. Here, we show that transcription termination is an important driving force for codon usage bias in eukaryotes. Using Neurospora crassa as a model organism, we demonstrated that introduction of rare codons results in premature transcription termination (PTT) within open reading frames and abolishment of full-length mRNA. PTT is a wide-spread phenomenon in Neurospora, and there is a strong negative correlation between codon usage bias and PTT events. Rare codons lead to the formation of putative poly(A) signals and PTT. A similar role for codon usage bias was also observed in mouse cells. Together, these results suggest that codon usage biases co-evolve with the transcription termination machinery to suppress premature termination of transcription and thus allow for optimal gene expression.

Transcription elongation rate affects nascent histone pre-mRNA folding and 3' end processing

In this study, Saldi et al. investigated how transcription elongation rate influences cotranscriptional pre-mRNA maturation. Their findings show that regulation of transcription speed can modulate pre-mRNA processing by changing nascent RNA structure and suggest a mechanism by which alternative processing could be controlled.

Transcription elongation rate influences cotranscriptional pre-mRNA maturation, but how such kinetic coupling works is poorly understood. The formation of nonadenylated histone mRNA 3′ ends requires recognition of an RNA structure by stem–loop-binding protein (SLBP). We report that slow transcription by mutant RNA polymerase II (Pol II) caused accumulation of polyadenylated histone mRNAs that extend past the stem–loop processing site. UV irradiation, which decelerates Pol II elongation, also induced long poly(A)⁺ histone transcripts. Inhibition of 3′ processing by slow Pol II correlates with failure to recruit SLBP to histone genes. Chemical probing of nascent RNA structure showed that the stem–loop fails to fold in transcripts made by slow Pol II, thereby explaining the absence of SLBP and failure to process 3′ ends. These results show that regulation of transcription speed can modulate pre-mRNA processing by changing nascent RNA structure and suggest a mechanism by which alternative processing could be controlled.