The changing paradigm of intron retention: regulation, ramifications and recipes

NSrp70 is a lymphocyte-essential splicing factor that controls thymocyte development

Alternative pre-mRNA splicing is a critical step to generate multiple transcripts, thereby dramatically enlarging the proteomic diversity. Thus, a common feature of most alternative splicing factor knockout models is lethality. However, little is known about lineage-specific alternative splicing regulators in a physiological setting. Here, we report that NSrp70 is selectively expressed in developing thymocytes, highest at the double-positive (DP) stage. Global splicing and transcriptional profiling revealed that NSrp70 regulates the cell cycle and survival of thymocytes by controlling the alternative processing of various RNA splicing factors, including the oncogenic splicing factor SRSF1. A conditional-knockout of Nsrp1 (NSrp70-cKO) using CD4Cre developed severe defects in T cell maturation to single-positive thymocytes, due to insufficient T cell receptor (TCR) signaling and uncontrolled cell growth and death. Mice displayed severe peripheral lymphopenia and could not optimally control tumor growth. This study establishes a model to address the function of lymphoid-lineage-specific alternative splicing factor NSrp70 in a thymic T cell developmental pathway.

Transcriptome analysis reveals Vernalization is independent of cold acclimation in Arabidopsis

Background

Through vernalization, plants achieve flowering competence by sensing prolonged cold exposure (constant exposure approximately 2-5 °C). During this process, plants initiate defense responses to endure cold conditions. Here, we conducted transcriptome analysis of Arabidopsis plants subjected to prolonged cold exposure (6 weeks) to explore the physiological dynamics of vernalization and uncover the relationship between vernalization and cold stress.

Results

Time-lag initiation of the two pathways and weighted gene co-expression network analysis (WGCNA) revealed that vernalization is independent of cold acclimation. Moreover, WGCNA revealed three major networks involving ethylene and jasmonic acid response, cold acclimation, and chromatin modification in response to prolonged cold exposure. Finally, throughout vernalization, the cold stress response is regulated via an alternative splicing-mediated mechanism.

Conclusion

These findings illustrate a comprehensive picture of cold stress- and vernalization-mediated global changes in Arabidopsis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07763-3.

Nuclear compartmentalization of TERT mRNA and TUG1 lncRNA is driven by intron retention

The spatial partitioning of the transcriptome in the cell is an important form of gene-expression regulation. Here, we address how intron retention influences the spatio-temporal dynamics of transcripts from two clinically relevant genes: TERT (Telomerase Reverse Transcriptase) pre-mRNA and TUG1 (Taurine-Upregulated Gene 1) lncRNA. Single molecule RNA FISH reveals that nuclear TERT transcripts uniformly and robustly retain specific introns. Our data suggest that the splicing of TERT retained introns occurs during mitosis. In contrast, TUG1 has a bimodal distribution of fully spliced cytoplasmic and intron-retained nuclear transcripts. We further test the functionality of intron-retention events using RNA-targeting thiomorpholino antisense oligonucleotides to block intron excision. We show that intron retention is the driving force for the nuclear compartmentalization of these RNAs. For both RNAs, altering this splicing-driven subcellular distribution has significant effects on cell viability. Together, these findings show that stable retention of specific introns can orchestrate spatial compartmentalization of these RNAs within the cell. This process reveals that modulating RNA localization via targeted intron retention can be utilized for RNA-based therapies.

Gene expression and alternative splicing dynamics are perturbed in female head transcriptomes following heterospecific copulation

BACKGROUND

Despite the growing interest in the female side of copulatory interactions, the roles played by differential expression and alternative splicing mechanisms of pre-RNA on tissues outside of the reproductive tract have remained largely unknown. Here we addressed these questions in the context of con- vs heterospecific matings between Drosophila mojavensis and its sister species, D. arizonae. We analyzed transcriptional responses in female heads using an integrated investigation of genome-wide patterns of gene expression, including differential expression (DE), alternative splicing (AS) and intron retention (IR).

