Multiplexed single-cell characterization of alternative polyadenylation regulators.

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity regulated by the cleavage and polyadenylation (CPA) machinery. To better understand how these proteins govern polyA site choice, we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a framework to detect perturbation-dependent changes in polyadenylation and characterize modules of co-regulated polyA sites. We find groups of intronic polyA sites regulated by distinct components of the nuclear RNA life cycle, including elongation, splicing, termination, and surveillance. We train and validate a deep neural network (APARENT-Perturb) for tandem polyA site usage, delineating a cis-regulatory code that predicts perturbation response and reveals interactions between regulatory complexes. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation.

Sites of transcription initiation drive mRNA isoform selection.

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

TDP-43 condensation properties specify its RNA-binding and regulatory repertoire.

Mutations causing amyotrophic lateral sclerosis (ALS) often affect the condensation properties of RNA-binding proteins (RBPs). However, the role of RBP condensation in the specificity and function of protein-RNA complexes remains unclear. We created a series of TDP-43 C-terminal domain (CTD) variants that exhibited a gradient of low to high condensation propensity, as observed in vitro and by nuclear mobility and foci formation. Notably, a capacity for condensation was required for efficient TDP-43 assembly on subsets of RNA-binding regions, which contain unusually long clusters of motifs of characteristic types and density. These "binding-region condensates" are promoted by homomeric CTD-driven interactions and required for efficient regulation of a subset of bound transcripts, including autoregulation of TDP-43 mRNA. We establish that RBP condensation can occur in a binding-region-specific manner to selectively modulate transcriptome-wide RNA regulation, which has implications for remodeling RNA networks in the context of signaling, disease, and evolution.

A Deep Neural Network for Predicting and Engineering Alternative Polyadenylation.

Alternative polyadenylation (APA) is a major driver of transcriptome diversity in human cells. Here, we use deep learning to predict APA from DNA sequence alone. We trained our model (APARENT, APA REgression NeT) on isoform expression data from over 3 million APA reporters. APARENT's predictions are highly accurate when tasked with inferring APA in synthetic and human 3'UTRs. Visualizing features learned across all network layers reveals that APARENT recognizes sequence motifs known to recruit APA regulators, discovers previously unknown sequence determinants of 3' end processing, and integrates these features into a comprehensive, interpretable, cis-regulatory code. We apply APARENT to forward engineer functional polyadenylation signals with precisely defined cleavage position and isoform usage and validate predictions experimentally. Finally, we use APARENT to quantify the impact of genetic variants on APA. Our approach detects pathogenic variants in a wide range of disease contexts, expanding our understanding of the genetic origins of disease.

Heteromeric RNP Assembly at LINEs Controls Lineage-Specific RNA Processing.

Long mammalian introns make it challenging for the RNA processing machinery to identify exons accurately. We find that LINE-derived sequences (LINEs) contribute to this selection by recruiting dozens of RNA-binding proteins (RBPs) to introns. This includes MATR3, which promotes binding of PTBP1 to multivalent binding sites within LINEs. Both RBPs repress splicing and 3' end processing within and around LINEs. Notably, repressive RBPs preferentially bind to evolutionarily young LINEs, which are located far from exons. These RBPs insulate the LINEs and the surrounding intronic regions from RNA processing. Upon evolutionary divergence, changes in RNA motifs within LINEs lead to gradual loss of their insulation. Hence, older LINEs are located closer to exons, are a common source of tissue-specific exons, and increasingly bind to RBPs that enhance RNA processing. Thus, LINEs are hubs for the assembly of repressive RBPs and also contribute to the evolution of new, lineage-specific transcripts in mammals. VIDEO ABSTRACT.

Nudt21 Controls Cell Fate by Connecting Alternative Polyadenylation to Chromatin Signaling.

Cell fate transitions involve rapid gene expression changes and global chromatin remodeling, yet the underlying regulatory pathways remain incompletely understood. Here, we identified the RNA-processing factor Nudt21 as a novel regulator of cell fate change using transcription-factor-induced reprogramming as a screening assay. Suppression of Nudt21 enhanced the generation of induced pluripotent stem cells, facilitated transdifferentiation into trophoblast stem cells, and impaired differentiation of myeloid precursors and embryonic stem cells, suggesting a broader role for Nudt21 in cell fate change. We show that Nudt21 directs differential polyadenylation of over 1,500 transcripts in cells acquiring pluripotency, although only a fraction changed protein levels. Remarkably, these proteins were strongly enriched for chromatin regulators, and their suppression neutralized the effect of Nudt21 during reprogramming. Collectively, our data uncover Nudt21 as a novel post-transcriptional regulator of cell fate and establish a direct, previously unappreciated link between alternative polyadenylation and chromatin signaling.

