# Generic eukaryotic core promoter prediction using structural features of DNA

Thomas Abeel, Yvan Saeys, Eric Bonnet, Pierre Rouzé, Yves van de Peer

► To cite this version:

Thomas Abeel, Yvan Saeys, Eric Bonnet, Pierre Rouzé, Yves van de Peer. Generic eukaryotic core promoter prediction using structural features of DNA. Genome Research, 2008, 18 (2), pp.310-323. 10.1101/gr.6991408 . cea-04497272

HAL Id: cea-04497272

https://cea.hal.science/cea-04497272

Submitted on 9 Mar 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Generic eukaryotic core promoter prediction using structural features of DNA

Thomas Abeel,¹,² Yvan Saeys,¹,² Eric Bonnet,¹,² Pierre Rouzé,¹,²,³ and Yves Van de Peer¹,²,<fnref n="4" />

¹Department of Plant Systems Biology, Flanders Institute for Biotechnology (VIB), 9052 Gent, Belgium; ²Department of Molecular Genetics, Ghent University, 9052 Gent, Belgium; ³Laboratoire Associé de l'INRA (France), Ghent University, 9052 Gent, Belgium

Despite many recent efforts, in silico identification of promoter regions is still in its infancy. However, the accurate identification and delineation of promoter regions is important for several reasons, such as improving genome annotation and devising experiments to study and understand transcriptional regulation. Current methods to identify the core region of promoters require large amounts of high-quality training data and often behave like black box models that output predictions that are difficult to interpret. Here, we present a novel approach for predicting promoters in whole-genome sequences by using large-scale structural properties of DNA. Our technique requires no training, is applicable to many eukaryotic genomes, and performs extremely well in comparison with the best available promoter prediction programs. Moreover, it is fast, simple in design, and has no size constraints, and the results are easily interpretable. We compared our approach with 14 current state-of-the-art implementations using human gene and transcription start site data and analyzed the ENCODE region in more detail. We also validated our method on 12 additional eukaryotic genomes, including vertebrates, invertebrates, plants, fungi, and protists.

[Supplemental material is available online at www.genome.org.]

Eukaryotic genomes are being sequenced at an ever-increasing pace. At the moment, nearly 50 complete genomes of eukaryotes are publicly available, and many more are in the pipeline to be sequenced in the next few years (Liolios et al. 2006). The proliferation of genome sequencing projects has driven the search for fast ways of sequence-based structural annotation, which involves the identification of genes and the modeling of their correct gene structure (Claverie et al. 1997; Mathé et al. 2002; Zhang 2002; Wang et al. 2004). Although great progress has been achieved in gene prediction, for instance by using comparative approaches (Wasserman et al. 2000; Liu et al. 2004; Jin et al. 2006; Wang and Zhang 2006), one of the more difficult tasks in the annotation of whole genomes remains the accurate identification and delineation of promoters (Fickett and Hatzigeorgiou 1997; Ohler 2000, 2001; Bajic et al. 2004, 2006a). Nevertheless, the prediction of the regions that control the transcriptional activation of genes is important for various reasons (Smale 2001; Butler and Kadonaga 2002; Bajic et al. 2004; Sonnenburg et al. 2006). On the one hand, promoter prediction can be used for the discovery of genes that are missed by gene predictors and/or for which experimental support (ESTs, cDNAs, etc.) is not available. On the other hand, the prediction of promoters is important for guiding further in silico searches and experimental work, for instance in narrowing down the regions that play the most important role in transcriptional regulation (Bajic et al. 2004, 2006a; Carninci et al. 2006; Solovyev et al. 2006).

The promoter is commonly referred to as the region upstream of a gene that contains the information permitting the proper activation or repression of the gene that it controls (Pe-

dersen et al. 1999; Smale and Kadonaga 2003). The promoter region itself is typically divided into three parts: (1) the core promoter, which is the region that is responsible for the actual binding of the transcription apparatus and which is typically situated ~35 bp upstream of the transcription start site (TSS); (2) the proximal promoter, a region containing several regulatory elements, which ranges up to a few hundred base pairs upstream of the TSS; and (3) the distal promoter, which can range several thousands of base pairs upstream of the TSS and contains additional regulatory elements called enhancers and silencers.

It has been known for quite some time that the properties of promoter regions are considerably different from those of other parts in the genome (Pedersen et al. 1998; Aerts et al. 2004; Florquin et al. 2005; Fukue et al. 2005; Tabach et al. 2007). Some features that have proven useful in the detection of promoters in vertebrate genomes are the so-called CpG islands close to the TSS (Delgado et al. 1998; Ioshikhes and Zhang 2000; Hannenhalli and Levy 2001), the presence of typical transcription factor binding sites (Solovyev and Shahmuradov 2003; Choi et al. 2004; Ohler 2006), and statistical properties of the core and proximal promoter (Down and Hubbard 2002; Bajic et al. 2006b; Fitzgerald et al. 2006). The similarities between orthologous promoters (Solovyev and Shahmuradov 2003; Jin et al. 2006) and information from mRNA transcripts (Liu and States 2002) have also been used to identify promoters. The more recent and sophisticated Promoter Prediction Programs (PPPs) look for these promoter-specific characteristics by using machine learning techniques such as discriminant analyses, Hidden Markov Models, and Artificial Neural Networks to predict and delineate promoters (for reviews, see Fickett and Hatzigeorgiou 1997; Rombauts et al. 2003; Bajic et al. 2004, 2006a; Sonnenburg et al. 2006). Programs and tools based on these techniques are difficult to train because they require a large amount of high-quality training data, preferably from an experimental setting (Munch and Krogh 2006).

<footnote n="4">⁴Corresponding author.
E-mail yves.vandepeer@psb.ugent.be; fax 32-(0)-9-33-13-809.
Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.6991408.</footnote>
# Generic eukaryotic core promoter prediction using structural features of DNA

Thomas Abeel,<sup>1,2</sup> Yvan Saeys,<sup>1,2</sup> Eric Bonnet,<sup>1,2</sup> Pierre Rouzé,<sup>1,2,3</sup>
and Yves Van de Peer<sup>1,2,4</sup>

<sup>1</sup>Department of Plant Systems Biology, Flanders Institute for Biotechnology (VIB), 9052 Gent, Belgium; <sup>2</sup>Department of Molecular Genetics, Ghent University, 9052 Gent, Belgium; <sup>3</sup>Laboratoire Associé de l'INRA (France), Ghent University, 9052 Gent, Belgium

Despite many recent efforts, in silico identification of promoter regions is still in its infancy. However, the accurate identification and delineation of promoter regions is important for several reasons, such as improving genome annotation and devising experiments to study and understand transcriptional regulation. Current methods to identify the core region of promoters require large amounts of high-quality training data and often behave like black box models that output predictions that are difficult to interpret. Here, we present a novel approach for predicting promoters in whole-genome sequences by using large-scale structural properties of DNA. Our technique requires no training, is applicable to many eukaryotic genomes, and performs extremely well in comparison with the best available promoter prediction programs. Moreover, it is fast, simple in design, and has no size constraints, and the results are easily interpretable. We compared our approach with 14 current state-of-the-art implementations using human gene and transcription start site data and analyzed the ENCODE region in more detail. We also validated our method on 12 additional eukaryotic genomes, including vertebrates, invertebrates, plants, fungi, and protists.

[Supplemental material is available online at www.genome.org.]

Eukaryotic genomes are being sequenced at an ever-increasing pace. At the moment, nearly 50 complete genomes of eukaryotes are publicly available, and many more are in the pipeline to be sequenced in the next few years (Liolios et al. 2006). The proliferation of genome sequencing projects has driven the search for fast ways of sequence-based structural annotation, which involves the identification of genes and the modeling of their correct gene structure (Claverie et al. 1997; Mathé et al. 2002; Zhang 2002; Wang et al. 2004). Although great progress has been achieved in gene prediction, for instance by using comparative approaches (Wasserman et al. 2000; Liu et al. 2004; Jin et al. 2006; Wang and Zhang 2006), one of the more difficult tasks in the annotation of whole genomes remains the accurate identification and delineation of promoters (Fickett and Hatzigeorgiou 1997; Ohler 2000, 2001; Bajic et al. 2004, 2006a). Nevertheless, the prediction of the regions that control the transcriptional activation of genes is important for various reasons (Smale 2001; Butler and Kadonaga 2002; Bajic et al. 2004; Sonnenburg et al. 2006). On the one hand, promoter prediction can be used for the discovery of genes that are missed by gene predictors and/or for which experimental support (ESTs, cDNAs, etc.) is not available. On the other hand, the prediction of promoters is important for guiding further in silico searches and experimental work, for instance in narrowing down the regions that play the most important role in transcriptional regulation (Bajic et al. 2004, 2006a; Carninci et al. 2006; Solovyev et al. 2006).

The promoter is commonly referred to as the region upstream of a gene that contains the information permitting the proper activation or repression of the gene that it controls (Pe-

dersen et al. 1999; Smale and Kadonaga 2003). The promoter region itself is typically divided into three parts: (1) the core promoter, which is the region that is responsible for the actual binding of the transcription apparatus and which is typically situated ~35 bp upstream of the transcription start site (TSS); (2) the proximal promoter, a region containing several regulatory elements, which ranges up to a few hundred base pairs upstream of the TSS; and (3) the distal promoter, which can range several thousands of base pairs upstream of the TSS and contains additional regulatory elements called enhancers and silencers.

It has been known for quite some time that the properties of promoter regions are considerably different from those of other parts in the genome (Pedersen et al. 1998; Aerts et al. 2004; Florquin et al. 2005; Fukue et al. 2005; Tabach et al. 2007). Some features that have proven useful in the detection of promoters in vertebrate genomes are the so-called CpG islands close to the TSS (Delgado et al. 1998; Ioshikhes and Zhang 2000; Hannenhalli and Levy 2001), the presence of typical transcription factor binding sites (Solovyev and Shahmuradov 2003; Choi et al. 2004; Ohler 2006), and statistical properties of the core and proximal promoter (Down and Hubbard 2002; Bajic et al. 2006b; Fitzgerald et al. 2006). The similarities between orthologous promoters (Solovyev and Shahmuradov 2003; Jin et al. 2006) and information from mRNA transcripts (Liu and States 2002) have also been used to identify promoters. The more recent and sophisticated Promoter Prediction Programs (PPPs) look for these promoter-specific characteristics by using machine learning techniques such as discriminant analyses, Hidden Markov Models, and Artificial Neural Networks to predict and delineate promoters (for reviews, see Fickett and Hatzigeorgiou 1997; Rombauts et al. 2003; Bajic et al. 2004, 2006a; Sonnenburg et al. 2006). Programs and tools based on these techniques are difficult to train because they require a large amount of high-quality training data, preferably from an experimental setting (Munch and Krogh 2006).

<footnote n="4">⁴Corresponding author.
E-mail yves.vandepeer@psb.ugent.be; fax 32-(0)-9-33-13-809.
Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.6991408.</footnote>
However, for most of the new genome projects there is only a
limited amount of such data available. Another problem is that
the outcome of programs based on these techniques is often dif-
ficult to interpret (Ratsch et al. 2006). Furthermore, all hitherto
available programs are species-specific; i.e., they are trained on
one species and are able to predict promoters only for that par-
ticular species. Another drawback of most PPPs is that they de-
pend on specific motifs expected to be present in the core pro-
moter. Indeed, some programs (Promoter2.0 [Knudsen 1999];
Eponine [Down and Hubbard 2002]; NNPP2.2 [Burden et al.
2004]) are based on the explicit presence of motifs such as the
TATA box, which are very common in certain species, such as
yeast (Struhl 1989), but much less common in mammals or
plants (Suzuki et al. 2001; Butler and Kadonaga 2002; Fukue et al.
2004; Florquin et al. 2005). This hampers the ability of a single
program to analyze different species with the same model, and
does not facilitate the discovery of different types of promoters.
Finally, there is the issue of scalability, as the best performing
programs are unable to process large datasets (Ohler et al. 2002;
Bajic et al. 2004, 2006a; Sonnenburg et al. 2006).

In light of all these caveats, we propose a simple technique
for identifying and delineating (core) promoters that is based on
the properties of long stretches of DNA. It has indeed been shown
that sequence properties such as GC content and more general
chemo-physical properties of the DNA, such as stabilizing energy
of Z-DNA (Ho et al. 1990), DNA denaturation values (Blake and
Delcourt 1998; Blake et al. 1999), protein-induced deformability
(Olson et al. 1998), and duplex-free energy (Sugimoto et al.
1996), among others (for review, see Florquin et al. 2005), can be
used to describe (core) promoters, and to discriminate between
(core) promoter sequences and non (core) promoter sequences
(Ohler et al. 2001; Florquin et al. 2005; Kanhere and Bansal 2005;
Uren et al. 2006; Wang and Benham 2006). Because all these
properties are calculated from conversion tables using di- or tri-
nucleotides, one may argue that these properties are in fact ex-
actly the same as the nucleotide sequence and do not offer any
additional information. However, several studies have shown
that this is not the case. Both Liao et al. (2000) and Baldi et al.
(1998) have analyzed the correlation between the different prop-
erties, and their main conclusion was that the properties are
largely independent. Moreover, Florquin et al. (2005) have clus-
tered promoters based on these structural properties. The genes
associated with the promoters in each cluster varied greatly for
the different properties, which again indicates that the different
properties contain complementary information. Bode et al.
(2006) have shown that it is very hard to identify scaffold/matrix
attachment regions (SMARs) from the sequence, as the important
scaffold proteins recognize structural features instead of specific
nucleotide sequences. Structural properties are known to have
long-range interactions (up to 10 kb), so they can exhibit prop-
erties that are not visible in the sequence (Merling et al. 2003;
Faiger et al. 2006). The Human Genomic Melting Map (Liu et al.
2007) shows a correlation between GC content and DNA dena-
turation, but due to the cooperative nature of DNA denaturation,
this correlation is weaker on scales <500 bp.

We applied these different properties to the problem of pro-
moter prediction and show that some properties give better per-
formance than others, which again indicates differences in in-
formation content. In particular, we present the Easy Promoter
Prediction Program (EP3), which uses GC content and large-scale
structural features of DNA to identify and delineate promoter
regions in whole-genome sequences. EP3 was applied to the

human genome and compared with other state-of-the-art PPPs
using two different evaluation techniques: One is commonly
used for PPPs, while the other is based on a novel scheme also
considering intergenic predictions that do not fall in a promoter.
We evaluated the different ways of encoding the genomic DNA,
and we also applied EP3 to a set of human noncoding RNA genes.
Finally, we evaluated EP3 on 12 other eukaryotic genomes.

Results and Discussion

Properties of the core promoter region

Figure 1 shows examples of numerical profiles for human pro-
moter sequences. These profiles were obtained by lining up all
the promoter sequences with the TSS at the same position. Next,
we converted all sequences using the different structural proper-
ties of DNA into numeric sequences and, for each position, plot-
ted the average over all numeric sequences. The X-axis is the
position relative to the TSS, and the Y-axis is the normalized
value of that property. It is clear that the human core promoter
adopts a very specific intrinsic structure that stretches over quite
long distances. Either a large peak is visible, representing highly
stable regions in the DNA, or a cleft is visible, which represents
highly unstable regions in the DNA. For instance, the GC profile
shows a peak, while the AT profile shows a cleft. When zoomed
in, it can be noticed that, within the broader region of stability,
there is a small unstable region that corresponds to the TSS, and
a second one that corresponds to the region where TBP binds (see
below). The promoter prediction technique we present will
identify those higher regions and are, thus, most likely to contain
a TSS.

To see how general these properties of promoter sequences
are, we created large-scale (2000 bp up- and downstream of the
TSS) datasets for 16 different species. For each of these datasets,
we analyzed the structural profiles. Over long-range distances, we
observe several types of profiles (Fig. 2). The first one applies
mostly to protists. In this type, the profile slowly decreases a few
hundred base pairs upstream of the TSS, shows a transition of a
strong peak and cleft on the TSS, and is flat again downstream of
the TSS, which may indicate a relatively small promoter area in
these organisms. The second type of profile is observed in mam-
mals and shows a broad elevated region around the TSS, with
again a strong local peak and drop on the TSS itself. Between
these two extremes, we observe a more gradual transition in
fungi and plants. The extremely high and low values that can be
seen on the TSS are discussed in more detail in the next section.

There seems to be a relationship between the genome size
and the size of the peaks or clefts: In prokaryotes (not shown) and
protists, these cover ~100 bp, while they gradually increase to
several hundreds of base pairs in fungi and up to ~1000 bp in
plants. In animals, smaller genomes such as that of *Drosophila*
have a peak ranging from a few hundred base pairs on each side
of the TSS; this gradually increases up to ~2000 bp in mammals.
In conclusion, the peak is most obvious in mammals, but is also
clearly visible in fish, insects, fungi, and plants; it seems consid-
erably smaller in protists.

In humans, as well as in other eukaryotes, three types of
RNA polymerases (RNAP) exist that are used to transcribe differ-
ent types of genes. All three of them manage to identify the TSS
in a whole-genome setting in their own way. Here, we focus on
humans because this species is among the best documented ones.
We retrieved the promoter sequences for several types of genes
However, for most of the new genome projects there is only a
limited amount of such data available. Another problem is that
the outcome of programs based on these techniques is often dif-
ficult to interpret (Ratsch et al. 2006). Furthermore, all hitherto
available programs are species-specific; i.e., they are trained on
one species and are able to predict promoters only for that par-
ticular species. Another drawback of most PPPs is that they de-
pend on specific motifs expected to be present in the core pro-
moter. Indeed, some programs (Promoter2.0 [Knudsen 1999];
Eponine [Down and Hubbard 2002]; NNPP2.2 [Burden et al.
2004]) are based on the explicit presence of motifs such as the
TATA box, which are very common in certain species, such as
yeast (Struhl 1989), but much less common in mammals or
plants (Suzuki et al. 2001; Butler and Kadonaga 2002; Fukue et al.
2004; Florquin et al. 2005). This hampers the ability of a single
program to analyze different species with the same model, and
does not facilitate the discovery of different types of promoters.
Finally, there is the issue of scalability, as the best performing
programs are unable to process large datasets (Ohler et al. 2002;
Bajic et al. 2004, 2006a; Sonnenburg et al. 2006).

In light of all these caveats, we propose a simple technique
for identifying and delineating (core) promoters that is based on
the properties of long stretches of DNA. It has indeed been shown
that sequence properties such as GC content and more general
chemo-physical properties of the DNA, such as stabilizing energy
of Z-DNA (Ho et al. 1990), DNA denaturation values (Blake and
Delcourt 1998; Blake et al. 1999), protein-induced deformability
(Olson et al. 1998), and duplex-free energy (Sugimoto et al.
1996), among others (for review, see Florquin et al. 2005), can be
used to describe (core) promoters, and to discriminate between
(core) promoter sequences and non (core) promoter sequences
(Ohler et al. 2001; Florquin et al. 2005; Kanhere and Bansal 2005;
Uren et al. 2006; Wang and Benham 2006). Because all these
properties are calculated from conversion tables using di- or tri-
nucleotides, one may argue that these properties are in fact ex-
actly the same as the nucleotide sequence and do not offer any
additional information. However, several studies have shown
that this is not the case. Both Liao et al. (2000) and Baldi et al.
(1998) have analyzed the correlation between the different prop-
erties, and their main conclusion was that the properties are
largely independent. Moreover, Florquin et al. (2005) have clus-
tered promoters based on these structural properties. The genes
associated with the promoters in each cluster varied greatly for
the different properties, which again indicates that the different
properties contain complementary information. Bode et al.
(2006) have shown that it is very hard to identify scaffold/matrix
attachment regions (SMARs) from the sequence, as the important
scaffold proteins recognize structural features instead of specific
nucleotide sequences. Structural properties are known to have
long-range interactions (up to 10 kb), so they can exhibit prop-
erties that are not visible in the sequence (Merling et al. 2003;
Faiger et al. 2006). The Human Genomic Melting Map (Liu et al.
2007) shows a correlation between GC content and DNA dena-
turation, but due to the cooperative nature of DNA denaturation,
this correlation is weaker on scales <500 bp.

We applied these different properties to the problem of pro-
moter prediction and show that some properties give better per-
formance than others, which again indicates differences in in-
formation content. In particular, we present the Easy Promoter
Prediction Program (EP3), which uses GC content and large-scale
structural features of DNA to identify and delineate promoter
regions in whole-genome sequences. EP3 was applied to the

