Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Why 30X

    Hello
    I'm in the process of organising sequencing of a number of tumours and matched blood samples. I have to related questions regarding coverage. I often hear of 30X coverage as if it is some kind of magic number but I'm curious to know where it has come from. Is it realistic to ask for 20X for my blood samples, since I am only interested in them as a control for my tumours, or is it some kind of logarithmic scale whereby 20X will only get me 10% of the information but spending lots of money on 40X will only add another 2%. My second question is about the tumours, which are likely to be a mixed population of cells. It seems to me that the heterogenous nature of the samples will mean 30X will not be enough, but what will be enough, is 50X OK, 60X if it's very mixed? I realise that these questions are a bit open ended and the answers will depend on what I want to do with the data, but I would appreciate any feedback or pointers to papers where this has been calculated. We're doing it on a HiSeq if error rates come into the calculations.
    Cheers
    Henry

  • #2
    I too would love to see answers to these same questions and several others. Is there a point of diminishing return on depth of coverage for cancer samples in detecting mutations? Can there also be a point at which your coverage is so high that the mutations one is detecting could actually be artifacts? Does the type of cancer (and whether the sample is from a cell line or actual tumor sample) have any bearing on the coverage one would use when sequencing for mutation detection?

    Comment


    • #3
      I've wondered about these too. One approach to answer it would be to take a genome at much higher coverage & then randomly sample the reads from the BAM file at different coverages & re-call variants. Obviously a bit compute intensive, but perhaps well worth it.

      I would guess sample purity is the main determinant for tumors. Once you have some tumor-specific mutations you can start estimating purity. Of course, the tumor itself is heterogeneous. This complicates things further. So one probably needs to state a goal clearly ("have 95% chance of detecting a mutation present in 80% of the tumor in a sample which is 55% tumor") and then compute out target coverage from there.

      Some tumors will be more heterogeneous than others; a first guess is that those with higher environmental mutation load OR those which have more DNA repair defects will be more heterogeneous.

      Comment


      • #4
        Originally posted by krobison View Post
        One approach to answer it would be to take a genome at much higher coverage & then randomly sample the reads from the BAM file at different coverages & re-call variants. Obviously a bit compute intensive, but perhaps well worth it.
        That is the approach someone else had suggested to me. So an N=2 makes me feel better about doing just that.

        Comment


        • #5
          That seems a pretty sensible idea to me too. Whether I can get it to work by Monday when I send my samples off is another thing. I think my suppliers are quite flexible in letting me do 20-30X now and then spend more money later to bump it up if need be.
          PS, what time do you get up Geneus. I'm in Europe so I'm supposed to be up and about at that time, but I make your first reply about 6.15am in your part of the world.

          Comment


          • #6
            Originally posted by henry.wood View Post
            PS, what time do you get up Geneus.
            Henry...let's just say early...and leave it at that.

            Comment


            • #7
              I've wondered about these too. One approach to answer it would be to take a genome at much higher coverage & then randomly sample the reads from the BAM file at different coverages & re-call variants. Obviously a bit compute intensive, but perhaps well worth it.
              What a coincidence (or maybe not) - I just started to consider doing this today (before reading this thread). I'll report back if I ever get around to it.

              Comment


              • #8
                for 'why 30x', check out the following url and corresponding reference. There's certainly more to say on the topic, but I'd thought I'd shoot this quick reply first:


                Bentley et al (2008) Accurate whole human genome sequencing using reversible terminator chemistry. Nature 456, 53-59
                doi:10.1038/nature07517

                Comment


                • #9
                  Originally posted by krobison View Post
                  I've wondered about these too. One approach to answer it would be to take a genome at much higher coverage & then randomly sample the reads from the BAM file at different coverages & re-call variants. Obviously a bit compute intensive, but perhaps well worth it.
                  One of the first WGS papers out of BGI did exactly this in the supplemental materials.

                  Comment


                  • #10
                    Originally posted by nilshomer View Post
                    One of the first WGS papers out of BGI did exactly this in the supplemental materials.
                    Do you have the details of which journal this appeared in?

                    Thanks.

                    Comment


                    • #11
                      Originally posted by Geneus View Post
                      Do you have the details of which journal this appeared in?

                      Thanks.
                      Hopefully you tried searching yourself, but see the section "Depth effect..." in http://www.nature.com/nature/journal/v456/n7218/full/nature07484.html".

                      Comment


                      • #12
                        Originally posted by nilshomer View Post
                        Hopefully you tried searching yourself
                        Thank you very much.

                        By the way, hope is not a strategy.

                        Comment


                        • #13
                          Thanks all for the input. It seems that I can save a bit of money on sequencing my control samples, since an extra 10X coverage going from 20-30X only catches an extra 1-2% of the SNPs. This will let me spend more on my tumours, which is what I'm really interested in. If it ever sees the light of day I'll try and remember to make a mention of the forum in the acknowledgments, even if it's just to get on the list of the "Greatest papers in the world".

                          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
                          8 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, Yesterday, 06:07 PM
                          0 responses
                          8 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