Unmasking alternative splicing inside protein-coding exons defines exitrons and their role in proteome plasticity

New dimensions in the molecular genetics of insect chemoreception

Chemoreception is the foundation of olfaction and taste, which in insects underlie the detection of humans to whom they spread disease and crops that they ravage. Recent advances have provided clear and in some cases surprising new insights into the molecular genetics of chemoreception. We describe mechanisms that govern the choice of a single Odorant receptor gene by an olfactory receptor neuron in Drosophila. We highlight genetic and epigenetic mechanisms by which chemoreceptor expression can be modulated. Exitrons, RNA editing, and pseudo-pseudogenes in chemosensory systems are described. We summarize key insights from the recent structural determinations of odorant and taste receptors. Finally, new molecular components of chemosensory systems, including long noncoding RNAs, are described.

The Exserohilum turcicum effector EtEC81 reprograms alternative splicing in maize and activates immunity

Some pathogen-derived effectors reprogram mRNA splicing in host plants to regulate plant immune responses. Whether effectors from Exserohilum turcicum, which causes northern corn leaf blight (NLB), interfere with RNA splicing remains unknown. We identify that the secreted protein EtEC81 (Exserohilum turcicum effector 81) modulates the alternative splicing (AS) of maize (Zea mays) pre-mRNAs and negatively regulates the pathogenicity of E. turcicum. EtEC81 physically interacts with MAIZE EtEC81-INTERACTING PROTEIN 1 (ZmEIP1), which associates with maize spliceosome components, modulates AS in host cells, and positively regulates defense responses against E. turcicum. Transcriptome analysis identifies 119 common events with altered AS in maize plants transiently overexpressing ZmEIP1 or EtEC81, suggesting that these factors cause the misregulation of cellular activities and thus induce immune responses. Together, our results suggest that the EtEC81 effector targets ZmEIP1 to reprogram pre-mRNA splicing in maize. These findings provide a mechanistic basis and potential target gene for preventing NLB.

Global dysregulation of circular RNAs in frontal cortex and whole blood from DM1 and DM2

Myotonic dystrophy type 1 (DM1) and type 2 (DM2) are autosomal dominant neuromuscular disorders associated with expansions of microsatellites, respectively, in DMPK and CNBP. Their pathogenesis is linked to the global aberrant alternative splicing (AAS) of many genes and marks mostly muscular and neuronal tissues, while blood is the least affected. Recent data in DM1 skeletal muscles indicated that abnormalities in RNA metabolism also include global upregulation of circular RNAs (circRNAs). CircRNAs are a heterogeneous group considered splicing errors and by-products of canonical splicing. To elucidate whether circRNA dysregulation is an inherent feature of the myotonic environment, we perform their analysis in the frontal cortex and whole blood of DM1 and DM2 patients. We find a global elevation of circRNAs in both tissues, and its magnitude is neither correlated with the differences in their parental gene expression nor is associated with AAS published earlier. Aberrantly spliced cassette exons of linear transcripts affected in DM1 and DM2 are not among the circularized exons, which unique genomic features prerequisite back-splicing. However, the blueprint of the AAS of linear RNAs is found in a variety of circRNA isoforms. The heterogeneity of circRNAs also originates from the utilization of exonic and intronic cryptic donors/acceptors in back splice junctions, and intron-containing circRNAs are more characteristic of the blood. Overall, this study reveals circRNA dysregulation in various tissues from DM1 and DM2; however, their levels do not correlate with the AAS in linear RNAs, suggesting a potential independent regulatory mechanism underlying circRNA upregulation in myotonic dystrophy.

Characterisation of LGP2 complex multitranscript system in humans: role in the innate immune response and evolution from non-human primates

Retinoic acid inducible gene I (RIG-I)-like receptors (RLRs), including RIG-I, MDA5 and LGP2, recognize viral RNA to mount an antiviral interferon (IFN) response RLRs share three different protein domains: C-terminal domain, DExD/H box RNA helicase domain, and an N-terminal domain with two tandem repeats (CARDs). LGP2 lacks tandem CARD and is not able to induce an IFN response. However, LGP2 positively enhances MDA5 and negatively regulates RIG-I signaling. In this study, we determined the LGP2 alternative transcripts in humans to further comprehend the mechanism of its regulation, their evolutionary origin, and the isoforms functionallity. The results showed new eight alternative transcripts in the samples tested. The presence of these transcripts demonstrated that the main mechanisms for the regulation of LGP2 expression are both by insertion of introns and by the loss of exons. The phylogenetic analysis of the comparison between sequences from exon 1 to exon 3 of humans and those previously described in non-human primates showed three well-differentiated groups (lineages) originating from gorillas, suggesting that the transspecies evolution has been maintained for 10 million years. The corresponding protein models (isoforms) were also established, obtaining four isoforms: one complete and three others lacking the C-terminal domain or this domain and the partial or total He2 Helicase domain, which would compromise the functionality of LGP2. In conclusion, this is the first study that elucidate the large genomic organization and complex transcriptional regulation of human LGP2, its pattern of sequence generation, and a mode of evolutionary inheritance across species.

At-RS31 orchestrates hierarchical cross-regulation of splicing factors and integrates alternative splicing with TOR-ABA pathways

Alternative splicing is essential for plants, enabling a single gene to produce multiple transcript variants to boost functional diversity and fine-tune responses to environmental and developmental cues. At-RS31, a plant-specific splicing factor in the Serine/Arginine (SR)-rich protein family, responds to light and the Target of Rapamycin (TOR) signaling pathway, yet its downstream targets and regulatory impact remain unknown.To identify At-RS31 targets, we applied individual-nucleotide resolution crosslinking and immunoprecipitation (iCLIP) and RNAcompete assays. Transcriptomic analyses of At-RS31 mutant and overexpressing plants further revealed its effects on alternative splicing.iCLIP identified 4,034 At-RS31 binding sites across 1,421 genes, enriched in CU-rich and CAGA RNA motifs. Comparative iCLIP and RNAcompete data indicate that the RS domain of At-RS31 may influence its binding specificity in planta, underscoring the value of combining in vivo and in vitro approaches. Transcriptomic analysis showed that At-RS31 modulates diverse splicing events, particularly intron retention and exitron splicing, and influences other splicing modulators, acting as a hierarchical regulator.By regulating stress-response genes and genes in both TOR and abscisic acid (ABA) signaling pathways, At-RS31 may help integrate these signals, balancing plant growth with environmental adaptability through alternative splicing.

RNA Isoform Diversity in Human Neurodegenerative Diseases

