Microglia specific alternative splicing alterations in multiple sclerosis.

Several aberrant alternative splicing (AS) events and their regulatory mechanisms are widely recognized in multiple sclerosis (MS). Yet the cell-type specific AS events have not been extensively examined. Here we assessed the diversity of AS events using web-based RNA-seq data of sorted CD15-CD11b+ microglia in white matter (WM) region from 10 patients with MS and 11 control subjects. The GSE111972 dataset was downloaded from GEO and ENA databases, aligned to the GRCh38 reference genome from ENSEMBL via STAR. rMATS was used to assess five types of AS events, alternative 3'SS (A3SS), alternative 5'SS (A5SS), skipped exon (SE), retained intron (RI) and mutually exclusive exons (MXE), followed by visualizing with rmats2sashimiplot and maser. Differential genes or transcripts were analyzed using the limma R package. Gene ontology (GO) analysis was performed with the clusterProfiler R package. 42,663 raw counts of AS events were identified and 132 significant AS events were retained based on the filtered criteria: 1) average coverage >10 and 2) delta percent spliced in (ΔPSI) >0.1. SE was the most common AS event (36.36%), followed by MXE events (32.58%), and RI (18.94%). Genes related to telomere maintenance and organization primarily underwent SE splicing, while genes associated with protein folding and mitochondrion organization were predominantly spliced in the MXE pattern. Conversely, genes experiencing RI were enriched in immune response and immunoglobulin production. In conclusion, we identified microglia-specific AS changes in the white matter of MS patients, which may shed light on novel pathological mechanisms underlying MS.

Single-cell long-read sequencing in human cerebral organoids uncovers cell-type-specific and autism-associated exons.

Dysregulation of alternative splicing has been repeatedly associated with neurodevelopmental disorders, but the extent of cell-type-specific splicing in human neural development remains largely uncharted. Here, single-cell long-read sequencing in induced pluripotent stem cell (iPSC)-derived cerebral organoids identifies over 31,000 uncatalogued isoforms and 4,531 cell-type-specific splicing events. Long reads uncover coordinated splicing and cell-type-specific intron retention events, which are challenging to study with short reads. Retained neuronal introns are enriched in RNA splicing regulators, showing shorter lengths, higher GC contents, and weaker 5' splice sites. We use this dataset to explore the biological processes underlying neurological disorders, focusing on autism. In comparison with prior transcriptomic data, we find that the splicing program in autistic brains is closer to the progenitor state than differentiated neurons. Furthermore, cell-type-specific exons harbor significantly more de novo mutations in autism probands than in siblings. Overall, these results highlight the importance of cell-type-specific splicing in autism and neuronal gene regulation.

Alzheimer's disease pathogenetic progression is associated with changes in regulated retained introns and editing of circular RNAs.

Introduction: The molecular changes leading to Alzheimer's disease (AD) progression are poorly understood. A decisive factor in the disease occurs when neurofibrillary tangles (NFT) composed of microtubule associated protein tau (MAPT) form in the entorhinal cortex and then spread throughout the brain. Methods: We therefore determined mRNA and circular RNA changes during AD progression, comparing Braak NFT stages I-VI. Total RNA was isolated from human brain (entorhinal and frontotemporal cortex). Poly(A)+ RNA was subjected to Nanopore sequencing, and total RNA was analyzed by standard Illumina sequencing. Circular RNAs were sequenced from RNase R treated and rRNA depleted total RNA. The sequences were analyzed using different bioinformatic tools, and expression constructs for circRNAs were analyzed in transfection experiments. Results: We detected 11,873 circRNAs of which 276 correlated with Braak NFT stages. Adenosine to inosine RNA editing increased about threefold in circRNAs during AD progression. Importantly, this correlation cannot be detected with mRNAs. CircMAN2A1 expression correlated with AD progression and transfection experiments indicated that RNA editing promoted its translation using start codons out of frame with linear mRNAs, which generates novel proteins. Discussion: Thus, we identified novel regulated retained introns that correlate with NFT Braak stages and provide evidence for a role of translated circRNAs in AD development.

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.

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.

