Question: minus integer copy number when using --purity for cnvkit
0
gravatar for xiaoqiaoliu
25 days 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 • 103 views
ADD COMMENTlink modified 22 days ago by Eric T.2.2k • written 25 days ago by xiaoqiaoliu0
0
gravatar for Kevin Blighe
24 days ago by
Kevin Blighe23k
Republic of Ireland
Kevin Blighe23k 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 24 days ago • written 24 days ago by Kevin Blighe23k

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

ADD REPLYlink written 23 days 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 23 days ago by Kevin Blighe23k
0
gravatar for Eric T.
22 days ago by
Eric T.2.2k
San Francisco, CA
Eric T.2.2k 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 22 days ago by Eric T.2.2k

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

ADD REPLYlink written 21 days 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: 731 users visited in the last hour