Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Allele frequency in sample below 10%

    Hello,

    In our pipeline we have a quality control step where 23 snps are called with GATK and with realtime PCR. Lately as the coverage depth increases, we start to see that TaqMan is calling a homozygous snp and the high throughput sequencing method a heterozygous snp. And this mostly in case where one allele is present below 10%.

    I have difficulty to imagine cases where the allele frequency in one individual is below 10%.

    Can someone help me to make sense of this?

    Thanks

  • #2
    What platform are you using, is it exome-capture or WGS, and how much was the DNA amplified? Biases can greatly reduce the rate of one allele.

    Of course, there are other possibilities like repeats, chimerism, and contamination that can cause odd allelic ratios.

    Comment


    • #3
      Thanks for the reply.

      We are using a Illumina HiSeq, with exaome capture, Agilent SureSelect version 5. The lab technicians aim for 18 picomolar and this particular run produced around 100 10^6 reads.

      I am not sure chimerism can explain this only 1 or 3 snps out of 23 show this genotype and this within a batch of 20 samples. I did not check if the region where the snps were called are repeat rich. Will check.

      Comment


      • #4
        This can also be caused by errors in the reads. It's a good idea to remove duplicates if your data was PCR-amplified; this can reduce the rate of false-positive variants. Also, requiring a variant to be seen both on plus- and minus-mapped reads, and requiring the base to be called with some minimum average or maximum quality, are also filters that can reduce false positives caused by errors or biases. And it's possible for called variants to be mapping artifacts, too; you may want to try a different aligner and see if you get the same results.

        Comment


        • #5
          The reads went to the full GATK/Picard pipeline, from FixMate and markDuplicate to the indel realignment and base recalibration. They passed all tests, no strand bias for example or anything else. They look valid... But as you suggest I will try with BWA as the aligner, currently we are using novoalign.

          Comment


          • #6
            We might have found the cause .... All samples came from the same lane and the lane was almost over-clustered. The genotype in the HTS pipeline is the same in all samples pointing to a contamination. So to conclude, we suspect a contamination but we don't know where exactly.

            Thanks a lot Brian for your helpful comments.

            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
            11 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