human genome and compared with other state-of-the-art PPPs
using two different evaluation techniques: One is commonly
used for PPPs, while the other is based on a novel scheme also
considering intergenic predictions that do not fall in a promoter.
We evaluated the different ways of encoding the genomic DNA,
and we also applied EP3 to a set of human noncoding RNA genes.
Finally, we evaluated EP3 on 12 other eukaryotic genomes.

Results and Discussion

Properties of the core promoter region

Figure 1 shows examples of numerical profiles for human pro-
moter sequences. These profiles were obtained by lining up all
the promoter sequences with the TSS at the same position. Next,
we converted all sequences using the different structural proper-
ties of DNA into numeric sequences and, for each position, plot-
ted the average over all numeric sequences. The X-axis is the
position relative to the TSS, and the Y-axis is the normalized
value of that property. It is clear that the human core promoter
adopts a very specific intrinsic structure that stretches over quite
long distances. Either a large peak is visible, representing highly
stable regions in the DNA, or a cleft is visible, which represents
highly unstable regions in the DNA. For instance, the GC profile
shows a peak, while the AT profile shows a cleft. When zoomed
in, it can be noticed that, within the broader region of stability,
there is a small unstable region that corresponds to the TSS, and
a second one that corresponds to the region where TBP binds (see
below). The promoter prediction technique we present will
identify those higher regions and are, thus, most likely to contain
a TSS.

To see how general these properties of promoter sequences
are, we created large-scale (2000 bp up- and downstream of the
TSS) datasets for 16 different species. For each of these datasets,
we analyzed the structural profiles. Over long-range distances, we
observe several types of profiles (Fig. 2). The first one applies
mostly to protists. In this type, the profile slowly decreases a few
hundred base pairs upstream of the TSS, shows a transition of a
strong peak and cleft on the TSS, and is flat again downstream of
the TSS, which may indicate a relatively small promoter area in
these organisms. The second type of profile is observed in mam-
mals and shows a broad elevated region around the TSS, with
again a strong local peak and drop on the TSS itself. Between
these two extremes, we observe a more gradual transition in
fungi and plants. The extremely high and low values that can be
seen on the TSS are discussed in more detail in the next section.

There seems to be a relationship between the genome size
and the size of the peaks or clefts: In prokaryotes (not shown) and
protists, these cover ~100 bp, while they gradually increase to
several hundreds of base pairs in fungi and up to ~1000 bp in
plants. In animals, smaller genomes such as that of *Drosophila*
have a peak ranging from a few hundred base pairs on each side
of the TSS; this gradually increases up to ~2000 bp in mammals.
In conclusion, the peak is most obvious in mammals, but is also
clearly visible in fish, insects, fungi, and plants; it seems consid-
erably smaller in protists.

In humans, as well as in other eukaryotes, three types of
RNA polymerases (RNAP) exist that are used to transcribe differ-
ent types of genes. All three of them manage to identify the TSS
in a whole-genome setting in their own way. Here, we focus on
humans because this species is among the best documented ones.
We retrieved the promoter sequences for several types of genes
<figure><img src="image_1.png" /><figcaption>**Figure 1.** Examples of numerical profiles representing structural properties of the DNA in human over long-range distances around the transcription start site (TSS) (2000 bp upstream and 2000 bp downstream). The numeric profiles contain the actual average values on each position and are not smoothed using a window. Structural properties can be divided into two broad categories: those for which the numerical profile shows a peak around the TSS (e.g., GC content, bendability, and DNA denaturation), and those for which the profile shows a cleft around the TSS (e.g., AT content, bDNA twist, and Duplex Stability-Free Energy). (X-axis) Position relative to the TSS; (Y-axis) normalized value for each of the different properties. The profiles are based on the properties listed in the first column of Table 2. The coordinates for the TSS were retrieved from Ensembl, and the sequences were extracted around this TSS.</figcaption></figure>

from Ensembl (see Methods) and constructed three different datasets for each of the promoter types: RNAP I, RNAP II, and RNAP III. For each of these datasets, we calculated a structural profile. Figure 3 shows the inverted base-stacking values over long-range distances around the TSS (2000 bp up- and downstream), as well as a zoomed-in view (top row) of a shorter region around the TSS (ranging from 200 bp upstream of the TSS to 50 bp downstream of the TSS) for sets of promoters recognized by the three different polymerases. Both on a large and a small scale, there are clear and important fluctuations in the profile around the TSS. One of the most striking features, at least for RNAP II promoters, is the presence of two clefts: one at the TSS and one near position –30. The latter one is related to the TATA binding protein that is important for transcription initiation.

There is, however, a clear difference between RNAP II pro-
moters on the one hand, and RNAP I and III promoters on the
other hand. RNAP I and III promoters lack the large-scale stable
region of RNAP II promoters and show local destabilization only
around the TSS. These local instabilities around the TSS may be
caused partly by the structural nature of the transcript originat-
ing from these genes. Drawing more definitive conclusions for
the RNAP I and III promoters would require much more data
than the few hundred sequences that are currently available in
Ensembl. The lack of large-scale features is interesting as it may
indicate that the region involved in regulation of these gene
types is in general quite small. This might be in line with the fact
that many genes transcribed by RNAP I and III promoters are
produced in very large quantities and need less regulation. For
instance, rRNA genes account for 80% of total steady-state cellu-

lar RNA and might thus require less sophisticated regulation (Ja-
cob 1995). Another observation is that the regulation of r-protein
genes is often post-transcriptional and therefore does not need
extensive transcriptional regulation (Amaldi et al. 1989; Presutti
et al. 1991). RNAP II promoters have an elevated region of nearly
2000 bp, which indicates a much larger area involved in tran-
scriptional regulation of protein-coding genes.

Relationship between structural profiles and known core
promoter elements

To investigate the relationship between the structural profile and
the occurrence of known core promoter elements, we analyzed
four elements (TATA, INR, BRE, and CpG islands) well known in
animal genomes and the relation between the presence or ab-
sence of the element in the promoter sequence and the presence
or absence of peaks and clefts in the structural profile. Again, we
limited ourselves to humans because to detect motifs we require
position weight matrices (PWM) or consensus sequences, which
are not available for all organisms.

The first element we considered is the TATA box, also
known as the Hogness box, which binds the TATA binding pro-
tein (TBP) and is often involved in transcription initiation (Smale
and Kadonaga 2003). The structural profile for the promoters
containing the TATA motif is shown in Figure 4 (top row). As can
be noticed, the numerical profile of promoters that do contain
the TATA box lies lower than those that do not. Furthermore, the
cleft around position –30 is much deeper in TATA-box-
containing promoters, implying that the DNA denatures more
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Examples of numerical profiles representing structural properties of the DNA in human over long-range distances around the transcription start site (TSS) (2000 bp upstream and 2000 bp downstream). The numeric profiles contain the actual average values on each position and are not smoothed using a window. Structural properties can be divided into two broad categories: those for which the numerical profile shows a peak around the TSS (e.g., GC content, bendability, and DNA denaturation), and those for which the profile shows a cleft around the TSS (e.g., AT content, bDNA twist, and Duplex Stability-Free Energy). (X-axis) Position relative to the TSS; (Y-axis) normalized value for each of the different properties. The profiles are based on the properties listed in the first column of Table 2. The coordinates for the TSS were retrieved from Ensembl, and the sequences were extracted around this TSS.</figcaption></figure>

from Ensembl (see Methods) and constructed three different datasets for each of the promoter types: RNAP I, RNAP II, and RNAP III. For each of these datasets, we calculated a structural profile. Figure 3 shows the inverted base-stacking values over long-range distances around the TSS (2000 bp up- and downstream), as well as a zoomed-in view (top row) of a shorter region around the TSS (ranging from 200 bp upstream of the TSS to 50 bp downstream of the TSS) for sets of promoters recognized by the three different polymerases. Both on a large and a small scale, there are clear and important fluctuations in the profile around the TSS. One of the most striking features, at least for RNAP II promoters, is the presence of two clefts: one at the TSS and one near position –30. The latter one is related to the TATA binding protein that is important for transcription initiation.

There is, however, a clear difference between RNAP II pro-
moters on the one hand, and RNAP I and III promoters on the
other hand. RNAP I and III promoters lack the large-scale stable
region of RNAP II promoters and show local destabilization only
around the TSS. These local instabilities around the TSS may be
caused partly by the structural nature of the transcript originat-
ing from these genes. Drawing more definitive conclusions for
the RNAP I and III promoters would require much more data
than the few hundred sequences that are currently available in
Ensembl. The lack of large-scale features is interesting as it may
indicate that the region involved in regulation of these gene
types is in general quite small. This might be in line with the fact
that many genes transcribed by RNAP I and III promoters are
produced in very large quantities and need less regulation. For
instance, rRNA genes account for 80% of total steady-state cellu-

lar RNA and might thus require less sophisticated regulation (Ja-
cob 1995). Another observation is that the regulation of r-protein
genes is often post-transcriptional and therefore does not need
extensive transcriptional regulation (Amaldi et al. 1989; Presutti
et al. 1991). RNAP II promoters have an elevated region of nearly
2000 bp, which indicates a much larger area involved in tran-
scriptional regulation of protein-coding genes.

Relationship between structural profiles and known core
promoter elements

To investigate the relationship between the structural profile and the occurrence of known core promoter elements, we analyzed four elements (TATA, INR, BRE, and CpG islands) well known in animal genomes and the relation between the presence or absence of the element in the promoter sequence and the presence or absence of peaks and clefts in the structural profile. Again, we limited ourselves to humans because to detect motifs we require position weight matrices (PWM) or consensus sequences, which are not available for all organisms.

The first element we considered is the TATA box, also
known as the Hogness box, which binds the TATA binding pro-
tein (TBP) and is often involved in transcription initiation (Smale
and Kadonaga 2003). The structural profile for the promoters
containing the TATA motif is shown in Figure 4 (top row). As can
be noticed, the numerical profile of promoters that do contain
the TATA box lies lower than those that do not. Furthermore, the
cleft around position –30 is much deeper in TATA-box-
containing promoters, implying that the DNA denatures more
<figure><img src="image_2.png" /><figcaption>**Figure 2.** Examples of numerical profiles representing different species using the same property (inverted base stacking value) over long-range distances around the transcription start site (TSS) (2000 bp upstream and 2000 bp downstream). (X-axis) Position relative to the TSS; (Y-axis) normalized value of the property. The profile is calculated with a window size of three.</figcaption></figure>

easily on that position, which conforms to the knowledge that the TBP is involved in the denaturation of DNA. Although it is clear that promoters containing the TATA motif have a much deeper cleft than those without, we must note that the cleft is visible as well. This indicates that the TBP, which is important for transcription to start, requires a region that denatures more easily, whether it contains the motif or not, to allow efficient binding (Comai et al. 1992; Cormack and Struhl 1992; White and Jackson 1992). Therefore, we can conclude that the TATA motif is responsible for the instability of this region, but other combinations of nucleotides also can result in a region that is unstable enough. Also, on a larger scale (Fig. 4, right side), we see that the

profile is slightly lower for promoters containing TATA than for non-TATA promoters.

A second well-known core promoter element is the Initiator element (INR) that is found mostly on position 0; i.e., the TSS itself. In Figure 4 (second row), the structural profiles for promoters with and without the INR motif are displayed. Overall, the two profiles are very similar, except for one remarkable difference visible on the small-scale plot: On position -1 there is a significant peak (up to -0.05) upstream of the TSS cleft in the INR-containing promoters, while this peak is missing in the promoters that lack INR. This seems to indicate that this peak, which increases the stability-instability transition, plays a role in the accurate transcription initiation because the INR motif is known to be important for accurate transcription initiation (Lo and Smale 1996). On a larger scale, we do not observe a significant difference between INR-containing promoters and those that do not.

The TFIIB recognition element (BRE) consists of two parts, BRE<sup>u</sup>, which is located upstream of the position of the TATA box (Lagrange et al. 1998), and BRE<sup>d</sup>, which is located downstream of the TATA motif (Deng and Roberts 2005). The profile of BRE promoters is plotted in the third row of Figure 4. It is clear that the profile of promoters that contain a BRE motif is higher than that of promoters lacking the BRE motif. This is not surprising because the BRE motif consists entirely of G and C nucleotides, which results in higher inverted base-stacking values. The local profile for BRE-containing promoters shows a rather stable area from position -40 to -20. These features fit well with the observation that TFIIB (GTF2B) interacts with the major groove upstream of the TATA box and with the minor groove downstream of the TATA box (Nikolov et al. 1995; Lagrange et al. 1998; Deng and Roberts 2005, 2006). Both interactions leave their tracks on the structural profile. The interaction with the major groove is visible as a peak at position -35, while the interaction with the minor groove is visible as a peak at position -25. These two peaks form a cleft at position -30 that will be used by TBP to bind. The BRE motif is more prevalent among promoters without the TATA motif, and the motif might be a functional substitute for the TATA box (Deng and Roberts 2006). On a larger scale, we see that the profile of promoters lacking the BRE motif has a much lower amplitude than that of promoters containing the BRE motif.

A last well-known characteristic of core promoters is the abundance of CpG islands near protein-coding genes in mam-
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Examples of numerical profiles representing different species using the same property (inverted base stacking value) over long-range distances around the transcription start site (TSS) (2000 bp upstream and 2000 bp downstream). (X-axis) Position relative to the TSS; (Y-axis) normalized value of the property. The profile is calculated with a window size of three.</figcaption></figure>

easily on that position, which conforms to the knowledge that the TBP is involved in the denaturation of DNA. Although it is clear that promoters containing the TATA motif have a much deeper cleft than those without, we must note that the cleft is visible as well. This indicates that the TBP, which is important for transcription to start, requires a region that denatures more easily, whether it contains the motif or not, to allow efficient binding (Comai et al. 1992; Cormack and Struhl 1992; White and Jackson 1992). Therefore, we can conclude that the TATA motif is responsible for the instability of this region, but other combinations of nucleotides also can result in a region that is unstable enough. Also, on a larger scale (Fig. 4, right side), we see that the

profile is slightly lower for promoters containing TATA than for non-TATA promoters.

A second well-known core promoter element is the Initiator element (INR) that is found mostly on position 0; i.e., the TSS itself. In Figure 4 (second row), the structural profiles for promoters with and without the INR motif are displayed. Overall, the two profiles are very similar, except for one remarkable difference visible on the small-scale plot: On position -1 there is a significant peak (up to -0.05) upstream of the TSS cleft in the INR-containing promoters, while this peak is missing in the promoters that lack INR. This seems to indicate that this peak, which increases the stability-instability transition, plays a role in the accurate transcription initiation because the INR motif is known to be important for accurate transcription initiation (Lo and Smale 1996). On a larger scale, we do not observe a significant difference between INR-containing promoters and those that do not.

The TFIIB recognition element (BRE) consists of two parts, BRE<sup>u</sup>, which is located upstream of the position of the TATA box (Lagrange et al. 1998), and BRE<sup>d</sup>, which is located downstream of the TATA motif (Deng and Roberts 2005). The profile of BRE promoters is plotted in the third row of Figure 4. It is clear that the profile of promoters that contain a BRE motif is higher than that of promoters lacking the BRE motif. This is not surprising because the BRE motif consists entirely of G and C nucleotides, which results in higher inverted base-stacking values. The local profile for BRE-containing promoters shows a rather stable area from position -40 to -20. These features fit well with the observation that TFIIB (GTF2B) interacts with the major groove upstream of the TATA box and with the minor groove downstream of the TATA box (Nikolov et al. 1995; Lagrange et al. 1998; Deng and Roberts 2005, 2006). Both interactions leave their tracks on the structural profile. The interaction with the major groove is visible as a peak at position -35, while the interaction with the minor groove is visible as a peak at position -25. These two peaks form a cleft at position -30 that will be used by TBP to bind. The BRE motif is more prevalent among promoters without the TATA motif, and the motif might be a functional substitute for the TATA box (Deng and Roberts 2006). On a larger scale, we see that the profile of promoters lacking the BRE motif has a much lower amplitude than that of promoters containing the BRE motif.

A last well-known characteristic of core promoters is the abundance of CpG islands near protein-coding genes in mam-
<figure><img src="image_3.png" /><figcaption>**Figure 3.** Structural profile of promoters for genes transcribed by different polymerases. (Top row) Numerical profiles over short-range distances around the TSS; (bottom row) profiles over long-range distances. The profile was calculated using the inverse of the base-stacking property. See text for details.</figcaption></figure>

mals (Bird 2002). CpG islands are regions that are mainly associated with mammalian promoters, but are also observed in plants (Rombauts et al. 2003). They are related to housekeeping genes (Baek et al. 2007), dynamic usage of TSSs (Kawaji et al. 2006), and bidirectional promoter activity (Trinklein et al. 2004). The CpG dinucleotide is underrepresented in the genome because the cytosine can be converted into thymine after methylation. Only DNA stretches that are under evolutionary stress, such as promoters, are preserved, and these regions are rich in CpG dinucleotides that are unmethylated. The unmethylated CpG-rich regions are often associated with promoters of protein-

coding genes (Bird 2002). These islands are typically found upstream of the core promoter, and, when present, these promoters usually lack a TATA box and often have alternative start sites. The local profile does not vary much between promoters having CpG islands and the others, with only an expected shift toward higher stability values due to GC richness. On a larger scale, as with BRE elements, the profile of promoters lacking CpG islands has a much lower amplitude than that of promoters with CpG islands. This may indicate that the existence of either feature is correlated with long-range regulation, which would require further investigation.

<figure><img src="image_4.png" /><figcaption>**Figure 4.** Structural profiles of RNAP II promoters containing known motifs or elements versus promoters for which the presence of motifs cannot be demonstrated. (Left two columns) Short-range distances around the TSS; (right two columns) long-range distances around the TSS. The profile was calculated using the inverse of the base-stacking property. (Inr) Initiator element, (BR) TFIIB recognition element.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Structural profile of promoters for genes transcribed by different polymerases. (Top row) Numerical profiles over short-range distances around the TSS; (bottom row) profiles over long-range distances. The profile was calculated using the inverse of the base-stacking property. See text for details.</figcaption></figure>

mals (Bird 2002). CpG islands are regions that are mainly associated with mammalian promoters, but are also observed in plants (Rombauts et al. 2003). They are related to housekeeping genes (Baek et al. 2007), dynamic usage of TSSs (Kawaji et al. 2006), and bidirectional promoter activity (Trinklein et al. 2004). The CpG dinucleotide is underrepresented in the genome because the cytosine can be converted into thymine after methylation. Only DNA stretches that are under evolutionary stress, such as promoters, are preserved, and these regions are rich in CpG dinucleotides that are unmethylated. The unmethylated CpG-rich regions are often associated with promoters of protein-

coding genes (Bird 2002). These islands are typically found upstream of the core promoter, and, when present, these promoters usually lack a TATA box and often have alternative start sites. The local profile does not vary much between promoters having CpG islands and the others, with only an expected shift toward higher stability values due to GC richness. On a larger scale, as with BRE elements, the profile of promoters lacking CpG islands has a much lower amplitude than that of promoters with CpG islands. This may indicate that the existence of either feature is correlated with long-range regulation, which would require further investigation.

<figure><img src="image_5.png" /><figcaption>**Figure 4.** Structural profiles of RNAP II promoters containing known motifs or elements versus promoters for which the presence of motifs cannot be demonstrated. (Left two columns) Short-range distances around the TSS; (right two columns) long-range distances around the TSS. The profile was calculated using the inverse of the base-stacking property. (Inr) Initiator element, (BR) TFIIB recognition element.</figcaption></figure>
## Promoter prediction and comparison to the state of the art

We have shown that in eukaryotes the region near the TSS has a very distinct structural profile, while this is absent in the remainder of intergenic sequences and in coding sequences (Fig. 1). The peak in this profile, for instance using base stacking values, is a well-positioned property of the eukaryotic promoter that stretches over several hundreds of base pairs with its top located on the TSS, which can be used to identify promoters. As computation of this profile for individual sequences involves smoothing over a few hundred base pairs, this feature is too coarse to predict the actual TSS; it is nevertheless very well suited to predict the core promoter region. Furthermore, due to the fact that it is such a large-scale property, there are very few false positives when scanning the genome. Using large-scale structural features is new to the field of promoter prediction, because previous studies usually analyzed no more than 200 bp around the TSS (Ohler and Niemann 2001; Fukue et al. 2004, 2005; Florquin et al. 2005; Kanhere and Bansal 2005).

