Splicing Analysis Research Articles

RNA splicing analysis using heterogeneous and large RNA-seq datasets

The ubiquity of RNA-seq has led to many methods that use RNA-seq data to analyze variations in RNA splicing. However, available methods are not well suited for handling heterogeneous and large datasets. Such datasets scale to thousands of samples across dozens of experimental conditions, exhibit increased variability compared to biological replicates, and involve thousands of unannotated splice variants resulting in increased transcriptome complexity. We describe here a suite of algorithms and tools implemented in the MAJIQ v2 package to address challenges in detection, quantification, and visualization of splicing variations from such datasets. Using both large scale synthetic data and GTEx v8 as benchmark datasets, we assess the advantages of MAJIQ v2 compared to existing methods. We then apply MAJIQ v2 package to analyze differential splicing across 2,335 samples from 13 brain subregions, demonstrating its ability to offer insights into brain subregion-specific splicing regulation.

Alternative splicing in shaping the molecular landscape of the cochlea.

The cochlea is a complex organ comprising diverse cell types with highly specialized morphology and function. Until now, the molecular underpinnings of its specializations have mostly been studied from a transcriptional perspective, but accumulating evidence points to post-transcriptional regulation as a major source of molecular diversity. Alternative splicing is one of the most prevalent and well-characterized post-transcriptional regulatory mechanisms. Many molecules important for hearing, such as cadherin 23 or harmonin, undergo alternative splicing to produce functionally distinct isoforms. Some isoforms are expressed specifically in the cochlea, while some show differential expression across the various cochlear cell types and anatomical regions. Clinical phenotypes that arise from mutations affecting specific splice variants testify to the functional relevance of these isoforms. All these clues point to an essential role for alternative splicing in shaping the unique molecular landscape of the cochlea. Although the regulatory mechanisms controlling alternative splicing in the cochlea are poorly characterized, there are animal models with defective splicing regulators that demonstrate the importance of RNA-binding proteins in maintaining cochlear function and cell survival. Recent technological breakthroughs offer exciting prospects for overcoming some of the long-standing hurdles that have complicated the analysis of alternative splicing in the cochlea. Efforts toward this end will help clarify how the remarkable diversity of the cochlear transcriptome is both established and maintained.

SpliceTools, a suite of downstream RNA splicing analysis tools to investigate mechanisms and impact of alternative splicing.

As a fundamental aspect of normal cell signaling and disease states, there is great interest in determining alternative splicing (AS) changes in physiologic, pathologic, and pharmacologic settings. High throughput RNA sequencing and specialized software to detect AS has greatly enhanced our ability to determine transcriptome-wide splicing changes. Despite the richness of this data, deriving meaning from sometimes thousands of AS events is a substantial bottleneck for most investigators. We present SpliceTools, a suite of data processing modules that arms investigators with the ability to quickly produce summary statistics, mechanistic insights, and functional significance of AS changes through command line or through an online user interface. Utilizing RNA-seq datasets for 186 RNA binding protein knockdowns, nonsense mediated RNA decay inhibition, and pharmacologic splicing inhibition, we illustrate the utility of SpliceTools to distinguish splicing disruption from regulated transcript isoform changes, we show the broad transcriptome footprint of the pharmacologic splicing inhibitor, indisulam, we illustrate the utility in uncovering mechanistic underpinnings of splicing inhibition, we identify predicted neo-epitopes in pharmacologic splicing inhibition, and we show the impact of splicing alterations induced by indisulam on cell cycle progression. Together, SpliceTools puts rapid and easy downstream analysis at the fingertips of any investigator studying AS.

Detection and targeting of splicing deregulation in pediatric acute myeloid leukemia stem cells

