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

Regulation potential of transcribed simple repeated sequences in developing neurons

Simple repeated sequences (SRSs), defined as tandem iterations of microsatellite- to satellite-sized DNA units, occupy a substantial part of the human genome. Some of these elements are known to be transcribed in the context of repeat expansion disorders. Mounting evidence suggests that the transcription of SRSs may also contribute to normal cellular functions. Here, we used genome-wide bioinformatics approaches to systematically examine SRS transcriptional activity in cells undergoing neuronal differentiation. We identified thousands of long noncoding RNAs containing >200-nucleotide-long SRSs (SRS-lncRNAs), with hundreds of these transcripts significantly upregulated in the neural lineage. We show that SRS-lncRNAs often originate from telomere-proximal regions and that they have a strong potential to form multivalent contacts with a wide range of RNA-binding proteins. Our analyses also uncovered a cluster of neurally upregulated SRS-lncRNAs encoded in a centromere-proximal part of chromosome 9, which underwent an evolutionarily recent segmental duplication. Using a newly established in vitro system for rapid neuronal differentiation of induced pluripotent stem cells, we demonstrate that at least some of the bioinformatically predicted SRS-lncRNAs, including those encoded in the segmentally duplicated part of chromosome 9, indeed increase their expression in developing neurons to readily detectable levels. These and other lines of evidence suggest that many SRSs may be expressed in a cell type and developmental stage-specific manner, providing a valuable resource for further studies focused on the functional consequences of SRS-lncRNAs in the normal development of the human brain, as well as in the context of neurodevelopmental disorders.

Integrating gene expression and splicing dynamics across dose-response oxidative modulators

Toxicological risk assessment increasingly utilizes transcriptomics to derive point of departure (POD) and modes of action (MOA) for chemicals. One essential biological process that allows a single gene to generate several different RNA isoforms is called alternative splicing. To comprehensively assess the role of splicing dysregulation in toxicological evaluation and elucidate its potential as a complementary endpoint, we performed RNA-seq on A549 cells treated with five oxidative stress modulators across a wide dose range. Differential gene expression (DGE) showed limited pathway enrichment except at high concentrations. However, alternative splicing analysis revealed variable intron retention events affecting diverse pathways for all chemicals in the absence of significant expression changes. For instance, diazinon elicited negligible gene expression changes but progressive increase in the number of intron retention events, suggesting splicing alterations precede expression responses. Benchmark dose modeling of intron retention data highlighted relevant pathways overlooked by expression analysis. Systematic integration of splicing datasets should be a useful addition to the toxicogenomic toolkit. Combining both modalities paint a more complete picture of transcriptomic dose-responses. Overall, evaluating intron retention dynamics afforded by toxicogenomics may provide biomarkers that can enhance chemical risk assessment and regulatory decision making. This work highlights splicing-aware toxicogenomics as a possible additional tool for examining cellular responses.

Systematic identification and characterization of exon-intron circRNAs

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

Edodin: A New Type of Toxin from Shiitake Mushroom (Lentinula edodes) That Inactivates Mammalian Ribosomes

Ribosome-inactivating proteins (RIPs) are a group of proteins with rRNA N-glycosylase activity that irreversibly inhibit protein synthesis and consequently cause cell death. Recently, an RIP called ledodin has been found in shiitake; it is cytotoxic, strongly inhibits protein synthesis, and shows rRNA N-glycosylase activity. In this work, we isolated and characterized a 50 kDa cytotoxic protein from shiitake that we named edodin. Edodin inhibits protein synthesis in a mammalian cell-free system, but not in insect-, yeast-, and bacteria-derived systems. It exhibits rRNA N-glycosylase and DNA-nicking activities, which relate it to plant RIPs. It was also shown to be toxic to HeLa and COLO 320 cells. Its structure is not related to other RIPs found in plants, bacteria, or fungi, but, instead, it presents the characteristic structure of the fold type I of pyridoxal phosphate-dependent enzymes. Homologous sequences have been found in other fungi of the class Agaricomycetes; thus, edodin could be a new type of toxin present in many fungi, some of them edible, which makes them of great interest in health, both for their involvement in food safety and for their potential biomedical and biotechnological applications.

LINC00460-FUS-MYC feedback loop drives breast cancer metastasis and doxorubicin resistance

Therapeutic resistance and metastasis largely contribute to mortality from breast cancer and therefore understanding the underlying mechanisms of such remains an urgent challenge. By cross-analysis of TCGA and GEO databases, LINC00460 was identified as an oncogenic long non-coding RNA, highly expressed in Doxorubicin resistant breast cancer. LINC00460 was further demonstrated to promote stem cell-like and epithelial-mesenchymal transition (EMT) characteristics in breast cancer cells. LINC00460 interacts with FUS protein with consequent enhanced stabilization, which further promotes MYC mRNA maturation. LINC00460 expression was transcriptionally enhanced by c-MYC protein, forming a positive feedback loop to promote metastasis and Doxorubicin resistance. LINC00460 depletion in Doxorubicin-resistant breast cancer cells restored sensitivity to Doxorubicin and increased the efficacy of c-MYC inhibitor therapy. Collectively, these findings implicate LINC00460 as a promising prognostic biomarker and potential therapeutic target to overcome Doxorubicin resistance in breast cancer.

TREX tetramer disruption alters RNA processing necessary for corticogenesis in THOC6 Intellectual Disability Syndrome