Distinct microglia alternative splicing in Alzheimer's disease.

Numerous alternative splicing (AS) events have been documented in Alzheimer's disease (AD). However, cell type-specific AS analysis is still lacking. We described AS events in the hippocampal microglia sorted by CD45 and CD11b from Aß precursor protein (APP) and non-transgenic (Ntg) mice. GSE171195 dataset was downloaded from GEO database, aligned to GRCm39 genome. Skipped exon (SE), alternative 3'SS (A3SS), retained intron (RI), alternative 5'SS (A5SS), and mutually exclusive exons (MXE) were evaluated using rMATS and maser. Differential expressed genes or transcripts were analyzed via limma. Gene ontology and correlation analyses were performed with clusterProfiler and ggcorrplot R packages. 36,340 raw counts of AS were identified, and 95 significant AS events were eventually selected with strict criteria: (1) average coverage >5; (2) delta percent spliced in >0.1. SE was the most common AS events (68.42%), followed by A3SS and RI. Autophagy genes were mainly spliced in SE events, actin depolymerization genes spliced in A3SS events, while synaptic plasticity related genes were mainly spliced in RI pattern. These significant AS events may be regulated by dysregulated splicing factors in AD. In conclusion, we revealed microglia specific AS events in AD, and our study provides novel pathological mechanisms in the pathogenesis of AD.

Ontology Specific Alternative Splicing Changes in Alzheimer's Disease.

Alternative splicing (AS) is a common phenomenon and correlates with aging and aging-related disorders including Alzheimer's disease (AD). We aimed to systematically characterize AS changes in the cerebral cortex of 9-month-old APP/PS1 mice. The GSE132177 dataset was downloaded from GEO and ENA databases, aligned to the GRCm39 reference genome from ENSEMBL via STAR. Alternative 3'SS (A3SS), alternative 5'SS (A5SS), skipped exon (SE), retained intron (RI), and mutually exclusive exons (MXE) AS events were evaluated using rMATS, rmats2sashimiplot, and maser. Differential genes or transcripts were analyzed using the limma R package. Gene ontology analysis was performed with the clusterProfiler R package. A total of 60,705 raw counts of AS were identified, and 113 significant AS events were finally selected in accordance with the selection criteria: 1) average coverage >10 and 2) delta percent spliced in (ΔPSI) >0.1. SE was the most abundant AS event (61.95%), and RI was the second most abundant AS type (13.27%), followed by A3SS (12.39%), thereafter A5SS and MXE comprised of 12.39%. Interestingly, genes that experienced SE were enriched in histone acetyltransferase (HAT) complex, while genes spliced by RI were enriched in autophagy and those which experienced A3SS were enriched in methyltransferase activity revealed by GO analysis. In conclusion, we revealed ontology specific AS changes in AD. Our analysis provides novel pathological mechanisms of AD.

Intron retention coupled with nonsense-mediated decay is involved in cellulase biosynthesis in cellulolytic fungi.

BACKGROUND: Knowledge on regulatory networks associated with cellulase biosynthesis is prerequisite for exploitation of such regulatory systems in enhancing cellulase production with low cost. The biological functions of intron retention (IR) and nonsense-mediated mRNA decay (NMD) in filamentous fungi is lack of study, let alone their roles in cellulase biosynthesis. RESULTS: We found that major cellulase genes (cel7a, cel7b, and cel3a) exhibited concomitant decrease in IR rates and increase in their gene expression in T. reesei under cellulase-producing condition (cellulose and lactose) that was accompanied with a more active NMD pathway, as compared to cellulase non-producing condition (glucose). In the presence of the NMD pathway inhibitor that successfully repressed the NMD pathway, the mRNA levels of cellulase genes were sharply down-regulated, but the rates of IR in these genes were significantly up-regulated. Consistently, the cellulase activities were severely inhibited. In addition, the NMD pathway inhibitor caused the downregulated mRNA levels of two important genes of the target of rapamycin (TOR) pathway, trfkbp12 and trTOR1. The absence of gene trfkbp12 made the cellulase production in T. reesei more sensitive to the NMD pathway inhibitor. CONCLUSIONS: All these findings suggest that the IR of cellulase genes regulates their own gene expression by coupling with the NMD pathway, which might involve the TOR pathway. Our results provide better understanding on intron retention, the NMD pathway, and cellulase production mechanism in filamentous fungi.

