Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
finding promoters using ensembl biomart narges Bioinformatics 3 12-03-2012 05:16 AM
Classify Bacterial and Archaeal fragments from 454-Metagenome Mithril Bioinformatics 4 10-11-2012 09:23 PM
classify RNA library Palgrave Bioinformatics 2 08-17-2011 06:13 AM
How classify sequences? Asgar Bioinformatics 0 12-07-2010 05:52 AM
Newbler --- classify many reads as Repeat NSTbioinformatics 454 Pyrosequencing 2 06-10-2010 04:43 AM

Thread Tools
Old 01-04-2013, 04:30 AM   #1
Location: Bilbao

Join Date: Apr 2012
Posts: 44
Default How to classify promoters according to CpG density

Hello everybody,

A research collaborator has asked us to classify promoters according to their CpG density. Data comes from a RRBS (Reduced Representation Bisulfite Sequencing) analysis in mouse (using mm9 sequence as template) and I don't have a clue on how to do this.

-Any of you knows about any tool, script or R-package that allows to do this ?
-Is there any repository where I can download a BED file where this kind of information is already available for organisms different than human?

We've been reading different papers that present such data but I haven't been able to understand how they manage to obtain such classification and then represent it graphically.

Any insight would be appreciated.

Thanks in advance.


PD. I've been given this set of conditions. I write them here just in case anybody (even the authors of some similar analysis) may recognize an analog issue reading them, so that they can suggest a tool, script or R package that perform this calculation:

"To determine promoter classes, we measured the GC content and the CpG ratio of observed to expected values in sliding 500-bp windows with a 5-bp offset in regions -900 bp to +400 bp relative to the TSS. Promoter classes were defined as follows: LCPs contain no 500-bp window with a CpG ratio >0.45; HCPs contain at least one 500-bp window with a CpG ratio >0.65 and GC content >55%; ICPs do not meet the previous criteria."

Last edited by Jluis; 01-04-2013 at 05:34 AM.
Jluis is offline   Reply With Quote
Old 01-06-2013, 04:33 PM   #2
Senior Member
Location: Australia

Join Date: Sep 2008
Posts: 136

The Bioconductor package Repitools has a function cpgDensityCalc which will run on any BSgenome organism available on Bioc, and whichever gene annotations you supply.
frozenlyse is offline   Reply With Quote
Old 01-07-2013, 12:28 AM   #3
Location: Bilbao

Join Date: Apr 2012
Posts: 44

Thank you frozenlyse, I'll take a look at it and see if it works in my particular analysis.

Jluis is offline   Reply With Quote

cpg island, promoters, rrbs

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 01:57 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO