Find Cpg Islands
2
2
Entering edit mode
7.7 years ago
niutster ▴ 100

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 • 15k views
ADD COMMENT
6
Entering edit mode
7.7 years ago

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.
ADD COMMENT
2
Entering edit mode
7.7 years ago
gammyknee ▴ 200

EMBOSS has a CpG island finder

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

ADD COMMENT
1
Entering edit mode

EMBOSS has a number of CpG related tools (see B.6.16. Applications in group Nucleic:cpg islands). The two CpG island finders: cpgplot and newcpgreport use slightly different criteria for calling a CpG island.

ADD REPLY

Login before adding your answer.

Traffic: 1194 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