2'-O-methylation at internal sites on mRNA promotes mRNA stability.

2'-O-methylation (Nm) is a prominent RNA modification well known in noncoding RNAs and more recently also found at many mRNA internal sites. However, their function and base-resolution stoichiometry remain underexplored. Here, we investigate the transcriptome-wide effect of internal site Nm on mRNA stability. Combining nanopore sequencing with our developed machine learning method, NanoNm, we identify thousands of Nm sites on mRNAs with a single-base resolution. We observe a positive effect of FBL-mediated Nm modification on mRNA stability and expression level. Elevated FBL expression in cancer cells is associated with increased expression levels for 2'-O-methylated mRNAs of cancer pathways, implying the role of FBL in post-transcriptional regulation. Lastly, we find that FBL-mediated 2'-O-methylation connects to widespread 3' UTR shortening, a mechanism that globally increases RNA stability. Collectively, we demonstrate that FBL-mediated Nm modifications at mRNA internal sites regulate gene expression by enhancing mRNA stability.

Alternative polyadenylation determines the functional landscape of inverted Alu repeats.

Inverted Alu repeats (IRAlus) are abundantly found in the transcriptome, especially in introns and 3' untranslated regions (UTRs). Yet, the biological significance of IRAlus embedded in 3' UTRs remains largely unknown. Here, we find that 3' UTR IRAlus silences genes involved in essential signaling pathways. We utilize J2 antibody to directly capture and map the double-stranded RNA structure of 3' UTR IRAlus in the transcriptome. Bioinformatic analysis reveals alternative polyadenylation as a major axis of IRAlus-mediated gene regulation. Notably, the expression of mouse double minute 2 (MDM2), an inhibitor of p53, is upregulated by the exclusion of IRAlus during UTR shortening, which is exploited to silence p53 during tumorigenesis. Moreover, the transcriptome-wide UTR lengthening in neural progenitor cells results in the global downregulation of genes associated with neurodegenerative diseases, including amyotrophic lateral sclerosis, via IRAlus inclusion. Our study establishes the functional landscape of 3' UTR IRAlus and its role in human pathophysiology.

Proximal termination generates a transcriptional state that determines the rate of establishment of Polycomb silencing.

The mechanisms and timescales controlling de novo establishment of chromatin-mediated transcriptional silencing by Polycomb repressive complex 2 (PRC2) are unclear. Here, we investigate PRC2 silencing at Arabidopsis FLOWERING LOCUS C (FLC), known to involve co-transcriptional RNA processing, histone demethylation activity, and PRC2 function, but so far not mechanistically connected. We develop and test a computational model describing proximal polyadenylation/termination mediated by the RNA-binding protein FCA that induces H3K4me1 removal by the histone demethylase FLD. H3K4me1 removal feeds back to reduce RNA polymerase II (RNA Pol II) processivity and thus enhance early termination, thereby repressing productive transcription. The model predicts that this transcription-coupled repression controls the level of transcriptional antagonism to PRC2 action. Thus, the effectiveness of this repression dictates the timescale for establishment of PRC2/H3K27me3 silencing. We experimentally validate these mechanistic model predictions, revealing that co-transcriptional processing sets the level of productive transcription at the locus, which then determines the rate of the ON-to-OFF switch to PRC2 silencing.

Regulation of alternative polyadenylation by the C2H2-zinc-finger protein Sp1.

Alternative polyadenylation (APA) enhances gene regulatory potential by increasing the diversity of mRNA transcripts. 3' UTR shortening through APA correlates with enhanced cellular proliferation and is a widespread phenomenon in tumor cells. Here, we show that the ubiquitously expressed transcription factor Sp1 binds RNA in vivo and is a common repressor of distal poly(A) site usage. RNA sequencing identified 2,344 genes (36% of the total mapped mRNA transcripts) with lengthened 3' UTRs upon Sp1 depletion. Sp1 preferentially binds the 3' UTRs of such lengthened transcripts and inhibits cleavage at distal sites by interacting with the subunits of the core cleavage and polyadenylation (CPA) machinery. The 3' UTR lengths of Sp1 target genes in breast cancer patient RNA-seq data correlate with Sp1 expression levels, implicating Sp1-mediated APA regulation in modulating tumorigenic properties. Taken together, our findings provide insights into the mechanism for dynamic APA regulation by unraveling a previously unknown function of the DNA-binding transcription factor Sp1.

