# Evidence of selection for an accessible nucleosomal array in human

Guénola Drillon, Benjamin Audit, Françoise Argoul, Alain Arneodo

► To cite this version:

Guénola Drillon, Benjamin Audit, Françoise Argoul, Alain Arneodo. Evidence of selection for an accessible nucleosomal array in human. BMC Genomics, 2016, 17 (1), pp.526 (1-20). 10.1186/s12864-016-2880-2 . hal-01391690

HAL Id: hal-01391690

https://hal.science/hal-01391690

Submitted on 3 Nov 2016

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Evidence of selection for an accessible nucleosomal array in human

Guénola Drillon<fnref n="1" />, Benjamin Audit<fnref n="1" />, Françoise Argoul<fnref n="1" />,<fnref n="2" /> and Alain Arneodo<fnref n="1" />,<fnref n="2" />*

**Abstract**

**Background:** Recently, a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix has been used to reveal some enrichment of nucleosome-inhibiting energy barriers (NIEBs) nearby ubiquitous human "master" replication origins. Here we use this model to predict the existence of about 1.6 millions NIEBs over the 22 human autosomes.

**Results:** We show that these high energy barriers of mean size 153 bp correspond to nucleosome-depleted regions (NDRs) in vitro, as expected, but also in vivo. On either side of these NIEBs, we observe, in vivo and in vitro, a similar compacted nucleosome ordering, suggesting an absence of chromatin remodeling. This nucleosomal ordering strongly correlates with oscillations of the GC content as well as with the interspecies and intraspecies mutation profiles along these regions. Comparison of these divergence rates reveals the existence of both positive and negative selections linked to nucleosome positioning around these intrinsic NDRs. Overall, these NIEBs and neighboring nucleosomes cover 37.5 % of the human genome where nucleosome occupancy is stably encoded in the DNA sequence. These 1 kb-sized regions of intrinsic nucleosome positioning are equally found in GC-rich and GC-poor isochores, in early and late replicating regions, in intergenic and genic regions but not at gene promoters.

**Conclusion:** The source of selection pressure on the NIEBs has yet to be resolved in future work. One possible scenario is that these widely distributed chromatin patterns have been selected in human to impair the condensation of the nucleosomal array into the 30 nm chromatin fiber, so as to facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner.

**Keywords:** Nucleosome depleted regions, Sequence evolution, Nucleosome ordering, Mutation rates, GC content, Human genome

**Background**

To better understand the function and evolution of
eukaryotic genomes, it is of prime importance to deci-
pher to which extent the different steps of DNA com-
paction, from the nucleosomal primary DNA structure
to the chromatin loop tertiary DNA structure, have
been encoded in the DNA sequence during evolution
[1, 2]. Nowadays, it is well recognized that the dynam-
ics of DNA folding and unfolding within living cells
plays an important role in the regulation of nuclear func-
tions [3]. Transcription, replication and the associated

repair mechanisms are known to underlie the increasing
complexity of neutral mutation patterns in mammalian
genomes [4]. In human, all substitution rates increase
monotonously with replication timing, as an indication of
some weakening in the fidelity of the replication machin-
ery during the S phase [5, 6]. Over evolutionary time-
scales, transcription and replication have been shown
to induce strand-specific asymmetry in mutations [1, 7].
Because chromatin structure and dynamics vary across
the genome, one can expect that chromatin influences
the course of sequence evolution over evolutionary time
and contributes to the observed regional variations in
the mutational landscape. Recent studies have revealed
rather contrasting situations depending on whether DNA
accessibility has a greater impact on mutagenic or repair
processes. Lower (resp. higher) substitution frequency is

<footnote>*Correspondence: alain.arneodo@u-bordeaux.fr</footnote>

<footnote n="1">¹Univ Lyon, Ens de Lyon, Univ Claude Bernard Lyon 1, CNRS, Laboratoire de Physique, F-69342 Lyon, France</footnote>

<footnote n="2">²LOMA, Université de Bordeaux, CNRS, UMR 5798, 51 Cours de le Libération, F-33405 Talence, France</footnote>
generally observed in open (resp. close) chromatin envi-
ronment [8]. At a fine scale, the nucleosome was early
recognized as being at the heart of the trade-off between
the necessity of compacting DNA in the cell nucleus and
the required accessibility to regulatory proteins [1, 3].
Recent studies of sequence evolution in various eukary-
otic organisms have confirmed the existence of some
correlation between sequence divergence and nucleosome
positioning. In yeast, the overall levels of divergence were
observed to be higher near the nucleosome cores than
in the linker regions [9, 10] and interpreted as the con-
sequence of the limited access of nucleosomal DNA to
DNA repair proteins [11]. Alternative explanations have
been proposed, invoking purifying selection in linker
regions in particular to maintain nucleosome-disfavoring
codons in genic regions and nucleosome-antipositioning
sequences at regulatory regions like gene promoters and
enhancers [9, 12–14]. A similar observation was reported
in medaka fish [15]. Independant analysis of the first avail-
able genome-wide nucleosome maps in human [16–18]
have confirmed that, in mammals like in yeast, single
nucleotide polymorphisms (SNPs) are enriched around
the position of the nucleosome cores and depleted in
the linker regions. By combining the analysis of inter-
species and intraspecies divergence rates, Prendergast and
Semple [18] have further provided evidence for selection
likely acting on particular base substitutions to maintain
some optimum GC composition in both core (high GC
content) and linker (low GC content) DNA. By showing
that these signatures of selection concern non-coding as
well as coding sequences, in bulk as well as in epigeneti-
cally modified nucleosomes, these authors have provided
some very promising understanding of the important role
played by the DNA sequence in nucleosome positioning.
These signatures of selection were derived from a rela-
tively small subset of well-positioned nucleosomes that
cover about 5 % of the human genome. If this coverage
is comparable to protein coding genes coverage, one may
wonder whether the nucleosomal 1D organization of the
chromatin fiber has a more widespread impact on the
divergence of regions traditionally considered as evolving
neutrally, including intronic and intergenic sequences.

The role of the genomic sequence on nucleosome posi-
tioning in vivo has become an important issue [19–22].
It was early recognized that the periodic occurence of
dinucleotides with a period of 10 – 11 bp favors the for-
mation of nucleosomes [23–26]. However, in *S. cerevisiae*,
Drosophila and mammals, this periodicity was shown to
account for no more than 20 % of in vivo nucleosome posi-
tioning above what is expected by chance [19, 27, 28] and
it is almost invisible in the human genome after masking
the recently integrated Alu-sequences [29]. An alterna-
tive to the tight histone binding obtained with favor-
able positioning sequences is the statistical positioning of

nucleosomes nearby nucleosome inhibitory energy barriers (NIEBs), that can be encoded either via unfavorable sequences that potentially resist to the structural distortions required by nucleosome formation or particular sequences that may recruit transcription factors, or/and other protein complexes such as chromatin regulators that may compete with the nucleosomes [19, 21, 22]. Due to the interplay between boundary confinement and nucleosome-nucleosome excluded volume interactions, some statistical short-range ordering can establish near a NIEB, and gradually disappear away from the barrier [19, 30–33]. In that context, a simple physical model of nucleosome assembly based on the computation of the free energy cost of bending a DNA fragment of a given sequence from its natural curvature to the final superheli-cal structure around the histone core, was recently shown to mimic in vitro nucleosome occupancy data remarkably well [19, 34, 35]. When compared to in vivo data in *S. cerevisiae* and *C. elegans*, this sequence-dependent ther-modynamical model performs as well as models based on statistical learning, suggesting that in these organisms, the in vivo nucleosome array organization is, to a large extent, controlled by the underlying genomic sequence, although it is also subject to the finite-range remodeling action of external factors.

A key issue to the understanding of chromatin-mediated
regulation of nuclear functions is the relative position-
ing of regulatory sites with respect to NIEBs encoded
in the DNA sequence. Active regulatory sequences are
commonly thought to lie in accessible chromatin regions.
This was verified in vivo in eukaryotic genomes but not
necessarily in vitro, as recently noticed in *Arabidopsis*
[36], Drosophila [37], mouse [37–40] and human [41, 42].
In *S. cerevisiae*, nucleosome-depleted regions (NDRs)
were mainly observed in vivo at transcription start
sites (TSS), transcription termination sites (TTS)
[31, 41, 43–45] and at active DNA replication origins
[46, 47]. Interestingly, most of these NDRs correspond to
intrinsic NDRs encoded in the DNA sequence, up to some
shape remodeling and phasing. If similar sequence-driven
NDR regulation of transcription and replication initiation
is likely to operate in different yeast species [48–50] and
possibly in *C. elegans* [51], the situation appears to be
different in mammals, where regulatory sequences were
shown to have higher-than-average intrinsic nucleo-
some occupancy. Unlike in yeast where genes are mostly
constitutively expressed and thus would be expected to
maintain low GC promoters, most human promoters
have a high GC content, that indeed corresponds to
nucleosome preference in vitro [41, 42] and reduced
frequency of rigid poly(dA:dT) sequences that impair
nucleosome formation and favor nucleosome disassembly
[23, 31, 52–55]. Recent experimental and bioinformatics
studies in the mouse genome [37, 40] have confirmed
that nucleosomes are intrinsically preferentially positioned at replication origins and removed when actived, as previously observed in *Arabidopsis* [36]. Altogether, these recent results suggest that, in higher eukaryotes, a high nucleosome affinity is directly programmed at regulatory sequences to intrinsically restrict access to regulatory information that will mostly be used in vivo in an epigenetically-controlled cell-type-dependent manner [41]. This probably explains that so far, only little attention has been paid to intrinsic NDRs in mammal chromatin.

Here we analyze in vitro and in vivo nucleosome occupancy data in human, using as a reference the set of NIEBs predicted by the physical model based on sequence-dependent DNA bending properties [19, 34, 35]. We show that these intrinsic NDRs, together with the flanking nucleosomes, are actually encoded in the GC content and are widely spread over the 22 human autosomes, covering about 37.5 % of the genome. When further investigating interspecies and intraspecies rates of divergence, we bring significant evidence of selection pressure to maintain both an optimal GC content at flanking nucleosomes and a robust GC content depletion in NDRs relative to the local bulk GC content. We comment on the observed deviations from neutral evolution as a possible indication of the selection of an open, accessible and dynamic 10 nm chromatin fiber to constitutively facilitate in each cell-type the epigenetic regulation of nuclear functions.

# Results

## Intrinsic NIEBs correspond to in vitro and in vivo NDRs

When running the physical model based on sequence-dependent DNA bending properties (Methods) over the 22 human autosomes, we predicted a total of 1,581,256 NIEBs of width ranging from 36 to 450 bp with a mean size of 153 bp (Additional file 1: Figure S1) comparable to the mean size of the NDRs observed in the nucleosome occupancy experimental data [56, 57]. As expected, the histogram of GC content computed in the NIEBs is significantly ($P \ll 10^{-3}$, see Methods) shifted towards lower values ($\overline{GC} = 26.5 \%$, Additional file 1: Figure S2B) as compared to the histogram obtained genome-wide ($\overline{GC} = 37.9 \%$, Additional file 1: Figure S2A) and to the histogram computed over the 300 bp flanking regions on the left and the right of these NIEBs ($\overline{GC} = 42.5 \%$, Additional file 1: Figure S2C). Hence NIEBs correspond to AT-rich regions likely containing rigid poly(dA:dT) tracks known to be important determinants of nucleosome depletion in vitro as well as in vivo. This is consistent with the observation that these regions also resist to DNase I cleavage and on average are less sensitive to DNase I than the genome average (normalized DNase I sensitivity raw tag density = 0.82, $P \ll 10^{-3}$).

When investigating nucleosome density inside these energetically unfavourable regions, we confirmed that

