Liqun Xi

A map of nucleosome positions in yeast at base-pair resolution.

The exact positions of nucleosomes along genomic DNA can influence many aspects of chromosome function, yet existing methods for mapping nucleosomes do not provide the necessary single base pair accuracy to determine these positions. Here we develop and apply a new approach for direct mapping of nucleosome centers based on chemical modification of engineered histones. The resulting map locates nucleosome positions genome-wide in unprecedented detail and accuracy. It reveals novel aspects of the in vivo nucleosome organization that are linked to transcription factor binding, RNA polymerase pausing, and the higher order structure of the chromatin fiber.The exact positions of nucleosomes along genomic DNA can influence many aspects of chromosome function. However, existing methods for mapping nucleosomes do not provide the necessary single-base-pair accuracy to determine these positions. Here we develop and apply a new approach for direct mapping of nucleosome centres on the basis of chemical modification of engineered histones. The resulting map locates nucleosome positions genome-wide in unprecedented detail and accuracy. It shows new aspects of the in vivo nucleosome organization that are linked to transcription factor binding, RNA polymerase pausing and the higher-order structure of the chromatin fibre.

Predicting nucleosome positioning using a duration Hidden Markov Model

BackgroundThe nucleosome is the fundamental packing unit of DNAs in eukaryotic cells. Its detailed positioning on the genome is closely related to chromosome functions. Increasing evidence has shown that genomic DNA sequence itself is highly predictive of nucleosome positioning genome-wide. Therefore a fast software tool for predicting nucleosome positioning can help understanding how a genomes nucleosome organization may facilitate genome function.ResultsWe present a duration Hidden Markov model for nucleosome positioning prediction by explicitly modeling the linker DNA length. The nucleosome and linker models trained from yeast data are re-scaled when making predictions for other species to adjust for differences in base composition. A software tool named NuPoP is developed in three formats for free download.ConclusionsSimulation studies show that modeling the linker length distribution and utilizing a base composition re-scaling method both improve the prediction of nucleosome positioning regarding sensitivity and false discovery rate. NuPoP provides a user-friendly software tool for predicting the nucleosome occupancy and the most probable nucleosome positioning map for genomic sequences of any size. When compared with two existing methods, NuPoP shows improved performance in sensitivity.

Preferentially quantized linker DNA lengths in Saccharomyces cerevisiae

The exact lengths of linker DNAs connecting adjacent nucleosomes specify the intrinsic three-dimensional structures of eukaryotic chromatin fibers. Some studies suggest that linker DNA lengths preferentially occur at certain quantized values, differing one from another by integral multiples of the DNA helical repeat, ∼10 bp; however, studies in the literature are inconsistent. Here, we investigate linker DNA length distributions in the yeast Saccharomyces cerevisiae genome, using two novel methods: a Fourier analysis of genomic dinucleotide periodicities adjacent to experimentally mapped nucleosomes and a duration hidden Markov model applied to experimentally defined dinucleosomes. Both methods reveal that linker DNA lengths in yeast are preferentially periodic at the DNA helical repeat (∼10 bp), obeying the forms 10n+5 bp (integer n). This 10 bp periodicity implies an ordered superhelical intrinsic structure for the average chromatin fiber in yeast.

Single-cell nucleosome mapping reveals the molecular basis of gene expression heterogeneity

Significance Nucleosomes limit access to DNA, which antagonizes gene expression and prevents recruitment of transcription factors that cannot bind DNA wrapped around the histone octamer. Numerous studies using large cell populations determined that active genes promoters tend to be nucleosome-depleted. We developed a method to examine nucleosome positioning in single cells and revealed significant heterogeneity of nucleosome configurations within a population. In an inactive gene loaded with nucleosomes, a small subpopulation of nucleosome-depleted cells exists that were engaged in transcription. Single-cell mapping revealed that even in apparently nucleosome-free regions, some cells were occupied by nucleosomes. These data reveal an underlying complexity of nucleosome positioning and its role in regulating gene expression. Nucleosomes, the basic unit of chromatin, have a critical role in the control of gene expression. Nucleosome positions have generally been determined by examining bulk populations of cells and then correlated with overall gene expression. Here, we describe a technique to determine nucleosome positioning in single cells by virtue of the ability of the nucleosome to protect DNA from GpC methylation. In the acid phosphatase inducible PHO5 gene, we find that there is significant cell-to-cell variation in nucleosome positions and shifts in nucleosome positioning correlate with changes in gene expression. However, nucleosome positioning is not absolute, and even with major shifts in gene expression, some cells fail to change nucleosome configuration. Mutations of the PHO5 promoter that introduce a poly(dA:dT) tract-stimulated gene expression under nonpermissive conditions led to shifts of positioned nucleosomes similar to induction of PHO5. By contrast, mutations that altered AA/TT/AT periodicity reduced gene expression upon PHO5 induction and stabilized nucleosomes in most cells, suggesting that enhanced nucleosome affinity for DNA antagonizes chromatin remodelers. Finally, we determined nucleosome positioning in two regions described as “fuzzy” or nucleosome-free when examined in a bulk assay. These regions consisted of distinct nucleosomes with a larger footprint for potential location and an increase population of cells lacking a nucleosome altogether. These data indicate an underlying complexity of nucleosome positioning that may contribute to the flexibility and heterogeneity of gene expression.

