Comprehensive Analysis Of Alternative Splicing Research Articles

Comprehensive analysis of alternative splicing across multiple transcriptomic cohorts reveals prognostic signatures in prostate cancer

BackgroundAlternative splicing (AS) plays a crucial role in transcriptomic diversity and is a hallmark of cancer that profoundly influences the development and progression of prostate cancer (PCa), a prevalent and potentially life-limiting cancer among men. Accumulating evidence has highlighted the association between AS dysregulation and the onset and progression of PCa. However, a comprehensive and integrative analysis of AS profiles at the event level, utilising data from multiple high-throughput cohorts and evaluating the prognosis of PCa progression, remains lacking and calls for thorough exploration.ResultsWe identified a differentially expressed retained intron event in ZWINT across three distinct cohorts, encompassing an original array-based dataset profiled by us previously and two RNA sequencing (RNA-seq) datasets. Subsequent in-depth analyses of these RNA-seq datasets revealed 141 altered events, of which 21 demonstrated a significant association with patients’ biochemical recurrence-free survival (BCRFS). We formulated an AS event-based prognostic signature, capturing six pivotal events in genes CYP4F12, NFATC4, PIGO, CYP3A5, ALS2CL, and FXYD3. This signature effectively differentiated high-risk patients diagnosed with PCa, who experienced shorter BCRFS, from their low-risk counterparts. Notably, the signature's predictive power surpassed traditional clinicopathological markers in forecasting 5-year BCRFS, demonstrating robust performance in both internal and external validation sets. Lastly, we constructed a novel nomogram that integrates patients’ Gleason scores with pathological tumour stages, demonstrating improved prognostication of BCRFS.ConclusionsPrediction of clinical progression remains elusive in PCa. This research uncovers novel splicing events associated with BCRFS, augmenting existing prognostic tools, thus potentially refining clinical decision-making.

Alternative splicing event associated with immunological features in bladder cancer.

Bladder cancer (BLCA) is the most prevalent urinary tumor with few treatments. Alternative splicing (AS) is closely related to tumor development and tumor immune microenvironment. However, the comprehensive analysis of AS and prognosis and immunological features in BLCA is still lacking. In this study, we downloaded RNA-Seq data and clinical information from The Cancer Genome Atlas (TCGA) database, and AS events were acquired from the TCGA Splice-seq. A total of eight prognostic AS events (C19orf57|47943|ES, ANK3|11845|AP, AK9|77203|AT, GRIK2|77096|AT, DYM|45472|ES, PTGER3|3415|AT, ACTG1|44120|RI, and TRMU|62711|AA) were identified by univariate analysis and least absolute shrinkage and selection operator (LASSO) regression analysis to construct a risk score model. The Kaplan-Meier analysis revealed that the high-risk group had a worse prognosis compared with the low-risk group. The area under the receiver operating characteristic (ROC) curves (AUCs) for this risk score model in 1, 3, and 5 years were 0.698, 0.742, and 0.772, respectively. One of the prognostic AS event-related genes, TRMU, was differentially expressed between tumor and normal tissues in BLCA. The single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT algorithm showed that both the risk score model and TRMU were significantly associated with tumor immune microenvironment and immune status (immune cells, immune-related pathway, and immune checkpoint) in BLCA patients. The TIMER database confirmed the relationship between the expression of TRMU and immune cells and checkpoint genes. Furthermore, Cytoscape software 3.8.0 was used to construct the regulatory network between AS and splicing factors (SFs). Our study demonstrated that AS events were powerful biomarkers to predict the prognosis and immune status in BLCA, which may be potential therapeutic targets in BLCA.

High-throughput total RNA sequencing in single cells using VASA-seq

Most methods for single-cell transcriptome sequencing amplify the termini of polyadenylated transcripts, capturing only a small fraction of the total cellular transcriptome. This precludes the detection of many long non-coding, short non-coding and non-polyadenylated protein-coding transcripts and hinders alternative splicing analysis. We, therefore, developed VASA-seq to detect the total transcriptome in single cells, which is enabled by fragmenting and tailing all RNA molecules subsequent to cell lysis. The method is compatible with both plate-based formats and droplet microfluidics. We applied VASA-seq to more than 30,000 single cells in the developing mouse embryo during gastrulation and early organogenesis. Analyzing the dynamics of the total single-cell transcriptome, we discovered cell type markers, many based on non-coding RNA, and performed in vivo cell cycle analysis via detection of non-polyadenylated histone genes. RNA velocity characterization was improved, accurately retracing blood maturation trajectories. Moreover, our VASA-seq data provide a comprehensive analysis of alternative splicing during mammalian development, which highlighted substantial rearrangements during blood development and heart morphogenesis.

