Question: CNV Segmentation density - non-zero?
0
gravatar for dompearcejankowski
18 months ago by
United Kingdom
dompearcejankowski10 wrote:

Hi All,

tl;dr my segment means look odd https://imgur.com/a/5yOTodf

I'm analysing Illumina CytoSNP 850k data, that has been pre-processed as a (I believe) combination of Genome Studio and BlueFuse.

Confusingly to me, when the distribution of log2 segment means (i.e. the "copy number") are visualised, it's bimodal around 0.

My expectation would be that most segments display, on average, no copy change and a single peak should be evident at ~0 (this is ignoring any smaller peaks that are assumed to be associated with actual CNV).

Is this a problem with the segmentation algorithm that has been used? Would it be prudent to perform the segmentation again, myself?

Thanks for any help! Dom

The plot itself can be found here - https://imgur.com/a/5yOTodf

cnv snp copy number cytosnp genome • 448 views
ADD COMMENTlink modified 18 months ago by markus.riester500 • written 18 months ago by dompearcejankowski10
2
gravatar for markus.riester
18 months ago by
markus.riester500 wrote:

The normalization algorithms assume that test and control cells contain the same amount of DNA. In your case, the test ploidy is significantly below 2, thus you get a right shift in log2-ratios. Have a look at tools like ABSOLUTE or ASCAT that assign integer copy numbers to these peaks by jointly calculating purity and ploidy. This question also has been asked a few times here, e.g. C: "Shifted" copy number ratios coming from cnvkit?, so you might find more detailed information in the archive.

ADD COMMENTlink written 18 months ago by markus.riester500

Thanks a lot, really appreciated

ADD REPLYlink written 18 months ago by dompearcejankowski10
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: 1885 users visited in the last hour