they correspond to NDRs as compared to genome average (Fig. 1a) and this not only for “Valouev” in vitro data (normalized nucleosome tag density = 0.15) but also for “Valouev” in vivo data (0.60) and “Schones” in vivo data (0.44) (Methods). Given the variability in sequencing depth in nucleosome MNase-seq data [42, 57], we tested whether the observed “mean” depletion reflects MNase-seq profile across individual NIEBs by comparing nucleosome tag densities within the barriers and in the two 300 bp flanking windows (Methods). We first excluded from our analysis the 43,364 NIEBs and flanking windows (2.7 % of our sample) lacking sufficient mappability to perform the test (Methods). Then, from this comparison, for each MNase-seq data set, we classified the 1,537,892 mappable NIEBs into 3 groups (Methods): (G1) NIEBs significantly depleted in nucleosome tags at the 5 % confidence level as compared to flanking windows; (G2) NIEBs significantly enriched in nucleosome tags at the 5 % confidence level as compared to flanking windows; and (G3) the remaining NIEBs. For “Valouev” in vitro data, we found that 76.7 % of the NIEBs are in G1 and correspond to regions where the experimental in vitro nucleosome tags density is in agreement with the NIEB predictions, and only 0.2 % are in G2 and have a significant nucleosome occupancy signal over the predicted NIEBs (Additional file 1: Figure S3A). Interestingly, similar distributions over the three groups were obtained when analyzing “Valouev” in vivo data (Additional file 1: Figure S3B), and also “Schones” in vivo data (Additional file 1: Figure S3C). Despite a significantly weaker sequencing depth in “Schones” in vivo data, we obtained similar percentages of NIEBs in G1 for the former (54.0 %) than for the latter (58.7 %), confirming that a majority of NIEBs correspond also to NDRs in vivo. Note that the corresponding percentages of NIEBs in (G2), namely 0.5 % in “Valouev” data and 0.9 % in “Schones” data, are both very small, as previously observed with in vitro data. Thus, independently of the initial strategy used to identify these regions of interest, the important point is that these AT-rich regions impair nucleosome formation in vitro and in vivo.

## In vivo nucleosome ordering near NIEBs is encoded in the DNA sequence

As expected from statistical positioning [19, 30], these NIEBs are flanked by a few (~ 2-3) rather well-positioned nucleosomes with a nucleosome repeat-length (NRL) ~ 160 bp robustly observed in vivo and slightly more pronounced in vitro (Fig. 1a, a'). This nucleosome ordering is in remarkable agreement with the normalized mean nucleosome density profiles predicted by the sequence-dependent physical model at low and high genome nucleosome coverages (Methods) with a robust NRL ~ 150 bp (Fig. 1b, b'). This indicates that not only the NDRs but the neighbouring nucleosome arrangement are programmed
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Normalized (with respect to genome average) mean nucleosome density on both sides of the 1,581,256 NIEBs predicted by the sequence-dependent physical model. **a** “Schones” in vivo (brown), “Valouev” in vivo (pink) and “Valouev” in vitro (purple) data (Methods). **b** Numerical mean profiles predicted by the physical model at low (dark green) and high (light green) genomic nucleosome coverages (Methods). **c** Mean GC content (blue), repeat-masked GC content (sky blue) and GC content at equilibrium (navy blue). **a', b', c'** are zooms on the profiles in (a,b,c) on the right-hand side of the in silico NIEBs; vertical blue lines correspond to local minima of the GC content (they will be reported in Figs. 4, 5, 6 and 7 for comparison). All profiles were computed at 1 bp resolution, except the ‘GC at equilibrium’ which was smoothed over 10 bp windows (its non-smoothed profile is shown as a background in (c)).</figcaption></figure>

