CNVkit cnv_ztest with low log2 values
Entering edit mode
2.4 years ago

Dear all, I am using CNVkit to analyse exome sequencing data (germline samples) in order to find new copy number variants. For the pipeline analysis I followed all the found suggestions for germline samples and here are the steps:

1) create a reference using all samples batch --normal ALL_BAM/*.bam --targets Exome_SureSelect_QXTV7_forCNVkit.bed --fasta hs37d5.fa --access access-5k-mappable.hg19.bed --output-reference reference.cnn

2) batch command using the previously created reference batch ALL_BAM/*.bam -r reference.cnn -d results/

3) add ci column

cd results/
for i in *.cnr ; do segmetrics -s `basename ${i%%.cnr}`.cn{s,r} --ci ; done

4) call command

for i in *segmetrics* ; do call $i --filter ci -m clonal --center mode -o call_cnvkit/`basename ${i%%segmetrics*}`.call.cns; done

5) cnv_ztest

for i in *.cnr; do $i -t -s call_cnvkit/`basename ${i%%.cnr}` -o cnv_ztest/`basename ${i%%.cnr}`.ztest.cnr; done

I am now looking at cnv_ztest results but I am not understanding log2 values reported in the file because they have quite negative values in all samples (like between -21 and -8). Is there something I am missing? Or probably is there an error in my pipeline? I also had a look at this thread but because I am using germline samples I didn't use the --drop-low-coverage as suggested. Do you have any idea about what is going wrong? Thanks in advance!


cnvkit germline exome sequencing cnv_ztest • 686 views
Entering edit mode

Hi Stephania,

is it for sure that your targeted enrichment kit was the same for both cases and controls?

Entering edit mode

Hi German, yes the target enrichment kit is the same for all samples.

Entering edit mode

Very low negative numbers show homozygous deletions. They are found as very low coverage. Somehow your tumor samples do not have coverage in these regions.

Having many recurrent homozygous deletions is very unlikely. Some genes such as CDKN2A are "famous" for that, but still, I would not expect more than 20% of samples to have it.


Login before adding your answer.

Traffic: 2164 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6