THOC6 variants are the genetic basis of autosomal recessive THOC6 Intellectual Disability Syndrome (TIDS). THOC6 is critical for mammalian Transcription Export complex (TREX) tetramer formation, which is composed of four six-subunit THO monomers. The TREX tetramer facilitates mammalian RNA processing, in addition to the nuclear mRNA export functions of the TREX dimer conserved through yeast. Human and mouse TIDS model systems revealed novel THOC6-dependent, species-specific TREX tetramer functions. Germline biallelic Thoc6 loss-of-function (LOF) variants result in mouse embryonic lethality. Biallelic THOC6 LOF variants reduce the binding affinity of ALYREF to THOC5 without affecting the protein expression of TREX members, implicating impaired TREX tetramer formation. Defects in RNA nuclear export functions were not detected in biallelic THOC6 LOF human neural cells. Instead, mis-splicing was detected in human and mouse neural tissue, revealing novel THOC6-mediated TREX coordination of mRNA processing. We demonstrate that THOC6 is required for key signaling pathways known to regulate the transition from proliferative to neurogenic divisions during human corticogenesis. Together, these findings implicate altered RNA processing in the developmental biology of TIDS neuropathology.

Alternative splicing and related RNA binding proteins in human health and disease

Alternative splicing (AS) serves as a pivotal mechanism in transcriptional regulation, engendering transcript diversity, and modifications in protein structure and functionality. Across varying tissues, developmental stages, or under specific conditions, AS gives rise to distinct splice isoforms. This implies that these isoforms possess unique temporal and spatial roles, thereby associating AS with standard biological activities and diseases. Among these, AS-related RNA-binding proteins (RBPs) play an instrumental role in regulating alternative splicing events. Under physiological conditions, the diversity of proteins mediated by AS influences the structure, function, interaction, and localization of proteins, thereby participating in the differentiation and development of an array of tissues and organs. Under pathological conditions, alterations in AS are linked with various diseases, particularly cancer. These changes can lead to modifications in gene splicing patterns, culminating in changes or loss of protein functionality. For instance, in cancer, abnormalities in AS and RBPs may result in aberrant expression of cancer-associated genes, thereby promoting the onset and progression of tumors. AS and RBPs are also associated with numerous neurodegenerative diseases and autoimmune diseases. Consequently, the study of AS across different tissues holds significant value. This review provides a detailed account of the recent advancements in the study of alternative splicing and AS-related RNA-binding proteins in tissue development and diseases, which aids in deepening the understanding of gene expression complexity and offers new insights and methodologies for precision medicine.

Genome-wide identification and expression profiles analysis of the authentic response regulator gene family in licorice

Introduction

As one of the traditional Chinese medicinal herbs that were most generally used, licorice attracts lots of interest due to its therapeutic potential. Authentic response regulators (ARRs) are key factors in cytokinin signal transduction and crucial for plant growth and stress response processes. Nevertheless, the characteristics and functions of the licorice ARR genes are still unknown.

Results

In present study, a systematic genome-wide identification and expression analysis of the licorice ARR gene family were conducted and 51 ARR members were identified. Collinearity analysis revealed the significant roles of segmental duplications in the expansion of licorice ARR genes. The cis-acting elements associated with development, stress and phytohormone responses were identified, implying their pivotal roles in diverse regulatory processes. RNA-seq and qRT-PCR results suggested that A-type, but not B-type ARRs were induced by zeatin. Additionally, ARRs participated in diverse abiotic stresses and phytohormones responses. Yeast one-hybrid assay demonstrated that GuARR1, GuARR2, GuARR11, GuARR12, GuARR10-1, GuARR10-2 and GuARR14 were able to bind to the promoter of GuARR8-3, and GuARR1, GuARR12 bound to the GuARR8-1 promoter. GuARR1, GuARR2, GuARR11 and GuARR10-2 bound to the GuARR6-2 promoter as well as GuARR12 and GuARR10-2 bound to the GuARR6-1 promoter.

Discussion

Collectively, these findings provide a basis for future ARR genes function investigations, shedding light on the potential medicinal properties and agricultural applications of licorice.

Extensive Alternative Splicing Patterns in Systemic Lupus Erythematosus Highlight Sexual Differences

Substantial evidence highlights divergences in immune responses between men and women. Women are more susceptible to autoimmunity, whereas men suffer from the more severe presentation of autoimmune disorders. The molecular mechanism of this sexual dimorphism remains elusive. Herein, we conducted a comprehensive analysis of sex differences in whole-blood gene expression focusing on alternative splicing (AS) events in systemic lupus erythematosus (SLE), which is a prototype sex-biased disease. This study included 79 SLE patients with active disease and 58 matched healthy controls who underwent whole-blood RNA sequencing. Sex differences in splicing events were widespread, existent in both SLE and a healthy state. However, we observed distinct gene sets and molecular pathways targeted by sex-dependent AS in SLE patients as compared to healthy subjects, as well as a notable sex dissimilarity in intron retention events. Sexually differential spliced genes specific to SLE patients were enriched for dynamic cellular processes including chromatin remodeling, stress and inflammatory responses. Remarkably, the extent of sexual differences in AS in the SLE patients and healthy individuals exceeded those in gene expression. Overall, this study reveals an unprecedent variation in sex-dependent splicing events in SLE and the healthy state, with potential implications for understanding the molecular basis of sexual dimorphism in autoimmunity.

Unleashing the potential of catalytic RNAs to combat mis-spliced transcripts

