BACKGROUND
coordinated regulation of gene expression is a fundamental process that depends on the binding of transcription factors to a gene's cis-regulatory sequences.
absolutely required for transcription initiation of metazoan protein-coding genes is the core promoter, the region of dna 35– <dig> bp upstream and downstream of the transcription start site   <cit> .
the core promoter contains sequence elements, referred to as "core promoter motifs," which interact with the basal transcription machinery, including rna polymerase ii and the tfiid complex  <cit> .
in recent years, it has become clear that the core promoter, rather than playing a passive role in the spatial and temporal regulation of gene expression, is an important active partner in these events  <cit> .
for instance, different promoter sequences are found preferentially associated with certain functional classes of genes, with genes expressed at particular developmental stages, and with genes expressed in the germ line versus the soma  <cit> .
various tissue-specific members of the tata box-binding protein  family, such as the tbp-related factors , bind preferentially to certain core promoters  <cit> .
there is also substantial evidence for preferred or specific promoter-enhancer interactions, whereby a distal cis-regulatory module  can stimulate activity from one promoter, but not another  <cit> .
a number of mechanisms have been demonstrated to restrict the activity of a crm to a particular promoter, including insulator elements  <cit> , insulator-bypass or promoter targeting elements  <cit> , short-range repression  <cit> , chromatin-mediated silencing  <cit> , and preferential interaction with promoters containing certain core promoter motifs  <cit> .
the relative prevalence of each of these mechanisms is unknown, as in most cases is a detailed understanding of how they function.
in particular, the molecular basis underlying core promoter preference has not been clearly defined.
the existence of crm-promoter specificity is all the more remarkable given that it is maintained despite the fact that there are sometimes other promoters closer to, or even interposed between, a crm and its target.
in fact, the latter may be a much more common scenario than typically credited, as it can occur not only with respect to the regulation of different genes, but also with respect to alternative promoters of the same gene.
in humans, it is estimated that upwards of 50% of all genes have at least one alternative promoter  <cit> , and there is growing evidence that alternative promoter usage plays important roles in both development and disease  <cit> .
it is unknown how frequently such alternative promoters are regulated by distinct crms, but the number could be large; kimura et al.
<cit>  suggest that over  <dig> sets of alternative promoters are regulated in a tissue-specific fashion.
except for the case of bidirectional promoters , few studies have focused specifically on promoters of neighboring genes or on alternative promoters, and little is known about the mechanisms that direct promoter usage choice.
baek et al.
<cit>  recently analyzed a subset of human promoters by dividing them into the four categories of cpg-island containing and non-containing single and alternative promoters, and observed differences in sequence properties, evolutionary conservation, biological roles, and degree of usage.
their data suggest that there may be differences among promoters depending on their relative position in the gene, with more upstream promoters being more highly expressed and more cpg-rich than the more downstream promoters.
interestingly, they found that the tata box and dpe core promoter motifs were more common in single than in alternative promoters.
however, a similar study by kimura et al.
<cit>  found little difference in the frequency of the tata box between the two groups, although they observed a large difference in the prevalence of cpg islands.
differences in the full set of promoters used and in how the promoters were grouped–the latter study did not look separately at the cpg-containing and non-containing promoters–may account for the discrepancies in the reported results.
a number of other sequence motifs, of unknown functional significance, were also seen to be differentially represented among the promoter classes  <cit> .
these studies suggest that there might be fundamental differences in the structure and function of single versus alternative promoters that could have broad implications for understanding how transcription is coordinated within the genome.
as a means to developing an estimate of how important the sequence of the promoter might be in dictating promoter usage choice and in mediating crm-promoter specificity, and as a prelude to experimental studies of the mechanisms of crm-promoter interactions, we undertook a global bioinformatics analysis of drosophila melanogaster promoters.
we found that there are marked differences in nucleotide composition and motif prevalence between single promoters and alternative promoters, and between the most 5' alternative promoters and more downstream alternative promoters.
we also observed that adjacent genes on the chromosome are more likely than expected to have promoters with a similar motif profile, and that this similarity in promoter configuration correlates with co-regulated gene expression.
our results suggest that promoter composition may play a larger-than-appreciated role in coordinating gene expression both between nearby genes and between multiple transcripts of the same gene.
RESULTS
in order to conduct a comprehensive genome-wide study of promoters, we looked at all protein coding genes annotated in the drosophila genome annotation release  <dig>   <cit> .
we considered all annotated tsss to be true tsss.
multiple tsss that were less than  <dig> bp apart were considered as sharing the same promoter .
all together, we obtained  <dig>  promoters from the genome, and separated these into three mutually exclusive classes : unique promoters , i.e., promoters for genes having only a single promoter; first alternative promoters , defined as the most 5' of a gene's multiple promoters ; and downstream alternative promoters , which are any alternative promoters 3' to the faps.
there were  <dig>   <dig>  and  <dig> promoters in ups, faps and daps respectively .
approximately 14% of genes had alternative promoters, with an average of  <dig>  promoters/gene  for those with more than one promoter.
although there are undoubtedly errors in the genome annotation with respect to tsss , we reasoned that these errors would be few relative to the total number of annotated genes and that therefore, while they might contribute noise to our analysis, they would not mask any clearly significant results.
however, to ensure the robustness of our results, we also generated two higher confidence sets of promoters: a "high quality" set of  <dig> promoters based on flybase transcript evidence annotations and a "cap-supported" set of  <dig> promoters in which all transcripts have their 5' ends accurately mapped based on cdna isolation using a 5' cap-dependent method .
we separated these smaller promoter sets into ups, faps, and daps, as we did with the full promoter set, and performed all analysis in parallel with the three sets .
in all cases, the results were essentially the same .
nucleotide distribution differs among the three promoter classes
as a starting point for comparisons among the three classes of promoters, we analyzed the distribution of the four nucleotides in each class by calculating the mean frequency of each nucleotide in a  <dig> bp sliding window across the tss region, ranging from - <dig> bp  to + <dig> bp .
promoters whose sequences overlapped that of an adjacent promoter of a different class were excluded from analysis.
similar to fitzgerald et al.
<cit> , who looked only at ups, we observed an at peak around - <dig> bp.
however, although this peak is clearly present for both ups and faps, it is absent for the daps, which have a markedly higher gc level between - <dig> bp and - <dig> bp.
although the nucleotide distributions of ups and faps are similar to one another from - <dig> bp to - <dig> bp, they begin to separate afterward: around the tsss, ups have more a but less c than the faps, and downstream of the tss, the ups are more gc rich than the faps .
the point at which the promoter classes begin to diverge is consistent with previous findings for human genes that the region beginning at approximately - <dig> bp comprises an extended promoter region with an important positive regulatory role  <cit> .
consistent with the trends we observed in the distribution plots, we found that the gc content differs significantly among all three classes of promoters with dap > fap > up   <dig> ) versus  <dig>  , kolmogorov-smirnov test p ≈ 0; dap versus fap:  <dig>   versus  <dig>  , kolmogorov-smirnov test p ≈  <dig> e-11; fap versus up:  <dig>   versus  <dig>  , kolmogorov-smirnov test p ≈  <dig> e-16).
note that the gc content was calculated after masking the coding regions in the promoters to avoid bias caused by the higher gc content of coding sequences  <cit> .
the nucleotide distributions indicate that all three promoter classes have distinct characteristics, which are most pronounced when comparing the daps to the others.
distribution of promoter motifs in the different promoter classes
the functional units of the core promoter are the "core promoter elements," the sequences that mediate binding of the general transcription factors to the promoter.
like most transcription factor binding sites, the sequences of these elements form a family of related subsequences known as a sequence "motif."
for brevity, and because in this study we do not explicitly evaluate binding but deal only with dna sequence, we refer here to both the sequence motifs themselves and to instances of these subsequences in promoter regions as "promoter motifs."
although the relationship between various core promoter motifs and tissue- and stage-specific gene regulation has been studied previously, as have possible associations between certain motifs and genes of particular functional classes  <cit> , there has been no systematic investigation of motif distributions in single versus alternative promoters, or among differently positioned alternative promoters.
therefore, we searched all  <dig>  promoters for the presence of the  <dig> promoter motifs identified by fitzgerald et al.
<cit>  in their analysis of drosophila ups that are overrepresented in core promoters or the extended promoter region up to - <dig> bp .
a full listing of the genomic coordinates of the mapped motifs and their sequences is provided .
for each of the  <dig> motifs, we looked at their relative distributions in ups, faps, and daps;  <dig> of them show a significant occurrence bias among the three promoter classes .
we see differences between unique and alternative promoters, between 5' promoters  and downstream promoters, and between first alternative and more downstream alternative promoters.
for example, tata/dmp <dig> is found in  <dig> % of ups,  <dig> -fold and  <dig> -fold higher than its occurrence in the two classes of alternative promoters, faps and daps, respectively.
ndm <dig> is significantly underrepresented in ups as compared to the two classes of alternative promoters.
dmv <dig>  dmv <dig>  dmv <dig>  ndm <dig> and e-box/ndm <dig> are underrepresented in the daps but do not differ significantly between the ups and faps.
inr/dmp <dig> is more prevalent in daps than in ups, whereas inr1/dmp <dig> and dpe1/dmp <dig> are less common in ups than in faps.
the presence of dmv <dig> and dre/ndm <dig> differs significantly among all three classes of promoters with fap > up > dap.
gaga/ndm1's occurrence is also significantly different among the three classes with dap > fap > up.
the occurrence biases we observed do not appear to correlate with the stage-specific motif usage noted by fitzgerald et al.
<cit> .
for instance, tata and gaga are both preferentially associated with adult-expressed genes, but while tata/dmp <dig> is overrepresented in ups, gaga/ndm <dig> is more prevalent in daps.
apercentage of the promoters in a class having the specific motif.
bcalculated using fisher's exact test.
holm's method was used for multiple hypothesis testing.
ns, not significant.
cpercentage of promoters without any known motifs.
as certain combinations of promoter motifs have been shown to preferentially co-occur–for example, tata and inr, and inr and dpe  <cit> –we also looked at the presence of motif combinations in the different promoter classes .
of the  <dig> significantly enriched motif combinations, six show significantly different occurrence among the three promoter classes .
overall, the pattern of motif combinations correlates with the occurrence of individual motifs.
for example, dmv <dig>  dre/ndm <dig> and e-box/ndm <dig> are significantly underrepresented in the daps, and the interactions between dmv <dig> and dre/ndm <dig> and between dre/ndm <dig> and e-box/ndm <dig> are also underrepresented in daps.
apercentage of the promoters in a class having the specific motif interaction.
bcalculated using fisher's exact test.
holm's method was used for multiple hypothesis testing.
ns, not significant.
one concern in this type of analysis is that the definition of, or method of searching for, the motifs may affect the results by altering the number and locations of the identified sequences.
we therefore repeated our analysis using two alternative methods for calculating the position requirement for the motifs, , and also by locating motif instances using patser  <cit>  and the position weight matrices defined by ohler et al.
<cit> .
all three methods gave qualitatively similar results with respect to the occurrence biases of individual motifs among the three different classes of promoters .
that is, although the different methods occasionally found different absolute numbers of motif instances , the motif occurrence biases among the respective promoter classes were the same.
human promoter classes also have distinct nucleotide compositions and motif preferences
we performed an analysis similar to that which we did for drosophila using  <dig> validated human promoters  <cit> .
nucleotide distribution plots show that both t and c levels vary among the three sets of promoters over the length of the promoter region, and a and g levels vary near the tsss .
the average gc content in human ups is significantly higher than that in faps and daps  versus  <dig>  , kolmogorov-smirnov test p =  <dig> e-06;  <dig>   versus  <dig>  , kolmogorov-smirnov test p = 0).
thus while the specific nucleotide compositions of the promoters are different between flies and humans, in both species there are pronounced differences in nucleotide preference among the three promoter classes.
to determine whether human promoters also display motif preferences among the three promoter classes, as we observed in drosophila, we mapped the locations of eight motifs previously identified as being overrepresented in human promoter sequences  <cit>  .
once again, we found significant differences in motif usage among the promoter classes .
in contrast to the fly promoter motifs, which collectively are distributed evenly among the promoter classes, most of the human promoter motifs have a higher occurrence frequency in ups, and the four motifs showing occurrence bias in different classes of promoters all differ significantly between unique and alternative promoters.
this compares well with the results of baek et al.
<cit>  showing a six-fold higher frequency of the tata box in cpg-poor ups.
our results suggest that not only does motif usage differ among human promoter classes, but that the motifs preferred in alternative promoters have not yet been identified in human.
consistent with this idea, we note that a much greater proportion of alternative promoters, as compared to unique promoters, lack any of the identified motifs we focused on in this study .
apercentage of the promoters in a class having the specific motif.
bcalculated using fisher's exact test.
holm's method was used for multiple hypothesis testing.
ns, not significant.
cpercentage of promoters without any known motifs.
promoter motif usage differs based on promoter position within a gene
the preceding analyses provide a general picture of promoter motif distributions among the different promoter classes, but do not tell us how motif profiles vary among the alternative promoters of the same gene.
however, this information could provide important insights into alternative promoter use and evolution.
for example, the presence of highly similar sets of motifs would suggest the possibility of extensive co-regulation of the alternative promoters.
conversely, very different motif compositions would be more consistent with differential regulation in which a distal regulatory element could interact with only one or a subset of the promoters.
landry et al.
<cit>  suggest a number of ways that alternative promoters may have evolved, one of which is through local sequence duplication.
this scenario would likely be reflected in higher-than-expected motif similarity, even given subsequent mutation away from the original sequence.
in order to evaluate the degree of relationship among the alternative promoters of individual genes, we compared which of the  <dig> promoter motifs are present between the first and second promoters of all drosophila genes with alternative promoters.
each pair of promoters was then assigned one of three levels of motif similarity: similar, intermediate, or different .
genes for which one or both of the two promoters contained no known motifs were excluded from the analysis.
in general, we find that the motif compositions between pairs of first and second alternative promoters are highly dissimilar.
however, in comparing alternative promoter pairs, we noticed that in some genes the two promoters are close enough in position that our motif mapping rules assign the same motif to both promoters.
without experimental data, we are unable to determine whether or not such motifs, which we refer to as "putatively shared motifs," are actually used by both promoters, or just one of the pair.
if we assume that putatively shared motifs are in fact used by both promoters, we find that although for the majority of the genes the first two promoters differ in motif composition, more promoter pairs are similar than we would expect at random .
if we consider that a putatively shared motif is used in only one of the two promoters, however, we find that although the degree of difference in promoter motif profiles between two alternative promoters of the same gene is still different from the random expectation, the statistical support for this conclusion is weaker and not significant in all three promoter datasets .
thus, while overall the motif composition between alternative promoters of the same gene is dissimilar, to what extent this dissimilarity is significantly less than what we would expect to see for a random pair of promoters depends on how the putatively shared motifs are actually used in the respective promoters.
the answer to this question must await detailed experimental investigation.
agenes with exactly two alternative promoters.
similar results were obtained for genes with multiple alternative promoters .
bpercentage of the alternative promoter pairs defined as similar, intermediate, and different .
cthe promoter set contains all promoters from drosophila genome annotation release  <dig> .
dthe promoter set contains only high quality promoters in fly genome; see methods.
ethe promoter set contains only cap-supported promoters in fly genome; see methods.
fmean of  <dig>  randomizations.
gempirical p-value calculated as the proportion of sampled randomizations where the percentage of the promoter pairs in each of the three similarity levels was no less than/no greater than the observed values
as mentioned above, promoter pairs that did not contain any of the  <dig> motifs in at least one promoter–approximately 84% of the potential promoter pairs –were omitted because we had no basis for determining how similar or different two such promoters were.
moreover, given the limited number of known promoter motifs, it is likely that for any of the promoter pairs there are additional but unidentified relevant sequence motifs.
in order to get around these twin difficulties, we compared how similar the full set of promoter pairs were to one another using d2z scores.
the d2z score is an alignment-free sequence comparison metric which compares k-mer word distributions between two sequences  <cit> .
the alternative promoter pairs of genes with exactly two promoters were more likely to have a higher d2z score than random expectation .
similar results were obtained for genes with multiple promoters .
this is not merely a result of the score distribution being driven by the subset of promoters with known motifs; breaking down the data into pairs for which both promoters contained known motifs, one promoter contained known motifs, or neither promoter contained known motifs revealed no differences in the score distribution for each subset .
these results indicate that for the full promoter data set, even without explicitly considering defined promoter motifs, alternative promoters of the same gene are more similar than random expectation, consistent with the results of our motif-based analysis of the smaller, motif-containing promoter subset.
similarities in motif composition between promoters of neighboring genes correlates with gene co-expression
not only is the analysis of alternative promoter pairs complicated by the possibility of shared motifs, but in most cases, we do not possess good data on promoter-specific gene expression.
on the other hand, questions about alternative promoter usage are paralleled by those with respect to neighboring genes: how do more distal regulatory elements select the proper promoter to activate when confronted by two or more promoters in relative proximity to one another?
unlike for alternative promoters, in the case of neighboring genes we have both unambiguous core promoter motif assignments and extensive gene expression data.
therefore, in addition to comparing promoter motif profiles among alternative promoters of the same gene, we looked at the motif composition between the promoters of neighboring genes.
so as not to confound our analysis with choices as to which alternative promoters to consider, we focused on neighboring genes with single promoters only , and compared the profiles of the  <dig> promoter motifs in the same way that we analyzed the motif composition of alternative promoter pairs.
although the motif profiles for the majority of neighboring genes are different, we found that the number of promoters with similar motifs is significantly higher than the random expectation .
there were a small number of neighboring unique promoters–most of which are bi-directional promoters–which, like we saw with some of the alternative promoters, could potentially share motifs.
however, removing the putatively shared motifs from one or the other of the neighboring promoters did not significantly change the result .
in other words, the promoters of neighboring genes are more closely related to one another than we would expect to see by chance alone.
apercentage of promoter pairs defined as similar, intermediate, and different .
bthe promoter set contains all promoters from drosophila genome annotation release  <dig> .
cthe promoter set contains only high quality promoters in fly genome; see methods.
dthe promoter set contains only cap-supported promoters in fly genome; see methods.
emean of  <dig>  randomizations.
fempirical p-value calculated as the proportion of sampled randomizations where the percentage of the promoter pairs in each of the three similarity levels was no less than/no greater than the observed values.
one possible reason for neighboring genes to have similar promoter organizations would be if the two genes were the result of a local sequence duplication event.
indeed, we observed higher sequence similarity in both transcribed regions and in the promoter regions from - <dig> bp to + <dig> bp  for neighboring unique genes that have similar motifs, compared to those with different motifs .
<dig> % of genes and  <dig> % of promoters have greater than 65% sequence alignment for promoters with similar motif profiles, versus  <dig> % and  <dig> %, respectively, for promoters with different profiles .
for neighboring unique promoters with similar motifs, but not for those with disparate motifs, promoter sequence similarity is highly correlated with gene sequence similarity , indicating that not just the known motifs but the sequences of the promoter region in general are highly related.
note that the majority of genes  with similar promoter motif profiles do not appear to be the result of gene duplication and that even then, promoter sequences have diverged much more rapidly than the gene sequences.
thus while potentially a contributing factor, gene duplication by itself cannot explain the unexpectedly high incidence of similar promoter motif profiles in neighboring genes.
we wondered whether the high number of neighboring genes with similar promoter make-ups might be indicative of their being coordinately regulated, e.g., by interaction of their promoters with the same enhancer.
we separated the neighboring gene pairs into two groups, those for which the promoters have a similar motif profile, and those whose promoters contain a different set of motifs.
we then compared the expression patterns of the genes belonging to the two different groups using gene expression data from  <dig> tissue-specific gene expression profiles contained in flyatlas  <cit> .
for each gene pair, we calculated the degree of co-expression over the  <dig> tissues and found that neighboring genes whose promoters have a similar motif composition are significantly more likely to be co-expressed than neighboring gene pairs with different motifs .
we also calculated how well correlated the level of expression was for genes of both groups.
for each tissue, the genes were ranked according to their mean expression level, and the pearson correlation coefficient of the expression ranks across all  <dig> tissues for any two genes was obtained.
neighboring genes whose promoters have similar motif profiles were more highly correlated in their expression level than neighboring genes having different motifs  =  <dig> ; one sided fisher's exact p =  <dig> e-04).
in order to be certain that these observed correlations between gene co-expression and promoter motif composition are not an artifact resulting from their high sequence similarity, we repeated our analysis using just those neighboring gene pairs with less than 65% gene sequence similarity.
in this group as well, gene pairs with similar promoter motif profiles are more likely to be co-expressed in same tissues and are more highly correlated in expression level than those with different motifs .
our data showing that approximately 18% of neighboring genes with unique promoters have similar promoter motif compositions, and that these genes tend to have correlated expression, is consistent with previous data on co-regulation of physically near-by drosophila genes  <cit> .
in particular, our results are reminiscent of the finding by spellman and rubin  <cit>  that roughly 20% of drosophila genes fall into clusters of adjacent genes with similar expression profiles.
we find that neighboring genes within one of these clusters are 63% more likely to have similar promoter motif profiles than those not within clusters .
thus at least some of the phenomenon described by spellman and rubin  <cit>  may be attributable to similarity in core promoter motifs among the genes in a co-expression cluster.
discussion
genome organization at the promoter level
we have systematically compared the nucleotide composition and promoter motif profiles of ups, faps, and daps throughout the drosophila genome.
our results demonstrate that the three types of promoters have distinct sequence and motif preferences.
intriguingly, although consistent with results from human promoters reported by baek et al.
<cit> , we observe clear differences based not only on whether a promoter is single or alternative, but also on the relative position of an alternative promoter among all of the alternative promoters.
eight of the  <dig> core and proximal promoter motifs we looked at occur differentially between faps and daps, and the daps have a distinct nucleotide profile.
thus, the genome appears to distinguish promoter types and positions.
our results suggest that each class may be subject to different modes of regulation or interact with differently constituted basal transcription complexes, and demonstrate that the genome has a complex organizational structure at the promoter level.
it is worth noting that there is no universally accepted method for accurately identifying genomic subsequences as being relevant instances of a particular defined sequence motif, and as such, a strictly in silico analysis will always be dependent on choices of method and parameters.
to guard against this, we used two different motif identification methods and three choices of range parameter for inclusion of a motif as part of the promoter.
we also used three different promoter sets , none of which is likely to be completely accurate and which will tend to variously include too many non-promoters  or too few real promoters .
reassuringly, we found that there were few substantive differences in our results when considering different methods and datasets; moreover, trends tended to be clearly preserved with the main differences appearing to stem from diminished statistical power when using the smaller data sets.
overall, we find our conclusions to be robust to choice of motif identification methods and search parameters, and promoter data sets.
our analysis of a more limited number of validated human promoters revealed that just as in drosophila, the different classes of promoters have distinct characteristics, and suggests that our observations point to a general metazoan organizational principle.
for example, similar to what we found in the fly genome, the tata motif is significantly overrepresented in human ups .
in fact, all four significant human motifs are overrepresented in ups, and there is a high percentage of alternative promoters that do not contain any known motifs.
thus promoter motifs in the human genome in general, and those specific to alternative promoters in particular, appear to be still awaiting identification.
one reason that these motifs have proven difficult to define may lie in the abundance of alternative promoters in the human genome; the differences in motif composition and gc content among ups, faps, and daps demonstrated by our results could contribute considerable noise to computational motif discovery attempts unless each promoter class is considered separately .
finding additional promoter motifs
more than half of the promoters in fly genome do not contain any of the  <dig> motifs we used for this study .
these promoters are problematic for purposes of promoter comparison, and some fraction of them may not in fact be true promoters, but rather might represent errors in the genome annotation that we used.
however, four lines of evidence suggest that the majority of these are genuine promoters.
one, we see a similar fraction of promoters lacking the known motifs when we use our more selective "high-quality" and "cap-supported" promoter sets.
two, analysis of the experimentally-verified drosophila promoters in the eukaryotic promoter database   <cit>  also reveals a high proportion without known motifs .
three, we see similar results overall for our analysis of human promoters, which relies exclusively on experimentally-verified promoter sequences.
four, our results using the d2z score demonstrate that the promoters without known motifs behave identically to those that have known motifs.
as a result, we believe that these represent true promoters and that a considerable number of promoter motifs remain to be identified.
indeed, a recent computational study has identified additional candidate drosophila upstream promoter motifs  <cit> , and it will be interesting to determine how these distribute with respect to the known motifs and to one another, and among the various promoter classes.
notably, all of the promoter motif discovery conducted to date has been performed on ups  <cit>  or on ups and faps jointly  <cit> .
as our data show that promoter motifs vary among the different single and alternative promoter classes, motifs specific for faps and daps may therefore be underrepresented among those that are known.
preliminary studies in our laboratory suggest that targeting motif discovery efforts to specific promoter subsets will be an effective strategy for identifying new motifs .
promoters and gene expression neighborhoods
the majority of adjacent promoters, either from neighboring genes or from alternative promoters of the same gene, have a highly dissimilar motif profile.
as promoter motifs have been implicated in helping to mediate specific enhancer-promoter interactions  <cit> , these differences are likely to represent one of the mechanisms used by the genome to prevent inappropriate gene activation by nearby crms.
nevertheless, neighboring genes are significantly more likely than expected to have highly similar promoters, with a strong correlation between motif similarity and strength of gene expression.
thus, a key role of the promoter may be in regulating levels of gene expression.
neighboring genes with similar promoters also show a concomitant increase in tissue co-expression, raising the possibility that they are either coordinately regulated by shared crms , or by individual crms that bind a similar complement of transcription factors  <cit> .
neighborhoods of co-regulated genes have been observed in many eukaryotes, including yeast, worm, fly, and human  <cit> .
a frequently proposed mechanism for this phenomenon is the presence of local chromatin domains; that is, that a local "open" chromatin configuration favorable for transcription–perhaps due to strong activation of one of the genes in the neighborhood–leads to spurious activation of other nearby genes, as opposed to regulated activation of each gene  <cit> .
in contrast, our data suggest that promoter sequences constitute a significant component of gene co-expression neighborhoods and point to a higher degree of genomic organization and regulation than seen for a chromatin-centric model.
note that the two mechanisms need not be mutually exclusive; for instance, the related promoters might all contribute to a strong local change in chromatin conformation that would affect even those neighboring genes with different promoter make-ups.
detailed experimental investigation along with careful mapping of chromatin modifications throughout the co-expression neighborhoods will be required to tease apart the various contributions of individual regulatory cis-trans interactions and epigenetic modifications.
however, our results suggest a much greater role for the promoter in mediating locally coordinated gene expression than heretofore appreciated.
CONCLUSIONS
our systematic investigation of drosophila promoters demonstrates that there are distinct sequence characteristics among unique, first alternative, and downstream alternative promoters, suggesting that different regulatory mechanisms may act preferentially on each class.
we also show that neighboring genes are unexpectedly likely to have similar promoter compositions, which correlates with an increased degree of gene coexpression and suggests a mechanism for the previously observed phenomenon of gene co-expression neighborhoods.
taken together, these data reveal a high degree of complex genome organization at the level of promoter sequences.

 coordinated regulation of gene expression is a fundamental process that depends on the binding of transcription factors to a gene's cis-regulatory sequences. absolutely required for transcription initiation of metazoan protein-coding genes is the core promoter, the region of dna 35– bp upstream and downstream of the transcription start site. the core promoter contains sequence elements, referred to as "core promoter motifs," which interact with the basal transcription machinery, including rna polymerase ii and the tfiid complex. in recent years, it has become clear that the core promoter, rather than playing a passive role in the spatial and temporal regulation of gene expression, is an important active partner in these events. for instance, different promoter sequences are found preferentially associated with certain functional classes of genes, with genes expressed at particular developmental stages, and with genes expressed in the germ line versus the soma. various tissue-specific members of the tata box-binding protein family, such as the tbp-related factors, bind preferentially to certain core promoters. there is also substantial evidence for preferred or specific promoter-enhancer interactions, whereby a distal cis-regulatory module can stimulate activity from one promoter, but not another. a number of mechanisms have been demonstrated to restrict the activity of a crm to a particular promoter, including insulator elements, insulator-bypass or promoter targeting elements, short-range repression, chromatin-mediated silencing, and preferential interaction with promoters containing certain core promoter motifs. the relative prevalence of each of these mechanisms is unknown, as in most cases is a detailed understanding of how they function. in particular, the molecular basis underlying core promoter preference has not been clearly defined. the existence of crm-promoter specificity is all the more remarkable given that it is maintained despite the fact that there are sometimes other promoters closer to, or even interposed between, a crm and its target. in fact, the latter may be a much more common scenario than typically credited, as it can occur not only with respect to the regulation of different genes, but also with respect to alternative promoters of the same gene. in humans, it is estimated that upwards of 50% of all genes have at least one alternative promoter, and there is growing evidence that alternative promoter usage plays important roles in both development and disease. it is unknown how frequently such alternative promoters are regulated by distinct crms, but the number could be large; kimura et al. suggest that over sets of alternative promoters are regulated in a tissue-specific fashion. except for the case of bidirectional promoters, few studies have focused specifically on promoters of neighboring genes or on alternative promoters, and little is known about the mechanisms that direct promoter usage choice. baek et al. recently analyzed a subset of human promoters by dividing them into the four categories of cpg-island containing and non-containing single and alternative promoters, and observed differences in sequence properties, evolutionary conservation, biological roles, and degree of usage. their data suggest that there may be differences among promoters depending on their relative position in the gene, with more upstream promoters being more highly expressed and more cpg-rich than the more downstream promoters. interestingly, they found that the tata box and dpe core promoter motifs were more common in single than in alternative promoters. however, a similar study by kimura et al. found little difference in the frequency of the tata box between the two groups, although they observed a large difference in the prevalence of cpg islands. differences in the full set of promoters used and in how the promoters were grouped–the latter study did not look separately at the cpg-containing and non-containing promoters–may account for the discrepancies in the reported results. a number of other sequence motifs, of unknown functional significance, were also seen to be differentially represented among the promoter classes. these studies suggest that there might be fundamental differences in the structure and function of single versus alternative promoters that could have broad implications for understanding how transcription is coordinated within the genome. as a means to developing an estimate of how important the sequence of the promoter might be in dictating promoter usage choice and in mediating crm-promoter specificity, and as a prelude to experimental studies of the mechanisms of crm-promoter interactions, we undertook a global bioinformatics analysis of drosophila melanogaster promoters. we found that there are marked differences in nucleotide composition and motif prevalence between single promoters and alternative promoters, and between the most 5' alternative promoters and more downstream alternative promoters. we also observed that adjacent genes on the chromosome are more likely than expected to have promoters with a similar motif profile, and that this similarity in promoter configuration correlates with co-regulated gene expression. our results suggest that promoter composition may play a larger-than-appreciated role in coordinating gene expression both between nearby genes and between multiple transcripts of the same gene.  in order to conduct a comprehensive genome-wide study of promoters, we looked at all protein coding genes annotated in the drosophila genome annotation release. we considered all annotated tsss to be true tsss. multiple tsss that were less than bp apart were considered as sharing the same promoter. all together, we obtained promoters from the genome, and separated these into three mutually exclusive classes: unique promoters, i.e., promoters for genes having only a single promoter; first alternative promoters, defined as the most 5' of a gene's multiple promoters; and downstream alternative promoters, which are any alternative promoters 3' to the faps. there were and promoters in ups, faps and daps respectively. approximately 14% of genes had alternative promoters, with an average of promoters/gene for those with more than one promoter. although there are undoubtedly errors in the genome annotation with respect to tsss, we reasoned that these errors would be few relative to the total number of annotated genes and that therefore, while they might contribute noise to our analysis, they would not mask any clearly significant results. however, to ensure the robustness of our results, we also generated two higher confidence sets of promoters: a "high quality" set of promoters based on flybase transcript evidence annotations and a "cap-supported" set of promoters in which all transcripts have their 5' ends accurately mapped based on cdna isolation using a 5' cap-dependent method. we separated these smaller promoter sets into ups, faps, and daps, as we did with the full promoter set, and performed all analysis in parallel with the three sets. in all cases, the results were essentially the same. nucleotide distribution differs among the three promoter classes as a starting point for comparisons among the three classes of promoters, we analyzed the distribution of the four nucleotides in each class by calculating the mean frequency of each nucleotide in a bp sliding window across the tss region, ranging from - bp to + bp. promoters whose sequences overlapped that of an adjacent promoter of a different class were excluded from analysis. similar to fitzgerald et al., who looked only at ups, we observed an at peak around - bp. however, although this peak is clearly present for both ups and faps, it is absent for the daps, which have a markedly higher gc level between - bp and - bp. although the nucleotide distributions of ups and faps are similar to one another from - bp to - bp, they begin to separate afterward: around the tsss, ups have more a but less c than the faps, and downstream of the tss, the ups are more gc rich than the faps. the point at which the promoter classes begin to diverge is consistent with previous findings for human genes that the region beginning at approximately - bp comprises an extended promoter region with an important positive regulatory role. consistent with the trends we observed in the distribution plots, we found that the gc content differs significantly among all three classes of promoters with dap > fap > up ) versus, kolmogorov-smirnov test p ≈ 0; dap versus fap: versus, kolmogorov-smirnov test p ≈ e-11; fap versus up: versus, kolmogorov-smirnov test p ≈ e-16). note that the gc content was calculated after masking the coding regions in the promoters to avoid bias caused by the higher gc content of coding sequences. the nucleotide distributions indicate that all three promoter classes have distinct characteristics, which are most pronounced when comparing the daps to the others. distribution of promoter motifs in the different promoter classes the functional units of the core promoter are the "core promoter elements," the sequences that mediate binding of the general transcription factors to the promoter. like most transcription factor binding sites, the sequences of these elements form a family of related subsequences known as a sequence "motif." for brevity, and because in this study we do not explicitly evaluate binding but deal only with dna sequence, we refer here to both the sequence motifs themselves and to instances of these subsequences in promoter regions as "promoter motifs." although the relationship between various core promoter motifs and tissue- and stage-specific gene regulation has been studied previously, as have possible associations between certain motifs and genes of particular functional classes, there has been no systematic investigation of motif distributions in single versus alternative promoters, or among differently positioned alternative promoters. therefore, we searched all promoters for the presence of the promoter motifs identified by fitzgerald et al. in their analysis of drosophila ups that are overrepresented in core promoters or the extended promoter region up to - bp. a full listing of the genomic coordinates of the mapped motifs and their sequences is provided. for each of the motifs, we looked at their relative distributions in ups, faps, and daps; of them show a significant occurrence bias among the three promoter classes. we see differences between unique and alternative promoters, between 5' promoters and downstream promoters, and between first alternative and more downstream alternative promoters. for example, tata/dmp is found in % of ups, -fold and -fold higher than its occurrence in the two classes of alternative promoters, faps and daps, respectively. ndm is significantly underrepresented in ups as compared to the two classes of alternative promoters. dmv dmv dmv ndm and e-box/ndm are underrepresented in the daps but do not differ significantly between the ups and faps. inr/dmp is more prevalent in daps than in ups, whereas inr1/dmp and dpe1/dmp are less common in ups than in faps. the presence of dmv and dre/ndm differs significantly among all three classes of promoters with fap > up > dap. gaga/ndm1's occurrence is also significantly different among the three classes with dap > fap > up. the occurrence biases we observed do not appear to correlate with the stage-specific motif usage noted by fitzgerald et al.. for instance, tata and gaga are both preferentially associated with adult-expressed genes, but while tata/dmp is overrepresented in ups, gaga/ndm is more prevalent in daps. apercentage of the promoters in a class having the specific motif. bcalculated using fisher's exact test. holm's method was used for multiple hypothesis testing. ns, not significant. cpercentage of promoters without any known motifs. as certain combinations of promoter motifs have been shown to preferentially co-occur–for example, tata and inr, and inr and dpe –we also looked at the presence of motif combinations in the different promoter classes. of the significantly enriched motif combinations, six show significantly different occurrence among the three promoter classes. overall, the pattern of motif combinations correlates with the occurrence of individual motifs. for example, dmv dre/ndm and e-box/ndm are significantly underrepresented in the daps, and the interactions between dmv and dre/ndm and between dre/ndm and e-box/ndm are also underrepresented in daps. apercentage of the promoters in a class having the specific motif interaction. bcalculated using fisher's exact test. holm's method was used for multiple hypothesis testing. ns, not significant. one concern in this type of analysis is that the definition of, or method of searching for, the motifs may affect the results by altering the number and locations of the identified sequences. we therefore repeated our analysis using two alternative methods for calculating the position requirement for the motifs,, and also by locating motif instances using patser and the position weight matrices defined by ohler et al.. all three methods gave qualitatively similar results with respect to the occurrence biases of individual motifs among the three different classes of promoters. that is, although the different methods occasionally found different absolute numbers of motif instances, the motif occurrence biases among the respective promoter classes were the same. human promoter classes also have distinct nucleotide compositions and motif preferences we performed an analysis similar to that which we did for drosophila using validated human promoters. nucleotide distribution plots show that both t and c levels vary among the three sets of promoters over the length of the promoter region, and a and g levels vary near the tsss. the average gc content in human ups is significantly higher than that in faps and daps versus, kolmogorov-smirnov test p = e-06; versus, kolmogorov-smirnov test p = 0). thus while the specific nucleotide compositions of the promoters are different between flies and humans, in both species there are pronounced differences in nucleotide preference among the three promoter classes. to determine whether human promoters also display motif preferences among the three promoter classes, as we observed in drosophila, we mapped the locations of eight motifs previously identified as being overrepresented in human promoter sequences. once again, we found significant differences in motif usage among the promoter classes. in contrast to the fly promoter motifs, which collectively are distributed evenly among the promoter classes, most of the human promoter motifs have a higher occurrence frequency in ups, and the four motifs showing occurrence bias in different classes of promoters all differ significantly between unique and alternative promoters. this compares well with the results of baek et al. showing a six-fold higher frequency of the tata box in cpg-poor ups. our results suggest that not only does motif usage differ among human promoter classes, but that the motifs preferred in alternative promoters have not yet been identified in human. consistent with this idea, we note that a much greater proportion of alternative promoters, as compared to unique promoters, lack any of the identified motifs we focused on in this study. apercentage of the promoters in a class having the specific motif. bcalculated using fisher's exact test. holm's method was used for multiple hypothesis testing. ns, not significant. cpercentage of promoters without any known motifs. promoter motif usage differs based on promoter position within a gene the preceding analyses provide a general picture of promoter motif distributions among the different promoter classes, but do not tell us how motif profiles vary among the alternative promoters of the same gene. however, this information could provide important insights into alternative promoter use and evolution. for example, the presence of highly similar sets of motifs would suggest the possibility of extensive co-regulation of the alternative promoters. conversely, very different motif compositions would be more consistent with differential regulation in which a distal regulatory element could interact with only one or a subset of the promoters. landry et al. suggest a number of ways that alternative promoters may have evolved, one of which is through local sequence duplication. this scenario would likely be reflected in higher-than-expected motif similarity, even given subsequent mutation away from the original sequence. in order to evaluate the degree of relationship among the alternative promoters of individual genes, we compared which of the promoter motifs are present between the first and second promoters of all drosophila genes with alternative promoters. each pair of promoters was then assigned one of three levels of motif similarity: similar, intermediate, or different. genes for which one or both of the two promoters contained no known motifs were excluded from the analysis. in general, we find that the motif compositions between pairs of first and second alternative promoters are highly dissimilar. however, in comparing alternative promoter pairs, we noticed that in some genes the two promoters are close enough in position that our motif mapping rules assign the same motif to both promoters. without experimental data, we are unable to determine whether or not such motifs, which we refer to as "putatively shared motifs," are actually used by both promoters, or just one of the pair. if we assume that putatively shared motifs are in fact used by both promoters, we find that although for the majority of the genes the first two promoters differ in motif composition, more promoter pairs are similar than we would expect at random. if we consider that a putatively shared motif is used in only one of the two promoters, however, we find that although the degree of difference in promoter motif profiles between two alternative promoters of the same gene is still different from the random expectation, the statistical support for this conclusion is weaker and not significant in all three promoter datasets. thus, while overall the motif composition between alternative promoters of the same gene is dissimilar, to what extent this dissimilarity is significantly less than what we would expect to see for a random pair of promoters depends on how the putatively shared motifs are actually used in the respective promoters. the answer to this question must await detailed experimental investigation. agenes with exactly two alternative promoters. similar results were obtained for genes with multiple alternative promoters. bpercentage of the alternative promoter pairs defined as similar, intermediate, and different. cthe promoter set contains all promoters from drosophila genome annotation release. dthe promoter set contains only high quality promoters in fly genome; see methods. ethe promoter set contains only cap-supported promoters in fly genome; see methods. fmean of randomizations. gempirical p-value calculated as the proportion of sampled randomizations where the percentage of the promoter pairs in each of the three similarity levels was no less than/no greater than the observed values as mentioned above, promoter pairs that did not contain any of the motifs in at least one promoter–approximately 84% of the potential promoter pairs –were omitted because we had no basis for determining how similar or different two such promoters were. moreover, given the limited number of known promoter motifs, it is likely that for any of the promoter pairs there are additional but unidentified relevant sequence motifs. in order to get around these twin difficulties, we compared how similar the full set of promoter pairs were to one another using d2z scores. the d2z score is an alignment-free sequence comparison metric which compares k-mer word distributions between two sequences. the alternative promoter pairs of genes with exactly two promoters were more likely to have a higher d2z score than random expectation. similar results were obtained for genes with multiple promoters. this is not merely a result of the score distribution being driven by the subset of promoters with known motifs; breaking down the data into pairs for which both promoters contained known motifs, one promoter contained known motifs, or neither promoter contained known motifs revealed no differences in the score distribution for each subset. these results indicate that for the full promoter data set, even without explicitly considering defined promoter motifs, alternative promoters of the same gene are more similar than random expectation, consistent with the results of our motif-based analysis of the smaller, motif-containing promoter subset. similarities in motif composition between promoters of neighboring genes correlates with gene co-expression not only is the analysis of alternative promoter pairs complicated by the possibility of shared motifs, but in most cases, we do not possess good data on promoter-specific gene expression. on the other hand, questions about alternative promoter usage are paralleled by those with respect to neighboring genes: how do more distal regulatory elements select the proper promoter to activate when confronted by two or more promoters in relative proximity to one another? unlike for alternative promoters, in the case of neighboring genes we have both unambiguous core promoter motif assignments and extensive gene expression data. therefore, in addition to comparing promoter motif profiles among alternative promoters of the same gene, we looked at the motif composition between the promoters of neighboring genes. so as not to confound our analysis with choices as to which alternative promoters to consider, we focused on neighboring genes with single promoters only, and compared the profiles of the promoter motifs in the same way that we analyzed the motif composition of alternative promoter pairs. although the motif profiles for the majority of neighboring genes are different, we found that the number of promoters with similar motifs is significantly higher than the random expectation. there were a small number of neighboring unique promoters–most of which are bi-directional promoters–which, like we saw with some of the alternative promoters, could potentially share motifs. however, removing the putatively shared motifs from one or the other of the neighboring promoters did not significantly change the result. in other words, the promoters of neighboring genes are more closely related to one another than we would expect to see by chance alone. apercentage of promoter pairs defined as similar, intermediate, and different. bthe promoter set contains all promoters from drosophila genome annotation release. cthe promoter set contains only high quality promoters in fly genome; see methods. dthe promoter set contains only cap-supported promoters in fly genome; see methods. emean of randomizations. fempirical p-value calculated as the proportion of sampled randomizations where the percentage of the promoter pairs in each of the three similarity levels was no less than/no greater than the observed values. one possible reason for neighboring genes to have similar promoter organizations would be if the two genes were the result of a local sequence duplication event. indeed, we observed higher sequence similarity in both transcribed regions and in the promoter regions from - bp to + bp for neighboring unique genes that have similar motifs, compared to those with different motifs. % of genes and % of promoters have greater than 65% sequence alignment for promoters with similar motif profiles, versus % and %, respectively, for promoters with different profiles. for neighboring unique promoters with similar motifs, but not for those with disparate motifs, promoter sequence similarity is highly correlated with gene sequence similarity, indicating that not just the known motifs but the sequences of the promoter region in general are highly related. note that the majority of genes with similar promoter motif profiles do not appear to be the result of gene duplication and that even then, promoter sequences have diverged much more rapidly than the gene sequences. thus while potentially a contributing factor, gene duplication by itself cannot explain the unexpectedly high incidence of similar promoter motif profiles in neighboring genes. we wondered whether the high number of neighboring genes with similar promoter make-ups might be indicative of their being coordinately regulated, e.g., by interaction of their promoters with the same enhancer. we separated the neighboring gene pairs into two groups, those for which the promoters have a similar motif profile, and those whose promoters contain a different set of motifs. we then compared the expression patterns of the genes belonging to the two different groups using gene expression data from tissue-specific gene expression profiles contained in flyatlas. for each gene pair, we calculated the degree of co-expression over the tissues and found that neighboring genes whose promoters have a similar motif composition are significantly more likely to be co-expressed than neighboring gene pairs with different motifs. we also calculated how well correlated the level of expression was for genes of both groups. for each tissue, the genes were ranked according to their mean expression level, and the pearson correlation coefficient of the expression ranks across all tissues for any two genes was obtained. neighboring genes whose promoters have similar motif profiles were more highly correlated in their expression level than neighboring genes having different motifs =; one sided fisher's exact p = e-04). in order to be certain that these observed correlations between gene co-expression and promoter motif composition are not an artifact resulting from their high sequence similarity, we repeated our analysis using just those neighboring gene pairs with less than 65% gene sequence similarity. in this group as well, gene pairs with similar promoter motif profiles are more likely to be co-expressed in same tissues and are more highly correlated in expression level than those with different motifs. our data showing that approximately 18% of neighboring genes with unique promoters have similar promoter motif compositions, and that these genes tend to have correlated expression, is consistent with previous data on co-regulation of physically near-by drosophila genes. in particular, our results are reminiscent of the finding by spellman and rubin that roughly 20% of drosophila genes fall into clusters of adjacent genes with similar expression profiles. we find that neighboring genes within one of these clusters are 63% more likely to have similar promoter motif profiles than those not within clusters. thus at least some of the phenomenon described by spellman and rubin may be attributable to similarity in core promoter motifs among the genes in a co-expression cluster. discussion genome organization at the promoter level we have systematically compared the nucleotide composition and promoter motif profiles of ups, faps, and daps throughout the drosophila genome. our results demonstrate that the three types of promoters have distinct sequence and motif preferences. intriguingly, although consistent with results from human promoters reported by baek et al., we observe clear differences based not only on whether a promoter is single or alternative, but also on the relative position of an alternative promoter among all of the alternative promoters. eight of the core and proximal promoter motifs we looked at occur differentially between faps and daps, and the daps have a distinct nucleotide profile. thus, the genome appears to distinguish promoter types and positions. our results suggest that each class may be subject to different modes of regulation or interact with differently constituted basal transcription complexes, and demonstrate that the genome has a complex organizational structure at the promoter level. it is worth noting that there is no universally accepted method for accurately identifying genomic subsequences as being relevant instances of a particular defined sequence motif, and as such, a strictly in silico analysis will always be dependent on choices of method and parameters. to guard against this, we used two different motif identification methods and three choices of range parameter for inclusion of a motif as part of the promoter. we also used three different promoter sets, none of which is likely to be completely accurate and which will tend to variously include too many non-promoters or too few real promoters. reassuringly, we found that there were few substantive differences in our results when considering different methods and datasets; moreover, trends tended to be clearly preserved with the main differences appearing to stem from diminished statistical power when using the smaller data sets. overall, we find our conclusions to be robust to choice of motif identification methods and search parameters, and promoter data sets. our analysis of a more limited number of validated human promoters revealed that just as in drosophila, the different classes of promoters have distinct characteristics, and suggests that our observations point to a general metazoan organizational principle. for example, similar to what we found in the fly genome, the tata motif is significantly overrepresented in human ups. in fact, all four significant human motifs are overrepresented in ups, and there is a high percentage of alternative promoters that do not contain any known motifs. thus promoter motifs in the human genome in general, and those specific to alternative promoters in particular, appear to be still awaiting identification. one reason that these motifs have proven difficult to define may lie in the abundance of alternative promoters in the human genome; the differences in motif composition and gc content among ups, faps, and daps demonstrated by our results could contribute considerable noise to computational motif discovery attempts unless each promoter class is considered separately. finding additional promoter motifs more than half of the promoters in fly genome do not contain any of the motifs we used for this study. these promoters are problematic for purposes of promoter comparison, and some fraction of them may not in fact be true promoters, but rather might represent errors in the genome annotation that we used. however, four lines of evidence suggest that the majority of these are genuine promoters. one, we see a similar fraction of promoters lacking the known motifs when we use our more selective "high-quality" and "cap-supported" promoter sets. two, analysis of the experimentally-verified drosophila promoters in the eukaryotic promoter database also reveals a high proportion without known motifs. three, we see similar results overall for our analysis of human promoters, which relies exclusively on experimentally-verified promoter sequences. four, our results using the d2z score demonstrate that the promoters without known motifs behave identically to those that have known motifs. as a result, we believe that these represent true promoters and that a considerable number of promoter motifs remain to be identified. indeed, a recent computational study has identified additional candidate drosophila upstream promoter motifs, and it will be interesting to determine how these distribute with respect to the known motifs and to one another, and among the various promoter classes. notably, all of the promoter motif discovery conducted to date has been performed on ups or on ups and faps jointly. as our data show that promoter motifs vary among the different single and alternative promoter classes, motifs specific for faps and daps may therefore be underrepresented among those that are known. preliminary studies in our laboratory suggest that targeting motif discovery efforts to specific promoter subsets will be an effective strategy for identifying new motifs. promoters and gene expression neighborhoods the majority of adjacent promoters, either from neighboring genes or from alternative promoters of the same gene, have a highly dissimilar motif profile. as promoter motifs have been implicated in helping to mediate specific enhancer-promoter interactions, these differences are likely to represent one of the mechanisms used by the genome to prevent inappropriate gene activation by nearby crms. nevertheless, neighboring genes are significantly more likely than expected to have highly similar promoters, with a strong correlation between motif similarity and strength of gene expression. thus, a key role of the promoter may be in regulating levels of gene expression. neighboring genes with similar promoters also show a concomitant increase in tissue co-expression, raising the possibility that they are either coordinately regulated by shared crms, or by individual crms that bind a similar complement of transcription factors. neighborhoods of co-regulated genes have been observed in many eukaryotes, including yeast, worm, fly, and human. a frequently proposed mechanism for this phenomenon is the presence of local chromatin domains; that is, that a local "open" chromatin configuration favorable for transcription–perhaps due to strong activation of one of the genes in the neighborhood–leads to spurious activation of other nearby genes, as opposed to regulated activation of each gene. in contrast, our data suggest that promoter sequences constitute a significant component of gene co-expression neighborhoods and point to a higher degree of genomic organization and regulation than seen for a chromatin-centric model. note that the two mechanisms need not be mutually exclusive; for instance, the related promoters might all contribute to a strong local change in chromatin conformation that would affect even those neighboring genes with different promoter make-ups. detailed experimental investigation along with careful mapping of chromatin modifications throughout the co-expression neighborhoods will be required to tease apart the various contributions of individual regulatory cis-trans interactions and epigenetic modifications. however, our results suggest a much greater role for the promoter in mediating locally coordinated gene expression than heretofore appreciated.  our systematic investigation of drosophila promoters demonstrates that there are distinct sequence characteristics among unique, first alternative, and downstream alternative promoters, suggesting that different regulatory mechanisms may act preferentially on each class. we also show that neighboring genes are unexpectedly likely to have similar promoter compositions, which correlates with an increased degree of gene coexpression and suggests a mechanism for the previously observed phenomenon of gene co-expression neighborhoods. taken together, these data reveal a high degree of complex genome organization at the level of promoter sequences. 
