Question: Why do I get a weird copy number in the tumour?
gravatar for w.jaratlerdsiri
4.1 years ago by
w.jaratlerdsiri10 wrote:

Weird copy number in Tumour

Please see the attached. I try to understand why Tumour 1 has most of the copy number equal to 1, but most of my past analyses (e.g. sequenza) and other works have most in the copy number =2 (same as Normal). Could you give me an insight of this interpretation?

I ran: batch export theta Sample_Tumor.cns Sample_Normal.cnr -o Sample.theta2.input 3./path/to/theta2/bin/RunTHetA Sample.theta2.input

Note I ran THeTA2 to get tumour purity and number of subclones. Is it something due to that "CNVkit only reports the estimated log2 copy ratio, and does not currently attempt a formal statistical test for estimating integer copy number"?

Thankyou, James

theta2 cnvkit • 1.7k views
ADD COMMENTlink modified 4.0 years ago by hywg20 • written 4.1 years ago by w.jaratlerdsiri10
gravatar for Eric T.
4.1 years ago by
Eric T.2.6k
San Francisco, CA
Eric T.2.6k wrote:

The procedure you ran is correct. THetA2 is not perfect, and it particularly struggles with samples that are noisy or have few large-scale CNVs. In this case I'd probably trust the result you got from Sequenza over THetA2.

I've updated CNVkit's "export theta" command to make it easy to run THetA2 with SNP allele frequencies to improve the results; you can try the latest CNVkit and THetA2 code from GitHub and see if that helps.

Otherwise, users have reported getting good results from PyClone, and it should be possible to run BubbleTree with CNVkit's segments if you coerce the data into the right form. Or you can proceed with Sequenza's purity estimate -- you only need a rough estimate of tumor purity to get useful results from CNVkit's call command.

ADD COMMENTlink modified 4.1 years ago • written 4.1 years ago by Eric T.2.6k

Hi Eric, may I ask how does the call command work? I have segment file output from sequenza. sequenza also gives a purity estimate. I want to scale the segment depth ratio by the tumor purity. The output of sequenaza segments is diff from cnvkit, so I am not sure how to use the call command. Thanks!

ADD REPLYlink written 2.4 years ago by Ming Tang2.6k
gravatar for hywg
4.0 years ago by
hywg20 wrote:

How do you run call on DNAnexus platform? I can run it on the Linux commandline, but found no way to run it on DNAnexus. Thanks in advance.

ADD COMMENTlink written 4.0 years ago by hywg20

The DNAnexus app doesn't expose all of CNVkit's functionality, only the basic pipeline. To run the call command you'll need to install CNVkit on your own Linux or OS X machine.

ADD REPLYlink written 4.0 years ago by Eric T.2.6k
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: 1889 users visited in the last hour