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
                    Recent Advances in Sequencing Analysis Tools
                    by seqadmin


                    The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
                    05-06-2024, 07:48 AM
                  • seqadmin
                    Essential Discoveries and Tools in Epitranscriptomics
                    by seqadmin




                    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
                    04-22-2024, 07:01 AM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, 05-14-2024, 07:03 AM
                  0 responses
                  26 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 05-10-2024, 06:35 AM
                  0 responses
                  46 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 05-09-2024, 02:46 PM
                  0 responses
                  59 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 05-07-2024, 06:57 AM
                  0 responses
                  47 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X