Single-nucleus RNA-sequencing (snRNA-seq) has revealed new levels of cellular organization and diversity within the human brain. However, full-length mRNA isoforms are not resolved in typical snRNA-seq analyses using short-read sequencing that cannot capture full-length transcripts. Here we combine standard 10x Genomics short-read snRNA-seq with targeted PacBio long-read snRNA-seq to examine isoforms of genes associated with neurological diseases at the single-cell level from prefrontal cortex samples of diseased and nondiseased human brain, assessing over 165,000 cells. Samples from 25 postmortem donors with Alzheimer's disease (AD), dementia with Lewy bodies (DLB), or Parkinson's disease (PD), along with age-matched controls, were compared. Analysis of the short-read libraries identified shared and distinct gene expression changes across the diseases. The same libraries were then assayed using enrichment probes to target 50 disease-related genes followed by long-read PacBio sequencing, enabling linkage between cell type and isoform expression. Vast mRNA isoform diversity was observed in all 50 targeted genes, even those that were not differentially expressed in the short-read data. We also developed an informatics method for detection of isoform structural differences in novel isoforms versus the reference annotation. These data expand available single-cell datasets of the human prefrontal cortical transcriptome with combined short- and long-read sequencing across AD, DLB, and PD, revealing increased mRNA isoform diversity that may contribute to disease features and could potentially represent therapeutic targets for neurodegenerative diseases.

Steering research on mRNA splicing in cancer towards clinical translation

Splicing factors are affected by recurrent somatic mutations and copy number variations in several types of haematologic and solid malignancies, which is often seen as prima facie evidence that splicing aberrations can drive cancer initiation and progression. However, numerous spliceosome components also 'moonlight' in DNA repair and other cellular processes, making their precise role in cancer difficult to pinpoint. Still, few would deny that dysregulated mRNA splicing is a pervasive feature of most cancers. Correctly interpreting these molecular fingerprints can reveal novel tumour vulnerabilities and untapped therapeutic opportunities. Yet multiple technological challenges, lingering misconceptions, and outstanding questions hinder clinical translation. To start with, the general landscape of splicing aberrations in cancer is not well defined, due to limitations of short-read RNA sequencing not adept at resolving complete mRNA isoforms, as well as the shallow read depth inherent in long-read RNA-sequencing, especially at single-cell level. Although individual cancer-associated isoforms are known to contribute to cancer progression, widespread splicing alterations could be an equally important and, perhaps, more readily actionable feature of human cancers. This is to say that in addition to 'repairing' mis-spliced transcripts, possible therapeutic avenues include exacerbating splicing aberration with small-molecule spliceosome inhibitors, targeting recurrent splicing aberrations with synthetic lethal approaches, and training the immune system to recognize splicing-derived neoantigens.

An anciently diverged family of RNA binding proteins maintain correct splicing of a class of ultra-long exons through cryptic splice site repression

Previously, we showed that the germ cell-specific nuclear protein RBMXL2 represses cryptic splicing patterns during meiosis and is required for male fertility (Ehrmann et al., 2019). Here, we show that in somatic cells the similar yet ubiquitously expressed RBMX protein has similar functions. RBMX regulates a distinct class of exons that exceed the median human exon size. RBMX protein-RNA interactions are enriched within ultra-long exons, particularly within genes involved in genome stability, and repress the selection of cryptic splice sites that would compromise gene function. The RBMX gene is silenced during male meiosis due to sex chromosome inactivation. To test whether RBMXL2 might replace the function of RBMX during meiosis we induced expression of RBMXL2 and the more distantly related RBMY protein in somatic cells, finding each could rescue aberrant patterns of RNA processing caused by RBMX depletion. The C-terminal disordered domain of RBMXL2 is sufficient to rescue proper splicing control after RBMX depletion. Our data indicate that RBMX and RBMXL2 have parallel roles in somatic tissues and the germline that must have been conserved for at least 200 million years of mammalian evolution. We propose RBMX family proteins are particularly important for the splicing inclusion of some ultra-long exons with increased intrinsic susceptibility to cryptic splice site selection.

Exitron splicing of odor receptor genes in Drosophila

Proper expression of odor receptor genes is critical for the function of olfactory systems. In this study, we identified exitrons (exonic introns) in four of the 39 Odorant receptor (Or) genes expressed in the Drosophila antenna. Exitrons are sequences that can be spliced out from within a protein-coding exon, thereby altering the encoded protein. We focused on Or88a, which encodes a pheromone receptor, and found that exitron splicing of Or88a is conserved across five Drosophila species over 20 My of evolution. The exitron was spliced out in 15% of Or88a transcripts. Removal of this exitron creates a non-coding RNA rather than an RNA that encodes a stable protein. Our results suggest the hypothesis that in the case of Or88a, exitron splicing could act in neuronal modulation by decreasing the level of functional Or transcripts. Activation of Or88a-expressing olfactory receptor neurons via either optogenetics or pheromone stimulation increased the level of exitron-spliced transcripts, with optogenetic activation leading to a 14-fold increase. A fifth Or can also undergo an alternative splicing event that eliminates most of the canonical open reading frame. Besides these cases of alternative splicing, we found alternative polyadenylation of four Ors, and exposure of Or67c to its ligand ethyl lactate in the antenna downregulated all of its 3' isoforms. Our study reveals mechanisms by which neuronal activity could be modulated via regulation of the levels of Or isoforms.

Importance of pre-mRNA splicing and its study tools in plants

Alternative splicing (AS) significantly enriches the diversity of transcriptomes and proteomes, playing a pivotal role in the physiology and development of eukaryotic organisms. With the continuous advancement of high-throughput sequencing technologies, an increasing number of novel transcript isoforms, along with factors related to splicing and their associated functions, are being unveiled. In this review, we succinctly summarize and compare the different splicing mechanisms across prokaryotes and eukaryotes. Furthermore, we provide an extensive overview of the recent progress in various studies on AS covering different developmental stages in diverse plant species and in response to various abiotic stresses. Additionally, we discuss modern techniques for studying the functions and quantification of AS transcripts, as well as their protein products. By integrating genetic studies, quantitative methods, and high-throughput omics techniques, we can discover novel transcript isoforms and functional splicing factors, thereby enhancing our understanding of the roles of various splicing modes in different plant species.

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

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

Responsive Alternative Splicing Events of Opisthopappus Species against Salt Stress

Salt stress profoundly affects plant growth, prompting intricate molecular responses, such as alternative splicing (AS), for environmental adaptation. However, the response of AS events to salt stress in Opisthopappus (Opisthopappus taihangensis and Opisthopappus longilobus) remains unclear, which is a Taihang Mountain cliff-dwelling species. Using RNA-seq data, differentially expressed genes (DEGs) were identified under time and concentration gradients of salt stress. Two types of AS, skipped exon (SE) and mutually exclusive exons (MXE), were found. Differentially alternative splicing (DAS) genes in both species were significantly enriched in "protein phosphorylation", "starch and sucrose metabolism", and "plant hormone signal transduction" pathways. Meanwhile, distinct GO terms and KEGG pathways of DAS occurred between two species. Only a small subset of DAS genes overlapped with DEGs under salt stress. Although both species likely adopted protein phosphorylation to enhance salt stress tolerance, they exhibited distinct responses. The results indicated that the salt stress mechanisms of both Opisthopappus species exhibited similarities and differences in response to salt stress, which suggested that adaptive divergence might have occurred between them. This study initially provides a comprehensive description of salt responsive AS events in Opisthopappus and conveys some insights into the molecular mechanisms behind species tolerance on the Taihang Mountains.