RESULTS

Our results indicated that early transcriptional responses were largely congruent between con- and heterospecific matings but are substantially perturbed over time. Conspecific matings induced functional pathways related to amino acid balance previously associated with the brain's physiology and female postmating behavior. Heterospecific matings often failed to activate regulation of some of these genes and induced expression of additional genes when compared with those of conspecifically-mated females. These mechanisms showed functional specializations with DE genes mostly linked to pathways of proteolysis and nutrient homeostasis, while AS genes were more related to photoreception and muscle assembly pathways. IR seems to play a more general role in DE regulation during the female postmating response.

CONCLUSIONS

We provide evidence showing that AS genes substantially perturbed by heterospecific matings in female heads evolve at slower evolutionary rates than the genome background. However, DE genes evolve at evolutionary rates similar, or even higher, than those of male reproductive genes, which highlights their potential role in sexual selection and the evolution of reproductive barriers.

The RNA binding protein FgRbp1 regulates specific pre-mRNA splicing via interacting with U2AF23 in Fusarium

Precursor messenger RNA (pre-mRNA) splicing is an essential and tightly regulated process in eukaryotic cells; however, the regulatory mechanisms for the splicing are not well understood. Here, we characterize a RNA binding protein named FgRbp1 in Fusarium graminearum, a fungal pathogen of cereal crops worldwide. Deletion of FgRbp1 leads to reduced splicing efficiency in 47% of the F. graminearum intron-containing gene transcripts that are involved in various cellular processes including vegetative growth, development, and virulence. The human ortholog RBM42 is able to fully rescue the growth defects of ΔFgRbp1. FgRbp1 binds to the motif CAAGR in its target mRNAs, and interacts with the splicing factor FgU2AF23, a highly conserved protein involved in 3' splice site recognition, leading to enhanced recruitment of FgU2AF23 to the target mRNAs. This study demonstrates that FgRbp1 is a splicing regulator and regulates the pre-mRNA splicing in a sequence-dependent manner in F. graminearum.

Regional identity of human neural stem cells determines oncogenic responses to histone H3.3 mutants

Point mutations within the histone H3.3 are frequent in aggressive childhood brain tumors known as pediatric high-grade gliomas (pHGGs). Intriguingly, distinct mutations arise in discrete anatomical regions: H3.3-G34R within the forebrain and H3.3-K27M preferentially within the hindbrain. The reasons for this contrasting etiology are unknown. By engineering human fetal neural stem cell cultures from distinct brain regions, we demonstrate here that cell-intrinsic regional identity provides differential responsiveness to each mutant that mirrors the origins of pHGGs. Focusing on H3.3-G34R, we find that the oncohistone supports proliferation of forebrain cells while inducing a cytostatic response in the hindbrain. Mechanistically, H3.3-G34R does not impose widespread transcriptional or epigenetic changes but instead impairs recruitment of ZMYND11, a transcriptional repressor of highly expressed genes. We therefore propose that H3.3-G34R promotes tumorigenesis by focally stabilizing the expression of key progenitor genes, thereby locking initiating forebrain cells into their pre-existing immature state.

Meta-Analyses of Splicing and Expression Quantitative Trait Loci Identified Susceptibility Genes of Glioma

BACKGROUND

The functions of most glioma risk alleles are unknown. Very few studies had evaluated expression quantitative trait loci (eQTL), and insights of susceptibility genes were limited due to scarcity of available brain tissues. Moreover, no prior study had examined the effect of glioma risk alleles on alternative RNA splicing.

OBJECTIVE

This study explored splicing quantitative trait loci (sQTL) as molecular QTL and improved the power of QTL mapping through meta-analyses of both cis eQTL and sQTL.

METHODS

