jSplice: a high-performance method for accurate prediction of alternative splicing events and its application to large-scale renal cancer transcriptome data

The Role of Alternative Splicing in Cancer: Regulatory Mechanism, Therapeutic Strategy, and Bioinformatics Application

[Formula: see text] Alternative splicing (AS) can generate distinct transcripts and subsequent isoforms that play differential functions from the same pre-mRNA. Recently, increasing numbers of studies have emerged, unmasking the association between AS and cancer. In this review, we arranged AS events that are closely related to cancer progression and presented promising treatments based on AS for cancer therapy. Obtaining proliferative capacity, acquiring invasive properties, gaining angiogenic features, shifting metabolic ability, and getting immune escape inclination are all splicing events involved in biological processes. Spliceosome-targeted and antisense oligonucleotide technologies are two novel strategies that are hopeful in tumor therapy. In addition, bioinformatics applications based on AS were summarized for better prediction and elucidation of regulatory routines mingled in. Together, we aimed to provide a better understanding of complicated AS events associated with cancer biology and reveal AS a promising target of cancer treatment in the future.

Structural and functional insights into CWC27/CWC22 heterodimer linking the exon junction complex to spliceosomes

Human CWC27 is an uncharacterized splicing factor and mutations in its gene are linked to retinal degeneration and other developmental defects. We identify the splicing factor CWC22 as the major CWC27 partner. Both CWC27 and CWC22 are present in published B^act spliceosome structures, but no interacting domains are visible. Here, the structure of a CWC27/CWC22 heterodimer bound to the exon junction complex (EJC) core component eIF4A3 is solved at 3Å-resolution. According to spliceosomal structures, the EJC is recruited in the C complex, once CWC27 has left. Our 3D structure of the eIF4A3/CWC22/CWC27 complex is compatible with the B^act spliceosome structure but not with that of the C complex, where a CWC27 loop would clash with the EJC core subunit Y14. A CWC27/CWC22 building block might thus form an intermediate landing platform for eIF4A3 onto the B^act complex prior to its conversion into C complex. Knock-down of either CWC27 or CWC22 in immortalized retinal pigment epithelial cells affects numerous common genes, indicating that these proteins cooperate, targeting the same pathways. As the most up-regulated genes encode factors involved in inflammation, our findings suggest a possible link to the retinal degeneration associated with CWC27 deficiencies.

psichomics: graphical application for alternative splicing quantification and analysis

Alternative pre-mRNA splicing generates functionally distinct transcripts from the same gene and is involved in the control of multiple cellular processes, with its dysregulation being associated with a variety of pathologies. The advent of next-generation sequencing has enabled global studies of alternative splicing in different physiological and disease contexts. However, current bioinformatics tools for alternative splicing analysis from RNA-seq data are not user-friendly, disregard available exon-exon junction quantification or have limited downstream analysis features. To overcome such limitations, we have developed psichomics, an R package with an intuitive graphical interface for alternative splicing quantification and downstream dimensionality reduction, differential splicing and gene expression and survival analyses based on The Cancer Genome Atlas, the Genotype-Tissue Expression project, the Sequence Read Archive project and user-provided data. These integrative analyses can also incorporate clinical and molecular sample-associated features. We successfully used psichomics in a laptop to reveal alternative splicing signatures specific to stage I breast cancer and associated novel putative prognostic factors.

Efficient and accurate detection of splice junctions from RNA-seq with Portcullis

Next-generation sequencing technologies enable rapid and cheap genome-wide transcriptome analysis, providing vital information about gene structure, transcript expression, and alternative splicing. Key to this is the accurate identification of exon-exon junctions from RNA sequenced (RNA-seq) reads. A number of RNA-seq aligners capable of splitting reads across these splice junctions (SJs) have been developed; however, it has been shown that while they correctly identify most genuine SJs available in a given sample, they also often produce large numbers of incorrect SJs. Here, we describe the extent of this problem using popular RNA-seq mapping tools and present a new method, called Portcullis, to rapidly filter false SJs derived from spliced alignments. We show that Portcullis distinguishes between genuine and false-positive junctions to a high degree of accuracy across different species, samples, expression levels, error profiles, and read lengths. Portcullis is portable, efficient, and, to our knowledge, currently the only SJ prediction tool that reliably scales for use with large RNA-seq datasets and large, highly fragmented genomes, while delivering accurate SJs.

