Tracking pre-mRNA maturation across subcellular compartments identifies developmental gene regulation through intron retention and nuclear anchoring

Shiba: A unified computational method for robust identification of differential RNA splicing across platforms

Alternative pre-mRNA splicing (AS) is a fundamental regulatory process that generates transcript diversity and cell type variation. We developed Shiba, a robust method integrating transcript assembly, splicing event identification, read counting, and statistical analysis, to efficiently quantify exon splicing levels across various types of RNA-seq datasets. Compared to existing pipelines, Shiba excels in capturing both annotated and unannotated or cryptic differential splicing events with superior accuracy, sensitivity, and reproducibility. Furthermore, Shiba's unique consideration of junction read imbalance and exon-body read coverage reduces false positives, essential for downstream functional analyses. We have further developed scShiba for single-cell/nucleus (sc/sn) RNA-seq data, enabling the exploration of splicing variations in heterogeneous cell populations. Both simulated and real data demonstrate Shiba's robustness across multiple sample sizes, including n=1 datasets and individual cell clusters from scRNA-seq. Application of Shiba on single replicates of RNA-seq identified new AS-NMD targets, and scShiba on snRNA-seq revealed intricate temporal AS regulation in dopaminergic neurons. Both Shiba and scShiba are provided in Docker/Singularity containers and Snakemake pipeline, enhancing accessibility and reproducibility. The comprehensive capabilities of Shiba and scShiba allow systematic and robust quantification of alternative splicing events, laying a solid foundation for mechanistic exploration of functional complexity in RNA splicing.

Mammalian RNAi represses pericentromeric lncRNAs to maintain genome stability

Mammalian pericentromeric tandem repeats produce long noncoding RNAs (lncRNAs) that are dysregulated in cancer and linked to genomic instability. Identifying the basic molecular characteristics of these lncRNAs and their regulation is important to understanding their biological function. Here, we determine that the Argonaute (Ago) proteins of the RNA interference (RNAi) pathway directly and uniformly repress bidirectional pericentromeric lncRNAs in a Dicer-dependent manner in mouse embryonic and adult stem cells. Ago-dependent and Dicer-dependent autoregulatory small RNAs were identified within pericentromeric lncRNA degradation intermediates. We develop an RNase H cleavage assay to determine the relative proportions and lengths of the pericentromeric lncRNA targets. We find that 5'-phosphate and non-polyadenylated bidirectional pericentromeric lncRNAs are expressed at similar proportions. These lncRNAs can span up to 9 repeats, with transcription from the reverse strand template yielding the longer products. Using pericentromeric repeat RNA reporters, we determine that Ago represses pericentromeric lncRNAs after S phase transcription. Upon loss of Ago, pericentromeric lncRNA dysregulation results in delayed cell cycle progression, a defective mitotic spindle assembly checkpoint (SAC) and genomic instability. These results show that an evolutionarily conserved Ago activity at pericentromeres contributes to mammalian genome stability.

The lncRNA Malat1 is trafficked to the cytoplasm as a localized mRNA encoding a small peptide in neurons

Synaptic function in neurons is modulated by local translation of mRNAs that are transported to distal portions of axons and dendrites. The metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) is broadly expressed across cell types, almost exclusively as a nuclear long noncoding RNA. We found that in differentiating neurons, a portion of Malat1 RNA redistributes to the cytoplasm. Depletion of Malat1 using antisense oligonucleotides (ASOs) stimulates the expression of particular pre- and postsynaptic proteins, implicating Malat1 in their regulation. Neuronal Malat1 is localized in puncta of both axons and dendrites that costain with Staufen1 protein, similar to neuronal RNA granules formed by locally translated mRNAs. Ribosome profiling of cultured mouse cortical neurons identified ribosome footprints within a 5' region of Malat1 containing short open reading frames. The upstream-most reading frame (M1) of the Malat1 locus was linked to the GFP-coding sequence in mouse embryonic stem cells. When these gene-edited cells were differentiated into glutamatergic neurons, the M1-GFP fusion protein was expressed. Antibody staining for the M1 peptide confirmed its presence in wild-type neurons and showed that M1 expression was enhanced by synaptic stimulation with KCl. Our results indicate that Malat1 serves as a cytoplasmic coding RNA in the brain that is both modulated by and modulates synaptic function.

