Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Why not variant call before pileup?

    (I'm very new to this particular field.)

    From what I've read, it seems all variant call workflows do: alignment -> pileup -> variant call -> filtering/etc, regardless of whether the data is from whole genome or exome sequencing, and I understand this approach is valid for some applications/diseases. But for cancer applications where a tumor (or tumors) may be heterogeneous and have multiple mutation profiles (e.g. within an exon), is it not more valid to do variant calls on each read (cluster), then do a 'pileup' on the calls?

    For example:

    Code:
    ref: ...AACGTG...
    
         ...AACGTG... 800x *clusters* had this sequence
         ...AACGAG... 100x *clusters* had this sequence
         ...ATCGTG... 100x *clusters* had this sequence
    
    The above data (assume 100% confidence in base call) will be concluded as:
         ...AACGTG... 90% wildtype
         ...ATCGAG... 10% mutant with two mutations
    	 
    ... when, in fact, it is two separate mutations at 10% each.
    If anyone can point me to any papers/etc that discuss this, it is much appreciated.
    Last edited by obk; 01-23-2014, 03:40 PM.

  • #2
    No, because if you call variants in a single read, you can not distinguish read errors and real variations.

    Comment


    • #3
      I did see a talk at PAG XXII where the person called variants from a pileup and then went back to the individual reads to fit the variants into haplotypes enforced by the reads. Of course, I can't recall the talk, or if it was even a new thing! But that would give you the results you want.
      Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

      Comment


      • #4
        Dear TiborNagy
        I understand I wouldn't want to make a call based on a single read, but in the simple example, you'd call the two mutations 100 times each, which I think would give me some confidence that they are not erroneous reads... I think what I'm wondering is: if you have enough confidence in the base calling technology (or have enough coverage per unique read (like in the example)), what is the difference between:
        1) pileup reads to get consensus read -> variant call -> filter: is it real? -> real SNVs
        2) variant call individual reads -> 'pileup' variant calls -> filter: is it real? -> real SNVs
        (this question may be specific to amplicon sequencing...)
        Thanks.

        Comment


        • #5
          You can do that. When we do genotyping of populations, we get reads along the lines of what describe (mixed haplotypes). So one way we look at it is to align reads, track the variants of each read, then filter. The one difference is that our reads are all in synch (a stack of 100 reads at position 100,000, then a stack of 100 reads at position 200,000, etc). You would have some reads that end in between variants, leading to a little more work interpreting that.
          Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

          Comment


          • #6
            Thanks for your comments SNPsaurus.
            Do you have strategies to do any quantitative analysis based on the stack of reads? If it's amplicon sequencing that you're doing, then I imagine it is difficult to account for PCR duplicates.

            Comment


            • #7
              PCR duplicates are an issue, since we can't use different start and stop locations as a way to distinguish independent events. We were mostly looking for the presence of haplotypes in the populations so the precise level wasn't a concern. I was impressed by the "call from pileup then use the reads for phasing" approach I saw at the meeting because it did allow the use of common pipelines up until the last step and I think using common tools is increasingly important.
              Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Recent Developments in Metagenomics
                by seqadmin





                Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
                09-23-2024, 06:35 AM
              • seqadmin
                Understanding Genetic Influence on Infectious Disease
                by seqadmin




                During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

                Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
                09-09-2024, 10:59 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 10-02-2024, 04:51 AM
              0 responses
              13 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-01-2024, 07:10 AM
              0 responses
              22 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 09-30-2024, 08:33 AM
              0 responses
              26 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 09-26-2024, 12:57 PM
              0 responses
              18 views
              0 likes
              Last Post seqadmin  
              Working...
              X