Human transcriptome can undergo RNA mis-splicing due to spliceopathies contributing to the increasing number of genetic diseases including muscular dystrophy (MD), Alzheimer disease (AD), Huntington disease (HD), myelodysplastic syndromes (MDS). Intron retention (IR) is a major inducer of spliceopathies where two or more introns remain in the final mature mRNA and account for many intronic expansion diseases. Potential removal of such introns for therapeutic purposes can be feasible when utilizing bioinformatics, catalytic RNAs, and nano-drug delivery systems. Overcoming delivery challenges of catalytic RNAs was discussed in this review as a future perspective highlighting the significance of utilizing synthetic biology in addition to high throughput deep sequencing and computational approaches for the treatment of mis-spliced transcripts.

Mechanisms of mRNA processing defects in inherited THOC6 intellectual disability syndrome

THOC6 is the genetic basis of autosomal recessive THOC6 Intellectual Disability Syndrome (TIDS). THOC6 facilitates the formation of the Transcription Export complex (TREX) tetramer, composed of four THO monomers. The TREX tetramer supports mammalian mRNA processing that is distinct from yeast TREX dimer functions. Human and mouse TIDS model systems allow novel THOC6-dependent TREX tetramer functions to be investigated. Biallelic loss-of-functon(LOF) THOC6 variants do not influence the expression and localization of TREX members in human cells, but our data suggests reduced binding affinity of ALYREF. Impairment of TREX nuclear export functions were not detected in cells with biallelic THOC6 LOF. Instead, mRNA mis-splicing was observed in human and mouse neural tissue, revealing novel insights into THOC6-mediated TREX coordination of mRNA processing. We demonstrate that THOC6 is required for regulation of key signaling pathways in human corticogenesis that dictate the transition from proliferative to neurogenic divisions that may inform TIDS neuropathology.

Improved enzymatic labeling of fluorescent in situ hybridization probes applied to the visualization of retained introns in cells

Fluorescence in situ hybridization (FISH) is a widely used tool for quantifying gene expression and determining the location of RNA molecules in cells. We present an improved method for FISH probe production that yields high-purity probes with a wide range of fluorophores using standard laboratory equipment at low cost. The method modifies an earlier protocol that uses terminal deoxynucleotidyl transferase to add fluorescently labeled nucleotides to synthetic deoxyoligonucleotides. In our protocol, amino-11-ddUTP is joined to an oligonucleotide pool prior to its conjugation to a fluorescent dye, thereby generating pools of probes ready for a variety of modifications. This order of reaction steps allows for high labeling efficiencies regardless of the GC content or terminal base of the oligonucleotides. The degree of labeling (DOL) for spectrally distinct fluorophores (Quasar, ATTO, and Alexa dyes) was mostly >90%, comparable with commercial probes. The ease and low cost of production allowed the generation of probe sets targeting a wide variety of RNA molecules. Using these probes, FISH assays in C2C12 cells showed the expected subcellular localization of mRNAs and pre-mRNAs for Polr2a (RNA polymerase II subunit 2a) and Gapdh, and of the long noncoding RNAs Malat1 and Neat1 Developing FISH probe sets for several transcripts containing retained introns, we found that retained introns in the Gabbr1 and Noc2l transcripts are present in subnuclear foci separate from their sites of synthesis and partially coincident with nuclear speckles. This labeling protocol should have many applications in RNA biology.

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.

Evidence for the role of transcription factors in the co-transcriptional regulation of intron retention

BACKGROUND

Alternative splicing is a widespread regulatory phenomenon that enables a single gene to produce multiple transcripts. Among the different types of alternative splicing, intron retention is one of the least explored despite its high prevalence in both plants and animals. The recent discovery that the majority of splicing is co-transcriptional has led to the finding that chromatin state affects alternative splicing. Therefore, it is plausible that transcription factors can regulate splicing outcomes.

RESULTS

We provide evidence for the hypothesis that transcription factors are involved in the regulation of intron retention by studying regions of open chromatin in retained and excised introns. Using deep learning models designed to distinguish between regions of open chromatin in retained introns and non-retained introns, we identified motifs enriched in IR events with significant hits to known human transcription factors. Our model predicts that the majority of transcription factors that affect intron retention come from the zinc finger family. We demonstrate the validity of these predictions using ChIP-seq data for multiple zinc finger transcription factors and find strong over-representation for their peaks in intron retention events.

CONCLUSIONS

This work opens up opportunities for further studies that elucidate the mechanisms by which transcription factors affect intron retention and other forms of splicing.

AVAILABILITY

Source code available at https://github.com/fahadahaf/chromir.

Comparative transcriptome responses of leaf and root tissues to salt stress in wheat strains with different salinity tolerances

Background: Salinity stress is a major adverse environmental factor that can limit crop yield and restrict normal land use. The selection of salt-tolerant strains and elucidation of the underlying mechanisms by plant breeding scientists are urgently needed to increase agricultural production in arid and semi-arid regions. Results: In this study, we selected the salt-tolerant wheat (Triticum aestivum) strain ST9644 as a model to study differences in expression patterns between salt-tolerant and salt-sensitive strains. High-throughput RNA sequencing resulted in more than 359.10 Gb of clean data from 54 samples, with an average of 6.65 Gb per sample. Compared to the IWGSC reference annotation, we identified 50,096 new genes, 32,923 of which have functional annotations. Comparisons of abundances between salt-tolerant and salt-sensitive strains revealed 3,755, 5,504, and 4,344 genes that were differentially expressed at 0, 6, and 24 h, respectively, in root tissue under salt stress. KEGG pathway analysis of these genes showed that they were enriched for phenylpropanoid biosynthesis (ko00940), cysteine and methionine metabolism (ko00270), and glutathione metabolism (ko00480). We also applied weighted gene co-expression network analysis (WGCNA) analysis to determine the time course of root tissue response to salt stress and found that the acute response lasts >6 h and ends before 12 h. We also identified key alternative splicing factors showing different splicing patterns in salt-sensitive and salt-tolerant strains; however, only few of them were differentially expressed in the two groups. Conclusion: Our results offer a better understanding of wheat salt tolerance and improve wheat breeding.

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.

