BACKGROUND
recent progress in high-throughput technologies has made it possible to collect a variety of "omics" data on transcripts and on the epigenetic behaviors of the histones that are often associated with these transcripts  <cit> .
cap analysis of gene expression  is a high-throughput method that enables large-scale identification of transcription start sites  of eukaryotic species.
this method measures gene expression levels simultaneously with tss identification by counting the sequenced 5' ends of full-length cdnas, termed cage tags  <cit> .
with the development of deep sequencing methods, more high-throughput, and high resolution "tag depth" measurements have become available   <cit> .
such recent whole-cell-level pictures of quantitative transcriptomes have revealed the complex transcriptional network of mammalian species  <cit> .
according to recent cage-based analyses of human tsss, the human "promotome" can be classified into two types of promoters by the degree of imprecision of their transcription initiation sites  <cit> .
one is the peak promoter, which initiates transcription strictly from a narrow genomic region , and the other is the broad promoter, which initiates transcription from wide-ranging positions   <cit> .
the peak promoters are suggested to be closely associated with the presence of the tata box  and with tissue-specific gene expression.
the broad promoters have been observed in the presence of cpg islands and drive relatively ubiquitous expression of the genes they control  <cit> .
the cpg-rich broad promoters are considered evolutionarily new and more likely to be controlled by epigenetic mechanisms, including dna methylation and sense-antisense regulation, than the peak promoters  <cit> .
these differences between broad and peak promoters raise questions of how these promoter types are associated with chromosomal structures and modifications and of how their difference confers cellular function.
in eukaryotic species, chromosomal dnas is packed into nucleosomes, each of which comprises approximately  <dig> base pairs wrapped around a histone protein octamer consisting of two copies of each of the four core histones, h2a, h2b, h <dig>  and h <dig>  <cit> .
two biologically important aspects of these histones are their positions and modifications, and it has been shown that these factors regulate transcription initiation  <cit> .
several methodologies have rapidly been developed for high-throughput identification of histone positions and modifications.
chip-chip identifies the histone-binding positions of genomic dna by using a combination of chromatin immunoprecipitation and tiling array  <cit> .
although chip-chip used to be a widely-used method, today, with the growing demand to develop high-throughput sequencing, the chip-seq method has been developed as a promising alternative to the tiling array-based approach in analyzing genome-wide nucleosome positioning  <cit> .
these methodologies have revealed several insights into the intertwining of gene expression with nucleosome position and histone modification.
for example, the degree of eviction of nucleosomes from the upstream regions of tsss is correlated with gene expression patterns in yeasts  <cit>  and humans  <cit> .
moreover, the methylated histone h <dig> at lysine  <dig>  and acetylated histone h <dig> at lysine  <dig> , located around tsss, are linked to gene activation  <cit> , whereas h3k27me <dig> and h3k9me <dig> are linked to gene repression  <cit> .
these modifications and related gene regulatory behaviors support the "histone code" hypothesis  <cit> , i.e.
that multiple histone modifications specify unique downstream functions.
however, the detailed mechanisms underlying transcriptional regulation by these histone behaviors are still obscure.
h3k9ac has recently been frequently observed around broad promoters  <cit> .
this implies that histone behavior is associated with promoter architecture, although this association has so far been found only in the case of h3k9ac, and the extent of such associations is unclear.
in this study, we systematically analyzed the relationships between histone behaviors and promoter architecture types by using information about  modified/unmodified histones;  their genomic positions relative to tsss;  their positional stabilities on the genome under two cellular conditions; and  gene expression.
the results showed that promoter architecture type and gene expression are tightly associated with the modification pattern and genomic positional stability of the histones forming nucleosomes.
they provide new insights into the epigenetic mechanisms of transcriptional regulation in terms of histone behavior.
RESULTS
promoter architecture and nucleosome positioning
we first focused on differences in nucleosome distribution around the two different types of transcription promoters .
we used human promoter positions for which information about the degree of transcription start imprecision had been obtained in a previous study  <cit> , as well as nucleosome positions defined as the genomic positions of histone h <dig> proteins in the resting condition in human cd4+ t-cells  <cit> .
we mapped them on human genomic sequences.
we then calculated the ratio of nucleosomes located at each genomic position relative to each peak and broad promoter.
we found that the nucleosome positions associated with broad promoters had markedly aligned and periodic patterns compared with those of peak promoters .
more strikingly, only in broad promoters, the first nucleosomes immediately downstream of the promoter were likely to be located in similar positions and those immediately upstream of the promoter were depleted .
this was contrary to our expectation; previous studies have reported that, in general, nucleosomes are distributed evenly around the promoter region  <cit> .
we had therefore expected that the nucleosome positions would be spread around the broad promoter and well aligned around the peak promoter, because tsss are widely spread in the broad promoter region but narrowly spread in the peak promoter region.
however, our results show that the broad promoter was specifically associated with a more aligned pattern of nucleosomes than the peak promoter.
h2a.z is a histone variant of h2a that is highly conserved among lower and higher eukaryotes.
enrichment of h2a.z around the promoter region has been also reported in yeast  <cit>  and humans  <cit> .
in terms of promoter architecture, we performed a similar analysis to the one of h <dig> shown in figure 1a of the positions of human nucleosomes harboring the histone variant h2a.z in human resting cd4+ t-cells  <cit> .
h2a.z was highly enriched around broad promoters but not peak promoters .
for example, the statistical significance of the enrichment was p <  <dig>  × 10- <dig>  for positions + <dig> to + <dig> with respect to the tss.
moreover, the distribution patterns of h2a.z were similar to those of h3; the positions of h2a.z were markedly aligned around broad promoters but not around peak promoters.
accessibility of transcription factor sp1
the two promoter architectures are associated with characteristic sequence contexts: the peak promoter is located close to a tata box and the broad promoter close to cpg islands  <cit> .
using the genomic positions of putative tata-box sites predicted by a position-specific weight matrix and the positions of cpg islands obtained from the ucsc genome browser database  <cit> , we confirmed that tata boxes were overrepresented in peak promoters and that broad promoters were highly associated with the presence of cpg islands .
it is possible that the aligned patterns of nucleosome positions around broad promoters are due to the accessibility of transcription factors to dna.
for instance, in the absence of the tata box, the ubiquitous transcription factor sp <dig> can recruit tata-binding proteins to initiate transcription  <cit> .
it has already been reported that consensus sp <dig> sites with high overall gc contents are overrepresented among broad promoters, and the positions of these sites for individual transcription units are less precise than those of tata boxes  <cit> .
consequently, we investigated the possibility that the nucleosomes around a broad promoter align in a more orderly fashion than those around the peak promoter because of the need to create a nucleosome-free region upstream of the tss to confer dna accessibility of transcription factor proteins.
we superimposed the distribution of putative sp <dig> sites  <cit>  around broad promoters onto that of the nucleosome positions , and we observed increased proportions of sp <dig> sites about - <dig> bp upstream of the broad promoter, where the nucleosome distribution was markedly depleted .
we conducted the same analysis for peak promoters.
the inverse relationship between sp <dig> site and nucleosome abundance around the broad promoter was much higher than that around the peak promoter, suggesting the plausibility of the dna accessibility model.
furthermore, we conducted a similar analysis for the binding sites of two other transcription factors, pu.
<dig> and maz, as a previous study  had analyzed the binding sites of these two factors in detail  <cit> .
the binding sites of both pu.
<dig> and maz were distributed on nucleosome-free regions around broad promoters, whereas no such trends were observed around peak promoters .
these results support the strong connection between the nucleosome-free region and the accessibility of transcription factors, which was specific to broad promoters.
positional stability of nucleosomes around broad promoters
if nucleosome positioning around broad promoters confers dna accessibility for the binding of transcription factors, then the nucleosome positions around broad promoters should be more stable throughout different cellular conditions than those around peak promoters, because broad promoters are usually associated with ubiquitously expressed gene   <cit>  and the genomic positions of transcription factor binding sites are fixed.
we analyzed the positional stability of nucleosomes located within positions + <dig> to + <dig> with respect to each promoter under "resting" and "activated" conditions of human cd4+ t-cells  <cit>  .
for each nucleosome position in the resting condition, we calculated the distance to the nearest nucleosome position in the activated condition in order to assess the positional stabilities of single nucleosomes under the two different cellular conditions.
the overall minimum distance was markedly shorter for broad promoter-associated nucleosomes than for peak promoter-associated ones .
in fact, the average absolute minimum distance in the case of the broad promoter  was significantly shorter than that for the peak promoter  .
these results demonstrated that the positions of nucleosomes around the broad promoters were more stable than those of nucleosomes around the peak promoters.
distribution of nucleosomes containing modified histones
it has been suggested that not only nucleosome position, but also nucleosomal histone modification, can regulate transcription  <cit> .
for instance, histone methylation is associated with either gene activation or repression, depending on the methylation site and state on the histone protein; in particular, methylation of histone h <dig>  in nucleosomes around the transcription promoter are well known to regulate gene expression  <cit> .
to investigate the differences in positional distribution of nucleosomes containing methylated histones around the two different types of promoter, we obtained nucleosome positions corresponding to each of three methylation types  in human cd4+ t cells from a previous study  <cit> , and we mapped these onto genomic sequences with the broad and peak promoter positions.
similar to the result for histone h <dig>  nucleosomes having h3k4me <dig>  - <dig>  and - <dig> were all highly enriched and well aligned around broad promoters, whereas they were depleted around peak promoter regions .
however, the alignment pattern of nucleosome positions differed depending on the type of methylation.
within the region downstream of the broad promoter, the first frequency peak of nucleosomes having h3k4me <dig> and  <dig> occurred in the + <dig> to + <dig> region , whereas those having h3k4me <dig> occurred in the + <dig> to + <dig> region .
for each methylation type, the difference in frequency of occurrence of nucleosomes with each type of modified histone in these regions between the peak and broad promoters was significant .
note that the values on the y-axes in figure  <dig> are not influenced by the absolute numbers of nucleosomes in each type of promoter, as they indicate the proportion of nucleosome-harboring tsss for each type of tss.
in addition to methylation, acetylation may control gene expression  <cit> .
we further analyzed nucleosome positioning corresponding to histone acetylation  in human cd4+ t cells and observed results similar to those for h3k4me <dig> .
for each of h <dig>  h2a.z, h3k4me <dig>  and h3k9ac, we estimated the abundance of nucleosomes associated with peak promoters relative to that of nucleosomes associated with broad promoters .
compared with nucleosomes carrying histone h <dig>  the relative abundances of nucleosomes carrying the modified histones or the histone variant were large, suggesting that the presence of histone modifications or a histone variant was highly associated with the broad promoter but not the peak promoter.
analysis of another genomic element that potentially influences histone behavior
methylation of cpg islands is tightly associated with the expression of downstream genes; a number of studies have therefore been conducted to analyze cpg islands at a genome-wide level  <cit> .
as described above, broad promoters are strongly associated with cpg islands .
therefore, it is possible that the enrichment of histone modifications and histone variants in the broad promoter region is derived merely from the effect of cpg islands and is independent of promoter architecture.
in fact, it has been shown that promoters with many cpg islands are more likely to harbor modified histones than promoters with fewer cpg islands  <cit> .
to address this issue, we analyzed the positions of nucleosomes having histone h <dig> and those having h3k4me <dig> around broad and peak promoters with and without cpg islands .
we found that, in the case where promoters were associated with cpg islands, nucleosomes with histone h3k4me <dig> were likely to be well aligned even around peak promoters.
however, broad promoter-associated nucleosomes were significantly more enriched than peak promoter-associated nucleosomes, especially in the region downstream of the promoter .
in contrast, when we focused only on promoters without cpg islands, nucleosomes having h3k4me <dig> were well aligned and enriched only around broad promoters ; the difference in the frequencies of downstream nucleosomes  potentially resulting from the difference in the alignment were significant .
broad promoters with cpg islands had an aligned pattern of nucleosomes carrying h <dig>  whereas no clear alignment was observed for peak promoters .
broad promoters without cpg islands still showed an aligned pattern of nucleosomes having h <dig> , whereas peak promoters had little alignment in the pattern .
these results show that the enrichment of nucleosomes having certain histones around a broad promoter is independent of the existence of cpg islands.
effect of histone modification on gene expression
to explore whether histone modification around the promoter affects gene expression, we analyzed the difference in expression levels of rnas transcribed from peak and broad promoters in terms of the existence of modified/unmodified histones in their surrounding regions.
we compared data sets of methylated/unmethylated histones and acetylated/unacetylated histones measured under resting conditions in human cd4+ t cells  <cit> .
gene expression data for resting cd4+ t cells were obtained from a previous study  <cit> ; we used only those genes for which the expression levels had been measured.
we classified promoters having at least one methylated/acetylated histone within the region from - <dig> to + <dig> as "promoters with methylated/acetylated histones" and all others as ones with unmethylated/unacetylated histones .
expression levels of genes 'associated with broad promoters that had methylated histones were significantly higher than those of genes associated with broad promoters with unmethylated histones .
conversely, the expression levels of genes associated with peak promoters having only unmethylated histones were as high as those of genes associated with peak promoters with methylated histones, and thus no significant difference was observed .
likewise, in the comparison between acetylated and unacetylated histones, the expression levels of genes associated with broad promoters that had acetylated histones were significantly higher than those of genes associated with broad promoters with no acetylated histones , but the expression levels of genes associated with peak promoters that had acetylated histones did not differ markedly from those of genes associated with peak promoters with only unacetylated histones .
these results suggest that the regulation of gene expression levels by histone modification is specific to broad promoter-associated genes.
discussion
we analyzed the global landscape of epigenetic relationships between histone modifications and transcription initiation by investigating genome-wide chip-seq data and deepcage data.
the results presented here show differences in the architecture of the broad and peak promoters that regulate gene expression.
especially, we revealed that the broad promoters were strongly associated with histones immediately downstream of the tss and they were frequently modified, presumably to regulate gene expression levels.
in previous studies, aligned patterns of nucleosome positions around tsss have been identified in yeasts and humans  <cit> .
however, we confirmed this alignment only for regions around tsss derived from broad promoters, not for those around tsss derived from peak promoters.
broad promoters have an aligned pattern of nucleosome positions around tsss and have large nucleosome-free regions immediately upstream of tsss.
studies in yeasts have validated the model of "open promoters," which have large, nucleosome-free regions immediately upstream of the tss and are often associated with tata-less promoters and poly -rich tracts, the sequences of which are unbendable and unstable for histone binding  <cit> .
the broad promoter characteristics that we found in humans are consistent with this model, because in humans the sequence patterns in cpg islands located upstream of tsss, in contrast to the yeast poly -rich tracts, have been shown to be unstable  <cit> .
our data indicate that the nucleosomes that are immediately downstream of tsss and associated with broad promoters are positioned in specific regions.
we suggest that broad promoters have these aligned patterns of nucleosome positions around tsss because the nucleosome position has a stronger impact on broad promotors than on peak promoters on the determination of tsss by transcription factors in the cell.
as an example of transcription factors that target broad promoters, we investigated the sp <dig> binding sites around tsss.
sp <dig> recognizes binding region of dna via its zinc finger domain whereas tbp recognizes tata box via its dna binding domain.
sp <dig> binding sites were enriched in the regions upstream of tsss corresponding to the nucleosome-free regions.
we observed similar tendencies for the binding sites of two transcription factors, pu.
<dig> and maz.
although biological experiments are necessary to investigate molecular mechanism behind this observation, we speculate that the nucleosome-free regions serve as "landing sites" for transcription factors, including sp <dig>  which have less precise binding motifs  than the tata box  <cit> .
in addition to histone h <dig>  we also analyzed the positions of the histone h2a variant h2a.z, which is enriched around tsss  <cit> , and we obtained similar results.
in contrast, peak promoters did not have aligned patterns of nucleosome positions.
one might suspect that the observation is due to high expression of genes associated with broad promoters, and low expression of those associated with peak promoters.
however even after we limited the analysis to broad and peak promoters both of which are associated with highly expressed genes, we still observed the preferences of h <dig> for broad promoters  compared to peak promoters .
although tsss for tata promoters are often fixed to single positions, our results suggest that such strictly controlled positions of tsss are not regulated by nucleosome position.
however, there is some evidence that the nucleosomes around tata promoters have regulatory roles in gene expression.
in yeasts, the tata promoter is one type of "covered promoter," and expression of the genes associated with such promoters is more likely to be inhibited by the presence of nucleosomes than expression of the genes associated with "open promoters," which are located in nucleosome-free regions  <cit> ; in covered promoters, nucleosomes often cover transcription factor binding sites to repress the expression of downstream genes.
it is also possible that, in humans, peak promoters associated with the tata box belong to one type of "covered promoter" where the expression of downstream genes is repressed by the presence of nucleosomes.
therefore, we speculate that transcription factor binding is controlled by nucleosome position in the case of peak promoters.
in our analysis of epigenetic control by histone modification, we uncovered an difference between broad and peak promoters.
h3k4me <dig>  - <dig>  and - <dig> and h3k9ac, which are associated with gene activation, were more highly enriched around tsss associated with broad promoters than around those associated with peak promoters.
thus broad promoters appeared to be under stronger epigenetic control than peak promoters.
we found a trend that further supported this hypothesis: the expression levels of genes associated with broad promoters that had modified histones had higher expression levels than genes associated with broad promoters without modified histones.
in contrast, peak promoters appeared to be under weaker epigenetic control, because far fewer of them harbored modified histones.
furthermore, there were no significant differences in the expression levels of genes associated with peak promoters that harbored or did not harbor modified histones.
it has been shown that promoters with many cpg islands are more likely to harbor modified histones than promoters with fewer cpg islands  <cit> .
however, even after we limited our analysis to promoters having cpg islands, number of broad promoters harboring h3k4me <dig> was still statistically higher than that of peak promoters.
even more remarkable differences were observed after we limited our analysis to promoters without cpg islands.
although these results may depend on the dataset of cpg islands we used, enrichment of h3k4me <dig> in downstream region  of broad promoters were still observed in the analysis using different dataset of cpg islands  <cit>  .
genes associated with broad promoters tend to be expressed ubiquitously, whereas those associated with peak promoters are likely to be expressed in specific tissues and may show low expression levels in most tissue types  <cit> .
therefore, if high levels of gene expression are directly associated with histone modifications around tsss, then we may observe spurious correlations between promoter type and histone modification.
in fact, h3k4me <dig> is known to upregulate the expression of downstream genes.
we therefore compared the distribution patterns of nucleosomes containing h3k4me <dig> around broad and peak promoters in cases where the downstream genes showed similar expression levels .
we found that the broad promoters also harbored more nucleosomes containing h3k4me <dig> in cases where the downstream genes showed similar expression levels ; the difference in the distributions of h3k4me <dig> around the broad and peak promoters was statistically significant , suggesting that promoter type was indeed associated with differences in epigenetic regulation by histone modifications.
peak promoters containing the tata box are regulated at their transcription initiation step, generally by the assembly of a pre-initiation complex with three additional components: the tata-associated factors, the so-called mediator complexes, and positive and negative cofactors.
we presume that peak promoters containing no tata box are regulated in a similar way.
this transcription system is widely used in various species, and our results suggest that it is unlikely to use epigenetic controls.
thus, broad and peak promoters have distinct systems to regulate gene expression.
throughout this work, we employed widely-accepted definition of peak promoters, i.e.
those which initiate transcription within the range of  <dig> bps.
changing this threshold to  <dig> bp did not have much effect on the distribution patterns of nucleosomes around broad and peak promoters as shown by pearson's correlation coefficients between histone distribution pattern around broad promoters  defined by >  <dig> bps threshold and that defined by >  <dig> bps threshold.
for h <dig> distribution patterns, correlation coefficients were  <dig>  and  <dig>  for broad and peak promoters, respectively.
for h3k4me <dig> distribution patterns, the correlation coefficients were  <dig>  for both broad and peak promoters.
these results suggest the robustness of the relationships between the imprecision of tss and patterns of histone distributions.
tata boxes are used in a wide range of organisms, including prokaryotes, and are thought to be part of an ancient transcriptional system.
in contrast, broad promoters are thought to be newly evolved  <cit>  and have incorporated histone modification systems.
our results showed that peak promoters, which are frequently associated with such ancient tata boxes, have not incorporated histone modification systems.
CONCLUSIONS
by using a computational approach, we discovered the general relationships between the two types of promoter architecture and histone behavior, including positioning and modification.
we first showed that the positions of histones around broad promoters were highly aligned and stable compared with those around peak promoters.
furthermore, we suggest that marked numbers of transcription initiations related to broad promoters are under the control of histone methylation and acetylation and are associated with gene expression level, whereas this is not the case with peak promoters.
these results indicate that the expression of genes associated with broad promoters, but not peak promoters, is highly associated with histone position and modification.
we believe that our study is a step in uncovering the general mechanisms underlying transcriptional systems and inferring how these systems have evolved.
this should eventually help us to understand the complexity of mammalian transcription.

 recent progress in high-throughput technologies has made it possible to collect a variety of "omics" data on transcripts and on the epigenetic behaviors of the histones that are often associated with these transcripts. cap analysis of gene expression is a high-throughput method that enables large-scale identification of transcription start sites of eukaryotic species. this method measures gene expression levels simultaneously with tss identification by counting the sequenced 5' ends of full-length cdnas, termed cage tags. with the development of deep sequencing methods, more high-throughput, and high resolution "tag depth" measurements have become available. such recent whole-cell-level pictures of quantitative transcriptomes have revealed the complex transcriptional network of mammalian species. according to recent cage-based analyses of human tsss, the human "promotome" can be classified into two types of promoters by the degree of imprecision of their transcription initiation sites. one is the peak promoter, which initiates transcription strictly from a narrow genomic region, and the other is the broad promoter, which initiates transcription from wide-ranging positions. the peak promoters are suggested to be closely associated with the presence of the tata box and with tissue-specific gene expression. the broad promoters have been observed in the presence of cpg islands and drive relatively ubiquitous expression of the genes they control. the cpg-rich broad promoters are considered evolutionarily new and more likely to be controlled by epigenetic mechanisms, including dna methylation and sense-antisense regulation, than the peak promoters. these differences between broad and peak promoters raise questions of how these promoter types are associated with chromosomal structures and modifications and of how their difference confers cellular function. in eukaryotic species, chromosomal dnas is packed into nucleosomes, each of which comprises approximately base pairs wrapped around a histone protein octamer consisting of two copies of each of the four core histones, h2a, h2b, h and h. two biologically important aspects of these histones are their positions and modifications, and it has been shown that these factors regulate transcription initiation. several methodologies have rapidly been developed for high-throughput identification of histone positions and modifications. chip-chip identifies the histone-binding positions of genomic dna by using a combination of chromatin immunoprecipitation and tiling array. although chip-chip used to be a widely-used method, today, with the growing demand to develop high-throughput sequencing, the chip-seq method has been developed as a promising alternative to the tiling array-based approach in analyzing genome-wide nucleosome positioning. these methodologies have revealed several insights into the intertwining of gene expression with nucleosome position and histone modification. for example, the degree of eviction of nucleosomes from the upstream regions of tsss is correlated with gene expression patterns in yeasts and humans. moreover, the methylated histone h at lysine and acetylated histone h at lysine, located around tsss, are linked to gene activation, whereas h3k27me and h3k9me are linked to gene repression. these modifications and related gene regulatory behaviors support the "histone code" hypothesis, i.e. that multiple histone modifications specify unique downstream functions. however, the detailed mechanisms underlying transcriptional regulation by these histone behaviors are still obscure. h3k9ac has recently been frequently observed around broad promoters. this implies that histone behavior is associated with promoter architecture, although this association has so far been found only in the case of h3k9ac, and the extent of such associations is unclear. in this study, we systematically analyzed the relationships between histone behaviors and promoter architecture types by using information about modified/unmodified histones; their genomic positions relative to tsss; their positional stabilities on the genome under two cellular conditions; and gene expression. the results showed that promoter architecture type and gene expression are tightly associated with the modification pattern and genomic positional stability of the histones forming nucleosomes. they provide new insights into the epigenetic mechanisms of transcriptional regulation in terms of histone behavior.  promoter architecture and nucleosome positioning we first focused on differences in nucleosome distribution around the two different types of transcription promoters. we used human promoter positions for which information about the degree of transcription start imprecision had been obtained in a previous study, as well as nucleosome positions defined as the genomic positions of histone h proteins in the resting condition in human cd4+ t-cells. we mapped them on human genomic sequences. we then calculated the ratio of nucleosomes located at each genomic position relative to each peak and broad promoter. we found that the nucleosome positions associated with broad promoters had markedly aligned and periodic patterns compared with those of peak promoters. more strikingly, only in broad promoters, the first nucleosomes immediately downstream of the promoter were likely to be located in similar positions and those immediately upstream of the promoter were depleted. this was contrary to our expectation; previous studies have reported that, in general, nucleosomes are distributed evenly around the promoter region. we had therefore expected that the nucleosome positions would be spread around the broad promoter and well aligned around the peak promoter, because tsss are widely spread in the broad promoter region but narrowly spread in the peak promoter region. however, our results show that the broad promoter was specifically associated with a more aligned pattern of nucleosomes than the peak promoter. h2a.z is a histone variant of h2a that is highly conserved among lower and higher eukaryotes. enrichment of h2a.z around the promoter region has been also reported in yeast and humans. in terms of promoter architecture, we performed a similar analysis to the one of h shown in figure 1a of the positions of human nucleosomes harboring the histone variant h2a.z in human resting cd4+ t-cells. h2a.z was highly enriched around broad promoters but not peak promoters. for example, the statistical significance of the enrichment was p < × 10- for positions + to + with respect to the tss. moreover, the distribution patterns of h2a.z were similar to those of h3; the positions of h2a.z were markedly aligned around broad promoters but not around peak promoters. accessibility of transcription factor sp1 the two promoter architectures are associated with characteristic sequence contexts: the peak promoter is located close to a tata box and the broad promoter close to cpg islands. using the genomic positions of putative tata-box sites predicted by a position-specific weight matrix and the positions of cpg islands obtained from the ucsc genome browser database, we confirmed that tata boxes were overrepresented in peak promoters and that broad promoters were highly associated with the presence of cpg islands. it is possible that the aligned patterns of nucleosome positions around broad promoters are due to the accessibility of transcription factors to dna. for instance, in the absence of the tata box, the ubiquitous transcription factor sp can recruit tata-binding proteins to initiate transcription. it has already been reported that consensus sp sites with high overall gc contents are overrepresented among broad promoters, and the positions of these sites for individual transcription units are less precise than those of tata boxes. consequently, we investigated the possibility that the nucleosomes around a broad promoter align in a more orderly fashion than those around the peak promoter because of the need to create a nucleosome-free region upstream of the tss to confer dna accessibility of transcription factor proteins. we superimposed the distribution of putative sp sites around broad promoters onto that of the nucleosome positions, and we observed increased proportions of sp sites about - bp upstream of the broad promoter, where the nucleosome distribution was markedly depleted. we conducted the same analysis for peak promoters. the inverse relationship between sp site and nucleosome abundance around the broad promoter was much higher than that around the peak promoter, suggesting the plausibility of the dna accessibility model. furthermore, we conducted a similar analysis for the binding sites of two other transcription factors, pu. and maz, as a previous study had analyzed the binding sites of these two factors in detail. the binding sites of both pu. and maz were distributed on nucleosome-free regions around broad promoters, whereas no such trends were observed around peak promoters. these results support the strong connection between the nucleosome-free region and the accessibility of transcription factors, which was specific to broad promoters. positional stability of nucleosomes around broad promoters if nucleosome positioning around broad promoters confers dna accessibility for the binding of transcription factors, then the nucleosome positions around broad promoters should be more stable throughout different cellular conditions than those around peak promoters, because broad promoters are usually associated with ubiquitously expressed gene and the genomic positions of transcription factor binding sites are fixed. we analyzed the positional stability of nucleosomes located within positions + to + with respect to each promoter under "resting" and "activated" conditions of human cd4+ t-cells. for each nucleosome position in the resting condition, we calculated the distance to the nearest nucleosome position in the activated condition in order to assess the positional stabilities of single nucleosomes under the two different cellular conditions. the overall minimum distance was markedly shorter for broad promoter-associated nucleosomes than for peak promoter-associated ones. in fact, the average absolute minimum distance in the case of the broad promoter was significantly shorter than that for the peak promoter. these results demonstrated that the positions of nucleosomes around the broad promoters were more stable than those of nucleosomes around the peak promoters. distribution of nucleosomes containing modified histones it has been suggested that not only nucleosome position, but also nucleosomal histone modification, can regulate transcription. for instance, histone methylation is associated with either gene activation or repression, depending on the methylation site and state on the histone protein; in particular, methylation of histone h in nucleosomes around the transcription promoter are well known to regulate gene expression. to investigate the differences in positional distribution of nucleosomes containing methylated histones around the two different types of promoter, we obtained nucleosome positions corresponding to each of three methylation types in human cd4+ t cells from a previous study, and we mapped these onto genomic sequences with the broad and peak promoter positions. similar to the result for histone h nucleosomes having h3k4me - and - were all highly enriched and well aligned around broad promoters, whereas they were depleted around peak promoter regions. however, the alignment pattern of nucleosome positions differed depending on the type of methylation. within the region downstream of the broad promoter, the first frequency peak of nucleosomes having h3k4me and occurred in the + to + region, whereas those having h3k4me occurred in the + to + region. for each methylation type, the difference in frequency of occurrence of nucleosomes with each type of modified histone in these regions between the peak and broad promoters was significant. note that the values on the y-axes in figure are not influenced by the absolute numbers of nucleosomes in each type of promoter, as they indicate the proportion of nucleosome-harboring tsss for each type of tss. in addition to methylation, acetylation may control gene expression. we further analyzed nucleosome positioning corresponding to histone acetylation in human cd4+ t cells and observed results similar to those for h3k4me. for each of h h2a.z, h3k4me and h3k9ac, we estimated the abundance of nucleosomes associated with peak promoters relative to that of nucleosomes associated with broad promoters. compared with nucleosomes carrying histone h the relative abundances of nucleosomes carrying the modified histones or the histone variant were large, suggesting that the presence of histone modifications or a histone variant was highly associated with the broad promoter but not the peak promoter. analysis of another genomic element that potentially influences histone behavior methylation of cpg islands is tightly associated with the expression of downstream genes; a number of studies have therefore been conducted to analyze cpg islands at a genome-wide level. as described above, broad promoters are strongly associated with cpg islands. therefore, it is possible that the enrichment of histone modifications and histone variants in the broad promoter region is derived merely from the effect of cpg islands and is independent of promoter architecture. in fact, it has been shown that promoters with many cpg islands are more likely to harbor modified histones than promoters with fewer cpg islands. to address this issue, we analyzed the positions of nucleosomes having histone h and those having h3k4me around broad and peak promoters with and without cpg islands. we found that, in the case where promoters were associated with cpg islands, nucleosomes with histone h3k4me were likely to be well aligned even around peak promoters. however, broad promoter-associated nucleosomes were significantly more enriched than peak promoter-associated nucleosomes, especially in the region downstream of the promoter. in contrast, when we focused only on promoters without cpg islands, nucleosomes having h3k4me were well aligned and enriched only around broad promoters; the difference in the frequencies of downstream nucleosomes potentially resulting from the difference in the alignment were significant. broad promoters with cpg islands had an aligned pattern of nucleosomes carrying h whereas no clear alignment was observed for peak promoters. broad promoters without cpg islands still showed an aligned pattern of nucleosomes having h, whereas peak promoters had little alignment in the pattern. these results show that the enrichment of nucleosomes having certain histones around a broad promoter is independent of the existence of cpg islands. effect of histone modification on gene expression to explore whether histone modification around the promoter affects gene expression, we analyzed the difference in expression levels of rnas transcribed from peak and broad promoters in terms of the existence of modified/unmodified histones in their surrounding regions. we compared data sets of methylated/unmethylated histones and acetylated/unacetylated histones measured under resting conditions in human cd4+ t cells. gene expression data for resting cd4+ t cells were obtained from a previous study; we used only those genes for which the expression levels had been measured. we classified promoters having at least one methylated/acetylated histone within the region from - to + as "promoters with methylated/acetylated histones" and all others as ones with unmethylated/unacetylated histones. expression levels of genes'associated with broad promoters that had methylated histones were significantly higher than those of genes associated with broad promoters with unmethylated histones. conversely, the expression levels of genes associated with peak promoters having only unmethylated histones were as high as those of genes associated with peak promoters with methylated histones, and thus no significant difference was observed. likewise, in the comparison between acetylated and unacetylated histones, the expression levels of genes associated with broad promoters that had acetylated histones were significantly higher than those of genes associated with broad promoters with no acetylated histones, but the expression levels of genes associated with peak promoters that had acetylated histones did not differ markedly from those of genes associated with peak promoters with only unacetylated histones. these results suggest that the regulation of gene expression levels by histone modification is specific to broad promoter-associated genes. discussion we analyzed the global landscape of epigenetic relationships between histone modifications and transcription initiation by investigating genome-wide chip-seq data and deepcage data. the results presented here show differences in the architecture of the broad and peak promoters that regulate gene expression. especially, we revealed that the broad promoters were strongly associated with histones immediately downstream of the tss and they were frequently modified, presumably to regulate gene expression levels. in previous studies, aligned patterns of nucleosome positions around tsss have been identified in yeasts and humans. however, we confirmed this alignment only for regions around tsss derived from broad promoters, not for those around tsss derived from peak promoters. broad promoters have an aligned pattern of nucleosome positions around tsss and have large nucleosome-free regions immediately upstream of tsss. studies in yeasts have validated the model of "open promoters," which have large, nucleosome-free regions immediately upstream of the tss and are often associated with tata-less promoters and poly -rich tracts, the sequences of which are unbendable and unstable for histone binding. the broad promoter characteristics that we found in humans are consistent with this model, because in humans the sequence patterns in cpg islands located upstream of tsss, in contrast to the yeast poly -rich tracts, have been shown to be unstable. our data indicate that the nucleosomes that are immediately downstream of tsss and associated with broad promoters are positioned in specific regions. we suggest that broad promoters have these aligned patterns of nucleosome positions around tsss because the nucleosome position has a stronger impact on broad promotors than on peak promoters on the determination of tsss by transcription factors in the cell. as an example of transcription factors that target broad promoters, we investigated the sp binding sites around tsss. sp recognizes binding region of dna via its zinc finger domain whereas tbp recognizes tata box via its dna binding domain. sp binding sites were enriched in the regions upstream of tsss corresponding to the nucleosome-free regions. we observed similar tendencies for the binding sites of two transcription factors, pu. and maz. although biological experiments are necessary to investigate molecular mechanism behind this observation, we speculate that the nucleosome-free regions serve as "landing sites" for transcription factors, including sp which have less precise binding motifs than the tata box. in addition to histone h we also analyzed the positions of the histone h2a variant h2a.z, which is enriched around tsss, and we obtained similar results. in contrast, peak promoters did not have aligned patterns of nucleosome positions. one might suspect that the observation is due to high expression of genes associated with broad promoters, and low expression of those associated with peak promoters. however even after we limited the analysis to broad and peak promoters both of which are associated with highly expressed genes, we still observed the preferences of h for broad promoters compared to peak promoters. although tsss for tata promoters are often fixed to single positions, our results suggest that such strictly controlled positions of tsss are not regulated by nucleosome position. however, there is some evidence that the nucleosomes around tata promoters have regulatory roles in gene expression. in yeasts, the tata promoter is one type of "covered promoter," and expression of the genes associated with such promoters is more likely to be inhibited by the presence of nucleosomes than expression of the genes associated with "open promoters," which are located in nucleosome-free regions; in covered promoters, nucleosomes often cover transcription factor binding sites to repress the expression of downstream genes. it is also possible that, in humans, peak promoters associated with the tata box belong to one type of "covered promoter" where the expression of downstream genes is repressed by the presence of nucleosomes. therefore, we speculate that transcription factor binding is controlled by nucleosome position in the case of peak promoters. in our analysis of epigenetic control by histone modification, we uncovered an difference between broad and peak promoters. h3k4me - and - and h3k9ac, which are associated with gene activation, were more highly enriched around tsss associated with broad promoters than around those associated with peak promoters. thus broad promoters appeared to be under stronger epigenetic control than peak promoters. we found a trend that further supported this hypothesis: the expression levels of genes associated with broad promoters that had modified histones had higher expression levels than genes associated with broad promoters without modified histones. in contrast, peak promoters appeared to be under weaker epigenetic control, because far fewer of them harbored modified histones. furthermore, there were no significant differences in the expression levels of genes associated with peak promoters that harbored or did not harbor modified histones. it has been shown that promoters with many cpg islands are more likely to harbor modified histones than promoters with fewer cpg islands. however, even after we limited our analysis to promoters having cpg islands, number of broad promoters harboring h3k4me was still statistically higher than that of peak promoters. even more remarkable differences were observed after we limited our analysis to promoters without cpg islands. although these results may depend on the dataset of cpg islands we used, enrichment of h3k4me in downstream region of broad promoters were still observed in the analysis using different dataset of cpg islands. genes associated with broad promoters tend to be expressed ubiquitously, whereas those associated with peak promoters are likely to be expressed in specific tissues and may show low expression levels in most tissue types. therefore, if high levels of gene expression are directly associated with histone modifications around tsss, then we may observe spurious correlations between promoter type and histone modification. in fact, h3k4me is known to upregulate the expression of downstream genes. we therefore compared the distribution patterns of nucleosomes containing h3k4me around broad and peak promoters in cases where the downstream genes showed similar expression levels. we found that the broad promoters also harbored more nucleosomes containing h3k4me in cases where the downstream genes showed similar expression levels; the difference in the distributions of h3k4me around the broad and peak promoters was statistically significant, suggesting that promoter type was indeed associated with differences in epigenetic regulation by histone modifications. peak promoters containing the tata box are regulated at their transcription initiation step, generally by the assembly of a pre-initiation complex with three additional components: the tata-associated factors, the so-called mediator complexes, and positive and negative cofactors. we presume that peak promoters containing no tata box are regulated in a similar way. this transcription system is widely used in various species, and our results suggest that it is unlikely to use epigenetic controls. thus, broad and peak promoters have distinct systems to regulate gene expression. throughout this work, we employed widely-accepted definition of peak promoters, i.e. those which initiate transcription within the range of bps. changing this threshold to bp did not have much effect on the distribution patterns of nucleosomes around broad and peak promoters as shown by pearson's correlation coefficients between histone distribution pattern around broad promoters defined by > bps threshold and that defined by > bps threshold. for h distribution patterns, correlation coefficients were and for broad and peak promoters, respectively. for h3k4me distribution patterns, the correlation coefficients were for both broad and peak promoters. these results suggest the robustness of the relationships between the imprecision of tss and patterns of histone distributions. tata boxes are used in a wide range of organisms, including prokaryotes, and are thought to be part of an ancient transcriptional system. in contrast, broad promoters are thought to be newly evolved and have incorporated histone modification systems. our results showed that peak promoters, which are frequently associated with such ancient tata boxes, have not incorporated histone modification systems.  by using a computational approach, we discovered the general relationships between the two types of promoter architecture and histone behavior, including positioning and modification. we first showed that the positions of histones around broad promoters were highly aligned and stable compared with those around peak promoters. furthermore, we suggest that marked numbers of transcription initiations related to broad promoters are under the control of histone methylation and acetylation and are associated with gene expression level, whereas this is not the case with peak promoters. these results indicate that the expression of genes associated with broad promoters, but not peak promoters, is highly associated with histone position and modification. we believe that our study is a step in uncovering the general mechanisms underlying transcriptional systems and inferring how these systems have evolved. this should eventually help us to understand the complexity of mammalian transcription. 
