Generating 5mCG or 5hmCG values per gene from BS/OXBS data
0
0
Entering edit mode
hasche • 0
@37f4c9de
Last seen 2.6 years ago
United States

Hi, I use the methread function to read the coverage files from Bismark. BS_21=methRead("2368_bismark_bt2_pe.deduplicated.bismark.cov.gz", sample.id="2368_BS",assembly="mm10", header=F,skip=0,sep="\t",context="CpG", pipeline="bismarkCoverage", mincov = 5) However, the resulting methylkit object seems to lose the methylation percentage values generated by Bismark.

Bismark: https://user-images.githubusercontent.com/64626735/161638426-bda3dd3c-4c06-4541-8797-74c7cd73e96f.png

MethyKitObject: https://user-images.githubusercontent.com/64626735/161638488-5652f1a1-92ee-4716-87c6-f103e093c960.png

While I have found out a quick way of incorporating the meth %s for each methylation call from BS-seq experiments; the same is tricky when it comes to 5hmC data that I derive using the adjustmethylC function (because then I have two meth %s; one from BS and other from OXBS).

Would a plain subtraction of percent.meth per CpG (derived from Bismark) (BS minus OXBS) be the ideal way to go?

N.B. As expected, I have way more methylation calls from BS than the OXBS experiment.

I think that this whole situation stems from the fact that after running methread; the percentage methylation values of the coverage output from Bismark are lost.

My end goal is to generate 5mC and 5hmc signals per gene.

Please advise.

methyl methylKit • 854 views
ADD COMMENT

Login before adding your answer.

Traffic: 930 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6