Identification of Splicing Variation Associated with Parental Behavior in the Burying Beetles (Nicrophorus orbicollis)

Alternative Splicing: New Insights from Global Analyses

BackgroundThe vast diversification of proteins in eukaryotic cells has been related with multiple transcript isoforms from a single gene that result in alternative splicing (AS) of primary transcripts. Analysis of RNA sequencing data from expressed sequence tags and next generation RNA sequencing has been crucial for AS identification and genome-wide AS studies. For the identification of AS events from the related legume species Phaseolus vulgaris and Glycine max, 157 and 88 publicly available RNA-seq libraries, respectively, were analyzed.ResultsWe identified 85,570 AS events from P. vulgaris in 72% of expressed genes and 134,316 AS events in 70% of expressed genes from G. max. These were categorized in seven AS event types with intron retention being the most abundant followed by alternative acceptor and alternative donor, representing ~75% of all AS events in both plants. Conservation of AS events in homologous genes between the two species was analyzed where an overrepresentation of AS affecting 5’UTR regions was observed for certain types of AS events. The conservation of AS events was experimentally validated for 8 selected genes, through RT-PCR analysis. The different types of AS events also varied by relative position in the genes. The results were consistent in both species.ConclusionsThe identification and analysis of AS events are first steps to understand their biological relevance. The results presented here from two related legume species reveal high conservation, over ~15–20 MY of divergence, and may point to the biological relevance of AS.

Splicing Factor SRSF1 Promotes Tumorigenesis Via Oncogenic Splice-Switching and Predicts Poor Prognosis in Multiple Myeloma

Background: Surgery, adjuvant chemotherapy, and radiotherapy remain the primary treatment options for soft tissue sarcomas (STSs). Identifying ways to improve the prognosis of patients with STS is a considerable challenge. Evidence shows that the dysregulation of alternative splicing events is involved in tumor pathogenesis and progression. The present study objectives were to identify survival-associated alternative splicing (AS) events that could serve as prognostic biomarkers and potentially serve as tumor-selective STS drug targets. Methods: STS-specific percent spliced in (PSI) values for splice events in 206 STS samples were downloaded from The Cancer Genome Atlas (TCGA) SpliceSeq®. Prognostic analysis was performed on seven types of AS events to determine their prognostic value in STS patients, for which prediction models were constructed with the formula Risk score = n ∑ PSIi*βi i To explore how the AS events function, mRNA level, protein level, copy number alteration, mutation, and methylation of several STS key genes were detected from TCGA RNA-seq data and relevant protein data, as well as Gene Expression Omnibus gene chip data. Additionally, Spearman’s rank correlation coefficients were calculated to evaluate any correlation between splicing factor expression and the PSI values of survival-associated AS events. Findings: Of the 40,184 AS events for 3,064 genes, 10,439 events were found to significantly correlate with patient survival rates. The area under the time-dependent receiver operating characteristic curve for the prognostic predictor of STS 2-year overall survival was 0·826. Notably, the splicing events of certain STS key genes were significantly associated with STS 2-year overall survival in the present study, including exon skip (ES) events in MDM2 and EWSR1 and alternate terminator events in CDKN2A and HMGA2 for dedifferentiated liposarcoma; ES in MDM2 and alternate promoter events in CDKN2A for leiomyosarcoma; and ES in EWSR1 for undifferentiated pleomorphic sarcoma. Moreover, splicing correlation networks between AS events and splicing factors revealed that almost all of the AS events associated with favorable prognosis were negatively correlated with the expression of splicing factors. Likewise, most of the AS events associated with poor prognosis were positively correlated with the expression of splicing factors. Interpretation: In-depth analysis of alternative RNA splicing could provide new insights into the mechanisms of STS oncogenesis, with the potential for novel avenues for this type of cancer therapy. Funding Statement: The Innovation Project of Guangxi Graduate Education (YCBZ2018038), Medical Excellence Award Funded by the Creative Research Development Grant from the First Affiliated Hospital of Guangxi Medical University, Guangxi Medical University Training Program for Distinguished Young Scholars, and Innovation Project of Guangxi Graduate Education (YCBZ20190xx). Declaration of Interests: The authors declare that they have no conflicts of interest.

