The authors confirm that all data underlying the findings are fully available without restriction. The data have been published previously and can be accessed through the following accession numbers GSE36063 (Brogaard, 2012), E-MEXP-1172 (Lee, 2007), GSE25133 (Valouev, 2011), and SRA030538.1 (Zhang, 2011), and the following URL <http://bioinfo.stats.northwestern.edu/~jzwang/Yigit2013NAR.html> (Yigit, 2013).

Introduction {#s1}
============

*In vivo*, promoters are characterized by a nucleosome free region (NFR) that is followed by a periodic phasing of well-positioned nucleosomes continuing into the gene body. In yeast, this phasing is absent *in vitro*, but can be restored by the addition of a whole cell extract (WCE) and ATP, presumably a result of ATP-dependent chromatin remodelers (CRs) [@pone.0110479-Zhang1]. The promoter NFR, however, is largely preserved *in vitro* because yeast promoters contain sequences that are inherently refractory to nucleosome formation, such as low G/C content [@pone.0110479-Tillo1] and poly-dA:dT tracts [@pone.0110479-Kaplan1].

Yeast promoters have a biased distribution of poly-A and poly-T elements flanking nucleosome free regions [@pone.0110479-Wu1], [@pone.0110479-Hampson1], which cannot be explained solely by the biased base content ([**Figure 1**](#pone-0110479-g001){ref-type="fig"}). This asymmetric poly-A/poly-T arrangement has no known function and is incongruous with the model that poly-dA:dT tracts simply exclude nucleosomes via a rigid DNA structure [@pone.0110479-Iyer1] since the DNA should resist bending equally in either orientation.

![Yeast promoters have a biased distribution of poly-As and poly-Ts.\
The observed and expected frequency of poly-A and poly-T (AAAAA/TTTTT) elements across yeast promoters is shown, with expected calculated given the base content of the region. A greater number of poly-Ts and poly-As occur than expected in the −115:−75 and −75:−35 regions, respectively (p\<10^−6^ by simulation; see methods).](pone.0110479.g001){#pone-0110479-g001}

Results and Discussion {#s2}
======================

Hypothesizing that the asymmetric arrangement of these elements in promoters may have evolved to maintain promoter NFRs through some effect on nucleosome occupancy, we identified all non-overlapping poly-A sequences of exactly length five (AAAAA) in the yeast genome and analyzed the nucleosome occupancy [@pone.0110479-Zhang1] surrounding these elements ([**Figure 2**](#pone-0110479-g002){ref-type="fig"}). *In vitro*, both poly-A and poly-T sequences are similarly depleted of nucleosomes in an approximately symmetric fashion, both in the presence and absence of a WCE. However, upon addition of ATP, which activates CRs present in the WCE, the sequence becomes further depleted, but in an asymmetric fashion; a nucleosome becomes well-positioned 5′ to the poly-A sequence, but not 3′, and the NFR is offset 5′ to the poly-A sequence, similar to the trend observed *in vivo* ([**Figure 2**](#pone-0110479-g002){ref-type="fig"}).

![Nucleosomes are arranged asymmetrically around poly-dA:dT tracts.\
Average nucleosome occupancy surrounding poly-A and poly-T sequences (AAAAA/TTTTT) for salt gradient dialysis (*in vitro*), WCE without ATP (WCE-ATP), WCE with ATP added (WCE+ATP), as well as *in vivo* occupancy [@pone.0110479-Zhang1]. The difference in occupancy between poly-As and poly-Ts is significant only for *in vivo* and WCE+ATP (by rank sum; see **Figure S1** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**).](pone.0110479.g002){#pone-0110479-g002}

We next asked how nucleosomes were positioned around the three possible distinct arrangements of poly-A sequences (poly-A/poly-A, poly-A/poly-T, poly-T/poly-A). *In vivo* [@pone.0110479-Lee1], when two poly-A elements are within ∼60 bp, a strong NFR that is offset 5′ to the poly-A sequences generally results ([**Figure 3A**](#pone-0110479-g003){ref-type="fig"}). The poly-A/poly-T arrangement is typically much less depleted between the two motifs and yields two NFRs; one 5′ to the poly-A and the other 3′ to the poly-T ([**Figure 3B**](#pone-0110479-g003){ref-type="fig"}). The poly-T/poly-A combination results in the most robust NFR ([**Figure 3C**](#pone-0110479-g003){ref-type="fig"}), which could explain why this arrangement is preferred in yeast promoters. Further, in all cases, nucleosomes tend to be more well-positioned 5′ to poly-A sequences (3′ to poly-T). *In vitro*, in the absence of WCE and ATP [@pone.0110479-Zhang1], there is little difference between the three possible poly-A/poly-T arrangements and, in general, nucleosomes are depleted symmetrically around each poly-dA:dT element (**Figure S2** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**). We note that the occupancy bias surrounding poly-dA:dT tracts in the presence of active chromatin remodelers is unlikely to result from differences in the nucleosome isolation/quantification procedures because the same procedures were used to generate all *in vitro* data [@pone.0110479-Zhang1], but the bias occurs only when WCE and ATP are both present ([**Figure 2**](#pone-0110479-g002){ref-type="fig"}). Further, the nucleosome occupancy bias surrounding poly-A/poly-T combinations is consistent between *in vivo* datasets that use different approaches for crosslinking (sulfhydryl [@pone.0110479-Brogaard1], formaldehyde [@pone.0110479-Zhang1], [@pone.0110479-Lee1]), cleavage (peroxide-mediate [@pone.0110479-Brogaard1], MNase [@pone.0110479-Zhang1], [@pone.0110479-Lee1]), and quantification (microarray [@pone.0110479-Lee1], sequencing [@pone.0110479-Zhang1], [@pone.0110479-Brogaard1]; see **Figure S3** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**).

![The different poly-A/poly-T arrangements result in vastly different nucleosome occupancy outcomes.\
*In vivo* nucleosome occupancy [@pone.0110479-Lee1] (heatmap) surrounding all instances of (**A**) poly-A/poly-A, (**B**) poly-A/poly-T, and (**C**) poly-T/poly-A combinations in the yeast genome separated by no more than 500 bp. Red and blue curves represent the outer motif edges of poly-Ts and poly-As, respectively. Note that the poly-T/poly-T combination is a mirror image of the poly-A/poly-A data.](pone.0110479.g003){#pone-0110479-g003}

We hypothesize that the CR-dependent asymmetric arrangement of nucleosomes surrounding poly-A elements reflects differences in the nucleosome translocation efficiency from upstream vs. downstream of poly-As. It is possible that such a difference could result from the different histone-DNA contacts of the two DNA-strands. However, mouse [@pone.0110479-Yigit1] and human [@pone.0110479-Valouev1], which have nucleosomes very similar to those of yeast (84% identical in histone fold domains, between mouse and yeast), display a trend opposite to yeast ([**Figure 4**](#pone-0110479-g004){ref-type="fig"}); poly-A/poly-T combinations tend to be more depleted than poly-T/poly-A combinations, two consecutive poly-As generally result in 3′-biased NFRs, and, overall, there appear to be a more robust nucleosome boundaries 3′ to poly-As (5\' to poly-Ts). This observation suggests that specific factors (e.g. CRs) are responsible for differentiating between poly-As and poly-Ts. For example, poly-A tracts could prevent binding of CRs such that they can move a nucleosome towards poly-A sequences, but once there, the CR binds the DNA less efficiently and so cannot move it away. Indeed, previous studies have hinted that the DNA sequence could influence the repositioning of nucleosomes by CRs *in vitro*, but the mechanism, *in vivo* relevance, and sequence determinants of this phenomenon remained unknown [@pone.0110479-Rippe1], [@pone.0110479-vanVugt1]. More detailed studies of nucleosome positioning in the presence or absence of different CRs will be needed to determine the specificities of these CRs.

![Mammalian nucleosome occupancy is also biased surrounding poly-As and poly-Ts, but the trend is opposite to yeast.\
*In vivo* nucleosome occupancy for (**A--C**) regions with available high-resolution nucleosome data from mouse Th1 cells [@pone.0110479-Yigit1] and (**D--F**) non-repetitive regions on chromosome 22, for human granulocytes [@pone.0110479-Valouev1] (heatmaps) surrounding all instances of (**A**, **D**) poly-A/poly-A, (**B**, **E**) poly-A/poly-T, and (**C**, **F**) poly-T/poly-A combinations. Gaussian smoothed between rows (SD = 10 and 50, for mouse and human, respectively). The distinct transitions from light to dark in the mouse data (**A-C**) result from using unsmoothed data, which corresponds roughly to nucleosome dyad occupancy (in this case the poly-A/poly-T bias was more obvious without smoothing). This distinct transition is presumably caused by the destabilization of nucleosomes as poly-dA:dT tracts are incorporated, and nucleosomes appear to be most unstable when the dyad is 69 bp from the proximal poly-dA:dT tract edge in human, mouse, and yeast (**Figure S4** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**).](pone.0110479.g004){#pone-0110479-g004}

Our data indicate that poly-A sequences form an asymmetric barrier to CR-mediated nucleosome transit, that this asymmetry is used in yeast promoter architecture, and that the same sequences are used differently in mammals. This phenomenon helps explain part of the discrepancy between *in vitro* and *in vivo* nucleosome occupancy and indicates that the DNA sequence may play a greater role in positioning nucleosomes in the cell than previously appreciated. More complex models of nucleosome occupancy that account for CR-mediated nucleosome transit may be needed to fully explain nucleosome occupancy in the dynamic environment of the cell.

Methods {#s3}
=======

Definition of poly-A/poly-T {#s3a}
---------------------------

For [**Figure 1**](#pone-0110479-g001){ref-type="fig"} and [**Figure 2**](#pone-0110479-g002){ref-type="fig"}, we defined a poly-A element as any instance of five As in a row in the genome, with poly-T defined similarly. For [**Figure 1**](#pone-0110479-g001){ref-type="fig"}, we calculated the expected occurrence of poly-As and poly-Ts by using the nucleotide frequency at every base pair in the region to calculate the proportion of promoters expected to contain a poly-A or poly-T sequence at any given position. For [**Figure 2**](#pone-0110479-g002){ref-type="fig"}, we only considered non-overlapping instances.

The poly-A/poly-T combinations in [**Figure 3**](#pone-0110479-g003){ref-type="fig"}, and **Figure S2** and **S3** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}** were derived by identifying all maximal poly-A and poly-T elements of at least 5 bp in the yeast genome and considering only those motif pairs whose (outer) edges lie within 500 bp and that have no additional poly-dA:dT tracts between them. [**Figure 4**](#pone-0110479-g004){ref-type="fig"} was created similarly, but only considering BAC-enriched regions for mouse data (i.e. regions for which high-resolution occupancy data are available) and only non-repetitive (by repeatmasker) regions of chromosome 22 for human. For **Figure S4** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**, only poly-A tracts of exactly length 5 were considered. In all cases, we used the NCBI v37 mouse genome, hg18 human genome, and R64 yeast genome.

Nucleosome occupancy normalization {#s3b}
----------------------------------

For the data displayed in [**Figure 2**](#pone-0110479-g002){ref-type="fig"}, [**Figure 4**](#pone-0110479-g004){ref-type="fig"} (**D--E**), and **Figure S2** and **S3** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}** (*in vitro* and *in vivo* yeast sequencing data [@pone.0110479-Zhang1], [@pone.0110479-Brogaard1], and MNase-digested chromatin from human granulocytes [@pone.0110479-Valouev1]), we smoothed the data within each locus (Gaussian, SD = 20 bp), while for the data in [**Figure 3**](#pone-0110479-g003){ref-type="fig"} (yeast microarray data [@pone.0110479-Lee1]) and [**Figure 4** (**A--C**)](#pone-0110479-g004){ref-type="fig"} (mouse Th1 sequencing data, representing the centres of 147 bp fragments isolated from native, MNase-digested chromatin [@pone.0110479-Yigit1]), we performed no such smoothing. Smoothing the data in this way makes it correspond more closely to nucleosome occupancy by distributing the dyad occupancy (nucleosome centre position) over the area covered by a nucleosome. We did not smooth the mouse data because doing so obscured the poly-A/poly-T bias. We noted that the sequencing data ([**Figure 2**](#pone-0110479-g002){ref-type="fig"}, [**Figure 4**](#pone-0110479-g004){ref-type="fig"}, and **Figure S2--S4** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**) displayed significant variation in the number of reads per locus, so, for these data, we scaled each locus so that they each had a comparable numbers of reads and threw out any loci containing fewer than 40 (yeast *in vitro*; **Figure S2** and **S4A** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**), 400 (mouse *in vivo*; [**Figure 4A-C**](#pone-0110479-g004){ref-type="fig"} and **Figure S4C** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**), or 100 reads (human *in vivo*, [**Figure 4D--F**](#pone-0110479-g004){ref-type="fig"}, **Figure S4B** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**). For [**Figure 3**](#pone-0110479-g003){ref-type="fig"}, [**Figure 4**](#pone-0110479-g004){ref-type="fig"}, and **Figure S2** and **S3** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"},** we also smoothed between loci to emphasize the overall occupancy trend (Gaussian, SD = 50, except for [**Figure 4** **A--C**](#pone-0110479-g004){ref-type="fig"}, for which we used SD = 10).

Significance of poly-As and poly-Ts in promoter regions {#s3c}
-------------------------------------------------------

To gauge the significance of the distribution of poly-As and poly-Ts in promoter regions, we generated "random-sequence promoters" (in equal proportion to the number of actual promoters analyzed) where, at every base position, that base had the same probability of being an A or T as the actual frequency of that base at that position. We repeated this procedure 10^6^ times and each time counted the number of occurrences of 5 Ts or 5 As in a row within the −115:−75 and −75:−35 regions (relative to the TSS), respectively, but we found no randomly generated set of promoters with as extreme an occurrence of poly-As and poly-Ts in these regions as observed *in vivo* (max simulated = 1653 and 480, actual = 4919 and 2449 for A5 and T5, respectively).

Significance of nucleosome bias surrounding poly-dA:dTs {#s3d}
-------------------------------------------------------

To gauge the significance of the biased distribution of nucleosomes surrounding poly-As and poly-Ts, we compared the distribution of normalized (as described above) reads surrounding these sequences within each experimental condition. We used the two-tailed (Mann-Whitney) rank sum test to gauge the significance of the difference in occupancy for poly-As compared to poly-Ts at equivalent positions relative to the poly-A/T. The result is plotted in **Figure S1** in **[File S1](#pone.0110479.s001){ref-type="supplementary-material"}**, along with the Bonferroni multiple hypothesis correction significance threshold.

Supporting Information {#s4}
======================

###### 

Supplementary figures and figure legends.

(DOC)

###### 

Click here for additional data file.

[^1]: **Competing Interests:**The authors have declared that no competing interests exist.

[^2]: Conceived and designed the experiments: CGD. Performed the experiments: CGD. Analyzed the data: CGD. Contributed reagents/materials/analysis tools: CGD TRH. Contributed to the writing of the manuscript: CGD TRH. Supervised the research: TRH.