Chemical map of Schizosaccharomyces pombe reveals species-specific features in nucleosome positioning

Significance This paper presents a high-resolution map of nucleosome positions of Schizosaccharomyces pombe. Comparison with the high-resolution map of Saccharomyces cerevisiae has provided important insights into nucleosome–DNA interaction and mechanistic variation in nucleosome positioning. The map shows a preponderance of linker lengths centered on 4/5 bp, placing adjacent nucleosomes on opposite faces of the DNA. The dinucleotide signature for nucleosome positioning is equally strong in exons as in introns. Unexpectedly, S. pombe nucleosomes have a preference for A/T residues surrounding the nucleosome dyad, and nucleosome occupancy is very mildly affected by poly (dA-dT) tracts. The preference for A/T residues around the dyad and its role in nucleosome phasing suggest a coevolution of genomes with the DNA binding preferences of nucleosomes across species. Using a recently developed chemical approach, we have generated a genome-wide map of nucleosomes in vivo in Schizosaccharomyces pombe (S. pombe) at base pair resolution. The shorter linker length previously identified in S. pombe is due to a preponderance of nucleosomes separated by ∼4/5 bp, placing nucleosomes on opposite faces of the DNA. The periodic dinucleotide feature thought to position nucleosomes is equally strong in exons as in introns, demonstrating that nucleosome positioning information can be superimposed on coding information. Unlike the case in Saccharomyces cerevisiae, A/T-rich sequences are enriched in S. pombe nucleosomes, particularly at ±20 bp around the dyad. This difference in nucleosome binding preference gives rise to a major distinction downstream of the transcription start site, where nucleosome phasing is highly predictable by A/T frequency in S. pombe but not in S. cerevisiae, suggesting that the genomes and DNA binding preferences of nucleosomes have coevolved in different species. The poly (dA-dT) tracts affect but do not deplete nucleosomes in S. pombe, and they prefer special rotational positions within the nucleosome, with longer tracts enriched in the 10- to 30-bp region from the dyad. S. pombe does not have a well-defined nucleosome-depleted region immediately upstream of most transcription start sites; instead, the −1 nucleosome is positioned with the expected spacing relative to the +1 nucleosome, and its occupancy is negatively correlated with gene expression. Although there is generally very good agreement between nucleosome maps generated by chemical cleavage and micrococcal nuclease digestion, the chemical map shows consistently higher nucleosome occupancy on DNA with high A/T content.

Archaeal nucleosome positioning in vivo and in vitro is directed by primary sequence motifs

