Question: Generation of QSEA TCGA 450K Human Methylation Calibration Dataset
gravatar for justinburgener
2.7 years ago by
justinburgener10 wrote:

Hi all,

This question is directed to those who may have experience with the QSEA R package, which proceeds from MEDIPS. 

I have MeDIP-seq data which I hope to begin analyzing very shortly, but unfortunately I'm not too familiar with processing genomic data. I've been able to following the tutorial for QSEA up until the normalization section, which suggests the use of a previously validated dataset. 450k TCGA data from lung is used in this example, and I such I presumed to use Head and Neck for mine. However I'm unable to figure out how to average the beta values within 500bp windows across all samples, as described within the tutorial, which looks like this:

I believe it may have been done by GRanges, or another R package which utilizes GRanges, but I am unable to generate the same object format. I've tried packages such as minfi, GRanges, but I lack the hands-on knowledge to further process the TCGA data after import.

I was wondering if anyone has any suggestions as to what R packages or software I can use to produce a similar R object.

Thank you,


ADD COMMENTlink modified 2.7 years ago by Matthias Lienhard140 • written 2.7 years ago by justinburgener10
Answer: Generation of QSEA TCGA 450K Human Methylation Calibration Dataset
gravatar for Matthias Lienhard
2.7 years ago by
Max Planck Institute for molecular Genetics, Berlin, Germany
Matthias Lienhard140 wrote:

Hi Justin,

I used two perl script to summarize the TCGA data into genome wide bins:

1) combines several sample files in a single table.

2) sumerizes bisulfite values in genome wide bins of fixed size by computing the average % methylation values.

I provide these scripts on my github repository:

The tables can then be averaged over normal and tumor samples in R directly.

For TCGA head and neck squamous cell carcinoma, I already prepared an 500 base overview table, in case you are interested, please contact me.

Best, Matthias

ADD COMMENTlink written 2.7 years ago by Matthias Lienhard140

Thank you for all the useful information! I'll take a look at the script and try to reproduce the tables to gain some more experience. I'll definitely contact you should I have difficulties and require the table you already prepared, very convenient! :) 

Thanks again,


ADD REPLYlink written 2.7 years ago by justinburgener10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 225 users visited in the last hour