A dynamic intron retention program enriched in RNA processing genes regulates gene expression during terminal erythropoiesis

Human erythrocytes' perplexing behaviour: erythrocytic microRNAs

Erythrocytes have the potential role in erythropoiesis and disease diagnosis. Thought to have lacked nucleic acid content, mammalian erythrocytes are nevertheless able to function for 120-140 days, metabolize heme, maintain oxidative stress, and so on. Mysteriously, erythrocytes proved as largest repositories of microRNAs (miRNAs) some of which are selectively retained and function in mature erythrocytes. They have unique expression patterns and have been found to be linked to specific conditions such as sickle cell anaemia, high-altitude hypoxia, chronic mountain sickness, cardiovascular and metabolic conditions as well as host-parasite interactions. They also have been implicated in cell storage-related damage and the regulation of its survival. However, the mechanism by which miRNAs function in the cell remains unclear. Investigations into the molecular mechanism of miRNAs in erythrocytes via extracellular vesicles have provided important clues in research studies on Plasmodium infection. Erythrocytes are also the primary source of circulating miRNAs but, how they affect the plasma/serum miRNAs profiles are still poorly understood. Erythrocyte-derived exosomal miRNAs, can interact with various body cell types, and have easy access to all regions, making them potentially crucial in various pathophysiological conditions. Which can also improve our understanding to identify potential treatment options and discovery related to non-invasive diagnostic markers. This article emphasizes the importance of erythrocytic miRNAs while focusing on the enigmatic behaviour of erythrocytes. It also sheds light on how this knowledge may be applied in the future to enhance the state of erythrocyte translational research from the standpoint of erythrocytic miRNAs.

Aberrant spliceosome activity via elevated intron retention and upregulation and phosphorylation of SF3B1 in chronic lymphocytic leukemia

Splicing factor 3b subunit 1 (SF3B1) is the largest subunit and core component of the spliceosome. Inhibition of SF3B1 was associated with an increase in broad intron retention (IR) on most transcripts, suggesting that IR can be used as a marker of spliceosome inhibition in chronic lymphocytic leukemia (CLL) cells. Furthermore, we separately analyzed exonic and intronic mapped reads on annotated RNA-sequencing transcripts obtained from B cells (n = 98 CLL patients) and healthy volunteers (n = 9). We measured intron/exon ratio to use that as a surrogate for alternative RNA splicing (ARS) and found that 66% of CLL-B cell transcripts had significant IR elevation compared with normal B cells (NBCs) and that correlated with mRNA downregulation and low expression levels. Transcripts with the highest IR levels belonged to biological pathways associated with gene expression and RNA splicing. A >2-fold increase of active pSF3B1 was observed in CLL-B cells compared with NBCs. Additionally, when the CLL-B cells were treated with macrolides (pladienolide-B), a significant decrease in pSF3B1, but not total SF3B1 protein, was observed. These findings suggest that IR/ARS is increased in CLL, which is associated with SF3B1 phosphorylation and susceptibility to SF3B1 inhibitors. These data provide additional support to the relevance of ARS in carcinogenesis and evidence of pSF3B1 participation in this process.

Accelerated DNA replication fork speed due to loss of R-loops in myelodysplastic syndromes with SF3B1 mutation

Myelodysplastic syndromes (MDS) with mutated SF3B1 gene present features including a favourable outcome distinct from MDS with mutations in other splicing factor genes SRSF2 or U2AF1. Molecular bases of these divergences are poorly understood. Here we find that SF3B1-mutated MDS show reduced R-loop formation predominating in gene bodies associated with intron retention reduction, not found in U2AF1- or SRSF2-mutated MDS. Compared to erythroblasts from SRSF2- or U2AF1-mutated patients, SF3B1-mutated erythroblasts exhibit augmented DNA synthesis, accelerated replication forks, and single-stranded DNA exposure upon differentiation. Importantly, histone deacetylase inhibition using vorinostat restores R-loop formation, slows down DNA replication forks and improves SF3B1-mutated erythroblast differentiation. In conclusion, loss of R-loops with associated DNA replication stress represents a hallmark of SF3B1-mutated MDS ineffective erythropoiesis, which could be used as a therapeutic target.

Minor introns impact on hematopoietic malignancies

In the intricate orchestration of the central dogma, pre-mRNA splicing plays a crucial role in the post-transcriptional process that transforms DNA into mature mRNA. Widely acknowledged as a pivotal RNA processing step, it significantly influences gene expression and alters the functionality of gene product proteins. Although U2-dependent spliceosomes efficiently manage the removal of over 99% of introns, a distinct subset of essential genes undergo splicing with a different intron type, denoted as minor introns, using U12-dependent spliceosomes. Mutations in spliceosome component genes are now recognized as prevalent genetic abnormalities in cancer patients, especially those with hematologic malignancies. Despite the relative rarity of minor introns, genes containing them are evolutionarily conserved and play crucial roles in functions such as the RAS-MAPK pathway. Disruptions in U12-type minor intron splicing caused by mutations in snRNA or its regulatory components significantly contribute to cancer progression. Notably, recurrent mutations associated with myelodysplastic syndrome (MDS) in the minor spliceosome component ZRSR2 underscore its significance. Examination of ZRSR2-mutated MDS cells has revealed that only a subset of minor spliceosome-dependent genes, such as LZTR1, consistently exhibit missplicing. Recent technological advancements have uncovered insights into minor introns, raising inquiries beyond current understanding. This review comprehensively explores the importance of minor intron regulation, the molecular implications of minor (U12-type) spliceosomal mutations and cis-regulatory regions, and the evolutionary progress of studies on minor, aiming to provide a sophisticated understanding of their intricate role in cancer biology.

ASAPA: a bioinformatic pipeline based on Iso-Seq that identifies the links among alternative splicing, alternative transcription initiation and alternative polyadenylation

BACKGROUND

Although the events associated with alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI) can be identified by many approaches based on isoform sequencing (Iso-Seq), these analyses are generally independent of each other and the links between these events are still rarely mentioned. However, an interdependency analysis can be achieved because the transcriptional start site, splice sites and polyA site could be simultaneously included in a long, full-length read from Iso-Seq.

RESULTS

We create ASAPA pipeline that enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data. We tested this pipeline using Arabidopsis data and found some interesting results: some adjacent introns tend to be simultaneously spliced or retained; coupling between AS and ATI or APA is limited to the initial or terminal intron; and ATI and APA are potentially linked in some special cases.

CONCLUSION

Our pipeline enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data, which is conducive to a better understanding of transcription landscape generation.

Intron retention is a mechanism of erythropoietin regulation in brain cell models

Intron retention is a mechanism of post-transcriptional gene regulation, including genes involved in erythropoiesis. Erythropoietin (EPO) is a hormone without evidence of intracellular vesicle storage that regulates erythropoiesis. We hypothesize that EPO uses intron retention as a mechanism of post-transcriptional regulation in response to hypoxia and ischemia. Cell models of hypoxia and ischemia for kidney, liver, and brain cells were examined for intron retention by real time quantitative PCR. EPO expression increased in most cells except for blood brain barrier and liver cells. The intron retained transcript ratio decreased in brain cells, except for Astrocytes, but showed no change in kidney or liver after 24 h of ischemia. The shift in intron ratio was maintained when using poly (A) enriched cDNA, suggesting that intron retention is not due to immature transcripts. The expression of EPO was elevated at variable time points amongst cell models with the intron ratio also changing over a time course of 2 to 16 h after ischemia. We conclude that intron retention is a mechanism regulating EPO expression in response to ischemia in a tissue specific manner.