Plant immunity suppressor SKRP encodes a novel RNA-binding protein that targets exon 3' end of unspliced RNA

The regulatory roles of RNA splicing in plant immunity are emerging but still largely obscure. We reported previously that Phytophthora pathogen effector Avr3c targets a soybean protein SKRP (serine/lysine/arginine-rich protein) to impair soybean basal immunity by regulating host pre-mRNA alternative splicing, while the biochemical nature of SKRP remains unknown. Here, by using Arabidopsis as a model, we studied the mechanism of SKRP in regulating pre-mRNA splicing and plant immunity. AtSKRP confers impaired plant immunity against Phytophthora capsici and associates with spliceosome component PRP8 and splicing factor SR45, which positively and negatively regulate plant immunity, respectively. Enhanced crosslinking and immunoprecipitation followed by high-throughput sequencing (eCLIP-seq) showed AtSKRP is a novel RNA-binding protein that targets exon 3' end of unspliced RNA. Such position-specific binding of SKRP is associated with its activity in suppressing intron retention, including at positive immune regulatory genes UBP25 and RAR1. In addition, we found AtSKRP self-interact and forms oligomer, and these properties are associated with its function in plant immunity. Overall, our findings reveal that the immune repressor SKRP is a spliceosome-associated protein that targets exon 3' end to regulate pre-mRNA splicing in Arabidopsis.

Identification of Intron Retention in the Slc16a3 Gene Transcript Encoding the Transporter MCT4 in the Brain of Aged and Alzheimer-Disease Model (APPswePS1dE9) Mice

The monocarboxylate transporter 4 (MCT4; Slc16a3) is expressed in the central nervous system, notably by astrocytes. It is implicated in lactate release and the regulation of glycolytic flux. Whether its expression varies during normal and/or pathological aging is unclear. As the presence of its mature transcript in the brain of young and old mice was determined, an unexpectedly longer RT-PCR fragment was detected in the mouse frontal cortex and hippocampus at 12 vs. 3 months of age. Cultured astrocytes expressed the expected 516 base pair (bp) fragment but treatment with IL-1β to mimic inflammation as can occur during aging led to the additional expression of a 928 bp fragment like that seen in aged mice. In contrast, cultured pericytes (a component of the blood-brain barrier) only exhibited the 516 bp fragment. Intriguingly, cultured endothelial cells constitutively expressed both fragments. When RT-PCR was performed on brain subregions of an Alzheimer mouse model (APPswePS1dE9), no fragment was detected at 3 months, while only the 928 bp fragment was present at 12 months. Sequencing of MCT4 RT-PCR products revealed the presence of a remaining intron between exon 2 and 3, giving rise to the longer fragment detected by RT-PCR. These results unravel the existence of intron retention for the MCT4 gene in the central nervous system. Such alternative splicing appears to increase with age in the brain and might be prominent in neurodegenerative diseases such as Alzheimer's disease. Hence, further studies in vitro and in vivo of intron 2 retention in the Slc16a3 gene transcript are required for adequate characterization concerning the biological roles of Slc16a3 isoforms in the context of aging and Alzheimer's disease pathology.

Cotranscriptional RNA processing and modification in plants

The activities of RNA polymerases shape the epigenetic landscape of genomes with profound consequences for genome integrity and gene expression. A fundamental event during the regulation of eukaryotic gene expression is the coordination between transcription and RNA processing. Most primary RNAs mature through various RNA processing and modification events to become fully functional. While pioneering results positioned RNA maturation steps after transcription ends, the coupling between the maturation of diverse RNA species and their transcription is becoming increasingly evident in plants. In this review, we discuss recent advances in our understanding of the crosstalk between RNA Polymerase II, IV, and V transcription and nascent RNA processing of both coding and noncoding RNAs.

Do not panic: An intron-centric guide to alternative splicing

This review is an attempt to establish concepts of splicing and alternative splicing giving proper relevance to introns, the key actors in this mechanism. It might also work as a guide for those who found their favorite gene undergoes alternative splicing and could benefit from gaining a theoretical framework to understand the possible impacts of this process. This is not a thorough review of all the work in the field, but rather a critical review of some of the most relevant work done to understand the underlying mechanisms of splicing and the key questions that remain unanswered such as: What is the physiological relevance of alternative splicing? What are the functions of the different outcomes? To what extent do different alternative splicing types contribute to the proteome? Intron retention is the most frequent alternative splicing event in plants and, although scientifically neglected, it is also common in animals. This is a heterogeneous type of alternative splicing that includes different sub-types with features that have distinctive consequences in the resulting transcripts. Remarkably, intron retention can be a dead end for a transcript, but it could also be a stable intermediate whose processing is resumed upon a particular signal or change in the cell status. New sequencing technologies combined with the study of intron lariats in different conditions might help to answer key questions and could help us to understand the actual relevance of introns in gene expression regulation.

NLRP3 leucine-rich repeats control induced and spontaneous inflammasome activation in cryopyrin-associated periodic syndrome

BACKGROUND

The cryopyrin-associated periodic syndromes (CAPS) comprise a group of rare autoinflammatory diseases caused by gain-of-function mutations in the NLRP3 gene. NLRP3 contains a leucine-rich repeats (LRR) domain with a highly conserved exonic organization that is subjected to extensive alternative splicing. Aberrant NLRP3 inflammasome assembly in CAPS causes chronic inflammation; however, the mechanisms regulating inflammasome function remain unclear.

OBJECTIVE

We aimed to elucidate the mechanisms regulating NLRP3-mediated autoinflammation in human disease, characterizing the role of LRR in inflammasome activation.

METHODS

We analyzed sequence read archive data to characterize the pattern of NLRP3 splicing in human monocytes and investigated the role of each LRR-coding exon in inflammasome regulation in genetically modified U937 cells representing CAPS and healthy conditions.

RESULTS

We detected a range of NLRP3 splice variants in human primary cells and monocytic cell lines, including 2 yet-undescribed splice variants. We observe that lipopolysaccharides affect the abundance of certain splice variants, suggesting that they may regulate NLRP3 activation by affecting alternative splicing. We showed that exons 4, 5, 7, and 9 are essential for inflammasome function, both in the context of wild-type NLRP3 activation by the agonist molecule nigericin and in a model of CAPS-mediated NLRP3 inflammasome assembly. Moreover, the SGT1-NLRP3 interaction is decreased in nonfunctional variants, suggesting that alternative splicing may regulate the recruitment of proteins that facilitate inflammasome assembly.

