SpliceNet: recovering splicing isoform-specific differential gene networks from RNA-Seq data of normal and diseased samples

Conventionally, overall gene expressions from microarrays are used to infer gene networks, but it is challenging to account splicing isoforms. High-throughput RNA Sequencing has made splice variant profiling practical. However, its true merit in quantifying splicing isoforms and isoform-specific exon expressions is not well explored in inferring gene networks. This study demonstrates SpliceNet, a method to infer isoform-specific co-expression networks from exon-level RNA-Seq data, using large dimensional trace. It goes beyond differentially expressed genes and infers splicing isoform network changes between normal and diseased samples. It eases the sample size bottleneck; evaluations on simulated data and lung cancer-specific ERBB2 and MAPK signaling pathways, with varying number of samples, evince the merit in handling high exon to sample size ratio datasets. Inferred network rewiring of well established Bcl-x and EGFR centered networks from lung adenocarcinoma expression data is in good agreement with literature. Gene level evaluations demonstrate a substantial performance of SpliceNet over canonical correlation analysis, a method that is currently applied to exon level RNA-Seq data. SpliceNet can also be applied to exon array data. SpliceNet is distributed as an R package available at http://www.jjwanglab.org/SpliceNet.     

Alternative Splicing in the Differentiation of Human Embryonic Stem Cells into Cardiac Precursors