RNA-Seq Analysis Reveals Localization-Associated Alternative Splicing across 13 Cell Lines.

Alternative splicing, a ubiquitous phenomenon in eukaryotes, is a regulatory mechanism for the biological diversity of individual genes. Most studies have focused on the effects of alternative splicing for protein synthesis. However, the transcriptome-wide influence of alternative splicing on RNA subcellular localization has rarely been studied. By analyzing RNA-seq data obtained from subcellular fractions across 13 human cell lines, we identified 8720 switching genes between the cytoplasm and the nucleus. Consistent with previous reports, intron retention was observed to be enriched in the nuclear transcript variants. Interestingly, we found that short and structurally stable introns were positively correlated with nuclear localization. Motif analysis reveals that fourteen RNA-binding protein (RBPs) are prone to be preferentially bound with such introns. To our knowledge, this is the first transcriptome-wide study to analyze and evaluate the effect of alternative splicing on RNA subcellular localization. Our findings reveal that alternative splicing plays a promising role in regulating RNA subcellular localization.

Molecular mechanisms regulating O-linked N-acetylglucosamine (O-GlcNAc)-processing enzymes.

The post-translational modification of proteins by O-linked N-acetylglucosamine (O-GlcNAc) dynamically programmes cellular physiology to maintain homoeostasis and tailor biochemical pathways to meet context-dependent cellular needs. Despite diverse roles of played by O-GlcNAc, only two enzymes act antagonistically to govern its cycling; O-GlcNAc transferase installs the monosaccharide on target proteins, and O-GlcNAc hydrolase removes it. The recent literature has exposed a network of mechanisms regulating these two enzymes to choreograph global, and target-specific, O-GlcNAc cycling in response to cellular stress and nutrient availability. Herein, we amalgamate these emerging mechanisms from a structural and molecular perspective to explore how the cell exerts fine control to regulate O-GlcNAcylation of diverse proteins in a selective fashion.

The nuclear-localized PPR protein OsNPPR1 is important for mitochondrial function and endosperm development in rice.

Pentatricopeptide repeat (PPR) proteins constitute one of the largest protein families in land plants. Recent studies revealed the functions of PPR proteins in organellar RNA metabolism and plant development, but the functions of most PPR proteins, especially PPRs localized in the nucleus, remain largely unknown. Here, we report the isolation and characterization of a rice mutant named floury and growth retardation1 (fgr1). fgr1 showed floury endosperm with loosely arranged starch grains, decreased starch and amylose contents, and retarded seedling growth. Map-based cloning showed that the mutant phenotype was caused by a single nucleotide substitution in the coding region of Os08g0290000. This gene encodes a nuclear-localized PPR protein, which we named OsNPPR1, that affected mitochondrial function. In vitro SELEX and RNA-EMSAs showed that OsNPPR1 was an RNA protein that bound to the CUCAC motif. Moreover, a number of retained intron (RI) events were detected in fgr1. Thus, OsNPPR1 was involved in regulation of mitochondrial development and/or functions that are important for endosperm development. Our results provide novel insights into coordinated interaction between nuclear-localized PPR proteins and mitochondrial function.

Gene expression patterns of chicken neuregulin 3 in association with copy number variation and frameshift deletion.