CONCLUSION

These findings demonstrate the contribution of the LRR domain in inflammasome function and suggest that navigating LRR exon usage within NLRP3 is sufficient to dampen inflammasome assembly in CAPS.

Exon Elongation Added Intrinsically Disordered Regions to the Encoded Proteins and Facilitated the Emergence of the Last Eukaryotic Common Ancestor

Most prokaryotic proteins consist of a single structural domain (SD) with little intrinsically disordered regions (IDRs) that by themselves do not adopt stable structures, whereas the typical eukaryotic protein comprises multiple SDs and IDRs. How eukaryotic proteins evolved to differ from prokaryotic proteins has not been fully elucidated. Here, we found that the longer the internal exons are, the more frequently they encode IDRs in eight eukaryotes including vertebrates, invertebrates, a fungus, and plants. Based on this observation, we propose the "small bang" model from the proteomic viewpoint: the protoeukaryotic genes had no introns and mostly encoded one SD each, but a majority of them were subsequently divided into multiple exons (step 1). Many exons unconstrained by SDs elongated to encode IDRs (step 2). The elongated exons encoding IDRs frequently facilitated the acquisition of multiple SDs to make the last common ancestor of eukaryotes (step 3). One prediction of the model is that long internal exons are mostly unconstrained exons. Analytical results of the eight eukaryotes are consistent with this prediction. In support of the model, we identified cases of internal exons that elongated after the rat-mouse divergence and discovered that the expanded sections are mostly in unconstrained exons and preferentially encode IDRs. The model also predicts that SDs followed by long internal exons tend to have other SDs downstream. This prediction was also verified in all the eukaryotic species analyzed. Our model accounts for the dichotomy between prokaryotic and eukaryotic proteins and proposes a selective advantage conferred by IDRs.

Transcriptomic and splicing changes underlying tomato responses to combined water and nutrient stress

Tomato is a horticultural crop of high economic and nutritional value. Suboptimal environmental conditions, such as limited water and nutrient availability, cause severe yield reductions. Thus, selection of genotypes requiring lower inputs is a goal for the tomato breeding sector. We screened 10 tomato varieties exposed to water deficit, low nitrate or a combination of both. Biometric, physiological and molecular analyses revealed different stress responses among genotypes, identifying T270 as severely affected, and T250 as tolerant to the stresses applied. Investigation of transcriptome changes caused by combined stress in roots and leaves of these two genotypes yielded a low number of differentially expressed genes (DEGs) in T250 compared to T270, suggesting that T250 tailors changes in gene expression to efficiently respond to combined stress. By contrast, the susceptible tomato activated approximately one thousand and two thousand genes in leaves and roots respectively, indicating a more generalized stress response in this genotype. In particular, developmental and stress-related genes were differentially expressed, such as hormone responsive factors and transcription factors. Analysis of differential alternative splicing (DAS) events showed that combined stress greatly affects the splicing landscape in both genotypes, highlighting the important role of AS in stress response mechanisms. In particular, several stress and growth-related genes as well as transcription and splicing factors were differentially spliced in both tissues. Taken together, these results reveal important insights into the transcriptional and post-transcriptional mechanisms regulating tomato adaptation to growth under reduced water and nitrogen inputs.

ScanExitronLR: characterization and quantification of exitron splicing events in long-read RNA-seq data

SUMMARY

Exitron splicing is a type of alternative splicing where coding sequences are spliced out. Recently, exitron splicing has been shown to increase proteome plasticity and play a role in cancer. Long-read RNA-seq is well suited for quantification and discovery of alternative splicing events; however, there are currently no tools available for the detection and annotation of exitrons in long-read RNA-seq data. Here, we present ScanExitronLR, an application for the characterization and quantification of exitron splicing events in long-reads. From a BAM alignment file, reference genome and reference gene annotation, ScanExitronLR outputs exitron events at the individual transcript level. Outputs of ScanExitronLR can be used in downstream analyses of differential exitron splicing. In addition, ScanExitronLR optionally reports exitron annotations such as truncation or frameshift type, nonsense-mediated decay status and Pfam domain interruptions. We demonstrate that ScanExitronLR performs better on noisy long-reads than currently published exitron detection algorithms designed for short-read data.

AVAILABILITY AND IMPLEMENTATION

ScanExitronLR is freely available at https://github.com/ylab-hi/ScanExitronLR and distributed as a pip package on the Python Package Index.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Landscape of exitrons in gastric cancer

BACKGROUND

Exitron is a new type of non-canonical alternative splicing. Accumulating evidence implies exitron may have pathological function and contribute to another source of anti-tumor immunogenicity in various cancers. Its role in gastric cancer remains poorly understood. Large-scale, multi-omics analysis could comprehensively characterize the landscape of exitrons in gastric cancer, reveal undiscovered mechanism and hopefully identify molecular biomarkers for predicting immunotherapy response.

METHODS

We collected datasets from five studies for analysis. RNA sequencing was used for exitron identification. Somatic mutations were detected by whole exome sequencing. Neopeptides were confirmed by proteome mass spectrometry.

FINDINGS

42174 gastric cancer-specific exitrons (GCSEs) were identified in 632 patients. GCSEs were clinically relevant to gender, age, Lauren type, tumor stage and prognosis. Tissue specificity test and pathogenic exitron prediction revealed their unique functional impact. GCSEs were mutually exclusive with mutations and demonstrated both unique and complementary function against TP53 mutation in gastric cancer. We further established splicing regulatory network to reveal upstream regulation of exitron splicing. We also evaluated the immunogenicity and diagnostic potential of GCSEs. Evidence of GCSEs-derived neopeptide expression was validated by whole proteome mass spectrometry. PD-1 and Siglecs were significantly increased in high neoantigen load patients. But exitron-related biomarkers failed to predict immunotherapy response, possibly due to small sample size and insufficient sequencing depth.

INTERPRETATION

The present study provided a comprehensive multidimensional landscape of gastric cancer exitrons and underscores insights into underexplored mechanism in gastric cancer pathology.

FUNDING

The Guangdong Provincial Key Laboratory of Precision Medicine for Gastroinstestinal Cancer (2020B121201004), the Guangdong Provincial Major Talents Project (No. 2019JC05Y361) and National Natural Science Foundation of China (grant number:82172960 and 81872013).

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

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

Alternative splicing as a key player in the fine-tuning of the immunity response in Arabidopsis

