Heads up:
We’re moving the GATK website, docs and forum to a new platform. Read the full story and breakdown of key changes on this blog.
Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
Attention:
We will be out of the office for a Broad Institute event from Dec 10th to Dec 11th 2019. We will be back to monitor the GATK forum on Dec 12th 2019. In the meantime we encourage you to help out other community members with their queries.
Thank you for your patience!

using MUTECT2 on tumor-only sample

BogdanBogdan Palo Alto, CAMember ✭✭

Dear all,

although this question has been asked a long time ago, if you do not mind asking it again, as I am looking for some updated workflows, strategies, ideas :

"what would be the acceptable pipeline for processing WES or WGS data from primary tumors that do not have matched germlines sequences ?"

thanks a lot,

-- bogdan

Tagged:

Best Answers

  • shleeshlee Cambridge ✭✭✭✭✭
    Accepted Answer

    Hi @Bogdan,

    For somatic SNV and indel calling, GATK4-Mutect2 (and FilterMutectCalls) is what we offer. See here for M2 documentation.

    As you know, SNVs and indels are only part of the tumor story. You should check out our somatic CNV calling pipeline to fill in the other half of the picture, e.g. with Tutorial#9143.

Answers

  • shleeshlee CambridgeMember, Broadie ✭✭✭✭✭
    Accepted Answer

    Hi @Bogdan,

    For somatic SNV and indel calling, GATK4-Mutect2 (and FilterMutectCalls) is what we offer. See here for M2 documentation.

    As you know, SNVs and indels are only part of the tumor story. You should check out our somatic CNV calling pipeline to fill in the other half of the picture, e.g. with Tutorial#9143.

  • BogdanBogdan Palo Alto, CAMember ✭✭
  • BogdanBogdan Palo Alto, CAMember ✭✭

    If I may add a question, please :

    after running MUTECT2 in TUMOR-ONLY mode, would it be a good practice to intersect the list of variants with 50-60 GERMLINE VCF files (that we have from other patients), in order to exclude as many germline variants as possible ?

  • shleeshlee CambridgeMember, Broadie ✭✭✭✭✭

    @Bogdan, presumably, your use of a PoN or known variants resources in your Mutect2 workflow will have already excluded germline variants.

  • BogdanBogdan Palo Alto, CAMember ✭✭

    Thank you. Please, when you have a minute, could I have the links to the webpages/information that describe the workflow, that describe the use of PoN and the exclusion of germline variants. Thanks ;)

  • shleeshlee CambridgeMember, Broadie ✭✭✭✭✭

    Hi @Bogdan,

    I've attached a PDF of the original GATK3-MuTect2 hands-on tutorial worksheet that I wrote back in January/February and that was presented by others in Leuven, Belgium (also back in February). Although the attached tutorial is about the previous implementation of M2, and does not use GATK4-M2, the concepts it discusses are still germane. The updated hands-on tutorial that uses GATK4-Mutect2 is written by another teammate. We can expect the team to post the recent UK workshop materials in the blog very shortly as the workshop ended yesterday.

  • BogdanBogdan Palo Alto, CAMember ✭✭

    many many thanks, the document is very very helpful ;)

  • BogdanBogdan Palo Alto, CAMember ✭✭

    Dear Lee, please, if I may add another question, namely where I could find a resource that would detail on how to do the selection of somatic variants that were obtained with MUTECT2 run in TUMOR-only mode. Thank you very much ..

  • BogdanBogdan Palo Alto, CAMember ✭✭

    many many thanks ;)

Sign In or Register to comment.