BackgroundSurgery, adjuvant chemotherapy, and radiotherapy are the primary treatment options for soft tissue sarcomas (STSs). However, identifying ways to improve the prognosis of patients with STS remains a considerable challenge. Evidence shows that the dysregulation of alternative splicing (AS) events is involved in tumor pathogenesis and progression. The present study objective was to identify survival-associated AS events that could serve as prognostic biomarkers and potentially serve as tumor-selective STS drug targets.MethodsSTS-specific ‘percent spliced in’ (PSI) values for splicing events in 206 STS samples were downloaded from The Cancer Genome Atlas SpliceSeq® database. Prognostic analyses were performed on seven types of AS events to determine their prognostic value in STS patients, for which prediction models were constructed with the risk score formula sumnolimits_{i}^{n} {PSIi; *;beta i}. Prediction models were also constructed to determine the prognostic value of AS events, and Spearman’s rank correlation coefficients were calculated to determine the degree of correlation between splicing factor expression and the PSI values.ResultsA total 10,439 events were found to significantly correlate with patient survival rates. The area under the time-dependent receiver operating characteristic curve for the prognostic predictor of STS overall survival was 0.826. Notably, the splicing events of certain STS key genes were significantly associated with STS 2-year overall survival in the present study, including exon skip (ES) events in MDM2 and EWSR1, alternate terminator events in CDKN2A and HMGA2 for dedifferentiated liposarcoma, ES in MDM2 and alternate promoter events in CDKN2A for leiomyosarcoma, and ES in EWSR1 for undifferentiated pleomorphic sarcoma. Moreover, splicing correlation networks between AS events and splicing factors revealed that almost all of the AS events showed negatively correlations with the expression of splicing factors.ConclusionAn in-depth analysis of alternative RNA splicing could provide new insights into the mechanisms of STS oncogenesis and the potential for novel approaches to this type of cancer therapy.

Dissecting the impact of regional identity and the oncogenic role of human-specific NOTCH2NL in an hESC model of H3.3G34R-mutant glioma

Alternative splicing is a molecular process that contributes greatly to the diversification of proteome and to gene functions. Understanding the mechanisms of stage-specific alternative splicing can provide a better understanding of the development of eukaryotes and the functions of different genes. Schistosoma japonicum is an infectious blood-dwelling trematode with a complex lifecycle that causes the tropical disease schistosomiasis. In this study, we analyzed the transcriptome of Schistosoma japonicum to discover alternative splicing events in this parasite, by applying RNA-seq to cDNA library of adults and schistosomula. Results were validated by RT-PCR and sequencing. We found 11,623 alternative splicing events among 7,099 protein encoding genes and average proportion of alternative splicing events per gene was 42.14%. We showed that exon skip is the most common type of alternative splicing events as found in high eukaryotes, whereas intron retention is the least common alternative splicing type. According to intron boundary analysis, the parasite possesses same intron boundaries as other organisms, namely the classic “GT-AG” rule. And in alternative spliced introns or exons, this rule is less strict. And we have attempted to detect alternative splicing events in genes encoding proteins with signal peptides and transmembrane helices, suggesting that alternative splicing could change subcellular locations of specific gene products. Our results indicate that alternative splicing is prevalent in this parasitic worm, and that the worm is close to its hosts. The revealed secretome involved in alternative splicing implies new perspective into understanding interaction between the parasite and its host.

Analysis of alternative splicing events by RNA sequencing in the ovaries of Xiang pig at estrous and diestrous