Plants, like animals, are constantly exposed to abiotic and biotic stresses, which often inhibit plant growth and development, and cause tissue damage, disease, and even plant death. Efficient and timely response to stress requires appropriate co- and posttranscriptional reprogramming of gene expression. Alternative pre-mRNA splicing provides an important layer of this regulation by controlling the level of factors involved in stress response and generating additional protein isoforms with specific features. Recent high-throughput studies have revealed that several defence genes undergo alternative splicing that is often affected by pathogen infection. Despite extensive work, the exact mechanisms underlying these relationships are still unclear, but the contribution of alternative protein isoforms to the defence response and the role of regulatory factors, including components of the splicing machinery, have been established. Modulation of gene expression in response to stress includes alternative splicing, chromatin remodelling, histone modifications, and nucleosome occupancy. How these processes affect plant immunity is mostly unknown, but these facets open new regulatory possibilities. Here we provide an overview of the current state of knowledge and recent findings regarding the growing importance of alternative splicing in plant response to biotic stress.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

BACKGROUND

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis.

RESULTS

We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage.

CONCLUSIONS

AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

The splicing factor 1-FLOWERING LOCUS M module spatially regulates temperature-dependent flowering by modulating FLOWERING LOCUS T and LEAFY expression

The AtSF1-FLM module spatially controls temperature-dependent flowering by negatively regulating the expression of FT and LFY in the leaf and shoot apex, respectively. Alternative splicing mediated by various splicing factors is important for the regulation of plant growth and development. Our recent reports have shown that a temperature-dependent interaction between Arabidopsis thaliana splicing factor 1 (AtSF1) and FLOWERING LOCUS M (FLM) pre-mRNA introns controls the differential production of FLM-β transcripts at different temperatures, eventually resulting in temperature-responsive flowering. However, the molecular and genetic interactions between the AtSF1-FLM module and floral activator genes remain unknown. Here, we aimed to identify the interactions among AtSF1, FLM, FLOWERING LOCUS T (FT), and LEAFY (LFY) by performing molecular and genetic analyses. FT and TWIN SISTER OF FT (TSF) expression in atsf1-2 mutants significantly increased in the morning and middle of the night at 16 and 23 °C, respectively, under long-day conditions. In addition, ft mutation suppressed the early flowering of atsf1-2 and atsf1-2 flm-3 mutants and masked the temperature response of atsf1-2 flm-3 mutants, suggesting that FT is a downstream target gene of the AtSF1-FLM module. LFY expression significantly increased in the diurnal samples of atsf1-2 mutants and in the shoot apex regions of atsf1-2 ft-10 mutants at different temperatures. The chromatin immunoprecipitation (ChIP) assay revealed that FLM directly binds to the genomic regions of LFY but not of APETALA1 (AP1). Moreover, lfy mutation suppressed the early flowering of flm-3 mutants, suggesting that LFY is another target of the AtSF1-FLM module. Our results reveal that the AtSF1-FLM module spatially modulates temperature-dependent flowering by regulating FT and LFY expressions.

A rust fungus effector directly binds plant pre-mRNA splice site to reprogram alternative splicing and suppress host immunity

Alternative splicing (AS) is a crucial post-transcriptional regulatory mechanism in plant resistance. However, whether and how plant pathogens target splicing in their host remains mostly unknown. For example, although infection by Puccinia striiformis f. sp. tritici (Pst), a pathogenic fungus that severely affects the yield of wheat worldwide, has been shown to significantly influence the levels of alternatively spliced transcripts in the host, the mechanisms that govern this process, and its functional consequence have not been examined. Here, we identified Pst_A23 as a new Pst arginine-rich effector that localizes to host nuclear speckles, nuclear regions enriched in splicing factors. We demonstrated that transient expression of Pst_A23 suppresses plant basal defence dependent on the Pst_A23 nuclear speckle localization and that this protein plays an important role in virulence, stable silencing of which improves wheat stripe rust resistance. Remarkably, RNA-Seq data revealed that AS patterns of 588 wheat genes are altered in Pst_A23-overexpressing lines compared to control plants. To further examine the direct relationship between Pst_A23 and AS, we confirmed direct binding between two RNA motifs predicted from these altered splicing sites and Pst_A23 in vitro. The two RNA motifs we chose occur in the cis-element of TaXa21-H and TaWRKY53, and we validated that Pst_A23 overexpression results in decreased functional transcripts of TaXa21-H and TaWRKY53 while silencing of TaXa21-H and TaWRKY53 impairs wheat resistance to Pst. Overall, this represents formal evidence that plant pathogens produce 'splicing' effectors, which regulate host pre-mRNA splicing by direct engagement of the splicing sites, thereby interfering with host immunity.

Investigating the disordered regions (MoRFs, SLiMs and LCRs) and functions of mimicry proteins/peptides in silico

Microbial mimicry of the host proteins/peptides can elicit host auto-reactive T- or B-cells resulting in autoimmune disease(s). Since intrinsically disordered protein regions (IDPRs) are involved in several host cell signaling and PPI networks, molecular mimicry of the IDPRs can help the pathogens in substituting their own proteins in the host cell-signaling and PPI networks and, ultimately hijacking the host cellular machinery. Thus, the present study was conducted to discern the structural disorder and intrinsically disordered protein regions (IDPRs) like, molecular recognition features (MoRFs), short linear motifs (SLiMs), and low complexity regions (LCRs) in the experimentally verified mimicry proteins and peptides (mimitopes) of bacteria, viruses and host. Also, functional characteristics of the mimicry proteins were studied in silico. Our results indicated that 78% of the bacterial host mimicry proteins and 45% of the bacterial host mimitopes were moderately/highly disordered while, 73% of the viral host mimicry proteins and 31% of the viral host mimitopes were moderately/highly disordered. Among the pathogens, 27% of the bacterial mimicry proteins and 13% of the bacterial mimitopes were moderately/highly disordered while, 53% of the viral mimicry proteins and 21% of the viral mimitopes were moderately/highly disordered. Though IDPR were frequent in host, bacterial and viral mimicry proteins, only a few mimitopes overlapped with the IDPRs like, MoRFs, SLiMs and LCRs. This suggests that most of the microbes cannot use molecular mimicry to modulate the host PPIs and hijack the host cell machinery. Functional analyses indicated that most of the pathogens exhibited mimicry with the host proteins involved in ion binding and signaling pathways. This is the first report on the disordered regions and functional aspects of experimentally proven host and microbial mimicry proteins.

An improved repertoire of splicing variants and their potential roles in Arabidopsis photomorphogenic development

Background

Light switches on the photomorphogenic development of young plant seedlings, allowing young seedlings to acquire photosynthetic capacities and gain survival fitness. Light regulates gene expression at all levels of the central dogma, including alternative splicing (AS) during the photomorphogenic development. However, accurate determination of full-length (FL) splicing variants has been greatly hampered by short-read RNA sequencing technologies.

Result