We first evaluated eQTLs and sQTLs of the CommonMind Consortium (CMC) and Genotype-Tissue Expression Project (GTEx) using genotyping, or whole-genome sequencing and RNA-seq data. Alternative splicing events were characterized using an annotation-free method that detected intron excision events. Then, we conducted meta-analyses by pooling the eQTL and sQTL results of CMC and GTEx using the inverse variance-weighted model. Afterward, we integrated QTL meta-analysis results (Q < 0.05) with the Glioma International Case Control Study (GICC) GWAS meta-analysis (case:12,496, control:18,190), using a summary statistics-based mendelian randomization (SMR) method.

RESULTS

Between CMC and GTEx, we combined the QTL data of 354 unique individuals of European ancestry. SMR analyses revealed 15 eQTLs in 11 loci and 32 sQTLs in 9 loci relevant to glioma risk. Two loci only harbored sQTLs (1q44 and 16p13.3). In seven loci, both eQTL and sQTL coexisted (2q33.3, 7p11.2, 11q23.3 15q24.2, 16p12.1, 20q13.33, and 22q13.1), but the target genes were different for five of these seven loci. Three eQTL loci (9p21.3, 20q13.33, and 22q13.1) and 4 sQTL loci (11q23.3, 16p13.3, 16q12.1, and 20q13.33) harbored multiple target genes. Eight target genes of sQTLs (C2orf80, SEC61G, TMEM25, PHLDB1, RP11-161M6.2, HEATR3, RTEL1-TNFRSF6B, and LIME1) had multiple alternatively spliced transcripts.

CONCLUSION

Our study revealed that the regulation of transcriptome by glioma risk alleles is complex, with the potential for eQTL and sQTL jointly affecting gliomagenesis in risk loci. QTLs of many loci involved multiple target genes, some of which were specific to alternative splicing. Therefore, quantitative trait loci that evaluate only total gene expression will miss many important target genes.

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.

Novel mutations in the WFS1 gene are associated with Wolfram syndrome and systemic inflammation

Mutations in the WFS1 gene, encoding wolframin (WFS1), cause endoplasmic reticulum (ER) stress and are associated with a rare autosomal-recessive disorder known as Wolfram syndrome (WS). WS is clinically characterized by childhood-onset diabetes mellitus, optic atrophy, deafness, diabetes insipidus and neurological signs. We identified two novel WFS1 mutations in a patient with WS, namely, c.316-1G > A (in intron 3) and c.757A > T (in exon 7). Both mutations, located in the N-terminal region of the protein, were predicted to generate a truncated and inactive form of WFS1. We found that although the WFS1 protein was not expressed in peripheral blood mononuclear cells (PBMCs) of the proband, no constitutive ER stress activation could be detected in those cells. In contrast, WS proband’s PBMCs produced very high levels of proinflammatory cytokines (i.e. TNF-α, IL-1β, and IL-6) in the absence of any stimulus. WFS1 silencing in PBMCs from control subjects by means of small RNA interference also induced a pronounced proinflammatory cytokine profile. The same cytokines were also significantly higher in sera from the WS patient as compared to matched healthy controls. Moreover, the chronic inflammatory state was associated with a dominance of proinflammatory T helper 17 (Th17)-type cells over regulatory T (Treg) lymphocytes in the WS PBMCs. The identification of a state of systemic chronic inflammation associated with WFS1 deficiency may pave the way to innovative and personalized therapeutic interventions in WS.

CRISPR/Cas9 Genome Editing of the Human Topoisomerase IIα Intron 19 5' Splice Site Circumvents Etoposide Resistance in Human Leukemia K562 Cells