BackgroundHistone wrapping of DNA into nucleosomes almost certainly evolved in the Archaea, and predates Eukaryotes. In Eukaryotes, nucleosome positioning plays a central role in regulating gene expression and is directed by primary sequence motifs that together form a nucleosome positioning code. The experiments reported were undertaken to determine if archaeal histone assembly conforms to the nucleosome positioning code.ResultsEukaryotic nucleosome positioning is favored and directed by phased helical repeats of AA/TT/AT/TA and CC/GG/CG/GC dinucleotides, and disfavored by longer AT-rich oligonucleotides. Deep sequencing of genomic DNA protected from micrococcal nuclease digestion by assembly into archaeal nucleosomes has established that archaeal nucleosome assembly is also directed and positioned by these sequence motifs, both in vivo in Methanothermobacter thermautotrophicus and Thermococcus kodakarensis and in vitro in reaction mixtures containing only one purified archaeal histone and genomic DNA. Archaeal nucleosomes assembled at the same locations in vivo and in vitro, with much reduced assembly immediately upstream of open reading frames and throughout the ribosomal rDNA operons. Providing further support for a common positioning code, archaeal histones assembled into nucleosomes on eukaryotic DNA and eukaryotic histones into nucleosomes on archaeal DNA at the same locations. T. kodakarensis has two histones, designated HTkA and HTkB, and strains with either but not both histones deleted grow normally but do exhibit transcriptome differences. Comparisons of the archaeal nucleosome profiles in the intergenic regions immediately upstream of genes that exhibited increased or decreased transcription in the absence of HTkA or HTkB revealed substantial differences but no consistent pattern of changes that would correlate directly with archaeal nucleosome positioning inhibiting or stimulating transcription.ConclusionsThe results obtained establish that an archaeal histone and a genome sequence together are sufficient to determine where archaeal nucleosomes preferentially assemble and where they avoid assembly. We confirm that the same nucleosome positioning code operates in Archaea as in Eukaryotes and presumably therefore evolved with the histone-fold mechanism of DNA binding and compaction early in the archaeal lineage, before the divergence of Eukaryotes.

A Chemical Approach to Mapping Nucleosomes at Base Pair Resolution in Yeast

Most eukaryotic DNA exists in DNA-protein complexes known as nucleosomes. The exact locations of nucleosomes along the genome play a critical role in chromosome functions and gene regulation. However, the current methods for nucleosome mapping do not provide the necessary accuracy to identify the precise nucleosome locations. Here we describe a new experimental approach that directly maps nucleosome center locations in vivo genome-wide at single base pair resolution.

High-resolution nucleosome mapping of targeted regions using BAC-based enrichment

We report a target enrichment method to map nucleosomes of large genomes at unprecedented coverage and resolution by deeply sequencing locus-specific mononucleosomal DNA enriched via hybridization with bacterial artificial chromosomes. We achieved ∼10 000-fold enrichment of specific loci, which enabled sequencing nucleosomes at up to ∼500-fold higher coverage than has been reported in a mammalian genome. We demonstrate the advantages of generating high-sequencing coverage for mapping the center of discrete nucleosomes, and we show the use of the method by mapping nucleosomes during T cell differentiation using nuclei from effector T-cells differentiated from clonal, isogenic, naïve, primary murine CD4 and CD8 T lymphocytes. The analysis reveals that discrete nucleosomes exhibit cell type-specific occupancy and positioning depending on differentiation status and transcription. This method is widely applicable to mapping many features of chromatin and discerning its landscape in large genomes at unprecedented resolution.

Estimating the population size with a behavioral response in capture-recapture experiment

A new estimating procedure is suggested to estimate the population size in a capture-recapture experiment. The capture intensities for first-capture and recapture are allowed to be different and time dependent but they are assumed to be proportional. It is shown that the information on the proportionality constant is crucial to the estimation of the population size. Sensitivity analysis with a misspecification of the proportionality constant is conducted. The method has also been extended to the case with an unknown proportionality. A real example is given.

Differential Nucleosome Occupancies across Oct4-Sox2 Binding Sites in Murine Embryonic Stem Cells.

The binding sequence for any transcription factor can be found millions of times within a genome, yet only a small fraction of these sequences encode functional transcription factor binding sites. One of the reasons for this dichotomy is that many other factors, such as nucleosomes, compete for binding. To study how the competition between nucleosomes and transcription factors helps determine a functional transcription factor site from a predicted transcription factor site, we compared experimentally-generated in vitro nucleosome occupancy with in vivo nucleosome occupancy and transcription factor binding in murine embryonic stem cells. Using a solution hybridization enrichment technique, we generated a high-resolution nucleosome map from targeted regions of the genome containing predicted sites and functional sites of Oct4/Sox2 regulation. We found that at Pax6 and Nes, which are bivalently poised in stem cells, functional Oct4 and Sox2 sites show high amounts of in vivo nucleosome displacement compared to in vitro. Oct4 and Sox2, which are active, show no significant displacement of in vivo nucleosomes at functional sites, similar to nonfunctional Oct4/Sox2 binding. This study highlights a complex interplay between Oct4 and Sox2 transcription factors and nucleosomes among different target genes, which may result in distinct patterns of stem cell gene regulation.