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.

Dynamic imaging of nascent RNA reveals general principles of transcription dynamics and stochastic splice site selection.

The activities of RNA polymerase and the spliceosome are responsible for the heterogeneity in the abundance and isoform composition of mRNA in human cells. However, the dynamics of these megadalton enzymatic complexes working in concert on endogenous genes have not been described. Here, we establish a quasi-genome-scale platform for observing synthesis and processing kinetics of single nascent RNA molecules in real time. We find that all observed genes show transcriptional bursting. We also observe large kinetic variation in intron removal for single introns in single cells, which is inconsistent with deterministic splice site selection. Transcriptome-wide footprinting of the U2AF complex, nascent RNA profiling, long-read sequencing, and lariat sequencing further reveal widespread stochastic recursive splicing within introns. We propose and validate a unified theoretical model to explain the general features of transcription and pervasive stochastic splice site selection.

Splice site m6A methylation prevents binding of U2AF35 to inhibit RNA splicing.

The N6-methyladenosine (m6A) RNA modification is used widely to alter the fate of mRNAs. Here we demonstrate that the C. elegans writer METT-10 (the ortholog of mouse METTL16) deposits an m6A mark on the 3' splice site (AG) of the S-adenosylmethionine (SAM) synthetase pre-mRNA, which inhibits its proper splicing and protein production. The mechanism is triggered by a rich diet and acts as an m6A-mediated switch to stop SAM production and regulate its homeostasis. Although the mammalian SAM synthetase pre-mRNA is not regulated via this mechanism, we show that splicing inhibition by 3' splice site m6A is conserved in mammals. The modification functions by physically preventing the essential splicing factor U2AF35 from recognizing the 3' splice site. We propose that use of splice-site m6A is an ancient mechanism for splicing regulation.

Spliceosome-targeted therapies trigger an antiviral immune response in triple-negative breast cancer.

Many oncogenic insults deregulate RNA splicing, often leading to hypersensitivity of tumors to spliceosome-targeted therapies (STTs). However, the mechanisms by which STTs selectively kill cancers remain largely unknown. Herein, we discover that mis-spliced RNA itself is a molecular trigger for tumor killing through viral mimicry. In MYC-driven triple-negative breast cancer, STTs cause widespread cytoplasmic accumulation of mis-spliced mRNAs, many of which form double-stranded structures. Double-stranded RNA (dsRNA)-binding proteins recognize these endogenous dsRNAs, triggering antiviral signaling and extrinsic apoptosis. In immune-competent models of breast cancer, STTs cause tumor cell-intrinsic antiviral signaling, downstream adaptive immune signaling, and tumor cell death. Furthermore, RNA mis-splicing in human breast cancers correlates with innate and adaptive immune signatures, especially in MYC-amplified tumors that are typically immune cold. These findings indicate that dsRNA-sensing pathways respond to global aberrations of RNA splicing in cancer and provoke the hypothesis that STTs may provide unexplored strategies to activate anti-tumor immune pathways.

Predicting Splicing from Primary Sequence with Deep Learning.

The splicing of pre-mRNAs into mature transcripts is remarkable for its precision, but the mechanisms by which the cellular machinery achieves such specificity are incompletely understood. Here, we describe a deep neural network that accurately predicts splice junctions from an arbitrary pre-mRNA transcript sequence, enabling precise prediction of noncoding genetic variants that cause cryptic splicing. Synonymous and intronic mutations with predicted splice-altering consequence validate at a high rate on RNA-seq and are strongly deleterious in the human population. De novo mutations with predicted splice-altering consequence are significantly enriched in patients with autism and intellectual disability compared to healthy controls and validate against RNA-seq in 21 out of 28 of these patients. We estimate that 9%-11% of pathogenic mutations in patients with rare genetic disorders are caused by this previously underappreciated class of disease variation.

Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing.

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

Dissecting the Causal Mechanism of X-Linked Dystonia-Parkinsonism by Integrating Genome and Transcriptome Assembly.

X-linked Dystonia-Parkinsonism (XDP) is a Mendelian neurodegenerative disease that is endemic to the Philippines and is associated with a founder haplotype. We integrated multiple genome and transcriptome assembly technologies to narrow the causal mutation to the TAF1 locus, which included a SINE-VNTR-Alu (SVA) retrotransposition into intron 32 of the gene. Transcriptome analyses identified decreased expression of the canonical cTAF1 transcript among XDP probands, and de novo assembly across multiple pluripotent stem-cell-derived neuronal lineages discovered aberrant TAF1 transcription that involved alternative splicing and intron retention (IR) in proximity to the SVA that was anti-correlated with overall TAF1 expression. CRISPR/Cas9 excision of the SVA rescued this XDP-specific transcriptional signature and normalized TAF1 expression in probands. These data suggest an SVA-mediated aberrant transcriptional mechanism associated with XDP and may provide a roadmap for layered technologies and integrated assembly-based analyses for other unsolved Mendelian disorders.