Stress-induced transcriptional readthrough into neighboring genes is linked to intron retention

Exposure to certain stresses leads to readthrough transcription. Using polyA-selected RNA-seq in mouse fibroblasts subjected to heat shock, oxidative, or osmotic stress, we found that readthrough transcription can proceed into proximal downstream genes, in a phenomenon previously termed "read-in." We found that read-in genes share distinctive genomic characteristics; they are GC-rich and extremely short , with genomic features conserved in human. Using ribosome profiling, we found that read-in genes show significantly reduced translation. Strikingly, read-in genes demonstrate marked intron retention, mostly in their first introns, which could not be explained solely by their short introns and GC-richness, features often associated with intron retention. Finally, we revealed H3K36me3 enrichment upstream to read-in genes. Moreover, demarcation of exon-intron junctions by H3K36me3 was absent in read-in first introns. Our data portray a relationship between read-in and intron retention, suggesting they may have co-evolved to facilitate reduced translation of read-in genes during stress.

Transcriptomic complexity of the human malaria parasite Plasmodium falciparum revealed by long-read sequencing

The Plasmodium falciparum human malaria parasite genome is incompletely annotated and does not accurately represent the transcriptomic diversity of this species. To address this need, we performed long-read transcriptomic sequencing. 5' capped mRNA was enriched from samples of total and nuclear-fractionated RNA from intra-erythrocytic stages and converted to cDNA library. The cDNA libraries were sequenced on PacBio and Nanopore long-read platforms. 12,495 novel isoforms were annotated from the data. Alternative 5' and 3' ends represent the majority of isoform events among the novel isoforms, with retained introns being the next most common event. The majority of alternative 5' ends correspond to genomic regions with features similar to those of the reference transcript 5' ends. However, a minority of alternative 5' ends showed markedly different features, including locations within protein-coding regions. Alternative 3' ends showed similar features to the reference transcript 3' ends, notably adenine-rich termination signals. Distinguishing features of retained introns could not be observed, except for a tendency towards shorter length and greater GC content compared with spliced introns. Expression of antisense and retained intron isoforms was detected at different intra-erythrocytic stages, suggesting developmental regulation of these isoform events. To gain insights into the possible functions of the novel isoforms, their protein-coding potential was assessed. Variants of P. falciparum proteins and novel proteins encoded by alternative open reading frames suggest that P. falciparum has a greater proteomic repertoire than the current annotation. We provide a catalog of annotated transcripts and encoded alternative proteins to support further studies on gene and protein regulation of this pathogen.

Integrative analysis in Pinus revealed long-term heat stress splicing memory

Due to the current climate change, many studies have described main drivers in abiotic stress. Recent findings suggest that alternative splicing (AS) has a critical role in controlling plant responses to high temperature. AS is a mechanism that allows organisms to create an assortment of RNA transcripts and proteins using a single gene. However, the most important roles of AS in stress could not be rigorously addressed because research has been focused on model species, covering only a narrow phylogenetic and lifecycle spectrum. Thus, AS degree of diversification among more dissimilar taxa in heat response is still largely unknown. To fill this gap, the present study employs a systems biology approach to examine how the AS landscape responds to and 'remembers' heat stress in conifers, a group which has received little attention even though their position can solve key evolutionary questions. Contrary to angiosperms, we found that potential intron retention may not be the most prevalent type of AS. Furthermore, our integrative analysis with metabolome and proteome data places splicing as the main source of variation during the response. Finally, we evaluated possible acquired long-term splicing memory in a diverse subset of events, and although this mechanism seems to be conserved in seed plants, AS dynamics are divergent. These discoveries reveal the particular way of remembering past temperature changes in long-lived plants and open the door to include species with unique features to determine the extent of conservation in gene expression regulation.

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.

Dynamic Transcriptional Landscape of Grass Carp (Ctenopharyngodon idella) Reveals Key Transcriptional Features Involved in Fish Development

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

A polymorphic transcriptional regulatory domain in the amyotrophic lateral sclerosis risk gene CFAP410 correlates with differential isoform expression

We describe the characterisation of a variable number tandem repeat (VNTR) domain within intron 1 of the amyotrophic lateral sclerosis (ALS) risk gene CFAP410 (Cilia and flagella associated protein 410) (previously known as C21orf2), providing insight into how this domain could support differential gene expression and thus be a modulator of ALS progression or risk. We demonstrated the VNTR was functional in a reporter gene assay in the HEK293 cell line, exhibiting both the properties of an activator domain and a transcriptional start site, and that the differential expression was directed by distinct repeat number in the VNTR. These properties embedded in the VNTR demonstrated the potential for this VNTR to modulate CFAP410 expression. We extrapolated these findings in silico by utilisation of tagging SNPs for the two most common VNTR alleles to establish a correlation with endogenous gene expression. Consistent with in vitro data, CFAP410 isoform expression was found to be variable in the brain. Furthermore, although the number of matched controls was low, there was evidence for one specific isoform being correlated with lower expression in those with ALS. To address if the genotype of the VNTR was associated with ALS risk, we characterised the variation of the CFAP410 VNTR in ALS cases and matched controls by PCR analysis of the VNTR length, defining eight alleles of the VNTR. No significant difference was observed between cases and controls, we noted, however, the cohort was unlikely to contain sufficient power to enable any firm conclusion to be drawn from this analysis. This data demonstrated that the VNTR domain has the potential to modulate CFAP410 expression as a regulatory element that could play a role in its tissue-specific and stimulus-inducible regulation that could impact the mechanism by which CFAP410 is involved in ALS.

