Question: minus integer copy number when using --purity for cnvkit
0
gravatar for xiaoqiaoliu
2.7 years ago by
xiaoqiaoliu0 wrote:

Dear all,

i tried to use estimated purity to adjust the allele specific copy number output from cnvkit, but it would generate minus copy number integer value.

below is my command.

cnvkit.py call sample.cns -v sample.vcf --purity 0.39 -m clonal -o sample_purity.call.cns

below is the output .cns file with minus integer copy number value

chromosome      start   end     gene    log2    baf     cn      cn1     cn2     depth   probes  weight
2       236578628       237490302       AGAP1,GBX2,ASB18,IQCA1,ACKR3    -9.96578        0.00387202      -1      0       -1      26.4653 99      49.9862

did i use cnvkit in wrong way?

thanks.

cnvkit • 1.0k views
ADD COMMENTlink modified 2.7 years ago by Eric T.2.6k • written 2.7 years ago by xiaoqiaoliu0
0
gravatar for Kevin Blighe
2.7 years ago by
Kevin Blighe71k
Republic of Ireland
Kevin Blighe71k wrote:

It just indicates loss in the sample, i.e., deletion, as evidenced also by the log2 and baf (B allele frequency) values.

The CNVkit manual states that, when cn1 and cn2 differ, it is indicative of allelic imbalance:

Allelic imbalance, including copy-number-neutral loss of heterozygosity (LOH), is then apparent when a segment’s “cn1” and “cn2” fields have different values.

[source: http://cnvkit.readthedocs.io/en/stable/pipeline.html]

There are many copy number programs out there that say that they can detect LOH, allelic imbalance, etc., and none of these has ever made its way to be a clear winner. So, we have many groups using many different programs. Thankfully, one output format has more or less becomes standardised, which is generally referred to as CN segment data:

sample   chr   start   end  num_probes  segment_mean / log2_ratio

You have this data, which is good.

I would not use the cn1 and cn2 values for downstream analyses.

Kevin

ADD COMMENTlink modified 2.7 years ago • written 2.7 years ago by Kevin Blighe71k

however, how can i get absolute copy number and LOH status for these segment with minus cn?

ADD REPLYlink written 2.7 years ago by xiaoqiaoliu0

You can output the absolute copy number with the call command.

In this case, it looks like the region is entirely deleted in tumour, i.e., no copies. The BAF is virtually nil / zero, which, coupled with the log2 ratio, implies complete deletion.

log2 of -9.96578 is equivalent to a tumour-to-normal ratio of:

2^(-9.96578) = 0.00100000296
ADD REPLYlink written 2.7 years ago by Kevin Blighe71k
0
gravatar for Eric T.
2.7 years ago by
Eric T.2.6k
San Francisco, CA
Eric T.2.6k wrote:

You are right that the cn and cn2 values there should be 0, not -1. It's a quirk of rounding -- you can interpret -1 the same as 0, meaning full deletion.

I thought I'd squashed this behavior in recent versions of CNVkit. Are you using the latest version, 0.9.3?

ADD COMMENTlink written 2.7 years ago by Eric T.2.6k

thanks for the reminding, Eric. i use v0.8.5. i should update to 0.9.3.

ADD REPLYlink written 2.7 years ago by xiaoqiaoliu0
Please log in to add an answer.

Help
Access

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