Question: best way to do joint genotyping
gravatar for bitpir
12 weeks ago by
bitpir60 wrote:

I have exome sequences from 6 families (~38 samples) of European descent with similar hereditary diseases (e.g. epilepsy). Each family is supposedly unrelated. Some samples are from healthy family members, others are infected. What's the best way to do joint genotyping for this type of analysis? Should I (i) do joint genotyping on all 38 samples, then split them into individual vcf afterwards (for variant analysis), or (ii) do joint genotyping only for samples in the same family (so 6 separate joint-genotypings)? Also, let's say my samples comprise of different ethnicities, is it still better to do joint-genotyping or it won't really improve the variant calling in this case? Thanks for your help!

vcf gatk joint-genotyping • 274 views
ADD COMMENTlink modified 6 weeks ago by Biostar ♦♦ 20 • written 12 weeks ago by bitpir60

All of them together. However, beware that you might loose singletons when using joint genotyping. R packages like VariantAnnotation can handle multi-sample-vcfs - depends on the tool though.

ADD REPLYlink written 11 weeks ago by JJ240

great! Thanks for the tip!

ADD REPLYlink written 11 weeks ago by bitpir60

I'd join genotype every sample together, and retain them as a single vcf unless retaining them in one file doesn't make sense. Given that you don't need to separate them to analyze them separately, I don't see why they should be in separate VCFs.

ADD REPLYlink written 12 weeks ago by Ram14k

Hi Ram, I don't need to analyze each of them separately, but was planning to analyze on per family basis as some genes could be dominant/recessive/compound hets, etc. Thanks for the advice!

ADD REPLYlink written 11 weeks ago by bitpir60

No matter how samples are analyzed, it doesn't matter how many VCFs they are stored in. Calling heavily unrelated samples together can confound your results, but you're the best judge of how phenotypically alike you deem your samples.

ADD REPLYlink written 11 weeks ago by Ram14k
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: 948 users visited in the last hour