In this study, we adopt PacBio isoform sequencing (Iso-seq) to overcome the limitation of the short-read RNA-seq technologies. Normalized cDNA libraries used for Iso-seq allows for comprehensive and effective identification of FL AS variants. Our analyses reveal more than 30,000 splicing variant models from approximately 16,500 gene loci and additionally identify approximately 700 previously unannotated genes. Among the variants, approximately 12,000 represent new gene models. Intron retention (IR) is the most frequently observed form of variants, and many IR-containing AS variants show evidence of engagement in translation. Our study reveals the formation of heterodimers of transcription factors composed of annotated and IR-containing AS variants. Moreover, transgenic plants overexpressing the IR forms of two B-BOX DOMAIN PROTEINs exhibits light-hypersensitive phenotypes, suggesting their regulatory roles in modulating optimal light responses.

Conclusions

This study provides an accurate and comprehensive portrait of full-length transcript isoforms and experimentally confirms the presence of de novo synthesized AS variants that impose regulatory functions in photomorphogenic development in Arabidopsis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02620-2.

How alternative splicing changes the properties of plant proteins

Most plant primary transcripts undergo alternative splicing (AS), and its impact on protein diversity is a subject of intensive investigation. Several studies have uncovered various mechanisms of how particular protein splice isoforms operate. However, the common principles behind the AS effects on protein function in plants have rarely been surveyed. Here, on the selected examples, we highlight diverse tissue expression patterns, subcellular localization, enzymatic activities, abilities to bind other molecules and other relevant features. We describe how the protein isoforms mutually interact to underline their intriguing roles in altering the functionality of protein complexes. Moreover, we also discuss the known cases when these interactions have been placed inside the autoregulatory loops. This review is particularly intended for plant cell and developmental biologists who would like to gain inspiration on how the splice variants encoded by their genes of interest may coordinately work.

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.

Mutually opposing activity of PIN7 splicing isoforms is required for auxin-mediated tropic responses in Arabidopsis thaliana

Advanced transcriptome sequencing has revealed that the majority of eukaryotic genes undergo alternative splicing (AS). Nonetheless, little effort has been dedicated to investigating the functional relevance of particular splicing events, even those in the key developmental and hormonal regulators. Combining approaches of genetics, biochemistry and advanced confocal microscopy, we describe the impact of alternative splicing on the PIN7 gene in the model plant Arabidopsis thaliana. PIN7 encodes a polarly localized transporter for the phytohormone auxin and produces two evolutionarily conserved transcripts, PIN7a and PIN7b. PIN7a and PIN7b, differing in a four amino acid stretch, exhibit almost identical expression patterns and subcellular localization. We reveal that they are closely associated and mutually influence each other's mobility within the plasma membrane. Phenotypic complementation tests indicate that the functional contribution of PIN7b per se is minor, but it markedly reduces the prominent PIN7a activity, which is required for correct seedling apical hook formation and auxin-mediated tropic responses. Our results establish alternative splicing of the PIN family as a conserved, functionally relevant mechanism, revealing an additional regulatory level of auxin-mediated plant development.

Pre-mRNA alternative splicing as a modulator for heat stress response in plants

The molecular responses of plants to the important abiotic stress, heat stress (HS), have been extensively studied at the transcriptional level. Alternative splicing (AS) is a post-transcriptional regulatory process in which an intron-containing gene can generate more than one mRNA variant. The impact of HS on the pre-mRNA splicing process has been reported in various eukaryotes but seldom discussed in-depth, especially in plants. Here, we review AS regulation in response to HS in different plant species. We discuss potential molecular mechanisms controlling heat-inducible AS regulation in plants and hypothesize that AS regulation participates in heat-priming establishment and HS memory maintenance. We propose that the pre-mRNA splicing variation is an important regulator of plant HS responses (HSRs).

Expression of a human cDNA in moss results in spliced mRNAs and fragmentary protein isoforms

Production of biopharmaceuticals relies on the expression of mammalian cDNAs in host organisms. Here we show that the expression of a human cDNA in the moss Physcomitrium patens generates the expected full-length and four additional transcripts due to unexpected splicing. This mRNA splicing results in non-functional protein isoforms, cellular misallocation of the proteins and low product yields. We integrated these results together with the results of our analysis of all 32,926 protein-encoding Physcomitrella genes and their 87,533 annotated transcripts in a web application, physCO, for automatized optimization. A thus optimized cDNA results in about twelve times more protein, which correctly localizes to the ER. An analysis of codon preferences of different production hosts suggests that similar effects occur also in non-plant hosts. We anticipate that the use of our methodology will prevent so far undetected mRNA heterosplicing resulting in maximized functional protein amounts for basic biology and biotechnology.

Alternative splicing in plant abiotic stress responses

Modifications of the cellular proteome pool upon stress allow plants to tolerate environmental changes. Alternative splicing is the most significant mechanism responsible for the production of multiple protein isoforms from a single gene. The spliceosome, a large ribonucleoprotein complex, together with several associated proteins, controls this pre-mRNA processing, adding an additional level of regulation to gene expression. Deep sequencing of transcriptomes revealed that this co- or post-transcriptional mechanism is highly induced by abiotic stress, and concerns vast numbers of stress-related genes. Confirming the importance of splicing in plant stress adaptation, key players of stress signaling have been shown to encode alternative transcripts, whereas mutants lacking splicing factors or associated components show a modified sensitivity and defective responses to abiotic stress. Here, we examine recent literature on alternative splicing and splicing alterations in response to environmental stresses, focusing on its role in stress adaptation and analyzing the future perspectives and directions for research.

Direct long-read RNA sequencing identifies a subset of questionable exitrons likely arising from reverse transcription artifacts

Resistance to CD19-directed immunotherapies in lymphoblastic leukemia has been attributed, among other factors, to several aberrant CD19 pre-mRNA splicing events, including recently reported excision of a cryptic intron embedded within CD19 exon 2. While "exitrons" are known to exist in hundreds of human transcripts, we discovered, using reporter assays and direct long-read RNA sequencing (dRNA-seq), that the CD19 exitron is an artifact of reverse transcription. Extending our analysis to publicly available datasets, we identified dozens of questionable exitrons, dubbed "falsitrons," that appear only in cDNA-seq, but never in dRNA-seq. Our results highlight the importance of dRNA-seq for transcript isoform validation.

Brief exposure to warm temperatures reduces intron retention in Kdm6b in a species with temperature-dependent sex determination

Animals with temperature-dependent sex determination (TSD) respond to thermal cues during early embryonic development to trigger gonadal differentiation. TSD has primarily been studied using constant temperature incubations, where embryos are exposed to constant male- or female-producing temperatures, and these studies have identified genes that display sex-specific expression in response to incubation temperature. Kdm6b, a histone demethylase gene, has received specific attention as it is among the initial genes to respond to incubation temperature and is necessary for testis development. Interestingly, Kdm6b retains an intron when eggs are incubated at a constant male-producing temperature, but the role of thermal variability in this developmental process is relatively understudied. Species with TSD regularly experience thermal cues that fluctuate between male- and female-producing temperatures throughout development but it is unclear how Kdm6b responds to such variable temperatures. In this study, we investigate temperature-sensitive splicing in Kdm6b by exposing embryos to male- and female-producing thermal conditions. We show a rapid decrease in levels of the intron retaining transcript of Kdm6b upon exposure to female-producing conditions. These results demonstrate that, under ecologically relevant conditions, temperature-sensitive splicing can differentially regulate genes critical to TSD.

