Question: minus integer copy number when using --purity for cnvkit
0
gravatar for xiaoqiaoliu
5 months 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 • 248 views
ADD COMMENTlink modified 4 months ago by Eric T.2.3k • written 5 months ago by xiaoqiaoliu0
0
gravatar for Kevin Blighe
5 months ago by
Kevin Blighe32k
Republic of Ireland
Kevin Blighe32k 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 5 months ago • written 5 months ago by Kevin Blighe32k

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

ADD REPLYlink written 4 months 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 4 months ago by Kevin Blighe32k
0
gravatar for Eric T.
4 months ago by
Eric T.2.3k
San Francisco, CA
Eric T.2.3k 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 4 months ago by Eric T.2.3k

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

ADD REPLYlink written 4 months 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: 1052 users visited in the last hour