Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Ion Torrent Ampliseq: duplicate removal / coverage / variants?

    Hi Folks

    I would love to get some opinions regarding duplicate removal during variant calling for the AmpliSeq Comprehensive Cancer Panel.

    I have noted some previous discussions on the topic on this forum, but I am experiencing some additional issues.

    We have been doing tumor vs. normal samples, and do the realignment and base quality score recalibration steps (no duplicate removal) using GATK. We are currently testing a variety of somatic variant callers (Strelka, VarScan, MuTect, JointSNVMix, SomaticSniper).

    Ion Torrent states the following:

    "Marking duplicate reads is not appropriate for Ion AmpliSeq data, because many independent reads are expected to share the same 5' alignment position and 3' adapter flow as each other. Marking duplicates on an Ion AmpliSeq run risks inappropriately flagging many reads that are in fact independent of one another."

    When doing visual inspection of (possibly) variant positions in IGV, a lot of the variants seem to have been called based on read depth that sure look like duplicates upon inspection. There also seems to be a trend for variant positions to be at the edges of these reads .

    However, if I do duplicate removal, my average coverage drops from 300x to 9x.

    The duplicate reads vs. true coverage depth obviously has quite serious implications for all variant calling statistics...

    Any advice / opinions / debate would be much appreciated.

    Best regards!

    Fourie

  • #2
    Well it's exactly like IonTorrent advise - you can't deduplicate amplicon data on any platform, neither can you do it with HaloPlex.

    This means it is impossible to remove PCR artefacts from the data. This is why I much prefer hybridisation capture for these kind of studies, as deduplication is required to keep that source of false positives under control.

    My best advice is if you're stuck with this system, run samples in duplicate - at least then if you have false positives from the amplification, they shouldn't be present in the other replicate.

    Comment


    • #3
      I have the same problem. Is there a way to rum Mutect without removing duplicates (or presumed duplicates)?

      Comment


      • #4
        You can use the VariantTools Biocondctor package, this gives you the
        number of unique in read positions for variant and reference. Using a GenomicRanges object generated from a vcf you can make it report specifically for the positions of interest. This can be done using the tally function.

        Comment


        • #5
          Here is the code:

          library(gmapR)
          library(VariantTools)
          library(VariantAnnotation)
          library(BiocParallel)

          biocParam <- MulticoreParam(workers = ncores)

          fastaFile <- rtracklayer::FastaFile(referencefile)
          gmapGenome <- GmapGenome(fastaFile, create=TRUE,directory = referencefolder)

          print(tallied)

          vcf = readVcf(vcffile)
          called <- as(unlist(vcf),"VRanges")


          tally.param <- TallyVariantsParam(gmapGenome, high_base_quality = 0L,minimum_mapq = 10L,
          which = unique(as(called,"GRanges")),ignore_duplicates = FALSE,read_pos_breaks = c(1,10,120,330),
          variant_strand = 1)


          tallied = tallyVariants(bam_file,tally.param,BPPARAM = biocParam)
          matched = called %in% tallied
          matched_tallied = tallied %in% called

          cur_called = called[matched]
          tallied = unique(tallied[matched_tallied])

          sampleNames(tallied) = sampleNames(cur_called)

          elementMetadata(tallied) = c(elementMetadata(cur_called),elementMetadata(tallied))

          print(tallied)
          Last edited by IonTom; 07-07-2014, 09:37 AM.

          Comment


          • #6
            Originally posted by nbahlis View Post
            I have the same problem. Is there a way to rum Mutect without removing duplicates (or presumed duplicates)?
            Yes, just skip the dedup step. It's just an upstream step to skip.

            Comment


            • #7
              Hi Fourie,
              I was wondering what you ended up doing in terms of variant calling with the Ion CCP panel data? Did you just end up using the torrent variant caller? If so, how did you do the downstream analysis for filtering the variants etc? Did you use Ion Reporter?
              I find it quite incredible that the Ion data seems to be quite incompatible with ANY of the variant callers like Mutect, Varscan to do paired tumor normal analysis.
              Thanks for the favor of a reply.
              Regards

              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
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              9 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              49 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              67 views
              0 likes
              Last Post seqadmin  
              Working...
              X