Timing is everything: advances in quantifying splicing kinetics

Splicing is a highly regulated process critical for proper pre-mRNA maturation and the maintenance of a healthy cellular environment. Splicing events are impacted by ongoing transcription, neighboring splicing events, and cis and trans regulatory factors on the respective pre-mRNA transcript. Within this complex regulatory environment, splicing kinetics have the potential to influence splicing outcomes but have historically been challenging to study in vivo. In this review, we highlight recent technological advancements that have enabled measurements of global splicing kinetics and of the variability of splicing kinetics at single introns. We demonstrate how identifying features that are correlated with splicing kinetics has increased our ability to form potential models for how splicing kinetics may be regulated in vivo.

Alternative splicing of a chromatin modifier alters the transcriptional regulatory programs of stem cell maintenance and neuronal differentiation

Development of embryonic stem cells (ESCs) into neurons requires intricate regulation of transcription, splicing, and translation, but how these processes interconnect is not understood. We found that polypyrimidine tract binding protein 1 (PTBP1) controls splicing of DPF2, a subunit of BRG1/BRM-associated factor (BAF) chromatin remodeling complexes. Dpf2 exon 7 splicing is inhibited by PTBP1 to produce the DPF2-S isoform early in development. During neuronal differentiation, loss of PTBP1 allows exon 7 inclusion and DPF2-L expression. Different cellular phenotypes and gene expression programs were induced by these alternative DPF2 isoforms. We identified chromatin binding sites enriched for each DPF2 isoform, as well as sites bound by both. In ESC, DPF2-S preferential sites were bound by pluripotency factors. In neuronal progenitors, DPF2-S sites were bound by nuclear factor I (NFI), while DPF2-L sites were bound by CCCTC-binding factor (CTCF). DPF2-S sites exhibited enhancer modifications, while DPF2-L sites showed promoter modifications. Thus, alternative splicing redirects BAF complex targeting to impact chromatin organization during neuronal development.

Systematic identification and characterization of exon-intron circRNAs

Exon-intron circRNAs (EIciRNAs) are a circRNA subclass with retained introns. Global features of EIciRNAs remain largely unexplored, mainly owing to the lack of bioinformatic tools. The regulation of intron retention (IR) in EIciRNAs and the associated functionality also require further investigation. We developed a framework, FEICP, which efficiently detected EIciRNAs from high-throughput sequencing (HTS) data. EIciRNAs are distinct from exonic circRNAs (EcircRNAs) in aspects such as with larger length, localization in the nucleus, high tissue specificity, and enrichment mostly in the brain. Deep learning analyses revealed that compared with regular introns, the retained introns of circRNAs (CIRs) are shorter in length, have weaker splice site strength, and have higher GC content. Compared with retained introns in linear RNAs (LIRs), CIRs are more likely to form secondary structures and show greater sequence conservation. CIRs are closer to the 5'-end, whereas LIRs are closer to the 3'-end of transcripts. EIciRNA-generating genes are more actively transcribed and associated with epigenetic marks of gene activation. Computational analyses and genome-wide CRISPR screening revealed that SRSF1 binds to CIRs and inhibits the biogenesis of most EIciRNAs. SRSF1 regulates the biogenesis of EIciLIMK1, which enhances the expression of LIMK1 in cis to boost neuronal differentiation, exemplifying EIciRNA physiological function. Overall, our study has developed the FEICP pipeline to identify EIciRNAs from HTS data, and reveals multiple features of CIRs and EIciRNAs. SRSF1 has been identified to regulate EIciRNA biogenesis. EIciRNAs and the regulation of EIciRNA biogenesis play critical roles in neuronal differentiation.

Emerging roles of RNA-binding proteins in fatty liver disease

