This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Find Cpg Islands

Hi,

I am looking for a tool to identify CpG islands, for a particular genome how do I find the position of these islands.

genome coordinates

2 answers

The UCSC browser contains a track for many genomes that provides the position of CpG islands.

For example, here is the table browser description for Human build hg19 http://genome.ucsc.edu/cgi-bin/hgTables?db=hg19&hgta_group=regulation&hgta_track=cpgIslandExt&hgta_table=cpgIslandExt&hgta_doSchema=describe+table+schema

CpG islands were predicted by searching the sequence one base at a time, scoring each dinucleotide (+17 for CG and -1 for others) and identifying maximally scoring segments. Each segment was then evaluated for the following criteria: GC content of 50% or greater, length greater than 200 bp, ratio greater than 0.6 of observed number of CG dinucleotides to the expected number on the basis of the number of Gs and Cs in the segment.

The CpG count is the number of CG dinucleotides in the island. The Percentage CpG is the ratio of CpG nucleotide bases (twice the CpG count) to the length. The ratio of observed to expected CpG is calculated according to the formula (cited in Gardiner-Garden et al. (1987)):

Obs/Exp CpG = Number of CpG * N / (Number of C * Number of G)   where N = length of sequence.

EMBOSS has a CpG island finder

http://emboss.sourceforge.net/apps/release/6.3/emboss/apps/cpgplot.html

Log in to answer this question.