Heat shock induces premature transcript termination and reconfigures the human transcriptome.

The heat shock (HS) response involves rapid induction of HS genes, whereas transcriptional repression is established more slowly at most other genes. Previous data suggested that such repression results from inhibition of RNA polymerase II (RNAPII) pause release, but here, we show that HS strongly affects other phases of the transcription cycle. Intriguingly, while elongation rates increase upon HS, processivity markedly decreases, so that RNAPII frequently fails to reach the end of genes. Indeed, HS results in widespread premature transcript termination at cryptic, intronic polyadenylation (IPA) sites near gene 5'-ends, likely via inhibition of U1 telescripting. This results in dramatic reconfiguration of the human transcriptome with production of new, previously unannotated, short mRNAs that accumulate in the nucleus. Together, these results shed new light on the basic transcription mechanisms induced by growth at elevated temperature and show that a genome-wide shift toward usage of IPA sites can occur under physiological conditions.

Developmentally regulated alternate 3' end cleavage of nascent transcripts controls dynamic changes in protein expression in an adult stem cell lineage.

Alternative polyadenylation (APA) generates transcript isoforms that differ in the position of the 3' cleavage site, resulting in the production of mRNA isoforms with different length 3' UTRs. Although widespread, the role of APA in the biology of cells, tissues, and organisms has been controversial. We identified >500 Drosophila genes that express mRNA isoforms with a long 3' UTR in proliferating spermatogonia but a short 3' UTR in differentiating spermatocytes due to APA. We show that the stage-specific choice of the 3' end cleavage site can be regulated by the arrangement of a canonical polyadenylation signal (PAS) near the distal cleavage site but a variant or no recognizable PAS near the proximal cleavage site. The emergence of transcripts with shorter 3' UTRs in differentiating cells correlated with changes in expression of the encoded proteins, either from off in spermatogonia to on in spermatocytes or vice versa. Polysome gradient fractionation revealed >250 genes where the long 3' UTR versus short 3' UTR mRNA isoforms migrated differently, consistent with dramatic stage-specific changes in translation state. Thus, the developmentally regulated choice of an alternative site at which to make the 3' end cut that terminates nascent transcripts can profoundly affect the suite of proteins expressed as cells advance through sequential steps in a differentiation lineage.

Overlapping Activities of ELAV/Hu Family RNA Binding Proteins Specify the Extended Neuronal 3' UTR Landscape in Drosophila.

The tissue-specific deployment of highly extended neural 3' UTR isoforms, generated by alternative polyadenylation (APA), is a broad and conserved feature of metazoan genomes. However, the factors and mechanisms that control neural APA isoforms are not well understood. Here, we show that three ELAV/Hu RNA binding proteins (Elav, Rbp9, and Fne) have similar capacities to induce a lengthened 3' UTR landscape in an ectopic setting. These factors promote accumulation of chromatin-associated, 3' UTR-extended, nascent transcripts, through inhibition of proximal polyadenylation site (PAS) usage. Notably, Elav represses an unannotated splice isoform of fne, switching the normally cytoplasmic Fne toward the nucleus in elav mutants. We use genomic profiling to reveal strong and broad loss of neural APA in elav/fne double mutant CNS, the first genetic background to largely abrogate this distinct APA signature. Overall, we demonstrate how regulatory interplay and functionally overlapping activities of neural ELAV/Hu RBPs drives the neural APA landscape.

ELAV and FNE Determine Neuronal Transcript Signatures through EXon-Activated Rescue.

The production of alternative RNA variants contributes to the tissue-specific regulation of gene expression. In the animal nervous system, a systematic shift toward distal sites of transcription termination produces transcript signatures that are crucial for neuron development and function. Here, we report that, in Drosophila, the highly conserved protein ELAV globally regulates all sites of neuronal 3' end processing and directly binds to proximal polyadenylation sites of target mRNAs in vivo. We uncover an endogenous strategy of functional gene rescue that safeguards neuronal RNA signatures in an ELAV loss-of-function context. When not directly repressed by ELAV, the transcript encoding the ELAV paralog FNE acquires a mini-exon, generating a new protein able to translocate to the nucleus and rescue ELAV-mediated alternative polyadenylation and alternative splicing. We propose that exon-activated functional rescue is a more widespread mechanism that ensures robustness of processes regulated by a hierarchy, rather than redundancy, of effectors.

A Cancer-Specific Ubiquitin Ligase Drives mRNA Alternative Polyadenylation by Ubiquitinating the mRNA 3' End Processing Complex.