FishExp: A comprehensive database and analysis platform for gene expression and alternative splicing of fish species

The publicly archived RNA-seq data has grown exponentially, while its valuable information has not yet been fully discovered and utilized, such as alternative splicing and its integration with gene expression. This is especially true for fish species which play important roles in ecology, research and the food industry. Furthermore, there is a lack of online platform to analyze users’ new data individually and jointly with existing data for the comprehensive analysis of alternative splicing and gene expression. Here, we present FishExp, a web-based data platform covering gene expression and alternative splicing in 26,081 RNA-seq experiments from 44 fishes. It allows users to query the data in a variety of ways, including gene identifier/symbol, functional term, and BLAST alignment. Moreover, users can customize experiments and tools to perform differential/specific expression and alternative splicing analysis, co-expression and cross-species analysis. In addition, functional enrichment is provided to confer biological significance. Notably, users are allowed to submit their own data and perform various analyses using the new data alone or alongside existing data in FishExp. Results of retrieval and analysis can be visualized on the gene-, transcript- and splicing event-level webpage in a highly interactive and intuitive manner. All data in FishExp can be downloaded for more in-depth analysis. The manually curated sample information, uniform data processing and various tools make it efficient for users to gain new insights from these large data sets, facilitating scientific hypothesis generation. FishExp is freely accessible at https://bioinfo.njau.edu.cn/fishExp.

Comprehensive Analysis of Alternative Splicing in Gastric Cancer Identifies Epithelial-Mesenchymal Transition Subtypes Associated with Survival.

This study presents a comprehensive analysis of alternative splicing in the context of patient classification, molecular mechanisms, and prognosis in gastric cancer.

Classification of glioma based on prognostic alternative splicing

BackgroundPreviously developed classifications of glioma have provided enormous advantages for the diagnosis and treatment of glioma. Although the role of alternative splicing (AS) in cancer, especially in glioma, has been validated, a comprehensive analysis of AS in glioma has not yet been conducted. In this study, we aimed at classifying glioma based on prognostic AS.MethodsUsing the TCGA glioblastoma (GBM) and low-grade glioma (LGG) datasets, we analyzed prognostic splicing events. Consensus clustering analysis was conducted to classified glioma samples and correlation analysis was conducted to characterize regulatory network of splicing factors and splicing events.ResultsWe analyzed prognostic splicing events and proposed novel splicing classifications across pan-glioma samples (labeled pST1–7) and across GBM samples (labeled ST1–3). Distinct splicing profiles between GBM and LGG were observed, and the primary discriminator for the pan-glioma splicing classification was tumor grade. Subtype-specific splicing events were identified; one example is AS of zinc finger proteins, which is involved in glioma prognosis. Furthermore, correlation analysis of splicing factors and splicing events identified SNRPB and CELF2 as hub splicing factors that upregulated and downregulated oncogenic AS, respectively.ConclusionA comprehensive analysis of AS in glioma was conducted in this study, shedding new light on glioma heterogeneity and providing new insights into glioma diagnosis and treatment.

Comprehensive characterization of the alternative splicing landscape in head and neck squamous cell carcinoma reveals novel events associated with tumorigenesis and the immune microenvironment.