A rampant and urgent global health issue of the 21st century is the emergence and progression of fatty liver disease (FLD), including alcoholic fatty liver disease and the more heterogenous metabolism-associated (or non-alcoholic) fatty liver disease (MAFLD/NAFLD) phenotypes. These conditions manifest as disease spectra, progressing from benign hepatic steatosis to symptomatic steatohepatitis, cirrhosis, and, ultimately, hepatocellular carcinoma. With numerous intricately regulated molecular pathways implicated in its pathophysiology, recent data have emphasized the critical roles of RNA-binding proteins (RBPs) in the onset and development of FLD. They regulate gene transcription and post-transcriptional processes, including pre-mRNA splicing, capping, and polyadenylation, as well as mature mRNA transport, stability, and translation. RBP dysfunction at every point along the mRNA life cycle has been associated with altered lipid metabolism and cellular stress response, resulting in hepatic inflammation and fibrosis. Here, we discuss the current understanding of the role of RBPs in the post-transcriptional processes associated with FLD and highlight the possible and emerging therapeutic strategies leveraging RBP function for FLD treatment. This article is categorized under: RNA in Disease and Development > RNA in Disease.

The lncRNA Malat1 is trafficked to the cytoplasm as a localized mRNA encoding a small peptide in neurons

Synaptic function is modulated by local translation of mRNAs that are transported to distal portions of axons and dendrites. The Metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) is broadly expressed across cell types, almost exclusively as a nuclear non-coding RNA. We found that in differentiating neurons, a portion of Malat1 RNA redistributes to the cytoplasm. Depletion of Malat1 from neurons stimulated expression of particular pre- and post- synaptic proteins, implicating Malat1 in their regulation. Neuronal Malat1 is localized to both axons and dendrites in puncta that co-stain with Staufen1 protein, similar to neuronal granules formed by locally translated mRNAs. Ribosome profiling of mouse cortical neurons identified ribosome footprints within a region of Malat1 containing short open reading frames. The upstream-most reading frame (M1) of the Malat1 locus was linked to the GFP coding sequence in mouse ES cells. When these gene-edited cells were differentiated into glutamatergic neurons, the M1-GFP fusion protein was expressed. Antibody staining for the M1 peptide confirmed its presence in wildtype neurons, and showed enhancement of M1 expression after synaptic stimulation with KCL. Our results indicate that Malat1 serves as a cytoplasmic coding RNA in the brain that is both modulated by and modulates synaptic function.

Subcellular RNA distribution and its change during human embryonic stem cell differentiation

The spatial localization of RNA within cells is closely related to its function and also involved in cell fate determination. However, the atlas of RNA distribution within cells and dynamic changes during the developmental process are largely unknown. In this study, five subcellular components, including cytoplasmic extract, membrane extract, soluble nuclear extract, chromatin-bound nuclear extract, and cytoskeletal extract, were isolated and the rules of subcellular RNA distribution in human embryonic stem cells (hESCs) and its change during hESC differentiation are summarized for the first time. The overall distribution patterns of coding and non-coding RNAs are revealed. Interestingly, some developmental genes are found to be transcribed but confined to the chromatin in undifferentiated hESC. Unexpectedly, alternative splicing and polyadenylation endow spatial heterogeneity among different isoforms of the same gene. Finally, the dynamic pattern of RNA distribution during hESC differentiation is characterized, which provides new clues for a comprehensive understanding hESC pluripotency and differentiation.

XIST directly regulates X-linked and autosomal genes in naive human pluripotent cells

X chromosome inactivation (XCI) serves as a paradigm for RNA-mediated regulation of gene expression, wherein the long non-coding RNA XIST spreads across the X chromosome in cis to mediate gene silencing chromosome-wide. In female naive human pluripotent stem cells (hPSCs), XIST is in a dispersed configuration, and XCI does not occur, raising questions about XIST's function. We found that XIST spreads across the X chromosome and induces dampening of X-linked gene expression in naive hPSCs. Surprisingly, XIST also targets specific autosomal regions, where it induces repressive chromatin changes and gene expression dampening. Thereby, XIST equalizes X-linked gene dosage between male and female cells while inducing differences in autosomes. The dispersed Xist configuration and autosomal localization also occur transiently during XCI initiation in mouse PSCs. Together, our study identifies XIST as the regulator of X chromosome dampening, uncovers an evolutionarily conserved trans-acting role of XIST/Xist, and reveals a correlation between XIST/Xist dispersal and autosomal targeting.

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

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

Systematic Analysis of Long Non-Coding RNAs in Inflammasome Activation in Monocytes/Macrophages

