SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
PubMed: Detection of low prevalence somatic mutations in solid tumors with ultra-deep Newsbot! Literature Watch 0 08-17-2012 12:00 AM
Somatic mutation profiling from 454 tumor-normal paired data pravee1216 Bioinformatics 0 05-12-2012 03:38 AM
Identifying Somatic mutation candidates question Phantom123 Bioinformatics 0 03-26-2012 06:47 AM
Mutation of FOXL2 in granulosa-cell tumors of the ovary. Chien-Yuan Chen Literature Watch 0 07-20-2009 03:01 PM

Reply
 
Thread Tools
Old 09-13-2012, 11:18 PM   #1
lethalfang
Member
 
Location: San Francisco, CA

Join Date: Aug 2011
Posts: 91
Default 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.
lethalfang is offline   Reply With Quote
Old 09-14-2012, 01:28 AM   #2
ymc
Senior Member
 
Location: Hong Kong

Join Date: Mar 2010
Posts: 498
Default

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?
ymc is offline   Reply With Quote
Old 09-14-2012, 02:16 AM   #3
adameur
Member
 
Location: Uppsala, Sweden

Join Date: Nov 2009
Posts: 23
Default

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..
adameur is offline   Reply With Quote
Old 09-14-2012, 10:07 AM   #4
lethalfang
Member
 
Location: San Francisco, CA

Join Date: Aug 2011
Posts: 91
Default

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.
lethalfang is offline   Reply With Quote
Old 09-15-2012, 05:36 AM   #5
adameur
Member
 
Location: Uppsala, Sweden

Join Date: Nov 2009
Posts: 23
Default

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..
adameur is offline   Reply With Quote
Old 09-17-2012, 07:39 PM   #6
lethalfang
Member
 
Location: San Francisco, CA

Join Date: Aug 2011
Posts: 91
Default

I'm trying this thing call Somatic Sniper now. I'll see how it goes:
http://gmt.genome.wustl.edu/somatic-...0.2/index.html
lethalfang is offline   Reply With Quote
Old 10-05-2013, 12:21 AM   #7
ymc
Senior Member
 
Location: Hong Kong

Join Date: Mar 2010
Posts: 498
Default

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.
ymc is offline   Reply With Quote
Old 10-05-2013, 10:36 AM   #8
lethalfang
Member
 
Location: San Francisco, CA

Join Date: Aug 2011
Posts: 91
Default

Quote:
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).
https://sites.google.com/site/strelk...variantcaller/
lethalfang is offline   Reply With Quote
Old 10-06-2013, 11:08 PM   #9
mikesh
Member
 
Location: California

Join Date: Jul 2012
Posts: 29
Default

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.
mikesh is offline   Reply With Quote
Reply

Tags
lifescope, ngs, snp, solid 5500, somatic mutations

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 08:21 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO