Population calling | joint calling | ensemble calling

[amizeranschi]

Description of feature

Hi!

I was trying out the joint_germline option, which produces a multi-sample VCF file from haplotypecaller. It would be helpful to have a similar feature for other callers like deepvariant, freebayes and strelka2.

It could also nice to have support for ensemble calling (from multiple callers), for some use cases where the number of samples is manageable.

One way to achieve this would be to use bcbio.variation.recall, which was created by Brad Chapman and used in bcbio-nextgen.

Relevant documention from bcbio-nextgen can be found here:

• population calling
• ensemble calling

For joint calling, another solution could be to have the variant callers output gVCF files, and then use GATK's GenotypeGVCFs or something similar on those. There was a discussion about this on the deepvariant GitHub page a few years ago: google/deepvariant#142 (comment).

[amizeranschi]

Past relevant feature request that I just noticed, but also ties into this:
#453

[FriederikeHanssen]

• Joint calling with mpileup  #1124
• FreeBayes
• Deepvariant
• Strelka2

[SofiaMAhmed]

Hi Thanks for this really useful tool. I used it recently but the joint calling option failed for freebayes. Have others reported this issue?

[FriederikeHanssen]

hi! There is no joint calling currently implemented for freebayes.

[SofiaMAhmed]

Thanks for the speedy response. Is this something that may be implemented in the near future? We are really keen to incorporate this pipeline for our routine work. My current work around is to use the sarek output bams then mark readgroups and run freebayes seperatley, which seemed like the best option, other than merging the outputted vcf's and recalling against the bams.
Cheers Soph

[FriederikeHanssen]

I don't think we have resources at the moment on the developer side. If you'd like to take a stab at it, we would be very happy to guide you through though :)

[amizeranschi]

I was looking into this a while ago, but couldn't figure out any practical way to do it. The bcbio-variation-recall package I mentioned above has long been discontinued. Outside of that, there was GLNexus which looked promising, but that has become discontinued as well a couple of years ago and there are currently some serious issues like this that were reported but never solved.