Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Sort of a noob here: looking for somatic mutation in tumors (NGS sequencing)

    As sort of a pilot study, we used Solid 5500 (Agilent SureSelect Exome) to sequence one
    tumor tissue and one matched normal tissue from the same patient.

    I used the vendor-provided LifeScope Genome Analysis Software to align the sequence and did the SNP callings. The initial SNP call has about 50,000 entries.

    For every SNP called in tumor, I aligned the Consensus Call of the same base position in the matched normal tissue, and then filtered out ones where the same base position in matched normal showed evidence of SNP.

    I also filtered out those in the dbSNP.

    We tested a small number of those potentially tumor-specific SNPs in Sanger sequencing, picking out ones with high coverage. 90% of them turned out to be sequencing error. For the other 10%, the same variant is also found in the matched normal.

    Well obviously, my noob approach didn't work.
    I picked out some SNPs to test because the bases had high coverage (although the reference vs. novel call is hardly 50-50, for most it's more like 80-20), although after the fact I do notice that many of them had pretty low Quality Score (QS = 20 or less).

    So anyway, not too many discussions about LifeScope, do you guys actually recommend some different software for such analysis? Or is there something I'm totally missing here?

    Thanks in advance.

  • #2
    I heard that the coverage for tumor should be at least 2x the normal cells, e.g. 100x vs 50x. Is this how u did it?

    why don't u give the standard bwa-GATK pipeline a chance also?

    Comment


    • #3
      Hi lethalfang

      I've been using Lifescope quite a lot for exome analyses and usually I detect around 50k-60k SNPs with 90% or so overlapping with dbSNP. So that seems to be about the same numbers as you find..

      I would be cautious to filter all dbSNP variants from the candidate list. As I understand dbSNP since can contain basically anything (including mutations). I usually use the dbSNPCommon instead for such filterings.

      Also, my experience is the variants with very high coverage can be problematic (could be due to repeats or similar) so maybe they are not the ideal candidates..

      Comment


      • #4
        Actually our average coverage was about 75X for both tumor and matched normal tissues.

        About 95% of the initial SNP calls overlap with dbSNP. I know there could be meaningful mutations in the dbSNP. I filtered it anyway because I thought there might be a greater proportion of meaningful mutations not in the dbSNP. But if I can't find anything in the other 5%, I'll go back to dig into those in the dbSNP.

        So adameur, what's your general experience with LifeScope?

        I also have a colleague who recommended using SomaticSniper.
        Has anyone used it? If so, how effective?

        Thanks.

        Comment


        • #5
          I'm not always impressed with Lifescope, but I think it works pretty well for variant calling.
          Around 50k SNPs with 95% overlapping with dbSNP sounds good to me..

          Comment


          • #6
            I'm trying this thing call Somatic Sniper now. I'll see how it goes:

            Comment


            • #7
              How do you filter SNPs reported by SomaticSniper? By default it reported too many SNPs for me (9572). The same dataset with Mutect gave me 42 PASS SNPs.

              Comment


              • #8
                Originally posted by ymc View Post
                How do you filter SNPs reported by SomaticSniper? By default it reported too many SNPs for me (9572). The same dataset with Mutect gave me 42 PASS SNPs.
                You can filter by quality score.
                Anyway, I've also used Strelka, and I found that to be easier, faster, and better than MuTect (for our data anyway).
                Access Google Sites with a personal Google account or Google Workspace account (for business use).

                Comment


                • #9
                  Generally speaking, if there is a sequencing error you cannot get rid of it if you'll sequence deeper. And 75x seems enough.
                  Some control, i.e. known mutation in tumor is desirable to fine-tune your pipeline. If you don't have one you can try intersect your SNPs with COSMIC database and check all the matches.

                  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