Acanthamoeba castellanii exhibits intron retention during encystment

Intron retention (IR) refers to the mechanism of alternative splicing in which an intron is not excised from the mature transcript. IR in the cosmopolitan free-living amoeba Acanthamoeba castellanii has not been studied. We performed an analysis of RNA sequencing data during encystment to identify genes that presented differentially retained introns during this process. We show that IR increases during cyst formation, indicating a potential mechanism of gene regulation that could help downregulate metabolism. We identify 69 introns from 67 genes that are differentially retained comparing the trophozoite stage and encystment after 24 and 48 h. These genes include several hypothetical proteins. We show different patterns of IR during encystment taking as examples a lipase, a peroxin-3 protein, an Fbox domain containing protein, a proteasome subunit, a polynucleotide adenylyltransferase, and a tetratricopeptide domain containing protein. A better understanding of IR in Acanthamoeba, and even other protists, could help elucidate changes in life cycle and combat disease such as Acanthamoeba keratitis in which the cyst is key for its persistence.

Genome-Wide Analysis of the Type-B Authentic Response Regulator Gene Family in Brassica napus

The type-B authentic response regulators (type-B ARRs) are positive regulators of cytokinin signaling and involved in plant growth and stress responses. In this study, we used bioinformatics, RNA-seq, and qPCR to study the phylogenetic and expression pattern of 35 type-B ARRs in Brassica napus. The BnARRs experienced gene expansion and loss during genome polyploidization and were classified into seven groups. Whole-genome duplication (WGD) and segmental duplication were the main forces driving type-B ARR expansion in B. napus. Several BnARRs with specific expression patterns during rapeseed development were identified, including BnARR12/14/18/23/33. Moreover, we found the type-B BnARRs were involved in rapeseed development and stress responses, through participating in cytokinin and ABA signaling pathways. This study revealed the origin, evolutionary history, and expression pattern of type-B ARRs in B. napus and will be helpful to the functional characterization of BnARRs.

Intronic Polyadenylation in Acquired Cancer Drug Resistance Circumvented by Utilizing CRISPR/Cas9 with Homology-Directed Repair: The Tale of Human DNA Topoisomerase IIα

Simple Summary

DNA topoisomerase IIα (170 kDa, TOP2α/170) resolves nucleic acid topological entanglements by generating transient double-strand DNA breaks. TOP2α inhibitors/poisons stabilize TOP2α-DNA covalent complexes resulting in persistent DNA damage and are frequently utilized to treat a variety of cancers. Acquired resistance to these chemotherapeutic agents is often associated with decreased TOP2α/170 expression levels. Studies have demonstrated that a reduction in TOP2α/170 results from a type of alternative polyadenylation designated intronic polyadenylation (IPA). As a consequence of IPA, variant TOP2α mRNA transcripts have been characterized that have resulted in the translation of C-terminal truncated TOP2α isoforms with altered biological activities. In this paper, an example is discussed where circumvention of acquired TOP2α-mediated drug resistance was achieved by utilizing CRISPR/Cas9 specific gene editing of an exon/intron boundary through homology directed repair (HDR) to reduce TOP2α IPA. These results illustrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant IPA.

Abstract

Intronic polyadenylation (IPA) plays a critical role in malignant transformation, development, progression, and cancer chemoresistance by contributing to transcriptome/proteome alterations. DNA topoisomerase IIα (170 kDa, TOP2α/170) is an established clinical target for anticancer agents whose efficacy is compromised by drug resistance often associated with a reduction of nuclear TOP2α/170 levels. In leukemia cell lines with acquired resistance to TOP2α-targeted drugs and reduced TOP2α/170 expression, variant TOP2α mRNA transcripts have been reported due to IPA that resulted in the translation of C-terminal truncated isoforms with altered nuclear-cytoplasmic distribution or heterodimerization with wild-type TOP2α/170. This review provides an overview of the various mechanisms regulating pre-mRNA processing and alternative polyadenylation, as well as the utilization of CRISPR/Cas9 specific gene editing through homology directed repair (HDR) to decrease IPA when splice sites are intrinsically weak or potentially mutated. The specific case of TOP2α exon 19/intron 19 splice site editing is discussed in etoposide-resistant human leukemia K562 cells as a tractable strategy to circumvent acquired TOP2α-mediated drug resistance. This example supports the importance of aberrant IPA in acquired drug resistance to TOP2α-targeted drugs. In addition, these results demonstrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant splicing/polyadenylation.

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.

Placental Transcription Profiling in 6-23 Weeks' Gestation Reveals Differential Transcript Usage in Early Development

The human placenta is a rapidly developing transient organ that is key to pregnancy success. Early development of the conceptus occurs in a low oxygen environment before oxygenated maternal blood begins to flow into the placenta at ~10-12 weeks' gestation. This process is likely to substantially affect overall placental gene expression. Transcript variability underlying gene expression has yet to be profiled. In this study, accurate transcript expression profiles were identified for 84 human placental chorionic villus tissue samples collected across 6-23 weeks' gestation. Differential gene expression (DGE), differential transcript expression (DTE) and differential transcript usage (DTU) between 6-10 weeks' and 11-23 weeks' gestation groups were assessed. In total, 229 genes had significant DTE yet no significant DGE. Integration of DGE and DTE analyses found that differential expression patterns of individual transcripts were commonly masked upon aggregation to the gene-level. Of the 611 genes that exhibited DTU, 534 had no significant DGE or DTE. The four most significant DTU genes ADAM10, VMP1, GPR126, and ASAH1, were associated with hypoxia-responsive pathways. Transcript usage is a likely regulatory mechanism in early placentation. Identification of functional roles will facilitate new insight in understanding the origins of pregnancy complications.

