Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • rare cancer mutation detection pipeline

    hello,

    what is the best analysis pipeline for rare (<1%) cancer mutation detection from targeted (bait hybridization enrichment) resequencing (illumina)?

    GATK?

    and who has the best bait development service?
    has someone experiences with MYbaits from http://www.mycroarray.com/index.html?

  • #2
    Not a "pipeline" as such but Illumina has a couple of variant callers (they are open source but no support is provided) that you may want to check out.

    For germline mutations ISAAC variant caller: https://github.com/sequencing/isaac_variant_caller

    For somatic mutations "Strelka": https://github.com/genome-vendor/strelka

    Comment


    • #3
      Originally posted by GenoMax View Post
      Not a "pipeline" as such but Illumina has a couple of variant callers (they are open source but no support is provided) that you may want to check out.

      For germline mutations ISAAC variant caller: https://github.com/sequencing/isaac_variant_caller

      For somatic mutations "Strelka": https://github.com/genome-vendor/strelka
      For germline, why would you use this ISAAC caller as oppose to the typical GATK?

      Comment


      • #4
        Originally posted by ymc View Post
        For germline, why would you use this ISAAC caller as oppose to the typical GATK?
        dietmar13 already had mentioned GATK in original post. ISAAC caller happens to be very fast for single samples.

        Comment


        • #5
          Originally posted by dietmar13 View Post
          what is the best analysis pipeline for rare (<1%) cancer mutation detection from targeted (bait hybridization enrichment) resequencing (illumina)?

          GATK?
          I think this is the worst application you could attempt to use GATK for.

          MuTect on the other hand might be more useful. I presume this is what Foundation Medicine use for their hybridisation (SureSelect) based panel from the reading of the 'competing financial interests' section of this paper:



          VarScan2 comes out pretty well in this comparison:



          I'm just more interested how you determine things that are present at <1% aren't just noise?

          Comment


          • #6
            thank you all for your suggestions

            @bukowski,

            we look at liquid biopsies from cancer patients and want identify tumor DNA among a large number of non-cancer cells...

            I will look into MuTect and VarScan2 first...

            best,

            dietmar

            Comment


            • #7
              Originally posted by dietmar13 View Post
              @bukowski,

              we look at liquid biopsies from cancer patients and want identify tumor DNA among a large number of non-cancer cells...

              I will look into MuTect and VarScan2 first...

              best,

              dietmar
              I think there should be specialized software for circulating tumor dna. I am looking for it too.

              Comment


              • #8
                Originally posted by Bukowski View Post
                I'm just more interested how you determine things that are present at <1% aren't just noise?
                That's an excellent point. Once you get to that level, there is a lot of noise.

                From my experience with different somatic variant callers, it's very difficult to get decent calls at below 5%. Often when they say "low frequency", they mean at least 5%.

                I had some non-paired samples with coverage above 10,000X and even at that level, it was very difficult to get any standard variant caller to detect anything below 1%. I found LoFreq (http://sourceforge.net/p/lofreq/wiki/Home/) to be a tool that got the job done fairly well. At least it's designed specifically for ultra low frequency mutations. Even then, I was looking for known mutations, so false positives weren't much of an issue.
                Last edited by id0; 12-17-2013, 04:35 PM.

                Comment


                • #9
                  First: I hope you have sequenced deeply, because, obviously, to see variants with less than 1%, you need coverage in the hundreds or better thousands, for both the tumour and the control.

                  If so, deepSNV might be what you need:

                  M Gerstung et al.: Reliable detection of subclonal single-nucleotide variants in tumour cell populations
                  Nature Communications (2011) 3: 811 doi:10.1038/ncomms1814

                  Comment


                  • #10
                    We carry out paired-end sequencing of short fragments, so the two reads confirm each others basecalls. With that, we can detect mutations at 1 in 5000 level. It still takes deep sequencing, but just to sample the complex population of DNA sufficiently to detect the rare allele, not (as with single-end sequencing) to try to overcome the base level of sequencing error.
                    Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

                    Comment


                    • #11
                      Originally posted by SNPsaurus View Post
                      We carry out paired-end sequencing of short fragments, so the two reads confirm each others basecalls. With that, we can detect mutations at 1 in 5000 level. It still takes deep sequencing, but just to sample the complex population of DNA sufficiently to detect the rare allele, not (as with single-end sequencing) to try to overcome the base level of sequencing error.
                      Is that custom code or are there tools available to match variants on read pairs?

                      Comment


                      • #12
                        Any overlapper program will do the brunt of the work. There is quite a bit of custom code to take those results for this kind of purpose. Some other groups are doing the same thing more or less... http://www.biomedcentral.com/1471-2164/14/96
                        Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

                        Comment


                        • #13
                          I am sorry, but I want to ask the question:"Why do you want to detect extremely low abundant mutations?"

                          Comment


                          • #14
                            Code:
                            "Why do you want to detect extremely low abundant mutations?"
                            because i have a mixture of many non-cancer cells with (probably) only few cancer cells and want detect the latter using targeted resequencing...

                            Comment


                            • #15
                              Originally posted by dietmar13 View Post
                              because i have a mixture of many non-cancer cells with (probably) only few cancer cells and want detect the latter using targeted resequencing...
                              Thanks for the reply.

                              Do you see only very rare mutations? Or you see other relatively more abundant mutations that are unlikely to be germline mutations. If there are more abundant somatic mutation candidates, how do you consider them?

                              Woody

                              Comment

                              Latest Articles

                              Collapse

                              • seqadmin
                                Strategies for Sequencing Challenging Samples
                                by seqadmin


                                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                                03-22-2024, 06:39 AM
                              • seqadmin
                                Techniques and Challenges in Conservation Genomics
                                by seqadmin



                                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                                Avian Conservation
                                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                                03-08-2024, 10:41 AM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, Yesterday, 06:37 PM
                              0 responses
                              12 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, Yesterday, 06:07 PM
                              0 responses
                              10 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-22-2024, 10:03 AM
                              0 responses
                              51 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-21-2024, 07:32 AM
                              0 responses
                              68 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X