Pediatric acute myeloid leukemia (pAML) is typified by high relapse rates and a relative paucity of somatic DNA mutations. Although seminal studies show that splicing factor mutations and mis-splicing fuel therapy-resistant leukemia stem cell (LSC) generation in adults, splicing deregulation has not been extensively studied in pAML. Herein, we describe single-cell proteogenomics analyses, transcriptome-wide analyses of FACS-purified hematopoietic stem and progenitor cells followed by differential splicing analyses, dual-fluorescence lentiviral splicing reporter assays, and the potential of a selective splicing modulator, Rebecsinib, in pAML. Using these methods, we discover transcriptomic splicing deregulation typified by differential exon usage. In addition, we discover downregulation of splicing regulator RBFOX2 and CD47 splice isoform upregulation. Importantly, splicing deregulation in pAML induces a therapeutic vulnerability to Rebecsinib in survival, self-renewal, and lentiviral splicing reporter assays. Taken together, the detection and targeting of splicing deregulation represent a potentially clinically tractable strategy for pAML therapy.

Farnesyl diphosphate synthase promotes cell proliferation by regulating gene expression and alternative splicing profiles in HeLa cells.

Farnesyl diphosphate synthase (FDPS), an essential enzyme involved in the mevalonate pathway, is implicated in various diseases, including multiple types of cancer. As an RNA-binding protein (RBP), FDPS is also involved in transcriptional and post-transcriptional regulation. However, to the best of our knowledge, transcriptome-wide targets of FDPS still remain unknown. In the present study, FDPS expression patterns in pan-cancer were analyzed. In addition, it was investigated how FDPS overexpression (FDPS-OE) regulates the transcriptome in HeLa cells. FDPS-OE increased the proliferation rate in HeLa cells by MTT assay. Using transcriptome-wide high throughput sequencing and bioinformatics analysis, it was found that FDPS upregulated the expression levels of genes enriched in cell proliferation and extracellular matrix organization, including the laminin subunit γ2, interferon-induced proteins with tetratricopeptide repeats 2 and matrix metallopeptidase 19 genes. According to alternative splicing (AS) analysis, FDPS modulated the splicing patterns of the bone morphogenic protein 1, semaphorin 4D, annexin A2 and sirtuin 2 genes, which are enriched in the cell cycle and DNA repair, and are related to cell proliferation. To corroborate the FDPS-regulated transcriptome findings, FDPS was overexpressed in human osteosarcoma cells. Differentially expressed genes and regulated AS genes in the cells were both validated by reverse transcription-quantitative PCR. The results suggested that, as an emerging RBP, FDPS may serve an important role in transcriptome profiles by altering gene expression and regulating AS. FDPS also affected the cell proliferation rate. These findings broaden the understanding of the molecular functions of FDPS, and the potential of FDPS as a target in therapy should be investigated.

Tip60's Novel RNA-Binding Function Modulates Alternative Splicing of Pre-mRNA Targets Implicated in Alzheimer's Disease.