BackgroundAlternative splicing (AS) is a process which generates several distinct mRNA isoforms from the same gene by splicing different portions out of the precursor transcript. Due to the (patho-)physiological importance of AS, a complete inventory of AS is of great interest. While this is in reach for human and mammalian model organisms, our knowledge of AS in plants has remained more incomplete. Experimental approaches for monitoring AS are either based on transcript sequencing or rely on hybridization to DNA microarrays. Among the microarray platforms facilitating the discovery of AS events, tiling arrays are well-suited for identifying intron retention, the most prevalent type of AS in plants. However, analyzing tiling array data is challenging, because of high noise levels and limited probe coverage.ResultsIn this work, we present a novel method to detect intron retentions (IR) and exon skips (ES) from tiling arrays. While statistical tests have typically been proposed for this purpose, our method instead utilizes support vector machines (SVMs) which are appreciated for their accuracy and robustness to noise. Existing EST and cDNA sequences served for supervised training and evaluation. Analyzing a large collection of publicly available microarray and sequence data for the model plant A. thaliana, we demonstrated that our method is more accurate than existing approaches. The method was applied in a genome-wide screen which resulted in the discovery of 1,355 IR events. A comparison of these IR events to the TAIR annotation and a large set of short-read RNA-seq data showed that 830 of the predicted IR events are novel and that 525 events (39%) overlap with either the TAIR annotation or the IR events inferred from the RNA-seq data.ConclusionsThe method developed in this work expands the scarce repertoire of analysis tools for the identification of alternative mRNA splicing from whole-genome tiling arrays. Our predictions are highly enriched with known AS events and complement the A. thaliana genome annotation with respect to AS. Since all predicted AS events can be precisely attributed to experimental conditions, our work provides a basis for follow-up studies focused on the elucidation of the regulatory mechanisms underlying tissue-specific and stress-dependent AS in plants.

Alternative splicing (AS) is a widely recognized post-transcriptional regulatory mechanism that plays a crucial role in plant evolution and environmental adaptation. In this study, five representative Panax species were systematically analyzed to examine the evolutionary conservation and functional characteristics of AS events. Results revealed an expansion in the number of AS events and associated genes across the Panax species, accompanied by a genome-wide shift in splicing types from a dominance of intron retention (IR) to an increase in exon skipping (ES), alternative donor (A5), and alternative acceptor (A3) events. Splicing preferences were also found to diverge among allotetraploid species, which exhibited more complex AS patterns. The genomic features of IR and ES events, such as GC content and length of the sequence involved in AS, were highly conserved among Panax species of different ploidy levels (diploid vs. allotetraploid). Genes harboring conserved IR events across all five species were identified, and functional annotation indicated that these genes are primarily involved in chromatin modification and RNA splicing-related processes. This study elucidates the dynamic remodeling of AS during the evolution of Panax and provides important insights into the evolutionary adaptive mechanisms of AS in plants.

Sacred lotus (<i>Nelumbo nucifera</i>) is a typical aquatic plant, belonging to basal eudicot plant, which is ideal for genome and genetic evolutionary study. Understanding lotus gene diversity is important for the study of molecular genetics and breeding. In this research, public RNA-seq data and the annotated reference genome were used to identify the genes in lotus. A total of 26,819 consensus and 1,081 novel genes were identified. Meanwhile, a comprehensive analysis of gene alternative splicing events was conducted, and a total of 19,983 “internal” alternative splicing (AS) events and 14,070 “complete” AS events were detected in 5,878 and 5,881 multi-exon expression genes, respectively. Observations made from the AS events show the predominance of intron retention (IR) subtype of AS events representing 33%. IR is followed by alternative acceptor (AltA), alternative donor (AltD) and exon skipping (ES), highlighting the universality of the intron definition model in plants. In addition, functional annotations of the gene with AS indicated its relationship to a number of biological processes such as cellular process and metabolic process, showing the key role for alternative splicing in influencing the growth and development of lotus. The results contribute to a better understanding of the current gene diversity in lotus, and provide an abundant resource for future functional genome analysis in lotus.