A Severe Dementia Syndrome Caused by Intron Retention and Cryptic Splice Site Activation in STUB1 and Exacerbated by TBP Repeat Expansions

Heterozygous pathogenic variants in the STIP1 homologous and U-box containing protein 1 (STUB1) gene have been identified as causes of autosomal dominant inherited spinocerebellar ataxia type 48 (SCA48). SCA48 is characterized by an ataxic movement disorder that is often, but not always, accompanied by a cognitive affective syndrome. We report a severe early onset dementia syndrome that mimics frontotemporal dementia and is caused by the intronic splice donor variant c.524+1G>A in STUB1. Impaired splicing was demonstrated by RNA analysis and in minigene assays of mutated and wild-type constructs of STUB1. The most striking consequence of this splicing impairment was retention of intron 3 in STUB1, which led to an in-frame insertion of 63 amino acids (aa) (p.Arg175_Glu176ins63) into the highly conserved coiled-coil domain of its encoded protein, C-terminus of HSP70-interacting protein (CHIP). To a lesser extent, activation of two cryptic splice sites in intron 3 was observed. The almost exclusively used one, c.524+86, was not predicted by in silico programs. Variant c.524+86 caused a frameshift (p.Arg175fs*93) that resulted in a truncated protein and presumably impairs the C-terminal U-box of CHIP, which normally functions as an E3 ubiquitin ligase. The cryptic splice site c.524+99 was rarely used and led to an in-frame insertion of 33 aa (p.Arg175_Glu176ins33) that resulted in disruption of the coiled-coil domain, as has been previously postulated for complete intron 3 retention. We additionally detected repeat expansions in the range of reduced penetrance in the TATA box-binding protein (TBP) gene by excluding other genes associated with dementia syndromes. The repeat expansion was heterozygous in one patient but compound heterozygous in the more severely affected patient. Therefore, we concluded that the observed severe dementia syndrome has a digenic background, making STUB1 and TBP important candidate genes responsible for early onset dementia syndromes.

Secondhand smoke affects reproductive functions by altering the mouse testis transcriptome, and leads to select intron retention in Pde1a

BACKGROUND

Human exposure to secondhand smoke (SHS) is known to result in adverse effects in multiple organ systems. However, the impact of SHS on the male reproductive system, particularly on the regulation of genes and molecular pathways that govern sperm production, maturation, and functions remains largely understudied.

OBJECTIVE

We investigated the effects of SHS on the testis transcriptome in a validated mouse model.

METHODS

Adult male mice were exposed to SHS (5 h/day, 5 days/week for 4 months) as compared to controls (clean air-exposed). RNA-seq analysis was performed on the testis of SHS-exposed mice and controls. Variant discovery and plink association analyses were also conducted to detect exposure-related transcript variants in SHS-treated mice.

RESULTS

Exposure of mice to SHS resulted in the aberrant expression of 131 testicular genes. Whilst approximately two thirds of the differentially expressed genes were protein-coding, the remaining (30.5%) comprised noncoding elements, mostly lncRNAs (19.1%). Variant discovery analysis identified a homozygous frameshift variant that is statistically significantly associated with SHS exposure (P = 7.744e-06) and is generated by retention of a short intron within Pde1a, a key regulator of spermatogenesis. Notably, this SHS-associated intron variant harbors an evolutionarily conserved, premature termination codon (PTC) that disrupts the open reading frame of Pde1a, presumably leading to its degradation via nonsense-mediated decay.

DISCUSSION

SHS alters the expression of genes involved in molecular pathways that are crucial for normal testis development and function. Preferential targeting of lncRNAs in the testis of SHS-exposed mice is especially significant considering their crucial role in the spatial and temporal modulation of spermatogenesis. Equally important is our discovery of a novel homozygous frameshift variant that is exclusively and significantly associated with SHS-exposure and is likely to represent a safeguard mechanism to regulate transcription of Pde1a and preserve normal testis function during harmful exposure to environmental agents.

Splice variants denote differences between a cancer stem cell side population of EWSR1‑ERG‑based Ewing sarcoma cells, its main population and EWSR1‑FLI‑based cells

Ewing sarcoma is a challenging cancer entity, which, besides the characteristic presence of a fusion gene, is driven by multiple alternative splicing events. So far, splice variants in Ewing sarcoma cells were mainly analyzed for EWSR1‑FLI1. The present study provided a comprehensive alternative splicing study on CADO‑ES1, an Ewing model cell line for an EWSR1‑ERG fusion gene. Based on a well‑-characterized RNA‑sequencing dataset with extensive control mechanisms across all levels of analysis, the differential spliced genes in Ewing cancer stem cells were ATP13A3 and EPB41, while the main population was defined by ACADVL, NOP58 and TSPAN3. All alternatively spliced genes were further characterized by their Gene Ontology (GO) terms and by their membership in known protein complexes. These results confirm and extend previous studies towards a systematic whole‑transcriptome analysis. A highlight is the striking segregation of GO terms associated with five basic splice events. This mechanistic insight, together with a coherent integration of all observations with prior knowledge, indicates that EWSR1‑ERG is truly a close twin to EWSR1‑FLI1, but still exhibits certain individuality. Thus, the present study provided a measure of variability in Ewing sarcoma, whose understanding is essential both for clinical procedures and basic mechanistic insight.