The NLRP3 inflammasome plays a pivotal role in regulating inflammation and immune responses. Its activation can lead to an inflammatory response and pyroptotic cell death. This is beneficial in the case of infections, but excessive activation can lead to chronic inflammation and tissue damage. Moreover, while most of the mammalian genome is transcribed as RNAs, only a small fraction codes for proteins. Among non-protein-coding RNAs, long non-coding RNAs (lncRNAs) have been shown to play key roles in regulating gene expression and cellular processes. They interact with DNA, RNAs, and proteins, and their dysregulation can provide insights into disease mechanisms, including NLRP3 inflammasome activation. Here, we systematically analyzed previously published RNA sequencing (RNA-seq) data of NLRP3 inflammasome activation in monocytes/macrophages to uncover inflammasome-regulated lncRNA genes. To uncover the functional importance of inflammasome-regulated lncRNA genes, one inflammasome-regulated lncRNA, ENSG00000273124, was knocked down in an in vitro model of macrophage polarization. The results indicate that silencing of ENSG00000273124 resulted in the up-regulation tumor necrosis factor (TNF), suggesting that this lncRNA might be involved in pro-inflammatory response in macrophages. To make our analyzed data more accessible, we developed the web database InflammasomeDB.

The brain's dark transcriptome: Sequencing RNA in distal compartments of neurons and glia

Transcriptomic approaches are powerful strategies to map the molecular diversity of cells in the brain. Single-cell genomic atlases have now been compiled for entire mammalian brains. However, complementary techniques are only just beginning to map the subcellular transcriptomes from distal cellular compartments. We review single-cell datasets alongside subtranscriptome data from the mammalian brain to explore the development of cellular and subcellular diversity. We discuss how single-cell RNA-seq misses transcripts localized away from cell bodies, which form the 'dark transcriptome' of the brain: a collection of subtranscriptomes in dendrites, axons, growth cones, synapses, and endfeet with important roles in brain development and function. Recent advances in subcellular transcriptome sequencing are beginning to reveal these elusive pools of RNA. We outline the success stories to date in uncovering the constituent subtranscriptomes of neurons and glia, as well as present the emerging toolkit that is accelerating the pace of subtranscriptome discovery.

Pre-mRNA splicing order is predetermined and maintains splicing fidelity across multi-intronic transcripts

Combinatorially, intron excision within a given nascent transcript could proceed down any of thousands of paths, each of which would expose different dynamic landscapes of cis-elements and contribute to alternative splicing. In this study, we found that post-transcriptional multi-intron splicing order in human cells is largely predetermined, with most genes spliced in one or a few predominant orders. Strikingly, these orders were conserved across cell types and stages of motor neuron differentiation. Introns flanking alternatively spliced exons were frequently excised last, after their neighboring introns. Perturbations to the spliceosomal U2 snRNA altered the preferred splicing order of many genes, and these alterations were associated with the retention of other introns in the same transcript. In one gene, early removal of specific introns was sufficient to induce delayed excision of three proximal introns, and this delay was caused by two distinct cis-regulatory mechanisms. Together, our results demonstrate that multi-intron splicing order in human cells is predetermined, is influenced by a component of the spliceosome and ensures splicing fidelity across long pre-mRNAs.

Bi-allelic SNAPC4 variants dysregulate global alternative splicing and lead to neuroregression and progressive spastic paraparesis

The vast majority of human genes encode multiple isoforms through alternative splicing, and the temporal and spatial regulation of those isoforms is critical for organismal development and function. The spliceosome, which regulates and executes splicing reactions, is primarily composed of small nuclear ribonucleoproteins (snRNPs) that consist of small nuclear RNAs (snRNAs) and protein subunits. snRNA gene transcription is initiated by the snRNA-activating protein complex (SNAPc). Here, we report ten individuals, from eight families, with bi-allelic, deleterious SNAPC4 variants. SNAPC4 encoded one of the five SNAPc subunits that is critical for DNA binding. Most affected individuals presented with delayed motor development and developmental regression after the first year of life, followed by progressive spasticity that led to gait alterations, paraparesis, and oromotor dysfunction. Most individuals had cerebral, cerebellar, or basal ganglia volume loss by brain MRI. In the available cells from affected individuals, SNAPC4 abundance was decreased compared to unaffected controls, suggesting that the bi-allelic variants affect SNAPC4 accumulation. The depletion of SNAPC4 levels in HeLa cell lines via genomic editing led to decreased snRNA expression and global dysregulation of alternative splicing. Analysis of available fibroblasts from affected individuals showed decreased snRNA expression and global dysregulation of alternative splicing compared to unaffected cells. Altogether, these data suggest that these bi-allelic SNAPC4 variants result in loss of function and underlie the neuroregression and progressive spasticity in these affected individuals.