as a whole, our data reveal an unexpected complexity of genomic organization at the promoter level with respect to both alternative and neighboring promoters. our results suggest that different regulatory mechanisms may apply to each of the three promoter classes, and provide a mechanism for "gene expression neighborhoods," local clusters of co-expressed genes. we have undertaken a systematic investigation of drosophila promoters. we divided promoters into three groups: unique promoters, first alternative promoters, and downstream alternative promoters. although many genes have multiple promoters, little attention has been focused on how these relate to one another; nor has much study been directed at relationships between promoters of adjacent genes. we also investigated the promoters of neighboring genes and found that a greater than expected number of adjacent genes have similar sequence motif profiles, which may allow the genes to be regulated in a coordinated fashion. in addition to its role as an assembly site for the basal transcriptional apparatus, the promoter plays a key part in mediating temporal and spatial aspects of gene expression through differential binding of transcription factors and selective interaction with distal enhancers. consistent with this, there is a positive correlation between similar promoter motifs and related gene expression profiles for these genes.
 the promoter is a critical necessary transcriptional cis-regulatory element. in addition to its role as an assembly site for the basal transcriptional apparatus, the promoter plays a key part in mediating temporal and spatial aspects of gene expression through differential binding of transcription factors and selective interaction with distal enhancers. although many genes have multiple promoters, little attention has been focused on how these relate to one another; nor has much study been directed at relationships between promoters of adjacent genes.  we have undertaken a systematic investigation of drosophila promoters. we divided promoters into three groups: unique promoters, first alternative promoters, and downstream alternative promoters. we observed distinct nucleotide distribution and sequence motif preferences among these three classes. we also investigated the promoters of neighboring genes and found that a greater than expected number of adjacent genes have similar sequence motif profiles, which may allow the genes to be regulated in a coordinated fashion. consistent with this, there is a positive correlation between similar promoter motifs and related gene expression profiles for these genes.  our results suggest that different regulatory mechanisms may apply to each of the three promoter classes, and provide a mechanism for "gene expression neighborhoods," local clusters of co-expressed genes. as a whole, our data reveal an unexpected complexity of genomic organization at the promoter level with respect to both alternative and neighboring promoters. 