Somatic mutation detection in HLA regions using hg38 alt-aware alignments
Answered
Hello,
I'm interested in identifying somatic mutations (using matched samples) both in the non-HLA and HLA regions.
For the HLA regions, one approach I've used, was to utilize known HLA-typing information about my samples to align to the specific HLA-alleles and then adapting the Mutect2 workflow. My post https://gatk.broadinstitute.org/hc/en-us/community/posts/4417333390619-ApplyBQSR-for-somatic-mutation-detection-in-HLA-regions referred to this approach.
Now, I would instead like to try a different approach, leveraging hg38 and the alt-aware alignments. I found the following (non-production) workflow: https://gatk.broadinstitute.org/hc/en-us/articles/360037498992#6, which has indication on how to call variants on HLA using alt-contigs. I'm wondering if anybody can point me to or can share an analogous workflow, but specific to somatic mutation calls using paired (tumor/normal) samples, aimed at HLA somatic mutation detection.
Thank you.
-
Hi Elisabetta Manduchi, I'm not aware of anything to point you to for using alt contigs to call somatic variants in the HLA, but a good method along the lines of your first approach is the Polysolver pipeline developed by Sachet Shukla, our friend at MD Anderson Cancer Center.
-
Thank you David Benjamin. Indeed my first approach was along the lines of Polysolver. I've also tried to use the Polysolver docker itself as is, which however did not yield any results in my case.
-
Where is polysolver ? the link seems to be broken
http://software.broadinstitute.org/cancer/cga/polysolver -
I think you may be able to get it from Docker, maybe here https://hub.docker.com/layers/sachet/polysolver/v4/images/sha256-4ca4f341da48058c32061137371dd3adb4ea816750a4adfa2a7cc40b730ac6af
-
Thank you Elisabetta Manduchi for the response.
Please sign in to leave a comment.
5 comments