This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to calling variants in RNAseq data from multiple samples?

Hi everybady! I would like to calling variants in RNSeq multiple samples. Let's say that I have two different group and I would like to know the differences among them. Each study group is formed by ~ bam files. Is the first time that I do variant calling, I've also read some post as Workflow Or Tutorial For Snp Calling? with a very clear workflow. I would like to know if I have to merge all the bam files for each group in a unique sam file after mark the duplicates. It should be correct? or should I use other specific pipeline for this analysis?

Thank you very very much in advance!

rna-seq gatk vc variant calling
• 4,240 views • link
updated by Samuel Brady Samuel Brady • written by Lila M Lila M

2 answers

See the GATK "best practise" for RNA Seq data here. Fair warning, that you're trying to do something that the data was never intended for. If you have SNPs in mind that you want to look at, I'd strongly recommend just looking at the alignments in IGV and eyeballing them. This protocol is a lot of work to get formal variant calls that won't offer you much more than eyeballing. Also if you're looking to find post transcriptional modifiers, then identifying them against noise is extremely difficult.

• 0 views • link
written by andrew.j.skelton73 andrew.j.skelton73

So the only way is to call variant in each sample separately and then eyeballing one by one? Does not exist other way to do variant calling in a whole population? Neither with other program?

Thank you!

• 0 views • link
written by Lila M Lila M

I'm not saying to call variants and then eyeball, I'm saying to look at your RNA Seq alignments in IGV and see in whatever bit you're interested in if something lights up.

You apply tools like the haplotype caller on a sample by sample basis, however if you have a decent number of samples you can do joint calling, although I don't recall that being benchmarked for RNA seq data.

If you're asking "What's my causal variant based on al my samples", that depends on your hypothesis. There are tools such as VEP to predict the effect of variants. If you want to do mendelian inheritance or something along those lines, take a look at GEMINI, MendelMD, or take a look at the tool list available at omicstools

• 0 views • link
written by andrew.j.skelton73 andrew.j.skelton73

and what about call variants independently and the merge the vcf file?

Thanks for all the information!

• 0 views • link
written by Lila M Lila M

I was thinking on it. So this approach should work, after following the GATK "best practise" for RNA Seq data, right?

• 0 views • link
written by Lila M Lila M

If you make a VCF for each run of the Haplotype Caller (i.e. each sample), then you can combine them to a multisample VCF using VCF-tools

• 0 views • link
written by andrew.j.skelton73 andrew.j.skelton73

yes, that's the idea! Thanks for the confirmation!

• 0 views • link
written by Lila M Lila M

We have successfully used UNCeqR to call mutations in RNA-Seq data. The first author if this paper is very responsive to questions if you need help with it. The github page is also very helpful.

• 0 views • link
written by Samuel Brady Samuel Brady

Log in to answer this question.

More posts like this