The severity of Alzheimer's disease (AD) progression involves a complex interplay of genetics, age, and environmental factors orchestrated by histone acetyltransferase (HAT)-mediated neuroepigenetic mechanisms. While disruption of Tip60 HAT action in neural gene control is implicated in AD, alternative mechanisms underlying Tip60 function remain unexplored. Here, we report a novel RNA binding function for Tip60 in addition to its HAT function. We show that Tip60 preferentially interacts with pre-mRNAs emanating from its chromatin neural gene targets in the Drosophila brain and this RNA binding function is conserved in human hippocampus and disrupted in Drosophila brains that model AD pathology and in AD patient hippocampus of either sex. Since RNA splicing occurs co-transcriptionally and alternative splicing (AS) defects are implicated in AD, we investigated whether Tip60-RNA targeting modulates splicing decisions and whether this function is altered in AD. Replicate multivariate analysis of transcript splicing (rMATS) analysis of RNA-Seq datasets from wild-type and AD fly brains revealed a multitude of mammalian-like AS defects. Strikingly, over half of these altered RNAs are identified as bona-fide Tip60-RNA targets that are enriched for in the AD-gene curated database, with some of these AS alterations prevented against by increasing Tip60 in the fly brain. Further, human orthologs of several Tip60-modulated splicing genes in Drosophila are well characterized aberrantly spliced genes in human AD brains, implicating disruption of Tip60's splicing function in AD pathogenesis. Our results support a novel RNA interaction and splicing regulatory function for Tip60 that may underly AS impairments that hallmark AD etiology.SIGNIFICANCE STATEMENT Alzheimer's disease (AD) has recently emerged as a hotbed for RNA alternative splicing (AS) defects that alter protein function in the brain yet causes remain unclear. Although recent findings suggest convergence of epigenetics with co-transcriptional AS, whether epigenetic dysregulation in AD pathology underlies AS defects remains unknown. Here, we identify a novel RNA interaction and splicing regulatory function for Tip60 histone acetyltransferase (HAT) that is disrupted in Drosophila brains modeling AD pathology and in human AD hippocampus. Importantly, mammalian orthologs of several Tip60-modulated splicing genes in Drosophila are well characterized aberrantly spliced genes in human AD brain. We propose that Tip60-mediated AS modulation is a conserved critical posttranscriptional step that may underlie AS defects now characterized as hallmarks of AD.

Alternative Splicing Changes in Human Aging Brain Samples

Alternative splicing affects differential isoform generation from genes and may also affect the protein products cellular function. Targeted splicing analysis of exon array analysis of young compared to middle and old post mortem brain samples from 1,234 samples (a total of 10 brain regions) aged from 16 to over 100 allowed me to detect significant alternative splicing changes upon human aging. This kind of data may allow us to better understand the underlying molecular mechanisms.

Revealing the History and Mystery of RNA-Seq.

Advances in RNA-sequencing technologies have led to the development of intriguing experimental setups, a massive accumulation of data, and high demand for tools to analyze it. To answer this demand, computational scientists have developed a myriad of data analysis pipelines, but it is less often considered what the most appropriate one is. The RNA-sequencing data analysis pipeline can be divided into three major parts: data pre-processing, followed by the main and downstream analyses. Here, we present an overview of the tools used in both the bulk RNA-seq and at the single-cell level, with a particular focus on alternative splicing and active RNA synthesis analysis. A crucial part of data pre-processing is quality control, which defines the necessity of the next steps; adapter removal, trimming, and filtering. After pre-processing, the data are finally analyzed using a variety of tools: differential gene expression, alternative splicing, and assessment of active synthesis, the latter requiring dedicated sample preparation. In brief, we describe the commonly used tools in the sample preparation and analysis of RNA-seq data.

One pipeline to predict them all? On the prediction of alternative splicing from RNA-Seq data