Switched alternative splicing events as attractive features in lung squamous cell carcinoma

Background

Alternative splicing (AS) plays important roles in transcriptome and proteome diversity. Its dysregulation has a close affiliation with oncogenic processes. This study aimed to evaluate AS-based biomarkers by machine learning algorithms for lung squamous cell carcinoma (LUSC) patients.

Method

The Cancer Genome Atlas (TCGA) database and TCGA SpliceSeq database were utilized. After data composition balancing, Boruta feature selection and Spearman correlation analysis were used for differentially expressed AS events. Random forests and a nested fivefold cross-validation were applied for lymph node metastasis (LNM) classifier building. Random survival forest combined with Cox regression model was performed for a prognostic model, based on which a nomogram was developed. Functional enrichment analysis and Spearman correlation analysis were also conducted to explore underlying mechanisms. The expression of some switch-involved AS events along with parent genes was verified by qRT-PCR with 20 pairs of normal and LUSC tissues.

Results

We found 16 pairs of splicing events from same parent genes which were strongly related to the splicing switch (intrapair correlation coefficient = − 1). Next, we built a reliable LNM classifier based on 13 AS events as well as a nice prognostic model, in which switched AS events behaved prominently. The qRT-PCR presented consistent results with previous bioinformatics analysis, and some AS events like ITIH5-10715-AT and QKI-78404-AT showed remarkable detection efficiency for LUSC.

Conclusion

AS events, especially switched ones from the same parent genes, could provide new insights into the molecular diagnosis and therapeutic drug design of LUSC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02429-2.

Genome-Wide CRISPR Screening to Identify Mammalian Factors that Regulate Intron Retention

Intron retention (IR) regulates gene expression to control fundamental biological processes like metabolism, differentiation, and cell cycle. Despite a wide variety of genes controlled by IR, few techniques are available to identify regulators of IR in an unbiased manner. Here, we describe a CRISPR knockout screening method that can be applied to uncover regulators of IR. This method uses GFP reporter constructs containing a retained intron from a gene of interest such that GFP signal is regulated by IR in the same fashion as the endogenous gene. The GFP levels are then used as a readout for genome-wide CRISPR screening. We have successfully used this approach to identify novel regulator of IR of the MAT2A transcript and propose that similar screens will be broadly applicable for the identification of novel factors that control IR of specific transcripts.

Family-Based Whole-Exome Analysis of Specific Language Impairment (SLI) Identifies Rare Variants in BUD13, a Component of the Retention and Splicing (RES) Complex

Specific language impairment (SLI) is a common neurodevelopmental disorder (NDD) that displays high heritability estimates. Genetic studies have identified several loci, but the molecular basis of SLI remains unclear. With the aim to better understand the genetic architecture of SLI, we performed whole-exome sequencing (WES) in a single family (ID: 489; n = 11). We identified co-segregating rare variants in three new genes: BUD13, APLP2, and NDRG2. To determine the significance of these genes in SLI, we Sanger sequenced all coding regions of each gene in unrelated individuals with SLI (n = 175). We observed 13 additional rare variants in 18 unrelated individuals. Variants in BUD13 reached genome-wide significance (p-value < 0.01) upon comparison with similar variants in the 1000 Genomes Project, providing gene level evidence that BUD13 is involved in SLI. Additionally, five BUD13 variants showed cohesive variant level evidence of likely pathogenicity. Bud13 is a component of the retention and splicing (RES) complex. Additional supportive evidence from studies of an animal model (loss-of-function mutations in BUD13 caused a profound neural phenotype) and individuals with an NDD phenotype (carrying a CNV spanning BUD13), indicates BUD13 could be a target for investigation of the neural basis of language.

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.

A Deep Exon Cryptic Splice Site Promotes Aberrant Intron Retention in a Von Willebrand Disease Patient

A translationally silent single nucleotide mutation in exon 44 (E44) of the von Willebrand factor (VWF) gene is associated with inefficient removal of intron 44 in a von Willebrand disease (VWD) patient. This intron retention (IR) event was previously attributed to reordered E44 secondary structure that sequesters the normal splice donor site. We propose an alternative mechanism: the mutation introduces a cryptic splice donor site that interferes with the function of the annotated site to favor IR. We evaluated both models using minigene splicing reporters engineered to vary in secondary structure and/or cryptic splice site content. Analysis of splicing efficiency in transfected K562 cells suggested that the mutation-generated cryptic splice site in E44 was sufficient to induce substantial IR. Mutations predicted to vary secondary structure at the annotated site also had modest effects on IR and shifted the balance of residual splicing between the cryptic site and annotated site, supporting competition among the sites. Further studies demonstrated that introduction of cryptic splice donor motifs at other positions in E44 did not promote IR, indicating that interference with the annotated site is context dependent. We conclude that mutant deep exon splice sites can interfere with proper splicing by inducing IR.

CYP11B1 variants influence skeletal maturation via alternative splicing

