Genome-wide profiling reveals functional interplay of DNA sequence composition, transcriptional activity, and nucleosome positioning in driving DNA supercoiling and helix destabilization in C. elegans.

DNA topology and alternative DNA structures are implicated in regulating diverse biological processes. Although biomechanical properties of these structures have been studied extensively in vitro, characterization in vivo, particularly in multicellular organisms, is limited. We devised new methods to map DNA supercoiling and single-stranded DNA in Caenorhabditis elegans embryos and diapause larvae. To map supercoiling, we quantified the incorporation of biotinylated psoralen into DNA using high-throughput sequencing. To map single-stranded DNA, we combined permanganate treatment with genome-wide sequencing of induced double-stranded breaks. We found high levels of negative supercoiling at transcription start sites (TSSs) in embryos. GC-rich regions flanked by a sharp GC-to-AT transition delineate boundaries of supercoil propagation. In contrast to TSSs in embryos, TSSs in diapause larvae showed dramatic reductions in negative supercoiling without concomitant attenuation of transcription, suggesting developmental-stage-specific regulation. To assess whether alternative DNA structures control chromosome architecture and gene expression, we examined DNA supercoiling in the context of X-Chromosome dosage compensation. We showed that the condensin dosage compensation complex creates negative supercoils locally at its highest-occupancy binding sites but found no evidence for large-scale supercoiling domains along X Chromosomes. In contrast to transcription-coupled negative supercoiling, single-strandedness, which is most pronounced at transcript end sites, is dependent on high AT content and symmetrically positioned nucleosomes. We propose that sharp transitions in sequence composition at functional genomic elements constitute a common regulatory code and that DNA structure and propagation of torsional stress at regulatory elements are critical parameters in shaping important developmental events.

Tripartite organization of centromeric chromatin in budding yeast.

The centromere is the genetic locus that organizes the proteinaceous kinetochore and is responsible for attachment of the chromosome to the spindle at mitosis and meiosis. In most eukaryotes, the centromere consists of highly repetitive DNA sequences that are occupied by nucleosomes containing the CenH3 histone variant, whereas in budding yeast, a ∼120-bp centromere DNA element (CDE) that is sufficient for centromere function is occupied by a single right-handed histone variant CenH3 (Cse4) nucleosome. However, these in vivo observations are inconsistent with in vitro evidence for left-handed octameric CenH3 nucleosomes. To help resolve these inconsistencies, we characterized yeast centromeric chromatin at single base-pair resolution. Intact particles containing both Cse4 and H2A are precisely protected from micrococcal nuclease over the entire CDE of all 16 yeast centromeres in both solubilized chromatin and the insoluble kinetochore. Small DNA-binding proteins protect CDEI and CDEIII and delimit the centromeric nucleosome to the ∼80-bp CDEII, only enough for a single DNA wrap. As expected for a tripartite organization of centromeric chromatin, loss of Cbf1 protein, which binds to CDEI, both reduces the size of the centromere-protected region and shifts its location toward CDEIII. Surprisingly, Cse4 overproduction caused genome-wide misincorporation of nonfunctional CenH3-containing nucleosomes that protect ∼135 base pairs and are preferentially enriched at sites of high nucleosome turnover. Our detection of two forms of CenH3 nucleosomes in the yeast genome, a singly wrapped particle at the functional centromere and octamer-sized particles on chromosome arms, reconcile seemingly conflicting in vivo and in vitro observations.

Distinct chromatin features characterize different classes of repeat sequences in Drosophila melanogaster.

BACKGROUND: Repeat sequences are abundant in eukaryotic genomes but many are excluded from genome assemblies. In Drosophila melanogaster classical studies of repeat content suggested variability between individuals, but they lacked the precision of modern high throughput sequencing technologies. Genome-wide profiling of chromatin features such as histone tail modifications and DNA-binding proteins relies on alignment to the reference genome and hence excludes highly repetitive sequences. RESULTS: By analyzing repeat libraries, sequence complexity and k-mer counts we determined the abundances of different D. melanogaster repeat classes in flies in two public datasets, DGRP and modENCODE. We found that larval DNA was depleted of all repeat classes relative to adult and embryonic DNA, as expected from the known depletion of repeat-rich pericentromeric regions during polytenization of larval tissues. By applying a method that is independent of alignment to the genome assembly, we found that satellite repeats associate with distinct H3 tail modifications, such as H3K9me2 and H3K9me3 for short repeats and H3K9me1 for 359 bp repeats. Short AT-rich repeats however are depleted of nucleosomes and hence all histone modifications and associated chromatin proteins. CONCLUSIONS: The total repeat content and association of repeat sequences with chromatin modifications can be determined despite repeats being excluded from genome assemblies, revealing unexpected distinctions in chromatin features based on sequence composition.

