Question: (Closed) GATK's Join Genotyping workflow on a large pedigree
gravatar for _r_am
5.9 years ago by
Baylor College of Medicine, Houston, TX
_r_am32k wrote:


I'm looking to use GATK's independent variant calling + joint genotyping workflow on a large pedigree with 10 different families. Is there any difference between the two approaches below:

  1. Run pipeline for all samples simultaneously
  2. Run pipeline for each family separately

I ask this because if the result of genotyping one family does not affect the results for other families, I might as well run them in parallel.

What do you say?


EDIT: I'm closing this question because I was mistaken in my assumption that GATK uses pedigree information during variant calling (HaplotypeCaller) or genotyping (GenotypeGVCFs).

pedigree joint genotyping gatk • 2.2k views
ADD COMMENTlink modified 5.9 years ago • written 5.9 years ago by _r_am32k
gravatar for Len Trigg
5.9 years ago by
Len Trigg1.5k
New Zealand
Len Trigg1.5k wrote:

I am no expert on GATK, but I believe that it does not incorporate pedigree information during the variant calling. Calling all the samples simultaneously will give somewhat better results for shared variants, although sensitivity to singletons may be reduced somewhat, so it depends a little on what your goals are. Both of these effects will be more pronounced in borderline cases (low coverage).

You might also want to consider using a variant caller such as Real Time Genomics that can directly incorporate pedigree information for multiple families during the calling (and give you a head start on de novo identification).


ADD COMMENTlink written 5.9 years ago by Len Trigg1.5k

Thank you for the inputs - I did not state my question clearly and I apologize for that. I was comparing running the entire pipeline (align+recalibrate+call+genotype) for all samples simultaneously versus per family. Does GATK incorporate pedigree information during genotyping?

And I've never heard of RTG - thank you for the pointer!

ADD REPLYlink written 5.9 years ago by _r_am32k

I guess I was mistaken - neither Variant Calling not Genotyping uses pedigree information

ADD REPLYlink written 5.9 years ago by _r_am32k
Please log in to add an answer.
The thread is closed. No new answers may be added.


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