With fasta/fastq input from an RNA-seq experiment, SNPs are found by KisSplice without using a reference. KisSplice provides only a local context around the SNP, but a reference transcriptome can be built from the RNAseq data using a full-lenth transcriptome assembler like Trinity. Then SNPs predicted by KisSplice can be positionned along the transcripts (with BLAT). Some SNPs that do not map on the transcripts of Trinity, called orphan SNPs, are harder to study but can still be of interest. We propose a method, KisSplice2RefTranscriptome, to predict a functional impact for the positioned SNPs, and intersect these results with condition-specific SNPs. Overall, starting from RNAseq data only, we obtain a list of condition-specific SNPs stratified by functional impact.
Version 1.3.1KisSplice2RefTranscriptome version 1.3.1 (2017-07-18)
- Improvements / Bugs fixed
- We know perform a test to check that the value given for parameter -K and the k-mer size used to run KisSplice are identical.
RequirementsKisSplice2RefTranscritpome is written in Python and is running on Linux platforms. A version of python 2.7 or more is needed.
How to get my input files ?
kissplice -s 1 -k 41 -r sample_example/reads1.fq -r sample_example/reads2.fq -r sample_example/readsN.fq
Trinity ORFs fileHere is a basic example command in order to run Trinity :
trinity --seqType fq --left reads_lefts.fq --right reads_right.fqThe assembled transcripts will be written in the Trinity.fasta file. You will then need to run Transdecoder (included in the Trinity package). Transdecoder will predict the ORFs of the transcripts that will be used for the prediction of the functional impact of SNPs. Here is a basic command in order to run Trandecoder on Trinity ouput :
Trandecoder.LongOrfs -t Trinity.fastaTrandecoder will output several files, KisSplice2RefTranscritpome only needs the .bed file.
Running BLATKisSplice2RefTranscriptome needs the bubbles to be postioned on the reference transcriptome. For now, only psl formated files are supported by K2rt. You can run BLAT to get this file. Here is a command in order to run BLAT :
blat --minIdentity=80 Trinity.fasta kissplice_output_type0_coherent.fa output.psl
Running KissDE (optional)You can intersect KisSplice2RefTranscriptomeresults with condition-specific SNPs by giving to K2rt the ouput of kissDE.
Using KisSplice2RefTranscriptomeYou have access to all options description with the help command :
Usage: kissplice2reftranscriptome [-h]
[-Q QC_MIN] [-C KC_MIN] [-K K_VAL]