NONO enhances mRNA processing of super-enhancer-associated GATA2 and HAND2 genes in neuroblastoma

High-risk neuroblastoma patients have poor survival rates and require better therapeutic options. High expression of a multifunctional DNA and RNA-binding protein, NONO, in neuroblastoma is associated with poor patient outcome; however, there is little understanding of the mechanism of NONO-dependent oncogenic gene regulatory activity in neuroblastoma. Here, we used cell imaging, biochemical and genome-wide molecular analysis to reveal complex NONO-dependent regulation of gene expression. NONO forms RNA- and DNA-tethered condensates throughout the nucleus and undergoes phase separation in vitro, modulated by nucleic acid binding. CLIP analyses show that NONO mainly binds to the 5' end of pre-mRNAs and modulates pre-mRNA processing, dependent on its RNA-binding activity. NONO regulates super-enhancer-associated genes, including HAND2 and GATA2. Abrogating NONO RNA binding, or phase separation activity, results in decreased expression of HAND2 and GATA2. Thus, future development of agents that target RNA-binding activity of NONO may have therapeutic potential in this cancer context.

ESPRESSO: Robust discovery and quantification of transcript isoforms from error-prone long-read RNA-seq data

Long-read RNA sequencing (RNA-seq) holds great potential for characterizing transcriptome variation and full-length transcript isoforms, but the relatively high error rate of current long-read sequencing platforms poses a major challenge. We present ESPRESSO, a computational tool for robust discovery and quantification of transcript isoforms from error-prone long reads. ESPRESSO jointly considers alignments of all long reads aligned to a gene and uses error profiles of individual reads to improve the identification of splice junctions and the discovery of their corresponding transcript isoforms. On both a synthetic spike-in RNA sample and human RNA samples, ESPRESSO outperforms multiple contemporary tools in not only transcript isoform discovery but also transcript isoform quantification. In total, we generated and analyzed ~1.1 billion nanopore RNA-seq reads covering 30 human tissue samples and three human cell lines. ESPRESSO and its companion dataset provide a useful resource for studying the RNA repertoire of eukaryotic transcriptomes.

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.

Kinetics of mRNA nuclear export regulate innate immune response gene expression

The abundance and stimulus-responsiveness of mature mRNA is thought to be determined by nuclear synthesis, processing, and cytoplasmic decay. However, the rate and efficiency of moving mRNA to the cytoplasm almost certainly contributes, but has rarely been measured. Here, we investigated mRNA export rates for innate immune genes. We generated high spatio-temporal resolution RNA-seq data from endotoxin-stimulated macrophages and parameterized a mathematical model to infer kinetic parameters with confidence intervals. We find that the effective chromatin-to-cytoplasm export rate is gene-specific, varying 100-fold: for some genes, less than 5% of synthesized transcripts arrive in the cytoplasm as mature mRNAs, while others show high export efficiency. Interestingly, effective export rates do not determine temporal gene responsiveness, but complement the wide range of mRNA decay rates; this ensures similar abundances of short- and long-lived mRNAs, which form successive innate immune response expression waves.

Spatiotemporal and genetic regulation of A-to-I editing throughout human brain development

Posttranscriptional RNA modifications by adenosine-to-inosine (A-to-I) editing are abundant in the brain, yet elucidating functional sites remains challenging. To bridge this gap, we investigate spatiotemporal and genetically regulated A-to-I editing sites across prenatal and postnatal stages of human brain development. More than 10,000 spatiotemporally regulated A-to-I sites were identified that occur predominately in 3' UTRs and introns, as well as 37 sites that recode amino acids in protein coding regions with precise changes in editing levels across development. Hyper-edited transcripts are also enriched in the aging brain and stabilize RNA secondary structures. These features are conserved in murine and non-human primate models of neurodevelopment. Finally, thousands of cis-editing quantitative trait loci (edQTLs) were identified with unique regulatory effects during prenatal and postnatal development. Collectively, this work offers a resolved atlas linking spatiotemporal variation in editing levels to genetic regulatory effects throughout distinct stages of brain maturation.