In applying EP3 to the human genome, we observed clear peaks near the start of most genes, which is exemplified in Figure 5 for a small region of the human genome. This figure shows a plot of the structural profile for 2 Mbp on the human chromosome 21 using DNA base-stacking values and a window size of 400 bp (see Methods). To demonstrate the predictive strength of these structural profiles, we have added the gene annotation retrieved from Ensembl to the predictions. This figure clearly shows that most genes do have a peak near the 5' end, and that there are virtually no peaks outside this region (false positives), indicating that the large-scale structural properties we observe

are distinct enough for EP3 to efficiently predict core promoter regions. It should be noted that, while this single property proves to work surprisingly well, future research may focus on combining multiple properties to describe the promoter region even better, as the overlap between the predictions obtained through different properties may vary significantly (Supplemental Table 1). For example, the overlap between Bendability and Duplex Stability-Free energy is only 1%, while the overlap between predictions from Duplex Stability-Free energy and DNA denaturation is ≥80%. The average overlap between a pair of properties is 44% (SD = 29%). Improving the models for the different structural properties may also increase the overall performance.

Next, we compared our program with a broad range of programs that are currently available for human promoter predictions (Table 1). We used both the Ensembl gene annotation and the CAGE transcription start data and tested GC content as the simplest feature and compared it with the other structural features. All programs were run using default settings unless indicated otherwise. For those programs that depend on a parameter, we included the best result in the table. For EP3, we included the values for the base-stacking property, as these gave the best results (Table 2). To rank the different programs, we used the F-measure of the program with a maximum allowed distance of 500 bp with the true TSS compared with the CAGE dataset. This is in contrast to other recent reviews and promoter prediction papers that used less strict validation (Bajic et al. 2004, 2006a; Sonnenburg et al. 2006; Xie et al. 2006). We included the less strict mismatch distances (1000 and 2000 bp) for reference with the other papers that allowed a distance of 1000 bp (Bajic et al. 2006a) or 2000 bp (Bajic et al. 2004; Sonnenburg et al. 2006).

When using larger values, more distant predictions will also be considered to be true positive (TP). This is not desirable, as predictions should be as close as possible to the TSS. The F-measure clearly decreases when the evaluation becomes stricter, indicating that all programs output predictions that are quite distant from the actual TSS.

For NNPP2.2, McPromoter, and Promoter2.0, we considered different parameter settings. We included only the setting that resulted in the highest F-measure. We considered two different techniques to calculate the number of TPs, false positives (FNs), and false negatives (FNs) (see Methods). The first one is based simply on gene annotation and was used before in several studies on promoter predictions (Bajic et al. 2004, 2006a; Sonnenburg et al. 2006; Xie et al. 2006). The second one is an alternative technique we developed that also includes intergenic predictions but requires transcriptional data, such as the CAGE dataset.

For the Ensembl annotation, we see that Dragon Gene Start Finder (GSF) (F = 0.53) and PromoterInspector (F = 0.49) perform better than EP3 (F = 44). However, when using the CAGE annotation, EP3 performs slightly

<figure><img src="image_5.png" /><figcaption>**Figure 5.** Structural profile (blue) of human chromosome 21 between position 32,000,000 bp and 33,000,000 bp. The profile, based on inverted base-stacking values, was made using a window size of 400 bp and using nonoverlapping windows. Experimentally annotated genes from Ensembl are indicated: (yellow bands) positive strand; (blue bands) negative strand. Peaks in the profile that exceed the classification threshold (T = 0.19, horizontal gray line) are predicted as promoter regions. (X-axis) Position in the genomic sequence, (Y-axis) inverse of the value for base-stacking energy. The labels on the plot show the gene names as they appear in Ensembl.</figcaption></figure>
## Promoter prediction and comparison to the state of the art

We have shown that in eukaryotes the region near the TSS has a very distinct structural profile, while this is absent in the remainder of intergenic sequences and in coding sequences (Fig. 1). The peak in this profile, for instance using base stacking values, is a well-positioned property of the eukaryotic promoter that stretches over several hundreds of base pairs with its top located on the TSS, which can be used to identify promoters. As computation of this profile for individual sequences involves smoothing over a few hundred base pairs, this feature is too coarse to predict the actual TSS; it is nevertheless very well suited to predict the core promoter region. Furthermore, due to the fact that it is such a large-scale property, there are very few false positives when scanning the genome. Using large-scale structural features is new to the field of promoter prediction, because previous studies usually analyzed no more than 200 bp around the TSS (Ohler and Niemann 2001; Fukue et al. 2004, 2005; Florquin et al. 2005; Kanhere and Bansal 2005).

In applying EP3 to the human genome, we observed clear peaks near the start of most genes, which is exemplified in Figure 5 for a small region of the human genome. This figure shows a plot of the structural profile for 2 Mbp on the human chromosome 21 using DNA base-stacking values and a window size of 400 bp (see Methods). To demonstrate the predictive strength of these structural profiles, we have added the gene annotation retrieved from Ensembl to the predictions. This figure clearly shows that most genes do have a peak near the 5' end, and that there are virtually no peaks outside this region (false positives), indicating that the large-scale structural properties we observe

are distinct enough for EP3 to efficiently predict core promoter regions. It should be noted that, while this single property proves to work surprisingly well, future research may focus on combining multiple properties to describe the promoter region even better, as the overlap between the predictions obtained through different properties may vary significantly (Supplemental Table 1). For example, the overlap between Bendability and Duplex Stability-Free energy is only 1%, while the overlap between predictions from Duplex Stability-Free energy and DNA denaturation is ≥80%. The average overlap between a pair of properties is 44% (SD = 29%). Improving the models for the different structural properties may also increase the overall performance.

Next, we compared our program with a broad range of programs that are currently available for human promoter predictions (Table 1). We used both the Ensembl gene annotation and the CAGE transcription start data and tested GC content as the simplest feature and compared it with the other structural features. All programs were run using default settings unless indicated otherwise. For those programs that depend on a parameter, we included the best result in the table. For EP3, we included the values for the base-stacking property, as these gave the best results (Table 2). To rank the different programs, we used the *F*-measure of the program with a maximum allowed distance of 500 bp with the true TSS compared with the CAGE dataset. This is in contrast to other recent reviews and promoter prediction papers that used less strict validation (Bajic et al. 2004, 2006a; Sonnenburg et al. 2006; Xie et al. 2006). We included the less strict mismatch distances (1000 and 2000 bp) for reference with the other papers that allowed a distance of 1000 bp (Bajic et al. 2006a) or 2000 bp (Bajic et al. 2004; Sonnenburg et al. 2006).

When using larger values, more distant predictions will also be considered to be true positive (TP). This is not desirable, as predictions should be as close as possible to the TSS. The *F*-measure clearly decreases when the evaluation becomes stricter, indicating that all programs output predictions that are quite distant from the actual TSS.

For NNPP2.2, McPromoter, and Promoter2.0, we considered different parameter settings. We included only the setting that resulted in the highest *F*-measure. We considered two different techniques to calculate the number of TPs, false positives (FPs), and false negatives (FNs) (see Methods). The first one is based simply on gene annotation and was used before in several studies on promoter predictions (Bajic et al. 2004, 2006a; Sonnenburg et al. 2006; Xie et al. 2006). The second one is an alternative technique we developed that also includes intergenic predictions but requires transcriptional data, such as the CAGE dataset.

For the Ensembl annotation, we see that Dragon Gene Start Finder (GSF) (*F* = 0.53) and PromoterInspector (*F* = 0.49) perform better than EP3 (*F* = 44). However, when using the CAGE annotation, EP3 performs slightly

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Structural profile (blue) of human chromosome 21 between position 32,000,000 bp and 33,000,000 bp. The profile, based on inverted base-stacking values, was made using a window size of 400 bp and using nonoverlapping windows. Experimentally annotated genes from Ensembl are indicated: (yellow bands) positive strand; (blue bands) negative strand. Peaks in the profile that exceed the classification threshold (*T* = 0.19, horizontal gray line) are predicted as promoter regions. (X-axis) Position in the genomic sequence, (Y-axis) inverse of the value for base-stacking energy. The labels on the plot show the gene names as they appear in Ensembl.</figcaption></figure>
<table><caption>**Table 1.** The performance of all current promoter prediction programs capable of analyzing the whole human genome</caption><thead><tr><th rowspan="3">Program</th><th rowspan="3">Reference</th><th colspan="10">Ensembl</th><th colspan="6">CAGE</th></tr><tr><th colspan="2">500</th><th colspan="2">1000</th><th colspan="2">2000</th><th colspan="2">500</th><th colspan="2">1000</th><th colspan="2">2000</th></tr><tr><th>Recall</th><th>Prec.</th><th>F</th><th>Recall</th><th>Prec.</th><th>F</th><th>Recall</th><th>Prec.</th><th>F</th><th>Recall</th><th>Prec.</th><th>F</th><th>Recall</th><th>Prec.</th><th>F</th></tr></thead><tbody><tr><td>EP3(Base stacking)</td><td>Bajic et al. (2003)</td><td>0.42</td><td>0.46</td><td>0.44</td><td>0.46</td><td>0.56</td><td>0.51</td><td>0.49</td><td>0.64</td><td>0.56</td><td>0.34</td><td>0.66</td><td><strong>0.45</strong></td><td>0.38</td><td>0.72</td><td>0.50</td><td>0.41</td><td>0.76</td><td>0.53</td></tr><tr><td>DragonGSF</td><td>Scherf et al. (2000)</td><td>0.45</td><td>0.63</td><td>0.53</td><td>0.50</td><td>0.74</td><td>0.60</td><td>0.54</td><td>0.80</td><td>0.64</td><td>0.31</td><td>0.75</td><td><strong>0.44</strong></td><td>0.37</td><td>0.81</td><td>0.51</td><td>0.41</td><td>0.85</td><td>0.55</td></tr><tr><td>PromoterInspector</td><td>Davuluri et al. (2001)</td><td>0.38</td><td>0.70</td><td>0.49</td><td>0.42</td><td>0.78</td><td>0.55</td><td>0.44</td><td>0.82</td><td>0.57</td><td>0.29</td><td>0.81</td><td><strong>0.43</strong></td><td>0.34</td><td>0.85</td><td>0.48</td><td>0.36</td><td>0.87</td><td>0.51</td></tr><tr><td>FirstEF</td><td>Down and Hubbard (2002)</td><td>0.58</td><td>0.34</td><td>0.43</td><td>0.64</td><td>0.39</td><td>0.48</td><td>0.69</td><td>0.43</td><td>0.53</td><td>0.41</td><td>0.42</td><td><strong>0.42</strong></td><td>0.48</td><td>0.48</td><td>0.48</td><td>0.54</td><td>0.53</td><td>0.53</td></tr><tr><td>Eponine</td><td>Gross and Brent (2006)</td><td>0.36</td><td>0.51</td><td>0.42</td><td>0.38</td><td>0.65</td><td>0.48</td><td>0.40</td><td>0.74</td><td>0.52</td><td>0.28</td><td>0.75</td><td><strong>0.41</strong></td><td>0.32</td><td>0.80</td><td>0.45</td><td>0.34</td><td>0.84</td><td>0.48</td></tr><tr><td>N-Scan</td><td>Ponger and Mouchiroud (2002)</td><td>0.55</td><td>0.51</td><td>0.53</td><td>0.60</td><td>0.55</td><td>0.57</td><td>0.63</td><td>0.58</td><td>0.60</td><td>0.33</td><td>0.45</td><td><strong>0.38</strong></td><td>0.39</td><td>0.50</td><td>0.44</td><td>0.43</td><td>0.54</td><td>0.48</td></tr><tr><td>CpgProD</td><td>Xie et al. (2006)</td><td>0.50</td><td>0.36</td><td>0.42</td><td>0.59</td><td>0.42</td><td>0.49</td><td>0.65</td><td>0.46</td><td>0.54</td><td>0.34</td><td>0.41</td><td><strong>0.37</strong></td><td>0.44</td><td>0.49</td><td>0.46</td><td>0.51</td><td>0.54</td><td>0.52</td></tr><tr><td>PromoterExplorer</td><td>OHler et al. (2000)</td><td>0.55</td><td>0.24</td><td>0.33</td><td>0.66</td><td>0.28</td><td>0.39</td><td>0.73</td><td>0.32</td><td>0.44</td><td>0.39</td><td>0.30</td><td><strong>0.34</strong></td><td>0.50</td><td>0.37</td><td>0.42</td><td>0.60</td><td>0.42</td><td>0.49</td></tr><tr><td>McPromoter (0.0)</td><td>Chen et al. (1997)</td><td>0.24</td><td>0.61</td><td>0.34</td><td>0.29</td><td>0.71</td><td>0.41</td><td>0.32</td><td>0.77</td><td>0.45</td><td>0.17</td><td>0.68</td><td><strong>0.28</strong></td><td>0.23</td><td>0.76</td><td>0.35</td><td>0.26</td><td>0.80</td><td>0.39</td></tr><tr><td>PromfD</td><td>Bajic et al. (2002)</td><td>0.55</td><td>0.14</td><td>0.22</td><td>0.61</td><td>0.16</td><td>0.25</td><td>0.68</td><td>0.18</td><td>0.28</td><td>0.44</td><td>0.16</td><td><strong>0.23</strong></td><td>0.50</td><td>0.18</td><td>0.27</td><td>0.60</td><td>0.21</td><td>0.32</td></tr><tr><td>DragonPF</td><td>Sonnenburg et al. (2006)</td><td>0.65</td><td>0.11</td><td>0.19</td><td>0.71</td><td>0.13</td><td>0.22</td><td>0.79</td><td>0.16</td><td>0.27</td><td>0.51</td><td>0.11</td><td><strong>0.18</strong></td><td>0.59</td><td>0.13</td><td>0.22</td><td>0.68</td><td>0.16</td><td>0.26</td></tr><tr><td>ARTS</td><td>Prestridge (1995)</td><td>0.77</td><td>0.08</td><td>0.14</td><td>0.82</td><td>0.09</td><td>0.16</td><td>0.88</td><td>0.10</td><td>0.18</td><td>0.63</td><td>0.06</td><td><strong>0.12</strong></td><td>0.71</td><td>0.08</td><td>0.14</td><td>0.80</td><td>0.09</td><td>0.16</td></tr><tr><td>PromoterScan</td><td>Knudsen (1999)</td><td>0.19</td><td>0.08</td><td>0.11</td><td>0.27</td><td>0.11</td><td>0.16</td><td>0.36</td><td>0.15</td><td>0.21</td><td>0.16</td><td>0.09</td><td><strong>0.12</strong></td><td>0.24</td><td>0.14</td><td>0.18</td><td>0.33</td><td>0.19</td><td>0.24</td></tr><tr><td>Promoter2.0 (medium)</td><td>Reese (2001)</td><td>0.68</td><td>0.03</td><td>0.06</td><td>0.91</td><td>0.04</td><td>0.08</td><td>0.99</td><td>0.04</td><td>0.08</td><td>0.63</td><td>0.04</td><td><strong>0.08</strong></td><td>0.87</td><td>0.06</td><td>0.11</td><td>0.97</td><td>0.07</td><td>0.12</td></tr><tr><td>NNPP 2.2(0.99)</td><td></td><td>0.03</td><td>0.02</td><td>0.02</td><td>0.06</td><td>0.04</td><td>0.05</td><td>0.10</td><td>0.08</td><td>0.09</td><td>0.03</td><td>0.03</td><td><strong>0.03</strong></td><td>0.05</td><td>0.06</td><td>0.05</td><td>0.09</td><td>0.11</td><td>0.10</td></tr></tbody></table>

The F-measure of the CAGE dataset with strictest maximum allowed mismatch was used to rank the programs (in bold). The FProm program (Solovyev et al. 2006) was omitted from the analyses because it is not for free for academic use. Prec., precision.
<table><caption>**Table 1.** The performance of all current promoter prediction programs capable of analyzing the whole human genome</caption><thead><tr><th rowspan="3">Program</th><th rowspan="3">Reference</th><th colspan="10">Ensembl</th><th colspan="6">CAGE</th></tr><tr><th colspan="2">500</th><th colspan="2">1000</th><th colspan="2">2000</th><th colspan="2">500</th><th colspan="2">1000</th><th colspan="2">2000</th></tr><tr><th>Recall</th><th>Prec.</th><th>F</th><th>Recall</th><th>Prec.</th><th>F</th><th>Recall</th><th>Prec.</th><th>F</th><th>Recall</th><th>Prec.</th><th>F</th><th>Recall</th><th>Prec.</th><th>F</th></tr></thead><tbody><tr><td>EP3(Base stacking)</td><td>Bajic et al. (2003)</td><td>0.42</td><td>0.46</td><td>0.44</td><td>0.46</td><td>0.56</td><td>0.51</td><td>0.49</td><td>0.64</td><td>0.56</td><td>0.34</td><td>0.66</td><td><strong>0.45</strong></td><td>0.38</td><td>0.72</td><td>0.50</td><td>0.41</td><td>0.76</td><td>0.53</td></tr><tr><td>DragonGSF</td><td>Scherf et al. (2000)</td><td>0.45</td><td>0.63</td><td>0.53</td><td>0.50</td><td>0.74</td><td>0.60</td><td>0.54</td><td>0.80</td><td>0.64</td><td>0.31</td><td>0.75</td><td><strong>0.44</strong></td><td>0.37</td><td>0.81</td><td>0.51</td><td>0.41</td><td>0.85</td><td>0.55</td></tr><tr><td>PromoterInspector</td><td>Davuluri et al. (2001)</td><td>0.38</td><td>0.70</td><td>0.49</td><td>0.42</td><td>0.78</td><td>0.55</td><td>0.44</td><td>0.82</td><td>0.57</td><td>0.29</td><td>0.81</td><td><strong>0.43</strong></td><td>0.34</td><td>0.85</td><td>0.48</td><td>0.36</td><td>0.87</td><td>0.51</td></tr><tr><td>FirstEF</td><td>Down and Hubbard (2002)</td><td>0.58</td><td>0.34</td><td>0.43</td><td>0.64</td><td>0.39</td><td>0.48</td><td>0.69</td><td>0.43</td><td>0.53</td><td>0.41</td><td>0.42</td><td><strong>0.42</strong></td><td>0.48</td><td>0.48</td><td>0.48</td><td>0.54</td><td>0.53</td><td>0.53</td></tr><tr><td>Eponine</td><td>Gross and Brent (2006)</td><td>0.36</td><td>0.51</td><td>0.42</td><td>0.38</td><td>0.65</td><td>0.48</td><td>0.40</td><td>0.74</td><td>0.52</td><td>0.28</td><td>0.75</td><td><strong>0.41</strong></td><td>0.32</td><td>0.80</td><td>0.45</td><td>0.34</td><td>0.84</td><td>0.48</td></tr><tr><td>N-Scan</td><td>Ponger and Mouchiroud (2002)</td><td>0.55</td><td>0.51</td><td>0.53</td><td>0.60</td><td>0.55</td><td>0.57</td><td>0.63</td><td>0.58</td><td>0.60</td><td>0.33</td><td>0.45</td><td><strong>0.38</strong></td><td>0.39</td><td>0.50</td><td>0.44</td><td>0.43</td><td>0.54</td><td>0.48</td></tr><tr><td>CpgProD</td><td>Xie et al. (2006)</td><td>0.50</td><td>0.36</td><td>0.42</td><td>0.59</td><td>0.42</td><td>0.49</td><td>0.65</td><td>0.46</td><td>0.54</td><td>0.34</td><td>0.41</td><td><strong>0.37</strong></td><td>0.44</td><td>0.49</td><td>0.46</td><td>0.51</td><td>0.54</td><td>0.52</td></tr><tr><td>PromoterExplorer</td><td>OHler et al. (2000)</td><td>0.55</td><td>0.24</td><td>0.33</td><td>0.66</td><td>0.28</td><td>0.39</td><td>0.73</td><td>0.32</td><td>0.44</td><td>0.39</td><td>0.30</td><td><strong>0.34</strong></td><td>0.50</td><td>0.37</td><td>0.42</td><td>0.60</td><td>0.42</td><td>0.49</td></tr><tr><td>McPromoter (0.0)</td><td>Chen et al. (1997)</td><td>0.24</td><td>0.61</td><td>0.34</td><td>0.29</td><td>0.71</td><td>0.41</td><td>0.32</td><td>0.77</td><td>0.45</td><td>0.17</td><td>0.68</td><td><strong>0.28</strong></td><td>0.23</td><td>0.76</td><td>0.35</td><td>0.26</td><td>0.80</td><td>0.39</td></tr><tr><td>PromF</td><td>Bajic et al. (2002)</td><td>0.55</td><td>0.14</td><td>0.22</td><td>0.61</td><td>0.16</td><td>0.25</td><td>0.68</td><td>0.18</td><td>0.28</td><td>0.44</td><td>0.16</td><td><strong>0.23</strong></td><td>0.50</td><td>0.18</td><td>0.27</td><td>0.60</td><td>0.21</td><td>0.32</td></tr><tr><td>DragonPF</td><td>Sonnenburg et al. (2006)</td><td>0.65</td><td>0.11</td><td>0.19</td><td>0.71</td><td>0.13</td><td>0.22</td><td>0.79</td><td>0.16</td><td>0.27</td><td>0.51</td><td>0.11</td><td><strong>0.18</strong></td><td>0.59</td><td>0.13</td><td>0.22</td><td>0.68</td><td>0.16</td><td>0.26</td></tr><tr><td>ARTS</td><td>Prestridge (1995)</td><td>0.77</td><td>0.08</td><td>0.14</td><td>0.82</td><td>0.09</td><td>0.16</td><td>0.88</td><td>0.10</td><td>0.18</td><td>0.63</td><td>0.06</td><td><strong>0.12</strong></td><td>0.71</td><td>0.08</td><td>0.14</td><td>0.80</td><td>0.09</td><td>0.16</td></tr><tr><td>PromoterScan</td><td>Knudsen (1999)</td><td>0.19</td><td>0.08</td><td>0.11</td><td>0.27</td><td>0.11</td><td>0.16</td><td>0.36</td><td>0.15</td><td>0.21</td><td>0.16</td><td>0.09</td><td><strong>0.12</strong></td><td>0.24</td><td>0.14</td><td>0.18</td><td>0.33</td><td>0.19</td><td>0.24</td></tr><tr><td>Promoter2.0 (medium)</td><td>Reese (2001)</td><td>0.68</td><td>0.03</td><td>0.06</td><td>0.91</td><td>0.04</td><td>0.08</td><td>0.99</td><td>0.04</td><td>0.08</td><td>0.63</td><td>0.04</td><td><strong>0.08</strong></td><td>0.87</td><td>0.06</td><td>0.11</td><td>0.97</td><td>0.07</td><td>0.12</td></tr><tr><td>NNPP 2.2(0.99)</td><td></td><td>0.03</td><td>0.02</td><td>0.02</td><td>0.06</td><td>0.04</td><td>0.05</td><td>0.10</td><td>0.08</td><td>0.09</td><td>0.03</td><td>0.03</td><td><strong>0.03</strong></td><td>0.05</td><td>0.06</td><td>0.05</td><td>0.09</td><td>0.11</td><td>0.10</td></tr></tbody></table>