G protein-coupled receptors (GPCRs) are the largest signaling family in the genome, serve an expansive array of functions, and are targets for approximately 50% of current therapeutics. In many tissues, such as airway smooth muscle (ASM), complex, unexpected, or paradoxical responses to agonists/antagonists occur without known mechanisms. We hypothesized that ASM express many more GPCRs than predicted, and that these undergo substantial alternative splicing, creating a highly diversified receptor milieu. Transcript arrays were designed detecting 434 GPCRs and their predicted splice variants. In this cell type, 353 GPCRs were detected (including 111 orphans), with expression levels varying by approximately 900-fold. Receptors used for treating airway disease were expressed lower than others with similar signaling properties, indicating potentially more effective targets. A disproportionate number of Class-A peptide-group receptors, and those coupling to G(q)/(11) or G(s) (vs. G(i)), was found. Importantly, 192 GPCRs had, on average, five different expressed receptor isoforms because of splicing events, including alternative splice donors and acceptors, novel introns, intron retentions, exon(s) skips, and novel exons, with the latter two events being most prevalent. The consequences of splicing were further investigated with the leukotriene B4 receptor, known for its aberrant responsiveness in lung. We found transcript expression of three variants because of alternative donor and acceptor splice sites, representing in-frame deletions of 38 and 100 aa, with protein expression of all three isoforms. Thus, alternative splicing, subject to conditional, temporal, and cell-type regulation, is a major mechanism that diversifies the GPCR superfamily, creating local recepteromes with specialized environments.

Alternative splicing (AS) can increase transcriptome diversity, protein diversity and protein yield, and is an important mechanism to regulate plant responses to stress. Oilseed rape (Brassica napus L.), one of the main oil crops in China, shows higher sensitivity to boron (B) deficiency than other species. Here, we demonstrated AS changes that largely increased the diversity of the mRNA expressed in response to B deficiency in B. napus. Each gene had two or more transcripts on average. A total of 33.3% genes in both Qingyou10 (QY10, B-efficient cultivar) and Westar10 (W10, B-inefficient cultivar) showed AS in both B conditions. The types of AS events were mainly intron retention, 3′ alternative splice site, 5′ alternative splice site and exon skipping. The tolerance ability of QY10 was higher than that of W10, possibly because there were far more differential alternative splicing (DAS) genes identified in QY10 at low B conditions than in W10. The number of genes with both DAS and differentially expressed (DE) was far lower than that of the genes that were either with DAS or DE in QY10 and W10, suggesting that the DAS and DE genes were independent. Four Serine/Arginine-rich (SR) splicing factors, BnaC06g14780D, BnaA01g14750D, BnaA06g15930D and BnaC01g41640D, underwent differentially alternative splicing in both cultivars. There existed gene–gene interactions between BnaC06g14780D and the genes associated with the function of B in oilseed rape at low B supply. This suggests that oilseed rape could regulate the alterative pre-mRNA splicing of SR protein related genes to increase the plant tolerance to B deficiency.

BackgroundGenome-wide computational analysis of alternative splicing (AS) in several flowering plants has revealed that pre-mRNAs from about 30% of genes undergo AS. Chlamydomonas, a simple unicellular green alga, is part of the lineage that includes land plants. However, it diverged from land plants about one billion years ago. Hence, it serves as a good model system to study alternative splicing in early photosynthetic eukaryotes, to obtain insights into the evolution of this process in plants, and to compare splicing in simple unicellular photosynthetic and non-photosynthetic eukaryotes. We performed a global analysis of alternative splicing in Chlamydomonas reinhardtii using its recently completed genome sequence and all available ESTs and cDNAs.ResultsOur analysis of AS using BLAT and a modified version of the Sircah tool revealed AS of 498 transcriptional units with 611 events, representing about 3% of the total number of genes. As in land plants, intron retention is the most prevalent form of AS. Retained introns and skipped exons tend to be shorter than their counterparts in constitutively spliced genes. The splice site signals in all types of AS events are weaker than those in constitutively spliced genes. Furthermore, in alternatively spliced genes, the prevalent splice form has a stronger splice site signal than the non-prevalent form. Analysis of constitutively spliced introns revealed an over-abundance of motifs with simple repetitive elements in comparison to introns involved in intron retention. In almost all cases, AS results in a truncated ORF, leading to a coding sequence that is around 50% shorter than the prevalent splice form. Using RT-PCR we verified AS of two genes and show that they produce more isoforms than indicated by EST data. All cDNA/EST alignments and splice graphs are provided in a website at http://combi.cs.colostate.edu/as/chlamy.ConclusionsThe extent of AS in Chlamydomonas that we observed is much smaller than observed in land plants, but is much higher than in simple unicellular heterotrophic eukaryotes. The percentage of different alternative splicing events is similar to flowering plants. Prevalence of constitutive and alternative splicing in Chlamydomonas, together with its simplicity, many available public resources, and well developed genetic and molecular tools for this organism make it an excellent model system to elucidate the mechanisms involved in regulated splicing in photosynthetic eukaryotes.