Alternative splicing (AS) has emerged as a key event in tumor development and microenvironment formation. However, comprehensive analysis of AS and its clinical significance in head and neck squamous cell carcinoma (HNSC) is urgently required. Methods: Genome-wide profiling of AS events using RNA-Seq data from The Cancer Genome Atlas (TCGA) program was performed in a cohort of 464 patients with HNSC. Cancer-associated AS events (CASEs) were identified between paired HNSC and adjacent normal tissues and evaluated in functional enrichment analysis. Splicing networks and prognostic models were constructed using bioinformatics tools. Unsupervised clustering of the CASEs identified was conducted and associations with clinical, molecular and immune features were analyzed. Results: We detected a total of 32,309 AS events and identified 473 CASEs in HNSC; among these, 91 were validated in an independent cohort (n = 15). Functional protein domains were frequently altered, especially by CASEs affecting cancer drivers, such as PCSK5. CASE parent genes were significantly enriched in pathways related to HNSC and the tumor immune microenvironment, such as the viral carcinogenesis (FDR < 0.001), Human Papillomavirus infection (FDR < 0.001), chemokine (FDR < 0.001) and T cell receptor (FDR < 0.001) signaling pathways. CASEs enriched in immune-related pathways were closely associated with immune cell infiltration and cytolytic activity. AS regulatory networks suggested a significant association between splicing factor (SF) expression and CASEs and might be regulated by SF methylation. Eighteen CASEs were identified as independent prognostic factors for overall and disease-free survival. Unsupervised clustering analysis revealed distinct correlations between AS-based clusters and prognosis, molecular characteristics and immune features. Immunogenic features and immune subgroups cooperatively depict the immune features of AS-based clusters. Conclusion: This comprehensive genome-wide analysis of the AS landscape in HNSC revealed novel AS events related to carcinogenesis and immune microenvironment, with implications for prognosis and therapeutic responses.

Comprehensive Analysis of Alternative Splicing in Gastric Cancer and Their Clinical Significance

Comprehensive Analysis of Alternative Splicing in Gastric Cancer and Their Clinical Significance

Comprehensive Analysis of Alternative Splicing Across Tumors from 8,705 Patients

Our comprehensive analysis of alternative splicing across 32 The Cancer Genome Atlas cancer types from 8,705 patients detects alternative splicing events and tumor variants by reanalyzing RNA and whole-exome sequencing data. Tumors have up to 30% more alternative splicing events than normal samples. Association analysis of somatic variants with alternative splicing events confirmed known trans associations with variants in SF3B1 and U2AF1 and identified additional trans-acting variants (e.g., TADA1, PPP2R1A). Many tumors have thousands of alternative splicing events not detectable in normal samples; on average, we identified ≈930 exon-exon junctions ("neojunctions") in tumors not typically found in GTEx normals. From ClinicalProteomic Tumor Analysis Consortium data available for breast and ovarian tumor samples, we confirmed≈1.7 neojunction- and ≈0.6 single nucleotide variant-derived peptides per tumor sample that are also predicted major histocompatibility complex-I binders ("putative neoantigens").

A Comprehensive Analysis of Alternative Splicing in Paleopolyploid Maize.

Identifying and characterizing alternative splicing (AS) enables our understanding of the biological role of transcript isoform diversity. This study describes the use of publicly available RNA-Seq data to identify and characterize the global diversity of AS isoforms in maize using the inbred lines B73 and Mo17, and a related species, sorghum. Identification and characterization of AS within maize tissues revealed that genes expressed in seed exhibit the largest differential AS relative to other tissues examined. Additionally, differences in AS between the two genotypes B73 and Mo17 are greatest within genes expressed in seed. We demonstrate that changes in the level of alternatively spliced transcripts (intron retention and exon skipping) do not solely reflect differences in total transcript abundance, and we present evidence that intron retention may act to fine-tune gene expression across seed development stages. Furthermore, we have identified temperature sensitive AS in maize and demonstrate that drought-induced changes in AS involve distinct sets of genes in reproductive and vegetative tissues. Examining our identified AS isoforms within B73 × Mo17 recombinant inbred lines (RILs) identified splicing QTL (sQTL). The 43.3% of cis-sQTL regulated junctions are actually identified as alternatively spliced junctions in our analysis, while 10 Mb windows on each side of 48.2% of trans-sQTLs overlap with splicing related genes. Using sorghum as an out-group enabled direct examination of loss or conservation of AS between homeologous genes representing the two subgenomes of maize. We identify several instances where AS isoforms that are conserved between one maize homeolog and its sorghum ortholog are absent from the second maize homeolog, suggesting that these AS isoforms may have been lost after the maize whole genome duplication event. This comprehensive analysis provides new insights into the complexity of AS in maize.

Alternative splicing enhances transcriptome complexity in desiccating seeds.