The role of alternative splicing in self-renewal, pluripotency and tissue lineage specification of human embryonic stem cells (hESCs) is largely unknown. To better define these regulatory cues, we modified the H9 hESC line to allow selection of pluripotent hESCs by neomycin resistance and cardiac progenitors by puromycin resistance. Exon-level microarray expression data from undifferentiated hESCs and cardiac and neural precursors were used to identify splice isoforms with cardiac-restricted or common cardiac/neural differentiation expression patterns. Splice events for these groups corresponded to the pathways of cytoskeletal remodeling, RNA splicing, muscle specification, and cell cycle checkpoint control as well as genes with serine/threonine kinase and helicase activity. Using a new program named AltAnalyze (http://www.AltAnalyze.org), we identified novel changes in protein domain and microRNA binding site architecture that were predicted to affect protein function and expression. These included an enrichment of splice isoforms that oppose cell-cycle arrest in hESCs and that promote calcium signaling and cardiac development in cardiac precursors. By combining genome-wide predictions of alternative splicing with new functional annotations, our data suggest potential mechanisms that may influence lineage commitment and hESC maintenance at the level of specific splice isoforms and microRNA regulation.     

Unusual Regulation of Splicing of the Cholinergic Locus in Caenorhabditis elegans

The essential neurotransmitter acetylcholine functions throughout the animal kingdom. In Caenorhabditis elegans, the acetylcholine biosynthetic enzyme [choline acetyltransferase (ChAT)] and vesicular transporter [vesicular acetylcholine transporter (VAChT)] are encoded by the cha-1 and unc-17 genes, respectively. These two genes compose a single complex locus in which the unc-17 gene is nested within the first intron of cha-1, and the two gene products arise from a common pre-messenger RNA (pre-mRNA) by alternative splicing. This genomic organization, known as the cholinergic gene locus (CGL), is conserved throughout the animal kingdom, suggesting that the structure is important for the regulation and function of these genes. However, very little is known about CGL regulation in any species. We now report the identification of an unusual type of splicing regulation in the CGL of C. elegans, mediated by two pairs of complementary sequence elements within the locus. We show that both pairs of elements are required for efficient splicing to the distal acceptor, and we also demonstrate that proper distal splicing depends more on sequence complementarity within each pair of elements than on the sequences themselves. We propose that these sequence elements are able to form stem-loop structures in the pre-mRNA; such structures would favor specific splicing alternatives and thus regulate CGL splicing. We have identified complementary elements at comparable locations in the genomes of representative species of other animal phyla; we suggest that this unusual regulatory mechanism may be a general feature of CGLs.     

MutPred Splice: machine learning-based prediction of exonic variants that disrupt splicing

We have developed a novel machine-learning approach, MutPred Splice, for the identification of coding region substitutions that disrupt pre-mRNA splicing. Applying MutPred Splice to human disease-causing exonic mutations suggests that 16% of mutations causing inherited disease and 10 to 14% of somatic mutations in cancer may disrupt pre-mRNA splicing. For inherited disease, the main mechanism responsible for the splicing defect is splice site loss, whereas for cancer the predominant mechanism of splicing disruption is predicted to be exon skipping via loss of exonic splicing enhancers or gain of exonic splicing silencer elements. MutPred Splice is available at http://mutdb.org/mutpredsplice.     

Using high-density DNA methylation arrays to profile copy number alterations

The integration of genomic and epigenomic data is an increasingly popular approach for studying the complex mechanisms driving cancer development. We have developed a method for evaluating both methylation and copy number from high-density DNA methylation arrays. Comparing copy number data from Infinium HumanMethylation450 BeadChips and SNP arrays, we demonstrate that Infinium arrays detect copy number alterations with the sensitivity of SNP platforms. These results show that high-density methylation arrays provide a robust and economic platform for detecting copy number and methylation changes in a single experiment. Our method is available in the ChAMP Bioconductor package: http://www.bioconductor.org/packages/2.13/bioc/html/ChAMP.html.     

RNAmotifs: prediction of multivalent RNA motifs that control alternative splicing

RNA-binding proteins (RBPs) regulate splicing according to position-dependent principles, which can be exploited for analysis of regulatory motifs. Here we present RNAmotifs, a method that evaluates the sequence around differentially regulated alternative exons to identify clusters of short and degenerate sequences, referred to as multivalent RNA motifs. We show that diverse RBPs share basic positional principles, but differ in their propensity to enhance or repress exon inclusion. We assess exons differentially spliced between brain and heart, identifying known and new regulatory motifs, and predict the expression pattern of RBPs that bind these motifs. RNAmotifs is available at https://bitbucket.org/rogrro/rna_motifs.     

FineSplice,enhanced splice junction detection and quantification: a novel pipeline based on the assessment of diverse RNA-Seq alignment solutions

Alternative splicing is the main mechanism governing protein diversity. The recent developments in RNA-Seq technology have enabled the study of the global impact and regulation of this biological process. However, the lack of standardized protocols constitutes a major bottleneck in the analysis of alternative splicing. This is particularly important for the identification of exon–exon junctions, which is a critical step in any analysis workflow. Here we performed a systematic benchmarking of alignment tools to dissect the impact of design and method on the mapping, detection and quantification of splice junctions from multi-exon reads. Accordingly, we devised a novel pipeline based on TopHat2 combined with a splice junction detection algorithm, which we have named FineSplice. FineSplice allows effective elimination of spurious junction hits arising from artefactual alignments, achieving up to 99% precision in both real and simulated data sets and yielding superior F₁ scores under most tested conditions. The proposed strategy conjugates an efficient mapping solution with a semi-supervised anomaly detection scheme to filter out false positives and allows reliable estimation of expressed junctions from the alignment output. Ultimately this provides more accurate information to identify meaningful splicing patterns. FineSplice is freely available at https://sourceforge.net/p/finesplice/.     

Interactome analysis brings splicing into focus

The spliceosome is a huge molecular machine that assembles dynamically onto its pre-mRNA substrates. A new study based on interactome analysis provides clues about how splicing-regulatory proteins modulate assembly of the spliceosome to either activate or repress splicing.Please see related Research article: http://www.genomebiology.com/2015/16/1/119/abstract     

Coordination of Recombination with Meiotic Progression in the Caenorhabditis elegans Germline by KIN-18, a TAO Kinase That Regulates the Timing of MPK-1 Signaling

Meiosis is a tightly regulated process requiring coordination of diverse events. A conserved ERK/MAPK-signaling cascade plays an essential role in the regulation of meiotic progression. The Thousand And One kinase (TAO) kinase is a MAPK kinase kinase, the meiotic role of which is unknown. We have analyzed the meiotic functions of KIN-18, the homolog of mammalian TAO kinases, in Caenorhabditis elegans. We found that KIN-18 is essential for normal meiotic progression; mutants exhibit accelerated meiotic recombination as detected both by analysis of recombination intermediates and by crossover outcome. In addition, ectopic germ-cell differentiation and enhanced levels of apoptosis were observed in kin-18 mutants. These defects correlate with ectopic activation of MPK-1 that includes premature, missing, and reoccurring MPK-1 activation. Late progression defects in kin-18 mutants are suppressed by inhibiting an upstream activator of MPK-1 signaling, KSR-2. However, the acceleration of recombination events observed in kin-18 mutants is largely MPK-1-independent. Our data suggest that KIN-18 coordinates meiotic progression by modulating the timing of MPK-1 activation and the progression of recombination events. The regulation of the timing of MPK-1 activation ensures the proper timing of apoptosis and is required for the formation of functional oocytes. Meiosis is a conserved process; thus, revealing that KIN-18 is a novel regulator of meiotic progression in C. elegans would help to elucidate TAO kinase’s role in germline development in higher eukaryotes.     

Systematic Analyses of rpm-1 Suppressors Reveal Roles for ESS-2 in mRNA Splicing in Caenorhabditis elegans

The PHR (Pam/Highwire/RPM-1) family of ubiquitin E3 ligases plays conserved roles in axon patterning and synaptic development. Genetic modifier analysis has greatly aided the discovery of the signal transduction cascades regulated by these proteins. In Caenorhabditis elegans, loss of function in rpm-1 causes axon overgrowth and aberrant presynaptic morphology, yet the mutant animals exhibit little behavioral deficits. Strikingly, rpm-1 mutations strongly synergize with loss of function in the presynaptic active zone assembly factors, syd-1 and syd-2, resulting in severe locomotor deficits. Here, we provide ultrastructural evidence that double mutants, between rpm-1 and syd-1 or syd-2, dramatically impair synapse formation. Taking advantage of the synthetic locomotor defects to select for genetic suppressors, previous studies have identified the DLK-1 MAP kinase cascade negatively regulated by RPM-1. We now report a comprehensive analysis of a large number of suppressor mutations of this screen. Our results highlight the functional specificity of the DLK-1 cascade in synaptogenesis. We also identified two previously uncharacterized genes. One encodes a novel protein, SUPR-1, that acts cell autonomously to antagonize RPM-1. The other affects a conserved protein ESS-2, the homolog of human ES2 or DGCR14. Loss of function in ess-2 suppresses rpm-1 only in the presence of a dlk-1 splice acceptor mutation. We show that ESS-2 acts to promote accurate mRNA splicing when the splice site is compromised. The human DGCR14/ES2 resides in a deleted chromosomal region implicated in DiGeorge syndrome, and its mutation has shown high probability as a risk factor for schizophrenia. Our findings provide the first functional evidence that this family of proteins regulate mRNA splicing in a context-specific manner.     

The Role of Exo1p Exonuclease in DNA End Resection to Generate Gene Conversion Tracts in Saccharomyces cerevisiae

The yeast Exo1p nuclease functions in multiple cellular roles: resection of DNA ends generated during recombination, telomere stability, DNA mismatch repair, and expansion of gaps formed during the repair of UV-induced DNA damage. In this study, we performed high-resolution mapping of spontaneous and UV-induced recombination events between homologs in exo1 strains, comparing the results with spontaneous and UV-induced recombination events in wild-type strains. One important comparison was the lengths of gene conversion tracts. Gene conversion events are usually interpreted as reflecting heteroduplex formation between interacting DNA molecules, followed by repair of mismatches within the heteroduplex. In most models of recombination, the length of the gene conversion tract is a function of the length of single-stranded DNA generated by end resection. Since the Exo1p has an important role in end resection, a reduction in the lengths of gene conversion tracts in exo1 strains was expected. In accordance with this expectation, gene conversion tract lengths associated with spontaneous crossovers in exo1 strains were reduced about twofold relative to wild type. For UV-induced events, conversion tract lengths associated with crossovers were also shorter for the exo1 strain than for the wild-type strain (3.2 and 7.6 kb, respectively). Unexpectedly, however, the lengths of conversion tracts that were unassociated with crossovers were longer in the exo1 strain than in the wild-type strain (6.2 and 4.8 kb, respectively). Alternative models of recombination in which the lengths of conversion tracts are determined by break-induced replication or oversynthesis during strand invasion are proposed to account for these observations.     

Remarkably Divergent Regions Punctuate the Genome Assembly of the Caenorhabditis elegans Hawaiian Strain CB4856

The Hawaiian strain (CB4856) of Caenorhabditis elegans is one of the most divergent from the canonical laboratory strain N2 and has been widely used in developmental, population, and evolutionary studies. To enhance the utility of the strain, we have generated a draft sequence of the CB4856 genome, exploiting a variety of resources and strategies. When compared against the N2 reference, the CB4856 genome has 327,050 single nucleotide variants (SNVs) and 79,529 insertion–deletion events that result in a total of 3.3 Mb of N2 sequence missing from CB4856 and 1.4 Mb of sequence present in CB4856 but not present in N2. As previously reported, the density of SNVs varies along the chromosomes, with the arms of chromosomes showing greater average variation than the centers. In addition, we find 61 regions totaling 2.8 Mb, distributed across all six chromosomes, which have a greatly elevated SNV density, ranging from 2 to 16% SNVs. A survey of other wild isolates show that the two alternative haplotypes for each region are widely distributed, suggesting they have been maintained by balancing selection over long evolutionary times. These divergent regions contain an abundance of genes from large rapidly evolving families encoding F-box, MATH, BATH, seven-transmembrane G-coupled receptors, and nuclear hormone receptors, suggesting that they provide selective advantages in natural environments. The draft sequence makes available a comprehensive catalog of sequence differences between the CB4856 and N2 strains that will facilitate the molecular dissection of their phenotypic differences. Our work also emphasizes the importance of going beyond simple alignment of reads to a reference genome when assessing differences between genomes.