how to combine methylation information from different platform?

Question: how to combine methylation information from different platform?

5.0 years ago by

liu4gre • 200

United States

liu4gre • 200 wrote:

Hi, guys,

I have a set of methylation data from 450K platform. And I want to test see how it compares to another set, which unfortunatly is from RRBS. I have the bed format file of the RRBS, which gives the CpG sites and methylation level. I am wondering how to combine 450K and RRBS, especially they don't have same CpG site.

Thanks.

rrbs methylation 450k • 3.1k views

ADD COMMENT • link •

modified 2.1 years ago by ibphuangchen • 10 • written 5.0 years ago by liu4gre • 200

5.0 years ago by

Charles Warden ♦ 6.6k

Duarte, CA

Charles Warden ♦ 6.6k wrote:

You can only compare signal that is commonly represented on both platforms.

If comparing 450k versus BS-Seq, you would obviously be ignoring a lot of information from the whole-genome BS-Seq. In your case, my understanding is that you may also have to work with differences in RRBS coverage between samples (even between technical replicates).

If it helps, there is a comparison between 450k signal and targeted BS-Seq data in Figure 4B of this paper (although that particular figure focuses on signal between islands rather than sites):

http://nar.oxfordjournals.org/content/41/11/e117.long

ADD COMMENT • link written 5.0 years ago by Charles Warden ♦ 6.6k

Thanks for the info.

It is true that 450K give less information than BS-Seq, such as RRBS. Actually I am only interested in those sites detected by 450K, as my sample is performed on this platform. The problem is my reference samples are performed on RRBS platform. I checked the CpG sites provided in these two samples, one is based on probe coordinates and another one is from the BED file. They don't match each other. So I am wondering what I should do?

Thanks.

ADD REPLY • link written 5.0 years ago by liu4gre • 200

You'll have to compare the files with a custom script, but the necessary information (chromsome, position, and beta / percentage methylation) are available in both cases. I believe the 450k annotation file (the .bpm file) provides both hg18 and hg19 coordinates (where the hg19 information is in CHR and MAPINFO). The 450k beta values are between 0 and 1 and the percentage methylation values in the .bed files (I assume from Bismark) will probably be between 0 and 100, so you need to change the scale by a factor of 100 (unless you are looking at differential methylation - then the chromosome and position are all you need to make something like the venn diagram in the figure that I was mentioning)

You can download that 450k annotation file from the Illumina website, and there is also a copy in the standalone version of COHCAP:

http://sourceforge.net/projects/cohcap/

ADD REPLY • link modified 5.0 years ago • written 5.0 years ago by Charles Warden ♦ 6.6k

Thank you very much for your help. Actually it is exact what I am asking. The coordinates from 450K (MAPINFO) are different from the coordinates provided in BBRS. How could I merge them? Or they are surely different CpG sites, meaning only the overlapped CpG sites can be extracted and compared (venn diagram)?

Thanks again.

ADD REPLY • link written 5.0 years ago by liu4gre • 200

Yes - you should only look for the overlap. The RRBS may just happen to not show enrichment in some areas covered by the 450k array, so this is my guess. You can visualize your alignment for a few cases to confirm this is in fact what is happening.

ADD REPLY • link written 5.0 years ago by Charles Warden ♦ 6.6k

Great. Thank you very much for the help.

ADD REPLY • link written 5.0 years ago by liu4gre • 200

Similar posts • Search »