BACKGROUND: Neuregulin 3 (NRG3) plays a key role in central nervous system development and is a strong candidate for human mental disorders. Thus, genetic variation in NRG3 may have some impact on a variety of phenotypes in non-mammalian vertebrates. Recently, genome-wide screening for short insertions and deletions in chicken (Gallus gallus) genomes has provided useful information about structural variation in functionally important genes. NRG3 is one such gene that has a putative frameshift deletion in exon 2, resulting in premature termination of translation. Our aims were to characterize the structure of chicken NRG3 and to compare expression patterns between NRG3 isoforms. RESULTS: Depending on the presence or absence of the 2-bp deletion in chicken NRG3, 3 breeds (red junglefowl [RJF], Boris Brown [BB], and Hinai-jidori [HJ]) were genotyped using flanking primers. In the commercial breeds (BB and HJ), approximately 45% of individuals had at least one exon 2 allele with the 2-bp deletion, whereas there was no deletion allele in RJF. The lack of a homozygous mutant indicated the existence of duplicated NRG3 segments in the chicken genome. Indeed, highly conserved elements consisting of exon 1, intron 1, exon 2, and part of intron 2 were found in the reference RJF genome, and quantitative PCR detected copy number variation (CNV) between breeds as well as between individuals. The copy number of conserved elements was significantly higher in chicks harboring the 2-bp deletion in exon 2. We identified 7 novel transcript variants using total mRNA isolated from the amygdala. Novel isoforms were found to lack the exon 2 cassette, which probably harbored the premature termination codon. The relative transcription levels of the newly identified isoforms were almost the same between chick groups with and without the 2-bp deletion, while chicks with the deletion showed significant suppression of the expression of previously reported isoforms. CONCLUSIONS: A putative frameshift deletion and CNV in chicken NRG3 are structural mutations that occurred before the establishment of commercial chicken lines. Our results further suggest that the putative frameshift deletion in exon 2 may potentially affect the expression level of particular isoforms of chicken NRG3.

Differential intron retention in Jumonji chromatin modifier genes is implicated in reptile temperature-dependent sex determination.

In many vertebrates, sex of offspring is determined by external environmental cues rather than by sex chromosomes. In reptiles, for instance, temperature-dependent sex determination (TSD) is common. Despite decades of work, the mechanism by which temperature is converted into a sex-determining signal remains mysterious. This is partly because it is difficult to distinguish the primary molecular events of TSD from the confounding downstream signatures of sexual differentiation. We use the Australian central bearded dragon, in which chromosomal sex determination is overridden at high temperatures to produce sex-reversed female offspring, as a unique model to identify TSD-specific features of the transcriptome. We show that an intron is retained in mature transcripts from each of two Jumonji family genes, JARID2 and JMJD3, in female dragons that have been sex-reversed by temperature but not in normal chromosomal females or males. JARID2 is a component of the master chromatin modifier Polycomb Repressive Complex 2, and the mammalian sex-determining factor SRY is directly regulated by an independent but closely related Jumonji family member. We propose that the perturbation of JARID2/JMJD3 function by intron retention alters the epigenetic landscape to override chromosomal sex-determining cues, triggering sex reversal at extreme temperatures. Sex reversal may then facilitate a transition from genetic sex determination to TSD, with JARID2/JMJD3 intron retention preserved as the decisive regulatory signal. Significantly, we also observe sex-associated differential retention of the equivalent introns in JARID2/JMJD3 transcripts expressed in embryonic gonads from TSD alligators and turtles, indicative of a reptile-wide mechanism controlling TSD.

In silico analysis of the sequence features responsible for alternatively spliced introns in the model green alga Chlamydomonas reinhardtii.