Inborn Errors of RNA Lariat Metabolism in Humans with Brainstem Viral Infection.

Viruses that are typically benign sometimes invade the brainstem in otherwise healthy children. We report bi-allelic DBR1 mutations in unrelated patients from different ethnicities, each of whom had brainstem infection due to herpes simplex virus 1 (HSV1), influenza virus, or norovirus. DBR1 encodes the only known RNA lariat debranching enzyme. We show that DBR1 expression is ubiquitous, but strongest in the spinal cord and brainstem. We also show that all DBR1 mutant alleles are severely hypomorphic, in terms of expression and function. The fibroblasts of DBR1-mutated patients contain higher RNA lariat levels than control cells, this difference becoming even more marked during HSV1 infection. Finally, we show that the patients' fibroblasts are highly susceptible to HSV1. RNA lariat accumulation and viral susceptibility are rescued by wild-type DBR1. Autosomal recessive, partial DBR1 deficiency underlies viral infection of the brainstem in humans through the disruption of tissue-specific and cell-intrinsic immunity to viruses.

The splicing regulators RBM5 and RBM10 are subunits of the U2 snRNP engaged with intron branch sites on chromatin.

Understanding the mechanisms of pre-mRNA splicing is limited by the technical challenges to examining spliceosomes in vivo. Here, we report the isolation of RNP complexes derived from precatalytic A or B-like spliceosomes solubilized from the chromatin pellet of mammalian cell nuclei. We found that these complexes contain U2 snRNP proteins and a portion of the U2 snRNA bound with protected RNA fragments that precisely map to intronic branch sites across the transcriptome. These U2 complexes also contained the splicing regulators RBM5 and RBM10. We found RBM5 and RBM10 bound to nearly all branch site complexes and not simply those at regulated exons. The deletion of a conserved RBM5/RBM10 peptide sequence, including a zinc finger motif, disrupted U2 interaction and rendered the proteins inactive for the repression of many alternative exons. We propose a model where RBM5 and RBM10 regulate splicing as components of the U2 snRNP complex following branch site base pairing.

Intron lariat spliceosomes convert lariats to true circles: implications for intron transposition.

Rare, full-length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envisioned and tested a hypothesis for their formation using Saccharomyces cerevisiae, documenting full-length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron lariat spliceosome (ILS) in which the 3'-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3' processing machinery) attacks the branch, joining the intron 3' end to the 5' splice site in a 3'-5' linked circle. Human U2 and U12 spliceosomes produce analogous full-length and processed circles. Postsplicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

Histone demethylase LSD1 is required for germinal center formation and BCL6-driven lymphomagenesis.

Germinal center (GC) B cells feature repression of many gene enhancers to establish their characteristic transcriptome. Here we show that conditional deletion of Lsd1 in GCs significantly impaired GC formation, associated with failure to repress immune synapse genes linked to GC exit, which are also direct targets of the transcriptional repressor BCL6. We found that BCL6 directly binds LSD1 and recruits it primarily to intergenic and intronic enhancers. Conditional deletion of Lsd1 suppressed GC hyperplasia caused by constitutive expression of BCL6 and significantly delayed BCL6-driven lymphomagenesis. Administration of catalytic inhibitors of LSD1 had little effect on GC formation or GC-derived lymphoma cells. Using a CRISPR-Cas9 domain screen, we found instead that the LSD1 Tower domain was critical for dependence on LSD1 in GC-derived B cells. These results indicate an essential role for LSD1 in the humoral immune response, where it modulates enhancer function by forming repression complexes with BCL6.

A functional link between lariat debranching enzyme and the intron-binding complex is defective in non-photosensitive trichothiodystrophy.

The pre-mRNA life cycle requires intron processing; yet, how intron-processing defects influence splicing and gene expression is unclear. Here, we find that TTDN1/MPLKIP, which is encoded by a gene implicated in non-photosensitive trichothiodystrophy (NP-TTD), functionally links intron lariat processing to spliceosomal function. The conserved TTDN1 C-terminal region directly binds lariat debranching enzyme DBR1, whereas its N-terminal intrinsically disordered region (IDR) binds the intron-binding complex (IBC). TTDN1 loss, or a mutated IDR, causes significant intron lariat accumulation, as well as splicing and gene expression defects, mirroring phenotypes observed in NP-TTD patient cells. A Ttdn1-deficient mouse model recapitulates intron-processing defects and certain neurodevelopmental phenotypes seen in NP-TTD. Fusing DBR1 to the TTDN1 IDR is sufficient to recruit DBR1 to the IBC and circumvents the functional requirement for TTDN1. Collectively, our findings link RNA lariat processing with splicing outcomes by revealing the molecular function of TTDN1.