Stimulus-specific remodeling of the neuronal transcriptome through nuclear intron-retaining transcripts

The nuclear envelope has long been considered primarily a physical barrier separating nuclear and cytosolic contents. More recently, nuclear compartmentalization has been shown to have additional regulatory functions in controlling gene expression. A sizeable proportion of protein-coding mRNAs is more prevalent in the nucleus than in the cytosol, suggesting regulated mRNA trafficking to the cytosol, but the mechanisms underlying controlled nuclear mRNA retention remain unclear. Here, we provide a comprehensive map of the subcellular localization of mRNAs in mature mouse cortical neurons, and reveal that transcripts retained in the nucleus comprise the majority of stable intron-retaining mRNAs. Systematically probing the fate of nuclear transcripts upon neuronal stimulation, we found opposite effects on sub-populations of transcripts: while some are targeted for degradation, others complete splicing to generate fully mature mRNAs that are exported to the cytosol and mediate rapid increases in protein levels. Finally, different forms of stimulation mobilize distinct groups of intron-retaining transcripts, with this selectivity arising from the activation of specific signaling pathways. Overall, our findings uncover a cue-specific control of intron retention as a major regulator of acute remodeling of the neuronal transcriptome.

A lifelong duty: how Xist maintains the inactive X chromosome

Female eutherians transcriptionally silence one X chromosome to balance gene dosage between the sexes. X-chromosome inactivation (XCI) is initiated by the lncRNA Xist, which assembles many proteins within the inactive X chromosome (Xi) to trigger gene silencing and heterochromatin formation. It is well established that gene silencing on the Xi is maintained through repressive epigenetic processes, including histone deacetylation and DNA methylation. Recent studies revealed a new mechanism where RNA-binding proteins that interact directly with the RNA contribute to the maintenance of Xist localization and gene silencing. In addition, a surprising plasticity of the Xi was uncovered with many genes becoming upregulated upon experimental deletion of Xist. Intriguingly, immune cells normally lose Xist from the Xi, suggesting that thisXist dependence is utilized in vivo to dynamically regulate gene expression from the Xi. These new studies expose fundamental regulatory mechanisms for the chromatin association of RNAs, highlight the need for studying the maintenance of XCI and Xist localization in a gene- and cell-type-specific manner, and are likely to have clinical impact.

Native RNA sequencing in fission yeast reveals frequent alternative splicing isoforms

The unicellular yeast Schizosaccharomyces pombe (fission yeast) retains many of the splicing features observed in humans and is thus an excellent model to study the basic mechanisms of splicing. Nearly half the genes contain introns, but the impact of alternative splicing in gene regulation and proteome diversification remains largely unexplored. Here we leverage Oxford Nanopore Technologies native RNA sequencing (dRNA), as well as ribosome profiling data, to uncover the full range of polyadenylated transcripts and translated open reading frames. We identify 332 alternative isoforms affecting the coding sequences of 262 different genes, 97 of which occur at frequencies higher than 20%, indicating that functional alternative splicing in S. pombe is more prevalent than previously suspected. Intron retention events make about 80% of the cases; these events may be involved in the regulation of gene expression and, in some cases, generate novel protein isoforms, as supported by ribosome profiling data in 18 of the intron retention isoforms. One example is the rpl22 gene, in which intron retention is associated with the translation of a protein of only 13 amino acids. We also find that lowly expressed transcripts tend to have longer poly(A) tails than highly expressed transcripts, highlighting an interdependence between poly(A) tail length and transcript expression level. Finally, we discover 214 novel transcripts that are not annotated, including 158 antisense transcripts, some of which also show translation evidence. The methodologies described in this work open new opportunities to study the regulation of splicing in a simple eukaryotic model.

Type I and II PRMTs inversely regulate post-transcriptional intron detention through Sm and CHTOP methylation