RNA-Seq has become the standard approach to quantify and compare gene expression and alternative splicing in different conditions. In many cases the limiting factor is not the sequencing itself but the bioinformatic analysis. A variety of software tools exist that predict alternative splicing patterns from RNA-Seq data, but surprisingly, a systematic comparison of the predictions obtained from different pipelines has not been performed. Here we compare results from frequently used bioinformatic tools using a high-quality RNA-Seq dataset. We show that there is little overlap in the splicing changes predicted by different tools and that GO-term analysis of the splicing changes predicted by the individual targets yields very different results. Validation of bioinformatic predictions by RT-PCR suggest a high number of false positives in the splicing changes predicated by each pipeline, which probably dominates GO-term analysis. The validation rate is strongly increased for targets predicted by several tools, offering a strategy to reduce false positives. Based on these results we offer some guidelines that may contribute to make alternative splicing predictions more reliable and may thus increase the impact of conclusions drawn from RNA-Seq studies. Furthermore, we created rmappet, a nextflow pipeline that performs alternative splicing analysis using rMATS and Whippet with subsequent overlapping of the results, enabling robust splicing analysis with only one command (https://github.com/didrikolofsson/rmappet/).

UHMK1 is a novel splicing regulatory kinase

The U2AF Homology Motif Kinase 1 (UHMK1) is the only kinase that contains the U2AF homology motif, a common protein interaction domain among splicing factors. Through this motif, UHMK1 interacts with the splicing factors SF1 and SF3B1, known to participate in the 3' splice site recognition during the early steps of spliceosome assembly. Although UHMK1 phosphorylates these splicing factors invitro, the involvement of UHMK1 in RNA processing has not previously been demonstrated. Here, we identify novel putative substrates of this kinase and evaluate UHMK1 contribution to overall gene expression and splicing, by integrating global phosphoproteomics, RNA-seq, and bioinformatics approaches. Upon UHMK1 modulation, 163 unique phosphosites were differentially phosphorylated in 117 proteins, of which 106 are novel potential substrates of this kinase. Gene Ontology analysis showed enrichment of terms previously associated with UHMK1 function, such as mRNA splicing, cell cycle, cell division, and microtubule organization. The majority of the annotated RNA-related proteins are components of the spliceosome but are also involved in several steps of gene expression. Comprehensive analysis of splicing showed that UHMK1 affected over 270 alternative splicing events. Moreover, splicing reporter assay further supported UHMK1 function on splicing. Overall, RNA-seq data demonstrated that UHMK1 knockdown had a minor impact on transcript expression and pointed to UHMK1 function in epithelial-mesenchymal transition. Functional assays demonstrated that UHMK1 modulation affects proliferation, colony formation, and migration. Taken together, our data implicate UHMK1 as a splicing regulatory kinase, connecting protein regulation through phosphorylation and gene expression in key cellular processes.

Integrated Analysis of Transcriptome Changes in Osteoarthritis: Gene Expression, Pathways and Alternative Splicing.

Osteoarthritis (OA) is the most prevalent joint disease characterized by the degeneration of articular cartilage and the remodeling of its underlying bones, resulting in pain and loss of function in the knees and hips. As far as we know, no curative treatments are available except for the joint replacement. The precise molecular mechanisms which are involved in the degradation of cartilage matrix and development of osteoarthritis are still unclear. By analyzing RNA-seq data, we found the molecular changes at the transcriptome level such as alternative splicing, gene expression, and molecular pathways in OA knees cartilage. Expression analysis have identified 457 differential expressed genes including 266 up-regulated genes such as TNFSF15, ST6GALNAC5, TGFBI, ASPM, and TYM, and 191 down-regulated genes such as ADM, JUN, IRE2, PIGA, and MAFF. Gene set enrichment analysis (GSEA) analysis identified down-regulated pathways related to translation, transcription, immunity, PI3K/AKT, and circadian as well as disturbed pathways related to extracellular matrix and collagen. Splicing analysis identified 442 differential alternative splicing events within 284 genes in osteoarthritis, including genes involved in extracellular matrix (ECM) and alternative splicing, and TIA1 was identified as a key regulator of these splicing events. These findings provide insights into disease etiology, and offer favorable information to support the development of more effective interventions in response to the global clinical challenge of osteoarthritis.

The subgenome Saccharum spontaneum contributes to sugar accumulation in sugarcane as revealed by full-length transcriptomic analysis

IntroductionModern sugarcane cultivars (Saccharum spp. hybrids) derived from crosses between S. officinarum and S. spontaneum, with high-sugar traits and excellent stress tolerance inherited respectively. However, the contribution of the S. spontaneum subgenome to sucrose accumulation is still unclear. ObjectiveTo compensate for the absence of a high-quality reference genome, a transcriptome analysis method is needed to analyze the molecular basis of differential sucrose accumulation in sugarcane hybrids and to find clues to the contribution of the S. spontaneum subgenome to sucrose accumulation. MethodsPacBio full-length sequencing was used to complement genome annotation, followed by the identification of differential genes between the high and low sugar groups using differential alternative splicing analysis and differential expression analysis. At the subgenomic level, the factors responsible for differential sucrose accumulation were investigated from the perspective of transcriptional and post-transcriptional regulation. ResultsA full-length transcriptome annotated at the subgenomic level was provided, complemented by 263,378 allele-defined transcript isoforms and 139,405 alternative splicing (AS) events. Differential alternative splicing (DA) analysis and differential expression (DE) analysis identified differential genes between high and low sugar groups and explained differential sucrose accumulation factors by the KEGG pathways. In some gene models, different or even opposite expression patterns of alleles from the same gene were observed, reflecting the potential evolution of these alleles toward novel functions in polyploid sugarcane. Among DA and DE genes in the sucrose source-sink complex pathway, we found some alleles encoding sucrose accumulation-related enzymes derived from the S. spontaneum subgenome were differentially expressed or had DA events between the two contrasting sugarcane hybrids. ConclusionFull-length transcriptomes annotated at the subgenomic level could better characterize sugarcane hybrids, and the S. spontaneum subgenome was found to contribute to sucrose accumulation.

A Potential Therapy Using Antisense Oligonucleotides to Treat Autosomal Recessive Polycystic Kidney Disease.

(1) Background: Autosomal recessive polycystic kidney disease (ARPKD) is a rare ciliopathy characterized by progressively enlarged kidneys with fusiform dilatation of the collecting ducts. Loss-of-function mutations in the PKHD1 gene, which encodes fibrocystin/polyductin, cause ARPKD; however, an efficient treatment method and drug for ARPKD have yet to be found. Antisense oligonucleotides (ASOs) are short special oligonucleotides which function to regulate gene expression and alter mRNA splicing. Several ASOs have been approved by the FDA for the treatment of genetic disorders, and many are progressing at present. We designed ASOs to verify whether ASOs mediate the correction of splicing further to treat ARPKD arising from splicing defects and explored them as a potential treatment option. (2) Methods: We screened 38 children with polycystic kidney disease for gene detection using whole-exome sequencing (WES) and targeted next-generation sequencing. Their clinical information was investigated and followed up. The PKHD1 variants were summarized and analyzed, and association analysis was carried out to analyze the relationship between genotype and phenotype. Various bioinformatics tools were used to predict pathogenicity. Hybrid minigene analysis was performed as part of the functional splicing analysis. Moreover, the de novo protein synthesis inhibitor cycloheximide was selected to verify the degraded pathway of abnormal pre-mRNAs. ASOs were designed to rescue aberrant splicing, and this was verified. (3) Results: Of the 11 patients with PKHD1 variants, all of them exhibited variable levels of complications of the liver and kidneys. We found that patients with truncating variants and variants in certain regions had a more severe phenotype. Two splicing variants of the PKHD1 genotypes were studied via the hybrid minigene assay: variants c.2141-3T>C and c.11174+5G>A. These cause aberrant splicing, and their strong pathogenicity was confirmed. We demonstrated that the abnormal pre-mRNAs produced from the variants escaped from the NMD pathway with the use of the de novo protein synthesis inhibitor cycloheximide. Moreover, we found that the splicing defects were rescued by using ASOs, which efficiently induced the exclusion of pseudoexons. (4) Conclusion: Patients with truncating variants and variants in certain regions had a more severe phenotype. ASOs are a potential drug for treating ARPKD patients harboring splicing mutations of the PKHD1 gene by correcting the splicing defects and increasing the expression of the normal PKHD1 gene.

Dravet syndrome and hemorrhagic shock and encephalopathy syndrome associated with an intronic deletion of SCN1A

ObjectiveHemorrhagic shock and encephalopathy syndrome (HSES) is a serious condition that requires intensive care and is associated with a high mortality rate. However, its pathogenesis remains unclear. In the present study, a genetic analysis was performed to determine the genetic background of patients with clinically suspected Dravet syndrome (DS) who developed HSES. MethodsWhole exome sequencing was performed, followed by minigene analysis of the intron variant detected by whole exome sequencing to confirm its effect on splicing. ResultsWhole exome sequencing revealed a novel 21-bp deletion in intron 3 of SCN1A NM_001165963.4 (NC_000002.11:g.166073675_166073695del). This deletion was not found in the patient’s parents and was proven to be de novo. Minigene analysis revealed an aberrant mRNA lacking 40 and 106 bp from the 5′ end of exon 4 of SCN1A. Therefore, we diagnosed this case as DS due to the deletion in intron 3 of SCN1A. ConclusionsWe report a case of DS with HSES caused by a 21-bp deletion in the intron of SCN1A that was confirmed by minigene analysis. The present case met Levin's criteria for HSES and the splicing analysis of SCN1A is an important finding. This study has important implications for understanding HSES pathogenesis.

Alternative polyadenylation transcriptome-wide association study identifies APA-linked susceptibility genes in brain disorders

Alternative polyadenylation (APA) plays an essential role in brain development; however, current transcriptome-wide association studies (TWAS) largely overlook APA in nominating susceptibility genes. Here, we performed a 3′ untranslated region (3′UTR) APA TWAS (3′aTWAS) for 11 brain disorders by combining their genome-wide association studies data with 17,300 RNA-seq samples across 2,937 individuals. We identified 354 3′aTWAS-significant genes, including known APA-linked risk genes, such as SNCA in Parkinson’s disease. Among these 354 genes, ~57% are not significant in traditional expression- and splicing-TWAS studies, since APA may regulate the translation, localization and protein-protein interaction of the target genes independent of mRNA level expression or splicing. Furthermore, we discovered ATXN3 as a 3′aTWAS-significant gene for amyotrophic lateral sclerosis, and its modulation substantially impacted pathological hallmarks of amyotrophic lateral sclerosis in vitro. Together, 3′aTWAS is a powerful strategy to nominate important APA-linked brain disorder susceptibility genes, most of which are largely overlooked by conventional expression and splicing analyses.

Deviating Alternative Splicing as a Molecular Subtype of Microsatellite Stable Colorectal Cancer.

Colorectal cancer (CRC) is the third most commonly diagnosed cancer worldwide. Biomarkers to aid in prognostication and treatment decisions are in high demand, and to facilitate their development, a better understanding of the underlying biology of the highly heterogeneous disease is needed. A genome-scale alternative splicing (AS) analysis using RNA-sequencing data from primary microsatellite stable (MSS) CRCs from 127 patients was performed. Splice variant-specific expression levels of individual cancer samples were compared with the total set of samples, and a metric for a tumor sample's global amount of deviating AS was developed. This metric varied considerably across the cohort and ranged from 6 to 282 deviating AS events per tumor sample. A threshold of 45 or more deviating events was set to distinguish cancers with high (n = 44) and low (n = 83) levels of deviating AS. Patients with high amounts of AS deviations had significantly shorter time to relapse compared with patients with fewer deviations (P = .04). Furthermore, differential gene expression analysis revealed nine known cancer-critical genes that are significantly upregulated in samples with high amounts of deviating AS. Validation of the results in an independent cohort of MSS CRCs showed the same tendency toward shorter progression-free survival among the high-deviation group. In both cohorts, enrichment for growth factors was identified among upregulated genes associated with this phenotype. There is a large variation in the amount of deviating AS among MSS CRCs, and we provide evidence that those with high amounts of deviations represent different cancer biology.

A spectrum of clinical severity of recessive titinopathies in prenatal.

Variants in TTN are associated with a broad range of clinical phenotypes, from dominant adult-onset dilated cardiomyopathy to recessive infantile-onset myopathy. However, few foetal cases have been reported for multiple reasons. Next-generation sequencing has facilitated the prenatal identification of a growing number of suspected titinopathy variants. We investigated six affected foetuses from three families, completed the intrauterine course of the serial phenotypic spectrum of TTN, and discussed the genotype-phenotype correlations from a broader perspective. The recognizable prenatal feature onset at the second trimester was started with reduced movement, then contracture 3-6 weeks later, followed with/without hydrops, finally at late pregnancy was accompanied with polyhydramnio (major) or oligohydramnios. Two cases with typical arthrogryposis-hydrops sequences identified a meta-only transcript variant c.36203-1G>T. Deleterious transcriptional consequences of the substitution were verified by minigene splicing analysis. Case 3 identified a homozygous splicing variant in the constitutively expressed Z-disc. It presented a milder phenotype than expected, which was presumably saved by the isoform of corons. A summary of the foetal-onset titinopathy cases implied that variants in TTN present with a series of signs and a spectrum of clinical severity, which followed the dosage/positional effect; the meta-only transcript allele involvement may be a prerequisite for the development of fatal hydrops.

Experimental and numerical investigation of core splicing configurations in an advanced composite sandwich structural system

In this study, an investigation is presented on the effects of the honeycomb core splice gaps on the performance of advanced composite sandwich structures consisting of aluminum honeycomb core and carbon-fiber/epoxy facesheets. The experimental and numerical investigation explores the effects of core thickness, facesheet thickness, splice gap width, and the presence of core splice adhesive on the structural response of the composite sandwich structure. The specimens fabricated with the splices exposed to primarily bending and shear to understand further how the core splice interacts with the overall structural response under such loading conditions. In addition, 3D finite element models investigated the failure modes in the specimens by examining the strain fields for core-spliced sections of the sandwich structure under the experimentally observed loading conditions. Various failure modes are observed in the experiments ranging from core shear failure, facesheet debonding, or fracture depending on the location of the core splice, core, and the facesheet properties. The results illustrate the importance of proper detailing and analysis of core splices in composite structures made from sandwich designs.

NR3C1 overexpression regulates the expression and alternative splicing of inflammation-associated genes involved in PTSD

NR3C1-encoding glucocorticoid receptors have dual roles as RNA-binding protein and transcription factor. Recent studies revealed that NR3C1 might play an important role in the pathogenesis of PTSD (Post-traumatic stress disorder). However, its molecular mechanism remained unclear. In the present study, a neuronal cell model was constructed by transfecting a NR3C1-overexpressing plasmid pIRES-hrGFP-1a-NR3C1 or empty vector into HT22 cells. The changes in global transcription levels and alternative splicing events in HT22 cells after NR3C1 overexpression were analyzed by RNA sequencing. Compared with the empty vector control, the expression of inflammatory factors were differentially regulated by NR3C1, including genes involved in chemokine signal pathway, PI3K-Akt signal pathway, cytokine receptor interaction, neural ligand-receptor interaction and so on. In addition, NR3C1 regulated the alternative splicing of many genes involved in immune response, axon formation, stress response and inflammation. This study was the first to perform a transcriptome analysis of differential gene expression and alternative splicing in a NR3C1-overexpressing HT22 cell model. Our results suggested that NR3C1 could manipulate the expression of inflammatory transcription factors and their alternative splicing patterns, subsequently affecting the expression of downstream targets, may be leading to the onset of PTSD. This study will provide new insights in the NR3C1-mediated gene regulation in relation to PTSD.

MARVEL: an integrated alternative splicing analysis platform for single-cell RNA sequencing data.

Alternative splicing is an important source of heterogeneity underlying gene expression between individual cells but remains an understudied area due to the paucity of computational tools to analyze splicing dynamics at single-cell resolution. Here, we present MARVEL, a comprehensive R package for single-cell splicing analysis applicable to RNA sequencing generated from the plate- and droplet-based methods. We performed extensive benchmarking of MARVEL against available tools and demonstrated its utility by analyzing multiple publicly available datasets in diverse cell types, including in disease. MARVEL enables systematic and integrated splicing and gene expression analysis of single cells to characterize the splicing landscape and reveal biological insights.