Before being dispersed in the environment, mature seeds need to be dehydrated. The survival of seeds after dispersal depends on their low hydration in combination with high desiccation tolerance. These characteristics are established during seed maturation. Some key seed maturation genes have been reported to be regulated by alternative splicing (AS). However, so far AS was described only for single genes and a comprehensive analysis of AS during seed maturation has been lacking. We investigated gene expression and AS during Arabidopsis thaliana seed development at a global level, before and after desiccation. Bioinformatics tools were developed to identify differentially spliced regions within genes. Our data suggest the importance and shows the peculiar features of AS during seed desiccation. We identified AS in 34% of genes that are expressed at both timepoints before and after desiccation. Most of these AS transcript variants had not been found before in other tissues. Among the AS genes some seed master regulators could be found. Interestingly, 6% of all expressed transcripts were not transcriptionally regulated during desiccation, but only modified by AS. We propose that AS should be more routinely taken into account in the analysis of transcriptomic data to prevent overlooking potentially important regulators.

Comprehensive analysis of alternative splicing in Digitalis purpurea by strand-specific RNA-Seq.

Digitalis purpurea (D. purpurea) is one of the most important medicinal plants and is well known in the treatment of heart failure because of the cardiac glycosides that are its main active compounds. However, in the absence of strand specific sequencing information, the post-transcriptional mechanism of gene regulation in D. purpurea thus far remains unknown. In this study, a strand-specific RNA-Seq library was constructed and sequenced using Illumina HiSeq platforms to characterize the transcriptome of D. purpurea with a focus on alternative splicing (AS) events and the effect of AS on protein domains. De novo RNA-Seq assembly resulted in 48,475 genes. Based on the assembled transcripts, we reported a list of 3,265 AS genes, including 5,408 AS events in D. purpurea. Interestingly, both glycosyltransferases and monooxygenase, which were involved in the biosynthesis of cardiac glycosides, are regulated by AS. A total of 2,422 AS events occurred in coding regions, and 959 AS events were located in the regions of 882 unique protein domains, which could affect protein function. This D. purpurea transcriptome study substantially increased the expressed sequence resource and presented a better understanding of post-transcriptional regulation to further facilitate the medicinal applications of D. purpurea for human health.

Extensive Alternative Splicing of the Repressor Element Silencing Transcription Factor Linked to Cancer

The repressor element silencing transcription factor (REST) is a coordinate transcriptional and epigenetic regulator which functions as a tumor suppressor or an oncogene depending on cellular context, and a truncated splice variant REST4 has been linked to various types of cancer. We performed a comprehensive analysis of alternative splicing (AS) of REST by rapid amplification of cDNA ends and PCR amplification of cDNAs from various tissues and cell lines with specific primers. We identified 8 novel alternative exons including an alternate last exon which doubles the REST gene boundary, along with numerous 5′/3′ splice sites and ends in the constitutive exons. With the combination of various splicing patterns (e.g. exon skipping and alternative usage of the first and last exons) that are predictive of altered REST activity, at least 45 alternatively spliced variants of coding and non-coding mRNA were expressed in a species- and cell-type/tissue-specific manner with individual differences. By examining the repertoire of REST pre-mRNA splicing in 27 patients with kidney, liver and lung cancer, we found that all patients without exception showed differential expression of various REST splice variants between paired tumor and adjacent normal tissues, with striking cell-type/tissue and individual differences. Moreover, we revealed that exon 3 skipping, which causes no frame shift but loss of a domain essential for nuclear translocation, was affected by pioglitazone, a highly selective activator of the peroxisome proliferator-activated receptor gamma (PPARγ) which contributes to cell differentiation and tumorigenesis besides its metabolic actions. Accordingly, this study demonstrates an extensive AS of REST pre-mRNA which redefines REST gene boundary and structure, along with a general but differential link between REST pre-mRNA splicing and various types of cancer. These findings advance our understanding of the complex, context-dependent regulation of REST gene expression and function, and provide potential biomarkers and therapeutic targets for cancer.

Comprehensive Analysis of Alternative Splicing and Functionality in Neuronal Differentiation of P19 Cells