The F-measure of the CAGE dataset with strictest maximum allowed mismatch was used to rank the programs (in bold). The FProm program (Solovyev et al. 2006) was omitted from the analyses because it is not for free for academic use. Prec., precision.
better (F = 0.45) than DragonGSF (F = 0.44) and PromoterInspector (F = 0.43). Although these differences seem small, it still concerns several hundreds of genes or TSSs. For instance, a 1% lower recall means 200 genes in the Ensembl dataset that are missed by that program. One percent less precision means 200 extra false predictions, so even the small differences are significant. In the case of the CAGE dataset, 1% difference corresponds to a difference of ~1250 TSSs. Combining this into the F-measure means that a program that has 1% lower value misses 1250 TSSs and predicts an additional 1250 wrong ones.

In conclusion, the best performing programs are EP3, DragonGSF, PromoterInspector, Eponine, and FirstEF, which have an F-measure > 0.4. Several programs (CpGProD, PromoterExplorer, N-Scan, and McPromoter) still perform quite well and have an F-measure > 0.25. The rest of the programs (PromFD, ARTS, DragonPF, PromoterScan, NNPP2.2, and Promoter2.0) do not perform very well on the complete human genome (F < 0.25). In all cases, these low F-measures are caused by a very low precision, which is obtained when the PPP outputs many FPs, a problem reminiscent of promoter prediction since the beginning.

The performance of PPPs on the Ensembl dataset is generally higher than that on the CAGE dataset, probably because of the different counting schemes for Ensembl data and the CAGE dataset, whereas the scheme for Ensembl data ignores all intergenic predictions. Some of the other programs may also perform better on the Ensembl data than EP3, because these programs have been trained on promoters of protein-coding genes and are therefore more gene-centric. For the top-performing programs, the balance between the recall (sensitivity) and precision (specificity) is mostly favoring the precision, with the exception of FirstEF, which is perfectly balanced on the CAGE dataset. The other programs often have high recall values, but at the cost of very low precision values.

## Performance on the ENCODE region

The ENCODE project aims to carefully annotate all functional elements in a small portion (1%) of the human genome. We used EP3 to make predictions on the 44 regions that cover ~30 Mb and compared the prediction with three different datasets. First, we compared our predictions to a set of known functional promoters (Cooper et al. 2006). The dataset is only partial for the ENCODE region because Cooper et al. tested only 642 putative promoters of the 921 they predicted. These partial data are insufficient to assess the precision of EP3. Of the 642 promoters tested, 387 were discovered to be functional. EP3 predicts 24% (recall) of the promoters that are marked as functional by Cooper et al. when using a maximum distance of 500 bp. This low recall rate is rather surprising, as the performance on the GENCODE and CAGE data is very good (see below). But when we look in detail, a much higher recall rate (40%) is obtained for the genes expressed in all 16 cell lines, which indicates that EP3 is biased toward broadly expressed genes. Of the 257 nonfunctional promoters, 18 (7%) are predicted by EP3.

Next, we compared the predictions of EP3 on the ENCODE region with the gene annotation from the GENCODE project (Harrow et al. 2006) and with the CAGE data from Riken. The GENCODE annotation was compared with the predictions using the classic method for calculating the performance, as it is a gene annotation similar to the one of Ensembl. We obtained a recall of 0.46, a precision of 0.72, and an F-measure of 0.56. The performance on the CAGE data was calculated with the novel method

presented here and has 0.61 recall, 0.87 precision, and 0.72 for the F-measure. Both performances were calculated using a maximum distance of only 500 bp. As expected from previous analyses of the ENCODE region (Bajic et al. 2006a), we see that the performance of our program is better on this region than on the rest of the genome, which indicates that some of the FP in the genome setting are actually missed genes or missed TSSs.

To test this last claim, we retrieved datasets from the ENCODE project for Affymetrix Transcribed Fragments, Yale Transcriptionally Active Regions (TARs), and novel TARs from the DART system (Rozowsky et al. 2007). We combined these three sets with the CAGE data from Riken (single CAGE tags included). This set is called the Evidence for Transcriptional Activity set (EFTA). When comparing the predictions of EP3, we found that of all predictions made by EP3, 87% have a hit with EFTA within 125 bp, 95% have a hit within 500 bp, and 98% have a hit within 2000 bp within EFTA. When excluding the single CAGE tags, the rates drop to 80%, 92%, and 97%, respectively. These numbers indicate that EP3 has a very strong specificity and that many of the so-called FPs discussed above are in fact associated with transcriptionally active regions.

Furthermore, we compared the predictions of EP3 with two sets of DNase hypersensitivity sites (DHSS) retrieved from the ENCODE project. For the first set (encodeNhgriDnaseHsMpssCd4, seven cell types), 50% of the DHSS are near a prediction of EP3, and for the second set (encodeRegulomeDnaseGM06990Sites, one cell type), 28% of the DHSS are near an EP3 prediction. For both sets, the recall rate is lower than that on the CAGE set, but this was to be expected, as the DNase dataset covers only a limited number of cell types.

## Performance on different eukaryotic genomes

In the previous sections, we demonstrated the performance of EP3 on the human genome. We have also tested its performance on a wide range of other eukaryotes, including animals (Mus musculus, Tetraodon nigroviridis, Drosophila melanogaster), fungi (Saccharomyces cerevisiae, Schizosaccharomyces pombe), algae (Ostreococcus tauri, Ostreococcus pacifica), higher plants (Arabidopsis thaliana, Oryza sativa, Populus trichocarpa), and a protist (Plasmodium falciparum). Only data for human and mouse are available from the CAGE technique; therefore, we limited the analyses for the other eukaryotes to the data available from Ensembl. Table 3 shows the result when we used EP3 to predict promoter regions in other eukaryotes. The F-measure ranges from 0.17 to 0.71 on the different species, with *P. falciparum* (F = 0.17) and *D. melanogaster* (F = 0.19) on the low end of the scale, and *O. pacifica* (F = 0.71) and *O. tauri* (F = 0.66) giving the best results. We evaluated the performance of EP3 only; the other programs are not suitable for all other genomes because they are specifically trained for a single species. EP3 obtains a good score for some species, while for other species the score is worse. The F-score is a bit higher for mouse than for human, which indicates that the program performs well for mammals. Within the green lineage (green algae and land plants), there seem to be two groups, based on the genome size. The performance for the two algae (first group) is excellent (F > 0.65), which is probably partly due to the very small genome size and the still large gene space (~8000 genes). Due to the window approach to assess TPs and FPs, most predictions will be a TP because of the small genome. The second group comprised of rice, *Arabidopsis*, and poplar has larger genomes, and the performance of EP3 is comparable to that on
better (F = 0.45) than DragonGSF (F = 0.44) and PromoterInspector (F = 0.43). Although these differences seem small, it still concerns several hundreds of genes or TSSs. For instance, a 1% lower recall means 200 genes in the Ensembl dataset that are missed by that program. One percent less precision means 200 extra false predictions, so even the small differences are significant. In the case of the CAGE dataset, 1% difference corresponds to a difference of ~1250 TSSs. Combining this into the F-measure means that a program that has 1% lower value misses 1250 TSSs and predicts an additional 1250 wrong ones.

In conclusion, the best performing programs are EP3, DragonGSF, PromoterInspector, Eponine, and FirstEF, which have an F-measure > 0.4. Several programs (CpGProD, PromoterExplorer, N-Scan, and McPromoter) still perform quite well and have an F-measure > 0.25. The rest of the programs (PromFD, ARTS, DragonPF, PromoterScan, NNPP2.2, and Promoter2.0) do not perform very well on the complete human genome (F < 0.25). In all cases, these low F-measures are caused by a very low precision, which is obtained when the PPP outputs many FPs, a problem reminiscent of promoter prediction since the beginning.

The performance of PPPs on the Ensembl dataset is generally higher than that on the CAGE dataset, probably because of the different counting schemes for Ensembl data and the CAGE dataset, whereas the scheme for Ensembl data ignores all intergenic predictions. Some of the other programs may also perform better on the Ensembl data than EP3, because these programs have been trained on promoters of protein-coding genes and are therefore more gene-centric. For the top-performing programs, the balance between the recall (sensitivity) and precision (specificity) is mostly favoring the precision, with the exception of FirstEF, which is perfectly balanced on the CAGE dataset. The other programs often have high recall values, but at the cost of very low precision values.

## Performance on the ENCODE region

The ENCODE project aims to carefully annotate all functional elements in a small portion (1%) of the human genome. We used EP3 to make predictions on the 44 regions that cover ~30 Mb and compared the prediction with three different datasets. First, we compared our predictions to a set of known functional promoters (Cooper et al. 2006). The dataset is only partial for the ENCODE region because Cooper et al. tested only 642 putative promoters of the 921 they predicted. These partial data are insufficient to assess the precision of EP3. Of the 642 promoters tested, 387 were discovered to be functional. EP3 predicts 24% (recall) of the promoters that are marked as functional by Cooper et al. when using a maximum distance of 500 bp. This low recall rate is rather surprising, as the performance on the GENCODE and CAGE data is very good (see below). But when we look in detail, a much higher recall rate (40%) is obtained for the genes expressed in all 16 cell lines, which indicates that EP3 is biased toward broadly expressed genes. Of the 257 nonfunctional promoters, 18 (7%) are predicted by EP3.

Next, we compared the predictions of EP3 on the ENCODE region with the gene annotation from the GENCODE project (Harrow et al. 2006) and with the CAGE data from Riken. The GENCODE annotation was compared with the predictions using the classic method for calculating the performance, as it is a gene annotation similar to the one of Ensembl. We obtained a recall of 0.46, a precision of 0.72, and an F-measure of 0.56. The performance on the CAGE data was calculated with the novel method

presented here and has 0.61 recall, 0.87 precision, and 0.72 for the F-measure. Both performances were calculated using a maximum distance of only 500 bp. As expected from previous analyses of the ENCODE region (Bajic et al. 2006a), we see that the performance of our program is better on this region than on the rest of the genome, which indicates that some of the FP in the genome setting are actually missed genes or missed TSSs.

To test this last claim, we retrieved datasets from the ENCODE project for Affymetrix Transcribed Fragments, Yale Transcriptionally Active Regions (TARs), and novel TARs from the DART system (Rozowsky et al. 2007). We combined these three sets with the CAGE data from Riken (single CAGE tags included). This set is called the Evidence for Transcriptional Activity set (EFTA). When comparing the predictions of EP3, we found that of all predictions made by EP3, 87% have a hit with EFTA within 125 bp, 95% have a hit within 500 bp, and 98% have a hit within 2000 bp within EFTA. When excluding the single CAGE tags, the rates drop to 80%, 92%, and 97%, respectively. These numbers indicate that EP3 has a very strong specificity and that many of the so-called FPs discussed above are in fact associated with transcriptionally active regions.

Furthermore, we compared the predictions of EP3 with two sets of DNase hypersensitivity sites (DHSS) retrieved from the ENCODE project. For the first set (encodeNhgriDnaseHsMpssCd4, seven cell types), 50% of the DHSS are near a prediction of EP3, and for the second set (encodeRegulomeDnaseGM06990Sites, one cell type), 28% of the DHSS are near an EP3 prediction. For both sets, the recall rate is lower than that on the CAGE set, but this was to be expected, as the DNase dataset covers only a limited number of cell types.

## Performance on different eukaryotic genomes

In the previous sections, we demonstrated the performance of EP3 on the human genome. We have also tested its performance on a wide range of other eukaryotes, including animals (Mus musculus, Tetraodon nigroviridis, Drosophila melanogaster), fungi (Saccharomyces cerevisiae, Schizosaccharomyces pombe), algae (Ostreococcus tauri, Ostreococcus pacifica), higher plants (Arabidopsis thaliana, Oryza sativa, Populus trichocarpa), and a protist (Plasmodium falciparum). Only data for human and mouse are available from the CAGE technique; therefore, we limited the analyses for the other eukaryotes to the data available from Ensembl. Table 3 shows the result when we used EP3 to predict promoter regions in other eukaryotes. The F-measure ranges from 0.17 to 0.71 on the different species, with *P. falciparum* (F = 0.17) and *D. melanogaster* (F = 0.19) on the low end of the scale, and *O. pacifica* (F = 0.71) and *O. tauri* (F = 0.66) giving the best results. We evaluated the performance of EP3 only; the other programs are not suitable for all other genomes because they are specifically trained for a single species. EP3 obtains a good score for some species, while for other species the score is worse. The F-score is a bit higher for mouse than for human, which indicates that the program performs well for mammals. Within the green lineage (green algae and land plants), there seem to be two groups, based on the genome size. The performance for the two algae (first group) is excellent (F > 0.65), which is probably partly due to the very small genome size and the still large gene space (~8000 genes). Due to the window approach to assess TPs and FPs, most predictions will be a TP because of the small genome. The second group comprised of rice, *Arabidopsis*, and poplar has larger genomes, and the performance of EP3 is comparable to that on
<table><caption>**Table 3.** Performance of EP3 on different eukaryotic genomes with a maximum allowed mismatch distance of 500 bp</caption><thead><tr><th>Species</th><th>F-measure</th><th>Size (Mb)<sup>a</sup></th></tr></thead><tbody><tr><td><i>P. falciparum</i></td><td>0.17</td><td>23</td></tr><tr><td><i>O. pacifica</i></td><td>0.71</td><td>13</td></tr><tr><td><i>O. tauri</i></td><td>0.66</td><td>13</td></tr><tr><td><i>A. thaliana</i></td><td>0.37</td><td>120</td></tr><tr><td><i>O. sativa</i></td><td>0.53</td><td>370</td></tr><tr><td><i>P. trichocarpa</i></td><td>0.46</td><td>300</td></tr><tr><td><i>S. cerevisiae</i></td><td>0.42</td><td>12</td></tr><tr><td><i>S. pombe</i></td><td>0.31</td><td>12</td></tr><tr><td><i>C. elegans</i></td><td>0.26</td><td>100</td></tr><tr><td><i>D. melanogaster</i></td><td>0.19</td><td>130</td></tr><tr><td><i>T. nigroviridis</i></td><td>0.23</td><td>220</td></tr><tr><td><i>M. musculus</i></td><td>0.46</td><td>2500</td></tr><tr><td><i>H. sapiens</i></td><td>0.44</td><td>3000</td></tr></tbody></table>

Note that the datasets that have been retrieved from Ensembl are not all of the same quality, even though the organisms presented were selected for having a good annotation.

mammals. The performance for the two yeasts is lower than that on the two *Ostreococcus* genomes, which have roughly the same genome size. This is probably due to the less obvious profile in yeasts compared with the one in algae (see Fig. 2). The performance of EP3 on *Drosophila* and *Plasmodium* is weak. In the case of the fruit fly, this is likely due to its very different structural profile, as observed in Figure 2, while for *Plasmodium* the low performance is probably caused by the rather flat profile observed in protists that makes it very difficult for EP3 to distinguish promoter regions from other parts of the genome. Therefore, for some species, for example *D. melanogaster*, specifically trained programs might perform better (Ohler 2006).

## Recognizing different promoter types

Besides genes that code for proteins, there are also genes that are transcribed but for which the RNA is not translated into proteins, so-called noncoding RNAs. These genes produce transcripts that function directly as structural, catalytic, or regulatory RNAs. Recent screens for such genes revealed a surprisingly large number of them, with prominent roles such as guiding the post-transcriptional regulation of protein-coding genes (Eddy 2001; Bartel 2004). Previous studies on promoter prediction focused mainly on a single type of promoter, most often the promoter of protein-coding genes, which are transcribed by RNAP II. Using the Ensembl annotation for humans, we show that our approach is also suited to predict other types of promoters. Although the program works best to identify and delineate promoters of protein-coding genes, it can also be used to detect promoters of snRNA, rRNA, miRNAs, snoRNA, and tRNA genes. Other types of noncoding genes such as scRNA and mitochondrial rRNA were not considered because of the lack of data in the Ensembl database. Table 4 shows the different recall (sensitivity) rates for the different types of genes in this study (snRNA, rRNA, miRNA, snoRNA, and tRNA). The precision cannot be calculated for these analyses because the program will always predict all types of promoters and it will never give predictions specific for a single type of promoter. Therefore, we focus only on how many known non-coding RNA promoters we can identify with our approach. Although EP3 can identify non-protein-coding genes, other programs have higher recall rates. However, the programs with high

recall rates are also the ones that performed worse when we applied them to humans. From Table 1, we see that the programs that have high recall and low precision have the lowest F-measure, those also being the programs that have the highest recall for the non-protein-coding genes. Although EP3 is thus also capable of predicting the promoters of noncoding genes, its performance is significantly lower, most probably because the peak in the profile is much smaller (see Fig. 2). Nevertheless, compared with the other top-performing programs on the whole genome (DragonGSF and PromoterInspector), EP3 has very similar recall rates. From this analysis, it is clear that general-purpose PPPs such as those listed in Table 1 are best suited for the prediction of protein-coding gene promoters. For miRNA and tRNA promoters, there are probably better approaches using specifically trained tools for identification of these types of promoters (Lowe and Eddy 1997; Zhou et al. 2007). Finally, EP3 is slightly biased toward GC-rich promoters because the structural feature is most outspoken in these promoters. This bias toward GC-rich promoters is also present in all top-performing PPPs from Table 1 (Scherf et al. 2001; Bajic et al. 2004) and indicates that CpG-island-associated housekeeping genes are favored in the predictions.