Gene Regulatory Network Perturbation by Genetic and Epigenetic Variation

Gene regulatory networks underlie biological function and cellular physiology. Alternative splicing (AS) is a fundamental step in gene regulatory networks and plays a key role in development and disease. In addition to the identification of aberrant AS events, an increasing number of studies are focusing on molecular determinants of AS, including genetic and epigenetic regulators. We review here recent efforts to identify various deregulated AS events as well as their molecular determinants that alter biological functions, and discuss clinical features of AS and their druggable potential.

Global isoform-specific transcript alterations and deregulated networks in clear cell renal cell carcinoma

Extensive genome-wide analyses of deregulated gene expression have now been performed for many types of cancer. However, most studies have focused on deregulation at the gene-level, which may overlook the alterations of specific transcripts for a given gene. Clear cell renal cell carcinoma (ccRCC) is one of the best-characterized and most pervasive renal cancers, and ccRCCs are well-documented to have aberrant RNA processing. In the present study, we examine the extent of aberrant isoform-specific RNA expression by reporting a comprehensive transcript-level analysis, using the new kallisto-sleuth-RATs pipeline, investigating coding and non-coding differential transcript expression in ccRCC. We analyzed 50 ccRCC tumors and their matched normal samples from The Cancer Genome Altas datasets. We identified 7,339 differentially expressed transcripts and 94 genes exhibiting differential transcript isoform usage in ccRCC. Additionally, transcript-level coexpression network analyses identified vasculature development and the tricarboxylic acid cycle as the most significantly deregulated networks correlating with ccRCC progression. These analyses uncovered several uncharacterized transcripts, including lncRNAs FGD5-AS1 and AL035661.1, as potential regulators of the tricarboxylic acid cycle associated with ccRCC progression. As ccRCC still presents treatment challenges, our results provide a new resource of potential therapeutics targets and highlight the importance of exploring alternative methodologies in transcriptome-wide studies.

Comparative Analysis and Classification of Cassette Exons and Constitutive Exons

Alternative splicing (AS) is a major engine that drives proteome diversity in mammalian genomes and is a widespread cause of human hereditary diseases. More than 95% of genes in the human genome are alternatively spliced, and the most common type of AS is the cassette exon. Recent discoveries have demonstrated that the cassette exon plays an important role in genetic diseases. To discover the formation mechanism of cassette exon events, we statistically analyze cassette exons and find that cassette exon events are strongly influenced by individual exons that are smaller in size and that have a lower GC content, more codon terminations, and weaker splice sites. We propose an improved random-forest-based hybrid method of distinguishing cassette exons from constitutive exons. Our method achieves a high accuracy in classifying cassette exons and constitutive exons and is verified to outperform previous approaches. It is anticipated that this study will facilitate a better understanding of the underlying mechanisms in cassette exons.

paraGSEA: a scalable approach for large-scale gene expression profiling

More studies have been conducted using gene expression similarity to identify functional connections among genes, diseases and drugs. Gene Set Enrichment Analysis (GSEA) is a powerful analytical method for interpreting gene expression data. However, due to its enormous computational overhead in the estimation of significance level step and multiple hypothesis testing step, the computation scalability and efficiency are poor on large-scale datasets. We proposed paraGSEA for efficient large-scale transcriptome data analysis. By optimization, the overall time complexity of paraGSEA is reduced from O(mn) to O(m+n), where m is the length of the gene sets and n is the length of the gene expression profiles, which contributes more than 100-fold increase in performance compared with other popular GSEA implementations such as GSEA-P, SAM-GS and GSEA2. By further parallelization, a near-linear speed-up is gained on both workstations and clusters in an efficient manner with high scalability and performance on large-scale datasets. The analysis time of whole LINCS phase I dataset (GSE92742) was reduced to nearly half hour on a 1000 node cluster on Tianhe-2, or within 120 hours on a 96-core workstation. The source code of paraGSEA is licensed under the GPLv3 and available at http://github.com/ysycloud/paraGSEA.