transcription promoters are fundamental genomic cis-elements controlling gene expression. here, we performed a systematic analysis of transcription promoters and gene expression, as well as of epigenetic histone behaviors, including genomic position, stability within the chromatin, and several modifications. eukaryotic transcription initiation is suggested to be associated with the genomic positions and modifications of nucleosomes. they can be classified into two types by the degree of imprecision of their transcription start sites: peak promoters, which initiate transcription from a narrow genomic region; and broad promoters, which initiate transcription from a wide-ranging region. for instance, it has been recently shown that histone with h3k acetylation is more likely to be distributed around broad promoters rather than peak promoters; it can thus be inferred that there is an association between histone h3k and promoter architecture. these results shed light on how epigenetic regulatory networks of histone modifications are associated with promoter architecture. we found that, in humans, broad promoters, but not peak promoters, generally had significant associations with nucleosome positioning and modification. specifically, around broad promoters histones were highly distributed and aligned in an orderly fashion. this feature was more evident with histones that were methylated or acetylated; moreover, the nucleosome positions around the broad promoters were more stable than those around the peak ones. more strikingly, the overall expression levels of genes associated with broad promoters with modified histones were significantly higher than the levels of genes associated with broad promoters with unmodified histones.
 transcription promoters are fundamental genomic cis-elements controlling gene expression. they can be classified into two types by the degree of imprecision of their transcription start sites: peak promoters, which initiate transcription from a narrow genomic region; and broad promoters, which initiate transcription from a wide-ranging region. eukaryotic transcription initiation is suggested to be associated with the genomic positions and modifications of nucleosomes. for instance, it has been recently shown that histone with h3k acetylation is more likely to be distributed around broad promoters rather than peak promoters; it can thus be inferred that there is an association between histone h3k and promoter architecture.  here, we performed a systematic analysis of transcription promoters and gene expression, as well as of epigenetic histone behaviors, including genomic position, stability within the chromatin, and several modifications. we found that, in humans, broad promoters, but not peak promoters, generally had significant associations with nucleosome positioning and modification. specifically, around broad promoters histones were highly distributed and aligned in an orderly fashion. this feature was more evident with histones that were methylated or acetylated; moreover, the nucleosome positions around the broad promoters were more stable than those around the peak ones. more strikingly, the overall expression levels of genes associated with broad promoters with modified histones were significantly higher than the levels of genes associated with broad promoters with unmodified histones.  these results shed light on how epigenetic regulatory networks of histone modifications are associated with promoter architecture. 