For most of our 25,000 genes, the removal of introns by pre-messenger RNA (pre-mRNA) splicing represents an essential step toward the production of functional messenger RNAs (mRNAs). Alternative splicing of a single pre-mRNA results in the production of different mRNAs. Although complex organisms use alternative splicing to expand protein function and phenotypic diversity, patterns of alternative splicing are often altered in cancer cells. Alternative splicing contributes to tumorigenesis by producing splice isoforms that can stimulate cell proliferation and cell migration or induce resistance to apoptosis and anticancer agents. Cancer-specific changes in splicing profiles can occur through mutations that are affecting splice sites and splicing control elements, and also by alterations in the expression of proteins that control splicing decisions. Recent progress in global approaches that interrogate splicing diversity should help to obtain specific splicing signatures for cancer types. The development of innovative approaches for annotating and reprogramming splicing events will more fully establish the essential contribution of alternative splicing to the biology of cancer and will hopefully provide novel targets and anticancer strategies. Metazoan genes are usually made up of several exons interrupted by introns. The introns are removed from the pre-mRNA by RNA splicing. In conjunction with other maturation steps, such as capping and polyadenylation, the spliced mRNA is then transported to the cytoplasm to be translated into a functional protein. The basic mechanism of splicing requires accurate recognition of each extremity of each intron by the spliceosome. Introns are identified by the binding of U1 snRNP to the 5' splice site and the U2AF65/U2AF35 complex to the 3' splice site. Following these interactions, other proteins and snRNPs are recruited to generate the complete spliceosomal complex needed to excise the intron. While many introns are constitutively removed by the spliceosome, other splice junctions are not used systematically, generating the phenomenon of alternative splicing. Alternative splicing is therefore the process by which a single species of pre-mRNA can be matured to produce different mRNA molecules (Fig.1). Depending on the number and types of alternative splicing events, a pre-mRNA can generate from two to several thousands different mRNAs leading to the production of a corresponding number of proteins. It is now believed that the expression of at least 70% of human genes is subjected to alternative splicing, implying an enormous contribution to proteomic diversity, and by extension, to the development and the evolution of complex animals. Defects in splicing have been associated with human diseases (Caceres and Kornblihtt, Trends Genet 18(4):186-93, 2002, Cartegni et al., Nat Rev Genet 3(4):285-98, 2002, Pagani and Baralle, Nat Rev Genet 5(5):389-96, 2004), including cancer (Brinkman, Clin Biochem 37(7):584-94, 2004, Venables, Bioessays 28(4):378-86, 2006, Srebrow and Kornblihtt, J Cell Sci 119(Pt 13):2635-2641, 2006, Revil et al., Bull Cancer 93(9):909-919, 2006, Venables, Transworld Res Network, 2006, Pajares et al., Lancet Oncol 8(4):349-57, 2007, Skotheim and Nees, Int J Biochem Cell Biol 39:1432-1449, 2007). Numerous studies have now confirmed the existence of specific differences in the alternative splicing profiles between normal and cancer tissues. Although there are a few cases where specific mutations are the primary cause for these changes, global alterations in alternative splicing in cancer cells may be primarily derived from changes in the expression of RNA-binding proteins that control splice site selection. Overall, these cancer-specific differences in alternative splicing offer an immense potential to improve the diagnosis and the prognosis of cancer. This review will focus on the functional impact of cancer-associated alternative splicing variants, the molecular determinants that alter the splicing decisions in cancer cells, and future therapeutic strategies.