Question: Variant detection in RNA seq datasets without a reference using kISSPLICE.
gravatar for koushikponnanna
7 weeks ago by
koushikponnanna0 wrote:

I am presently working on RNA seq paired end data from few Drosophila tissues. As I dont have an established genome reference I need to rely on reference free variant calling tools and I find KISSPLICE suitable. I have outputs generated for my datasets. I have few questions related to the outputs generated.

  1. The paired end inputs as indicated in manual needs to be loaded separately as -r read.1 and -r read.2. A single output that is generated contains coverage as C1 and C2 for lower and upper path, do I have to sort the outputs to get information for my sample or the output that is generated is final output that I can rely on for that sample? Find the log file I have attached with this mail for one sample. I carried the same analysis by merging 2 fastq files into one file. I got slightly different outputs.

  2. Are there any graphical representations options that I can apply to the outputs? I have a total of 6 paired end samples . What is the best possible method to quantify and compare these samples? source of samples are testis tissue from different species.

thank you

kissplice snp rna-seq de novo • 276 views
ADD COMMENTlink modified 2 hours ago by vincent.lacroix50 • written 7 weeks ago by koushikponnanna0

from few Drosophila tissues. As I dont have an established genome reference

Am I missing something here, because I thought Drosophila was pretty well characterised, and there was a decent reference - Here for example.

ADD REPLYlink modified 7 weeks ago • written 7 weeks ago by andrew.j.skelton733.9k

I am aware of the well characterized reference for few species. But the problem is those are quiet diverged, and using them to call variants would not serve my purpose. That is the reason I am looking to call variants from raw reads.

ADD REPLYlink written 7 weeks ago by koushikponnanna0
gravatar for vincent.lacroix
2 hours ago by
vincent.lacroix50 wrote:

Dear Koushik Ponnanna, Thanks for your interest in KisSplice. Below are some answers to your questions: 1- When you load the reads from a pair separately, C1 will correspond to the read count of read 1, and C2 will correspond to the read count of read 2. In principle, C1+C2 should give you the read count obtained with the merged files. I think the difference you see is due to the fact that you used the raw data for the merged files, and the trimmed data for the separated files. 2- If you have replicates, you may want to use KissDE to find variants associated to a species, as described here: If you do not have replicates, you cannot use KissDE. But you can still run all other steps of our pipeline. This will produce an .xls file which you can exploit for your downstream analysis. For instance, you can select SNPs which affect the protein sequence. The final table contains counts, and estimated allele frequencies. You can then plot those in R, as in Supp Figure 6 of our NAR paper. If you need a visual output for specific genes, you will probably need to assemble the full-length transcripts with Trinity (or Oases), map your reads to the assembled transcripts, and visualise your transcripts of interest through IGV. Best regards, Vincent

ADD COMMENTlink written 2 hours ago by vincent.lacroix50
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1547 users visited in the last hour