Question: SciClone: Tumor purity estimation
gravatar for Xiaolong Zhang
2.8 years ago by
Xiaolong Zhang0 wrote:

I am trying to find a mutation-based tumor purity estimation method to due with my WES data. I was noticed that SciClone could cluster mutations of muti-region sampled tumors (remove SNVs located in CNV regions) as well as estimate tumor purities.

However, I found that the purity estimating function is closed ("doPurityScaling = FALSE") by anthor in the script of "sciClone" function. Anyhow, I edited this script to turn on this function and estimated the tumor purities of different samples. Surprisingly, the purity scores are highly unstable (range from 0.0193 to 1.6606) and many purity values were more than 1.

How dose it happens? If the purity estimating function is unusable in SciClone? If not, do you have any good method to estimate tumor purity with mutations VAF (not based on CNVs)? Thank you!

purity R genome • 1.5k views
ADD COMMENTlink modified 2.7 years ago by Chris Miller21k • written 2.8 years ago by Xiaolong Zhang0
gravatar for Chris Miller
2.7 years ago by
Chris Miller21k
Washington University in St. Louis, MO
Chris Miller21k wrote:

As you've found, there was originally some experimental (and fairly naive) purity estimation written into sciClone, but it did not work as well as I'd like. The only thing it's doing is taking the mean VAF of the founding clone (generally cluster 1) and doubling it.

In cases where your clusters aren't well constructed (due to a small number of points, or due to CN-altered points that were not removed), that purity estimation can be pretty far off.

ADD COMMENTlink written 2.7 years ago by Chris Miller21k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1152 users visited in the last hour