BackgroundAlternative splicing, which produces multiple mRNAs from a single gene, occurs in most human genes and contributes to protein diversity. Many alternative isoforms are expressed in a spatio-temporal manner, and function in diverse processes, including in the neural system.Methodology/Principal FindingsThe purpose of the present study was to comprehensively investigate neural-splicing using P19 cells. GeneChip Exon Array analysis was performed using total RNAs purified from cells during neuronal cell differentiation. To efficiently and readily extract the alternative exon candidates, 9 filtering conditions were prepared, yielding 262 candidate exons (236 genes). Semiquantitative RT-PCR results in 30 randomly selected candidates suggested that 87% of the candidates were differentially alternatively spliced in neuronal cells compared to undifferentiated cells. Gene ontology and pathway analyses suggested that many of the candidate genes were associated with neural events. Together with 66 genes whose functions in neural cells or organs were reported previously, 47 candidate genes were found to be linked to 189 events in the gene-level profile of neural differentiation. By text-mining for the alternative isoform, distinct functions of the isoforms of 9 candidate genes indicated by the result of Exon Array were confirmed.Conclusions/SignificanceAlternative exons were successfully extracted. Results from the informatics analyses suggested that neural events were primarily governed by genes whose expression was increased and whose transcripts were differentially alternatively spliced in the neuronal cells. In addition to known functions in neural cells or organs, the uninvestigated alternative splicing events of 11 genes among 47 candidate genes suggested that cell cycle events are also potentially important. These genes may help researchers to differentiate the roles of alternative splicing in cell differentiation and cell proliferation.

Structural organization of the sea urchin DNA (cytosine-5)-methyltransferase gene and characterization of five alternative spliced transcripts

Sea urchin DNA (cytosine-5)-methyltransferase ( Dnmt1) that is responsible for maintenance of DNA methylation patterns clearly shares similarity with various Dnmt1s identified in vertebrates. In this study, we determined the structure of the sea urchin Dnmt1 gene by screening a genomic library of the sea urchin Paracentrotus lividus with the complementary DNA (cDNA) as probe. Analysis of the positive clones demonstrated that the Dnmt1 gene consists of 34 exons and 33 introns spanning a distance of 35 kb. All exon-intron junction sequences agree with the GT/AG consensus with the exception of the 3′ acceptor site of intron 8 where CT replaces AG consensus. The differences in the total number of exons between sea urchin and mouse genes reside mainly in the N-terminal region of the protein (exons 5–7 of the sea urchin, exons 5–12 of the mouse) where there is very low similarity in the amino acid sequence. By reverse transcription-polymerase chain reaction using oligonucleotides spanning different regions of the cDNA we carried out a comprehensive analysis of alternative splicing of the Dnmt1 messenger RNA (mRNA) in sea urchin embryos at different stages of development. We demonstrated the presence of at least five alternative spliced mRNAs that are regulated during development and are translated in truncated or deleted proteins.

Alternative splicing of transcripts for the alpha 3 chain of mouse collagen VI: identification of an abundant isoform lacking domains N7-N10 in mouse and human.

Three distinct alpha chains form the collagen VI monomer, the alpha 3(VI) chain being much larger than the alpha 1(VI) and alpha 2(VI) chains. The alpha 3(VI) chain has 10 von Willebrand Factor type A domains of approximately 200 amino acids at the N-terminus (N1-N10) compared with only one such domain in the alpha 1(VI) and alpha 2(VI) chains. Domains N10, N9, N7 and N3 of the alpha 3(VI) chain are subject to alternative splicing in chick and/or human tissues, indicating the possibility of isoforms that have different functions depending on which N-terminal domains are included or excluded. In this study we have PCR amplified and sequenced mouse alpha 3(VI) cDNA encoding the N2-N10 domains. By reverse transcription-PCR using oligonucleotides spanning different regions of the cDNA we have undertaken a comprehensive analysis of alternative splicing of the alpha 3(VI) mRNA in embryonic and adult mouse tissues. We demonstrate that domains N10, N9 and N7 are also subject to alternative splicing in mouse tissues and in addition identify an abundant novel variant transcript that lacks all four N-terminal domains (N7-N10) in mouse tissues and human cells. We also identify less abundant transcripts that lack a large part of the N3 domain, and transcripts lacking the entire N5 domain. Using specific RNase protection assays we show that the shorter transcripts containing domains (N8+N7+N6), (N8+N6) and N6 are present at higher levels than transcripts containing the N10 and/or N9 domains, with tissue-specific variation in the levels of variant transcripts. These studies demonstrate a larger range of collagen VI protein variants than previously described.