Question: best way to do joint genotyping
0
gravatar for bitpir
8 months ago by
bitpir90
bitpir90 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 • 520 views
ADD COMMENTlink modified 7 months ago by Biostar ♦♦ 20 • written 8 months ago by bitpir90
1

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 8 months ago by JJ350

great! Thanks for the tip!

ADD REPLYlink written 8 months ago by bitpir90

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 8 months ago by RamRS18k

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 8 months ago by bitpir90

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 8 months ago by RamRS18k
Please log in to add an answer.

Help
Access

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