Question: how to interpret a CNV scatter plot
0
gravatar for sara_beygi87
25 days ago by
sara_beygi870 wrote:

Dear all,

I needed help with the interpretation of these copy number to scatter plots from 2 patients. I'm interested in copy number changes in the CCR4 gene on chr3. The data is derived from a custom target panel. I realize from the plots that both patients have copy number losses in CCR4. However, I need help understanding 1) why patient C has 2 orange segments in the CCR4 region (would each indicate a different allele?) 2) whether we can comment on the bi-allelic or mono-allelic nature of the losses in each patient

Appreciate your thoughts!

CNV-scatter-plots-of-Chr-3

sequencing next-gen genome • 118 views
ADD COMMENTlink modified 24 days ago by Medhat8.8k • written 25 days ago by sara_beygi870

scatter plots from chr 3 in patients B and C

ADD REPLYlink written 25 days ago by sara_beygi870
3
gravatar for Amitm
24 days ago by
Amitm2.0k
UK
Amitm2.0k wrote:

Hi, I am assuming that this is from sequencing-based approach. PatientB has a very strong signal for a deletion segment spanning ~30thMb to 60th Mb region, which includes your gene. PatientC appears to have a deletion signal but which appears to be focal around CCR4 only. The orange line is the 'segmentation' outcome which takes into account how many probes/markers in surrounding area were different to, say the base trend i.e. copy ratio '0' (no change in copy number). Because the deletion signal in patientB is much large (~30thMb to 60th Mb), the 'segmentation' process binned that part of the region into a separate segment (i.e. a region which is in a different copy number state), and hence the separate orange line at copy ratio '-1'. In patientC, the deletion signal is not widespread and hence the 'segmentation' process hasn't created a separate 'segment' (orange line).

In terms of the fig. and whats available from the plot, it could be said that pat.B has mono-allelic loss (as the separate orange line is at '-1'). You could say the same for pat.C as well, albeit with somewhat lower confidence. Hope this helps.

ADD COMMENTlink written 24 days ago by Amitm2.0k

Hi Amit, thanks so much, this was very helpful. Sorry I didn't ask my question more clearly, but in my first question I was wondering why there are 2 separate orange segments in the CCR4 region in patient C (one around -0.6, the other one around -.1.8)? also thanks for your comment regarding mono vs bi-allelic, what would make you say a deletion is bi-allelic, if the segmentation line goes far down beyond -5? Thanks a lot again, much appreciated!

ADD REPLYlink written 24 days ago by sara_beygi870
1

No worries. To understand better the outcome in pat.C, I would recommend to look into the 'segmentation' output from which these plots are made. It should be in a format similar to

chr<TAB>Start<TAB>End<TAB>seg.mean<TAB>marker support

That should tell what those two small segments in pat.C are like. The genomic window of the plot is quite large and its hard say by looking at the plot only. Also, note that the 'segmentation' process is an algorithm and not everything that comes out of it can be necessarily explained biologically. To ascertain if there could be bi-allelic loss, given that this is sequencing data (?), you could open the BAM file in IGV browser. For a bi-allelic loss (meaning both alleles gone), there would be hardly any reads mapping in the loci. But it doesn't seem like so from the above plot. To add to all this, there is also aspect of the proportion of 'diseased' cell vs. 'non-diseased' cell in the sample. In an ideal situation with 100% of cells being 'diseased', bi-allelic should mean no reads/ signal at all.

ADD REPLYlink written 24 days ago by Amitm2.0k

Got it, thank you very much! I will check the segmentation output and yes you are correct, this is sequencing data.

ADD REPLYlink written 24 days ago by sara_beygi870
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: 1465 users visited in the last hour