Alternatively spliced introns are the ones that are usually spliced but can be occasionally retained in a transcript isoform. They are the most frequently used alternative splice form in plants (~50% of alternative splicing events). Chlamydomonas reinhardtii, a unicellular alga, is a good model to understand alternative splicing (AS) in plants from an evolutionary perspective as it diverged from land plants a billion years ago. Using over 7 million cDNA sequences from both pyrosequencing and Sanger sequencing, we found that a much higher percentage of genes (~20% of multi-exon genes) undergo AS than previously reported (3-5%). We found a full component of SR and SR-like proteins possibly involved in AS. The most prevalent type of AS event (40%) was retention of introns, most of which were supported by multiple cDNA evidence (72%) while only 20% of them have coding capacity. By comparing retained and constitutive introns, we identified sequence features potentially responsible for the retention of introns, in the framework of an "intron definition" model for splicing. We find that retained introns tend to have a weaker 5' splice site, more Gs in their poly-pyrimidine tract and a lesser conservation of nucleotide 'C' at position -3 of the 3' splice site. In addition, the sequence motifs found in the potential branch-point region differed between retained and constitutive introns. Furthermore, the enrichment of G-triplets and C-triplets among the first and last 50 nt of the introns significantly differ between constitutive and retained introns. These could serve as intronic splicing enhancers. All the alternative splice forms can be accessed at http://bioinfolab.miamioh.edu/cgi-bin/PASA_r20140417/cgi-bin/status_report.cgi?db=Chre_AS .

Preferential binding of a stable G3BP ribonucleoprotein complex to intron-retaining transcripts in mouse brain and modulation of their expression in the cerebellum.

Neuronal granules play an important role in the localization and transport of translationally silenced messenger ribonucleoproteins in neurons. Among the factors associated with these granules, the RNA-binding protein G3BP1 (stress-granules assembly factor) is involved in neuronal plasticity and is induced in Alzheimer's disease. We immunopurified a stable complex containing G3BP1 from mouse brain and performed high-throughput sequencing and cross-linking immunoprecipitation to identify the associated RNAs. The G3BP-complex contained the deubiquitinating protease USP10, CtBP1 and the RNA-binding proteins Caprin-1, G3BP2a and splicing factor proline and glutamine rich, or PSF. The G3BP-complex binds preferentially to transcripts that retain introns, and to non-coding sequences like 3'-untranslated region and long non-coding RNAs. Specific transcripts with retained introns appear to be enriched in the cerebellum compared to the rest of the brain and G3BP1 depletion decreased this intron retention in the cerebellum of G3BP1 knockout mice. Among the enriched transcripts, we found an overrepresentation of genes involved in synaptic transmission, especially glutamate-related neuronal transmission. Notably, G3BP1 seems to repress the expression of the mature Grm5 (metabotropic glutamate receptor 5) transcript, by promoting the retention of an intron in the immature transcript in the cerebellum. Our results suggest that G3BP is involved in a new functional mechanism to regulate non-coding RNAs including intron-retaining transcripts, and thus have broad implications for neuronal gene regulation, where intron retention is widespread.

Fungal alternative splicing is associated with multicellular complexity and virulence: a genome-wide multi-species study.

Alternative splicing (AS) is a cellular process that increases a cell's coding capacity from a limited set of genes. Although AS is common in higher plants and animals, its prevalence in other eukaryotes is mostly unknown. In fungi the involvement of AS in gene expression and its effect on multi-cellularity and virulence is of great medical and economic interest. We present a genome-wide comparative study of AS in 23 informative fungi of different taxa, based on alignments of public transcript sequences. Random sampling of expressed sequence tags allows for robust and comparable estimations of AS rates. We find that a greater fraction of fungal genes than previously expected is associated with AS. We estimate that on average, 6.4% of the annotated genes are affected by AS, with Cryptococcus neoformans showing an extraordinary rate of 18%. The investigated Basidiomycota show higher average AS rates (8.6%) than the Ascomycota (6.0%), although not significant. We find that multi-cellular complexity and younger evolutionary age associate with higher AS rates. Furthermore, AS affects genes involved in pathogenic lifestyle, particularly in functions of stress response and dimorphic switching. Together, our analysis strongly supports the view that AS is a rather common phenomenon in fungi and associates with higher multi-cellular complexity.