## Conclusion

The evaluation of PPPs in a whole-genome context is crucial to understand the true performance of the program. Evaluation on a small test set, such as the Eukaryotic Promoter Database (Schmid et al. 2006), does not provide sufficient insight into the real performance of the program when used in actual genome annotation projects. The recall and precision values we found in our analysis are lower than those reported in the original papers, where in most cases the evaluation was done on a (much) smaller dataset. Even the evaluation of a complete chromosome is not sufficient, as there are huge differences in nucleotide content and gene density. If possible, it is advisable to use transcription data, such as the CAGE data, to assess the performance. Transcription data are superior to the gene annotation and its associated way of counting TPs; the gene annotation may not give a complete picture on the performance because it completely ignores intergenic predictions. In short, one should assess a promoter predictor on the whole genome, preferably validating with TSS data.

<table><caption>**Table 4.** Recall (sensitivity) for the different gene types for the different programs</caption><thead><tr><th rowspan="2">Program</th><th colspan="5">Recall (%, 2000 bp maximum distance)</th></tr><tr><th>mRNA</th><th>miRNA</th><th>snRNA</th><th>snoRNA</th><th>rRNA</th></tr></thead><tbody><tr><td>ARTS</td><td>91</td><td>62</td><td>44</td><td>47</td><td>59</td></tr><tr><td>CpgProD</td><td>69</td><td>27</td><td>14</td><td>17</td><td>14</td></tr><tr><td>DragonGSF</td><td>59</td><td>15</td><td>5</td><td>14</td><td>4</td></tr><tr><td>DragonPF</td><td>82</td><td>55</td><td>37</td><td>43</td><td>37</td></tr><tr><td>Eponine</td><td>44</td><td>11</td><td>2</td><td>7</td><td>8</td></tr><tr><td>FirstEF</td><td>74</td><td>32</td><td>10</td><td>20</td><td>16</td></tr><tr><td>McPromoter (0.0)</td><td>35</td><td>10</td><td>4</td><td>10</td><td>7</td></tr><tr><td>McPromoter (-0.05)</td><td>87</td><td>80</td><td>51</td><td>65</td><td>64</td></tr><tr><td>NNPP2.2 (0.99)</td><td>10</td><td>9</td><td>9</td><td>7</td><td>7</td></tr><tr><td>PromoterExplorer</td><td>77</td><td>40</td><td>17</td><td>26</td><td>20</td></tr><tr><td>Promoter2.0 (high)</td><td>61</td><td>52</td><td>70</td><td>62</td><td>65</td></tr><tr><td>Promoter2.0 (medium)</td><td>99</td><td>96</td><td>95</td><td>94</td><td>94</td></tr><tr><td>EP3</td><td>53</td><td>19</td><td>2</td><td>7</td><td>12</td></tr></tbody></table>

For each program, the sensitivity percentages are shown when using a maximum mismatch of 2000 bp.

<footnote><sup>a</sup>Approximate genome size in megabases (Mb).</footnote>
<table><caption>**Table 3.** Performance of EP3 on different eukaryotic genomes with a maximum allowed mismatch distance of 500 bp</caption><thead><tr><th>Species</th><th>F-measure</th><th>Size (Mb)<sup>a</sup></th></tr></thead><tbody><tr><td><i>P. falciparum</i></td><td>0.17</td><td>23</td></tr><tr><td><i>O. pacifica</i></td><td>0.71</td><td>13</td></tr><tr><td><i>O. tauri</i></td><td>0.66</td><td>13</td></tr><tr><td><i>A. thaliana</i></td><td>0.37</td><td>120</td></tr><tr><td><i>O. sativa</i></td><td>0.53</td><td>370</td></tr><tr><td><i>P. trichocarpa</i></td><td>0.46</td><td>300</td></tr><tr><td><i>S. cerevisiae</i></td><td>0.42</td><td>12</td></tr><tr><td><i>S. pombe</i></td><td>0.31</td><td>12</td></tr><tr><td><i>C. elegans</i></td><td>0.26</td><td>100</td></tr><tr><td><i>D. melanogaster</i></td><td>0.19</td><td>130</td></tr><tr><td><i>T. nigroviridis</i></td><td>0.23</td><td>220</td></tr><tr><td><i>M. musculus</i></td><td>0.46</td><td>2500</td></tr><tr><td><i>H. sapiens</i></td><td>0.44</td><td>3000</td></tr></tbody></table>

Note that the datasets that have been retrieved from Ensembl are not all of the same quality, even though the organisms presented were selected for having a good annotation.

mammals. The performance for the two yeasts is lower than that on the two *Ostreococcus* genomes, which have roughly the same genome size. This is probably due to the less obvious profile in yeasts compared with the one in algae (see Fig. 2). The performance of EP3 on *Drosophila* and *Plasmodium* is weak. In the case of the fruit fly, this is likely due to its very different structural profile, as observed in Figure 2, while for *Plasmodium* the low performance is probably caused by the rather flat profile observed in protists that makes it very difficult for EP3 to distinguish promoter regions from other parts of the genome. Therefore, for some species, for example *D. melanogaster*, specifically trained programs might perform better (Ohler 2006).

## Recognizing different promoter types

Besides genes that code for proteins, there are also genes that are transcribed but for which the RNA is not translated into proteins, so-called noncoding RNAs. These genes produce transcripts that function directly as structural, catalytic, or regulatory RNAs. Recent screens for such genes revealed a surprisingly large number of them, with prominent roles such as guiding the post-transcriptional regulation of protein-coding genes (Eddy 2001; Bartel 2004). Previous studies on promoter prediction focused mainly on a single type of promoter, most often the promoter of protein-coding genes, which are transcribed by RNAP II. Using the Ensembl annotation for humans, we show that our approach is also suited to predict other types of promoters. Although the program works best to identify and delineate promoters of protein-coding genes, it can also be used to detect promoters of snRNA, rRNA, miRNAs, snoRNA, and tRNA genes. Other types of noncoding genes such as scRNA and mitochondrial rRNA were not considered because of the lack of data in the Ensembl database. Table 4 shows the different recall (sensitivity) rates for the different types of genes in this study (snRNA, rRNA, miRNA, snoRNA, and tRNA). The precision cannot be calculated for these analyses because the program will always predict all types of promoters and it will never give predictions specific for a single type of promoter. Therefore, we focus only on how many known non-coding RNA promoters we can identify with our approach. Although EP3 can identify non-protein-coding genes, other programs have higher recall rates. However, the programs with high

recall rates are also the ones that performed worse when we applied them to humans. From Table 1, we see that the programs that have high recall and low precision have the lowest F-measure, those also being the programs that have the highest recall for the non-protein-coding genes. Although EP3 is thus also capable of predicting the promoters of noncoding genes, its performance is significantly lower, most probably because the peak in the profile is much smaller (see Fig. 2). Nevertheless, compared with the other top-performing programs on the whole genome (DragonGSF and PromoterInspector), EP3 has very similar recall rates. From this analysis, it is clear that general-purpose PPPs such as those listed in Table 1 are best suited for the prediction of protein-coding gene promoters. For miRNA and tRNA promoters, there are probably better approaches using specifically trained tools for identification of these types of promoters (Lowe and Eddy 1997; Zhou et al. 2007). Finally, EP3 is slightly biased toward GC-rich promoters because the structural feature is most outspoken in these promoters. This bias toward GC-rich promoters is also present in all top-performing PPPs from Table 1 (Scherf et al. 2001; Bajic et al. 2004) and indicates that CpG-island-associated housekeeping genes are favored in the predictions.

## Conclusion

The evaluation of PPPs in a whole-genome context is crucial to understand the true performance of the program. Evaluation on a small test set, such as the Eukaryotic Promoter Database (Schmid et al. 2006), does not provide sufficient insight into the real performance of the program when used in actual genome annotation projects. The recall and precision values we found in our analysis are lower than those reported in the original papers, where in most cases the evaluation was done on a (much) smaller dataset. Even the evaluation of a complete chromosome is not sufficient, as there are huge differences in nucleotide content and gene density. If possible, it is advisable to use transcription data, such as the CAGE data, to assess the performance. Transcription data are superior to the gene annotation and its associated way of counting TPs; the gene annotation may not give a complete picture on the performance because it completely ignores intergenic predictions. In short, one should assess a promoter predictor on the whole genome, preferably validating with TSS data.

<table><caption>**Table 4.** Recall (sensitivity) for the different gene types for the different programs</caption><thead><tr><th rowspan="2">Program</th><th colspan="5">Recall (%, 2000 bp maximum distance)</th></tr><tr><th>mRNA</th><th>miRNA</th><th>snRNA</th><th>snoRNA</th><th>rRNA</th></tr></thead><tbody><tr><td>ARTS</td><td>91</td><td>62</td><td>44</td><td>47</td><td>59</td></tr><tr><td>CpgProD</td><td>69</td><td>27</td><td>14</td><td>17</td><td>14</td></tr><tr><td>DragonGSF</td><td>59</td><td>15</td><td>5</td><td>14</td><td>4</td></tr><tr><td>DragonPF</td><td>82</td><td>55</td><td>37</td><td>43</td><td>37</td></tr><tr><td>Eponine</td><td>44</td><td>11</td><td>2</td><td>7</td><td>8</td></tr><tr><td>FirstEF</td><td>74</td><td>32</td><td>10</td><td>20</td><td>16</td></tr><tr><td>McPromoter (0.0)</td><td>35</td><td>10</td><td>4</td><td>10</td><td>7</td></tr><tr><td>McPromoter (-0.05)</td><td>87</td><td>80</td><td>51</td><td>65</td><td>64</td></tr><tr><td>NNPP2.2 (0.99)</td><td>10</td><td>9</td><td>9</td><td>7</td><td>7</td></tr><tr><td>PromoterExplorer</td><td>77</td><td>40</td><td>17</td><td>26</td><td>20</td></tr><tr><td>Promoter2.0 (high)</td><td>61</td><td>52</td><td>70</td><td>62</td><td>65</td></tr><tr><td>Promoter2.0 (medium)</td><td>99</td><td>96</td><td>95</td><td>94</td><td>94</td></tr><tr><td>EP3</td><td>53</td><td>19</td><td>2</td><td>7</td><td>12</td></tr></tbody></table>

For each program, the sensitivity percentages are shown when using a maximum mismatch of 2000 bp.

<footnote><sup>a</sup>Approximate genome size in megabases (Mb).</footnote>
While EP3 does not outperform its peers by much, the program has several additional advantages compared with other PPPs. EP3 requires no training or parameter tuning, unlike other programs that need extensive amounts of experimentally determined data for the training of their model (Ohler et al. 2000; Scherf et al. 2000; Davuluri et al. 2001; Down and Hubbard 2002; Bajic et al. 2003). When working on a genomic scale, speed and memory requirements also are of importance. EP3 is very fast (for instance, it takes <1 h to annotate the complete human genome), requires little memory, and can thus be run on a home computer; in contrast, some programs require a computer cluster of 80 machines for nearly a week to process the human genome (data not shown). Besides performing very well, especially in light of its simplicity, EP3 can handle many eukaryotic genomes without modifications, as was shown here.

Up to now, most PPPs have used biologically driven small-scale features to build a model for promoter recognition (Fickett and Hatzigeorgiou 1997; Hannenhalli and Levy 2001; Werner 2003). These features are different for each species, and thus the current programs need to be retrained for each species. Here, we have presented a simple approach that employs more global structural features of the DNA sequence in promoter and non-promoter regions. Our technique does not use any complex machine learning algorithms, for which it is often impossible to infer any new knowledge from the model itself. On the contrary, our method requires no training whatsoever and can be applied to any eukaryotic genome. Numerical profiles representing sequence-dependent properties investigated on a large scale and a very large number of genes show a remarkable feature of promoters. To the best of our knowledge, this is the first time that this feature is described for multiple eukaryotes. The feature associated with promoters is an extended region where the DNA is more stable. This large-scale feature is not present in other types of genomic sequence. While it is most outspoken in vertebrate promoters, it is also present in other eukaryotic promoters, so it seems to be universal and largely independent of the presence or absence of binding sites. As a result, EP3 is capable of performing very well on several eukaryotic species, without the need for any training, which is a unique achievement among promoter prediction software, allowing true "ab initio" promoter prediction. Finally, analysis of the well documented ENCODE region shows that most predictions are indeed associated with transcriptionally active regions.

## Methods

### Datasets

