This is a test version of Biostars. For the public version, visit https://www.biostars.org.
CNVkit: Overlapping Segments?

After running cnvkit.py segment, I've noticed that some segments appear to overlap with each other. For instance,

 18  18520842  60987290  BCL2  0.83552 424 131.823
 18  60987231  78016748  BCL2  0.04287 175 54.1546

The end of the first segment (60,987,290) overlaps with the beginning of the second segment (60,987,231). This is expected? I was expecting that the segments would be non-overlapping and the reason why is an issue is when overlapping mutations onto these segments. The mutation may be associated with multiple copy number predictions which should not be the case.

Thanks,

cnvkit
• 3,893 views • link
written by fongchunchan • fongchunchan

2 answers

In general the segments should not overlap. Diagnostics:

  • Which version of CNVkit are you using? The current is v0.7.11.
  • Is your sequencing protocol hybrid capture? Did you use the baited regions or "capture targets" BED file?
  • Did you use the "--split" argument with the batch or `target?
  • Which segmentation method did you use, the default CBS or another one?

(The note in the gainloss docs refers to a different situation where there is a breakpoint in the middle of the gene, so the "left" half of the gene belongs to one segment and the right half belongs to another.)

• 1 views • link
written by Eric T. Eric T.
  • 0.7.11
  • Yes it's a hybrid capture and yes I used a capture targets bed file.
  • I am not using the --split argument
  • Using the default segmentation method.
• 1 views • link
written by fongchunchan • fongchunchan

Thanks. My guess is that your target BED file contains overlapping intervals. Using the --split option will merge the overlapping intervals and then divide them evenly, which should solve your problem.

• 1 views • link
written by Eric T. Eric T.

Thanks. Using the --split option seems to have solved the overlapping intervals.

• 1 views • link
written by fongchunchan • fongchunchan

+1 for the detailed answer , I guess the OP should have a workaround now and if these answers solves the problem suing the --split parameter and changing the segmentation method, it would be nice to get the follow up reply from him. It keeps up with the interest of the forum.

• 1 views • link
written by ivivek_ngs ivivek_ngs

See this pdf:

https://media.readthedocs.org/pdf/cnvkit/v0.4.1/cnvkit.pdf and find there:

2.3.2 gainloss

"If segments are given, the log2 ratio value reported for each gene will be the value of the segment covering the gene. Where more than one segment overlaps the gene, i.e. if the gene contains a breakpoint, each segment’s value will be reported as a separate row for the same gene. If a large-scale CNA covers multiple genes, each of those genes will be listed individually".

• 1 views • link
written by natasha.sernova natasha.sernova

Log in to answer this question.

More posts like this

No similar posts found.