Protein arginine methyltransferases (PRMTs) are required for the regulation of RNA processing factors. Type I PRMT enzymes catalyze mono- and asymmetric dimethylation; Type II enzymes catalyze mono- and symmetric dimethylation. To understand the specific mechanisms of PRMT activity in splicing regulation, we inhibited Type I and II PRMTs and probed their transcriptomic consequences. Using the newly developed Splicing Kinetics and Transcript Elongation Rates by Sequencing (SKaTER-seq) method, analysis of co-transcriptional splicing demonstrated that PRMT inhibition resulted in altered splicing rates. Surprisingly, co-transcriptional splicing kinetics did not correlate with final changes in splicing of polyadenylated RNA. This was particularly true for retained introns (RI). By using actinomycin D to inhibit ongoing transcription, we determined that PRMTs post-transcriptionally regulate RI. Subsequent proteomic analysis of both PRMT-inhibited chromatin and chromatin-associated polyadenylated RNA identified altered binding of many proteins, including the Type I substrate, CHTOP, and the Type II substrate, SmB. Targeted mutagenesis of all methylarginine sites in SmD3, SmB, and SmD1 recapitulated splicing changes seen with Type II PRMT inhibition, without disrupting snRNP assembly. Similarly, mutagenesis of all methylarginine sites in CHTOP recapitulated the splicing changes seen with Type I PRMT inhibition. Examination of subcellular fractions further revealed that RI were enriched in the nucleoplasm and chromatin. Taken together, these data demonstrate that, through Sm and CHTOP arginine methylation, PRMTs regulate the post-transcriptional processing of nuclear, detained introns.

Towards resolution of the intron retention paradox in breast cancer

BACKGROUND

After many years of neglect in the field of alternative splicing, the importance of intron retention (IR) in cancer has come into focus following landmark discoveries of aberrant IR patterns in cancer. Many solid and liquid tumours are associated with drastic increases in IR, and such patterns have been pursued as both biomarkers and therapeutic targets. Paradoxically, breast cancer (BrCa) is the only tumour type in which IR is reduced compared to adjacent normal breast tissue.

METHODS

In this study, we have conducted a pan-cancer analysis of IR with emphasis on BrCa and its subtypes. We explored mechanisms that could cause aberrant and pathological IR and clarified why normal breast tissue has unusually high IR.

RESULTS

Strikingly, we found that aberrantly decreasing IR in BrCa can be largely attributed to normal breast tissue having the highest occurrence of IR events compared to other healthy tissues. Our analyses suggest that low numbers of IR events in breast tumours are associated with poor prognosis, particularly in the luminal B subtype. Interestingly, we found that IR frequencies negatively correlate with cell proliferation in BrCa cells, i.e. rapidly dividing tumour cells have the lowest number of IR events. Aberrant RNA-binding protein expression and changes in tissue composition are among the causes of aberrantly decreasing IR in BrCa.

CONCLUSIONS

Our results suggest that IR should be considered for therapeutic manipulation in BrCa patients with aberrantly low IR levels and that further work is needed to understand the cause and impact of high IR in other tumour types.

Post-Transcriptional Control of Mating-Type Gene Expression during Gametogenesis in Saccharomyces cerevisiae

Gametogenesis in diploid cells of the budding yeast Saccharomyces cerevisiae produces four haploid meiotic products called spores. Spores are dormant until nutrients trigger germination, when they bud asexually or mate to return to the diploid state. Each sporulating diploid produces a mix of spores of two haploid mating types, a and α. In asexually dividing haploids, the mating types result from distinct, mutually exclusive gene expression programs responsible for production of mating pheromones and the receptors to sense them, all of which are silent in diploids. It was assumed that spores only transcribe haploid- and mating-type-specific genes upon germination. We find that dormant spores of each mating type harbor transcripts representing all these genes, with the exception of Mata1, which we found to be enriched in a spores. Mata1 transcripts, from a rare yeast gene with two introns, were mostly unspliced. If the retained introns reflect tethering to the MATa locus, this could provide a mechanism for biased inheritance. Translation of pheromones and receptors were repressed at least until germination. We find antisense transcripts to many mating genes that may be responsible. These findings add to the growing number of examples of post-transcriptional regulation of gene expression during gametogenesis.