The CAGE datasets for human have been retrieved from the FANTOM3 project (http://fantom.gsc.riken.go.jp/). The database was compiled by Carninci et al. (2006) and was obtained through the CAGE technique, which identifies all possible TSSs with very high accuracy (Shiraki et al. 2003). This dataset covers the whole human genome and can thus be used for validation purposes. Only tag clusters with at least two mapped tags on the same genomic location were considered to be real TSSs, which should filter out most FPs from the CAGE technique. This filtering resulted in 123,400 unique start sites for human. The whole genome sequences for human (hg17) and mouse (mm5) were retrieved from the UCSC Genome Bioinformatics Site (http://genome.ucsc.edu/) (Kuhn et al. 2007).

We retrieved the following genes in human using the BioMart tool at the Ensembl website, release 37 (http://www.

ensembl.org): (1) 20,297 protein-coding, (2) 1382 small nuclear RNA (snRNA), (3) 330 ribosomal RNA (rRNA), (4) 326 microRNA (miRNA), and (5) 642 small nucleolar RNA (snoRNA) (Hubbard et al. 2007). The dataset for transfer RNA (tRNA) was downloaded from the tRNAscan-SE homepage (http://lowelab.ucsc.edu/GtRNAbd/Hsapi/) (Lowe and Eddy 1997). The datasets with annotations for different species (Table 3) were retrieved from two different sources: The protein-coding genes for human, mouse, Tetraodon, fruit fly, and yeast were retrieved using BioMart, while the genomic sequences were retrieved from the Ensembl ftp server (ftp://ftp.ensembl.org/pub/). The sequences and annotation for the five plant species were retrieved from in-house data (http://bioinformatics.psb.ugent.be). The genome and annotation for *P. falciparum* were retrieved from the Sanger database (http://www.sanger.ac.uk/).

The data for the ENCODE analysis have been downloaded from the UCSC ENCODE repository (http://genome.ucsc.edu/ENCODE/). The data for the functional promoters in the encode region were retrieved from Supplemental Table D (Cooper et al. 2006). The data for the novel transcriptionally active regions were retrieved from the DART database (http:// dart.gersteinlab.org/ENCODE/GT/).

### Calculating structural profiles

Profiles from DNA sequences are calculated as follows: First, the nucleotide sequence is converted into a sequence of numbers (i.e., a numerical profile). This is done by replacing each di-nucleotide (or trinucleotide, depending on the physico-chemical feature used) with its corresponding structural value, which is obtained from experimentally validated conversion tables. Florquin et al. (2005) provide references to the protocols to obtain these conversion tables. These contain values, for example, for stabilizing energy of Z-DNA (Ho et al. 1990), DNA denaturation (Blake and Delcourt 1998; Blake et al. 1999), protein-induced deformability (Olson et al. 1998), and duplex-free energy (Sugimoto et al. 1996). Next, we average over several values, and the number of values we use for computing this average is called the window size, in our case 400 bp. The value 400 was chosen because it gave the best performance in a whole-genome context (see Results). If a value in the profile exceeds a certain threshold, this points to a putative promoter region (see below). For properties that show a cleft instead of a peak (see below), we took the inverse value to make sure we also obtained peaks. This was done to be certain we could use the same approach for all structural features.

### Prediction algorithm

The algorithm we use (EP3) to identify promoters has two internal parameters, one for the length of the window and one for the deviation from the average. Both have been determined empirically. The optimal length of the window is 400 bp. We determined this by trying values between 50 and 2000 and seeing which one performed best. Smaller windows result in more predictions, but the number of false predictions increases faster than the rate of true predictions. Larger values result in fewer predictions but higher precision. The *F*-measure indicates that 400 bp is optimal for all organisms. This value appears to be independent of the size of the genome because the organisms we analyzed ranged in size from 12 Mbp to 3 Gbp, and 400 bp was always the optimal value. Because the overall profile of a genome can vary significantly between different species due to different GC content, the threshold should be determined dynamically for each genome. As a basis, we use the average and standard deviation of the values of a property for the whole genome. For bigger ge-
While EP3 does not outperform its peers by much, the program has several additional advantages compared with other PPPs. EP3 requires no training or parameter tuning, unlike other programs that need extensive amounts of experimentally determined data for the training of their model (Ohler et al. 2000; Scherf et al. 2000; Davuluri et al. 2001; Down and Hubbard 2002; Bajic et al. 2003). When working on a genomic scale, speed and memory requirements also are of importance. EP3 is very fast (for instance, it takes <1 h to annotate the complete human genome), requires little memory, and can thus be run on a home computer; in contrast, some programs require a computer cluster of 80 machines for nearly a week to process the human genome (data not shown). Besides performing very well, especially in light of its simplicity, EP3 can handle many eukaryotic genomes without modifications, as was shown here.

Up to now, most PPPs have used biologically driven small-scale features to build a model for promoter recognition (Fickett and Hatzigeorgiou 1997; Hannenhalli and Levy 2001; Werner 2003). These features are different for each species, and thus the current programs need to be retrained for each species. Here, we have presented a simple approach that employs more global structural features of the DNA sequence in promoter and non-promoter regions. Our technique does not use any complex machine learning algorithms, for which it is often impossible to infer any new knowledge from the model itself. On the contrary, our method requires no training whatsoever and can be applied to any eukaryotic genome. Numerical profiles representing sequence-dependent properties investigated on a large scale and a very large number of genes show a remarkable feature of promoters. To the best of our knowledge, this is the first time that this feature is described for multiple eukaryotes. The feature associated with promoters is an extended region where the DNA is more stable. This large-scale feature is not present in other types of genomic sequence. While it is most outspoken in vertebrate promoters, it is also present in other eukaryotic promoters, so it seems to be universal and largely independent of the presence or absence of binding sites. As a result, EP3 is capable of performing very well on several eukaryotic species, without the need for any training, which is a unique achievement among promoter prediction software, allowing true "ab initio" promoter prediction. Finally, analysis of the well documented ENCODE region shows that most predictions are indeed associated with transcriptionally active regions.

## Methods

### Datasets

The CAGE datasets for human have been retrieved from the FANTOM3 project (http://fantom.gsc.riken.go.jp/). The database was compiled by Carninci et al. (2006) and was obtained through the CAGE technique, which identifies all possible TSSs with very high accuracy (Shiraki et al. 2003). This dataset covers the whole human genome and can thus be used for validation purposes. Only tag clusters with at least two mapped tags on the same genomic location were considered to be real TSSs, which should filter out most FPs from the CAGE technique. This filtering resulted in 123,400 unique start sites for human. The whole genome sequences for human (hg17) and mouse (mm5) were retrieved from the UCSC Genome Bioinformatics Site (http://genome.ucsc.edu/) (Kuhn et al. 2007).

We retrieved the following genes in human using the BioMart tool at the Ensembl website, release 37 (http://www.

ensembl.org): (1) 20,297 protein-coding, (2) 1382 small nuclear RNA (snRNA), (3) 330 ribosomal RNA (rRNA), (4) 326 microRNA (miRNA), and (5) 642 small nucleolar RNA (snoRNA) (Hubbard et al. 2007). The dataset for transfer RNA (tRNA) was downloaded from the tRNAscan-SE homepage (http://lowelab.ucsc.edu/GtRNAbd/Hsapi/) (Lowe and Eddy 1997). The datasets with annotations for different species (Table 3) were retrieved from two different sources: The protein-coding genes for human, mouse, Tetraodon, fruit fly, and yeast were retrieved using BioMart, while the genomic sequences were retrieved from the Ensembl ftp server (ftp://ftp.ensembl.org/pub/). The sequences and annotation for the five plant species were retrieved from in-house data (http://bioinformatics.psb.ugent.be). The genome and annotation for *P. falciparum* were retrieved from the Sanger database (http://www.sanger.ac.uk/).

The data for the ENCODE analysis have been downloaded from the UCSC ENCODE repository (http://genome.ucsc.edu/ENCODE/). The data for the functional promoters in the encode region were retrieved from Supplemental Table D (Cooper et al. 2006). The data for the novel transcriptionally active regions were retrieved from the DART database (http://dart.gersteinlab.org/ENCODE/GT/).

### Calculating structural profiles

Profiles from DNA sequences are calculated as follows: First, the nucleotide sequence is converted into a sequence of numbers (i.e., a numerical profile). This is done by replacing each di-nucleotide (or trinucleotide, depending on the physico-chemical feature used) with its corresponding structural value, which is obtained from experimentally validated conversion tables. Florquin et al. (2005) provide references to the protocols to obtain these conversion tables. These contain values, for example, for stabilizing energy of Z-DNA (Ho et al. 1990), DNA denaturation (Blake and Delcourt 1998; Blake et al. 1999), protein-induced deformability (Olson et al. 1998), and duplex-free energy (Sugimoto et al. 1996). Next, we average over several values, and the number of values we use for computing this average is called the window size, in our case 400 bp. The value 400 was chosen because it gave the best performance in a whole-genome context (see Results). If a value in the profile exceeds a certain threshold, this points to a putative promoter region (see below). For properties that show a cleft instead of a peak (see below), we took the inverse value to make sure we also obtained peaks. This was done to be certain we could use the same approach for all structural features.

### Prediction algorithm

The algorithm we use (EP3) to identify promoters has two internal parameters, one for the length of the window and one for the deviation from the average. Both have been determined empirically. The optimal length of the window is 400 bp. We determined this by trying values between 50 and 2000 and seeing which one performed best. Smaller windows result in more predictions, but the number of false predictions increases faster than the rate of true predictions. Larger values result in fewer predictions but higher precision. The *F*-measure indicates that 400 bp is optimal for all organisms. This value appears to be independent of the size of the genome because the organisms we analyzed ranged in size from 12 Mbp to 3 Gbp, and 400 bp was always the optimal value. Because the overall profile of a genome can vary significantly between different species due to different GC content, the threshold should be determined dynamically for each genome. As a basis, we use the average and standard deviation of the values of a property for the whole genome. For bigger ge-
nomes (human and mouse), this deviation should be larger than for other species. Empirically, we determined that the optimal threshold is the average plus one standard deviation for all species, except for human and mouse where this is the average plus three standard deviations. Future research with more available genomes, especially ones that fill in the gap between the small genomes (≤300 Mbp) and the large ones (≥2.5 Gbp) may allow us to define a function that can calculate the optimal multiplier for the standard deviation from the genome size. When the window size (400 bp) and the threshold (depends on genome size) are known, the prediction algorithm is straightforward. The program calculates the profile, and each time a value in the profile is above the threshold for this genome, the location is predicted to be a putative promoter region.

## Evaluating predictions

A measure for the performance of a PPP is the harmonic mean of the recall (sensitivity) and the precision (specificity), known as the *F*-measure (van Rijsbergen 1979). The higher this value, the better the program is able to correctly predict promoters. The recall or sensitivity is the number of predicted promoters (TP) divided by the total number of promoters (TP+FN). The precision or specificity is the number of correct predictions (TP) divided by the total number of predictions (TP+FP). To assess the number of TPs, FPs, and FNs, we need to define the maximum allowed mismatch between the prediction and the true TSS. In previous studies, this value was set to 2000 bp (Bajic et al. 2004; Sonnenburg et al. 2006; Xie et al. 2006), although more recently, for the ENCODE project and for two PPPs, a maximum mismatch of 1000 bp also was evaluated (Bajic et al. 2006a). In Tables 1–3, we display the performance of the programs for three different maximum mismatch values: 2000 bp, 1000 bp, and 500 bp. The smaller the mismatch allowed, the stricter the evaluation. We ranked the performance of the different programs according to the smallest window size (500 bp). In the classic way to evaluate a PPP (Bajic et al. 2004), one will start from gene annotations, e.g., as compiled at Ensembl. A TP is then defined as a TSS that has a prediction within the maximum allowed distance from the true TSS as annotated in the database. A FP is a prediction that lies inside the gene but not within the maximum allowed distance from the TSS. A FN is a true TSS from the database that has no prediction within the maximum allowed mismatch. All predictions that fall within an intergenic region and further from a TSS than the maximum allowed mismatch are ignored. Counting a prediction inside a gene as a FP will be incorrect in some cases, as Carninci et al. (2006) showed that transcription can start in the 5' untranslated region (UTR) of a gene, as well as in an exon or in the 3' UTR. Therefore, not taking into account any of the intergenic predictions or only considering the TSSs located at the 5' UTR gives a biased view of a promoter predictor. Using a gene annotation and ignoring the intergenic predictions will not give a real estimate of the performance of the prediction program for two reasons: (1) predictions inside a gene may in fact be true TSSs; TPs might be underestimated, and consequently FPs are then overestimated; (2) intergenic predictions outside the TSS region are ignored because it is impossible to know whether they are real FPs or whether they point to genes missed by the annotation process. This may in turn lead to underestimation or overestimation of FPs. When using gene annotation with all its uncertainties, it is impossible to know the actual performance of the program in a genomic setting, especially in light of the evidence that much more of the genome is transcribed than just the protein-coding regions (Kapranov et al. 2007; Peters et al. 2007; Pongavic et al. 2007).

To address some of these caveats, we developed an alternative way to assess the performance of a promoter prediction program that does not exhibit these shortcomings. Recently available whole-genome TSS data allow for a new way to count and evaluate. We used a dataset of TSSs for the whole human genome that were characterized experimentally using high-throughput cap-analysis of gene expression (CAGE) (for review, see Carninci 2006). In contrast to previous studies (Bajic et al. 2004, 2006a; Sonnenburg et al. 2006), this dataset allows counting the number of TPs and FPs in an objective way. The dataset contains locations where transcription starts. A TP is a prediction that is within the maximum allowed mismatch from a true TSS, a FN is a true TSS that has no prediction, and a FP is a prediction that is not near a true TSS from the dataset.

## Availability

A user-friendly online implementation of our promoter prediction software EP3 is available at http://bioinformatics.psb.ugent.be/. The program allows the user to process FASTA-formatted sequence files. The online implementation is also available as a standalone application.

## Acknowledgments

We thank Pierre Hilson for fruitful discussions and three anonymous referees for their constructive comments that greatly helped improve the original version of the manuscript. This work was supported by a grant from the Research Foundation–Flanders (3G031805). T.A. thanks the Institute for the Promotion of Innovation by Science and Technology in Flanders for a pre-doctoral fellowship. Y.S. is a Postdoctoral Fellow of the Research Foundation–Flanders.

## References

* Aerts, S., Thijs, G., Dabrowski, M., Moreau, Y., and De Moor, B. 2004. Comprehensive analysis of the base composition around the transcription start site in Metazoa. *BMC Genomics* **5**: 34. doi: 10.1186/1471-2164-5-34.

* Amaldi, F., Bozzoni, I., Beccari, E., and Pierandrei-Amaldi, P. 1989. Expression of ribosomal protein genes and regulation of ribosome biosynthesis in *Xenopus* development. *Trends Biochem. Sci.* **14**: 175-178.

* Baek, D., Davis, C., Ewing, B., Gordon, D., and Green, P. 2007. Characterization and predictive discovery of evolutionarily conserved mammalian alternative promoters. *Genome Res.* **17**: 145-155.

* Bajic, V.B., Seah, S.H., Chong, A., Zhang, G., Koh, J.L.Y., and Brusic, V. 2002. Dragon Promoter Finder: Recognition of vertebrate RNA polymerase II promoters. *Bioinformatics* **18**: 198-199.

* Bajic, V.B., Seah, S.H., Chong, A., Krishnan, S.P.T., Koh, J.L.Y., and Brusic, V. 2003. Computer model for recognition of functional transcription start sites in RNA polymerase II promoters of vertebrates. *J. Mol. Graph. Model.* **21**: 323-332.

* Bajic, V.B., Tan, S.L., Suzuki, Y., and Sugano, S. 2004. Promoter prediction analysis on the whole human genome. *Nat. Biotechnol.* **22**: 1467-1473.

* Bajic, V.B., Brent, M.R., Brown, R.H., Frankish, A., Harrow, J., Ohler, U., Solovyev, V.V., and Tan, S.L. 2006a. Performance assessment of promoter predictions on ENCODE regions in the EGASP experiment. *Genome Biol.* **7**: S3.1-3.13. doi: 10.1186/gb-2006-7-s1-s3.

* Bajic, V.B., Tan, S.L., Christoffels, A., Schonbach, C., Lipovich, L., Yang, L., Hofmann, O., Kruger, A., Hide, W., Kai, C., et al. 2006b. Mice and men: Their promoter properties. *PLoS Genet.* **2**: e54. doi: 10.1371/journal.pgen.00020054.

* Baldi, P., Chauvin, Y., Brunak, S., Gorodkin, J., and Pedersen, A.G. 1998. Computational applications of DNA structural scales. *Proc. Int. Conf. Intell. Syst. Mol. Biol.* **6**: 35-42.

* Bartel, D.P. 2004. MicroRNAs: Genomics, biogenesis, mechanism, and function. *Cell* **116**: 281-297.

* Bird, A. 2002. DNA methylation patterns and epigenetic memory. *Genes & Dev.* **16**: 6-21.
nomes (human and mouse), this deviation should be larger than for other species. Empirically, we determined that the optimal threshold is the average plus one standard deviation for all species, except for human and mouse where this is the average plus three standard deviations. Future research with more available genomes, especially ones that fill in the gap between the small genomes (≤300 Mbp) and the large ones (≥2.5 Gbp) may allow us to define a function that can calculate the optimal multiplier for the standard deviation from the genome size. When the window size (400 bp) and the threshold (depends on genome size) are known, the prediction algorithm is straightforward. The program calculates the profile, and each time a value in the profile is above the threshold for this genome, the location is predicted to be a putative promoter region.

## Evaluating predictions

A measure for the performance of a PPP is the harmonic mean of the recall (sensitivity) and the precision (specificity), known as the *F*-measure (van Rijsbergen 1979). The higher this value, the better the program is able to correctly predict promoters. The recall or sensitivity is the number of predicted promoters (TP) divided by the total number of promoters (TP+FN). The precision or specificity is the number of correct predictions (TP) divided by the total number of predictions (TP+FP). To assess the number of TPs, FPs, and FNs, we need to define the maximum allowed mismatch between the prediction and the true TSS. In previous studies, this value was set to 2000 bp (Bajic et al. 2004; Sonnenburg et al. 2006; Xie et al. 2006), although more recently, for the ENCODE project and for two PPPs, a maximum mismatch of 1000 bp also was evaluated (Bajic et al. 2006a). In Tables 1–3, we display the performance of the programs for three different maximum mismatch values: 2000 bp, 1000 bp, and 500 bp. The smaller the mismatch allowed, the stricter the evaluation. We ranked the performance of the different programs according to the smallest window size (500 bp). In the classic way to evaluate a PPP (Bajic et al. 2004), one will start from gene annotations, e.g., as compiled at Ensembl. A TP is then defined as a TSS that has a prediction within the maximum allowed distance from the true TSS as annotated in the database. A FP is a prediction that lies inside the gene but not within the maximum allowed distance from the TSS. A FN is a true TSS from the database that has no prediction within the maximum allowed mismatch. All predictions that fall within an intergenic region and further from a TSS than the maximum allowed mismatch are ignored. Counting a prediction inside a gene as a FP will be incorrect in some cases, as Carninci et al. (2006) showed that transcription can start in the 5' untranslated region (UTR) of a gene, as well as in an exon or in the 3' UTR. Therefore, not taking into account any of the intergenic predictions or only considering the TSSs located at the 5' UTR gives a biased view of a promoter predictor. Using a gene annotation and ignoring the intergenic predictions will not give a real estimate of the performance of the prediction program for two reasons: (1) predictions inside a gene may in fact be true TSSs; TPs might be underestimated, and consequently FPs are then overestimated; (2) intergenic predictions outside the TSS region are ignored because it is impossible to know whether they are real FPs or whether they point to genes missed by the annotation process. This may in turn lead to underestimation or overestimation of FPs. When using gene annotation with all its uncertainties, it is impossible to know the actual performance of the program in a genomic setting, especially in light of the evidence that much more of the genome is transcribed than just the protein-coding regions (Kapranov et al. 2007; Peters et al. 2007; Prijavic et al. 2007).

To address some of these caveats, we developed an alternative way to assess the performance of a promoter prediction program that does not exhibit these shortcomings. Recently available whole-genome TSS data allow for a new way to count and evaluate. We used a dataset of TSSs for the whole human genome that were characterized experimentally using high-throughput cap-analysis of gene expression (CAGE) (for review, see Carninci 2006). In contrast to previous studies (Bajic et al. 2004, 2006a; Sonnenburg et al. 2006), this dataset allows counting the number of TPs and FPs in an objective way. The dataset contains locations where transcription starts. A TP is a prediction that is within the maximum allowed mismatch from a true TSS, a FN is a true TSS that has no prediction, and a FP is a prediction that is not near a true TSS from the dataset.

## Availability

A user-friendly online implementation of our promoter prediction software EP3 is available at http://bioinformatics.psb.ugent.be/. The program allows the user to process FASTA-formatted sequence files. The online implementation is also available as a standalone application.

## Acknowledgments

We thank Pierre Hilson for fruitful discussions and three anonymous referees for their constructive comments that greatly helped improve the original version of the manuscript. This work was supported by a grant from the Research Foundation–Flanders (3G031805). T.A. thanks the Institute for the Promotion of Innovation by Science and Technology in Flanders for a pre-doctoral fellowship. Y.S. is a Postdoctoral Fellow of the Research Foundation–Flanders.

## References

* Aerts, S., Thijs, G., Dabrowski, M., Moreau, Y., and De Moor, B. 2004. Comprehensive analysis of the base composition around the transcription start site in Metazoa. *BMC Genomics* **5**: 34. doi: 10.1186/1471-2164-5-34.

* Amaldi, F., Bozzoni, I., Beccari, E., and Pierandrei-Amaldi, P. 1989. Expression of ribosomal protein genes and regulation of ribosome biosynthesis in *Xenopus* development. *Trends Biochem. Sci.* **14**: 175-178.

* Baek, D., Davis, C., Ewing, B., Gordon, D., and Green, P. 2007. Characterization and predictive discovery of evolutionarily conserved mammalian alternative promoters. *Genome Res.* **17**: 145-155.

* Bajic, V.B., Seah, S.H., Chong, A., Zhang, G., Koh, J.L.Y., and Brusic, V. 2002. Dragon Promoter Finder: Recognition of vertebrate RNA polymerase II promoters. *Bioinformatics* **18**: 198-199.

* Bajic, V.B., Seah, S.H., Chong, A., Krishnan, S.P.T., Koh, J.L.Y., and Brusic, V. 2003. Computer model for recognition of functional transcription start sites in RNA polymerase II promoters of vertebrates. *J. Mol. Graph. Model.* **21**: 323-332.

* Bajic, V.B., Tan, S.L., Suzuki, Y., and Sugano, S. 2004. Promoter prediction analysis on the whole human genome. *Nat. Biotechnol.* **22**: 1467-1473.

* Bajic, V.B., Brent, M.R., Brown, R.H., Frankish, A., Harrow, J., Ohler, U., Solovyev, V.V., and Tan, S.L. 2006a. Performance assessment of promoter predictions on ENCODE regions in the EGASP experiment. *Genome Biol.* **7**: S3.1-3.13. doi: 10.1186/gb-2006-7-s1-s3.

* Bajic, V.B., Tan, S.L., Christoffels, A., Schonbach, C., Lipovich, L., Yang, L., Hofmann, O., Kruger, A., Hide, W., Kai, C., et al. 2006b. Mice and men: Their promoter properties. *PLoS Genet.* **2**: e54. doi: 10.1371/journal.pgen.0002005.

* Baldi, P., Chauvin, Y., Brunak, S., Gorodkin, J., and Pedersen, A.G. 1998. Computational applications of DNA structural scales. *Proc. Int. Conf. Intell. Syst. Mol. Biol.* **6**: 35-42.

* Bartel, D.P. 2004. MicroRNAs: Genomics, biogenesis, mechanism, and function. *Cell* **116**: 281-297.

* Bird, A. 2002. DNA methylation patterns and epigenetic memory. *Genes & Dev.* **16**: 6-21.
Blake, R.D. and Delcourt, S.G. 1998. Thermal stability of DNA. *Nucleic Acids Res.* **26**: 3323–3332.

Blake, R.D., Bizzaro, J.W., Blake, J.D., Day, G.R., Delcourt, S.G., Knowles, J., Marx, K.A., and SantaLucia Jr., J. 1999. Statistical mechanical simulation of polymeric DNA melting with MELTSIM. *Bioinformatics* **15**: 370–375.

Bode, J., Winkelmann, S., Gotze, S., Spiker, S., Tsutsui, K., Bi, C., Prashanth, A.K., and Benham, C. 2006. Correlations between scaffold/matrix attachment region (S/MAR) binding activity and DNA duplex destabilization energy. *J. Mol. Biol.* **358**: 597–613.

Breslauer, K.J., Frank, R., Blocker, H., and Marky, L.A. 1986. Predicting DNA duplex stability from the base sequence. *Proc. Natl. Acad. Sci.* **83**: 3746–3750.

Brukner, I., Sanchez, R., Suck, D., and Pongor, S. 1995. Trinucleotide models for DNA bending propensity: Comparison of models based on DNasel digestion and nucleosome packaging data. *J. Biomol. Struct. Dyn.* **13**: 309–317.

Burden, S., Lin, Y.X., and Zhang, R. 2004. Improving promoter prediction for the NNPP2.2 algorithm: A case study using *Escherichia coli* DNA sequences. *Bioinformatics* **21**: 601–607.

Butler, J.E.F. and Kadonaga, J.T. 2002. The RNA polymerase II core promoter: A key component in the regulation of gene expression. *Genes & Dev.* **16**: 2583–2592.

Carninci, P. 2006. Tagging mammalian transcription complexity. *Trends Genet.* **22**: 501–510.

Carninci, P., Sandelin, A., Lenhard, B., Katayama, S., Shimokawa, K., Ponjavic, J., Semple, C.A.M., Taylor, M.S., Engström, P.G., Frith, M.C., et al. 2006. Genome-wide analysis of mammalian promoter architecture and evolution. *Nat. Genet.* **38**: 626–635.

Chen, Q.K., Hertz, G.Z., and Stormo, G.D. 1997. PromFD 1.0: A computer program that predicts eukaryotic pol II promoters using strings and IMD matrices. *Comput. Appl. Biosci.* **13**: 29–35.

Choi, C.H., Kalosakas, G., Rasmussen, K.O., Hiromura, M., Bishop, A.R., and Usheva, A. 2004. DNA dynamically directs its own transcription initiation. *Nucleic Acids Res.* **32**: 1584–1590.

Claverie, J.M., Poirot, O., and Lopez, F. 1997. The difficulty of identifying genes in anonymous vertebrate sequences. *Comput. Chem.* **21**: 203–214.

Comai, L., Tanese, N., and Tjian, R. 1992. The TATA-binding protein and associated factors are integral components of the RNA polymerase I transcription factor, SL1. *Cell* **68**: 965–976.

Cooper, S.J., Trinklein, N.D., Anton, E.D., Nguyen, L., and Myers, R.M. 2006. Comprehensive analysis of transcriptional promoter structure and function in 1% of the human genome. *Genome Res.* **16**: 1–10.

Cormack, B.P. and Struhl, K. 1992. The TATA-binding protein is required for transcription by all three nuclear RNA polymerases in yeast cells. *Cell* **69**: 685–696.

Davuluri, R.V., Grosse, I., and Zhang, M.Q. 2001. Computational identification of promoters and first exons in the human genome. *Nat. Genet.* **29**: 412–417.

Delgado, S., Gómez, M., Bird, A., and Antequera, F. 1998. Initiation of DNA replication at CpG islands in mammalian chromosomes. *EMBO J.* **17**: 2426–2435.

Deng, W. and Roberts, S.G. 2005. A core promoter element downstream of the TATA box that is recognized by TFIIB. *Genes & Dev.* **19**: 2418–2423.

Deng, W. and Roberts, S.G. 2006. Core promoter elements recognized by transcription factor IIB. *Biochem. Soc. Trans.* **34**: 1051–1053.

Down, T.A. and Hubbard, T.J.P. 2002. Computational detection and location of transcription start sites in mammalian genomic DNA. *Genome Res.* **12**: 458–461.

Eddy, S.R. 2001. Non-coding RNA genes and the modern RNA world. *Nat. Rev. Genet.* **2**: 919–929.

El Hassan, M.A. and Calladine, C.R. 1996. Propeller-twisting of base-pairs and the conformational mobility of dinucleotide steps in DNA. *J. Mol. Biol.* **259**: 95–103.

Faiger, H., Ivanchenko, M., Cohen, I., and Haran, T.E. 2006. TBP flanking sequences: Asymmetry of binding, long-range effects and consensus sequences. *Nucleic Acids Res.* **34**: 104–119.

Fickett, J.W. and Hatzigeorgiou, A.G. 1997. Eukaryotic promoter recognition. *Genome Res.* **7**: 861–878.

Fitzgerald, P., Sturgill, D., Shyakhtenko, A., Oliver, B., and Vinson, C. 2006. Comparative genomics of *Drosophila* and human core promoters. *Genome Biol.* **7**: R53. doi: 10.1186/gb-2006-7-7-r53.

Florquin, K., Saeks, Y., Degroeve, S., Rouzé, P., and Van de Peer, Y. 2005. Large-scale structural analysis of the core promoter in mammalian and plant genomes. *Nucleic Acids Res.* **33**: 4255–4264.

Fukue, Y., Sumida, N., Nishikawa, J.-i., and Ohyama, T. 2004. Core promoter elements of eukaryotic genes have a highly distinctive mechanical property. *Nucleic Acids Res.* **32**: 5834–5840.

Fukue, Y., Sumida, N., Tanase, J.-i., and Ohyama, T. 2005. A highly

distinctive mechanical property found in the majority of human promoters and its transcriptional relevance. *Nucleic Acids Res.* **33**: 3821–3827.

Gorin, A.A., Zhurkin, V.B., and Olson, W.K. 1995. B-DNA twisting correlates with base-pair morphology. *J. Mol. Biol.* **247**: 34–48.

Greenbaum, J.A., Pang, B., and Tullius, T.D. 2007. Construction of a genome-scale structural map at single-nucleotide resolution. *Genome Res.* **17**: 947–953.

Gross, S.S. and Brent, M.R. 2006. Using multiple alignments to improve gene prediction. *J. Comput. Biol.* **13**: 379–393.

Hannenhalli, S. and Levy, S. 2001. Promoter prediction in the human genome. *Bioinformatics* **17**: S90–S96.

Harrow, J., Denoeud, F., Frankish, A., Reymond, A., Chen, C.K., Crast, J., Lagarde, J., Gilbert, J.G., Storey, R., Swarbreck, D., et al. 2006. Gencode: Producing a reference annotation for ENCODE. *Genome Biol.* **7**: 1–9.

Ho, P.S., Zhou, G.W., and Clark, L.B. 1990. Polarized electronic spectra of Z-DNA single crystals. *Biopolymers* **30**: 151–163.

Hubbard, T.J., Aken, B.L., Beal, K., Ballester, B., Caccamo, M., Chen, Y., Clarke, L., Coates, G., Cunningham, F., Cutts, T., et al. 2007. Ensembl 2007. *Nucleic Acids Res.* **35**: D610–D617.

Ioshikhes, I.P. and Zhang, M.Q. 2000. Large-scale human promoter mapping using CpG islands. *Nat. Genet.* **26**: 61–63.

Ivanov, V.I. and Minchenkova, L.E. 1994. The A-form of DNA: In search of the biological role. *Mol. Biol. (Mosk.)* **28**: 1258–1271.

Jacob, S.T. 1995. Regulation of ribosomal gene transcription. *Biochem. J.* **306**: 617–626.

Jin, V.X., Singer, G.A.C., Agosto-Pérez, F.J., Liyanarachchi, S., and Davuluri, R.V. 2006. Genome-wide analysis of core promoter elements from conserved human and mouse orthologous pairs. *BMC Bioinformatics* **7**: 114. doi: 10.1186/1471-2105-7-114.

Kanhere, A. and Bansal, M. 2005. Structural properties of promoters: Similarities and differences between prokaryotes and eukaryotes. *Nucleic Acids Res.* **33**: 3165–3175.

Kapranov, P., Willingham, A.T., and Gingeras, T.R. 2007. Genome-wide transcription and the implications for genomic organization. *Nat. Rev. Genet.* **8**: 413–423.

Kawaji, H., Frith, M.C., Katayama, S., Sandelin, A., Kai, C., Kawai, J., Carninci, P., and Hayashizaki, Y. 2006. Dynamic usage of transcription start sites within core promoters. *Genome Biol.* **7**: R118. doi: 10.1186/gb-2006-7-12-r118.

Knudsen, S. 1999. Promoter2.0: For the recognition of PolII promoter sequences. *Bioinformatics* **15**: 356–361.

Kuhn, R.M., Karolchik, D., Zweig, A.S., Trumbower, H., Thomas, D.J., Thakkapallayil, A., Sugnet, C.W., Stanke, M., Smith, K.E., Siepel, A., et al. 2007. The UCSC genome browser database: Update 2007. *Nucleic Acids Res.* **35**: D668–D673.

Lagrange, T., Kapanidis, A.N., Tang, H., Reinberg, D., and Ebright, R.H. 1998. New core promoter element in RNA polymerase II-dependent transcription: Sequence-specific DNA binding by transcription factor IIB. *Genes & Dev.* **12**: 34–44.

Liao, G-c., Rehm, E.J., and Rubin, G.M. 2000. Insertion site preferences of the P transposable element in *Drosophila melanogaster*. *Proc. Natl. Acad. Sci.* **97**: 3347–3351.

Liolios, K., Tavernarakis, N., Hugenholtz, P., and Kyrpides, N.C. 2006. The Genomes On Line Database (GOLD) v.2: A monitor of genome projects worldwide. *Nucleic Acids Res.* **34**: D332–D334.

Liu, R. and States, D.J. 2002. Consensus promoter identification in the human genome utilizing expressed gene markers and gene modeling. *Genome Res.* **12**: 462–469.

Liu, Y., Liu, X.S., Wei, L., Altman, R.B., and Batzoglou, S. 2004. Eukaryotic regulatory element conservation analysis and identification using comparative genomics. *Genome Res.* **14**: 451–458.

Liu, F., Tostesen, E., Sundet, J.K., Jenssen, T.K., Bock, C., Jerstad, G.I., Thilly, W.G., and Hovig, E. 2007. The human genomic melting map. *PLoS Comput. Biol.* **3**: e93. doi: 10.1371/journal.pcbi.0030093.

Lo, K. and Smale, S.T. 1996. Generality of a functional initiator consensus sequence. *Gene* **182**: 13–22.

Lowe, T.M. and Eddy, S.R. 1997. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. *Nucleic Acids Res.* **25**: 955–964.

Mathé, C., Sagot, M.-F., Schiex, T., and Rouzé, P. 2002. Current methods of gene prediction, their strengths and weaknesses. *Nucleic Acids Res.* **30**: 4103–4117.

Merling, A., Sagaydakova, N., and Haran, T.E. 2003. A-tract polarity dominate the curvature in flanking sequences. *Biochemistry* **42**: 4978–4984.

Munch, K. and Krogh, A. 2006. Automatic generation of gene finders for eukaryotic species. *BMC Bioinformatics* **7**: 263. doi: 10.1186/1471-2105-7-263.
Blake, R.D. and Delcourt, S.G. 1998. Thermal stability of DNA. *Nucleic Acids Res.* **26**: 3323-3332.

Blake, R.D., Bizzaro, J.W., Blake, J.D., Day, G.R., Delcourt, S.G., Knowles, J., Marx, K.A., and SantaLucia Jr., J. 1999. Statistical mechanical simulation of polymeric DNA melting with MELTSIM. *Bioinformatics* **15**: 370-375.

Bode, J., Winkelmann, S., Gotze, S., Spiker, S., Tsutsui, K., Bi, C., Prashanth, A.K., and Benham, C. 2006. Correlations between scaffold/matrix attachment region (S/MAR) binding activity and DNA duplex destabilization energy. *J. Mol. Biol.* **358**: 597-613.

Breslauer, K.J., Frank, R., Blocker, H., and Marky, L.A. 1986. Predicting DNA duplex stability from the base sequence. *Proc. Natl. Acad. Sci.* **83**: 3746-3750.

Brukner, I., Sanchez, R., Suck, D., and Pongor, S. 1995. Trinucleotide models for DNA bending propensity: Comparison of models based on DNasel digestion and nucleosome packaging data. *J. Biomol. Struct. Dyn.* **13**: 309-317.

Burden, S., Lin, Y.X., and Zhang, R. 2004. Improving promoter prediction for the NNPP2.2 algorithm: A case study using *Escherichia coli* DNA sequences. *Bioinformatics* **21**: 601-607.

Butler, J.E.F. and Kadonaga, J.T. 2002. The RNA polymerase II core promoter: A key component in the regulation of gene expression. *Genes & Dev.* **16**: 2583-2592.

Carninci, P. 2006. Tagging mammalian transcription complexity. *Trends Genet.* **22**: 501-510.

Carninci, P., Sandelin, A., Lenhard, B., Katayama, S., Shimokawa, K., Ponjavic, J., Semple, C.A.M., Taylor, M.S., Engström, P.G., Frith, M.C., et al. 2006. Genome-wide analysis of mammalian promoter architecture and evolution. *Nat. Genet.* **38**: 626-635.

Chen, Q.K., Hertz, G.Z., and Stormo, G.D. 1997. PromFD 1.0: A computer program that predicts eukaryotic pol II promoters using strings and IMD matrices. *Comput. Appl. Biosci.* **13**: 29-35.

Choi, C.H., Kalosakas, G., Rasmussen, K.O., Hiromura, M., Bishop, A.R., and Usheva, A. 2004. DNA dynamically directs its own transcription initiation. *Nucleic Acids Res.* **32**: 1584-1590.

Claverie, J.M., Poirot, O., and Lopez, F. 1997. The difficulty of identifying genes in anonymous vertebrate sequences. *Comput. Chem.* **21**: 203-214.

Comai, L., Tanese, N., and Tjian, R. 1992. The TATA-binding protein and associated factors are integral components of the RNA polymerase I transcription factor, SL1. *Cell* **68**: 965-976.

Cooper, S.J., Trinklein, N.D., Anton, E.D., Nguyen, L., and Myers, R.M. 2006. Comprehensive analysis of transcriptional promoter structure and function in 1% of the human genome. *Genome Res.* **16**: 1-10.

Cormack, B.P. and Struhl, K. 1992. The TATA-binding protein is required for transcription by all three nuclear RNA polymerases in yeast cells. *Cell* **69**: 685-696.

Davuluri, R.V., Grosse, I., and Zhang, M.Q. 2001. Computational identification of promoters and first exons in the human genome. *Nat. Genet.* **29**: 412-417.

Delgado, S., Gómez, M., Bird, A., and Antequera, F. 1998. Initiation of DNA replication at CpG islands in mammalian chromosomes. *EMBO J.* **17**: 2426-2435.

Deng, W. and Roberts, S.G. 2005. A core promoter element downstream of the TATA box that is recognized by TFIIB. *Genes & Dev.* **19**: 2418-2423.

Deng, W. and Roberts, S.G. 2006. Core promoter elements recognized by transcription factor IIB. *Biochem. Soc. Trans.* **34**: 1051-1053.

Down, T.A. and Hubbard, T.J.P. 2002. Computational detection and location of transcription start sites in mammalian genomic DNA. *Genome Res.* **12**: 458-461.

Eddy, S.R. 2001. Non-coding RNA genes and the modern RNA world. *Nat. Rev. Genet.* **2**: 919-929.

El Hassan, M.A. and Calladine, C.R. 1996. Propeller-twisting of base-pairs and the conformational mobility of dinucleotide steps in DNA. *J. Mol. Biol.* **259**: 95-103.

Faiger, H., Ivanchenko, M., Cohen, I., and Haran, T.E. 2006. TBP flanking sequences: Asymmetry of binding, long-range effects and consensus sequences. *Nucleic Acids Res.* **34**: 104-119.

Fickett, J.W. and Hatzigeorgiou, A.G. 1997. Eukaryotic promoter recognition. *Genome Res.* **7**: 861-878.

Fitzgerald, P., Sturgill, D., Shyakhtenko, A., Oliver, B., and Vinson, C. 2006. Comparative genomics of *Drosophila* and human core promoters. *Genome Biol.* **7**: R53. doi: 10.1186/gb-2006-7-7-r53.

Florquin, K., Saeks, Y., Degroeve, S., Rouzé, P., and Van de Peer, Y. 2005. Large-scale structural analysis of the core promoter in mammalian and plant genomes. *Nucleic Acids Res.* **33**: 4255-4264.

Fukue, Y., Sumida, N., Nishikawa, J.-i., and Ohyama, T. 2004. Core promoter elements of eukaryotic genes have a highly distinctive mechanical property. *Nucleic Acids Res.* **32**: 5834-5840.

Fukue, Y., Sumida, N., Tanase, J.-i., and Ohyama, T. 2005. A highly

distinctive mechanical property found in the majority of human promoters and its transcriptional relevance. *Nucleic Acids Res.* **33**: 3821-3827.

Gorin, A.A., Zhurkin, V.B., and Olson, W.K. 1995. B-DNA twisting correlates with base-pair morphology. *J. Mol. Biol.* **247**: 34-48.

Greenbaum, J.A., Pang, B., and Tullius, T.D. 2007. Construction of a genome-scale structural map at single-nucleotide resolution. *Genome Res.* **17**: 947-953.

Gross, S.S. and Brent, M.R. 2006. Using multiple alignments to improve gene prediction. *J. Comput. Biol.* **13**: 379-393.

Hannenhalli, S. and Levy, S. 2001. Promoter prediction in the human genome. *Bioinformatics* **17**: S90-S96.

Harrow, J., Denoeud, F., Frankish, A., Reymond, A., Chen, C.K., Crast, J., Lagarde, J., Gilbert, J.G., Storey, R., Swarbreck, D., et al. 2006. Gencode: Producing a reference annotation for ENCODE. *Genome Biol.* **7**: 1-9.

Ho, P.S., Zhou, G.W., and Clark, L.B. 1990. Polarized electronic spectra of Z-DNA single crystals. *Biopolymers* **30**: 151-163.

Hubbard, T.J., Aken, B.L., Beal, K., Ballester, B., Caccamo, M., Chen, Y., Clarke, L., Coates, G., Cunningham, F., Cutts, T., et al. 2007. Ensembl 2007. *Nucleic Acids Res.* **35**: D610-D617.

Ioshikhes, I.P. and Zhang, M.Q. 2000. Large-scale human promoter mapping using CpG islands. *Nat. Genet.* **26**: 61-63.

Ivanov, V.I. and Minchenkova, L.E. 1994. The A-form of DNA: In search of the biological role. *Mol. Biol. (Mosk.)* **28**: 1258-1271.

Jacob, S.T. 1995. Regulation of ribosomal gene transcription. *Biochem. J.* **306**: 617-626.

Jin, V.X., Singer, G.A.C., Agosto-Pérez, F.J., Liyanarachchi, S., and Davuluri, R.V. 2006. Genome-wide analysis of core promoter elements from conserved human and mouse orthologous pairs. *BMC Bioinformatics* **7**: 114. doi: 10.1186/1471-2105-7-114.

Kanhere, A. and Bansal, M. 2005. Structural properties of promoters: Similarities and differences between prokaryotes and eukaryotes. *Nucleic Acids Res.* **33**: 3165-3175.

Kapranov, P., Willingham, A.T., and Gingeras, T.R. 2007. Genome-wide transcription and the implications for genomic organization. *Nat. Rev. Genet.* **8**: 413-423.

Kawaji, H., Frith, M.C., Katayama, S., Sandelin, A., Kai, C., Kawai, J., Carninci, P., and Hayashizaki, Y. 2006. Dynamic usage of transcription start sites within core promoters. *Genome Biol.* **7**: R118. doi: 10.1186/gb-2006-7-12-r118.

Knudsen, S. 1999. Promoter2.0: For the recognition of PolII promoter sequences. *Bioinformatics* **15**: 356-361.

Kuhn, R.M., Karolchik, D., Zweig, A.S., Trumbower, H., Thomas, D.J., Thakkapallayil, A., Sugnet, C.W., Stanke, M., Smith, K.E., Siepel, A., et al. 2007. The UCSC genome browser database: Update 2007. *Nucleic Acids Res.* **35**: D668-D673.

Lagrange, T., Kapanidis, A.N., Tang, H., Reinberg, D., and Ebright, R.H. 1998. New core promoter element in RNA polymerase II-dependent transcription: Sequence-specific DNA binding by transcription factor IIB. *Genes & Dev.* **12**: 34-44.

Liao, G-c., Rehm, E.J., and Rubin, G.M. 2000. Insertion site preferences of the P transposable element in *Drosophila melanogaster*. *Proc. Natl. Acad. Sci.* **97**: 3347-3351.

Liolios, K., Tavernarakis, N., Hugenholtz, P., and Kyrpides, N.C. 2006. The Genomes On Line Database (GOLD) v.2: A monitor of genome projects worldwide. *Nucleic Acids Res.* **34**: D332-D334.

Liu, R. and States, D.J. 2002. Consensus promoter identification in the human genome utilizing expressed gene markers and gene modeling. *Genome Res.* **12**: 462-469.

Liu, Y., Liu, X.S., Wei, L., Altman, R.B., and Batzoglou, S. 2004. Eukaryotic regulatory element conservation analysis and identification using comparative genomics. *Genome Res.* **14**: 451-458.

Liu, F., Tostesen, E., Sundet, J.K., Jenssen, T.K., Bock, C., Jerstad, G.I., Thilly, W.G., and Hovig, E. 2007. The human genomic melting map. *PLoS Comput. Biol.* **3**: e93. doi: 10.1371/journal.pcbi.0030093.

Lo, K. and Smale, S.T. 1996. Generality of a functional initiator consensus sequence. *Gene* **182**: 13-22.

Lowe, T.M. and Eddy, S.R. 1997. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. *Nucleic Acids Res.* **25**: 955-964.

Mathé, C., Sagot, M.-F., Schiex, T., and Rouzé, P. 2002. Current methods of gene prediction, their strengths and weaknesses. *Nucleic Acids Res.* **30**: 4103-4117.

Merling, A., Sagaydakova, N., and Haran, T.E. 2003. A-tract polarity dominate the curvature in flanking sequences. *Biochemistry* **42**: 4978-4984.

Munch, K. and Krogh, A. 2006. Automatic generation of gene finders for eukaryotic species. *BMC Bioinformatics* **7**: 263. doi: 10.1186/1471-2105-7-263.
Nikolov, D.B., Chen, H., Halay, E.D., Usheva, A.A., Hisatake, K., Lee, D.K., Roeder, R.G., and Burley, S.K. 1995. Crystal structure of a TFIIIB-TBP-TATA-element ternary complex. *Nature* **377**: 119–128.

Ohler, U. 2000. Promoter prediction on a genomic scale—The Adh experience. *Genome Res.* **10**: 539–542.

Ohler, U. 2001. *Computational promoter recognition in eukaryotic genomic DNA*. Technische Fakultät der Universität Erlangen, Nürnberg.

Ohler, U. 2006. Identification of core promoter modules in *Drosophila* and their application in accurate transcription start site prediction. *Nucleic Acids Res.* **34**: 5943–5950.

Ohler, U. and Niemann, H. 2001. Identification and analysis of eukaryotic promoters: Recent computational approaches. *Trends Genet.* **17**: 56–60.

Ohler, U., Stemmer, G., Harbeck, S., and Niemann, H. 2000. Stochastic segment models of eukaryotic promoter regions. *Pac. Symp. Biocomput.* **5**: 377–388.

Ohler, U., Niemann, H., Liao, G.-c., and Rubin, G.M. 2001. Joint modeling of DNA sequence and physical properties to improve eukaryotic promoter recognition. *Bioinformatics* **17**: S199–S206.

Ohler, U., Liao, G.-c., Niemann, H., and Rubin, G.M. 2002. Computational analysis of core promoters in the *Drosophila* genome. *Genome Biol.* **3**: doi: 10.1186/gb-2002-3-12-research0087.

Olson, W.K., Gorin, A.A., Lu, X.-J., Hock, L.M., and Zhurkin, V.B. 1998. DNA sequence-dependent deformability deduced from protein-DNA crystal complexes. *Proc. Natl. Acad. Sci.* **95**: 11163–11168.

Ornstein, R.L., Rein, R., Breen, D.L., and MacElroy, R.D. 1978. Optimized potential function for calculation of nucleic-acid interaction energies. 1. Base stacking. *Biopolymers* **17**: 2341–2360.

Pedersen, A.G., Baldi, P., Chauvin, Y., and Brunak, S.O. 1998. DNA structure in human RNA polymerase II promoters. *J. Mol. Biol.* **281**: 663–673.

Pedersen, A.G., Baldi, P., Chauvin, Y., and Brunak, S.O. 1999. The biology of eukaryotic promoter prediction—A review. *Comput. Chem.* **23**: 191–207.

Peters, B.A., St Croix, B., Sjöblom, T., Cummins, J.M., Silliman, N., Ptak, J., Saha, S., Kinzler, K.W., Hatzis, C., and Velculescu, V.E. 2007. Large-scale identification of novel transcripts in the human genome. *Genome Res.* **17**: 287–292.

Ponger, L. and Mouchiroud, D. 2002. CpGProD: Identifying CpG islands associated with transcription start sites in large genomic mammalian sequences. *Bioinformatics* **18**: 631–633.

Ponjavic, J., Ponting, C.P., and Lunter, G. 2007. Functionality or transcriptional noise? Evidence for selection within long noncoding RNAs. *Genome Res.* **17**: 556–565.

Prestridge, D.S. 1995. Predicting Pol II promoter sequences using transcription factor binding sites. *J. Mol. Biol.* **249**: 923–932.

Presutti, C., Ciafre, S.A., and Bozzoni, I. 1991. The ribosomal protein L2 in *S. cerevisiae* controls the level of accumulation of its own mRNA. *EMBO J.* **10**: 2215–2221.

Ratsch, G., Sonnenburg, S., and Schafer, C. 2006. Learning interpretable SVMs for biological sequence classification. *BMC Bioinformatics* **7**: S9. doi: 10.1186/1471-2105-7-S1-S9.

Reese, M.G. 2001. Application of a time-delay neural network to promoter annotation in the *Drosophila melanogaster* genome. *Comput. Chem.* **26**: 51–56.

Rombauts, S., Florquin, K., Lescot, M., Marchal, K., Rouzé, P., and Van de Peer, Y. 2003. Computational approaches to identify promoters and cis-regulatory elements in plant genomes. *Plant Physiol.* **132**: 1162–1176.

Rozowsky, J.S., Newburger, D., Sayward, F., Wu, J., Jordan, G., Korbel, J.O., Nagalakshmi, U., Yang, J., Zheng, D., Guigo, R., et al. 2007. The DART classification of unannotated transcription within the ENCODE regions: Associating transcription with known and novel loci. *Genome Res.* **17**: 732–745.

Satchwell, S.C., Drew, H.R., and Travers, A.A. 1986. Sequence periodicities in chicken nucleosome core DNA. *J. Mol. Biol.* **191**: 659–675.

Scherf, M., Klingenhoff, A., and Werner, T. 2000. Highly specific localization of promoter regions in large genomic sequences by PromoterInspector: A novel context analysis approach. *J. Mol. Biol.* **297**: 599–606.

Scherf, M., Klingenhoff, A., Frech, K., Quandt, K., Schneider, R., Grote, K., Frisch, M., Gailus-Durner, V., Seidel, A., Brack-Werner, R., et al. 2001. First pass annotation of promoters on human chromosome 22. *Genome Res.* **11**: 333–340.

Schmid, C.D., Perier, R., Praz, V., and Bucher, P. 2006. EPD in its twentieth year: Towards complete promoter coverage of selected model organisms. *Nucleic Acids Res.* **34**: D82–D85.

Shiraki, T., Kondo, S., Katayama, S., Waki, K., Kasukawa, T., Kawaji, H., Kodzius, R., Watahiki, A., Nakamura, M., Arakawa, T., et al. 2003. Cap analysis gene expression for high-throughput analysis of transcriptional starting point and identification of promoter usage. *Proc. Natl. Acad. Sci.* **100**: 15776–15781.

Sivolob, A.V. and Khrapunov, S.N. 1995. Translational positioning of nucleosomes on DNA: The role of sequence-dependent isotropic DNA bending stiffness. *J. Mol. Biol.* **247**: 918–931.

Smale, S.T. 2001. Core promoters: Active contributors to combinatorial gene regulation. *Genes & Dev.* **15**: 2503–2508.

Smale, S.T. and Kadonaga, J.T. 2003. The RNA polymerase II core promoter. *Annu. Rev. Biochem.* **72**: 449–479.

Solovyev, V.V. and Shahmuradov, I.A. 2003. PromH: Promoters identification using orthologous genomic sequences. *Nucleic Acids Res.* **31**: 3540–3545.

Solovyev, V., Kosarev, P., Seledsov, I., and Vorobyev, D. 2006. Automatic annotation of eukaryotic genes, pseudogenes and promoters. *Genome Biol.* **7**: 11–12.

Sonnenburg, S.O., Zien, A., and Rätsch, G. 2006. ARTS: Accurate recognition of transcription starts in human. *Bioinformatics* **22**: 472–480.

Struhl, K. 1989. Molecular mechanisms of transcriptional regulation in yeast. *Annu. Rev. Biochem.* **58**: 1051–1077.

Sugimoto, N., Nakano, S., Yoneyama, M., and Honda, K. 1996. Improved thermodynamic parameters and helix initiation factor to predict stability of DNA duplexes. *Nucleic Acids Res.* **24**: 4501–4505.

Suzuki, Y., Tsunoda, T., Sese, J., Taira, H., Mizushima-Sugano, J., Hata, H., Ota, T., Isogai, T., Tanaka, T., Nakamura, Y., et al. 2001. Identification and characterization of the potential promoter regions of 1031 kinds of human genes. *Genome Res.* **11**: 677–684.

Tabach, Y., Brosh, R., Buganim, Y., Reiner, A., Zuk, O., Yitzhaky, A., Koudritsky, M., Rotter, V., and Domany. E. 2007. Wide-scale analysis of human functional transcription factor binding reveals a strong bias towards the transcription start site. *PLoS ONE* **2**: e807. doi: 10.1371/journal.pone.0000807.

Trinklein, N.D., Aldred, S.F., Hartman, S.J., Schroeder, D.I., Otillar, R.P., and Myers, R.M. 2004. An abundance of bidirectional promoters in the human genome. *Genome Res.* **14**: 62–66.

Uren, P., Cameron-Jones, M., and Sale, A. 2006. Promoter prediction using physic-chemical properties of DNA. *Lect. Notes Comput. Sci.* **4216**: 21–31.

van Rijssbergen, C.J. 1979. *Information retrieval, 2nd edition*. Butterworths, London.

Wang, H. and Benham, C.J. 2006. Promoter prediction and annotation of microbial genomes based on DNA sequence and structural responses to superhelical stress. *BMC Bioinformatics* **7**: 248. doi: 10.1186/1471-2105-7-248.

Wang, G. and Zhang, W. 2006. A steganalysis-based approach to comprehensive identification and characterization of functional regulatory elements. *Genome Biol.* **7**: R49. doi: 10.1186/gb-2006-7-6-r49.

Wang, Z., Chen, Y., and Li, Y. 2004. A brief review of computational gene prediction methods. *Genom. Proteom. Bioinformatics* **2**: 216–221.

Wasserman, W.W., Palumbo, M., Thompson, W., Fickett, J.W., and Lawrence, C.E. 2000. Human-mouse genome comparisons to locate regulatory sites. *Nat. Genet.* **26**: 225–228.

Werner, T. 2003. The state of the art in mammalian promoter recognition. *Brief. Bioinform.* **4**: 22–30.

White, R.J. and Jackson, S.P. 1992. The TATA-binding protein: A central role in transcription by RNA polymerases I, II and III. *Trends Genet.* **8**: 284–288.

Xie, X., Wu, S., Lam, K.M., and Yan, H. 2006. PromoterExplorer: An effective promoter identification method based on the AdaBoost algorithm. *Bioinformatics* **22**: 2722–2728.

Zhang, M.Q. 2002. Computational prediction of eukaryotic protein-coding genes. *Nat. Rev. Genet.* **3**: 698–709.

Zhou, X., Ruan, J., Wang, G., and Zhang, W. 2007. Characterization and identification of microRNA core promoters in four model species. *PLoS Comput. Biol.* **3**: e37. doi: 10.1371/journal.pcbi.0030037.

Received August 3, 2007; accepted in revised form November 14, 2007.
Nikolov, D.B., Chen, H., Halay, E.D., Usheva, A.A., Hisatake, K., Lee, D.K., Roeder, R.G., and Burley, S.K. 1995. Crystal structure of a TFIIIB-TBP-TATA-element ternary complex. *Nature* **377**: 119–128.

Ohler, U. 2000. Promoter prediction on a genomic scale—The Adh experience. *Genome Res.* **10**: 539–542.

Ohler, U. 2001. *Computational promoter recognition in eukaryotic genomic DNA*. Technische Fakultät der Universität Erlangen, Nürnberg.

Ohler, U. 2006. Identification of core promoter modules in *Drosophila* and their application in accurate transcription start site prediction. *Nucleic Acids Res.* **34**: 5943–5950.

Ohler, U. and Niemann, H. 2001. Identification and analysis of eukaryotic promoters: Recent computational approaches. *Trends Genet.* **17**: 56–60.

Ohler, U., Stemmer, G., Harbeck, S., and Niemann, H. 2000. Stochastic segment models of eukaryotic promoter regions. *Pac. Symp. Biocomput.* **5**: 377–388.

Ohler, U., Niemann, H., Liao, G.-c., and Rubin, G.M. 2001. Joint modeling of DNA sequence and physical properties to improve eukaryotic promoter recognition. *Bioinformatics* **17**: S199–S206.

Ohler, U., Liao, G.-c., Niemann, H., and Rubin, G.M. 2002. Computational analysis of core promoters in the *Drosophila* genome. *Genome Biol.* **3**: doi: 10.1186/gb-2002-3-12-research0087.

Olson, W.K., Gorin, A.A., Lu, X.-J., Hock, L.M., and Zhurkin, V.B. 1998. DNA sequence-dependent deformability deduced from protein-DNA crystal complexes. *Proc. Natl. Acad. Sci.* **95**: 11163–11168.

Ornstein, R.L., Rein, R., Breen, D.L., and MacElroy, R.D. 1978. Optimized potential function for calculation of nucleic-acid interaction energies. 1. Base stacking. *Biopolymers* **17**: 2341–2360.

Pedersen, A.G., Baldi, P., Chauvin, Y., and Brunak, S.O. 1998. DNA structure in human RNA polymerase II promoters. *J. Mol. Biol.* **281**: 663–673.

Pedersen, A.G., Baldi, P., Chauvin, Y., and Brunak, S.O. 1999. The biology of eukaryotic promoter prediction—A review. *Comput. Chem.* **23**: 191–207.

Peters, B.A., St Croix, B., Sjöblom, T., Cummins, J.M., Silliman, N., Ptak, J., Saha, S., Kinzler, K.W., Hatzis, C., and Velculescu, V.E. 2007. Large-scale identification of novel transcripts in the human genome. *Genome Res.* **17**: 287–292.

Ponger, L. and Mouchiroud, D. 2002. CpGProD: Identifying CpG islands associated with transcription start sites in large genomic mammalian sequences. *Bioinformatics* **18**: 631–633.

Ponjavic, J., Ponting, C.P., and Lunter, G. 2007. Functionality or transcriptional noise? Evidence for selection within long noncoding RNAs. *Genome Res.* **17**: 556–565.

Pretridge, D.S. 1995. Predicting Pol II promoter sequences using transcription factor binding sites. *J. Mol. Biol.* **249**: 923–932.

Presutti, C., Ciafre, S.A., and Bozzoni, I. 1991. The ribosomal protein L2 in *S. cerevisiae* controls the level of accumulation of its own mRNA. *EMBO J.* **10**: 2215–2221.

Ratsch, G., Sonnenburg, S., and Schafer, C. 2006. Learning interpretable SVMs for biological sequence classification. *BMC Bioinformatics* **7**: S9. doi: 10.1186/1471-2105-7-S1-S9.

Reese, M.G. 2001. Application of a time-delay neural network to promoter annotation in the *Drosophila melanogaster* genome. *Comput. Chem.* **26**: 51–56.

Rombauts, S., Florquin, K., Lescot, M., Marchal, K., Rouzé, P., and Van de Peer, Y. 2003. Computational approaches to identify promoters and cis-regulatory elements in plant genomes. *Plant Physiol.* **132**: 1162–1176.

Rozowsky, J.S., Newburger, D., Sayward, F., Wu, J., Jordan, G., Korbel, J.O., Nagalakshmi, U., Yang, J., Zheng, D., Guigo, R., et al. 2007. The DART classification of unannotated transcription within the ENCODE regions: Associating transcription with known and novel loci. *Genome Res.* **17**: 732–745.

Satchwell, S.C., Drew, H.R., and Travers, A.A. 1986. Sequence periodicities in chicken nucleosome core DNA. *J. Mol. Biol.* **191**: 659–675.

Scherf, M., Klingenhoff, A., and Werner, T. 2000. Highly specific localization of promoter regions in large genomic sequences by PromoterInspector: A novel context analysis approach. *J. Mol. Biol.* **297**: 599–606.

Scherf, M., Klingenhoff, A., Frech, K., Quandt, K., Schneider, R., Grote, K., Frisch, M., Gailus-Durner, V., Seidel, A., Brack-Werner, R., et al. 2001. First pass annotation of promoters on human chromosome 22. *Genome Res.* **11**: 333–340.

Schmid, C.D., Perier, R., Praz, V., and Bucher, P. 2006. EPD in its twentieth year: Towards complete promoter coverage of selected model organisms. *Nucleic Acids Res.* **34**: D82–D85.

Shiraki, T., Kondo, S., Katayama, S., Waki, K., Kasukawa, T., Kawaji, H., Kodzius, R., Watahiki, A., Nakamura, M., Arakawa, T., et al. 2003. Cap analysis gene expression for high-throughput analysis of transcriptional starting point and identification of promoter usage. *Proc. Natl. Acad. Sci.* **100**: 15776–15781.

Sivolob, A.V. and Khrapunov, S.N. 1995. Translational positioning of nucleosomes on DNA: The role of sequence-dependent isotropic DNA bending stiffness. *J. Mol. Biol.* **247**: 918–931.

Smale, S.T. 2001. Core promoters: Active contributors to combinatorial gene regulation. *Genes & Dev.* **15**: 2503–2508.

Smale, S.T. and Kadonaga, J.T. 2003. The RNA polymerase II core promoter. *Annu. Rev. Biochem.* **72**: 449–479.

Solovyev, V.V. and Shahmuradov, I.A. 2003. PromH: Promoters identification using orthologous genomic sequences. *Nucleic Acids Res.* **31**: 3540–3545.

Solovyev, V., Kosarev, P., Seledsov, I., and Vorobyev, D. 2006. Automatic annotation of eukaryotic genes, pseudogenes and promoters. *Genome Biol.* **7**: 11–12.

Sonnenburg, S.O., Zien, A., and Rätsch, G. 2006. ARTS: Accurate recognition of transcription starts in human. *Bioinformatics* **22**: 472–480.

Struhl, K. 1989. Molecular mechanisms of transcriptional regulation in yeast. *Annu. Rev. Biochem.* **58**: 1051–1077.

Sugimoto, N., Nakano, S., Yoneyama, M., and Honda, K. 1996. Improved thermodynamic parameters and helix initiation factor to predict stability of DNA duplexes. *Nucleic Acids Res.* **24**: 4501–4505.

Suzuki, Y., Tsunoda, T., Sese, J., Taira, H., Mizushima-Sugano, J., Hata, H., Ota, T., Isogai, T., Tanaka, T., Nakamura, Y., et al. 2001. Identification and characterization of the potential promoter regions of 1031 kinds of human genes. *Genome Res.* **11**: 677–684.

Tabach, Y., Brosh, R., Buganim, Y., Reiner, A., Zuk, O., Yitzhaky, A., Koudritsky, M., Rotter, V., and Domany. E. 2007. Wide-scale analysis of human functional transcription factor binding reveals a strong bias towards the transcription start site. *PLoS ONE* **2**: e807. doi: 10.1371/journal.pone.0000807.

Trinklein, N.D., Aldred, S.F., Hartman, S.J., Schroeder, D.I., Otillar, R.P., and Myers, R.M. 2004. An abundance of bidirectional promoters in the human genome. *Genome Res.* **14**: 62–66.

Uren, P., Cameron-Jones, M., and Sale, A. 2006. Promoter prediction using physic-chemical properties of DNA. *Lect. Notes Comput. Sci.* **4216**: 21–31.

van Rijssbergen, C.J. 1979. *Information retrieval, 2nd edition*. Butterworths, London.

Wang, H. and Benham, C.J. 2006. Promoter prediction and annotation of microbial genomes based on DNA sequence and structural responses to superhelical stress. *BMC Bioinformatics* **7**: 248. doi: 10.1186/1471-2105-7-248.

Wang, G. and Zhang, W. 2006. A steganalysis-based approach to comprehensive identification and characterization of functional regulatory elements. *Genome Biol.* **7**: R49. doi: 10.1186/gb-2006-7-6-r49.

Wang, Z., Chen, Y., and Li, Y. 2004. A brief review of computational gene prediction methods. *Genom. Proteom. Bioinformatics* **2**: 216–221.

Wasserman, W.W., Palumbo, M., Thompson, W., Fickett, J.W., and Lawrence, C.E. 2000. Human-mouse genome comparisons to locate regulatory sites. *Nat. Genet.* **26**: 225–228.

Werner, T. 2003. The state of the art in mammalian promoter recognition. *Brief. Bioinform.* **4**: 22–30.

White, R.J. and Jackson, S.P. 1992. The TATA-binding protein: A central role in transcription by RNA polymerases I, II and III. *Trends Genet.* **8**: 284–288.

Xie, X., Wu, S., Lam, K.M., and Yan, H. 2006. PromoterExplorer: An effective promoter identification method based on the AdaBoost algorithm. *Bioinformatics* **22**: 2722–2728.

Zhang, M.Q. 2002. Computational prediction of eukaryotic protein-coding genes. *Nat. Rev. Genet.* **3**: 698–709.

Zhou, X., Ruan, J., Wang, G., and Zhang, W. 2007. Characterization and identification of microRNA core promoters in four model species. *PLoS Comput. Biol.* **3**: e37. doi: 10.1371/journal.pcbi.0003007.

Received August 3, 2007; accepted in revised form November 14, 2007.
<figure><img src="image_6.png" /></figure>

GENOME
RESEARCH

# Generic eukaryotic core promoter prediction using structural features of DNA

Thomas Abeel, Yvan Saeys, Eric Bonnet, et al.

*Genome Res.* 2008 18: 310-323 originally published online December 20, 2007
Access the most recent version at doi:10.1101/gr.6991408

**Supplemental Material** http://genome.cshlp.org/content/suppl/2007/12/20/gr.6991408.DC1

**References** This article cites 112 articles, 26 of which can be accessed free at: http://genome.cshlp.org/content/18/2/310.full.html#ref-list-1

License

**Email Alerting Service** Receive free email alerts when new articles cite this article - sign up in the box at the top right corner of the article or click here.

<figure><img src="image_7.png" /></figure>

To subscribe to *Genome Research* go to:
[https://genome.cshlp.org/subscriptions](https://genome.cshlp.org/subscriptions)
<figure><img src="image_7.png" /></figure>

GENOME
RESEARCH

# Generic eukaryotic core promoter prediction using structural features of DNA

Thomas Abeel, Yvan Saeys, Eric Bonnet, et al.

*Genome Res.* 2008 18: 310-323 originally published online December 20, 2007
Access the most recent version at doi:10.1101/gr.6991408

**Supplemental Material** http://genome.cshlp.org/content/suppl/2007/12/20/gr.6991408.DC1

**References** This article cites 112 articles, 26 of which can be accessed free at: http://genome.cshlp.org/content/18/2/310.full.html#ref-list-1

License

**Email Alerting Service** Receive free email alerts when new articles cite this article - sign up in the box at the top right corner of the article or click here.

<figure><img src="image_8.png" /></figure>

To subscribe to *Genome Research* go to:
[https://genome.cshlp.org/subscriptions](https://genome.cshlp.org/subscriptions)