Gene Architecture Facilitates Intron-Mediated Enhancement of Transcription

Introns impact several vital aspects of eukaryotic organisms like proteomic plasticity, genomic stability, stress response and gene expression. A role for introns in the regulation of gene expression at the level of transcription has been known for more than thirty years. The molecular basis underlying the phenomenon, however, is still not entirely clear. An important clue came from studies performed in budding yeast that indicate that the presence of an intron within a gene results in formation of a multi-looped gene architecture. When looping is defective, these interactions are abolished, and there is no enhancement of transcription despite normal splicing. In this review, we highlight several potential mechanisms through which looping interactions may enhance transcription. The promoter-5′ splice site interaction can facilitate initiation of transcription, the terminator-3′ splice site interaction can enable efficient termination of transcription, while the promoter-terminator interaction can enhance promoter directionality and expedite reinitiation of transcription. Like yeast, mammalian genes also exhibit an intragenic interaction of the promoter with the gene body, especially exons. Such promoter-exon interactions may be responsible for splicing-dependent transcriptional regulation. Thus, the splicing-facilitated changes in gene architecture may play a critical role in regulation of transcription in yeast as well as in higher eukaryotes.

A pan-cancer transcriptome analysis of exitron splicing identifies novel cancer driver genes and neoepitopes

Exitron splicing (EIS) creates a cryptic intron (called an exitron) within a protein-coding exon to increase proteome diversity. EIS is poorly characterized, but emerging evidence suggests a role for EIS in cancer. Through a systematic investigation of EIS across 33 cancers from 9,599 tumor transcriptomes, we discovered that EIS affected 63% of human coding genes and that 95% of those events were tumor specific. Notably, we observed a mutually exclusive pattern between EIS and somatic mutations in their affected genes. Functionally, we discovered that EIS altered known and novel cancer driver genes for causing gain- or loss-of-function, which promotes tumor progression. Importantly, we identified EIS-derived neoepitopes that bind to major histocompatibility complex (MHC) class I or II. Analysis of clinical data from a clear cell renal cell carcinoma cohort revealed an association between EIS-derived neoantigen load and checkpoint inhibitor response. Our findings establish the importance of considering EIS alterations when nominating cancer driver events and neoantigens.

An Intron of Invertebrate Microphthalmia Transcription Factor Gene Is Evolved from a Longer Ancestral Sequence

Introns are highly variable in number and size. Sequence simulation is an effective method to elucidate intron evolution patterns. Previously, we have reported that introns are more likely to evolve through mutation-and-deletion (MD) rather than through mutation-and-insertion (MI). In the present study, we further studied evolution models by allowing insertion in the MD model and by allowing deletion in the MI model at various frequencies. It was found that all deletion-biased models with proper parameter settings could generate sequences with attributes matchable to 16 invertebrate introns from the microphthalmia transcription factor gene, whereas all insertion-biased models with any parameter settings failed to generate such sequences. We conclude that the examined invertebrate introns may have evolved from a longer ancestral sequence in a deletion-biased pattern. The constructed models are useful for studying the evolution of introns from other genes and/or from other taxonomic groups. (C++ scripts of all deletion- and insertion-biased models are available upon request.)

Targeting alternative splicing by RNAi: from the differential impact on splice variants to triggering artificial pre-mRNA splicing

Alternative splicing generates multiple transcript and protein isoforms from a single gene and controls transcript intracellular localization and stability by coupling to mRNA export and nonsense-mediated mRNA decay (NMD). RNA interference (RNAi) is a potent mechanism to modulate gene expression. However, its interactions with alternative splicing are poorly understood. We used artificial microRNAs (amiRNAs, also termed shRNAmiR) to knockdown all splice variants of selected target genes in Arabidopsis thaliana. We found that splice variants, which vary by their protein-coding capacity, subcellular localization and sensitivity to NMD, are affected differentially by an amiRNA, although all of them contain the target site. Particular transcript isoforms escape amiRNA-mediated degradation due to their nuclear localization. The nuclear and NMD-sensitive isoforms mask RNAi action in alternatively spliced genes. Interestingly, Arabidopsis SPL genes, which undergo alternative splicing and are targets of miR156, are regulated in the same manner. Moreover, similar results were obtained in mammalian cells using siRNAs, indicating cross-kingdom conservation of these interactions among RNAi and splicing isoforms. Furthermore, we report that amiRNA can trigger artificial alternative splicing, thus expanding the RNAi functional repertoire. Our findings unveil novel interactions between different post-transcriptional processes in defining transcript fates and regulating gene expression.

Differential nucleosome occupancy modulates alternative splicing in Arabidopsis thaliana

Alternative splicing (AS) is a major gene regulatory mechanism in plants. Recent evidence supports co-transcriptional splicing in plants, hence the chromatin state can impact AS. However, how dynamic changes in the chromatin state such as nucleosome occupancy influence the cold-induced AS remains poorly understood. Here, we generated transcriptome (RNA-Seq) and nucleosome positioning (MNase-Seq) data for Arabidopsis thaliana to understand how nucleosome positioning modulates cold-induced AS. Our results show that characteristic nucleosome occupancy levels are strongly associated with the type and abundance of various AS events under normal and cold temperature conditions in Arabidopsis. Intriguingly, exitrons, alternatively spliced internal regions of protein-coding exons, exhibit distinctive nucleosome positioning pattern compared to other alternatively spliced regions. Likewise, nucleosome patterns differ between exitrons and retained introns, pointing to their distinct regulation. Collectively, our data show that characteristic changes in nucleosome positioning modulate AS in plants in response to cold.

Viral Perturbation of Alternative Splicing of a Host Transcript Benefits Infection

Pathogens disturb alternative splicing patterns of infected eukaryotic hosts. However, in plants it is unknown if this is incidental to infection or represents a pathogen-induced remodeling of host gene expression needed to support infection. Here, we compared changes in transcription and protein accumulation with changes in transcript splicing patterns in maize (Zea mays) infected with the globally important pathogen sugarcane mosaic virus (SCMV). Our results suggested that changes in alternative splicing play a major role in determining virus-induced proteomic changes. Focusing on maize phytoene synthase1 (ZmPSY1), which encodes the key regulatory enzyme in carotenoid biosynthesis, we found that although SCMV infection decreases total ZmPSY1 transcript accumulation, the proportion of splice variant T001 increases by later infection stages so that ZmPSY1 protein levels are maintained. We determined that ZmPSY1 has two leaf-specific transcripts, T001 and T003, distinguished by differences between the respective 3'-untranslated regions (UTRs). The shorter 3'-UTR of T001 makes it the more efficient mRNA. Nonsense ZmPSY1 mutants or virus-induced silencing of ZmPSY1 expression suppressed SCMV accumulation, attenuated symptoms, and decreased chloroplast damage. Thus, ZmPSY1 acts as a proviral host factor that is required for virus accumulation and pathogenesis. Taken together, our findings reveal that SCMV infection-modulated alternative splicing ensures that ZmPSY1 synthesis is sustained during infection, which supports efficient virus infection.