Alternative polyadenylation (APA) contributes to transcriptome complexity by generating mRNA isoforms with varying 3' UTR lengths. APA leading to 3' UTR shortening (3' US) is a common feature of most cancer cells; however, the molecular mechanisms are not understood. Here, we describe a widespread mechanism promoting 3' US in cancer through ubiquitination of the mRNA 3' end processing complex protein, PCF11, by the cancer-specific MAGE-A11-HUWE1 ubiquitin ligase. MAGE-A11 is normally expressed only in the male germline but is frequently re-activated in cancers. MAGE-A11 is necessary for cancer cell viability and is sufficient to drive tumorigenesis. Screening for targets of MAGE-A11 revealed that it ubiquitinates PCF11, resulting in loss of CFIm25 from the mRNA 3' end processing complex. This leads to APA of many transcripts affecting core oncogenic and tumor suppressors, including cyclin D2 and PTEN. These findings provide insights into the molecular mechanisms driving APA in cancer and suggest therapeutic strategies.

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

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 quantitative trait methylation mapping in human cancers provides clues into the molecular mechanisms of APA.

Genetic variants are involved in the orchestration of alternative polyadenylation (APA) events, while the role of DNA methylation in regulating APA remains unclear. We generated a comprehensive atlas of APA quantitative trait methylation sites (apaQTMs) across 21 different types of cancer (1,612 to 60,219 acting in cis and 4,448 to 142,349 in trans). Potential causal apaQTMs in non-cancer samples were also identified. Mechanistically, we observed a strong enrichment of cis-apaQTMs near polyadenylation sites (PASs) and both cis- and trans-apaQTMs in proximity to transcription factor (TF) binding regions. Through the integration of ChIP-signals and RNA-seq data from cell lines, we have identified several regulators of APA events, acting either directly or indirectly, implicating novel functions of some important genes, such as TCF7L2, which is known for its involvement in type 2 diabetes and cancers. Furthermore, we have identified a vast number of QTMs that share the same putative causal CpG sites with five different cancer types, underscoring the roles of QTMs, including apaQTMs, in the process of tumorigenesis. DNA methylation is extensively involved in the regulation of APA events in human cancers. In an attempt to elucidate the potential underlying molecular mechanisms of APA by DNA methylation, our study paves the way for subsequent experimental validations into the intricate biological functions of DNA methylation in APA regulation and the pathogenesis of human cancers. To present a comprehensive catalog of apaQTM patterns, we introduce the Pancan-apaQTM database, available at https://pancan-apaqtm-zju.shinyapps.io/pancanaQTM/.

Co-transcriptional Loading of RNA Export Factors Shapes the Human Transcriptome.

During gene expression, RNA export factors are mainly known for driving nucleo-cytoplasmic transport. While early studies suggested that the exon junction complex (EJC) provides a binding platform for them, subsequent work proposed that they are only recruited by the cap binding complex to the 5' end of RNAs, as part of TREX. Using iCLIP, we show that the export receptor Nxf1 and two TREX subunits, Alyref and Chtop, are recruited to the whole mRNA co-transcriptionally via splicing but before 3' end processing. Consequently, Alyref alters splicing decisions and Chtop regulates alternative polyadenylation. Alyref is recruited to the 5' end of RNAs by CBC, and our data reveal subsequent binding to RNAs near EJCs. We demonstrate that eIF4A3 stimulates Alyref deposition not only on spliced RNAs close to EJC sites but also on single-exon transcripts. Our study reveals mechanistic insights into the co-transcriptional recruitment of mRNA export factors and how this shapes the human transcriptome.

Selective Roles of Vertebrate PCF11 in Premature and Full-Length Transcript Termination.

The pervasive nature of RNA polymerase II (Pol II) transcription requires efficient termination. A key player in this process is the cleavage and polyadenylation (CPA) factor PCF11, which directly binds to the Pol II C-terminal domain and dismantles elongating Pol II from DNA in vitro. We demonstrate that PCF11-mediated termination is essential for vertebrate development. A range of genomic analyses, including mNET-seq, 3' mRNA-seq, chromatin RNA-seq, and ChIP-seq, reveals that PCF11 enhances transcription termination and stimulates early polyadenylation genome-wide. PCF11 binds preferentially between closely spaced genes, where it prevents transcriptional interference and consequent gene downregulation. Notably, PCF11 is sub-stoichiometric to the CPA complex. Low levels of PCF11 are maintained by an auto-regulatory mechanism involving premature termination of its own transcript and are important for normal development. Both in human cell culture and during zebrafish development, PCF11 selectively attenuates the expression of other transcriptional regulators by premature CPA and termination.

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.