iFLAS: positive-unlabeled learning facilitates full-length transcriptome-based identification and functional exploration of alternatively spliced isoforms in maize

The advent of full-length transcriptome sequencing technologies has accelerated the discovery of novel splicing isoforms. However, existing alternative splicing (AS) tools are either tailored for short-read RNA-Seq data or designed for human and animal studies. The disparities in AS patterns between plants and animals still pose a challenge to the reliable identification and functional exploration of novel isoforms in plants. Here, we developed integrated full-length alternative splicing analysis (iFLAS), a plant-optimized AS toolkit that introduced a semi-supervised machine learning method known as positive-unlabeled (PU) learning to accurately identify novel isoforms. iFLAS also enables the investigation of AS functions from various perspectives, such as differential AS, poly(A) tail length, and allele-specific AS (ASAS) analyses. By applying iFLAS to three full-length transcriptome sequencing datasets, we systematically identified and functionally characterized maize (Zea mays) AS patterns. We found intron retention not only introduces premature termination codons, resulting in lower expression levels of isoforms, but may also regulate the length of 3'UTR and poly(A) tail, thereby affecting the functional differentiation of isoforms. Moreover, we observed distinct ASAS patterns in two genes within heterosis offspring, highlighting their potential value in breeding. These results underscore the broad applicability of iFLAS in plant full-length transcriptome-based AS research.

Replicative aging in yeast involves dynamic intron retention patterns associated with mRNA processing/export and protein ubiquitination

