awesome-alternative-splicing

March 31, 2018 ยท View on GitHub

What is this?

This is a resource for software and information about alternative splicing. Contributions welcome...

Types of alternative splicing

  • Skipped exon or cassette exon (SE): An exon can be retained or spliced out of the primary transcript.
  • Mutually exclusion exons (MXE): One of two exons is retained in mRNA splicing, but not both.
  • Alternative 5' splice site (A5SS): An alternative 5' splice junction (donor site) is used that changes the 3' boundary of the upstream exon.
  • Alternative 3' splice site (A3SS): An alternative 3' splice junction (acceptor site) is used that changes the 5' boundary of the downstream exon.
  • Retained intron (RI): A sequence may be spliced out as an intron or retained, and there are no flanking introns.

Types of alternative splicing

Source: http://rnaseq-mats.sourceforge.net/rmats3.0.9/

Software

  • ALEXA-Seq - alternative expression analysis by massively parallel sequencing.
  • AltAnalyze - analyze alternative splicing from single-cell and RNA-Seq data.
  • Cufflinks - assemble and quantify transcripts.
  • DEXSeq - identify differential exon usage.
  • flotilla - reproduce machine learning analysis of gene expression and alternative splicing data.
  • GMAP and GSNAP - detect complex variants and splicing in short reads, SNP-tolerant.
  • G-Mo.R-Se - maps splice junctions to genome.
  • HMMSplicer - discovery canonical and non-canonical splice junctions in short read datasets.
  • JunctionSeq - identify differential splice junctions.
  • MapSplice - map RNA-seq data to reference genome for splice junction discovery.
  • MISO - determine alternative splicing expression.
  • MMES - statistically determine alternative splicing.
  • outrigger - calculate alternative splicing scores of RNA-Seq data based on junction reads and a de novo, custom annotation created with a graph database, especially made for single-cell analyses.
  • rMATS - RNA-Seq Multavariate Analysis of Transcript Splicing. Reading rMATS output
  • rmats2sashimiplot - visualize rMATS output using sashimi plots.
  • SAW - identify splicing events from RNA-Seq data.
  • Scripture - reconstruct transcript isoforms.
  • SingleSplice - detect biological variation in alternative splicing within a population of single cells.
  • SpliceMap - discover and align splice junctions for RNA-Seq reads.
  • SpliceR - detect alternative splicing and predict coding potential.
  • SplicingCompass - detect differential splicing using RNA-Seq data.
  • SplitSeek - predict splice events from RNA-Seq data.
  • STAR - identify alternative splicing.
  • SUPPA - identify alternative splicing.
  • TopHat - map splice junctions for RNA-Seq reads.

Databases

  • ASIP - Alternative Splicing in Plants.
  • ASG - Alternative Splicing Gallery for human genes.
  • ASPicDB - Alternative Splicing PredICtion DataBase.
  • ENSEMBL - Human and mouse genome annotations.
  • FAST DB/Easana - Friendly Alternative Splicing and Transcripts Database.
  • Hollywood exon annotation database - A website for querying a relational database of constitutive and alternative human exons, by using biological and descriptive features.
  • HS3D - Data set of Homo Sapiens Exon, Intron and Splice regions extracted from GenBank Rel.123.
  • H-DBAS - Human-transcriptome DataBase for Alternative Splicing.
  • MAASE - Convenient access, identification, and annotation of alternative splicing events (ASEs), designed specifically with experimentalists in mind.
  • Pro-Splicer - Alternative splicing database based on protein, mRNA, and EST Sequences.
  • SpliceNest - Visualizing splicing of genes from EST Data for human, mouse, Drodophila and Arabidopsis.

Review of RNA-Seq splicing tools