We performed genome-wide association study meta-analysis to identify genetic determinants of skeletal age (SA) deviating in multiple growth disorders. The joint meta-analysis (N = 4557) in two multiethnic cohorts of school-aged children identified one locus, CYP11B1 (expression confined to the adrenal gland), robustly associated with SA (rs6471570-A; β = 0.14; P = 6.2 × 10^-12). rs6410 (a synonymous variant in the first exon of CYP11B1 in high LD with rs6471570), was prioritized for functional follow-up being second most significant and the one closest to the first intron-exon boundary. In 208 adrenal RNA-seq samples from GTEx, C-allele of rs6410 was associated with intron 3 retention (P = 8.11 × 10^-40), exon 4 inclusion (P = 4.29 × 10^-34), and decreased exon 3 and 5 splicing (P = 7.85 × 10^-43), replicated using RT-PCR in 15 adrenal samples. As CYP11B1 encodes 11-β-hydroxylase, involved in adrenal glucocorticoid and mineralocorticoid biosynthesis, our findings highlight the role of adrenal steroidogenesis in SA in healthy children, suggesting alternative splicing as a likely underlying mechanism.

Modulation of alternative splicing during early infection of human primary B lymphocytes with Epstein-Barr virus (EBV): a novel function for the viral EBNA-LP protein

Epstein-Barr virus (EBV) is a human herpesvirus associated with human cancers worldwide. Ex vivo, the virus efficiently infects resting human B lymphocytes and induces their continuous proliferation. This process is accompanied by a global reprogramming of cellular gene transcription. However, very little is known on the impact of EBV infection on the regulation of alternative splicing, a pivotal mechanism that plays an essential role in cell fate determination and is often deregulated in cancer. In this study, we have developed a systematic time-resolved analysis of cellular mRNA splice variant expression during EBV infection of resting B lymphocytes. Our results reveal that major modifications of alternative splice variant expression appear as early as day 1 post-infection and suggest that splicing regulation provides—besides transcription—an additional mechanism of gene expression regulation at the onset of B cell activation and proliferation. We also report a role for the viral proteins, EBNA2 and EBNA-LP, in the modulation of specific alternative splicing events and reveal a previously unknown function for EBNA-LP—together with the RBM4 splicing factor—in the alternative splicing regulation of two important modulators of cell proliferation and apoptosis respectively, NUMB and BCL-X.

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.

Alternative Splicing: A New Cause and Potential Therapeutic Target in Autoimmune Disease

Alternative splicing (AS) is a complex coordinated transcriptional regulatory mechanism. It affects nearly 95% of all protein-coding genes and occurs in nearly all human organs. Aberrant alternative splicing can lead to various neurological diseases and cancers and is responsible for aging, infection, inflammation, immune and metabolic disorders, and so on. Though aberrant alternative splicing events and their regulatory mechanisms are widely recognized, the association between autoimmune disease and alternative splicing has not been extensively examined. Autoimmune diseases are characterized by the loss of tolerance of the immune system towards self-antigens and organ-specific or systemic inflammation and subsequent tissue damage. In the present review, we summarized the most recent reports on splicing events that occur in the immunopathogenesis of systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) and attempted to clarify the role that splicing events play in regulating autoimmune disease progression. We also identified the changes that occur in splicing factor expression. The foregoing information might improve our understanding of autoimmune diseases and help develop new diagnostic and therapeutic tools for them.

Extensive Changes in Transcription Dynamics Reflected on Alternative Splicing Events in Systemic Lupus Erythematosus Patients

In addition to increasing the complexity of the transcriptional output, alternative RNA splicing can lead to the reduction of mRNA translation or the production of non-functional or malfunctional proteins, thus representing a vital component of the gene regulation process. Herein, we set out to detect and characterize alternative splicing events that occur in whole-blood samples of patients with Systemic Lupus Erythematosus (SLE) as compared to healthy counterparts. Through the implementation of a computational pipeline on published RNA-sequencing data, we identified extensive changes in the transcription dynamics affecting a large number of genes. We found a predominance of intron retention events, with the majority introducing premature stop codons, suggestive of gene repression, in both inactive and active SLE patient samples. Alternative splicing affected a distinct set of genes from the ones detected as differentially expressed in the same comparisons, while alternatively spliced genes tended to reside in genome areas associated with increased gene co-expression. Functional analysis of genes affected by alternative splicing pointed towards particular functions related to metabolism and histone acetylation as of potential interest. Together, our findings underline the importance of incorporating alternative splicing analyses in the context of molecular characterization of complex diseases such as SLE.

Different Patterns of mRNA Nuclear Retention during Meiotic Prophase in Larch Microsporocytes

Recent studies show a crucial role of post-transcriptional processes in the regulation of gene expression. Our research has shown that mRNA retention in the nucleus plays a significant role in such regulation. We studied larch microsporocytes during meiotic prophase, characterized by pulsatile transcriptional activity. After each pulse, the transcriptional activity is silenced, but the transcripts synthesized at this time are not exported immediately to the cytoplasm but are retained in the cell nucleus and especially in Cajal bodies, where non-fully-spliced transcripts with retained introns are accumulated. Analysis of the transcriptome of these cells and detailed analysis of the nuclear retention and transport dynamics of several mRNAs revealed two main patterns of nuclear accumulation and transport. The majority of studied transcripts followed the first one, consisting of a more extended retention period and slow release to the cytoplasm. We have shown this in detail for the pre-mRNA and mRNA encoding RNA pol II subunit 10. In this pre-mRNA, a second (retained) intron is posttranscriptionally spliced at a precisely defined time. Fully mature mRNA is then released into the cytoplasm, where the RNA pol II complexes are produced. These proteins are necessary for transcription in the next pulse to occur.mRNAs encoding translation factors and SERRATE followed the second pattern, in which the retention period was shorter and transcripts were rapidly transferred to the cytoplasm. The presence of such a mechanism in various cell types from a diverse range of organisms suggests that it is an evolutionarily conserved mechanism of gene regulation.

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.