Saccharomyces cerevisiae (baker's yeast) has yielded relevant insights into some of the basic mechanisms of organismal aging. Among these are genomic instability, oxidative stress, caloric restriction and mitochondrial dysfunction. Several genes are known to have an impact on the aging process, with corresponding mutants exhibiting short- or long-lived phenotypes. Research dedicated to unraveling the underlying cellular mechanisms can support the identification of conserved mechanisms of aging in other species. One of the hitherto less studied fields in yeast aging is how the organism regulates its gene expression at the transcriptional level. To our knowledge, we present the first investigation into alternative splicing, particularly intron retention, during replicative aging of S. cerevisiae. This was achieved by utilizing the IRFinder algorithm on a previously published RNA-seq data set by Janssens et al. (2015). In the present work, 44 differentially retained introns in 43 genes were identified during replicative aging. We found that genes with altered intron retention do not display significant changes in overall transcript levels. It was possible to functionally assign distinct groups of these genes to the cellular processes of mRNA processing and export (e.g., YRA1) in early and middle-aged yeast, and protein ubiquitination (e.g., UBC5) in older cells. In summary, our work uncovers a previously unexplored layer of the transcriptional program of yeast aging and, more generally, expands the knowledge on the occurrence of alternative splicing in baker's yeast.

Exitrons: offering new roles to retained introns-the novel regulators of protein diversity and utility

Exitrons are exonic introns. This subclass of intron retention alternative splicing does not contain a Pre-Terminating stop Codon. Therefore, when retained, they are always a part of a protein. Intron retention is a frequent phenomenon predominantly found in plants, which results in either the degradation of the transcripts or can serve as a stable intermediate to be processed upon induction by specific signals or the cell status. Interestingly, exitrons have coding ability and may confer additional attributes to the proteins that retain them. Therefore, exitron-containing and exitron-spliced isoforms will be a driving force for creating protein diversity in the proteome of an organism. This review establishes a basic understanding of exitron, discussing its genesis, key features, identification methods and functions. We also try to depict its other potential roles. The present review also aims to provide a fundamental background to those who found such exitronic sequences in their gene(s) and to speculate the future course of studies.

Noncoding rules of survival: epigenetic regulation of normal and malignant hematopoiesis

Hematopoiesis is an essential process for organismal development and homeostasis. Epigenetic regulation of gene expression is critical for stem cell self-renewal and differentiation in normal hematopoiesis. Increasing evidence shows that disrupting the balance between self-renewal and cell fate decisions can give rise to hematological diseases such as bone marrow failure and leukemia. Consequently, next-generation sequencing studies have identified various aberrations in histone modifications, DNA methylation, RNA splicing, and RNA modifications in hematologic diseases. Favorable outcomes after targeting epigenetic regulators during disease states have further emphasized their importance in hematological malignancy. However, these targeted therapies are only effective in some patients, suggesting that further research is needed to decipher the complexity of epigenetic regulation during hematopoiesis. In this review, an update on the impact of the epigenome on normal hematopoiesis, disease initiation and progression, and current therapeutic advancements will be discussed.

The transcriptomic landscape of normal and ineffective erythropoiesis at single-cell resolution

The anemias of myelodysplastic syndrome (MDS) and Diamond Blackfan anemia (DBA) are generally macrocytic and always reflect ineffective erythropoiesis yet result from diverse genetic mutations. To delineate shared mechanisms that lead to cell death, we studied the fate of single erythroid marrow cells from individuals with DBA or MDS-5q. We defined an unhealthy (vs healthy) differentiation trajectory using transcriptional pseudotime and cell surface proteins. The pseudotime trajectories diverge immediately after cells upregulate transferrin receptor (CD71), import iron, and initiate heme synthesis, although cell death occurs much later. Cells destined to die express high levels of heme-responsive genes, including ribosomal protein and globin genes, whereas surviving cells downregulate heme synthesis and upregulate DNA damage response, hypoxia, and HIF1 pathways. Surprisingly, 24% ± 12% of cells from control subjects follow the unhealthy trajectory, implying that heme might serve as a rheostat directing cells to live or die. When heme synthesis was inhibited with succinylacetone, more DBA cells followed the healthy trajectory and survived. We also noted high numbers of messages with retained introns that increased as erythroid cells matured, confirmed the rapid cycling of colony forming unit-erythroid, and demonstrated that cell cycle timing is an invariant property of differentiation stage. Including unspliced RNA in pseudotime determinations allowed us to reliably align independent data sets and accurately query stage-specific transcriptomic changes. MDS-5q (unlike DBA) results from somatic mutation, so many normal (unmutated) erythroid cells persist. By independently tracking erythroid differentiation of cells with and without chromosome 5q deletions, we gained insight into why 5q+ cells cannot expand to prevent anemia.

Posttranslational splicing modifications as a key mechanism in cytarabine resistance in acute myeloid leukemia

Despite the approval of several drugs for AML, cytarabine is still widely used as a therapeutic approach. However, 85% of patients show resistance and only 10% overcome the disease. Using RNA-seq and phosphoproteomics, we show that RNA splicing and serine-arginine-rich (SR) proteins phosphorylation were altered during cytarabine resistance. Moreover, phosphorylation of SR proteins at diagnosis were significantly lower in responder than non-responder patients, pointing to their utility to predict response. These changes correlated with altered transcriptomic profiles of SR protein target genes. Notably, splicing inhibitors were therapeutically effective in treating sensitive and resistant AML cells as monotherapy or combination with other approved drugs. H3B-8800 and venetoclax combination showed the best efficacy in vitro, demonstrating synergistic effects in patient samples and no toxicity in healthy hematopoietic progenitors. Our results establish that RNA splicing inhibition, alone or combined with venetoclax, could be useful for the treatment of newly diagnosed or relapsed/refractory AML.

Post-transcriptional regulation of erythropoiesis

Erythropoiesis is a complex, precise, and lifelong process that is essential for maintaining normal body functions. Its strict regulation is necessary to prevent a variety of blood diseases. Normal erythropoiesis is precisely regulated by an intricate network that involves transcription levels, signal transduction, and various epigenetic modifications. In recent years, research on post-transcriptional levels in erythropoiesis has expanded significantly. The dynamic regulation of splicing transitions is responsible for changes in protein isoform expression that add new functions beneficial for erythropoiesis. RNA-binding proteins adapt the translation of transcripts to the protein requirements of the cell, yielding mRNA with dynamic translation efficiency. Noncoding RNAs, such as microRNAs and lncRNAs, are indispensable for changing the translational efficiency and/or stability of targeted mRNAs to maintain the normal expression of genes related to erythropoiesis. N6-methyladenosine-dependent regulation of mRNA translation plays an important role in maintaining the expression programs of erythroid-related genes and promoting erythroid lineage determination. This review aims to describe our current understanding of the role of post-transcriptional regulation in erythropoiesis and erythroid-associated diseases, and to shed light on the physiological and pathological implications of the post-transcriptional regulation machinery in erythropoiesis. These may help to further enrich our understanding of the regulatory network of erythropoiesis and provide new strategies for the diagnosis and treatment of erythroid-related diseases.

Biophysical modeling with variational autoencoders for bimodal, single-cell RNA sequencing data

We motivate and present biVI, which combines the variational autoencoder framework of scVI with biophysically motivated, bivariate models for nascent and mature RNA distributions. While previous approaches to integrate bimodal data via the variational autoencoder framework ignore the causal relationship between measurements, biVI models the biophysical processes that give rise to observations. We demonstrate through simulated benchmarking that biVI captures cell type structure in a low-dimensional space and accurately recapitulates parameter values and copy number distributions. On biological data, biVI provides a scalable route for identifying the biophysical mechanisms underlying gene expression. This analytical approach outlines a generalizable strategy for treating multimodal datasets generated by high-throughput, single-cell genomic assays.

Abnormalities in intron retention characterize patients with systemic lupus erythematosus

Regulation of intron retention (IR), a form of alternative splicing, is a newly recognized checkpoint in gene expression. Since there are numerous abnormalities in gene expression in the prototypic autoimmune disease systemic lupus erythematosus (SLE), we sought to determine whether IR was intact in patients with this disease. We, therefore, studied global gene expression and IR patterns of lymphocytes in SLE patients. We analyzed RNA-seq data from peripheral blood T cell samples from 14 patients suffering from systemic lupus erythematosus (SLE) and 4 healthy controls and a second, independent data set of RNA-seq data from B cells from16 SLE patients and 4 healthy controls. We identified intron retention levels from 26,372 well annotated genes as well as differential gene expression and tested for differences between cases and controls using unbiased hierarchical clustering and principal component analysis. We followed with gene-disease enrichment analysis and gene-ontology enrichment analysis. Finally, we then tested for significant differences in intron retention between cases and controls both globally and with respect to specific genes. Overall decreased IR was found in T cells from one cohort and B cells from another cohort of patients with SLE and was associated with increased expression of numerous genes, including those encoding spliceosome components. Different introns within the same gene displayed both up- and down-regulated retention profiles indicating a complex regulatory mechanism. These results indicate that decreased IR in immune cells is characteristic of patients with active SLE and may contribute to the abnormal expression of specific genes in this autoimmune disease.

SUMOylation regulates pre-mRNA splicing to overcome DNA damage in fungi

Pathogenic fungi are subject to DNA damage stress derived from host immune responses during infection. Small ubiquitin-like modifier (SUMO) modification and precursor (pre)-mRNA splicing are both involved in DNA damage response (DDR). However, the mechanisms of how SUMOylation and splicing coordinated in DDR remain largely unknown. Combining with biochemical analysis, RNA-Seq method, and biological analysis, we report that SUMO pathway participates in DDR and virulence in Fusarium graminearum, a causal agent of Fusarium head blight of cereal crops world-wide. Interestingly, a key transcription factor FgSR is SUMOylated upon DNA damage stress. SUMOylation regulates FgSR nuclear-cytoplasmic partitioning and its phosphorylation by FgMec1, and promotes its interaction with chromatin remodeling complex SWI/SNF for activating the expression of DDR-related genes. Moreover, the SWI/SNF complex was found to further recruit splicing-related NineTeen Complex, subsequently modulates pre-mRNA splicing during DDR. Our findings reveal a novel function of SUMOylation in DDR by regulating a transcription factor to orchestrate gene expression and pre-mRNA splicing to overcome DNA damage during the infection of F. graminearum, which advances the understanding of the delicate regulation of DDR by SUMOylation in pathogenic fungi, and extends the knowledge of cooperation of SUMOylation and pre-mRNA splicing in DDR in eukaryotes.

RNA splicing factors in normal hematopoiesis and hematologic malignancies: novel therapeutic targets and strategies

RNA splicing, a crucial transesterification-based process by which noncoding regions are removed from premature RNA to create mature mRNA, regulates various cellular functions, such as proliferation, survival, and differentiation. Clinical and functional studies over the past 10 y have confirmed that mutations in RNA splicing factors are among the most recurrent genetic abnormalities in hematologic neoplasms, including myeloid malignancies, chronic lymphocytic leukemia, mantle cell lymphoma, and clonal hematopoiesis. These findings indicate an important role for splicing factor mutations in the development of clonal hematopoietic disorders. Mutations in core or accessory components of the RNA spliceosome complex alter splicing sites in a manner of change of function. These changes can result in the dysregulation of cancer-associated gene expression and the generation of novel mRNA transcripts, some of which are not only critical to disease development but may be also serving as potential therapeutic targets. Furthermore, multiple studies have revealed that hematopoietic cells bearing mutations in splicing factors depend on the expression of the residual wild-type allele for survival, and these cells are more sensitive to reduced expression of wild-type splicing factors or chemical perturbations of the splicing machinery. These findings suggest a promising possibility for developing novel therapeutic opportunities in tumor cells based on mutations in splicing factors. Here, we combine current knowledge of the mechanistic and functional effects of frequently mutated splicing factors in normal hematopoiesis and the effects of their mutations in hematologic malignancies. Moreover, we discuss the development of potential therapeutic opportunities based on these mutations.

LAFITE Reveals the Complexity of Transcript Isoforms in Subcellular Fractions

Characterization of the subcellular distribution of RNA is essential for understanding the molecular basis of biological processes. Here, the subcellular nanopore direct RNA-sequencing (DRS) of four lung cancer cell lines (A549, H1975, H358, and HCC4006) is performed, coupled with a computational pipeline, Low-abundance Aware Full-length Isoform clusTEr (LAFITE), to comprehensively analyze the full-length cytoplasmic and nuclear transcriptome. Using additional DRS and orthogonal data sets, it is shown that LAFITE outperforms current methods for detecting full-length transcripts, particularly for low-abundance isoforms that are usually overlooked due to poor read coverage. Experimental validation of six novel isoforms exclusively identified by LAFITE further confirms the reliability of this pipeline. By applying LAFITE to subcellular DRS data, the complexity of the nuclear transcriptome is revealed in terms of isoform diversity, 3'-UTR usage, m6A modification patterns, and intron retention. Overall, LAFITE provides enhanced full-length isoform identification and enables a high-resolution view of the RNA landscape at the isoform level.

Reversal of splicing infidelity is a pre-activation step in B cell differentiation

Introduction

B cell activation and differentiation is central to the adaptive immune response. Changes in exon usage can have major impacts on cellular signaling and differentiation but have not been systematically explored in differentiating B cells.

Methods

We analyzed exon usage and intron retention in RNA-Seq data from subsets of human B cells at various stages of differentiation, and in an in vitro laboratory model of B cell activation and differentiation (Epstein Barr virus infection).

Results

Blood naïve B cells were found to have an unusual splicing profile, with unannotated splicing events in over 30% of expressed genes. Splicing changed substantially upon naïve B cell entry into secondary lymphoid tissue and before activation, involving significant increases in exon commitment and reductions in intron retention. These changes preferentially involved short introns with weak splice sites and were likely mediated by an overall increase in splicing efficiency induced by the lymphoid environment. The majority of transcripts affected by splicing changes showed restoration of encoded conserved protein domains and/or reduced targeting to the nonsense-mediated decay pathway. Affected genes were enriched in functionally important immune cell activation pathways such as antigen-mediated signaling, cell cycle control and mRNA processing and splicing.

Discussion

Functional observations from donor B cell subsets in progressive states of differentiation and from timecourse experiments using the in vitro model suggest that these widespread changes in mRNA splicing play a role in preparing naïve B cells for the decisive step of antigen-mediated activation and differentiation.

Identification and characterization of RBM12 as a novel regulator of fetal hemoglobin expression

The fetal-to-adult hemoglobin transition is clinically relevant because reactivation of fetal hemoglobin (HbF) significantly reduces morbidity and mortality associated with sickle cell disease (SCD) and β-thalassemia. Most studies on the developmental regulation of the globin genes, including genome-wide genetics screens, have focused on DNA binding proteins, including BCL11A and ZBTB7A/LRF and their cofactors. Our understanding of RNA binding proteins (RBPs) in this process is much more limited. Two RBPs, LIN28B and IGF2BP1, are known posttranscriptional regulators of HbF production, but a global view of RBPs is still lacking. Here, we carried out a CRISPR/Cas9-based screen targeting RBPs harboring RNA methyltransferase and/or RNA recognition motif (RRM) domains and identified RNA binding motif 12 (RBM12) as a novel HbF suppressor. Depletion of RBM12 induced HbF expression and attenuated cell sickling in erythroid cells derived from patients with SCD with minimal detrimental effects on cell maturation. Transcriptome and proteome profiling revealed that RBM12 functions independently of major known HbF regulators. Enhanced cross-linking and immunoprecipitation followed by high-throughput sequencing revealed strong preferential binding of RBM12 to 5' untranslated regions of transcripts, narrowing down the mechanism of RBM12 action. Notably, we pinpointed the first of 5 RRM domains as essential, and, in conjunction with a linker domain, sufficient for RBM12-mediated HbF regulation. Our characterization of RBM12 as a negative regulator of HbF points to an additional regulatory layer of the fetal-to-adult hemoglobin switch and broadens the pool of potential therapeutic targets for SCD and β-thalassemia.

Retained introns in long RNA-seq reads are not reliably detected in sample-matched short reads

BACKGROUND

There is growing interest in retained introns in a variety of disease contexts including cancer and aging. Many software tools have been developed to detect retained introns from short RNA-seq reads, but reliable detection is complicated by overlapping genes and transcripts as well as the presence of unprocessed or partially processed RNAs.

RESULTS

We compared introns detected by 8 tools using short RNA-seq reads with introns observed in long RNA-seq reads from the same biological specimens. We found significant disagreement among tools (Fleiss' [Formula: see text]) such that 47.7% of all detected intron retentions were not called by more than one tool. We also observed poor performance of all tools, with none achieving an F1-score greater than 0.26, and qualitatively different behaviors between general-purpose alternative splicing detection tools and tools confined to retained intron detection.

CONCLUSIONS

Short-read tools detect intron retention with poor recall and precision, calling into question the completeness and validity of a large percentage of putatively retained introns called by commonly used methods.

Physiological intron retaining transcripts in the cytoplasm abound during human motor neurogenesis

Intron retention (IR) is now recognized as a dominant splicing event during motor neuron (MN) development; however, the role and regulation of intron-retaining transcripts (IRTs) localized to the cytoplasm remain particularly understudied. Here we show that IR is a physiological process that is spatiotemporally regulated during MN lineage restriction and that IRTs in the cytoplasm are detected in as many as 13% (n = 2297) of the genes expressed during this process. We identify a major class of cytoplasmic IRTs that are not associated with reduced expression of their own genes but instead show a high capacity for RNA-binding protein and miRNA occupancy. Finally, we show that ALS-causing VCP mutations lead to a selective increase in cytoplasmic abundance of this particular class of IRTs, which in turn temporally coincides with an increase in the nuclear expression level of predicted miRNA target genes. Altogether, our study identifies a previously unrecognized class of cytoplasmic intronic sequences with potential regulatory function beyond gene expression.

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.

RNA velocity unraveled

We perform a thorough analysis of RNA velocity methods, with a view towards understanding the suitability of the various assumptions underlying popular implementations. In addition to providing a self-contained exposition of the underlying mathematics, we undertake simulations and perform controlled experiments on biological datasets to assess workflow sensitivity to parameter choices and underlying biology. Finally, we argue for a more rigorous approach to RNA velocity, and present a framework for Markovian analysis that points to directions for improvement and mitigation of current problems.

Single-cell sequencing data are snapshots of biological processes, making it challenging to infer dynamic relationships between cell types. RNA velocity attempts to bypass this challenge by treating the unspliced RNA content as a proxy for spliced RNA content in the near future, and using this “extrapolation” to build directional relationships. However, the method, as implemented in several software packages, is not yet reliable enough to be actionable, in part due to the large number of arbitrary, user-set hyperparameters, as well as fundamental incompatibilities between the biophysics of transcription in the living cell and the models used throughout the velocity workflows. In this study, we review these issues, and use existing results from the fields of stochastic modeling and fluorescence transcriptomics to develop an alternative theoretical framework. We show that our framework can facilitate the development and inference of physically consistent models for sequencing data, as well as the unification of single-cell analyses to self-consistently treat variation due to cell type dynamics and identities, the stochasticity inherent to single-molecule processes, and the uncertainty introduced by sequencing experiments.

Dwarf and High Tillering1 represses rice tillering through mediating the splicing of D14 pre-mRNA

Strigolactones (SLs) constitute a class of plant hormones that regulate many aspects of plant development, including repressing tillering in rice (Oryza sativa). However, how SL pathways are regulated is still poorly understood. Here, we describe a rice mutant dwarf and high tillering1 (dht1), which exhibits pleiotropic phenotypes (such as dwarfism and increased tiller numbers) similar to those of mutants defective in SL signaling. We show that DHT1 encodes a monocotyledon-specific hnRNP-like protein that acts as a previously unrecognized intron splicing factor for many precursor mRNAs (pre-mRNAs), including for the SL receptor gene D14. We find that the dht1 (DHT1I232F) mutant protein is impaired in its stability and RNA binding activity, causing defective splicing of D14 pre-mRNA and reduced D14 expression, and consequently leading to the SL signaling-defective phenotypes. Overall, our findings deepen our understanding of the functional diversification of hnRNP-like proteins and establish a connection between posttranscriptional splicing and SL signaling in the regulation of plant development.

Intron retention is a stress response in sensor genes and is restored by Japanese herbal medicines: A basis for future clinical applications

Intron retention (IR) is a regulatory mechanism that can retard protein production by acting at the level of mRNA processing. We recently demonstrated that IR occurs at the pre-symptomatic state during the aging process of a mouse model of aging, providing a promising biomarker for that state, and can be restored to the normal state by juzentaihoto (JTT), a Japanese herbal medicine (Kampo) (Okada et al. 2021). Here we characterized the genes that accumulate retained introns, examined the biological significance of increased IR in these genes for the host, and determined whether drugs other than JTT can have this effect. By analyzing RNA-sequencing data generated from the hippocampus of the 19-week-old SAMP8 mouse, a model for studying age-related depression and Alzheimer's disease, we showed that genes with increased IR are generally involved in multiple metabolic pathways and have pivotal roles in sensing homeostasis. We thus propose that IR is a stress response and works to fine-tune the expression of many downstream target genes, leading to lower levels of their translation under stress conditions. Interestingly, Kampo medicines, as well as other organic compounds, restored splicing of a specific set of retained introns in these sensor genes in accordance with the physiological recovery conditions of the host, which corresponds with the recovery of transcripts represented by differentially expressed genes. Thus, analysis of IR genes may have broad applicability in evaluating the pre-symptomatic state based on the extent of IR of selective sensor genes, opening a promising early diagnosis of any diseases and a strategy for evaluating efficacies of several drugs based on the extent of IR restoration of these sensor genes.

Dysregulated minor intron splicing in cancer

Pre‐mRNA splicing is now widely recognized as a cotranscriptional and post‐transcriptional mechanism essential for regulating gene expression and modifying gene product function. Mutations in genes encoding core spliceosomal proteins and accessory regulatory splicing factors are now considered among the most recurrent genetic abnormalities in patients with cancer, particularly hematologic malignancies. These include mutations in the major (U2‐type) and minor (U12‐type) spliceosomes, which remove >99% and ~0.35% of introns, respectively. Growing evidence indicates that aberrant splicing of evolutionarily conserved U12‐type minor introns plays a crucial role in cancer as the minor spliceosome component, ZRSR2, is subject to recurrent, leukemia‐associated mutations, and intronic mutations have been shown to disrupt the splicing of minor introns. Here, we review the importance of minor intron regulation, the molecular effects of the minor (U12‐type) spliceosomal mutations and cis‐regulatory regions, and the development of minor intron studies for better understanding of cancer biology.

Meta-Analysis Suggests That Intron Retention Can Affect Quantification of Transposable Elements from RNA-Seq Data

Simple Summary

Transposable elements (TEs) are repetitive sequences comprising more than one third of the human genome with the original ability to change their location within the genome. Owing to their repetitive nature, the quantification of TEs results often challenging. RNA-seq is a useful tool for genome-wide TEs quantification, nevertheless it also presents technical issues, including low reads mappability and erroneous quantification derived from the transcription of TEs fragments embedded in canonical transcripts. Fragments derived from TEs are found within the introns of most genes, which led to the hypothesis that intron retention (IR) can affect the unbiased quantification of TEs expression. Performing meta-analysis of public RNA-seq datasets, here we observe that IR can indeed impact the quantification of TEs by increasing the number of reads mapped on intronic TE copies. Our work highlights a correlation between IR and TEs expression measurement by RNA-seq that should be taken into account to achieve reliable TEs quantification, especially in samples characterized by extensive IR, because differential IR might be confused with differential TEs expression.

Abstract

Transposable elements (TEs), also known as “jumping genes”, are repetitive sequences with the capability of changing their location within the genome. They are key players in many different biological processes in health and disease. Therefore, a reliable quantification of their expression as transcriptional units is crucial to distinguish between their independent expression and the transcription of their sequences as part of canonical transcripts. TEs quantification faces difficulties of different types, the most important one being low reads mappability due to their repetitive nature preventing an unambiguous mapping of reads originating from their sequences. A large fraction of TEs fragments localizes within introns, which led to the hypothesis that intron retention (IR) can be an additional source of bias, potentially affecting accurate TEs quantification. IR occurs when introns, normally removed from the mature transcript by the splicing machinery, are maintained in mature transcripts. IR is a widespread mechanism affecting many different genes with cell type-specific patterns. We hypothesized that, in an RNA-seq experiment, reads derived from retained introns can introduce a bias in the detection of overlapping, independent TEs RNA expression. In this study we performed meta-analysis using public RNA-seq data from lymphoblastoid cell lines and show that IR can impact TEs quantification using established tools with default parameters. Reads mapped on intronic TEs were indeed associated to the expression of TEs and influence their correct quantification as independent transcriptional units. We confirmed these results using additional independent datasets, demonstrating that this bias does not appear in samples where IR is not present and that differential TEs expression does not impact on IR quantification. We concluded that IR causes the over-quantification of intronic TEs and differential IR might be confused with differential TEs expression. Our results should be taken into account for a correct quantification of TEs expression from RNA-seq data, especially in samples in which IR is abundant.

The RNA-binding protein SRSF3 has an essential role in megakaryocyte maturation and platelet production

RNA processing is increasingly recognized as a critical control point in the regulation of different hematopoietic lineages including megakaryocytes responsible for the production of platelets. Platelets are anucleate cytoplasts that contain a rich repertoire of RNAs encoding proteins with essential platelet functions derived from the parent megakaryocyte. It is largely unknown how RNA binding proteins contribute to the development and functions of megakaryocytes and platelets. We show that serine-arginine-rich splicing factor 3 (SRSF3) is essential for megakaryocyte maturation and generation of functional platelets. Megakaryocyte-specific deletion of Srsf3 in mice led to macrothrombocytopenia characterized by megakaryocyte maturation arrest, dramatically reduced platelet counts, and abnormally large functionally compromised platelets. SRSF3 deficient megakaryocytes failed to reprogram their transcriptome during maturation and to load platelets with RNAs required for normal platelet function. SRSF3 depletion led to nuclear accumulation of megakaryocyte mRNAs, demonstrating that SRSF3 deploys similar RNA regulatory mechanisms in megakaryocytes as in other cell types. Our study further suggests that SRSF3 plays a role in sorting cytoplasmic megakaryocyte RNAs into platelets and demonstrates how SRSF3-mediated RNA processing forms a central part of megakaryocyte gene regulation. Understanding SRSF3 functions in megakaryocytes and platelets provides key insights into normal thrombopoiesis and platelet pathologies as SRSF3 RNA targets in megakaryocytes are associated with platelet diseases.

Characterization of Alternative Splicing Events in Porcine Skeletal Muscles with Different Intramuscular Fat Contents

Meat quality is one of the most important economic traits in pig breeding and production. Intramuscular fat (IMF) is a major factor that improves meat quality. To better understand the alternative splicing (AS) events underlying meat quality, long-read isoform sequencing (Iso-seq) was used to identify differential (D)AS events between the longissimus thoracis (LT) and semitendinosus (ST), which differ in IMF content, together with short-read RNA-seq. Through Iso-seq analysis, we identified a total of 56,789 novel transcripts covering protein-coding genes, lncRNA, and fusion transcripts that were not previously annotated in pigs. We also identified 456,965 AS events, among which 3930 were DAS events, corresponding to 2364 unique genes. Through integrative analysis of Iso-seq and RNA-seq, we identified 1174 differentially expressed genes (DEGs), among which 122 were DAS genes, i.e., DE-DAS genes. There are 12 overlapped pathways between the top 20 DEGs and DE-DAS genes, as revealed by KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, indicating that DE-DAS genes play important roles in the differential phenotype of LT and ST. Further analysis showed that upregulated DE-DAS genes are more important than downregulated ones in IMF deposition. Fatty acid degradation and the PPAR (peroxisome proliferator-activated receptor) signaling pathway were found to be the most important pathways regulating the differential fat deposition of the two muscles. The results update the existing porcine genome annotations and provide data for the in-depth exploration of the mechanisms underlying meat quality and IMF deposition.

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.

Uncovering the impacts of alternative splicing on the proteome with current omics techniques

The high‐throughput sequencing of cellular RNAs has underscored a broad effect of isoform diversification through alternative splicing on the transcriptome. Moreover, the differential production of transcript isoforms from gene loci has been recognized as a critical mechanism in cell differentiation, organismal development, and disease. Yet, the extent of the impact of alternative splicing on protein production and cellular function remains a matter of debate. Multiple experimental and computational approaches have been developed in recent years to address this question. These studies have unveiled how molecular changes at different steps in the RNA processing pathway can lead to differences in protein production and have functional effects. New and emerging experimental technologies open exciting new opportunities to develop new methods to fully establish the connection between messenger RNA expression and protein production and to further investigate how RNA variation impacts the proteome and cell function.

This article is categorized under:

RNA Processing > Splicing Regulation/Alternative Splicing

Translation > Regulation

RNA Evolution and Genomics > Computational Analyses of RNA

Differential fates of introns in gene expression due to global alternative splicing

The discovery of introns over four decades ago revealed a new vision of genes and their interrupted arrangement. Throughout the years, it has appeared that introns play essential roles in the regulation of gene expression. Unique processing of excised introns through the formation of lariats suggests a widespread role for these molecules in the structure and function of cells. In addition to rapid destruction, these lariats may linger on in the nucleus or may even be exported to the cytoplasm, where they remain stable circular RNAs (circRNAs). Alternative splicing (AS) is a source of diversity in mature transcripts harboring retained introns (RI-mRNAs). Such RNAs may contain one or more entire retained intron(s) (RIs), but they may also have intron fragments resulting from sequential excision of smaller subfragments via recursive splicing (RS), which is characteristic of long introns. There are many potential fates of RI-mRNAs, including their downregulation via nuclear and cytoplasmic surveillance systems and the generation of new protein isoforms with potentially different functions. Various reports have linked the presence of such unprocessed transcripts in mammals to important roles in normal development and in disease-related conditions. In certain human neurological-neuromuscular disorders, including myotonic dystrophy type 2 (DM2), frontotemporal dementia/amyotrophic lateral sclerosis (FTD/ALS) and Duchenne muscular dystrophy (DMD), peculiar processing of long introns has been identified and is associated with their pathogenic effects. In this review, we discuss different mechanisms involved in the processing of introns during AS and the functions of these large sections of the genome in our biology.

Dynamic intron retention modulates gene expression in the monocytic differentiation pathway

Monocytes play a crucial role in maintaining homeostasis and mediating a successful innate immune response. They also act as central players in diverse pathological conditions, thus making them an attractive therapeutic target. Within the bone marrow, monocytes arise from a committed precursor termed cMoP (Common Monocyte Progenitor). However, molecular mechanisms that regulate the differentiation of cMoP to various monocytic subsets remain unclear. Herein, we purified murine myeloid precursors for deep poly-A enriched RNA sequencing to understand the role of alternative splicing in the development and differentiation of monocytes under homeostasis. Our analyses revealed intron retention to be the major alternative splicing mechanism involved in the monocyte differentiation cascade, especially in the differentiation of Ly6C^hi monocytes to Ly6C^lo monocytes. Furthermore, we found that the key genes regulated by intron retention in the differentiation of murine Ly6C^hi to Ly6C^lo monocytes were also conserved in humans. Our data highlight the unique role of intron retention in the regulation of the monocytic differentiation pathway.

Using the Zebrafish as a Genetic Model to Study Erythropoiesis

Vertebrates generate mature red blood cells (RBCs) via a highly regulated, multistep process called erythropoiesis. Erythropoiesis involves synthesis of heme and hemoglobin, clearance of the nuclei and other organelles, and remodeling of the plasma membrane, and these processes are exquisitely coordinated by specific regulatory factors including transcriptional factors and signaling molecules. Defects in erythropoiesis can lead to blood disorders such as congenital dyserythropoietic anemias, Diamond–Blackfan anemias, sideroblastic anemias, myelodysplastic syndrome, and porphyria. The molecular mechanisms of erythropoiesis are highly conserved between fish and mammals, and the zebrafish (Danio rerio) has provided a powerful genetic model for studying erythropoiesis. Studies in zebrafish have yielded important insights into RBC development and established a number of models for human blood diseases. Here, we focus on latest discoveries of the molecular processes and mechanisms regulating zebrafish erythropoiesis and summarize newly established zebrafish models of human anemias.

The Expanding Role of Alternative Splicing in Vascular Smooth Muscle Cell Plasticity

Vascular smooth muscle cells (VSMCs) display extraordinary phenotypic plasticity. This allows them to differentiate or dedifferentiate, depending on environmental cues. The ability to ‘switch’ between a quiescent contractile phenotype to a highly proliferative synthetic state renders VSMCs as primary mediators of vascular repair and remodelling. When their plasticity is pathological, it can lead to cardiovascular diseases such as atherosclerosis and restenosis. Coinciding with significant technological and conceptual innovations in RNA biology, there has been a growing focus on the role of alternative splicing in VSMC gene expression regulation. Herein, we review how alternative splicing and its regulatory factors are involved in generating protein diversity and altering gene expression levels in VSMC plasticity. Moreover, we explore how recent advancements in the development of splicing-modulating therapies may be applied to VSMC-related pathologies.

A dynamic intron retention program regulates the expression of several hundred genes during pollen meiosis

Intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during early stages of pollen development. To improve our understanding of the gene regulatory mechanisms that drive developmental processes, we performed a genome-wide study of alternative splicing and isoform switching during five key stages of pollen development in field mustard, Brassica rapa. Surprisingly, for several hundred genes (12.3% of the genes analysed), isoform switching results in stage-specific expression of intron-retaining transcripts at the meiotic stage of pollen development. In such cases, we report temporally regulated switching between expression of a canonical, translatable isoform and an intron-retaining transcript that is predicted to produce a truncated and presumably inactive protein. The results suggest a new pervasive mechanism underlying modulation of protein levels in a plant developmental program. The effect is not based on gene expression induction but on the type of transcript produced. We conclude that intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during meiosis, especially genes related to ribosome biogenesis, mRNA transport and nuclear envelope architecture. We also propose that stage-specific expression of a non-functional isoform of Brassica rapa BrSDG8, a non-redundant member of histone methyltransferase gene family, linked to alternative splicing regulation, may contribute to the intron retention observed.

Single-cell transcriptomics reveal temporal dynamics of critical regulators of germ cell fate during mouse sex determination

Despite the importance of germ cell (GC) differentiation for sexual reproduction, the gene networks underlying their fate remain unclear. Here, we comprehensively characterize the gene expression dynamics during sex determination based on single-cell RNA sequencing of 14 914 XX and XY mouse GCs between embryonic days (E) 9.0 and 16.5. We found that XX and XY GCs diverge transcriptionally as early as E11.5 with upregulation of genes downstream of the bone morphogenic protein (BMP) and nodal/Activin pathways in XY and XX GCs, respectively. We also identified a sex-specific upregulation of genes associated with negative regulation of mRNA processing and an increase in intron retention consistent with a reduction in mRNA splicing in XY testicular GCs by E13.5. Using computational gene regulation network inference analysis, we identified sex-specific, sequential waves of putative key regulator genes during GC differentiation and revealed that the meiotic genes are regulated by positive and negative master modules acting in an antagonistic fashion. Finally, we found that rare adrenal GCs enter meiosis similarly to ovarian GCs but display altered expression of master genes controlling the female and male genetic programs, indicating that the somatic environment is important for GC function. Our data are available on a web platform and provide a molecular roadmap of GC sex determination at single-cell resolution, which will serve as a valuable resource for future studies of gonad development, function, and disease.

Identification of altered cell signaling pathways using proteomic profiling in stable and progressive chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) is characterized by significant biologic and clinical heterogeneity. This study was designed to explore CLL B-cells' proteomic profile in order to identify biologic processes affected at an early stage and during disease evolution as stable or progressive. Purified B cells from 11 untreated CLL patients were tested at two time points by liquid chromatography-tandem mass spectrometry. Patients included in the study evolved to either progressive (n = 6) or stable disease (n = 5). First, at an early stage of the disease (Binet stage A), based on the relative abundance levels of 389 differentially expressed proteins (DEPs), samples were separated into stable and progressive clusters with the main differentiating factor being the RNA splicing pathway. Next, in order to test how the DEPs affect RNA splicing, a RNA-Seq study was conducted showing 4217 differentially spliced genes between the two clusters. Distinct longitudinal evolutions were observed with predominantly proteomic modifications in the stable CLL group and spliced genes in the progressive CLL group. Splicing events were shown to be six times more frequent in the progressive CLL group. The main aberrant biologic processes controlled by DEPs and spliced genes in the progressive group were cytoskeletal organization, Wnt/β-catenin signaling, and mitochondrial and inositol phosphate metabolism with a downstream impact on CLL B-cell survival and migration. This study suggests that proteomic profiles at the early stage of CLL can discriminate progressive from stable disease and that RNA splicing dysregulation underlies CLL evolution, which opens new perspectives in terms of biomarkers and therapy.

Coordinated changes in gene expression kinetics underlie both mouse and human erythroid maturation

BACKGROUND

Single-cell technologies are transforming biomedical research, including the recent demonstration that unspliced pre-mRNA present in single-cell RNA-Seq permits prediction of future expression states. Here we apply this RNA velocity concept to an extended timecourse dataset covering mouse gastrulation and early organogenesis.

RESULTS

Intriguingly, RNA velocity correctly identifies epiblast cells as the starting point, but several trajectory predictions at later stages are inconsistent with both real-time ordering and existing knowledge. The most striking discrepancy concerns red blood cell maturation, with velocity-inferred trajectories opposing the true differentiation path. Investigating the underlying causes reveals a group of genes with a coordinated step-change in transcription, thus violating the assumptions behind current velocity analysis suites, which do not accommodate time-dependent changes in expression dynamics. Using scRNA-Seq analysis of chimeric mouse embryos lacking the major erythroid regulator Gata1, we show that genes with the step-changes in expression dynamics during erythroid differentiation fail to be upregulated in the mutant cells, thus underscoring the coordination of modulating transcription rate along a differentiation trajectory. In addition to the expected block in erythroid maturation, the Gata1-chimera dataset reveals induction of PU.1 and expansion of megakaryocyte progenitors. Finally, we show that erythropoiesis in human fetal liver is similarly characterized by a coordinated step-change in gene expression.

CONCLUSIONS

By identifying a limitation of the current velocity framework coupled with in vivo analysis of mutant cells, we reveal a coordinated step-change in gene expression kinetics during erythropoiesis, with likely implications for many other differentiation processes.

Splicing regulation in hematopoiesis

PURPOSE OF REVIEW

Splicing mutations are among the most recurrent genetic perturbations in hematological malignancies, highlighting an important impact of splicing regulation in hematopoietic development. However, compared to our understanding of splicing factor mutations in hematological malignancies, studies of splicing components and alternative splicing in normal hematopoiesis have been less well investigated. Here, we outline the most recent findings on splicing regulation in normal hematopoiesis and discuss the important questions in the field.

RECENT FINDINGS

Recent studies have highlighted the critical role of splicing regulation in hematopoiesis, including characterization of splicing components in normal hematopoiesis, investigation of transcriptional alterations on splicing, and identification of stage-specific alternative splicing events during hematopoietic development.

SUMMARY

These interesting findings provide insights on hematopoietic regulation at a co-transcriptional level. More high-throughput RNA ribonucleic acid (RNA) sequencing and functional genomic screens are needed to advance our knowledge of critical alternative splicing patterns in shaping hematopoiesis.

HNRNPM controls circRNA biogenesis and splicing fidelity to sustain cancer cell fitness

High spliceosome activity is a dependency for cancer cells, making them more vulnerable to perturbation of the splicing machinery compared to normal cells. To identify splicing factors important for prostate cancer (PCa) fitness, we performed pooled shRNA screens in vitro and in vivo. Our screens identified heterogeneous nuclear ribonucleoprotein M (HNRNPM) as a regulator of PCa cell growth. RNA- and eCLIP-sequencing identified HNRNPM binding to transcripts of key homeostatic genes. HNRNPM binding to its targets prevents aberrant exon inclusion and backsplicing events. In both linear and circular mis-spliced transcripts, HNRNPM preferentially binds to GU-rich elements in long flanking proximal introns. Mimicry of HNRNPM-dependent linear-splicing events using splice-switching-antisense-oligonucleotides was sufficient to inhibit PCa cell growth. This suggests that PCa dependence on HNRNPM is likely a result of mis-splicing of key homeostatic coding and non-coding genes. Our results have further been confirmed in other solid tumors. Taken together, our data reveal a role for HNRNPM in supporting cancer cell fitness. Inhibition of HNRNPM activity is therefore a potential therapeutic strategy in suppressing growth of PCa and other solid tumors.

Integrative functional genomic analysis of intron retention in human and mouse brain with Alzheimer's disease

Intron retention (IR) has been implicated in the pathogenesis of complex diseases such as cancers; its association with Alzheimer's disease (AD) remains unexplored. We performed genome-wide analysis of IR through integrating genetic, transcriptomic, and proteomic data of AD subjects and mouse models from the Accelerating Medicines Partnership-Alzheimer's Disease project. We identified 4535 and 4086 IR events in 2173 human and 1736 mouse genes, respectively. Quantitation of IR enabled the identification of differentially expressed genes that conventional exon-level approaches did not reveal. There were significant correlations of intron expression within innate immune genes, like HMBOX1, with AD in humans. Peptides with a high probability of translation from intron-retained mRNAs were identified using mass spectrometry. Further, we established AD-specific intron expression Quantitative Trait Loci, and identified splicing-related genes that may regulate IR. Our analysis provides a novel resource for the search for new AD biomarkers and pathological mechanisms.

Intron retention as a new pre-symptomatic marker of aging and its recovery to the normal state by a traditional Japanese multi-herbal medicine

Intron retention (IR) is an important regulatory mechanism that affects gene expression and protein functions. Using klotho mice at the pre-symptomatic state, we discovered that retained-introns accumulated in several organs including the liver and that among these retained introns in the liver a subset was recovered to the normal state by a Japanese traditional herbal medicine. This is the first report of IR recovery by a medicine. IR-recovered genes fell into two categories: those involved in liver-specific metabolism and in splicing. Metabolome analysis of the liver showed that the klotho mice were under starvation stress. In addition, our differentially expressed gene analysis showed that liver metabolism was actually recovered by the herbal medicine at the transcriptional level. By analogy with the widespread accumulation of intron-retained pre-mRNAs induced by heat shock stress, we propose a model in which retained-introns in klotho mice were induced by an aging stress and in which this medicine-related IR recovery is indicative of the actual recovery of liver-specific metabolic function to the healthy state. Accumulation of retained-introns was also observed at the pre-symptomatic state of aging in wild-type mice and may be an excellent marker for this state in general.

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

Steps of mRNA maturation are important gene regulatory events that occur in distinct cellular locations. However, transcriptomic analyses often lose information on the subcellular distribution of processed and unprocessed transcripts. We generated extensive RNA-seq data sets to track mRNA maturation across subcellular locations in mouse embryonic stem cells, neuronal progenitor cells, and postmitotic neurons. We find disparate patterns of RNA enrichment between the cytoplasmic, nucleoplasmic, and chromatin fractions, with some genes maintaining more polyadenylated RNA in chromatin than in the cytoplasm. We bioinformatically defined four regulatory groups for intron retention, including complete cotranscriptional splicing, complete intron retention in the cytoplasmic RNA, and two intron groups present in nuclear and chromatin transcripts but fully excised in cytoplasm. We found that introns switch their regulatory group between cell types, including neuronally excised introns repressed by polypyrimidine track binding protein 1 (PTBP1). Transcripts for the neuronal gamma-aminobutyric acid (GABA) B receptor, 1 (Gabbr1) are highly expressed in mESCs but are absent from the cytoplasm. Instead, incompletely spliced Gabbr1 RNA remains sequestered on chromatin, where it is bound by PTBP1, similar to certain long noncoding RNAs. Upon neuronal differentiation, Gabbr1 RNA becomes fully processed and exported for translation. Thus, splicing repression and chromatin anchoring of RNA combine to allow posttranscriptional regulation of Gabbr1 over development. For this and other genes, polyadenylated RNA abundance does not indicate functional gene expression. Our data sets provide a rich resource for analyzing many other aspects of mRNA maturation in subcellular locations and across development.

Nuclear mechanisms of gene expression control: pre-mRNA splicing as a life or death decision

Thousands of genes produce polyadenylated mRNAs that still contain one or more introns. These transcripts are known as retained intron RNAs (RI-RNAs). In the past 10 years, RI-RNAs have been linked to post-transcriptional alternative splicing in a variety of developmental contexts, but they can also be dead-end products fated for RNA decay. Here we discuss the role of intron retention in shaping gene expression programs, as well as recent evidence suggesting that the biogenesis and fate of RI-RNAs is regulated by nuclear organization. We discuss the possibility that proximity of RNA to nuclear speckles - biomolecular condensates that are highly enriched in splicing factors and other RNA binding proteins - is associated with choices ranging from efficient co-transcriptional splicing, export and stability to regulated post-transcriptional splicing and possible vulnerability to decay.

A splicing factor switch controls hematopoietic lineage specification of pluripotent stem cells

Alternative splicing (AS) leads to transcriptome diversity in eukaryotic cells and is one of the key regulators driving cellular differentiation. Although AS is of crucial importance for normal hematopoiesis and hematopoietic malignancies, its role in early hematopoietic development is still largely unknown. Here, by using high-throughput transcriptomic analyses, we show that pervasive and dynamic AS takes place during hematopoietic development of human pluripotent stem cells (hPSCs). We identify a splicing factor switch that occurs during the differentiation of mesodermal cells to endothelial progenitor cells (EPCs). Perturbation of this switch selectively impairs the emergence of EPCs and hemogenic endothelial progenitor cells (HEPs). Mechanistically, an EPC-induced alternative spliced isoform of NUMB dictates EPC specification by controlling NOTCH signaling. Furthermore, we demonstrate that the splicing factor SRSF2 regulates splicing of the EPC-induced NUMB isoform, and the SRSF2-NUMB-NOTCH splicing axis regulates EPC generation. The identification of this splicing factor switch provides a new molecular mechanism to control cell fate and lineage specification.

Lithium increases mitochondrial respiration in iPSC-derived neural precursor cells from lithium responders

Lithium (Li), valproate (VPA) and lamotrigine (LTG) are commonly used to treat bipolar disorder (BD). While their clinical efficacy is well established, the mechanisms of action at the molecular level are still incompletely understood. Here we investigated the molecular effects of Li, LTG and VPA treatment in induced pluripotent stem cell (iPSC)-derived neural precursor cells (NPCs) generated from 3 healthy controls (CTRL), 3 affective disorder Li responsive patients (Li-R) and 3 Li non-treated patients (Li-N) after 6 h and 1 week of exposure. Differential expression (DE) analysis after 6 h of treatment revealed a transcriptional signature that was associated with all three drugs and most significantly enriched for ribosome and oxidative phosphorylation (OXPHOS) pathways. In addition to the shared DE genes, we found that Li exposure was associated with 554 genes uniquely regulated in Li-R NPCs and enriched for spliceosome, OXPHOS and thermogenesis pathways. In-depth analysis of the treatment-associated transcripts uncovered a significant decrease in intron retention rate, suggesting that the beneficial influence of these drugs might partly be related to splicing. We examined the mitochondrial respiratory function of the NPCs by exploring the drugs' effects on oxygen consumption rate (OCR) and glycolytic rate (ECAR). Li improved OCR levels only in Li-R NPCs by enhancing maximal respiration and reserve capacity, while VPA enhanced maximal respiration and reserve capacity in Li-N NPCs. Overall, our findings further support the involvement of mitochondrial functions in the molecular mechanisms of mood stabilizers and suggest novel mechanisms related to the spliceosome, which warrant further investigation.