Epigenome characterization at single base-pair resolution.

We have combined standard micrococcal nuclease (MNase) digestion of nuclei with a modified protocol for constructing paired-end DNA sequencing libraries to map both nucleosomes and subnucleosome-sized particles at single base-pair resolution throughout the budding yeast genome. We found that partially unwrapped nucleosomes and subnucleosome-sized particles can occupy the same position within a cell population, suggesting dynamic behavior. By varying the time of MNase digestion, we have been able to observe changes that reflect differential sensitivity of particles, including the eviction of nucleosomes. To characterize DNA-binding features of transcription factors, we plotted the length of each fragment versus its position in the genome, which defined the minimal protected region of each factor. This process led to the precise mapping of protected and exposed regions at and around binding sites, and also determination of the degree to which they are flanked by phased nucleosomes and subnucleosome-sized particles. Our protocol and mapping method provide a general strategy for epigenome characterization, including nucleosome phasing and dynamics, ATP-dependent nucleosome remodelers, and transcription factors, from a single-sequenced sample.

Functional screening in human HSPCs identifies optimized protein-based enhancers of Homology Directed Repair.

Homology Directed Repair (HDR) enables precise genome editing, but the implementation of HDR-based therapies is hindered by limited efficiency in comparison to methods that exploit alternative DNA repair routes, such as Non-Homologous End Joining (NHEJ). In this study, we develop a functional, pooled screening platform to identify protein-based reagents that improve HDR in human hematopoietic stem and progenitor cells (HSPCs). We leverage this screening platform to explore sequence diversity at the binding interface of the NHEJ inhibitor i53 and its target, 53BP1, identifying optimized variants that enable new intermolecular bonds and robustly increase HDR. We show that these variants specifically reduce insertion-deletion outcomes without increasing off-target editing, synergize with a DNAPK inhibitor molecule, and can be applied at manufacturing scale to increase the fraction of cells bearing repaired alleles. This screening platform can enable the discovery of future gene editing reagents that improve HDR outcomes.

The budding yeast Centromere DNA Element II wraps a stable Cse4 hemisome in either orientation in vivo.

In budding yeast, a single cenH3 (Cse4) nucleosome occupies the ∼120-bp functional centromere, however conflicting structural models for the particle have been proposed. To resolve this controversy, we have applied H4S47C-anchored cleavage mapping, which reveals the precise position of histone H4 in every nucleosome in the genome. We find that cleavage patterns at centromeres are unique within the genome and are incompatible with symmetrical structures, including octameric nucleosomes and (Cse4/H4)2 tetrasomes. Centromere cleavage patterns are compatible with a precisely positioned core structure, one in which each of the 16 yeast centromeres is occupied by oppositely oriented Cse4/H4/H2A/H2B hemisomes in two rotational phases within the population. Centromere-specific hemisomes are also inferred from distances observed between closely-spaced H4 cleavages, as predicted from structural modeling. Our results indicate that the orientation and rotational position of the stable hemisome at each yeast centromere is not specified by the functional centromere sequence. DOI: http://dx.doi.org/10.7554/eLife.01861.001.

Modeling disordered regions in proteins using Rosetta.

Protein structure prediction methods such as Rosetta search for the lowest energy conformation of the polypeptide chain. However, the experimentally observed native state is at a minimum of the free energy, rather than the energy. The neglect of the missing configurational entropy contribution to the free energy can be partially justified by the assumption that the entropies of alternative folded states, while very much less than unfolded states, are not too different from one another, and hence can be to a first approximation neglected when searching for the lowest free energy state. The shortcomings of current structure prediction methods may be due in part to the breakdown of this assumption. Particularly problematic are proteins with significant disordered regions which do not populate single low energy conformations even in the native state. We describe two approaches within the Rosetta structure modeling methodology for treating such regions. The first does not require advance knowledge of the regions likely to be disordered; instead these are identified by minimizing a simple free energy function used previously to model protein folding landscapes and transition states. In this model, residues can be either completely ordered or completely disordered; they are considered disordered if the gain in entropy outweighs the loss of favorable energetic interactions with the rest of the protein chain. The second approach requires identification in advance of the disordered regions either from sequence alone using for example the DISOPRED server or from experimental data such as NMR chemical shifts. During Rosetta structure prediction calculations the disordered regions make only unfavorable repulsive contributions to the total energy. We find that the second approach has greater practical utility and illustrate this with examples from de novo structure prediction, NMR structure calculation, and comparative modeling.