in the DNA sequence. Indeed, the local GC content (Fig. 1c, c') provides an excellent prediction of the experimental mean nucleosome occupancy (i.e. the probability that a given locus belongs to a DNA fragment involved in the nucleosome complex) profiles and this not only for in vitro data, as previously pointed out in various organisms, including *S. cerevisiae* [19, 53, 54, 58], *C. elegans* [19, 51, 53, 58] and human [41, 42], but also for in vivo data. Let us point out that a quite similar mean GC content profile is obtained when considering regions of the chimpanzee genome homologous to the human

intrinsic NDRs (data not shown, Methods) suggesting that the nucleosome ordering associated with NIEBs is robust across primate evolution. Note that the intrinsic mean nucleosome spacing predicted by the physical model (Fig. 1b, b') likely encoded in the local GC-content (Fig. 1c, c') and consistently observed in vitro and in vivo (Fig. 1a, a'), namely NRL ~ 150 – 160 bp, is significantly smaller than the in vivo genome average NRL ~ 203 bp, the average heterochromatin NRL ≳ 205 bp and more surprisingly the average NRL ~ 178 bp in euchromatin around active promoters and enhancers estimated
in CD4⁺ [42]. The highly compacted nucleosome arrangement predicted and observed on both sides of the NIEBs with a rather short DNA linker size ~ 10 – 20 bp, is likely to affect and probably to impair the condensation of the nucleosomal array into the 30 nm chromatin fiber [59–63].

We reproduced this analysis of nucleosome occupancy nearby NIEBs according to the G1, G2 or G3 group classification. For each set of experimental data, taking into account the mappability of the considererd regions (Additional file 1: Figure S4), we plotted the mean nucleosome density (Additional file 1: Figure S5) and the mean GC content profile (Additional file 1: Figure S6). What-ever the group, the GC profile displays a robust oscillatory pattern similar to the one previously observed on the genome-wide mean GC profile (Fig. 1c, c'). As expected, a compact nucleosome arrangement like the one observed genome-wide (Fig. 1a, a') was found around the G1 energy barriers that indeed constitute a large majority of predicted NDRs in the “Valouev” in vitro and also in the “Valouev” and “Schones” in vivo data. More importantly, when averaging the nucleosome tag densities around the set of 43,364 NIEBs lacking sufficient mappability to perform the test (Additional file 1: Figure S4), we recovered a consistent well ordered but less pronounced nucleosome ordering flanking the predicted NDRs in the in vitro and also in the two in vivo nucleosome data sets (Additional file 1: Figure S5). The additional observation that no significant departure from the compact nucleosome ordering picture was observed around the G3 energy barriers resulted in the following overall estimates: among the 1,581,256 NIEBs, 1,578,827 (99.8 %) in “Valouev” in vitro data, 1,573,887 (99.5 %) in “Valouev” in vivo data and 1,566,659 (99.1 %) in “Schones” in vivo data, correspond to experimentally observed NDRs bordered by rather well-positioned and strongly compacted

neighbouring nucleosomes. Thus, only an extreme minority (< 1 %) of the predicted NIEBs, namely the ones that belong to the group G2, were found to have a nucleosome occupancy higher than their two flanking 300 bp windows.

## NIEBs and intrinsic flanking nucleosomes are widely distributed over human autosomes

To characterize the spatial distribution of these NIEBs along human chromosomes, we performed a statistical analysis of the border to border interdistances between successive NIEBs. The obtained histogram (Fig. 2) displays an exponential tail with a characteristic interdistance compatible with the mean interdistance $\bar{d} = 1.54$ kb as the signature of a Poisson-like distribution. Interestingly, for interdistances $d < 1$ kb, this histogram switches to a quantized distribution with peaks equally separated by a remarkably stable and robust distance 153 bp, quite similar to the characteristic DNA length 147 bp involved in the nucleosome complex. We generated two-dimensional (2D) maps of nucleosome density that reveal a strikingly organized nucleosome ordering along inter-NIEB regions (Fig. 3). A crystal-like nucleosome organization with a well-defined number of regularly spaced nucleosomes is clearly observed in both “Schones” (Fig. 3a) and “Valouev” (Fig. 3b) in vivo data. As the interdistance $d$ increases, the 1,581,005 inter-NIEB regions clusters into classes displaying the same number $n$ of nucleosomes, from $n = 1$ (145,463; 9.2 %), 2 (147,922; 9.4 %), 3 (117,678; 7.4 %), 4 (116,399; 7.4 %), to 5 (108,324; 6.9 %) for $d < 0.8$ kb. For larger $d$, the periodic nucleosome positioning remains visible, but appears fuzzy as the signature of some weakening of the statistical positioning away enough from the two confining bordering NIEBs [19, 30, 32]. Nevertheless, a close inspection of the in vivo experimental 2D maps shows that if the central nucleosomes have no

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Histogram of (1,581,005) border to border inter-distances between successive NIEBs. The inset corresponds to a log-lin representation of the tail of this histogram, which puts into light a Poisson-like exponential behavior with a mean interdistance $\bar{d} = 1,541$ bp (red straight line). The green vertical lines mark interdistances $d = 117 + 153 k$ (bp), $k = 0, 1, 2, 3$ and 4</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Heat map of nucleosome density and GC content in between NIEBs. The 1,581,005 inter-NIEB regions were centered at 0 and ordered vertically from the smallest (top) to the largest (bottom). **a** “Schones” in vivo data; **b** “Valouev” in vivo data; **c** “Valouev” in vitro data; **d** in silico nucleosome density computed at high nucleosome coverage (Methods); **e** GC content. Only the central part of inter-NIEB regions larger than 1.6 kb were shown to enlighten the absence of nucleosome positioning in these regions. Each horizontal line represents the mean nucleosome coverage over 200 inter-NIEB regions coded from white (0) to full color (1) (Methods)</figcaption></figure>

preferential positioning, the first two nucleosomes flanking the NIEBs are well positioned and highly compacted with a very short NRL ≃ 153 bp and this even for the largest inter-NIEB regions. This remarkable crystal-like nucleosomal chromatin architecture reminds us of the intragenic nucleosome ordering previously observed in S. cerevisiae [32], except that contrary to its vanishing in in vitro data in yeast, this periodic ordering is still observed in “Valouev” in vitro data in human (Fig. 3c) and in the physical model at low nucleosome coverage (Fig. 1b, b') as the footprint of some underlying positioning sequences. This conclusion is corroborated by the quite similar 2D maps predicted by our physical model at low (in vitro) and high (in vivo) nucleosome coverages (Methods) that remarkably reproduce (Fig. 3d) the nucleosomal patterns observed in vivo (Fig. 3a, b) and in vitro (Fig. 3c). Indeed, when representing the corresponding heat map obtained for the GC content (Fig. 3e), we confirmed that nucleosome occupancy nearby NIEBs is intrinsically encoded in the local GC content, in the crystal-like nucleosome organization observed in small (d < 0.8 kb) inter-NIEB regions, as well as in the first two flanking nucleosomes for larger (d > 0.8 kb) inter-NIEB regions. Altogether, relative to the overall 2,681,301,120 bp of sequenced DNA in the 22 human autosomes, the 1,581,256 NIEBs cover 241,129,337 bp (9.0 %), the 635,786 inter-NIEB regions with crystal-like nucleosome ordering (n = 1 to 5) cover 254,831,369 (9.5 %) and the well-positioned first two nucleosomes flanking the 945,219 widely spaced (d > 800 bp) successive NIEBs cover 510,418,260 bp (19.0 %). This means that in 37.5 % of the human genome, the

GC content provides a remarkable prediction of nucleosome occupancy. Thus, in more than a third of the human genome, in vivo nucleosome positioning is intrinsically encoded in the DNA sequence without evidence of local chromatin remodeling. These findings also suggest that nucleosome positioning could be at the basis of local GC content variation as opposed to a mutational basis such as bias gene conversion [64–66].

With regard to the large scale organization of mammalian genomes [4, 67], we investigated the density of NIEBs in 100 kb non-overlapping windows that we further classified according to their GC content, mean replication timing (MRT) data, DNase I sensitivity raw tag density and meiotic recombination rate respectively (Table 1). We used MRT [68] and DNase I sensitivity [69] data from the lymphoblastoid cell line Gm06990 as surrogates of corresponding data in CD4⁺ (Methods). As compared to the genome-wide mean density = 0.59 NIEB/kb, the densities obtained in the GC-poor light isochore regions L1 (0.54 NIEB/kb) and L2 (0.63 NIEB/kb) are not so different from the ones obtained in the GC-rich heavy isochore regions H1 (0.65 NIEB/kb), H2 (0.56 NIEB/kb) and H3 (0.44 NIEB/kb) with a significant enrichment in L2 and H1 relative to some significant depletion in L1 and H3 (P ≪ 10⁻³, Table 1). This non-monotonous variation of the NIEB density versus the GC content contrasts with the anticorrelation observed with the MRT, the early replicating windows containing relatively more NIEBs (0.65 NIEB/kb) than mid S phase windows (0.59 NIEB/kb) and late replicating windows (0.54 NIEB/kb) (P ≪ 10⁻³, Table 1). Consistently with previous observation that early
<table><caption>**Table 1** Large-scale distribution of the predicted 1,581,256 NIEBs along the 22 human autosomes</caption><thead><tr><th></th><th>(% genome)</th><th>NIEB/kb</th><th>cov (%)</th></tr></thead><tbody><tr><td>All NIEBs</td><td>100</td><td>0.59</td><td>9.0</td></tr><tr><td>L1 (GC &lt; 38 %)</td><td>33.3</td><td>0.54</td><td>9.5</td></tr><tr><td>L2 (38 ≤ GC &lt; 42 %)</td><td>32.3</td><td>0.63</td><td>9.6</td></tr><tr><td>H1 (42 ≤ GC &lt; 47 %)</td><td>21.4</td><td>0.65</td><td>9.0</td></tr><tr><td>H2 (47 ≤ GC &lt; 52 %)</td><td>9.2</td><td>0.56</td><td>6.9</td></tr><tr><td>H3 (52 % ≤ GC)</td><td>3.7</td><td>0.44</td><td>5.0</td></tr><tr><td>Early (MRT &lt; 0.36)</td><td>30.0</td><td>0.65</td><td>9.1</td></tr><tr><td>Medium (0.36 ≤ MRT &lt; 0.69)</td><td>40.0</td><td>0.59</td><td>8.9</td></tr><tr><td>Late (0.69 ≤ MRT)</td><td>30.0</td><td>0.54</td><td>9.1</td></tr><tr><td>Low DNase (&lt; 14.3 reads/kb)</td><td>30.0</td><td>0.53</td><td>9.0</td></tr><tr><td>Medium DNase</td><td>40.0</td><td>0.60</td><td>9.3</td></tr><tr><td>High DNase (&gt; 29.3 reads/kb)</td><td>30.0</td><td>0.63</td><td>8.7</td></tr><tr><td>Low Recomb (&lt; 0.378 cM/Mb)</td><td>30.0</td><td>0.61</td><td>9.7</td></tr><tr><td>Medium Recomb</td><td>40.0</td><td>0.59</td><td>9.0</td></tr><tr><td>High Recomb (&gt; 1.681 cM/Mb)</td><td>30.0</td><td>0.58</td><td>8.4</td></tr><tr><td>Genes</td><td>47.1</td><td>0.61</td><td>9.2</td></tr><tr><td>Intergenes</td><td>51.6</td><td>0.57</td><td>8.9</td></tr></tbody></table>

NIEB density was computed in 100 kb non-overlapping windows that were classified according to their GC content, MRT, DNase I sensitivity raw tag density and meiotic recombination rate (Methods). The last two rows correspond to mean NIEB densities computed in 23,329 genes (Methods) and complementary intergenic regions (Methods). The first column defines the region, the second column its corresponding genome coverage (%), the third column the mean NIEB density (NIEB/kb) and the fourth column, the corresponding coverage of this region by the set of NIEBs (%). For each row, we estimated the standard error of the mean NIEB density to be SEM < 0.01

replicating regions are more sensitive to DNase I cleavage than late replicating regions [70, 71], we found that the NIEB density indeed correlates with DNase I tag density, the less sensitive 100 kb windows containing relatively less NIEBs (0.53 NIEB/kb) than the most sensitive (0.63 NIEB/kb) ($P \ll 10^{-3}$, Table 1). No significant correlation was observed with the meiotic recombination rate with only a slight enrichment (0.61 NIEB/kb) in windows with low cross-over rate (Table 1). Since GC-rich, early replicating open chromatin regions accessible to DNase I digestion are known to be highly genic regions [72], we further investigated the NIEB density in genic and intergenic regions. Human genes contain slightly more NIEBs (0.61 NIEB/kb) than intergenic regions (0.57 NIEB/kb) ($P \ll 10^{-3}$, Table 1). But as observed genome-wide, the density of NIEBs in the intergenic regions displays quite comparable anti-correlation with MRT and correlation with DNase I tag density (Additional file 1: Table S1), indicating that these trends do not simply reflect the difference in gene density in early replicating transcriptionally active open chromatin-regions relative to late replicating gene desert heterochromatin regions. Interestingly, as compared to genome average, the more (resp. less

accessible euchromatin (resp. heterochromatin) regions were found to be enriched (resp. depleted) in crystal-like nucleosomal patterns bordered by two NIEBs and this in genic as well as in intergenic regions (Additional file 1: Table S2). Fuzzy nucleosome occupancy in larger ($d > 0.8$ kb) inter-NIEB regions were more homogeneously observed along human autosomes without significant correlation with GC content, MRT, DNase I tag density and meiotic recombination rate (Additional file 1: Table S2). Altogether these results show that the predicted 1,581,256 NIEBs are widely distributed along human autosomes, with a density at 100 kb resolution that robustly ranges from 0.54 to 0.65 NIEB/kb, meaning that on average, one can expect to find a NIEB every 1.5–1.8 kb along human chromosomes and this whatever the genomic (nucleotide composition, gene, intergene) and epigenetic (MRT, DNase I sensitivity) context. When adding to these NIEBs, the flanking nucleosomes that were also shown to be encoded in the local GC content (Fig. 3), this is indeed a ~ 1 kb long constitutive nucleosomosomal chromatin pattern that is widely distributed along human chromosomes covering overall more than a third of the genome.

## Complex patterns of divergence around NIEBs

To detect possible evidence of selection linked to the nucleosomal organization of the chromatin fiber, we first used single nucleotide substitution rates tabulated in the human lineage since its divergence from chimpanzee using macaque and orangutan as outgroups [6] (Methods). We found that with respect to the mean substitution rate (0.50 %) obtained when averaging over the 22 human autosomes, the level of divergence in the NIEBs is significantly lower (0.46 %) ($P \ll 10^{-3}$, Fig. 4a), revealing that these regions unfavourable to nucleosome formation were highly conserved, at least during recent evolution. As previously observed in yeast [9–13], medaka fish [15] and human [16–18], we recovered that the overall levels of divergence are higher than genome average in the well-positioned nucleosomes flanking the NIEBs and lower in linker regions but not as low as in the NIEBs. Interestingly, quite similar substitution rate patterns were obtained in genes and in intergenic regions (Fig. 4b), around a mean substitution rate that is larger (0.53 %) in intergenic regions than in genes (0.47 %) ($P \ll 10^{-3}$), consistent with previous reports [5, 6] that higher substitution frequency is generally observed in late replicating, gene poor, heterochromatin regions. We computed the GC compositional profile toward which the sequences are evolving as the result of the so-computed substitution rate matrix (Methods). The resulting profile of GC composition at equilibrium (Fig. 1c, c') is in remarkable agreement with the “Valouev” in vitro and “Valouev” and “Schones” in vivo experimental nucleosome occupancy data (Fig. 1a, a'), with a NRL ~ 155-160 bp. This further
<figure><img src="image_6.png" /></figure>

supports the interpretation of local GC variation as a
physical basis accommodating for a specific nucleosomal
positioning pattern. When comparing this equilibrium
GC profile to the GC profiles obtained from the native
and repeat-masked sequences, some additional oscilla-
tions are present in the core of the flanking nucleosomes
in the native GC profile, suggesting the presence nearby
the NIEBs of sequences that do not reflect long-term evo-
lutionary compositional patterns, e.g. repeated sequences
that have been inserted recently in the human genome.

Next we examined the relative contributions of individ-
ual base changes. Strong peaks of S (C/G) → W (A/T)
substitution rates were observed in the NIEBs (Fig. 4c, e)
with matching low rates of W → S substitution rates in
the same regions (Fig. 4d, f). As a consequence, there is a
clear preference for substitutions increasing the AT con-
tent of these NIEBs. A similar preference was observed in
the linker regions between the successive flanking nucleos-
osomes but this preference progressively vanishes when
moving away from the NIEBs (Fig. 4c–f). Different sub-
stitution patterns were obtained at the well-positioned
flanking nucleosomes with high rates of W → S substitu-
tions (Fig. 4d, f) coming along with low rates of S → W
substitutions (Fig. 4c, e). This combination of divergence

patterns contributes to increasing the GC content at the
flanking nucleosome positions, consistent with previous
analysis of well-positioned nucleosome in “Schones” in
vivo data by Prendergast and Semple [18]. Note that its
importance progressively damps out when considering
nucleosomes further away from the NIEBs as expected
from statistical ordering near an excluding energy
barrier [19, 30].

Using single nucleotide polymorphism (SNP) data from
several fully sequenced human genomes [73], we fur-
ther investigated the rates of intraspecies divergence
(Methods). The broad variation in total polymorphism
density in and around the set of NIEBs is somehow dif-
ferent from the overall pattern of interspecies divergence
(Fig. 4a). Among the dominating intraspecies divergence
rates, we recovered in the recent human lineage strong
peaks of S → W substitution rates in the NIEBs and to a
lesser extent in the linker regions (Fig. 4c, e) that correlate
with low rates of W → S substitutions in the same regions
(Fig. 4d). Also both interspecies and intraspecies T → A
and A → T substitution rates display some increase in
the NIEBs (Fig. 4g). As far as the well-positioned flanking
nucleosomes are concerned, we consistently found high
rates of A → G and T → C intraspecies substitutions
(Fig. 4d) with matching lower C → A and G → T substitution rates (Fig. 4e). But the comparison of inter- and intraspecies divergence rates also revealed striking differences for certain classes of substitutions. In particular the T → G and A → C intraspecies substitution rate profiles were found rather flat over the flanking nucleosome regions with some relative increase (and not decrease) in the NIEBs (Fig. 4f). These observed differences in rates of interspecies and intraspecies divergence strongly suggest the presence of some underlying modes of selection associated with nucleosome ordering nearby the NIEBs.

**The genomic sequences associated with NIEBs and flanking nucleosomes are not evolving neutrally**

To investigate the possible existence of selection in the regions of interest, we used McDonald-Kreitman (MK) test and its variants [74] which amounts to compare rates of fixed interspecies divergence and intraspecies polymorphisms. Positive (resp. negative) selection is expected to lead to an excess (resp. depletion) of interspecies divergence relative to intraspecies divergence. To avoid to be biased by confounding factors like altered rates of mutation and repair [18] that potentially can induce differences

in substitution rates in some regions, we not only computed the rates of inter- and intraspecies divergence at the sites of interest but also the rates of inter- and intraspecies divergence at distant loci with no preferential nucleosome positioning as a “proxy” of the rate of neutral divergence (Methods). Practically, we used the average rates of inter- and intraspecies divergence computed in the central regions of the 945,219 largest ($d > 0.8$ kb) inter-NIEB regions, where nucleosome density was predicted (Fig. 3d) and observed experimentally (Fig. 3a–c), to be flat, as reflecting the background mutation rates (Methods). Thus, if the NIEBs and flanking nucleosomes were evolving neutrally, we would expect that the ratio of interspecies and intraspecies base changes corrected for “proxy” neutral rates, be robustly $S_{X→Y} = 1$ in the regions of interest as further away from the NIEBs. If instead these regions were under selective pressure, we would expect to see some deviation $S_{X→Y} > 1$ (positive selection) or $S_{X→Y} < 1$ (negative selection) across the ~1 kb regions encompassing the NIEBs and the flanking (on the right and on the left) well-positioned nucleosomes.

As shown in Fig. 5, we obtained convincing signatures of selection in all $W → S$ and $S → W$ substitutions, but

<figure><img src="image_7.png" /><figcaption>**Fig. 5** Evidence of selection in and around the 1,581,256 NIEBs. Ratios $S_{X→Y}$ of background corrected inter- and intraspecies divergence rates plotted against the position from the closest NIEB border (negative distances correspond to loci inside the NIEBs). The pannels correspond to the substitution rates: (a) G → A and C → T; (b) A → G and T → C; (c) C → A and G → T; (d) T → G and A → C; (e) T → A and A → T; (f) C → G and G → C. In each panel the first (resp. second) substitution is represented in dark green (resp. light green). Curves were smoothed over 10 bp windows. The vertical blue lines have the same meaning as in Fig. 1a'–c'. The horizontal dark green (resp. light green) dashed lines, mark the 95 % confidence intervals for $S_{X→Y}$ ratios under the hypothesis of neutral evolution (Methods); the probability for $S_{X→Y}$ to be above (resp. below) these limits is 0.025</figcaption></figure>
the observed patterns of selection look rather complex, with opposite forces of selection in the NIEBs and linkers regions, as compared to the flanking nucleosomes. For instance, a relative excess of the rates of S → W interspecies divergence compared to intraspecies divergence rates was found in the NIEBs and to a lesser extent in the linker regions as a indication of positive selection (Fig. 5a, c). The maximal $S_{X→Y}$ value observed over 10000 bootstrap simulations under the hypothesis of neutral evolution (Methods) is 1.07 (resp. 1.12) for G → A and C → T (resp. C → A and G → T) (Additional file 1: Figure S7). $S_{X→Y}$ values in NIEBs and linker regions are larger that these maximal values underlining the high significance ($P \ll 10^{-4}$) of the positive selection acting on these regions. On the contrary, in the flanking nucleosomes some depletion of the same S → W interspecies substitution relative to intraspecies divergence rates was observed as indicative of purifying selection that progressively vanishes when considering nucleosomes further away from the NIEBs ($P \lesssim 0.025$ on the first nucleosome) (Fig. 5a, c). Interestingly, the opposite patterns of selection were found for W → S substitutions. Negative selection was observed in the NIEBs ($P \ll 10^{-4}$) and

linker regions ($P < 0.025$) and positive selection in the neighbouring nucleosome positions ($P < 0.025$ on the first nucleosome) (Fig. 5b, d). These individual complex selection patterns actually cooperate to the overall maintenance of lower GC composition in NIEBs and linker regions and of high GC composition in the flanking nucleosomes. Note that the selection patterns found around the dyads of the 5,474,320 well-positioned first two flanking nucleosomes in the present study are in perfect agreement with those previously reported by Prendergast and Semple [18] from a similar analysis of 817,774 well-positioned nucleosomes identified in “Schones” in vivo nucleosome occupancy data of which ~ 30% are located within 300 bp of a NIEB border.

As compared to the rather flat selection score profile $S_{X→Y} \simeq 1$ indicating the absence of selection in C → G and G → C substitutions in the flanking nucleosomes and a very weak positive selection in the NIEBs (Fig. 5f), some signature of purifying negative selection was observed in T → A and A → T in the NIEBs ($P \ll 10^{-4}$) and in the linkers ($P < 0.025$) (Fig. 5e). This reflects the presence of highly conserved poly(A) and poly(T) patterns (Fig. 6a) that

<figure><img src="image_8.png" /></figure>
are known to impair nucleosome formation [43, 55, 58] and favour nucleosome disassembly [52] by increasing the DNA wrapping free energy cost around the histone octamer. Interestingly, besides being globally AT-rich, the NIEBs are delimited on each side by these rigid DNA fragments that indeed contribute to delineate their edges. The additional estimation of context-dependent substitution rates proved that $tTt \to tAt$ and $aAa \to aTa$ substitutions are highly anti-selected (Fig. 6b) on the opposite to $sTs \to sAs$ and $sAs \to sTs$ (Fig. 6c), where $s=(c,g)$, confirming that these poly(dA:dT) patterns are evolutionary conserved. Besides this negative selection to disrupt poly(A) and poly(T), we also found some positive selection in $aGa \to aAa$ and $tCt \to tTt$ (Fig. 6d) and also in $aCa \to aAa$ and $tGt \to tTt$ (data not shown) to create these nucleosome inhibitory patterns at the NIEB edges and in the linkers, that slightly exceeds the one observed in $\neg aG\neg a \to \neg aA\neg a$ and $\neg tC\neg t \to \neg tT\neg t$ (Fig. 6e), where $\neg a=(c,t,g)$ and $\neg t=(a,c,g)$. These results provide some clear understanding of the localized excess of positive selection in $S \to W$ substitutions at the NIEB edges and in the linker (Fig. 5a, c) concomitant with a localized excess of negative selection in $W \to S$ substitutions (Fig. 5b, d).

To avoid possible bias in the MK test-inspired analysis induced by spurious low frequency polymorphisms [75], we took into account in our study only SNPs with a minor allele frequency > 1%, which represent 15% of the original sampling set only. However, when considering every SNP, we recovered very similar complex selection patterns with only little differences (Additional file 1: Figure S8), confirming that these signatures of positive and negative selections were maintained even in presence of slightly deleterious mutations.

As a control study, we reproduced our analysis of inter- and intraspecies substitution rates after removing repeated sequences using RepeatMasker (Methods). This had little effect on the broad selection patterns previously observed (Additional file 1: Figure S9). Note that our demonstration is based on the observation that similar oscillating mean GC content profiles are found in chimpanzee and human nearby homologous intrinsic NDRs suggesting a robust nucleosome ordering in these regions (no nucleosome repositioning [76]) across primate evolution.

## Selection maintains optimal GC profile for nucleosome positioning nearby NIEBs

Even though NIEBs were shown to be rather homogeneously distributed in GC-rich and GC-poor isochores (Table 1), one may wonder whether the local GC composition has some effect on the divergence rates observed in the NIEBs as well as in the flanking nucleosomes and linker regions. As shown in Fig. 7a, b, when using as a reference the average GC content computed in the central

regions of the 945,219 largest ($d > 0.8$ kb) inter-NIEB regions, we clearly evidenced significant differences in the GC content profiles obtained for different classes of local reference GC content. The higher the local GC content, the largest the relative GC composition depletion inside NIEBs and linker regions, and in contrast the lowest the relative GC composition excess at the flanking nucleosome positions. In particular, for local GC ≥ 46%, the GC composition inside the flanking well-positioned nucleosomes shows almost no difference with the reference GC content characteristics of background ill-positioned nucleosomes. The computation of the corresponding GC profiles at equilibrium (Methods) confirmed the existence of an optimal GC content at human nucleosome cores (Fig. 7c) consistent with previous study of mononucleotide and 5mer frequencies associated with the set of well-positioned nucleosomes identified in “Schones” in vivo data [18]. In contrast, the equilibrium GC content in NIEBs still depends on the reference GC content so that the difference in equilibrium GC content between NIEBs and inter-NIEB central regions is constant (Fig. 7d). Hence, the evolution of NIEB GC content is relative to its local GC environment.

To explicitly quantify how selection patterns were affected by the local GC content, we re-examined the modes and strengths of selection according to which class of reference GC content the considered inter-NIEB region belongs to (Additional File 1: Figure S10). We observed that the excess of $W \to S$ and in particular of $A \to G$ and $T \to C$ interspecies substitution rates compared to intraspecies rates at flanking nucleosome positions is stronger in a low GC environment and almost vanishing in a high GC environment (Additional file 1: Figure S10B,D). More significantly, the purifying selection observed at these flanking nucleosomes for $S \to W$ is stronger in low GC local environment and almost absent in high GC environment (Additional file 1: Figure S10A, C), as the signature of selection mainly acting to counterbalance the local high mutation rates towards high AT content. As far as the NIEBs and linker regions are concerned, no clear effect on the reference GC composition was observed on the positive and negative selection scores (Additional file 1: Figure S10) strongly suggesting the robustness of these selection patterns as regards to the local compositional environment. Altogether, these results confirm that the NIEBs and flanking nucleosomes correspond to chromatin patterns that have been selected during evolution according to complex selection patterns involving positive and purifying selection. The strengths and directions of selection depend on the local background GC content and combine to maintain a high difference in GC composition (Fig. 7c, d) between the lowest GC composition in the NIEBs and the highest GC composition in the closest flanking nucleosomes.
<figure><img src="image_9.png" /></figure>

**Fig. 7** a Repeat-masked GC content in and around the predicted NIEBs bordering the 945,219 largest (d > 0.8 kb) inter-NIEB regions. These barriers were classified according to the repeat-masked GC content computed in the central part of these large inter-NIEB regions: GC < 0.38 (blue), 0.38 ≤ GC < 0.46 (green) and 0.46 ≤ GC (red). b Repeat-masked GC content normalized by its local repeat-masked reference value. c Same as (a) for GC content at equilibrium. d Equilibrium GC content normalized by its local equilibrium reference value. In (c, d) curves were smoothed over 30 bp windows. The vertical blue lines have the same meaning as in Fig. 1a'–c'

## NIEBs at gene transcription start and termination sites

Since the intrinsic chromatin pattern made of a NIEB and its flanking nucleosomes was found almost equally distributed in genes and in intergenic regions (Table 1, Additional file 1: Table S1), we repeated our analysis by focusing on the 790,285 totally intergenic chromatin patterns. Quite similar patterns of selection were obtained (Additional file 1: Figure S11) demonstrating that the different modes and strengths of selection evidenced in this study (Fig. 5) are generally acting over the entire human autosomes shaping the genome landscape not only relative to functional regulatory sequences, genes and exonic sequences submitted to coding selective pressure, but also in Mb long heterochromatin intergenic regions that were not known to be constrained to maintain adequate nucleosome organization. When focusing our analysis at a smaller scale around human gene TSS, we found some significant depletion of NIEBs coverage in a [−500 bp, 500 bp] window (Additional file 1: Figure S12A). Indeed as originally observed in vitro by Valouev et al. [42] (Additional file 1: Figure S13B), the sequence-dependent physical model confirmed a nucleosome positioning preference on these regulatory sequences and this

for expressed as well as non-expressed genes (Additional file 1: Figure S13A) [77]. As compared to genome average density, a notable excess in nucleosome density is actually programmed to extend over a few kbs upstream of the TSS in the promoter region and also downstream in the gene body. This excess is predicted to be slightly more important for expressed than for non-expressed genes (Additional file 1: Figure S13A), in agreement with the in vitro data (Additional file 1: Figure S13B). These results are quite consistent with the emerging idea that in human, gene promoters, and more generally regulatory sequences, are basally protected by a high nucleosome occupancy [41]. When comparing the intrinsic (Additional file 1: Figure S13A) and in vitro (Additional file 1: Figure S13B) nucleosome density profiles around human gene TSS to the ones observed in vivo by Valouev et al. [42] (Additional file 1: Figure S13C) and Schones et al. [57] (Additional file 1: Figure S13D), we recovered for active genes, a situation previously observed in other organisms such as yeast and fly, with a NDR of ~ 150 bp overlapping the TSS and flanked on both sides by an array of well-positioned and phased nucleosomes [77]. These NDRs and strong nucleosome ordering can be explained as resulting from the
binding of Pol II acting as an “exluding” energy barrier that contributes to the statistical packing of neighbouring nucleosomes [42, 57]. This likely explains part of the ~ 70 % of well-positioned nucleosomes found in “Schones” in vivo data away (> 300 bp) from NIEB borders. This pronounced nucleosome depletion and nearby nucleosome positioning and phasing are actually lost in inactive promoters (Additional file 1: Figure S13C,D) that display rather flat nucleosome density profiles. This suggests that the transition from an inactive to an active promoter involves, in a Pol II-dependent manner, nucleosome sliding and/or nucleosome eviction [28, 42, 57]. Because genes regulation in higher eukaryotes is typically cell-type specific, it thus seems more advantageous to keep the promoters and related regulatory sites intrinsically occupied by nucleosomes unless they need to be used. Note that, as expected by the presence in most protein-coding genes of a poly(A) signal at their 3' end, we observed a remarkable enrichment in NIEBs at TTS (Additional file 1: Figure S12B). As predicted by the sequence-dependent model (Additional file 1: Figure S13A'), and observed in vitro (Additional file 1: Figure S13B'), an intrinsic NDR is programmed at gene TTS with well-positioned flanking nucleosomes, a chromatin pattern which appears to be significantly remodeled in vivo (Additional file 1: Figure S13C, D').

## Mutational pattern at NIEBs cannot be explained by biased gene conversion

Several studies have identified biased gene conversion (BGC) as a neutral process leading to strong variation of W → S mutation rates along human chromosomes via the biased repair of A:C and G:T mismatch through meiotic recombination [64–66]. Thus one may wonder to which extent BGC may underlie some of the base changes observed in our study. Here, we have identified large variations that affect all type of W → S and S → W substitution rates. These variations are not specific of GC-rich, nucleosome-poor regions of the human genome which are known to experience higher recombination rates [78]. We have shown that the selected chromatin patterns involving intrinsic NDRs and flanking nucleosomes are equally distributed in high, medium and low meiotic recombination rate regions (Table 1 and Additional file 1: Table S1). Moreover we found that, independently of their size, double strand break (DSB) hotspots in human individuals [79] are highly depleted in NIEBs (Additional file 1: Figure S14A, B) and consistently that the NIEBs are void of DSB hotspots (Additional file 1: Figure S14C). Altogether, these results demonstrate that the complex selection patterns revealed in this study, their dependence on the local GC content and the highly conserved AT-rich intrinsic NDRs cannot be explained by the neutral BGC process.

## Discussion

Nucleosome organization along the so-called 10 nm chromatin fiber definitely conditions the next steps of DNA compaction into successive higher order structures including condensation into the 30 nm fiber and the formation of chromatin loops, up to a full extent of condensation in metaphase chromosomes [3]. The existence of a (~ 1 kb) chromatin pattern made of an intrinsic NDRs flanked by highly compacted and well-positioned nucleosomes (NRL ~ 150-160 bp) encoded in the sequence and almost homogeneously Poisson-like-distributed along human chromosomes with a mean interdistance ~ 1.5 kb that does not depend much on the GC content of the region, the gene density and the MRT, and that overall cover about 37.5 % of the genome, raises the question of which chromatin structure has been selected during the evolution and if so, to favour or facilitate which function? According to geometrical modelling of the constitutive 30 nm chromatin fiber [59–63], a small nucleosome spacing with a rather short DNA linker size ~ 10-20 bp is likely to impair the condensation of the nucleosomal array into the chromatin fiber, leaving a well-organized chromatin structure open and accessible. But it is well documented that in somatic cells, there are regions of open accessible euchromatin, but also regions of highly repressive heterochromatin. We have reproduced our substitution and selection analysis after classifying our chromatin patterns according to the epigenetic chromatin state they belong to in the Gm12878 differentiated cell type. Following the multivariate analysis of epigenetic marks at 100 kb resolution [72, 80], we considered on the one hand chromatin patterns in gene rich, high GC, early constant replication timing regions (CTRs) enriched in open chromatin marks (e.g. H3K4me3, H3K27ac) and on the other hand chromatin patterns in gene desert, low GC, late CTRs enriched in the repressive heterochromatin mark H3K9me3 associated with the heterochromatin protein 1 (HP1). We found quite similar selection patterns as observed genome-wide (Fig. 5) for both classes of chromatin patterns (Additional file 1: Figure S15A–F), confirming that they likely operate uniformly and robustly along human chromosomes. These results suggest that an open and accessible basal 10 nm chromatin fiber has been selected in human to intrinsically facilitate the epigenetic regulation of nuclear functions in a cell-type-specific manner. This interpretation is strongly supported by the additional observation of quite similar selection patterns in pluripotent H1hesc cell line (Additional file 1: Figure S15A’-F’), which is known to have a highly dynamic and accessible chromatin refractory to both HP1 and polycomb heterochromatin spreading [72]. These results are consistent with the recent visualization at nanoscale resolution of the organization of the nucleosomes in intact nuclei and in single cells in human and mouse [81]. In pluripotent as
well as in differentiated cell types, the reported exper-
imental results argue against the existence of a highly
ordered secondary structure such as the 30 nm chromatin
fiber.

Repetitive mobile elements constitute nearly half of the human genome [67]. During this study, in multiple occasions, we have noticed that the results obtained with the native DNA sequence were slightly different from the ones obtained with the repeat-masked sequence, as exemplified by the native and masked GC content profiles around the NIEBs (Fig. 1c, c'). As compared to the rather smooth two bump masked GC profile that perfectly matches in vitro and in vivo nucleosome occupancy data, the native GC profile displays some striking robust oscillatory internal patterns suggesting the presence of some repeated sequences nearby a non negligible subset of the predicted NIEBs. When systematically investigating the principal families of interdispersed repeats SINES (Alu, MIR) and LINES (Line1, Line2), we found that a lot of (52 %) Alu elements were inserted flanking a NIEB. The Alu retro-transposons have been inserting in primate genomes for the past 65 million years and have reached over one million copies covering ~ 11 % of the human genome [82]. A typical dimeric Alu element is about 300 bp long and is composed of two distinct GC- and CpG-rich monomers, separated by a short A-rich region. Importantly, the 3' end of an Alu element has a longer poly(A) track that plays a critical role in its amplification mechanism [82, 83]. Interestingly, we found that the orientation of the Alu element is strongly dependent on which NIEB side they were inserted. They are mainly sense at the NIEB 5' end and antisense at the NIEB 3' end (Additional file 1: Figure S16A), so that the body of the Alu element is external to the NIEB. The remarkable phasing of the Alu elements at NIEB 5' end (resp. 3' end) results from the matching of the poly(A) tail of the sense (resp. antisense) Alu with the poly(A) (resp. poly(T)) tracks that were shown to define the edges of some of the predicted NIEBs (Fig. 6a) leading to an asymmetric distribution of poly(A) and poly(T) tracks at the NIEB borders (Additional file 1: Figure S16B). However, a majority (61 %) of NIEBs are free on either side of Alu elements and do not result from mechanisms underlying Alu integration [82]. Indeed, for these Alu-free NIEBs, we observed a symmetric enrichment of poly(A) and poly(T) tracks at NIEB borders (Additional file 1: Figure S16C). Hence, NIEBs likely preexists Alu insertion and provide adequate chromatin patterns favorable to Alu integration. A systematic analysis of the spatial distribution of the Alus in relation to NIEBs and related chromatin patterns, should shed a new light on Alu integration and evolution as well as on their role in the chromatin-mediated regulation of gene expression and of the replication spatiotemporal program.

**Conclusions**

Pioneering studies in yeast have shown that rates of diver-
gence are higher in nucleosome cores than in linker DNA
[9, 10]. Various interpretations of these fluctuations in
divergence rates have been proposed, including the pos-
sible consequence of the limited access of DNA wrapped
around the histones to DNA repair proteins [11] and/or
the existence of purifying selection in linker regions to
maintain correct positioning of nucleosomes [9, 12, 13].
AT-rich linker regions relative to GC-rich nucleosomal
sequences were further shown to evolve under compen-
satory dynamics that maintain heterogenous level of AT
content through spatially coupled W(A/T)-losing and W-
gaining substitutions [48]. The recent analysis of well-
positioned nucleosomes in CD4⁺ human T-cells [18] has
revealed the existence of similar substitution patterns in
the linker regions as a part of a complex combination
of positive and negative selection acting to maintain a
favourable base compositions in both linker and core
regions of human nucleosomes. But in regard to the
nucleosomal organization of the chromatin fiber, these
signatures of selection were obtained from a relatively
small subset of well-positioned nucleosomes that cover
less than 5 % of the human genome [18]. In this study,
we have developed a similar analysis of modes of selec-
tion associated with nucleosome occupancy along human
chromosomes, taking as reference chromatin sequences
no longer the ones favouring nucleosome positioning but
instead the ones hindering nucleosome formation as pre-
dicted by a physical model based on sequence-dependent
DNA bending properties. We have shown that these
1,581,256 NIEBs of mean size 153 bp together with their
first two flanking nucleosomes are actually encoded in
the GC content and are widely distributed along the 22
human autosomes in GC-rich and GC-poor isochores,
in early and late replicating regions, in intergenic and
genic regions (Table 1). By comparing rates of inter- and
intraspecies divergences, we have brought evidence of
complex patterns of positive and negative selection that
depend on the local GC content to maintain some optimal
difference in GC composition between the AT-rich intrin-
sic NDRs and the GC-rich well-positioned first flanking
nucleosomes. This set of selected chromatin patterns
(~1 kb) made of well-positioned nucleosomes (~ two on
each side) flanking an intrinsic NDR overall cover about
37.5 % of the human genome, confirming that nucleo-
some positioning and antipositioning are likely to impact
the divergence of many regions traditionally considered
as evolving neutrally, such as intronic and intergenic
DNA. We hope that the generalization of our study to
other (non primate) mammalian and eukaryotic genomes
will provide important clue to understanding genome
evolution and epigenetic regulation in both health
and disease.
## Methods

All this work was done on the 22 autosomes of the hg18 human assembly. For data coming from hg19 assembly, we used liftOver to convert coordinates from hg19 to hg18.

### Physical model

We used the physical model based on sequence-dependent DNA bending properties described in [19, 34, 35]. For the prediction of NIEBs, we started with the set of 1,747,548 NIEBs detected over the 22 human autosomes by the physical model (see Section 3.4 in [77]), from which we discarded the largest ones ($l > 450$ bp) that were shown to have non negligible in vivo nucleosome occupancy. All our analysis was based on the 1,581,256 resulting NIEBs. Note that when we studied the inter-NIEB regions (Fig. 3), we did not take into account the 229 non-fully sequenced inter-NIEBs. For the in silico nucleosome density, we used the same parameters than those detailed in Section 3.2 and 3.3 in [77] for modeling low (resp. high) nucleosome coverage observed in vitro (resp. in vivo). We used liftOver to determine the coordinates of the chimpanzee regions homologous to human NIEBs (liftOver from hg18 to panTro4).

### Genome-wide DNase I sensitivity data

DNase I sensitivity raw tag density data from Gm06990 cells [69] were downloaded from the UCSC FTP site (http://hgdownload.cse.ucsc.edu/goldenPath/hg18/encode DCC/wgEncodeUwDGF/) as a bed file (wgEncodeUw DGFSignalGm06990-bedGraph4.gz).

### Genome-wide nucleosome positioning data

We used genome-wide maps of nucleosome positioning in vivo in resting human CD4⁺ T cells and in vitro obtained by Valouev et al. [42] using micrococcal nuclease (MNase) digestion of chromatin and massive DNA sequencing technique (MNase-seq). Data files were downloaded from Gene Expression Omnibus (GEO, www.ncbi.nlm.nih.gov/geo/) under accession number GSE25133. For every read on the positive strand, we count a score of ‘1 read’ at the corresponding dyad locus, i.e. 70 bp further. We do the same for the negative reads, we count ‘1 read’ 70 bp before the read match.

We also used for comparison, the in vivo genome-wide map of nucleosome positioning in resting human CD4⁺ T cells obtained by Schones et al. [57] using MNase-seq technique. We downloaded the tag coordinate bed files of the resting nucleosomes from dir.nhlbi.nih.gov/papers/lmi/epigenomes/hgtcellnucleosomes.aspx. As for “Valouev” nucleosome data, for every read ‘+’ (resp. ‘-’), we count a score of ‘1 read’ at the corresponding dyad locus, i.e. 70 bp further (resp. before).

For the heat maps of nucleosome density (Fig. 3a-c), we draw, for each NIEB, a horizontal segment around all loci

where some tags were observed, of length equal to twice the number of tags, and then each line corresponds to the average over 200 NIEBs.

Both “Valouev” and “Schones” nucleosome data take into account only uniquely mapped reads (of length ~ 25 bp). To avoid to be biased by the fact that some 25 bp sequences are repeated in the human genome, and therefore cannot be uniquely mapped, we considered, in our analysis, only sites with a mappability equal to 1 [84].

Schones et al. [57] well-positioned nucleosome data as used in [18] were obtained from the weighted sequence model of Reynolds et al. [85].

### Statistical analysis

To test whether the observed mean nucleosome depletion observed in the predicted set of NIEBs actually reflects MNase-seq profile across individual NIEBs, we compared nucleosome tag densities within the barriers and in the two 300 bp flanking windows (Additional file 1: Figure S3). To do so, we used a Z-test. For each NIEB, we computed the mean number $m_1$ of tags per locus within the NIEB and the mean number $m_2$ of tags per locus within the two 300 bp flanking windows as well as their associated standard deviation, $\sigma_1$ and $\sigma_2$ respectively. As explained above, we took into account only the $n_1$ (resp. $n_2$) mappable sites of these two regions. To perform the Z-test, $n_1$ and $n_2$ have to be > 30. Thus, 43,364 NIEBs (2.7 % of our sample) were excluded from our analysis for lack of sufficient mappability. For the remaining 1,537,892 mappable NIEBs and flanking regions, we computed a z-score as:

$$z = \frac{m_2 - m_1}{\sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}} \quad (1)$$

Under the null hypothesis (there is as many tags inside NIEBs as besides them), this Z-test can be approximated by a normal distribution. Hence, there is 5 % of chance to get a z-score > 1.645 (and 5 % of chance to get a z-score < -1.645). For each MNase-seq data set, we classified the 1,537,892 mappable NIEBs into 3 groups: (G1) NIEBs significantly depleted in nucleosome tags at the 5 % confidence level as compared to flanking windows (z-score > 1.645); (G2) NIEBs significantly enriched in nucleosome tags at the 5 % confidence level as compared to flanking windows (z-score < -1.645); and (G3) the remaining NIEBs (-1.645 ≤ z-score ≤ 1.645).

When comparing the mean value of two distributions (of GC content, DNase I HS site coverage, nucleosome coverage, NIEB density or substitution rate) obtained for two classes of region (NIEB vs. genome wide, NIEB vs. flanking regions, different classes of 100 kb windows, gene vs. intergene, ...), we used the one-sided Welch’s unequal variance t-test, where the t-test distribution was approximated by a normal distribution.
### Interspecies divergence rates

We used the point-mutation data obtained by Chen et al. [6]. Nucleotide substitutions were tabulated in the human lineage since its divergence from chimpanzee using both the orangutan and macaque as outgroups. To minimize the effects of alignment artifacts, only isolated substitutions defined as those flanked by sites that are identical in the four species were tabulated. Sequences were divided into CpG and non-CpG sites. CpG sites were defined as the sites having the following human/chimpanzee/orangutan/macaque pattern: NG/CG/CG/CG or CG/NG/CG/CG or CN/CG/CG/CG or CG/CN/CG/CG, where N is any nucleotide. CpG sites were excluded for the computation of C → W/S substitution rates. Substitution rates were calculated by dividing the number of substitution events of the appropriate type by the number of potentially mutable sites that meet the same criteria. Since the divergence between these four catarrhini species is small, possible multiple substitutions were ignored (using two outgroups instead of one lowers the amount of multiple mutations, in particular at CpG sites, and preferentially eliminates sites that are not ancestral to human and chimpanzee).

### GC composition at equilibrium

We computed the GC profile towards which the sequences are evolving as the result of the observed substitution profiles without taking into account neighbor-dependency, including CpG substitutions.

### Intraspecies divergence rates

We used the SNP data from The 1000 Genomes Project [73] (http://www.ncbi.nlm.nih.gov/variation/tools/1000genomes/). First, we dowloaded the ASN.1 flat files from ftp://ftp.ncbi.nlm.nih.gov/snp/organisms/human_9606_b142_GRCh37p13/ASN1_flat/ and we did a lift-Over, from hg19 to hg18, of the SNP coordinates. Then, the ancestral base at any position was assumed to be any observed allele (including that observed in the reference genome) that matched the corresponding base in the ancestral genome (the one defined above for the definition of the interspecies divergence rates), or if this latter is not defined, the corresponding base in the chimp genome. Sites where the ancestral base could not be determined were excluded. As for interspecies rates, sequences were divided into CpG and non-CpG sites and CpG sites were excluded for the computation of C → W/S substitution rates.

Note that because MK test-inspired analyses were shown to be potentially skewed by the presence of slightly deleterious mutations (disproportionately segregate at low frequencies) [75], all this work was performed by taking into account only alleles observed in more than 1 % of the genomes, which represent only 15 % of the total amount of

SNPs described in the files. However, this had little effect on the broad patterns seen in our analysis (Additional file 1: Figure S8 as compared to Fig. 5).

### Rates of selection

The rate of base changes at each position from the NIEB border was measured by dividing the observed number of base changes by the total number of matching ancestral bases at each position [18]:

$$d\text{Inter}_{\chi \to \gamma} = \frac{\text{InterDiff}_{\chi \to \gamma}}{\text{ancestralBase}_{\chi}}, \quad (2)$$

$$d\text{Intra}_{\chi \to \gamma} = \frac{\text{IntraDiff}_{\chi \to \gamma}}{\text{ancestralBase}_{\chi}}. \quad (3)$$

Equations 2 and 3 calculate the rates of base change at given positions from the NIEB border. $\chi$ and $\gamma$ correspond to the bases before and after the specific change, respectively, $\chi$ being the ancestral base and $\gamma$ being the base observed in the human lineage. $\text{InterDiff}_{\chi \to \gamma}$ and $\text{IntraDiff}_{\chi \to \gamma}$ are the total number of relevant interspecies and intraspecies changes observed at the position of interest relative to the NIEB border. $\text{ancestralBase}_{\chi}$ is the number of corresponding ancestral bases observed at the same position.

Positions relative to the NIEB border where interspecies rates of change showed unusual deviations from intraspecies rates were identified by first correcting each value of $d\text{Inter}_{\chi \to \gamma}$ and $d\text{Intra}_{\chi \to \gamma}$ for mean/background rates of divergence. Mean/background rates of change, backgroundInter and backgroundIntra, were estimated by averaging over the central positions in the middle (±300 bp) of the largest inter-NIEBs ($l > 800$ bp). Those regions look like having rather flat inter- and intraspecies divergence rates (Fig. 4). (Note that for the Additional file 1: Figure S10, associated to the Fig. 7, three different backgrounds were computed for each base changes $\chi \to \gamma$: one for each class of local GC content.) The rate of interspecies base changes observed across all NIEBs and flanking regions at each position was then divided by the corresponding rate of intraspecies change to provide an indication of selection:

$$S_{\chi \to \gamma} = \frac{d\text{Inter}_{\chi \to \gamma}/\text{backgroundInter}_{\chi \to \gamma}}{d\text{Intra}_{\chi \to \gamma}/\text{backgroundIntra}_{\chi \to \gamma}}. \quad (4)$$

The ratio of background corrected inter- and intraspecies divergence rates was calculated for each base change and at each position from the NIEB border. $S_{\chi \to \gamma} > 1$ indicates an excess of interspecies change and positive selection relative to mean rates, and values $S_{\chi \to \gamma} < 1$ indicates an excess of intraspecies changes and negative selection.

In order to determine the $S_{\chi \to \gamma}$ values clearly distinguishing a hypothesis of selection from a hypothesis of
neutral drift, we determined the expected $S_{X→Y}$ distributions under the neutral hypothesis as the distributions obtained for the middle of the largest inter-NIEBs ($l > 800$ bp). Given the 10 bp resolution adopted in this work, we randomly selected a 10 bp window in each of the large inter-NIEBs to compute one realization of $S_{X→Y}$ values (Eq. (4)) under the neutral hypothesis. Repeating this operation 10000 times, we obtained the expected $S_{X→Y}$ histograms (Additional file 1: Figure S7). We derived from these distributions 95 % confidence intervals for $S_{X→Y}$ under the neutral hypothesis. Note that since large inter-NIEBs are only a fraction of the total number of regions averaged in Fig. 5 closed to and within the the NIEBs, these confidence intervals are conservative estimates. As expected, these intervals nicely bound the $S_{X→Y}$ fluctuations away from the NIEBs (Fig. 5) illustrating that they can be used to visually assert the expected $S_{X→Y}$ fluctuations without performing an explicit simulation in Fig. 6 and Additional file 1: Figures S8–S11 and S15). Note that since each 10 bp non-overlapping window along an observed S ratio curve (Fig. 5) is an independent observation of a S ratio value, it is justified to compute the P-value of these values using the expected $S_{X→Y}$ histograms (Additional file 1: Figure S7) without correcting for multiple tests.

### Meiotic recombination

Crossover rate data of the human genome were retrieved from the International HapMap Project (http://www.hapmap.org) [86]. The crossover rate for a given 100 kb window (Table 1) was computed by taking the two most distant markers belonging to this window and dividing their distance in the genetic map in centimorgans (CM) by their distance in the sequence map in megabases (Mb). Note that only the 100 kb windows with at least 5 markers were considered. We used the genome-wide recombination initiation maps of Pratto et al. [79] (Additional file 1: Figure S14). We downloaded the double strand break (DSB) hotspots in human individuals data (the GSE59836 peak data file) from http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM1447330.

### Mean replication timing data

Mean replication timing (MRT) was derived from RepliSeq data [6, 68] for 100 kb non-overlapping windows in hg18 coordinates. The MRT profile of Gm06990 was obtained from the authors [87].

### Gene expression data

For the Additional file 1: Figure S13, expression data were retrieved from the Genome Browser of the University of California Santa Cruz (UCSC). We downloaded expression values from the release 2 of Caltech RNA-seq track (ENCODE project at UCSC, http://

hgdownload.cse.ucsc.edu/goldenPath/hg18/encodeDCC/wgEncodeCalthechRnaSeq/).

Expression for one gene is given in reads per kb of exon model per million mapped reads (RPKM) [88]. RPKM is defined as:

$$R = \frac{10^9 C}{NL}, \qquad (5)$$

where C is the number of mappable reads that fall into gene exons (union of exons for genes with alternative splicing), N is the total number of mappable reads in the experiment, and L is the total length of the exons in base pairs. We associated 17 872 genes with a valid RPKM value in Gm12878.

To define the genes and the intergenes regions in Table 1, as well as the totally intergenic 100 kb windows, we determined the genome coordinates of each gene (labeled by its RefSeq identifier) by using RefSeq Genes track. For genes associated with more than one splicing variant, we merged exons coordinates by taking their union. We obtained a table of 23 329 genes. The intergenes coordinates correspond to intervals between genes from which we remove 2 kb from each side (to be sure to be in an 'intergenic' region). For the Additional file 1: Figure S12, we used the transcription start site (TSS) associated to these genes (placed at the beginning of the first exon).

### Euchromatin and heterochromatin regions

Data for the transcriptionally active, early replicating euchromatin regions (C1+C2 in Gm12878 and EC1+EC2 in H1hesc) and transcriptionally inactive, late replicating heterochromatin regions (C3+C4 in Gm12878 and EC3+EC4 in H1hesc), used in Additional file 1: Figure S15, were obtained from the authors [72, 89].

### Repeat sequence annotations

We retrieved repeat sequence annotations from the UCSC Genome Browser (rmsk table) obtained using Repeat-Masker [90].

### Additional file

Additional file 1: Supplementary Online Material. Additional data file 1 contains supporting Figures S1 to S16 and supporting Tables S1 and S2 (PDF). (PDF 1802 kb)

### Abbreviations

BGC, biased gene conversion; DSB, double strand break; MRT, mean replication timing; MK, McDonald-Kreitman; NDR, nucleosome-depleted region; NIEB, nucleosome-inhibiting energy barrier; NRL, nucleosome repeat-length; RPKM, reads per kb of exon model per million mapped reads; SNP, single nucleotide polymorphisms; TSS, transcription start site; TTS, transcription termination site

### Acknowledgements

We would like to thank C.-L. Chen for providing assistance with human/chimpanzee mutational data. We are very grateful to L. Duret and C. Vaillant for helpful discussions.
**Funding**

This work was supported by the Institut National de la Santé et de la Recherche Médicale (Plan Cancer 2012 01-84862). The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

**Availability of data and materials**

The NIEBs predictions supporting the conclusions of this article can be obtained from the authors.

**Authors' contributions**

All authors designed the study. GD performed the analysis. BA and FA contributed to elaborate analysis tools. BA, FA and AA supervised the study. AA wrote the manuscript with contributions from all co-authors. All authors read and approved the final manuscript.

**Competing interests**

All authors declare that they have no competing interests.

**Consent for publication**

Not applicable.

**Ethics approval and consent to participate**

Not applicable.

Received: 14 December 2015 Accepted: 4 July 2016

Published online: 29 July 2016

**References**

1. Arneodo A, Vaillant C, Audit B, Argoul F, d'Aubenton-Carafa Y, Thermes C. Multi-scale coding of genomic information: From DNA sequence to genome structure and function. Phys Rep. 2011;498:45-188.

2. Parker SCJ, Tullius TD. DNA shape, genetic codes, and evolution. Curr Opin Struct Biol. 2011;21(3):342-7. doi:10.1016/j.sbi.2011.03.002.

3. Alberts B, Bray D, Lewis J, Raff M, Roberts K, Watson JD. Molecular Biology of the Cell, 4th ed. New-York: Garland Publishing; 2002.

4. Graur D, Li WH. Fundamentals of Molecular Evolution. Sunderland: Sinauer Associates; 1999.

5. Stamatoyannopoulos JA, Adzhubei I, Thurman RE, Kryukov GV, Mirkin SM, Sunyaev SR. Human mutation rate associated with DNA replication timing. Nat Genet. 2009;41(4):393-5. doi:10.1038/ng.363.

6. Chen CL, Rappailles A, Duquenne L, Huvet M, Guilbaud G, Farinelli L, Audit B, d'Aubenton-Carafa Y, Arneodo A, Hyrien O, Thermes C. Impact of replication timing on non-CpG and CpG substitution rates in mammalian genomes. Genome Res. 2010;20:447-57. doi:10.1101/gr.098947.109.

7. Chen CL, Duquenne L, Audit B, Guilbaud G, Rappailles A, Baker A, Huvet M, d'Aubenton-Carafa Y, Hyrien O, Arneodo A, Thermes C. Replication-associated mutational asymmetry in the human genome. Mol Biol Evol. 2011;28(8):2327-337. doi:10.1093/molbev/msr056.

8. Makova KD, Hardison RC. The effects of chromatin organization on variation in mutation rates in the genome. Nat Rev Genet. 2015;16(4): 213-23. doi:10.1038/nrg3890.

9. Warnecke T, Batada NN, Hurst LD. The impact of the nucleosome code on protein-coding sequence evolution in yeast. PLoS Genet. 2008;4(11): 1000250. doi:10.1371/journal.pgen.1000250.

10. Washietl S, Machné R, Goldman N. Evolutionary footprints of nucleosome positions in yeast. Trends Genet. 2008;24(12):583-7. doi:10.1016/j.tig.2008.09.003.

11. Chen X, Chen Z, Chen H, Su Z, Yang J, Lin F, Shi S, He X. Nucleosomes suppress spontaneous mutations base-specifically in eukaryotes. Science. 2012;335(6073):1235-1238. doi:10.1126/science.1217580.

12. Babbitt GA, Kim Y. Inferring natural selection on fine-scale chromatin organization in yeast. Mol Biol Evol. 2008;25(8):1714-1727. doi:10.1093/molbev/msn127.

13. Babbitt GA, Cotter CR. Functional conservation of nucleosome formation selectively biases presumably neutral molecular variation in yeast genomes. Genome Biol Evol. 2011;3:15-22. doi:10.1093/gbe/evq081.

14. Langley SA, Karpen GH, Langley CH. Nucleosomes shape DNA polymorphism and divergence. PLoS Genet. 2014;10(7):1004457. doi:10.1371/journal.pgen.1004457.

15. Sasaki S, Mello CC, Shimada A, Nakatani Y, Hashimoto SI, Ogawa M, Matsushima K, Gu SG, Kasahara M, Ahsan B, Sasaki A, Saito T, Suzuki Y, Sugano S, Kohara Y, Takeda H, Fire A, Morishita S. Chromatin-associated periodicity in genetic variation downstream of transcriptional start sites. Science. 2009;323(5912):401-4. doi:10.1126/science.1163183.

16. Ying H, Epps J, Williams R, Huttley G. Evidence that localized variation in primate sequence divergence arises from an influence of nucleosome placement on DNA repair. Mol Biol Evol. 2010;27(3):637-49. doi:10.1093/molbev/msp253.

17. Tolstorukov MY, Volfovsky N, Stephens RM, Park PJ. Impact of chromatin structure on sequence variability in the human genome. Nat Struct Mol Biol. 2011;18(4):510-5. doi:10.1038/nsmb.2012.

18. Prendergast JGD, Semple CAM. Widespread signatures of recent selection linked to nucleosome positioning in the human lineage. Genome Res. 2011;21(11):1777-1787. doi:10.1101/gr.122275.111.

19. Chevereau G, Arneodo A, Vaillant C. Influence of the genomic sequence on the primary structure of chromatin. Front Life Sci. 2011;5(1-2):29-68. doi:10.1080/21553769.2012.708882.

20. Travers AA, Vaillant C, Arneodo A, Muskhelishvili G. DNA structure, nucleosome placement and chromatin remodelling: a perspective. Biochem Soc Trans. 2012;40(2):335-40. doi:10.1042/BST20110757.

21. Struhl K, Segal E. Determinants of nucleosome positioning. Nat Struct Mol Biol. 2013;20(3):267-73. doi:10.1038/nsmb.2506.

22. Hughes AL, Rando OJ. Mechanisms underlying nucleosome positioning in vivo. Annu Rev Biophys. 2014;43:41-63. doi:10.1146/annurev-biophys-051013-023114.

23. Satchwell SC, Drew HR, Travers AA. Sequence periodicities in chicken nucleosome core DNA. J Mol Biol. 1986;191:659-75.

24. Ioshikhes I, Bolshoy A, Derenshteyn K, Borodovsky M, Trifonov EN. Nucleosome DNA sequence pattern revealed by multiple alignment of experimentally mapped sequences. J Mol Biol. 1996;262(2):129-39. doi:10.1006/jmbi.1996.0503.

25. Lowary PT, Widom J. New DNA sequence rules for high affinity binding to histone octamer and sequence-directed nucleosome positioning. J Mol Biol. 1998;276(1):19-42. doi:10.1006/jmbi.1997.1494.

26. Segal E, Fondufe-Mittendorf Y, Chen L, Thäström A, Field Y, Moore IK, Wang J-PZ, Widom J. A genomic code for nucleosome positioning. Nature. 2006;442(7104):772-8. doi:10.1038/nature04979.

27. Peckham HE, Thurman RE, Fu Y, Stamatoyannopoulos JA, Noble WS, Struhl K, Weng Z. Nucleosome positioning signals in genomic DNA. Genome Res. 2007;17(8):1170-1177. doi:10.1101/gr.6101007.

28. Zhang Y, Moqtaderi Z, Rattner BP, Euskirchen G, Snyder M, Kadonaga JT, Liu XS, Struhl K. Intrinsic histone-DNA interactions are not the major determinant of nucleosome positions in vivo. Nat Struct Mol Biol. 2009;16(8):847-52. doi:10.1038/nsmb.1636.

29. Bettecken T, Frenkel ZM, Trifonov EN. Human nucleosomes: special role of CG dinucleotides and Alu-nucleosomes. BMC Genomics. 2011;12(1): 273. doi:10.1186/1471-2164-12-273.

30. Kornberg RD, Stryer L. Statistical distributions of nucleosomes: nonrandom locations by a stochastic mechanism. Nucleic Acids Res. 1988;16(14A):6677-690.

31. Mavrich TN, Ioshikhes IP, Venters BJ, Jiang C, Tomsho LP, Qi J, Schuster SC, Albert I, Pugh BF. A barrier nucleosome model for statistical positioning of nucleosomes throughout the yeast genome. Genome Res. 2008;18(7):1073-1083. doi:10.1101/gr.078261.108.

32. Vaillant C, Palmeira L, Chevereau G, Audit B, d'Aubenton-Carafa Y, Thermes C, Arneodo A. A novel strategy of transcription regulation by intragenic nucleosome ordering. Genome Res. 2010;20(1):59-67. doi:10.1101/gr.096644.109.

33. Möbius W, Gerland U. Quantitative test of the barrier nucleosome model for statistical positioning of nucleosomes up- and downstream of transcription start sites. PLoS Comput Biol. 2010;6(8):1000891. doi:10.1371/journal.pcbi.1000891.

34. Vaillant C, Audit B, Arneodo A. Experiments confirm the influence of genome long-range correlations on nucleosome positioning. Phys Rev Lett. 2007;99(21):218103. doi:10.1103/PhysRevLett.99.218103.

35. Milani P, Chevereau G, Vaillant C, Audit B, Haftek-Terreau Z, Marilley M, Bouvet P, Argoul F, Arneodo A. Nucleosome positioning by genomic excluding-energy barriers. Proc Natl Acad Sci USA. 2009;106(52): 22257-2262. doi:10.1073/pnas.0909511106.
36. Stroud H, Otero S, Desvoyes B, Ramírez-Parra E, Jacobsen SE, Gutierrez C. Genome-wide analysis of histone H3.1 and H3.3 variants in *Arabidopsis thaliana*. Proc Natl Acad Sci USA. 2012;109(14):5370–375. doi:10.1073/pnas.1203145109.

37. Cayrou C, Coulombe P, Puy A, Rialle S, Kaplan N, Segal E, Méchali M. New insights into replication origin characteristics in metazoans. Cell Cycle. 2012;11(4):658–67. doi:10.4161/cc.11.4.19097.

38. Teif VB, Vainshtein Y, Caudron-Herger M, Mallm JP, Marth C, Höfer T, Rippe K. Genome-wide nucleosome positioning during embryonic stem cell development. Nat Struct Mol Biol. 2012;19(11):1185–1192. doi:10.1038/nsmb.2419.

39. Fenouil R, Cauchy P, Koch F, Descostes N, Cabeza JZ, Innocenti C, Ferrier P, Spicuglia S, Gut M, Gut I, Araújo JC. CpG islands and GC content dictate nucleosome depletion in a transcription-independent manner at mammalian promoters. Genome Res. 2012;22(12):2399–408. doi:10.1101/gr.138776.112.

40. Lombraña R, Almeida R, Revuelta I, Madeira S, Herranz G, Saiz N, Bastolla U, Gómez M. High-resolution analysis of DNA synthesis start sites and nucleosome architecture at efficient mammalian replication origins. EMBO J. 2013;32(19):2631–644. doi:10.1038/emboj.2013.195.

41. Tillo D, Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Field Y, Lieb JD, Widom J, Segal E, Hughes TR. High nucleosome occupancy is encoded at human regulatory sequences. PLoS One. 2010;5(2):9129. doi:10.1371/journal.pone.0009129.

42. Valouev A, Johnson SM, Boyd SD, Smith CL, Fire AZ, Sidow A. Determinants of nucleosome organization in primary human cells. Nature. 2011;474(7352):516–20. doi:10.1038/nature10002.

43. 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–30. doi:10.1126/science.1112178.

44. 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. doi:10.1038/ng1217.

45. Field Y, Kaplan N, Fondufe-Mittendorf Y, Moore IK, Sharon E, Lubling Y, Widom J, Segal E. Distinct modes of regulation by chromatin encoded through nucleosome positioning signals. PLoS Comput Biol. 2008;4(11): 1000216. doi:10.1371/journal.pcbi.1000216.

46. Eaton ML, Galani K, Kang S, Bell SP, MacAlpine DM. Conserved nucleosome positioning defines replication origins. Genes Dev. 2010;24(8):748–53. doi:10.1101/gad.1913210.

47. Berbenetz NM, Nislow C, Brown GW. Diversity of eukaryotic DNA replication origins revealed by genome-wide analysis of chromatin structure. PLoS Genet. 2010;6(9):1001092. doi:10.1371/journal.pgen.1001092.

48. Kenigsberg E, Bar A, Segal E, Tanay A. Widespread compensatory evolution conserves DNA-encoded nucleosome organization in yeast. PLoS Comput Biol. 2010;6(12):1001039. doi:10.1371/journal.pcbi.1001039.

49. Tsankov A, Yanagisawa Y, Rhind N, Regev A, Rando OJ. Evolutionary divergence of intrinsic and trans-regulated nucleosome positioning sequences reveals plastic rules for chromatin organization. Genome Res. 2011;21(11):1851–1862. doi:10.1101/gr.122267.111.

50. Tsui K, Dubuis S, Gebbia M, Morse RH, Barkai N, Tirosh I, Nislow C. Evolution of nucleosome occupancy: conservation of global properties and divergence of gene-specific patterns. Mol Cell Biol. 2011;31(21): 4348–355. doi:10.1128/MCB.05276-11.

51. Valouev A, Ichikawa J, Tonthat T, Stuart J, Ranade S, Peckham H, Zeng K, Malek JA, Costa G, McKernan K, Sidow A, Fire A, Johnson SM. A high-resolution, nucleosome position map of *C. elegans* reveals a lack of universal sequence-dictated positioning. Genome Res. 2008;18(7): 1051–1063. doi:10.1101/gr.076463.108.

52. Iyer V, Struhl K. Poly(dA:dT) is a ubiquitous promoter element that stimulates transcription via its intrinsic DNA structure. EMBO J. 1995;14(11):2570–579.

53. Miele V, Vaillant C, d'Aubenton-Carafa Y, Thermes C, Grange T. DNA physical properties determine nucleosome occupancy from yeast to fly. Nucleic Acids Res. 2008;36(11):3746–756. doi:10.1093/nar/gkn262.

54. Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Tillo D, Field Y, LeProust EM, Hughes TR, Lieb JD, Widom J, Segal E. The DNA-encoded nucleosome organization of a eukaryotic genome. Nature. 2009;458(7236):362–6. doi:10.1038/nature07667.

55. Segal E, Widom J. Poly(dA:dT) tracts: major determinants of nucleosome organization. Curr Opin Struct Biol. 2009;19(1):65–71. doi:10.1016/j.sbi.2009.01.004.

56. Ozsolak F, Song JS, Liu XS, Fisher DE. High-throughput mapping of the chromatin structure of human promoters. Nat Biotechnol. 2007;25(2): 244–8. doi:10.1038/nbt1279.

57. 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–98. doi:10.1016/j.cell.2008.02.022.

58. Tillo D, Hughes TR. G+C content dominates intrinsic nucleosome occupancy. BMC Bioinforma. 2009;10:442. doi:10.1186/1471-2105-10-442.

59. Schiessel H. Topical Review: The physics of chromatin. J Phys Condens Matter. 2003;15:699–774.

60. Lesne A, Victor JM. Chromatin fiber functional organization: some plausible models. Eur Phys J E. 2006;19(3):279–90. doi:10.1140/epje/i2005-10050-6.

61. Kepper N, Foethke D, Stehr R, Wedemann G, Rippe K. Nucleosome geometry and internucleosomal interactions control the chromatin fiber conformation. Biophys J. 2008;95(8):3692–705. doi:10.1529/biophysj.107.121079.

62. Diesinger PM, Heermann DW. Depletion effects massively change chromatin properties and influence genome folding. Biophys J. 2009;97(8):2146–153. doi:16/j.bpj.2009.06.057.

63. Beshnova DA, Cherstvy AG, Vainshtein Y, Teif VB. Regulation of the nucleosome repeat length in vivo by the DNA sequence, protein concentrations and long-range interactions. PLoS Comput Biol. 2014;10(7):1003698. doi:10.1371/journal.pcbi.1003698.

64. Eyre-Walker A, Hurst LD. The evolution of isochores. Nat Rev Genet. 2001;2(7):549–55. doi:10.1038/35080577.

65. Duret L, Eyre-Walker A, Galtier N. A new perspective on isochore evolution. Gene. 2006;385:71–4. doi:10.1016/j.gene.2006.04.030.

66. Duret L, Arndt PF. The impact of recombination on nucleotide substitutions in the human genome. PLoS Genet. 2008;4(5):1000071. doi:10.1371/journal.pgen.1000071.

67. Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, Funke R, Gage D, Harris K, Heaford A, Howland J, Kann L, Lehoczky J, LeVine R, McEwan P, McKernan K, Meldrim J, Mesirov JP, Miranda C, Morris W, Naylor J, Raymond C, Rossetti M, Santos R, Sheridan A, Sougnez C, Stange-Thomann N, Stojanovic N, Subramanian A, Wyman D, Rogers J, Sulston J, Ainscough R, Beck S, Bentley D, Burton J, Clee C, Carter N, Coulson A, Deadman R, Deloukas P, Dunham A, Dunham I, Durbin R, French L, Grafham D, Gregory S, Hubbard T, Humphray S, Hunt A, Jones M, Lloyd C, McMurray A, Matthews L, Mercer S, Milne S, Mullikin JC, Mungall A, Plumb R, Ross M, Shownkeen R, Sims S, Waterston RH, Wilson RK, Hillier LW, McPherson JD, Marra MA, Mardis ER, Fulton LA, Chinwalla AT, Pepin KH, Gish WR, Chissoe SL, Wendl MC, Delehaunty KD, Miner TL, Delehaunty A, Kramer JB, Cook LL, Fulton RS, Johnson DL, Minx PJ, Clifton SW, Hawkins T, Branscomb E, Predki P, Richardson P, Wenning S, Slezak T, Doggett N, Cheng JF, Olsen A, Lucas S, Elkin C, Uberbacher E, Frazier M, Gibbs RA, Muzny DM, Scherer SE, Bouck JB, Sodergren EJ, Worley KC, Rives CM, Gorrell JH, Metzker ML, Naylor SL, Kucherlapati RS, Nelson DL, Weinstock GM, Sakaki Y, Fujiyama A, Hattori M, Yada T, Toyoda A, Itoh T, Kawagoe C, Watanabe H, Totoki Y, Taylor T, Weissenbach J, Heilig R, Saurin W, Artiguenave F, Brottier P, Bruls T, Pelletier E, Robert C, Wincker P, Smith DR, Doucette-Stamm L, Rubenfield M, Weinstock K, Lee HM, Dubois J, Rosenthal A, Platzer M, Nyakatura G, Taudien S, Rump A, Yang H, Yu J, Wang J, Huang G, Gu J, Hood L, Rowen L, Madan A, Qin S, Davis RW, Federspiel NA, Abola AP, Proctor MJ, Myers RM, Schmutz J, Dickson M, Grimwood J, Cox DR, Olson MV, Kaul R, Raymond C, Shimizu N, Kawasaki K, Minoshima S, Evans GA, Athanasiou M, Schultz R, Roe BA, Chen F, Pan H, Ramser J, Lehrach H, Reinhardt R, McCombie WR, de la Bastide M, Dedhia N, Blöcker H, Hornischer K, Nordsiek G, Agarwala R, Aravind L, Bailey JA, Bateman A, Batzoglou S, Birney E, Bork P, Brown DG, Burge CB, Cerutti L, Chen HC, Church D, Clamp M, Copley RR, Doerks T, Eddy SR, Eichler EE, Furey TS, Galagan J, Gilbert JG, Harmon C, Hayashizaki Y, Haussler D, Hermjakob H, Hokamp K, Jang W, Johnson LS, Jones TA, Kasif S, Kasprzyk A, Kennedy S, Kent WJ, Kitts P, Koonin EV, Korf I, Kulp D, Lancet D, Lowe TM, McLysaght A, Mikkelsen T, Moran JV, Mulder N, Pollara VJ, Ponting CP, Schuler G, Schultz J, Slater G, Smit AF, Stupka E, Szustakowski J, Thierry-Mieg D, Thierry-Mieg D.
J, Wagner L, Wallis J, Wheeler R, Williams A, Wolf YI, Wolfe KH, Yang SP, Yeh RF, Collins F, Guyer MS, Peterson J, Felsenfeld A, Wetterstrand KA, Patrinos A, Morgan MJ, de Jong P, Catanese JJ, Osoegawa K, Shizuya H, Choi S, Chen YJ, Szustakowki J, Consortium IHGS. Initial sequencing and analysis of the human genome. Nature. 2001;409(6822):860–921.

68. Hansen RS, Thomas S, Sandstrom R, Canfield TK, Thurman RE, Weaver M, Dorschner MO, Gartler SM, Stamatoyannopoulos JA. Sequencing newly replicated DNA reveals widespread plasticity in human replication timing. Proc Natl Acad Sci USA. 2010;107(1):139–44.
doi:10.1073/pnas.0912402107.

69. Sabo PJ, Kuehn MS, Thurman R, Johnson BE, Johnson EM, Cao H, Yu M, Rosenzweig E, Goldy J, Haydock A, Weaver M, Shafer A, Lee K, Neri F, Humbert R, Singer MA, Richmond TA, Dorschner MO, McArthur M, Hawrylycz M, Green RD, Navas PA, Noble WS, Stamatoyannopoulos JA. Genome-scale mapping of DNase I sensitivity in vivo using tiling DNA microarrays. Nat Methods. 2006;3(7):511–8. doi:10.1038/nmeth890.

70. Audit B, Zaghloul L, Vaillant C, Chevereau G, d'Aubenton-Carafa Y, Thermes C, Arneodo A. Open chromatin encoded in DNA sequence is the signature of "master" replication origins in human cells. Nucleic Acids Res. 2009;37(18):6064–075. doi:10.1093/nar/gkp631.

71. Gindin Y, Valenzuela MS, Aladjem MI, Meltzer PS, Bilke S. A chromatin structure-based model accurately predicts DNA replication timing in human cells. Mol Syst Biol. 2014;10:722.

72. Julienne H, Audit B, Arneodo A. Embryonic stem cell specific "master" replication origins at the heart of the loss of pluripotency. PLoS Comput Biol. 2015;11(2):1003969. doi:10.1371/journal.pcbi.1003969.

73. The 1000 Genomes Project Consortium, Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, Gibbs RA, Hurles ME, McVean GA. A map of human genome variation from population-scale sequencing. Nature. 2010;467(7319):1061–1073. doi:10.1038/nature09534.

74. McDonald JH, Kreitman M. Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991;351(6328):652–4. doi:10.1038/351652a0.

75. Charlesworth J, Eyre-Walker A. The McDonald-Kreitman test and slightly deleterious mutations. Mol Biol Evol. 2008;25(6):1007–1015.
doi:10.1093/molbev/msn005.

76. Warnecke T, Becker EA, Facciotti MT, Nislow C, Lehner B. Conserved substitution patterns around nucleosome footprints in eukaryotes and archaea derive from frequent nucleosome repositioning through evolution. PLoS Comput Biol. 2013;9(11):1003373.
doi:10.1371/journal.pcbi.1003373.

77. Drillon G, Audit B, Argoul F, Arneodo A. Ubiquitous human 'master' origins of replication are encoded in the DNA sequence via a local enrichment in nucleosome excluding energy barriers. J Phys Condens Matter. 2015;27(6):064102. doi:10.1088/0953-8984/27/6/064102.

78. Duret L, Galtier N. Biased gene conversion and the evolution of mammalian genomic landscapes. Annu Rev Genomics Hum Genet. 2009;10:285–311. doi:10.1146/annurev-genom-082908-150001.

79. Pratto F, Brick K, Khil P, Smagulova F, Petukhova GV, Camerini-Otero RD. DNA recombination. Recombination initiation maps of individual human genomes. Science. 2014;346(6211):1256442.
doi:10.1126/science.1256442.

80. Julienne H, Zoufir A, Audit B, Arneodo A. Human genome replication proceeds through four chromatin states. PLoS Comput Biol. 2013;9(10): 1003233. doi:10.1371/journal.pcbi.1003233.

81. Ricci MA, Manzo C, García-Parajo MF, Lakadamyali M, Cosma MP. Chromatin fibers are formed by heterogeneous groups of nucleosomes in vivo. Cell. 2015;160(6):1145–1158. doi:10.1016/j.cell.2015.01.054.

82. Deininger P. Alu elements: know the SINEs. Genome Biol. 2011;12(12):236.
doi:10.1186/gb-2011-12-12-236.

83. Mighell AJ, Markham AF, Robinson PA. Alu sequences. FEBS Lett. 1997;417(1):1–5.

84. Derrien T, Estellé J, Sola SM, Knowles DG, Raineri E, Guigó R, Ribeca P. Fast computation and applications of genome mappability. PLoS One. 2012;7(1):30377. doi:10.1371/journal.pone.0030377.

85. Reynolds SM, Bilmes JA, Noble WS. Learning a weighted sequence model of the nucleosome core and linker yields more accurate predictions in *Saccharomyces cerevisiae* and *Homo sapiens*. PLoS Comput Biol. 2010;6(7):1000834. doi:10.1371/journal.pcbi.1000834.

86. Frazer KA, Eskin E, Kang HM, Bogue MA, Hinds DA, Beilharz EJ, Gupta RV, Montgomery J, Morenzoni MM, Nilsen GB, Pethiyagoda CL, Stuve LL, Johnson FM, Daly MJ, Wade CM, Cox DR. A sequence-based

variation map of 8.27 million SNPs in inbred mouse strains. Nature.
2007;448(7157):1050–1053. doi:10.1038/nature06067.

87. Baker A, Audit B, Chen CL, Moindrot B, Leleu A, Guilbaud G, Rappailles A, Vaillant C, Goldar A, Mongelard F, d'Aubenton-Carafa Y, Hyrien O, Thermes C, Arneodo A. Replication fork polarity gradients revealed by megabase-sized U-shaped replication timing domains in human cell lines. PLoS Comput Biol. 2012;8(4):1002443. doi:10.1371/journal.pcbi.1002443.

88. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8. doi:10.1038/nmeth.1226.

89. Boulos RE, Julienne H, Baker A, Chen CL, Petryk N, Kahli M, d'Aubenton-Carafa Y, Goldar A, Jensen P, Hyrien O, Thermes C, Arneodo A, Audit B. From the chromatin interaction network to the organization of the human genome into replication N/U-domains. New J Phys. 2014;16:115014.

90. Smit AFA, Hubley R, Green P. RepeatMasker 3.01996–2010. http://www.repeatmasker.org.

Submit your next manuscript to BioMed Central
and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central