Recognition of transcription terminators during retrotransposition: How to keep a group II intron quiet.

Chung et al. recently presented the structure of a primitive group IIC intron with its DNA target, which reveals the structural requirements that this class of intron uses to recognize a transcription terminator stem loop at the DNA level for insertion during retrotransposition.

PTBP1-activated co-transcriptional splicing controls epigenetic status of pluripotent stem cells.

Many spliceosomal introns are excised from nascent transcripts emerging from RNA polymerase II (RNA Pol II). The extent of cell-type-specific regulation and possible functions of such co-transcriptional events remain poorly understood. We examined the role of the RNA-binding protein PTBP1 in this process using an acute depletion approach followed by the analysis of chromatin- and RNA Pol II-associated transcripts. We show that PTBP1 activates the co-transcriptional excision of hundreds of introns, a surprising effect given that this protein is known to promote intron retention. Importantly, some co-transcriptionally activated introns fail to complete their splicing without PTBP1. In a striking example, retention of a PTBP1-dependent intron triggers nonsense-mediated decay of transcripts encoding DNA methyltransferase DNMT3B. We provide evidence that this regulation facilitates the natural decline in DNMT3B levels in developing neurons and protects differentiation-specific genes from ectopic methylation. Thus, PTBP1-activated co-transcriptional splicing is a widespread phenomenon mediating epigenetic control of cellular identity.

Structural basis of pre-tRNA intron removal by human tRNA splicing endonuclease.

Removal of the intron from precursor-tRNA (pre-tRNA) is essential in all three kingdoms of life. In humans, this process is mediated by the tRNA splicing endonuclease (TSEN) comprising four subunits: TSEN2, TSEN15, TSEN34, and TSEN54. Here, we report the cryo-EM structures of human TSEN bound to full-length pre-tRNA in the pre-catalytic and post-catalytic states at average resolutions of 2.94 and 2.88 Å, respectively. Human TSEN features an extended surface groove that holds the L-shaped pre-tRNA. The mature domain of pre-tRNA is recognized by conserved structural elements of TSEN34, TSEN54, and TSEN2. Such recognition orients the anticodon stem of pre-tRNA and places the 3'-splice site and 5'-splice site into the catalytic centers of TSEN34 and TSEN2, respectively. The bulk of the intron sequences makes no direct interaction with TSEN, explaining why pre-tRNAs of varying introns can be accommodated and cleaved. Our structures reveal the molecular ruler mechanism of pre-tRNA cleavage by TSEN.

Minor intron splicing is critical for survival of lethal prostate cancer.

The evolutionarily conserved minor spliceosome (MiS) is required for protein expression of ∼714 minor intron-containing genes (MIGs) crucial for cell-cycle regulation, DNA repair, and MAP-kinase signaling. We explored the role of MIGs and MiS in cancer, taking prostate cancer (PCa) as an exemplar. Both androgen receptor signaling and elevated levels of U6atac, a MiS small nuclear RNA, regulate MiS activity, which is highest in advanced metastatic PCa. siU6atac-mediated MiS inhibition in PCa in vitro model systems resulted in aberrant minor intron splicing leading to cell-cycle G1 arrest. Small interfering RNA knocking down U6atac was ∼50% more efficient in lowering tumor burden in models of advanced therapy-resistant PCa compared with standard antiandrogen therapy. In lethal PCa, siU6atac disrupted the splicing of a crucial lineage dependency factor, the RE1-silencing factor (REST). Taken together, we have nominated MiS as a vulnerability for lethal PCa and potentially other cancers.

FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns.

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

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.

Role of the proline-rich disordered domain of DROSHA in intronic microRNA processing.

DROSHA serves as a gatekeeper of the microRNA (miRNA) pathway by processing primary transcripts (pri-miRNAs). While the functions of structured domains of DROSHA have been well documented, the contribution of N-terminal proline-rich disordered domain (PRD) remains elusive. Here we show that the PRD promotes the processing of miRNA hairpins located within introns. We identified a DROSHA isoform (p140) lacking the PRD, which is produced by proteolytic cleavage. Small RNA sequencing revealed that p140 is significantly impaired in the maturation of intronic miRNAs. Consistently, our minigene constructs demonstrated that PRD enhances the processing of intronic hairpins, but not those in exons. Splice site mutations did not affect the PRD's enhancing effect on intronic constructs, suggesting that the PRD acts independently of splicing reaction by interacting with sequences residing within introns. The N-terminal regions from zebrafish and Xenopus DROSHA can replace the human counterpart, indicating functional conservation despite poor sequence alignment. Moreover, we found that rapidly evolving intronic miRNAs are generally more dependent on PRD than conserved ones, suggesting a role of PRD in miRNA evolution. Our study reveals a new layer of miRNA regulation mediated by a low-complexity disordered domain that senses the genomic contexts of miRNA loci.