An essential function of DNA topoisomerase IIα (TOP2α; 170 kDa, TOP2α/170) is to resolve DNA topologic entanglements during chromosome disjunction by introducing transient DNA double-stranded breaks. TOP2α/170 is an important target for DNA damage-stabilizing anticancer drugs, whose clinical efficacy is compromised by drug resistance often associated with decreased TOP2α/170 expression. We recently demonstrated that an etoposide-resistant K562 clonal subline, K/VP.5, with reduced levels of TOP2α/170, expresses high levels of a novel C-terminal truncated TOP2α isoform (90 kDa, TOP2α/90). TOP2α/90, the translation product of a TOP2α mRNA that retains a processed intron 19 (I19), heterodimerizes with TOP2α/170 and is a resistance determinant through a dominant-negative effect on drug activity. We hypothesized that genome editing to enhance I19 removal would provide a tractable strategy to circumvent acquired TOP2α-mediated drug resistance. To enhance I19 removal in K/VP.5 cells, CRISPR/Cas9 was used to make changes (GAG//GTAA AC →GAG//GTAA GT ) in the TOP2α gene's suboptimal exon 19/intron 19 5' splice site (E19/I19 5' SS). Gene-edited clones were identified by quantitative polymerase chain reaction and verified by sequencing. Characterization of a clone with all TOP2α alleles edited revealed improved I19 removal, decreased TOP2α/90 mRNA/protein, and increased TOP2α/170 mRNA/protein. Sensitivity to etoposide-induced DNA damage (γH2AX, Comet assays) and growth inhibition was restored to levels comparable to those in parental K562 cells. Together, the results indicate that our gene-editing strategy for optimizing the TOP2α E19/I19 5' SS in K/VP.5 cells circumvents resistance to etoposide and other TOP2α-targeted drugs. SIGNIFICANCE STATEMENT: Results presented here indicate that CRISPR/Cas9 gene editing of a suboptimal exon 19/intron 19 5' splice site in the DNA topoisomerase IIα (TOP2α) gene results in circumvention of acquired drug resistance to etoposide and other TOP2α-targeted drugs in a clonal K562 cell line by enhancing removal of intron 19 and thereby decreasing formation of a truncated TOP2α 90 kDa isoform and increasing expression of full-length TOP2α 170 kDa in these resistant cells. Results demonstrate the importance of RNA processing in acquired drug resistance to TOP2α-targeted drugs.

Molecular characterization and functional analysis of Eimeria tenella citrate synthase

Chicken coccidiosis, caused by an obligate intracellular protozoan parasite of the genus Eimeria, is a major parasitic disease in the intensively reared poultry industry. Due to the widespread use of anticoccidial drugs, resistance has become an inevitable problem. In our previous study, Eimeria tenella citrate synthase (EtCS) was found to be up-expressed in two drug-resistant strains (diclazuril-resistant and maduramycin-resistant strains) compared to drug-sensitive strain by RNA sequence. In this study, we cloned and expressed EtCS and obtain its polyclonal antibodies. Quantitative real-time polymerase chain (qPCR) reactions and Western blots were used to analyze the transcription and translation levels of EtCS in sensitive and three drug-resistant strains. Compared with the sensitive strain, the transcription of EtCS was both significantly upregulated in diclazuril-resistant and maduramycin-resistant strains, but was not significantly different in salinomycin-resistant strain. No significant difference was seen in translation level in the three drug-resistant strains. Indirect immunofluorescence indicated that EtCS was mainly located in the cytoplasm of sporozoites except for posterior refractile bodies and in the cytoplasm and surface of merozoites. Anti-rEtCS antibody has inhibitory effects on E. tenella sporozoite invasion of DF-1 cells and the inhibition rate is more than 83%. Binding of the protein to chicken macrophage (HD11) cells was confirmed by immunofluorescence assays. When macrophages were treated with rEtCS, secretion of nitric oxide and cell proliferation of the macrophages were substantially reduced. These results showed that EtCS may be related to host cell invasion of E. tenella and involve in the development of E.tenella resistance to some drugs.

Possible role played by the SINE2 element in gene regulation, as demonstrated by differential processing and polyadenylation in avirulent strains of E. histolytica

Entamoeba histolytica represents a useful model in parasitic organisms due to its complex genomic organization and survival mechanisms. To counteract pathogenic organisms, it is necessary to characterize their molecular biology to design new strategies to combat them. In this report, we investigated a less-known genetic element, short interspersed nuclear element 2 (SINE2), that is present in this ameba and is highly transcribed and polyadenylated. In this study, we show that in two different nonvirulent strains of E. histolytica, SINE2 is differentially processed into two transcript fragments, that is, a full-length 560-nt fragment and a shorter 393-nt fragment bearing an approximately 18-nt polyadenylation tail. Sequence analysis of the SINE2 transcript showed that a Musashi-like protein may bind to it. Also, two putative Musashi-like sequences were identified on the transcript. Semiquantitative expression analysis of the two Musashi-like proteins identified in the E. histolytica genome (XP_648918 and XP_649094) showed that XP_64094 is overexpressed in the nonvirulent strains tested. The information available in the literature and the results presented in this report indicate that SINE2 may affect other genes, as observed with the epigenetic silencing of the G3 strain, by an antisense mechanism or via RNA-protein interactions that may ultimately be involved in the phenotype of nonvirulent strains of E. histolytica.

A Novel Splice-Site Mutation in the ELN Gene Suggests an Alternative Mechanism for Vascular Elastinopathies

The ELN gene encodes elastin, a fundamental protein of the extracellular matrix that confers elasticity to different tissues including blood vessels. The formation of elastin fibers is a complex process involving monomer coacervation and subsequent crosslinking. Mutations in exons 1-29 of the ELN gene have been linked to supravalvular aortic stenosis (SVAS) whereas mutations in exons 30-33 are associated with autosomal dominant cutis laxa (ADCL). This striking segregation has led to the hypothesis that distinct molecular mechanisms underlie both diseases. SVAS is believed to arise through haploinsufficiency while ADCL is hypothesized to be caused by a dominant negative effect. Here, we describe a patient with SVAS harboring a novel splice-site mutation in the last exon of ELN. The location of this mutation is not consistent with current knowledge of SVAS, since all mutations reported in the C-terminus have been found in ADCL patients, and a thorough evaluation did not reveal significant skin involvement in this case. RT-PCR analysis of skin tissue showed that C-terminal mutations in the region can lead to the production of aberrant transcripts through intron retention and activation of cryptic splice sites and suggest that disruption of the very last exon can lead to functional haploinsufficiency potentially related to SVAS.

Intron retention is a robust marker of intertumoral heterogeneity in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer with a 5-year survival rate of <8%. Unsupervised clustering of 76 PDAC patients based on intron retention (IR) events resulted in two clusters of tumors (IR-1 and IR-2). While gene expression-based clusters are not predictive of patient outcome in this cohort, the clusters we developed based on intron retention were associated with differences in progression-free interval. IR levels are lower and clinical outcome is worse in IR-1 compared with IR-2. Oncogenes were significantly enriched in the set of 262 differentially retained introns between the two IR clusters. Higher IR levels in IR-2 correlate with higher gene expression, consistent with detention of intron-containing transcripts in the nucleus in IR-2. Out of 258 genes encoding RNA-binding proteins (RBP) that were differentially expressed between IR-1 and IR-2, the motifs for seven RBPs were significantly enriched in the 262-intron set, and the expression of 25 RBPs were highly correlated with retention levels of 139 introns. Network analysis suggested that retention of introns in IR-2 could result from disruption of an RBP protein-protein interaction network previously linked to efficient intron removal. Finally, IR-based clusters developed for the majority of the 20 cancer types surveyed had two clusters with asymmetrical distributions of IR events like PDAC, with one cluster containing mostly intron loss events. Taken together, our findings suggest IR may be an important biomarker for subclassifying tumors.

Ctcf haploinsufficiency mediates intron retention in a tissue-specific manner

CTCF is a master regulator of gene transcription and chromatin organisation with occupancy at thousands of DNA target sites genome-wide. While CTCF is essential for cell survival, CTCF haploinsufficiency is associated with tumour development and hypermethylation. Increasing evidence demonstrates CTCF as a key player in several mechanisms regulating alternative splicing (AS), however, the genome-wide impact of Ctcf dosage on AS has not been investigated.

We examined the effect of Ctcf haploinsufficiency on gene expression and AS in five tissues from Ctcf hemizygous (Ctcf^+/-) mice. Reduced Ctcf levels caused distinct tissue-specific differences in gene expression and AS in all tissues. An increase in intron retention (IR) was observed in Ctcf^+/- liver and kidney. In liver, this specifically impacted genes associated with cytoskeletal organisation, splicing and metabolism. Strikingly, most differentially retained introns were short, with a high GC content and enriched in Ctcf binding sites in their proximal upstream genomic region. This study provides new insights into the effects of CTCF haploinsufficiency on organ transcriptomes and the role of CTCF in AS regulation.

Intron Retention as a Mode for RNA-Seq Data Analysis

Intron retention (IR) is an alternative splicing mode whereby introns, rather than being spliced out as usual, are retained in mature mRNAs. It was previously considered a consequence of mis-splicing and received very limited attention. Only recently has IR become of interest for transcriptomic data analysis owing to its recognized roles in gene expression regulation and associations with complex diseases. In this article, we first review the function of IR in regulating gene expression in a number of biological processes, such as neuron differentiation and activation of CD4⁺ T cells. Next, we briefly review its association with diseases, such as Alzheimer's disease and cancers. Then, we describe state-of-the-art methods for IR detection, including RNA-seq analysis tools IRFinder and iREAD, highlighting their underlying principles and discussing their advantages and limitations. Finally, we discuss the challenges for IR detection and potential ways in which IR detection methods could be improved.

Direct RNA Sequencing for the Study of Synthesis, Processing, and Degradation of Modified Transcripts

It has been known for a few decades that transcripts can be marked by dozens of different modifications. Yet, we are just at the beginning of charting these marks and understanding their functional impact. High-quality methods were developed for the profiling of some of these marks, and approaches to finely study their impact on specific phases of the RNA life-cycle are available, including RNA metabolic labeling. Thanks to these improvements, the most abundant marks, including N6-methyladenosine, are emerging as important determinants of the fate of marked RNAs. However, we still lack approaches to directly study how the set of marks for a given RNA molecule shape its fate. In this perspective, we first review current leading approaches in the field. Then, we propose an experimental and computational setup, based on direct RNA sequencing and mathematical modeling, to decipher the functional consequences of RNA modifications on the fate of individual RNA molecules and isoforms.

Systematic Analysis of the DNA Methylase and Demethylase Gene Families in Rapeseed (Brassica napus L.) and Their Expression Variations After Salt and Heat stresses

DNA methylation is a process through which methyl groups are added to the DNA molecule, thereby modifying the activity of a DNA segment without changing the sequence. Increasing evidence has shown that DNA methylation is involved in various aspects of plant growth and development via a number of key processes including genomic imprinting and repression of transposable elements. DNA methylase and demethylase are two crucial enzymes that play significant roles in dynamically maintaining genome DNA methylation status in plants. In this work, 22 DNA methylase genes and six DNA demethylase genes were identified in rapeseed (Brassica napus L.) genome. These DNA methylase and DNA demethylase genes can be classified into four (BnaCMTs, BnaMET1s, BnaDRMs and BnaDNMT2s) and three (BnaDMEs, BnaDML3s and BnaROS1s) subfamilies, respectively. Further analysis of gene structure and conserved domains showed that each sub-class is highly conserved between rapeseed and Arabidopsis. Expression analysis conducted by RNA-seq as well as qRT-PCR suggested that these DNA methylation/demethylation-related genes may be involved in the heat/salt stress responses in rapeseed. Taken together, our findings may provide valuable information for future functional characterization of these two types of epigenetic regulatory enzymes in polyploid species such as rapeseed, as well as for analyzing their evolutionary relationships within the plant kingdom.