- Research article
- Open Access
Nucleosome-coupled expression differences in closely-related species
© Guan et al; licensee BioMed Central Ltd. 2011
- Received: 8 November 2010
- Accepted: 26 September 2011
- Published: 26 September 2011
Genome-wide nucleosome occupancy is negatively related to the average level of transcription factor motif binding based on studies in yeast and several other model organisms. The degree to which nucleosome-motif interactions relate to phenotypic changes across species is, however, unknown.
We address this challenge by generating nucleosome positioning and cell cycle expression data for Saccharomyces bayanus and show that differences in nucleosome occupancy reflect cell cycle expression divergence between two yeast species, S. bayanus and S. cerevisiae. Specifically, genes with nucleosome-depleted MBP1 motifs upstream of their coding sequence show periodic expression during the cell cycle, whereas genes with nucleosome-shielded motifs do not. In addition, conserved cell cycle regulatory motifs across these two species are more nucleosome-depleted compared to those that are not conserved, suggesting that the degree of conservation of regulatory sites varies, and is reflected by nucleosome occupancy patterns. Finally, many changes in cell cycle gene expression patterns across species can be correlated to changes in nucleosome occupancy on motifs (rather than to the presence or absence of motifs).
Our observations suggest that alteration of nucleosome occupancy is a previously uncharacterized feature related to the divergence of cell cycle expression between species.
- Nucleosome Position
- Nucleosome Occupancy
- Gene Expression Divergence
- Upstream Regulatory Region
- Cell Cycle Regulate Gene
An organism's DNA contains numerous regulatory sequences that are used to modulate gene expression; yet DNA sequence alone does not explain why some regulatory sequences are functional while others are not. Because most genomic DNA (80% on average) is tightly packaged into nucleosomes , alternating nucleosome occupancy has been proposed as an important strategy to regulate gene expression since its initial discovery [2, 3]. Indeed, higher expression levels are commonly associated with nucleosome depletion at promoters and other genomic locations, e.g. rDNA [1, 4–6]. It has also been demonstrated that nucleosome occupancy affects the accessibility of DNA sequence motifs to transcriptional regulators; as a consequence different DNA sequences can display different nucleosome occupancy levels [1, 4, 7]. Further, motifs recognized and bound by active transcription factors are more likely to be nucleosome-depleted than those recognized by inactive ones [1, 8–13]. Differential occupancy on many motifs has been observed in certain environmental conditions [14, 15] and following environmental stresses . However, it remains controversial whether changes of nucleosome occupancy  or their initial positioning  determines levels of gene expression.
Most previous studies have focused on measurements of average transcription levels and average nucleosome occupancy over regulatory regions. The one-to-one connection between the occupancy of individual motifs and the resulting effect on gene expression has been tested only for a small number of genes. A recent study demonstrated that nucleosome depletion at two cell cycle-regulated promoters, CLN2pr and HOpr, ensures periodic expression pattern of genes involved in cell-cycle progression . These experiments clearly linked a specific expression pattern (cell-cycle periodicity) to nucleosome occupancy. The generality of this phenomenon for genes containing cell cycle regulating motifs remains to be tested through genome-wide experiments.
An average correlation between expression level and nucleosome occupancy at promoters across species has been reported , but it is not, however, clear how motif-specific nucleosome occupancy patterns affect the expression of individual genes across different species. To address this question, we sought an analysis approach that transcends the average expression level and targets the response at a specific class of motifs under specific conditions. Although predictions of nucleosome occupancy often assume that nucleosome positions are identical on conserved DNA sequences , experimental data is needed to test this assumption to better understand how nucleosome occupancy on motifs relates to phenotypic evolution. Such comparison across species can provide insight that augments ongoing efforts to define the relative contributions of cis and trans acting factors in phenotype divergence.
In this study, we determined the genome-wide nucleosome positions in the yeast S. bayanus, and compared these findings to patterns of gene expression during the cell cycle of S. cerevisiae and S. bayanus, two closely related sensu stricto yeast species. We show that changes in nucleosome occupancy on motifs are correlated with phenotypic divergence between species. In particular, our results show that nucleosomes provide a conspicuous genome-wide signature of MBP1 cell-cycle motif recognition in these two yeasts and this signature distinguishes which motifs result in periodic, cyclic expression patterns of the downstream genes. Although averaged expression level has previously been negatively linked to nucleosome occupancy at promoters [1, 4–6], our data provide a high-resolution, genome-wide demonstration of how the interplay between nucleosome occupancy and motif content is related to a specific expression pattern (i.e. in the cell cycle) of individual genes. Conserved transcription factor binding sites are more likely to be nucleosome-depleted, suggesting that patterns of nucleosome occupancy may reflect conservation of regulatory circuits across species. Finally, our cross-species comparison of transcription factor binding sites and nucleosome occupancy patterns reveals that changes in nucleosome-motif interactions are correlated to expression divergence, i.e. despite their conserved presence, motifs that become nucleosome-occupied during evolution no longer regulate downstream gene expression.
Global nucleosome occupancy in S. bayanus
We compared the closely related Saccharomyces species S. bayanus and S. cerevisiae to investigate how the interplay between nucleosomes and transcription factor binding motifs may affect downstream gene expression divergence. The two species are separated by 20 million years of evolution, a practical distance allowing us to investigate expression alterations between related genomes while still allowing for relatively unambiguous ortholog assignments. To map the genome-wide profile of nucleosome occupancy of S. bayanus, we used high-throughput, short read sequencing to detect nucleosome positioning as described in [1, 15] (data deposited in the short read archive GSE24356). Applying a wavelet model , we identified 47,777 well-defined nucleosomes (Additional file 1, Test S1), consistent with published work , with 97.95% of the S. bayanus genome assembly sequence covered (Additional file 2, Figure S1, Additional file 2, Figure S2).
As a quality control, we examined the S. bayanus genome for several features known to be present. Previous observations in S. cerevisiae and S. bayanus showed stereotypical nucleosome positioning relative to the start codon [1, 18]. Consistent with previous observations of nucleosome depletion at start codons (ATG) in S. cerevisiae, in S. bayanus we observe nucleosome depletion centered at ~150 bp upstream of the start codon (Additional file 2, Figure S2B) and at ~90 bp downstream of the stop codon (Additional file 2, Figure S2C). In addition, we found nucleosome occupancy in S. bayanus peaks immediately downstream of the start codon, as observed in previous publications [18, 22]. This confirms the technical quality of the sequencing data and provides an independent demonstration that the overall nucleosome occupancy pattern is conserved between S. bayanus and S. cerevisiae.
Nucleosome depletion on cell cycle regulatory motifs predicts periodic expression for downstream genes in S. bayanus and S. cerevisiae
Previous studies have reported that, on average, nucleosomes are depleted over transcription factor binding motifs [1, 8], leading to the suggestion that the active transcription factors are correlated to depleted nucleosomes on their target sites [1, 8–12]. We wanted to test if this correlation holds true at the level of individual genes, i.e. if, at a given motif upstream of several genes that vary in their gene expression levels, the level of nucleosome occupancy is correlated to these differences in gene expression. To test this possibility, we analyzed cell cycle gene expression data from S. bayanus cultures synchronized with alpha-factor. Genes were ranked by their cell cycle expression periodicity as determined by Fourier transform (Additional file 2, Figure S3). The top motifs enriched in cell cycle-regulated genes were determined by a mutual information based algorithm, FIRE , with the consensus sequences [AGT]ACGCG[AT][ACG]A. This motif in S. bayanus maps to the S. cerevisiae cell cycle regulatory motifs for the transcription factor Swi4p (e = 2.01e-09) based on the motif comparative tool STAMP  and the AlignACE database . The SBF (Swi4p-Swi6p) complex, in concert with MBF (Mbp1-Swi6p), is the principal transcriptional regulator of the yeast cell cycle .
To study the generality of this motif-nucleosome interaction on individual gene expression, we explored whether the motif-nucleosome interplay holds true in the related species, S. cerevisiae. Using cell cycle expression data from S. cerevisiae cells synchronized with alpha factor under similar conditions to our S. bayanus experiment , we identified genes displaying cell cycle periodicity (Additional file 2, Figure S3). The most significant motif we identified, with the consensus [AGT][AT]CGCGT[CT][AGT], corresponded to the MBP1 motif (e = 4.71e-07 using STAMP ). Mbp1p is a core member of the MBF complex with a consensus sequence very similar to that of SWI4 (e = 7.58e-05, STAMP), the motif we identified in S. bayanus. Consistent with our observations in S. bayanus, nucleosome depletion ( data from) at MBP1 motifs is correlated with periodic expression during the cell cycle (p = 6.27E-29, Mann-Whitney U test). Genes downstream of these depleted sites show peak expression at G1/S, the cell cycle stage regulated by SBF-MBF complexes  (Figure 1B). In contrast, genes with nucleosome-occupied motifs did not show a periodic pattern of expression. These results highlight that our observations can be generalized across species that have diverged by at least 20 million years.
Motif preservation for orthologous genes is correlated with nucleosome depletion in S. cerevisiae and S. bayanus
Correlation between changes in gene expression and changes of nucleosome-motif interactions across species
To investigate whether there is a correlation between nucleosome and motif interaction to cell-cycle expression conservation, we ranked the orthologous cell cycle regulated genes with conserved MBP1 motifs according to whether they are cyclically expressed in S. bayanus. Nucleosome depletion is a clear signal for those genes that remain periodically expressed in S. bayanus (Figure 4B, p = 7.2348E-05, Mann-Whitley U test). The conservation of periodic expression pattern is not an artifact of the existence of non-expressed genes: when we restrict the analysis to expressed genes (defined as the top 90% genes in expression level ), we observe the same pattern (Additional file 2, Figure S6). Similar results were obtained when testing the motif and gene expression conservation patterns in the opposite direction, i.e., from S. bayanus to S. cerevisiae (Additional file 2, Figure S7), supporting the generality of our observation and supporting the idea that motif gain/loss and changes in nucleosome-motif interactions contribute to alterations in gene expression patterns. This result suggests that alteration of nucleosome-motif interactions in the upstream regulatory region of genes is a biologically relevant phenomenon relating to gene expression divergence that is distinct from the loss of specific motifs in the regulatory region.
Orthologous motif analysis on expression-nucleosome changes across species
We also performed the analysis in the reverse direction, considering the conservation of nucleosome occupancy from S. bayanus to S. cerevisiae. These orthologous genes in S. cerevisiae were categorized into those that exhibit a cyclic pattern during the cell cycle and those that do not. The motifs upstream of periodically expressed genes show greater nucleosome depletion than genes which are not regulated by the cell cycle (Figure 5C, p = 0.05, Mann-Whitney U test), suggesting that for some classes of genes nucleosome depletion alone can be correlated to expression changes independent of the movement of transcription factor binding motifs.
Our observations provide two important perspectives on the current understanding of the control of gene expression. First our study provides strong evidence that cell-cycle regulation of gene expression is correlated not only to transcription factor binding motifs, but also to nucleosome occupancy at these sites. This genome-wide observation is consistent with a recent observation that nucleosome-depletion could ensure cyclic expression in two cell cycle-regulated promoters, CLN2pr and HOpr. It is important to note that the nucleosome data presented here are from an asynchronous culture, and therefore the reported occupancy reflects the integration of the nucleosome occupancy signal throughout the cell cycle. Nucleosome occupancy has been observed to fluctuate during the cell cycle in concert with periodic gene expression , and further studies at the level of individual genes and groups of genes will help to elucidate this cause-and-effect relationship. Indeed, deciphering specific cause-effect relationships between transcription factor binding and nucleosome depletion remains an important challenge.
Second, we show that changes in gene expression across species can be correlated to either loss of the regulatory motif or a change in nucleosome occupancy on conserved motifs. Compared to previous studies that focused on average expression level , we demonstrate that such correlation is significant for cell cycle expression regulation. Future work on nucleosome occupancy change on motifs between species and the resulting changes in gene expression may help, in part, to explain the disconnect between regulatory motif divergence and gene expression divergence. Specifically, it has been reported that most of the differences between species in transcription factor (TF)-binding motifs in yeasts and mammals have a very limited ability to predict gene expression divergence . Our results demonstrate that alteration of gene expression patterns can be related to a combination of motif turnover and nucleosome occupancy. Therefore, studies of transcription factor binding motifs may benefit from considering nucleosome occupancy data in parallel.
Our observations on the evolution of cell cycle gene expression suggest that phenotype divergence is strongly related to changes of nucleosome-motif interaction. However, despite this general trend, not all expression pattern changes observed in our analysis are explained by changes in nucleosome occupancy or by alteration of motifs. This observation likely reflects the real underlying biological sophistication of gene expression control. As such, our results contribute a novel phenomenon that should be considered in studies aimed at understanding the divergence of gene expression during evolution and its impact on evolutionary change. The datasets provided with this study comprise a rich resource for additional analyses (e.g. those that incorporate additional genome-wide, high-resolution data) to address these and related questions.
Through analyzing nucleosome occupancy pattern and cell cycle expression divergence between two yeast species S. bayanus and S. cerevisiae, we found that many changes in cell cycle gene expression patterns across species can be correlated to changes in nucleosome occupancy on motifs (rather than to the presence or absence of motifs). This observation suggests that alteration of nucleosome occupancy is a previously uncharacterized feature related to the divergence of cell cycle expression between species.
1. Nucleosome data preparation
Nucleosome sequencing data preparation
S, bayanus, also referred to as S. bayanus var uvarum, (623-6C) was obtained from Mark Johnston at the University of Colorado, Denver. In brief, yeast cells were grown in YPD at 27°C and harvested at mid-log phase (OD = 0.6). Nucleosomal DNA was isolated using protocols modified from  where mono-nucleosome sized DNA fragments were gel-purified and sequenced on the Illumina platform to produce 36 nt reads.
Mapping nucleosome reads
S. bayanus reference genome  was obtained from SGD (http://www.yeastgenome.org). We aligned nucleosome short reads (in total 11,841,061 reads) to this reference genome using Maq 0.7.1 (Mapping and Assembly with Quality) . This resulted in the number of matched reads at each chromosome location for sense and antisense strands. Nucleosome positions were defined according to . This nucleosome signal was further normalized by the ratio of the log value to the log of the median of the genome. The genome median is 1.
Modeling nucleosome positions for visualization
To visualize approximate nucleosome positions in S. bayanus, we applied NPS (Nucleosome Positioning from Sequencing)  to predict the positions of nucleosomes based on a wavelet model. Using FindPeaks 4.0, we converted the nucleosome signals based on the triangle distribution Additional file 2, Figure S2).
Given the estimated nucleosome positions, we calculated the percentage of nucleosome coverage of each contig by finding the ratio of nucleosomes within each contig (in bps) to the total length of each contig (in bps). In general, contigs with a reasonable size are well-covered, while very short contigs are likely to be poorly covered and also do not contain open reading frames. We therefore ranked the contigs by their coverage Additional file 2, Figure S1), and found 20% coverage (as determined by the nucleosome center regions by NPS ) to be a reasonable cutoff for well-covered contigs. To quantify the well-covered regions, we calculated the percentage of all contigs with greater than 20% nucleosome coverage (in bps) against the size of the entire genome. The open reading frames in these well-covered contigs were included in the following analyses. Two genes with unusually high numbers of mapped reads (> 100 times the genome average) were excluded from the analysis, and this resulted in the final 4935 genes included in the study.
2. Cell cycle data analysis
Cell cycle data and normalization
We acquired S. cerevisiae cell-cycle data from the alpha factor synchronization study described in Pramila et al. 2006  and S. bayanus alpha factor synchronization cell-cycle data described in  (available at GSE16544). Genes with more than half of the values missing were removed from the analyses and other missing values were imputed with KNNimpute , with K = 10, Euclidean distance. Genes with duplicate measurements were averaged. For each gene in the cell cycle data, the expression values were centered so that the average over the time course equals to 0. The orthologous gene matches of S. bayanus and S. cerevisiae were obtained from .
Identification of cell cycle frequency and the motifs that drive the cell cycle
For S. bayanus, we directly applied Fourier transform  to identify the major cell cycle frequency (Additional file 2, Figure S3A). Fourier transform decomposes the time-course expression into different frequencies. For S. bayanus, the frequency with maximum amplitude over all genes is a noise signal because the canonical cell cycle regulated genes such as MBP1 do not show maximum amplitude at this frequency. The frequency with the secondary maximum amplitude is the cell cycle frequency. Genes were ranked according to their differences from the cell cycle frequency.
For S. cerevisiae, a set of cell-cycle regulated genes were identified in . We identified the major cell cycle frequency based on the expression pattern of this set of genes (Additional file 2, Figure S3B) using Fourier transform. All genes were ranked according to their difference from the cell cycle frequency.
For both datasets, this ranking was used to identify the motifs most enriched in the cell cycle-regulated sets using FIRE . Motifs identified by FIRE  were mapped to the AlignACE database by STAMP . Often each motif could be mapped to a set of related motifs, for example, STB1, SWI4 and MBP1, in which case we labeled the motif with the best match. The positions of these motifs were then scanned throughout the whole genome using FIRE .
For each gene, we applied Fourier transform and identified the frequency with the maximum amplitude. The difference d of this frequency and the cell cycle frequency is calculated. The distribution of d g is plotted for all genes g, and the inflection point of this distribution is identified with the difference value of d inflection . Cell cycle-regulated genes were defined as those whose d g is smaller than d inflection , which indicates their pattern of expression is similar to cell cycle frequency. Conservation of cell cycle-regulation is defined as in both species, the periodicity difference is smaller than that of the inflection point.
The nucleosome signal around motifs and the cell cycle expression data were visualized in MeV4.6 by matching genes.
3. ChIP-chip data
We obtained Mbp1p binding sites data from , and mapped these sites to the S. cerevisiae May 2006 genome assembly, so as to match the nucleosome data we used. The genome-wide MBP1 positions were scanned using FIRE  based on the consensus defined in the AlignACE database.
4. Statistical analysis of changes in nucleosome occupancy
Calculating statistical significance between nucleosome occupancy and cell cycle periodicity
For each gene, we define its closeness to cell cycle periodicity based on the difference of its major frequency from the cell cycle frequency. For all the matches between a motif and its downstream gene, we used the Mann-Whitney U test to calculate the significance of the similarity of the distribution between the nucleosome signal and the gene's closeness to cell cycle periodicity.
Calculating statistical significance of nucleosome occupancy pattern on motifs
We used the following method to calculate the statistical significance when comparing nucleosome occupancy between conserved and non-conserved motifs, or between motifs driving periodic expression and non-functional motifs. For all occurrences of a motif, we identified the trough (lowest) position of the average nucleosome signal on these occurrences. The individual values on this trough position were recorded. The p value was calculated using Welch's t test.
For analysis of motifs that are conserved during evolution between S. cerevisiae and S. bayanus, we used alignment from . Only motifs that occur within the orthologous regions are considered in such analysis.
We thank James Broach and Amy Caudy for discussion and suggestions on the manuscript, Lance Parsons for technical support and Lance Parsons for technical support. OGT is supported by the NIH grants R01 GM071966 and R01HG005998, NSF grant IIS-0513552 and NSF CAREER award DBI-0546275. CN is supported by the Canadian Institutes for Health Research MOP-86705. YG, MJD and OGT were partially supported by the NIGMS Center of Excellence P50 GM071508, and by donations from the A. V. Davis Foundation and Princeton University for funding of QCB301, Experimental Project Laboratory. MJD is a Rita Allen Foundation Scholar and supported by NIH-P41-RR011823.
All nucleosome and expression data associated with this study is available at http://function.princeton.edu/cm/data/cellcycle_nucleosome/
Cell cycle expression data are available at GSE16544
- Lee W, Tillo D, Bray N, Morse RH, Davis RW, Hughes TR, Nislow C: A high-resolution atlas of nucleosome occupancy in yeast. Nat Genet. 2007, 39 (10): 1235-1244. 10.1038/ng2117.View ArticlePubMedGoogle Scholar
- Olins AL, Olins DE: Spheroid chromatin units (v bodies). Science. 1974, 183 (4122): 330-332. 10.1126/science.183.4122.330.View ArticlePubMedGoogle Scholar
- Kornberg RD: Chromatin structure: a repeating unit of histones and DNA. Science. 1974, 184 (139): 868-871. 10.1126/science.184.4139.868.View ArticlePubMedGoogle Scholar
- Bernstein BE, Liu CL, Humphrey EL, Perlstein EO, Schreiber SL: Global nucleosome occupancy in yeast. Genome Biol. 2004, 5 (9): R62-10.1186/gb-2004-5-9-r62.View ArticlePubMedPubMed CentralGoogle Scholar
- Ozsolak F, Song JS, Liu XS, Fisher DE: High-throughput mapping of the chromatin structure of human promoters. Nat Biotechnol. 2007, 25 (2): 244-248. 10.1038/nbt1279.View ArticlePubMedGoogle Scholar
- Schones DE, Cui K, Cuddapah S, Roh TY, Barski A, Wang Z, Wei G, Zhao K: Dynamic regulation of nucleosome positioning in the human genome. Cell. 2008, 132 (5): 887-898. 10.1016/j.cell.2008.02.022.View ArticlePubMedGoogle Scholar
- Segal E, Fondufe-Mittendorf Y, Chen L, Thastrom A, Field Y, Moore IK, Wang JP, Widom J: A genomic code for nucleosome positioning. Nature. 2006, 442 (7104): 772-778. 10.1038/nature04979.View ArticlePubMedPubMed CentralGoogle Scholar
- Yuan GC, Liu YJ, Dion MF, Slack MD, Wu LF, Altschuler SJ, Rando OJ: Genome-scale identification of nucleosome positions in S. cerevisiae. Science. 2005, 309 (5734): 626-630. 10.1126/science.1112178.View ArticlePubMedGoogle Scholar
- Owen-Hughes T, Workman JL: Experimental analysis of chromatin function in transcription control. Crit Rev Eukaryot Gene Expr. 1994, 4 (4): 403-441.PubMedGoogle Scholar
- Liu X, Lee CK, Granek JA, Clarke ND, Lieb JD: Whole-genome comparison of Leu3 binding in vitro and in vivo reveals the importance of nucleosome occupancy in target site selection. Genome Res. 2006, 16 (12): 1517-1528. 10.1101/gr.5655606.View ArticlePubMedPubMed CentralGoogle Scholar
- Morohashi N, Nakajima K, Kurihara D, Mukai Y, Mitchell AP, Shimizu M: A nucleosome positioned by alpha2/Mcm1 prevents Hap1 activator binding in vivo. Biochem Biophys Res Commun. 2007, 364 (3): 583-588. 10.1016/j.bbrc.2007.10.037.View ArticlePubMedPubMed CentralGoogle Scholar
- Sekinger EA, Moqtaderi Z, Struhl K: Intrinsic histone-DNA interactions and low nucleosome density are important for preferential accessibility of promoter regions in yeast. Mol Cell. 2005, 18 (6): 735-748. 10.1016/j.molcel.2005.05.003.View ArticlePubMedGoogle Scholar
- Badis G, Chan ET, van Bakel H, Pena-Castillo L, Tillo D, Tsui K, Carlson CD, Gossett AJ, Hasinoff MJ, Warren CL, et al: A library of yeast transcription factor motifs reveals a widespread function for Rsc3 in targeting nucleosome exclusion at promoters. Mol Cell. 2008, 32 (6): 878-887. 10.1016/j.molcel.2008.11.020.View ArticlePubMedPubMed CentralGoogle Scholar
- Zawadzki KA, Morozov AV, Broach JR: Chromatin-dependent transcription factor accessibility rather than nucleosome remodeling predominates during global transcriptional restructuring in Saccharomyces cerevisiae. Mol Biol Cell. 2009, 20 (15): 3503-3513. 10.1091/mbc.E09-02-0111.View ArticlePubMedPubMed CentralGoogle Scholar
- Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Tillo D, Field Y, LeProust EM, Hughes TR, Lieb JD, Widom J, et al: The DNA-encoded nucleosome organization of a eukaryotic genome. Nature. 2009, 458 (7236): 362-366. 10.1038/nature07667.View ArticlePubMedGoogle Scholar
- Shivaswamy S, Bhinge A, Zhao Y, Jones S, Hirst M, Iyer VR: Dynamic remodeling of individual nucleosomes across a eukaryotic genome in response to transcriptional perturbation. PLoS Biol. 2008, 6 (3): e65-10.1371/journal.pbio.0060065.View ArticlePubMedPubMed CentralGoogle Scholar
- Bai L, Charvin G, Siggia ED, Cross FR: Nucleosome-depleted regions in cell-cycle-regulated promoters ensure reliable gene expression in every cell cycle. Dev Cell. 2010, 18 (4): 544-555. 10.1016/j.devcel.2010.02.007.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsankov AM, Thompson DA, Socha A, Regev A, Rando OJ: The role of nucleosome positioning in the evolution of gene regulation. PLoS Biol. 2010, 8 (7): e1000414-10.1371/journal.pbio.1000414.View ArticlePubMedPubMed CentralGoogle Scholar
- Washietl S, Machne R, Goldman N: Evolutionary footprints of nucleosome positions in yeast. Trends Genet. 2008, 24 (12): 583-587. 10.1016/j.tig.2008.09.003.View ArticlePubMedGoogle Scholar
- Zhang Y, Shin H, Song JS, Lei Y, Liu XS: Identifying positioned nucleosomes with epigenetic marks in human from ChIP-Seq. BMC Genomics. 2008, 9: 537-10.1186/1471-2164-9-537.View ArticlePubMedPubMed CentralGoogle Scholar
- Cliften P, Sudarsanam P, Desikan A, Fulton L, Fulton B, Majors J, Waterston R, Cohen BA, Johnston M: Finding functional features in Saccharomyces genomes by phylogenetic footprinting. Science. 2003, 301 (5629): 71-76. 10.1126/science.1084337.View ArticlePubMedGoogle Scholar
- Ioshikhes IP, Albert I, Zanton SJ, Pugh BF: Nucleosome positions predicted through comparative genomics. Nat Genet. 2006, 38 (10): 1210-1215. 10.1038/ng1878.View ArticlePubMedGoogle Scholar
- Elemento O, Slonim N, Tavazoie S: A universal framework for regulatory element discovery across all genomes and data types. Mol Cell. 2007, 28 (2): 337-350. 10.1016/j.molcel.2007.09.027.View ArticlePubMedPubMed CentralGoogle Scholar
- Mahony S, Benos PV: STAMP: a web tool for exploring DNA-binding motif similarities. Nucleic Acids Res. 2007, W253-258. 35 Web ServerGoogle Scholar
- Hughes JD, Estep PW, Tavazoie S, Church GM: Computational identification of cis-regulatory elements associated with groups of functionally related genes in Saccharomyces cerevisiae. J Mol Biol. 2000, 296 (5): 1205-1214. 10.1006/jmbi.2000.3519.View ArticlePubMedGoogle Scholar
- Koch C, Moll T, Neuberg M, Ahorn H, Nasmyth K: A role for the transcription factors Mbp1 and Swi4 in progression from G1 to S phase. Science. 1993, 261 (5128): 1551-1557. 10.1126/science.8372350.View ArticlePubMedGoogle Scholar
- Pramila T, Wu W, Miles S, Noble WS, Breeden LL: The Forkhead transcription factor Hcm1 regulates chromosome segregation genes and fills the S-phase gap in the transcriptional circuitry of the cell cycle. Genes Dev. 2006, 20 (16): 2266-2278. 10.1101/gad.1450606.View ArticlePubMedPubMed CentralGoogle Scholar
- Mahony S, Auron PE, Benos PV: DNA familial binding profiles made easy: comparison of various motif alignment and clustering strategies. PLoS Comput Biol. 2007, 3 (3): e61-10.1371/journal.pcbi.0030061.View ArticlePubMedPubMed CentralGoogle Scholar
- Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, et al: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431 (7004): 99-104. 10.1038/nature02800.View ArticlePubMedPubMed CentralGoogle Scholar
- Doniger SW, Fay JC: Frequent gain and loss of functional transcription factor binding sites. PLoS Comput Biol. 2007, 3 (5): e99-10.1371/journal.pcbi.0030099.View ArticlePubMedPubMed CentralGoogle Scholar
- Martin OC, DeSevo CG, Guo BZ, Koshland DE, Dunham MJ, Zheng Y: Telomere behavior in a hybrid yeast. Cell Res. 2009, 19 (7): 910-912. 10.1038/cr.2009.65.View ArticlePubMedGoogle Scholar
- Hogan GJ, Lee CK, Lieb JD: Cell cycle-specified fluctuation of nucleosome occupancy at gene promoters. PLoS Genet. 2006, 2 (9): e158-10.1371/journal.pgen.0020158.View ArticlePubMedPubMed CentralGoogle Scholar
- Tirosh I, Weinberger A, Bezalel D, Kaganovich M, Barkai N: On the relation between promoter divergence and gene expression evolution. Mol Syst Biol. 2008, 4: 159-View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Ruan J, Durbin R: Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res. 2008, 18 (11): 1851-1858. 10.1101/gr.078212.108.View ArticlePubMedPubMed CentralGoogle Scholar
- Guan Y, Dunham M, Caudy A, Troyanskaya O: Systematic planning of genome-scale experiments in poorly studied species. PLoS Comput Biol. 2010, 6 (3): e1000698-10.1371/journal.pcbi.1000698.View ArticlePubMedPubMed CentralGoogle Scholar
- Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, Botstein D, Altman RB: Missing value estimation methods for DNA microarrays. Bioinformatics. 2001, 17 (6): 520-525. 10.1093/bioinformatics/17.6.520.View ArticlePubMedGoogle Scholar
- Byrne KP, Wolfe KH: Visualizing syntenic relationships among the hemiascomycetes with the Yeast Gene Order Browser. Nucleic Acids Res. 2006, D452-455. 34 DatabaseGoogle Scholar
- Zhu G, Spellman PT, Volpe T, Brown PO, Botstein D, Davis TN, Futcher B: Two yeast forkhead genes regulate the cell cycle and pseudohyphal growth. Nature. 2000, 406 (6791): 90-94. 10.1038/35017581.View ArticlePubMedGoogle Scholar
- Clifton P, Sudarsanam P, Desikan A, Fulton L, Fulton B, Majors J, Waterston R, Cohen AB, Johnston M: Finding functional features in Saccharomyces Genome by phylogenetic footprinting. Science. 2003, 301 (5629): 71-6. 10.1126/science.1084337.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.