ToolPerforms split-read alignmentTranscript reconstruction (assembly)Expression Analysis (any)Gene expression analysisTranscript specific expression analysisExon junction expressionQuantitative alternative expression analysisExpression level sensitivityOutputMinimum read length required or recommendedVisualization toolPerforms comparisons between conditions (ex. tumor vs normal)Relevant comparison to Alexa-seqData type supportedCitation
Alexa-SeqNNYYYYYJunction"Expression and structure information for junctions and genes; UCSC track infoextensive alternative expression visualization / statistics / graphs"No minimum (tested on 36bp-100bp reads)"Extensiveincludes custom graphs and links to UCSC browser"Y-
CufflinksNYYYYNYTranscript"Transcript information and expression statisticsBEDGTF"75bpUCSC browserYY
ScriptureNYYNN*NNTranscript*"Transcript structure information and non-parsimonious expression statisticsBED"75bpUCSC browserNNIllumina
SpliceMapNNYYNYNJunction"Alignments. SAMBEDWig"50bpUCSC browserNN
TopHatYNNNNYNJunction"Alignments. SAMBEDWig"75bpNNN
MMESYNNNNYNJunctionIdentified splice junctions and p-values25bpNNNIllumina"Wang et al.
G-Mo.R-SeYYNNNNNTranscript*Transcript structure information. GFF25bpGrape Genome BrowserNNIllumina"Denoeud et al.
SplitSeekYNNNNYNJunctionAlignments. BED50bpUCSC browserNNSOLiD only"Ameur et al.2010
GSNAPYYNNNNNN/AAlignments. SAM and FASTAMinimum 14bp (tested on 36bp reads)UCSC browserNNIlluminasodium bisulfite-treated DNA sequencing (for analysis of methylation status)
ToolPerforms split-read alignmentTranscript reconstruction (assembly)Expression Analysis (any)Gene expression analysisTranscript specific expression analysisExon or Junction expressionQuantitative alternative expression analysisExpression level sensitivityOutputMinimum read length required or recommendedAbility to identify rearrangements / indelsJunction IdentificationImplementationPublic toolOpen Source
Alexa-seqNNYYYYYJunction"Expression and structure information for junctions and genes; UCSC track infoextensive alternative expression visualization / statistics / graphs"No minimum (tested on 36bp-100bp reads)NDatabasePerl/R/UnixY
CufflinksNYYYYNYTranscript"Transcript information and expression statisticsBEDGTF"75bpNPredicted from dataC++
ScriptureNYYNN*NNTranscript*"Transcript structure information and non-parsimonious expression statisticsBED"75bpNPredicted from dataJavaY
SpliceMapNNYYNYNJunction"Alignments. SAMBEDWig"50bpNPredicted from dataPython
TopHatYNNNNYNJunction"Alignments. SAMBEDWig"75bpNPredicted from dataC++/Python
MMESYNNNNYNJunctionIdentified splice junctions and p-values25bpNPredicted from dataPublished algorithm onlyNN/A
G-Mo.R-SeYYNNNNNTranscript*Transcript structure information. GFF25bpNPredicted from dataPerlY**Y
SplitSeekYNNNNYNJunctionAlignments. BED50bpYPredicted from dataPerlYY (GPL)
GSNAPYYNNNNNN/AAlignments. SAM and FASTAMinimum 14bp (tested on 36bp reads)YDatabase or predicted form data"Source code in Cutility programs in Perl"Y
ToolPerforms split-read alignmentTranscript reconstruction (assembly)Expression Analysis (any)Gene expression analysisTranscript specific expression analysisExon or Junction expressionQuantitative alternative expression analysisExpression level sensitivityOutputMinimum read length required or recommendedAbility to identify rearrangements / indelsJunction IdentificationImplementationPublic toolOpen Source
Alexa-seqNNYYYYYJunction"Expression and structure information for junctions and genes; UCSC track infoextensive alternative expression visualization / statistics / graphs"No minimum (tested on 36bp-100bp reads)NDatabasePerl/R/UnixY
CufflinksNYYYYNYTranscript"Transcript information and expression statisticsBEDGTF"75bpNPredicted from dataC++
ScriptureNYYNN*NNTranscript*"Transcript structure information and non-parsimonious expression statisticsBED"75bpNPredicted from dataJavaY
SpliceMapNNYYNYNJunction"Alignments. SAMBEDWig"50bpNPredicted from dataPython
TopHatYNNNNYNJunction"Alignments. SAMBEDWig"75bpNPredicted from dataC++/Python
MMESYNNNNYNJunctionIdentified splice junctions and p-values25bpNPredicted from dataPublished algorithm onlyNN/A
G-Mo.R-SeYYNNNNNTranscript*Transcript structure information. GFF25bpNPredicted from dataPerlY**Y
SplitSeekYNNNNYNJunctionAlignments. BED50bpYPredicted from dataPerlYY (GPL)
GSNAPYYNNNNNN/AAlignments. SAM and FASTAMinimum 14bp (tested on 36bp reads)YDatabase or predicted form data"Source code in Cutility programs in Perl"Y

* Transcript expression values are non-parsimonious

** Not supported (current version not stable)