Intron retention and its impact on gene expression and protein diversity: A review and a practical guide

Intron retention (IR) occurs when a complete and unspliced intron remains in mature mRNA. An increasing body of literature has demonstrated a major role for IR in numerous biological functions, including several that impact human health and disease. Although experimental technologies used to study other forms of mRNA splicing can also be used to investigate IR, a specialized downstream computational analysis is optimal for IR discovery and analysis. Here we provide a review of IR and its biological implications, as well as a practical guide for how to detect and analyze it. Several methods, including long read third generation direct RNA sequencing, are described. We have developed an R package, FakIR, to facilitate the execution of the bioinformatic tasks recommended in this review and a tutorial on how to fit them to users aims. Additionally, we provide guidelines and experimental protocols to validate IR discovery and to evaluate the potential impact of IR on gene expression and protein output. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Splicing Regulation/Alternative Splicing RNA Methods > RNA Analyses in vitro and In Silico.

Holding on to Junk Bonds: Intron Retention in Cancer and Therapy

Intron retention (IR) in cancer was for a long time overlooked by the scientific community, as it was previously considered to be an artifact of a dysfunctional spliceosome. Technological advancements made in the last decade offer unique opportunities to explore the role of IR as a widespread phenomenon that contributes to the transcriptional diversity of many cancers. Numerous studies in cancer have shed light on dysregulation of cellular mechanisms that lead to aberrant and pathologic IR. IR is not merely a mechanism of gene regulation, but rather it can mediate cancer pathogenesis and therapeutic resistance in various human diseases. The burden of IR in cancer is governed by perturbations to mechanisms known to regulate this phenomenon and include epigenetic variation, mutations within the gene body, and splicing factor dysregulation. This review summarizes possible causes for aberrant IR and discusses the role of IR in therapy or as a consequence of disease treatment. As neoepitopes originating from retained introns can be presented on the cancer cell surface, the development of personalized cancer vaccines based on IR-derived neoepitopes should be considered. Ultimately, a deeper comprehension about the origins and consequences of aberrant IR may aid in the development of such personalized cancer vaccines.

Analysis of Plasmodium vivax schizont transcriptomes from field isolates reveals heterogeneity of expression of genes involved in host-parasite interactions

Plasmodium vivax gene regulation remains difficult to study due to the lack of a robust in vitro culture method, low parasite densities in peripheral circulation and asynchronous parasite development. We adapted an RNA-seq protocol “DAFT-seq” to sequence the transcriptome of four P. vivax field isolates that were cultured for a short period ex vivo before using a density gradient for schizont enrichment. Transcription was detected from 78% of the PvP01 reference genome, despite being schizont-enriched samples. This extensive data was used to define thousands of 5′ and 3′ untranslated regions, some of which overlapped with neighbouring transcripts, and to improve the gene models of 352 genes, including identifying 20 novel gene transcripts. This dataset has also significantly increased the known amount of heterogeneity between P. vivax schizont transcriptomes from individual patients. The majority of genes found to be differentially expressed between the isolates lack Plasmodium falciparum homologs and are predicted to be involved in host-parasite interactions, with an enrichment in reticulocyte binding proteins, merozoite surface proteins and exported proteins with unknown function. An improved understanding of the diversity within P. vivax transcriptomes will be essential for the prioritisation of novel vaccine targets.

Phytophthora Effectors Modulate Genome-wide Alternative Splicing of Host mRNAs to Reprogram Plant Immunity

Alternative splicing (AS) of pre-mRNAs increases transcriptome and proteome diversity, regulates gene expression through multiple mechanisms, and plays important roles in plant development and stress responses. However, the prevalence of genome-wide plant AS changes during infection and the mechanisms by which pathogens modulate AS remain poorly understood. Here, we examined the global AS changes in tomato leaves infected with Phytophthora infestans, the infamous Irish famine pathogen. We show that more than 2000 genes exhibiting significant changes in AS are not differentially expressed, indicating that AS is a distinct layer of transcriptome reprogramming during plant-pathogen interactions. Furthermore, our results show that P. infestans subverts host immunity by repressing the AS of positive regulators of plant immunity and promoting the AS of susceptibility factors. To study the underlying mechanism, we established a luminescence-based AS reporter system in Nicotiana benthamiana to screen pathogen effectors modulating plant AS. We identified nine splicing regulatory effectors (SREs) from 87 P. infestans effectors. Further studies revealed that SRE3 physically binds U1-70K to manipulate the plant AS machinery and subsequently modulates AS-mediated plant immunity. Our study not only unveils genome-wide plant AS reprogramming during infection but also establishes a novel AS screening tool to identify SREs from a wide range of plant pathogens, providing opportunities to understand the splicing regulatory mechanisms through which pathogens subvert plant immunity.

Alternative Splicing During the Chlamydomonas reinhardtii Cell Cycle

Genome-wide analysis of transcriptome data in Chlamydomonas reinhardtii shows periodic patterns in gene expression levels when cultures are grown under alternating light and dark cycles so that G1 of the cell cycle occurs in the light phase and S/M/G0 occurs during the dark phase. However, alternative splicing, a process that enables a greater protein diversity from a limited set of genes, remains largely unexplored by previous transcriptome based studies in C. reinhardtii. In this study, we used existing longitudinal RNA-seq data obtained during the light-dark cycle to investigate the changes in the alternative splicing pattern and found that 3277 genes (19.75% of 17,746 genes) undergo alternative splicing. These splicing events include Alternative 5′ (Alt 5′), Alternative 3′ (Alt 3′) and Exon skipping (ES) events that are referred as alternative site selection (ASS) events and Intron retention (IR) events. By clustering analysis, we identified a subset of events (26 ASS events and 10 IR events) that show periodic changes in the splicing pattern during the cell cycle. About two-thirds of these 36 genes either introduce a pre-termination codon (PTC) or introduce insertions or deletions into functional domains of the proteins, which implicate splicing in altering gene function. These findings suggest that alternative splicing is also regulated during the Chlamydomonas cell cycle, although not as extensively as changes in gene expression. The longitudinal changes in the alternative splicing pattern during the cell cycle captured by this study provides an important resource to investigate alternative splicing in genes of interest during the cell cycle in Chlamydomonas reinhardtii and other eukaryotes.