Schizosaccharomyces pombe genome-wide nucleosome mapping reveals positioning mechanisms distinct from those of Saccharomyces cerevisiae

Compromised 2-start zigzag chromatin folding in immature mouse retina cells driven by irregularly spaced nucleosomes with short DNA linkers

The formation of condensed heterochromatin is critical for establishing cell-specific transcriptional programs. To reveal structural transitions underlying heterochromatin formation in maturing mouse rod photoreceptors, we apply cryo-EM tomography, AI-assisted deep denoising, and molecular modeling. We find that chromatin isolated from immature retina cells contains many closely apposed nucleosomes with extremely short or absent nucleosome linkers, which are inconsistent with the typical two-start zigzag chromatin folding. In mature retina cells, the fraction of short-linker nucleosomes is much lower, supporting stronger chromatin compaction. By Cryo-EM-assisted nucleosome interaction capture we observe that chromatin in immature retina is enriched with i±1 interactions while chromatin in mature retina contains predominantly i±2 interactions typical of the two-start zigzag. By mesoscale modeling and computational simulation, we clarify that the unusually short linkers typical of immature retina are sufficient to inhibit the two-start zigzag and chromatin compaction by the interference of very short linkers with linker DNA stems. We propose that this short linker composition renders nucleosome arrays more open in immature retina and that, as the linker DNA length increases in mature retina, chromatin fibers become globally condensed via tight zigzag folding. This mechanism may be broadly utilized to introduce higher chromatin folding entropy for epigenomic plasticity.

Phosphorylation of HP1/Swi6 relieves competition with Suv39/Clr4 on nucleosomes and enables H3K9 trimethyl spreading

Heterochromatin formation in Schizosaccharomyces pombe requires the spreading of histone 3 (H3) Lysine 9 (K9) methylation (me) from nucleation centers by the H3K9 methylase, Suv39/Clr4, and the reader protein, HP1/Swi6. To accomplish this, Suv39/Clr4 and HP1/Swi6 have to associate with nucleosomes both nonspecifically, binding DNA and octamer surfaces and specifically, via recognition of methylated H3K9 by their respective chromodomains. However, how both proteins avoid competition for the same nucleosomes in this process is unclear. Here, we show that phosphorylation tunes the nucleosome affinity of HP1/Swi6 such that it preferentially partitions onto Suv39/Clr4's trimethyl product rather than its unmethylated substrates. Preferential partitioning enables efficient conversion from di-to trimethylation on nucleosomes in vitro and H3K9me3 spreading in vivo. Together, our data suggests that phosphorylation of HP1/Swi6 creates a regime that relieves competition with the "read-write" mechanism of Suv39/Clr4 for productive heterochromatin spreading.

Evolution of promoter-proximal pausing enabled a new layer of transcription control

Promoter-proximal pausing of RNA polymerase II (Pol II) is a key regulatory step during transcription. Despite the central role of pausing in gene regulation, we do not understand the evolutionary processes that led to the emergence of Pol II pausing or its transition to a rate-limiting step actively controlled by transcription factors. Here we analyzed transcription in species across the tree of life. Unicellular eukaryotes display a slow acceleration of Pol II near transcription start sites that transitioned to a longer-lived, focused pause in metazoans. This event coincided with the evolution of new subunits in the NELF and 7SK complexes. Depletion of NELF in mammals shifted the promoter-proximal buildup of Pol II from the pause site into the early gene body and compromised transcriptional activation for a set of heat shock genes. Our work details the evolutionary history of Pol II pausing and sheds light on how new transcriptional regulatory mechanisms evolve.

Interpretable deep residual network uncovers nucleosome positioning and associated features

Nucleosomes represent elementary building units of eukaryotic chromosomes and consist of DNA wrapped around a histone octamer flanked by linker DNA segments. Nucleosomes are central in epigenetic pathways and their genomic positioning is associated with regulation of gene expression, DNA replication, DNA methylation and DNA repair, among other functions. Building on prior discoveries that DNA sequences noticeably affect nucleosome positioning, our objective is to identify nucleosome positions and related features across entire genome. Here, we introduce an interpretable framework based on the concepts of deep residual networks (NuPoSe). Trained on high-coverage human experimental MNase-seq data, NuPoSe is able to learn sequence and structural patterns associated with nucleosome organization in human genome. NuPoSe can be also applied to unseen data from different organisms and cell types. Our findings point to 43 informative features, most of them constitute tri-nucleotides, di-nucleotides and one tetra-nucleotide. Most features are significantly associated with the nucleosomal structural characteristics, namely, periodicity of nucleosomal DNA and its location with respect to a histone octamer. Importantly, we show that features derived from the 27 bp linker DNA flanking nucleosomes contribute up to 10% to the quality of the prediction model. This, along with the comprehensive training sets, deep-learning architecture, and feature selection method, may contribute to the NuPoSe's 80-89% classification accuracy on different independent datasets.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

Asymmetrical nucleosomal DNA signatures regulate transcriptional directionality

Despite the symmetrical structure of nucleosomes, in vitro studies have shown that transcription proceeds with different efficiency depending on the orientation of the DNA sequence around them. However, it is unclear whether this functional asymmetry is present in vivo and whether it could regulate transcriptional directionality. Here, we report that the proximal and distal halves of nucleosomal DNA contribute differentially to nucleosome stability in the genome. In +1 nucleosomes, this asymmetry facilitates or hinders transcription depending on the orientation of its underlying DNA, and this difference is associated with an asymmetrical interaction between DNA and histones. These properties are encoded in the DNA signature of +1 nucleosomes, since its incorporation in the two orientations into downstream nucleosomes renders them asymmetrically accessible to MNase and inverts the balance between sense and antisense transcription. Altogether, our results show that nucleosomal DNA endows nucleosomes with asymmetrical properties that modulate the directionality of transcription.

Topoisomerase 1 facilitates nucleosome reassembly at stress genes during recovery

Chromatin remodeling is essential to allow full development of alternative gene expression programs in response to environmental changes. In fission yeast, oxidative stress triggers massive transcriptional changes including the activation of hundreds of genes, with the participation of histone modifying complexes and chromatin remodelers. DNA transcription is associated to alterations in DNA topology, and DNA topoisomerases facilitate elongation along gene bodies. Here, we test whether the DNA topoisomerase Top1 participates in the RNA polymerase II-dependent activation of the cellular response to oxidative stress. Cells lacking Top1 are resistant to H2O2 stress. The transcriptome of Δtop1 strain was not greatly affected in the absence of stress, but activation of the anti-stress gene expression program was more sustained than in wild-type cells. Top1 associated to stress open reading frames. While the nucleosomes of stress genes are partially and transiently evicted during stress, the chromatin configuration remains open for longer times in cells lacking Top1, facilitating RNA polymerase II progression. We propose that, by removing DNA tension arising from transcription, Top1 facilitates nucleosome reassembly and works in synergy with the chromatin remodeler Hrp1 as opposing forces to transcription and to Snf22 / Hrp3 opening remodelers.

HapX-mediated H2B deub1 and SreA-mediated H2A.Z deposition coordinate in fungal iron resistance

Plant pathogens are challenged by host-derived iron starvation or excess during infection, but the mechanism through which pathogens counteract iron stress is unclear. Here, we found that Fusarium graminearum encounters iron excess during the colonization of wheat heads. Deletion of heme activator protein X (FgHapX), siderophore transcription factor A (FgSreA) or both attenuated virulence. Further, we found that FgHapX activates iron storage under iron excess by promoting histone H2B deubiquitination (H2B deub1) at the promoter of the responsible gene. Meanwhile, FgSreA is shown to inhibit genes mediating iron acquisition during iron excess by facilitating the deposition of histone variant H2A.Z and histone 3 lysine 27 trimethylation (H3K27 me3) at the first nucleosome after the transcription start site. In addition, the monothiol glutaredoxin FgGrx4 is responsible for iron sensing and control of the transcriptional activity of FgHapX and FgSreA via modulation of their enrichment at target genes and recruitment of epigenetic regulators, respectively. Taken together, our findings elucidated the molecular mechanisms for adaptation to iron excess mediated by FgHapX and FgSreA during infection in F. graminearum and provide novel insights into regulation of iron homeostasis at the chromatin level in eukaryotes.

The structure of the NuA4-Tip60 complex reveals the mechanism and importance of long-range chromatin modification

Histone acetylation regulates most DNA transactions and is dynamically controlled by highly conserved enzymes. The only essential histone acetyltransferase (HAT) in yeast, Esa1, is part of the 1-MDa NuA4 complex, which plays pivotal roles in both transcription and DNA-damage repair. NuA4 has the unique capacity to acetylate histone targets located several nucleosomes away from its recruitment site. Neither the molecular mechanism of this activity nor its physiological importance are known. Here we report the structure of the Pichia pastoris NuA4 complex, with its core resolved at 3.4-Å resolution. Three subunits, Epl1, Eaf1 and Swc4, intertwine to form a stable platform that coordinates all other modules. The HAT module is firmly anchored into the core while retaining the ability to stretch out over a long distance. We provide structural, biochemical and genetic evidence that an unfolded linker region of the Epl1 subunit is critical for this long-range activity. Specifically, shortening the Epl1 linker causes severe growth defects and reduced H4 acetylation levels over broad chromatin regions in fission yeast. Our work lays the foundations for a mechanistic understanding of NuA4's regulatory role and elucidates how its essential long-range activity is attained.

Epigenetic switching with asymmetric bridging interactions

Gene expression states are often stably sustained in cis despite massively disruptive events like DNA replication. This is achieved by on-going enzymatic activity that maintains parts of the DNA in either heterochromatic (packed) or euchromatic (free) states, each of which is stabilized by both positive feedback and bridging interactions between individual nucleosomes. In contrast to condensed matter, however, the dynamics is not only governed by equilibrium binding interactions but is also mediated by enzymes that recognize and act on specific amino acid tails of the nucleosomes. The mechanical result is that some nucleosomes can bind to one another and form tightly packed polymer configurations, whereas others remain unbound and form free, noncompact polymer configurations. Here, we study the consequences of such an asymmetric interaction pattern on the dynamics of epigenetic switching. We develop a 3D polymer model and show that traits associated with epigenetic switching, such as bistability and epigenetic memory, are permitted by such a model. We find, however, that the experimentally observed burst-like nature of some epigenetic switches is difficult to reproduce by this biologically motivated interaction. Instead, the behavior seen in experiments can be explained by introducing partial confinement, which particularly affects the euchromatic regions of the chromosome.

Opposing Roles of FACT for Euchromatin and Heterochromatin in Yeast

DNA is stored in the nucleus of a cell in a folded state; however, only the necessary genetic information is extracted from the required group of genes. The key to extracting genetic information is chromatin ambivalence. Depending on the chromosomal region, chromatin is characterized into low-density "euchromatin" and high-density "heterochromatin", with various factors being involved in its regulation. Here, we focus on chromatin regulation and gene expression by the yeast FACT complex, which functions in both euchromatin and heterochromatin. FACT is known as a histone H2A/H2B chaperone and was initially reported as an elongation factor associated with RNA polymerase II. In budding yeast, FACT activates promoter chromatin by interacting with the transcriptional activators SBF/MBF via the regulation of G1/S cell cycle genes. In fission yeast, FACT plays an important role in the formation of higher-order chromatin structures and transcriptional repression by binding to Swi6, an HP1 family protein, at heterochromatin. This FACT property, which refers to the alternate chromatin-regulation depending on the binding partner, is an interesting phenomenon. Further analysis of nucleosome regulation within heterochromatin is expected in future studies.

Histone H3 core domain in chromatin with different DNA linker lengths studied by 1H-Detected solid-state NMR spectroscopy

Chromatin, a dynamic protein-DNA complex that regulates eukaryotic genome accessibility and essential functions, is composed of nucleosomes connected by linker DNA with each nucleosome consisting of DNA wrapped around an octamer of histones H2A, H2B, H3 and H4. Magic angle spinning solid-state nuclear magnetic resonance (NMR) spectroscopy can yield unique insights into histone structure and dynamics in condensed nucleosomes and nucleosome arrays representative of chromatin at physiological concentrations. Recently we used J-coupling-based solid-state NMR methods to investigate with residue-specific resolution the conformational dynamics of histone H3 N-terminal tails in 16-mer nucleosome arrays containing 15, 30 or 60 bp DNA linkers. Here, we probe the H3 core domain in the 16-mer arrays as a function of DNA linker length via dipolar coupling-based 1H-detected solid-state NMR techniques. Specifically, we established nearly complete assignments of backbone chemical shifts for H3 core residues in arrays with 15-60 bp DNA linkers reconstituted with 2H,13C,15N-labeled H3. Overall, these chemical shifts were similar irrespective of the DNA linker length indicating no major changes in H3 core conformation. Notably, however, multiple residues at the H3-nucleosomal DNA interface in arrays with 15 bp DNA linkers exhibited relatively pronounced differences in chemical shifts and line broadening compared to arrays with 30 and 60 bp linkers. These findings are consistent with increased heterogeneity in nucleosome packing and structural strain within arrays containing short DNA linkers that likely leads to side-chains of these interfacial residues experiencing alternate conformations or shifts in their rotamer populations relative to arrays with the longer DNA linkers.

Prediction of nucleosome dynamic interval based on long-short-term memory network (LSTM)

Nucleosome localization is a dynamic process and consists of nucleosome dynamic intervals (NDIs). We preprocessed nucleosome sequence data as time series data (TSD) and developed a long short-term memory network (LSTM) model for training time series data (TSD; LSTM-TSD model) using iterative training and feature learning that predicts NDIs with high accuracy. Sn, Sp, Acc, and MCC of the obtained LSTM model is 91.88%, 92.72%, 92.30%, and 84.61%, respectively. LSTM model could precisely predict the NDIs of yeast 16 chromosome. The NDIs contain 90.29% of nucleosome core DNA and 91.20% of nucleosome central sites, indicating that NDIs have high confidence. We found that the binding sites of transcriptional proteins and other proteins are outside NDIs, not in NDIs. These results are important for analysis of nucleosome localization and gene transcriptional regulation.

Local chromatin context regulates the genetic requirements of the heterochromatin spreading reaction

Heterochromatin spreading, the expansion of repressive chromatin structure from sequence-specific nucleation sites, is critical for stable gene silencing. Spreading re-establishes gene-poor constitutive heterochromatin across cell cycles but can also invade gene-rich euchromatin de novo to steer cell fate decisions. How chromatin context (i.e. euchromatic, heterochromatic) or different nucleation pathways influence heterochromatin spreading remains poorly understood. Previously, we developed a single-cell sensor in fission yeast that can separately record heterochromatic gene silencing at nucleation sequences and distal sites. Here we couple our quantitative assay to a genetic screen to identify genes encoding nuclear factors linked to the regulation of heterochromatin nucleation and the distal spreading of gene silencing. We find that mechanisms underlying gene silencing distal to a nucleation site differ by chromatin context. For example, Clr6 histone deacetylase complexes containing the Fkh2 transcription factor are specifically required for heterochromatin spreading at constitutive sites. Fkh2 recruits Clr6 to nucleation-distal chromatin sites in such contexts. In addition, we find that a number of chromatin remodeling complexes antagonize nucleation-distal gene silencing. Our results separate the regulation of heterochromatic gene silencing at nucleation versus distal sites and show that it is controlled by context-dependent mechanisms. The results of our genetic analysis constitute a broad community resource that will support further analysis of the mechanisms underlying the spread of epigenetic silencing along chromatin.

Repressive structures, or heterochromatin, are seeded at specific genome sequences and then “spread” to silence nearby chromosomal regions. While much is known about the factors that seed heterochromatin, the genetic requirements for spreading are less clear. We devised a fission yeast single-cell method to examine how gene silencing is propagated by the heterochromatin spreading process specifically. Here we use this platform to ask if specific genes are required for the spreading process and whether the same or different genes direct spreading from different chromosomal seeding sites. We find a significant number of genes that specifically promote or antagonize the heterochromatin spreading process. However, different genes are required to enact spreading from different seeding sites. These results have potential implications for cell fate specification, where genes are newly silenced by heterochromatin spreading from diverse chromosomal sites. In a central finding, we show that the Clr6 protein complex, which removes chromatin marks linked to active genes, associates with the Forkhead 2 transcription factor to promote spreading of silencing structures from seeding sites at numerous chromosomal loci. In contrast, we show that proteins that remodel chromatin antagonize the spreading of gene silencing.

GC content strongly influences the role of poly(dA) in the intrinsic nucleosome positioning in Saccharomyces cerevisiae

The nucleosome is the basic structural element of genomic DNA packaging and plays a role in transcription, replication, and recombination. Poly(dA) tracts are considered major sequence determinants of nucleosome positioning, although their role is not well understood. Here, we show that the homopolymeric character and the low GC content of poly(dA)s play different roles in nucleosome formation. We found that the inherent low GC content of poly(dA) alone can account for the deep and anisotropic nucleosome depletion at structurally and functionally important regions of promoters and origins of replication. We also show that the level of nucleosome occupancy at poly(dA) is strongly related to the local nucleotide background and its high frequency of occurrence in Saccharomyces cerevisiae does not appear merely to be associated with its intrinsic nucleosome-excluding properties. In addition, we show that the GC content alone can predict more than 60% of the in vitro nucleosome map, providing further evidence that the intrinsic nucleosome positioning is more greatly determined by GC content than poly(dA) stretches. Our results are consistent with a model in which poly(dA) stretches act at two distinct levels: first, by its low GC content, which intrinsically contributes to hinder nucleosome formation, and second, by its contiguous runs of dA that selectively drive the recruitment of non-histone proteins with structural and functional roles.

Long RNA-Mediated Chromatin Regulation in Fission Yeast and Mammals

As part of a complex network of genome control, long regulatory RNAs exert significant influences on chromatin dynamics. Understanding how this occurs could illuminate new avenues for disease treatment and lead to new hypotheses that would advance gene regulatory research. Recent studies using the model fission yeast Schizosaccharomyces pombe (S. pombe) and powerful parallel sequencing technologies have provided many insights in this area. This review will give an overview of key findings in S. pombe that relate long RNAs to multiple levels of chromatin regulation: histone modifications, gene neighborhood regulation in cis and higher-order chromosomal ordering. Moreover, we discuss parallels recently found in mammals to help bridge the knowledge gap between the study systems.

The local integration preference of the Tf1 retrotransposon in Schizosaccharomyces pombe

Transposons are mobile DNAs that can move to different locations in host genomes. The integration site selection of transposons is critical for both themselves and host cells. Studies on the integration of retrotransposons and retroviruses have focused more on the global preference than on the local preference. The local preferences of retrotransposons are usually weak and of large diversity. Here, we analyzed hundreds of thousands of independent integration events of the Tf1 retrotransposon in Schizosaccharomyces pombe. The consensus sequence at the Tf1 integration sites shows a palindromic pattern, which can be divided into four sections, each of them contains one or more CGnTA units with a period of 10 base pairs, indicating interaction with subunits of the integrase oligomer in the pre-integration complex. Moreover, the analysis on the nucleosome occupancy flanking Tf1 target sites shows that Tf1 integration favors regions with one entire nucleosome depletion.

The Active Mechanism of Nucleosome Depletion by Poly(dA:dT) Tracts In Vivo

Poly(dA:dT) tracts cause nucleosome depletion in many species, e.g., at promoters and replication origins. Their intrinsic biophysical sequence properties make them stiff and unfavorable for nucleosome assembly, as probed by in vitro nucleosome reconstitution. The mere correlation between nucleosome depletion over poly(dA:dT) tracts in in vitro reconstituted and in in vivo chromatin inspired an intrinsic nucleosome exclusion mechanism in vivo that is based only on DNA and histone properties. However, we compile here published and new evidence that this correlation does not reflect mechanistic causation. (1) Nucleosome depletion over poly(dA:dT) in vivo is not universal, e.g., very weak in S. pombe. (2) The energy penalty for incorporating poly(dA:dT) tracts into nucleosomes is modest (<10%) relative to ATP hydrolysis energy abundantly invested by chromatin remodelers. (3) Nucleosome depletion over poly(dA:dT) is much stronger in vivo than in vitro if monitored without MNase and (4) actively maintained in vivo. (5) S. cerevisiae promoters evolved a strand-biased poly(dA) versus poly(dT) distribution. (6) Nucleosome depletion over poly(dA) is directional in vivo. (7) The ATP dependent chromatin remodeler RSC preferentially and directionally displaces nucleosomes towards 5′ of poly(dA). Especially distribution strand bias and displacement directionality would not be expected for an intrinsic mechanism. Together, this argues for an in vivo mechanism where active and species-specific read out of intrinsic sequence properties, e.g., by remodelers, shapes nucleosome organization.

Conformational Dynamics of Histone H3 Tails in Chromatin

Chromatin is a supramolecular DNA-protein complex that compacts eukaryotic genomes and regulates their accessibility and functions. Dynamically disordered histone H3 N-terminal tails are among key chromatin regulatory components. Here, we used high-resolution-magic-angle-spinning NMR measurements of backbone amide ¹⁵N spin relaxation rates to investigate, with residue-specific detail, the dynamics and interactions of H3 tails in recombinant ¹³C,¹⁵N-enriched nucleosome arrays containing 15, 30, or 60 bp linker DNA between the nucleosome repeats. These measurements were compared to analogous data available for mononucleosomes devoid of linker DNA or containing two 20 bp DNA overhangs. The H3 tail dynamics in nucleosome arrays were found to be considerably attenuated compared with nucleosomes with or without linker DNA due to transient electrostatic interactions with the linker DNA segments and the structured chromatin environment. Remarkably, however, the H3 tail dynamics were not modulated by the specific linker DNA length within the 15-60 bp range investigated here.

Coupling between Sequence-Mediated Nucleosome Organization and Genome Evolution

The nucleosome is a major modulator of DNA accessibility to other cellular factors. Nucleosome positioning has a critical importance in regulating cell processes such as transcription, replication, recombination or DNA repair. The DNA sequence has an influence on the position of nucleosomes on genomes, although other factors are also implicated, such as ATP-dependent remodelers or competition of the nucleosome with DNA binding proteins. Different sequence motifs can promote or inhibit the nucleosome formation, thus influencing the accessibility to the DNA. Sequence-encoded nucleosome positioning having functional consequences on cell processes can then be selected or counter-selected during evolution. We review the interplay between sequence evolution and nucleosome positioning evolution. We first focus on the different ways to encode nucleosome positions in the DNA sequence, and to which extent these mechanisms are responsible of genome-wide nucleosome positioning in vivo. Then, we discuss the findings about selection of sequences for their nucleosomal properties. Finally, we illustrate how the nucleosome can directly influence sequence evolution through its interactions with DNA damage and repair mechanisms. This review aims to provide an overview of the mutual influence of sequence evolution and nucleosome positioning evolution, possibly leading to complex evolutionary dynamics.

Pickle Recruits Retinoblastoma Related 1 to Control Lateral Root Formation in Arabidopsis

Lateral root (LR) formation is an example of a plant post-embryonic organogenesis event. LRs are issued from non-dividing cells entering consecutive steps of formative divisions, proliferation and elongation. The chromatin remodeling protein PICKLE (PKL) negatively regulates auxin-mediated LR formation through a mechanism that is not yet known. Here we show that PKL interacts with RETINOBLASTOMA-RELATED 1 (RBR1) to repress the LATERAL ORGAN BOUNDARIES-DOMAIN 16 (LBD16) promoter activity. Since LBD16 function is required for the formative division of LR founder cells, repression mediated by the PKL–RBR1 complex negatively regulates formative division and LR formation. Inhibition of LR formation by PKL–RBR1 is counteracted by auxin, indicating that, in addition to auxin-mediated transcriptional responses, the fine-tuned process of LR formation is also controlled at the chromatin level in an auxin-signaling dependent manner.

Establishment of heterochromatin in domain-size-dependent bursts

Methylation of histone H3K9 is a hallmark of epigenetic silencing in eukaryotes. Nucleosome modifications often rely on positive feedback where enzymes are recruited by modified nucleosomes. A combination of local and global feedbacks has been proposed to account for some dynamic properties of heterochromatin, but the range at which the global feedbacks operate and the exact mode of heterochromatin propagation are not known. We investigated these questions in fission yeast. Guided by mathematical modeling, we incrementally increased the size of the mating-type region and profiled heterochromatin establishment over time. We observed exponential decays in the proportion of cells with active reporters, with rates that decreased with domain size. Establishment periods varied from a few generations in wild type to >200 generations in the longest region examined, and highly correlated silencing of two reporters located outside the nucleation center was observed. On a chromatin level, this indicates that individual regions are silenced in sudden bursts. Mathematical modeling accounts for these bursts if heterochromatic nucleosomes facilitate a deacetylation or methylation reaction at long range, in a distance-independent manner. A likely effector of three-dimensional interactions is the evolutionarily conserved Swi6^HP1 H3K9me reader, indicating the bursting behavior might be a general mode of heterochromatin propagation.

Repression of a large number of genes requires interplay between homologous recombination and HIRA

During homologous recombination, Dbl2 protein is required for localisation of Fbh1, an F-box helicase that efficiently dismantles Rad51-DNA filaments. RNA-seq analysis of dbl2Δ transcriptome showed that the dbl2 deletion results in upregulation of more than 500 loci in Schizosaccharomyces pombe. Compared with the loci with no change in expression, the misregulated loci in dbl2Δ are closer to long terminal and long tandem repeats. Furthermore, the misregulated loci overlap with antisense transcripts, retrotransposons, meiotic genes and genes located in subtelomeric regions. A comparison of the expression profiles revealed that Dbl2 represses the same type of genes as the HIRA histone chaperone complex. Although dbl2 deletion does not alleviate centromeric or telomeric silencing, it suppresses the silencing defect at the outer centromere caused by deletion of hip1 and slm9 genes encoding subunits of the HIRA complex. Moreover, our analyses revealed that cells lacking dbl2 show a slight increase of nucleosomes at transcription start sites and increased levels of methylated histone H3 (H3K9me2) at centromeres, subtelomeres, rDNA regions and long terminal repeats. Finally, we show that other proteins involved in homologous recombination, such as Fbh1, Rad51, Mus81 and Rad54, participate in the same gene repression pathway.

Nucleosome Positioning and Spacing: From Mechanism to Function

Eukaryotes associate their genomes with histone proteins, forming nucleosome particles. Nucleosomes regulate and protect the genetic information. They often assemble into evenly spaced arrays of nucleosomes. These regular nucleosome arrays cover significant portions of the genome, in particular over genes. The presence of these evenly spaced nucleosome arrays is highly conserved throughout the entire eukaryotic domain. Here, we review the mechanisms behind the establishment of this primary structure of chromatin with special emphasis on the biogenesis of evenly spaced nucleosome arrays. We highlight the roles that transcription, nucleosome remodelers, DNA sequence, and histone density play towards the formation of evenly spaced nucleosome arrays and summarize our current understanding of their cellular functions. We end with key unanswered questions that remain to be explored to obtain an in-depth understanding of the biogenesis and function of the nucleosome landscape.

Phosphorylation of the Transcription Factor Atf1 at Multiple Sites by the MAP Kinase Sty1 Controls Homologous Recombination and Transcription

Transcription factors are often the downstream effectors of signaling cascades. In fission yeast, the transcription factor Atf1 is phosphorylated by the MAP kinase Sty1 under several environmental stressors to promote transcription initiation of stress genes. However, Sty1 and Atf1 have also been involved in other cellular processes such as homologous recombination at hotspots, ste11 gene expression during mating and meiosis, or regulation of fbp1 gene transcription under glucose starvation conditions. Using different phospho-mutants of Atf1, we have investigated the role of Atf1 phosphorylation by Sty1 in those biological processes. An Atf1 mutant lacking the canonical MAP kinase phosphorylation sites cannot activate fbp1 transcription when glucose is depleted, but it is still able to induce recombination at ade6.M26 and to induce ste11 after nitrogen depletion; in these last cases, Sty1 is still required, suggesting that additional non-canonical sites are activating the transcription factor. In all cases, an Atf1 phosphomimetic mutant bypasses the requirement of the Sty1 kinase in these diverse biological processes, highlighting the essential role of the DNA binding factor Atf1 on chromatin remodeling and cell adaptation to nutritional changes. We propose that post-translational modifications of Atf1 by Sty1, either at canonical or non-canonical sites, are sufficient to activate some of the functions of Atf1, those involving chromatin remodeling and transcription initiation. However, in the case of fbp1 where Atf1 acts synergistically with other transcription factors, elimination of the canonical sites is sufficient to hamper some of the interactions required in this complex scenario and to impair transcription initiation.

Q-Nuc: a bioinformatics pipeline for the quantitative analysis of nucleosomal profiles

Nucleosomal profiling is an effective method to determine the positioning and occupancy of nucleosomes, which is essential to understand their roles in genomic processes. However, the positional randomness across the genome and its relationship with nucleosome occupancy remains poorly understood. Here we present a computational method that segments the profile into nucleosomal domains and quantifies their randomness and relative occupancy level. Applying this method to published data, we find on average ~ 3-fold differences in the degree of positional randomness between regions typically considered "well-ordered", as well as an unexpected predominance of only two types of domains of positional randomness in yeast cells. Further, we find that occupancy levels between domains actually differ maximally by ~ 2-3-fold in both cells, which has not been described before. We also developed a procedure by which one can estimate the sequencing depth that is required to identify nucleosomal positions even when regional positional randomness is high. Overall, we have developed a pipeline to quantitatively characterize domain-level features of nucleosome randomness and occupancy genome-wide, enabling the identification of otherwise unknown features in nucleosomal organization.

YeasTSS: an integrative web database of yeast transcription start sites

The transcription initiation landscape of eukaryotic genes is complex and highly dynamic. In eukaryotes, genes can generate multiple transcript variants that differ in 5' boundaries due to usages of alternative transcription start sites (TSSs), and the abundance of transcript isoforms are highly variable. Due to a large number and complexity of the TSSs, it is not feasible to depict details of transcript initiation landscape of all genes using text-format genome annotation files. Therefore, it is necessary to provide data visualization of TSSs to represent quantitative TSS maps and the core promoters (CPs). In addition, the selection and activity of TSSs are influenced by various factors, such as transcription factors, chromatin remodeling and histone modifications. Thus, integration and visualization of functional genomic data related to these features could provide a better understanding of the gene promoter architecture and regulatory mechanism of transcription initiation. Yeast species play important roles for the research and human society, yet no database provides visualization and integration of functional genomic data in yeast. Here, we generated quantitative TSS maps for 12 important yeast species, inferred their CPs and built a public database, YeasTSS (www.yeastss.org). YeasTSS was designed as a central portal for visualization and integration of the TSS maps, CPs and functional genomic data related to transcription initiation in yeast. YeasTSS is expected to benefit the research community and public education for improving genome annotation, studies of promoter structure, regulated control of transcription initiation and inferring gene regulatory network.

Set1/COMPASS repels heterochromatin invasion at euchromatic sites by disrupting Suv39/Clr4 activity and nucleosome stability

Protection of euchromatin from invasion by gene-repressive heterochromatin is critical for cellular health and viability. In addition to constitutive loci such as pericentromeres and subtelomeres, heterochromatin can be found interspersed in gene-rich euchromatin, where it regulates gene expression pertinent to cell fate. While heterochromatin and euchromatin are globally poised for mutual antagonism, the mechanisms underlying precise spatial encoding of heterochromatin containment within euchromatic sites remain opaque. We investigated ectopic heterochromatin invasion by manipulating the fission yeast mating type locus boundary using a single-cell spreading reporter system. We found that heterochromatin repulsion is locally encoded by Set1/COMPASS on certain actively transcribed genes and that this protective role is most prominent at heterochromatin islands, small domains interspersed in euchromatin that regulate cell fate specifiers. Sensitivity to invasion by heterochromatin, surprisingly, is not dependent on Set1 altering overall gene expression levels. Rather, the gene-protective effect is strictly dependent on Set1's catalytic activity. H3K4 methylation, the Set1 product, antagonizes spreading in two ways: directly inhibiting catalysis by Suv39/Clr4 and locally disrupting nucleosome stability. Taken together, these results describe a mechanism for spatial encoding of euchromatic signals that repel heterochromatin invasion.

HMMRATAC: a Hidden Markov ModeleR for ATAC-seq

ATAC-seq has been widely adopted to identify accessible chromatin regions across the genome. However, current data analysis still utilizes approaches initially designed for ChIP-seq or DNase-seq, without considering the transposase digested DNA fragments that contain additional nucleosome positioning information. We present the first dedicated ATAC-seq analysis tool, a semi-supervised machine learning approach named HMMRATAC. HMMRATAC splits a single ATAC-seq dataset into nucleosome-free and nucleosome-enriched signals, learns the unique chromatin structure around accessible regions, and then predicts accessible regions across the entire genome. We show that HMMRATAC outperforms the popular peak-calling algorithms on published human ATAC-seq datasets. We find that single-end sequenced or size-selected ATAC-seq datasets result in a loss of sensitivity compared to paired-end datasets without size-selection.

Diverse DNA Sequence Motifs Activate Meiotic Recombination Hotspots Through a Common Chromatin Remodeling Pathway

In meiosis, multiple different DNA sequence motifs help to position homologous recombination at hotspots in the genome. How do the seemingly disparate cis-acting regulatory modules each promote locally the activity of the basal recombination machinery? We defined molecular mechanisms of action for five different hotspot-activating DNA motifs (M26, CCAAT, Oligo-C, 4095, 4156) located independently at the same site within the ade6 locus of the fission yeast Schizosaccharomyces pombe Each motif promoted meiotic recombination (i.e., is active) within this context, and this activity required the respective binding proteins (transcription factors Atf1, Pcr1, Php2, Php3, Php5, Rst2). High-resolution analyses of chromatin structure by nucleosome scanning assays revealed that each motif triggers the displacement of nucleosomes surrounding the hotspot motif in meiosis. This chromatin remodeling required the respective sequence-specific binding proteins, was constitutive for two motifs, and was enhanced meiotically for three others. Hotspot activity of each motif strongly required the ATP-dependent chromatin remodeling enzyme Snf22 (Snf2/Swi2), with lesser dependence on Gcn5, Mst2, and Hrp3. These findings support a model in which most meiotic recombination hotspots are positioned by the binding of transcription factors to their respective DNA sites. The functional redundancy of multiple, sequence-specific protein-DNA complexes converges upon shared chromatin remodeling pathways that help provide the basal recombination machinery (Spo11/Rec12 complex) access to its DNA substrates within chromatin.

The histone variant H2A.Z in gene regulation

The histone variant H2A.Z is involved in several processes such as transcriptional control, DNA repair, regulation of centromeric heterochromatin and, not surprisingly, is implicated in diseases such as cancer. Here, we review the recent developments on H2A.Z focusing on its role in transcriptional activation and repression. H2A.Z, as a replication-independent histone, has been studied in several model organisms and inducible mammalian model systems. Its loading machinery and several modifying enzymes have been recently identified, and some of the long-standing discrepancies in transcriptional activation and/or repression are about to be resolved. The buffering functions of H2A.Z, as supported by genome-wide localization and analyzed in several dynamic systems, are an excellent example of transcriptional control. Posttranslational modifications such as acetylation and ubiquitination of H2A.Z, as well as its specific binding partners, are in our view central players in the control of gene expression. Understanding the key-mechanisms in either turnover or stabilization of H2A.Z-containing nucleosomes as well as defining the H2A.Z interactome will pave the way for therapeutic applications in the future.

Nucleosome dynamics of human iPSC during neural differentiation

Nucleosome positioning is important for neurodevelopment, and genes mediating chromatin remodelling are strongly associated with human neurodevelopmental disorders. To investigate changes in nucleosome positioning during neural differentiation, we generate genome‐wide nucleosome maps from an undifferentiated human‐induced pluripotent stem cell (hiPSC) line and after its differentiation to the neural progenitor cell (NPC) stage. We find that nearly 3% of nucleosomes are highly positioned in NPC, but significantly, there are eightfold fewer positioned nucleosomes in pluripotent cells, indicating increased positioning during cell differentiation. Positioned nucleosomes do not strongly correlate with active chromatin marks or gene transcription. Unexpectedly, we find a small population of nucleosomes that occupy similar positions in pluripotent and neural progenitor cells and are found at binding sites of the key gene regulators NRSF/REST and CTCF. Remarkably, the presence of these nucleosomes appears to be independent of the associated regulatory complexes. Together, these results present a scenario in human cells, where positioned nucleosomes are sparse and dynamic, but may act to alter gene expression at a distance via the structural conformation at sites of chromatin regulation.

Nucleosome positioning and spacing: from genome-wide maps to single arrays

The positioning of nucleosomes relative to DNA and their neighboring nucleosomes represents a fundamental layer of chromatin organization. Changes in nucleosome positioning and spacing affect the accessibility of DNA to regulatory factors and the formation of higher order chromatin structures. Sequencing of mononucleosomal fragments allowed mapping nucleosome positions on a genome-wide level in many organisms. This revealed that successions of evenly spaced and well-positioned nucleosomes-so called phased nucleosome arrays-occur at the 5' end of many active genes and in the vicinity of transcription factor and other protein binding sites. Phased arrays arise from the interplay of barrier elements on the DNA, which position adjacent nucleosomes, and the nucleosome spacing activity of ATP-dependent chromatin remodelers. A shortcoming of classic mononucleosomal mapping experiments is that they only reveal nucleosome spacing and array regularity at select sites in the genome with well-positioned nucleosomes. However, new technological approaches elucidate nucleosome array structure throughout the genome and with single-cell resolution. In the future, it will be interesting to see whether changes in nucleosome array regularity and spacing contribute to the formation of higher order chromatin structures and the spatial organization of the genome in vivo.

Comprehensive profiling of the fission yeast transcription start site activity during stress and media response

Fission yeast, Schizosaccharomyces pombe, is an attractive model organism for transcriptional and chromatin biology research. Such research is contingent on accurate annotation of transcription start sites (TSSs). However, comprehensive genome-wide maps of TSSs and their usage across commonly applied laboratory conditions and treatments for S. pombe are lacking. To this end, we profiled TSS activity genome-wide in S. pombe cultures exposed to heat shock, nitrogen starvation, hydrogen peroxide and two commonly applied media, YES and EMM2, using Cap Analysis of Gene Expression (CAGE). CAGE-based annotation of TSSs is substantially more accurate than existing PomBase annotation; on average, CAGE TSSs fall 50-75 bp downstream of PomBase TSSs and co-localize with nucleosome boundaries. In contrast to higher eukaryotes, dispersed TSS distributions are not common in S. pombe. Our data recapitulate known S. pombe stress expression response patterns and identify stress- and media-responsive alternative TSSs. Notably, alteration of growth medium induces changes of similar magnitude as some stressors. We show a link between nucleosome occupancy and genetic variation, and that the proximal promoter region is genetically diverse between S. pombe strains. Our detailed TSS map constitutes a central resource for S. pombe gene regulation research.

Shelterin and subtelomeric DNA sequences control nucleosome maintenance and genome stability

Telomeres and the shelterin complex cap and protect the ends of chromosomes. Telomeres are flanked by the subtelomeric sequences that have also been implicated in telomere regulation, although their role is not well defined. Here, we show that, in Schizosaccharomyces pombe, the telomere-associated sequences (TAS) present on most subtelomeres are hyper-recombinogenic, have metastable nucleosomes, and unusual low levels of H3K9 methylation. Ccq1, a subunit of shelterin, protects TAS from nucleosome loss by recruiting the heterochromatic repressor complexes CLRC and SHREC, thereby linking nucleosome stability to gene silencing. Nucleosome instability at TAS is independent of telomeric repeats and can be transmitted to an intrachromosomal locus containing an ectopic TAS fragment, indicating that this is an intrinsic property of the underlying DNA sequence. When telomerase recruitment is compromised in cells lacking Ccq1, DNA sequences present in the TAS promote recombination between chromosomal ends, independent of nucleosome abundance, implying an active function of these sequences in telomere maintenance. We propose that Ccq1 and fragile subtelomeres co-evolved to regulate telomere plasticity by controlling nucleosome occupancy and genome stability.

Cryo-ET reveals the macromolecular reorganization of S. pombe mitotic chromosomes in vivo

Chromosomes condense during mitosis in most eukaryotes. This transformation involves rearrangements at the nucleosome level and has consequences for transcription. Here, we use cryo-electron tomography (cryo-ET) to determine the 3D arrangement of nuclear macromolecular complexes, including nucleosomes, in frozen-hydrated Schizosaccharomyces pombe cells. Using 3D classification analysis, we did not find evidence that nucleosomes resembling the crystal structure are abundant. This observation and those from other groups support the notion that a subset of fission yeast nucleosomes may be partially unwrapped in vivo. In both interphase and mitotic cells, there is also no evidence of monolithic structures the size of Hi-C domains. The chromatin is mingled with two features: pockets, which are positions free of macromolecular complexes; and "megacomplexes," which are multimegadalton globular complexes like preribosomes. Mitotic chromatin is more crowded than interphase chromatin in subtle ways. Nearest-neighbor distance analyses show that mitotic chromatin is more compacted at the oligonucleosome than the dinucleosome level. Like interphase, mitotic chromosomes contain megacomplexes and pockets. This uneven chromosome condensation helps explain a longstanding enigma of mitosis: a subset of genes is up-regulated.

The INO80 complex activates the transcription of S-phase genes in a cell cycle-regulated manner

Chromatin structure is an essential factor in the proper regulation of DNA repair, DNA replication and transcription. The INO80 complex and the SWR complex have been shown to play a fundamental role in transcription regulation through remodeling chromatin at specific genes and loci. Here, we report that the Schizosaccharomyces pombe INO80 complex physically interacts with the mlui-binding factor (MBF) complex. Furthermore, we are able to detect the INO80 complex in MBF-regulated promoters. Binding of INO80 to these genes is cell cycle regulated, with a maximum binding preceding their transcription and accumulation of their mRNAs. In fact, the INO80 complex is required to fully and timely activate the transcription of these genes. We also show that the accumulation of acetylated H2A.Z at the +1 nucleosome is cell cycle regulated. Cells in which H2A.Z acetylation is abolished still have some cell cycle-regulated transcription of MBF-dependent genes, although to a much lesser extent.

Local enrichment of HP1alpha at telomeres alters their structure and regulation of telomere protection

Enhanced telomere maintenance is evident in malignant cancers. While telomeres are thought to be inherently heterochromatic, detailed mechanisms of how epigenetic modifications impact telomere protection and structures are largely unknown in human cancers. Here we develop a molecular tethering approach to experimentally enrich heterochromatin protein HP1α specifically at telomeres. This results in increased deposition of H3K9me3 at cancer cell telomeres. Telomere extension by telomerase is attenuated, and damage-induced foci at telomeres are reduced, indicating augmentation of telomere stability. Super-resolution STORM imaging shows an unexpected increase in irregularity of telomeric structure. Telomere-tethered chromo shadow domain (CSD) mutant I165A of HP1α abrogates both the inhibition of telomere extension and the irregularity of telomeric structure, suggesting the involvement of at least one HP1α-ligand in mediating these effects. This work presents an approach to specifically manipulate the epigenetic status locally at telomeres to uncover insights into molecular mechanisms underlying telomere structural dynamics.

Long noncoding RNA repertoire and targeting by nuclear exosome, cytoplasmic exonuclease, and RNAi in fission yeast

Long noncoding RNAs (lncRNAs), which are longer than 200 nucleotides but often unstable, contribute a substantial and diverse portion to pervasive noncoding transcriptomes. Most lncRNAs are poorly annotated and understood, although several play important roles in gene regulation and diseases. Here we systematically uncover and analyze lncRNAs in Schizosaccharomyces pombe. Based on RNA-seq data from twelve RNA-processing mutants and nine physiological conditions, we identify 5775 novel lncRNAs, nearly 4× the previously annotated lncRNAs. The expression of most lncRNAs becomes strongly induced under the genetic and physiological perturbations, most notably during late meiosis. Most lncRNAs are cryptic and suppressed by three RNA-processing pathways: the nuclear exosome, cytoplasmic exonuclease, and RNAi. Double-mutant analyses reveal substantial coordination and redundancy among these pathways. We classify lncRNAs by their dominant pathway into cryptic unstable transcripts (CUTs), Xrn1-sensitive unstable transcripts (XUTs), and Dicer-sensitive unstable transcripts (DUTs). XUTs and DUTs are enriched for antisense lncRNAs, while CUTs are often bidirectional and actively translated. The cytoplasmic exonuclease, along with RNAi, dampens the expression of thousands of lncRNAs and mRNAs that become induced during meiosis. Antisense lncRNA expression mostly negatively correlates with sense mRNA expression in the physiological, but not the genetic conditions. Intergenic and bidirectional lncRNAs emerge from nucleosome-depleted regions, upstream of positioned nucleosomes. Our results highlight both similarities and differences to lncRNA regulation in budding yeast. This broad survey of the lncRNA repertoire and characteristics in S. pombe, and the interwoven regulatory pathways that target lncRNAs, provides a rich framework for their further functional analyses.

Application of MNase-Seq in the Global Mapping of Nucleosome Positioning in Plants

The precise positioning of nucleosomes along the underlying DNA is critical for a variety of biological processes, especially in regulating transcription. The interplay between nucleosomes and transcription factors for accessing the underlying DNA sequences is one of the key determinants that affect transcriptional regulation. Moreover, nucleosomes with various packing statuses confer distinct functions in regulating gene expressions in response to various internal or external signals. Therefore, global mapping of nucleosome positions is one informative way to elucidate the relationship between patterns of nucleosome positioning/occupancy and transcriptional regulations. MNase digestion coupled with high-throughput sequencing (MNase-seq) has been utilized widely for global mapping of nucleosome positioning in eukaryotes that have a sequenced genome. We have developed a robust MNase-seq procedure in plants. It mainly includes plant nuclei isolation, treatment of purified nuclei with MNase, gel recovery of MNase-trimmed mononucleosomal DNA with an approximate size of 150 bp, MNase-seq library preparation followed by Illumina sequencing, and data analysis. MNase-seq has already been successfully applied to identify genome-wide nucleosome positioning in model plants, rice, and Arabidopsis thaliana.

Major Determinants of Nucleosome Positioning

The compact structure of the nucleosome limits DNA accessibility and inhibits the binding of most sequence-specific proteins. Nucleosomes are not randomly located on the DNA but positioned with respect to the DNA sequence, suggesting models in which critical binding sites are either exposed in the linker, resulting in activation, or buried inside a nucleosome, resulting in repression. The mechanisms determining nucleosome positioning are therefore of paramount importance for understanding gene regulation and other events that occur in chromatin, such as transcription, replication, and repair. Here, we review our current understanding of the major determinants of nucleosome positioning: DNA sequence, nonhistone DNA-binding proteins, chromatin-remodeling enzymes, and transcription. We outline the major challenges for the future: elucidating the precise mechanisms of chromatin opening and promoter activation, identifying the complexes that occupy promoters, and understanding the multiscale problem of chromatin fiber organization.

Regulation of tRNA gene transcription by the chromatin structure and nucleosome dynamics

The short, non-coding genes transcribed by the RNA polymerase (pol) III, necessary for survival of a cell, need to be repressed under the stress conditions in vivo. The pol III-transcribed genes have adopted several novel chromatin-based regulatory mechanisms to their advantage. In the budding yeast, the sub-nucleosomal size tRNA genes are found in the nucleosome-free regions, flanked by positioned nucleosomes at both the ends. With their chromosomes-wide distribution, all tRNA genes have a different chromatin context. A single nucleosome dynamics controls the accessibility of the genes for transcription. This dynamics operates under the influence of several chromatin modifiers in a gene-specific manner, giving the scope for differential regulation of even the isogenes within a tRNA gene family. The chromatin structure around the pol III-transcribed genes provides a context conducive for steady-state transcription as well as gene-specific transcriptional regulation upon signaling from the environmental cues. This article is part of a Special Issue entitled: SI: Regulation of tRNA synthesis and modification in physiological conditions and disease edited by Dr. Boguta Magdalena.

Tracking the evolution of 3D gene organization demonstrates its connection to phenotypic divergence

It has recently been shown that the organization of genes in eukaryotic genomes, and specifically in 3D, is strongly related to gene expression and function and partially conserved between organisms. However, previous studies of 3D genomic organization analyzed each organism independently from others. Here, we propose an approach for unified inter-organismal analysis of gene organization based on a network representation of Hi-C data. We define and detect four classes of spatially co-evolving orthologous modules (SCOMs), i.e. gene families that co-evolve in their 3D organization, based on patterns of divergence and conservation of distances. We demonstrate our methodology on Hi-C data from Saccharomyces cerevisiae and Schizosaccharomyces pombe, and identify, among others, modules relating to RNA splicing machinery and chromatin silencing by small RNA which are central to S. pombe's lifestyle. Our results emphasize the importance of 3D genomic organization in eukaryotes and suggest that the evolutionary mechanisms that shape gene organization affect the organism fitness and phenotypes. The proposed algorithms can be utilized in future studies of genome evolution and comparative analysis of spatial genomic organization in different tissues, conditions and single cells.

Capturing Structural Heterogeneity in Chromatin Fibers

Chromatin fiber organization is implicated in processes such as transcription, DNA repair and chromosome segregation, but how nucleosomes interact to form higher-order structure remains poorly understood. We solved two crystal structures of tetranucleosomes with approximately 11-bp DNA linker length at 5.8 and 6.7 Å resolution. Minimal intramolecular nucleosome-nucleosome interactions result in a fiber model resembling a flat ribbon that is compatible with a two-start helical architecture, and that exposes histone and DNA surfaces to the environment. The differences in the two structures combined with electron microscopy reveal heterogeneous structural states, and we used site-specific chemical crosslinking to assess the diversity of nucleosome-nucleosome interactions through identification of structure-sensitive crosslink sites that provide a means to characterize fibers in solution. The chromatin fiber architectures observed here provide a basis for understanding heterogeneous chromatin higher-order structures as they occur in a genomic context.

Nucleosomal organization and DNA base composition patterns

Nucleosomes are the basic units of chromatin. They compact the genome inside the nucleus and regulate the access of proteins to DNA. In the yeast genome, most nucleosomes occupy well-defined positions, which are maintained under many different physiological situations and genetic backgrounds. Although several short sequence elements have been described that favor or reduce the affinity between histones and DNA, the extent to which the DNA sequence affects nucleosome positioning in the genomic context remains unclear. Recent analyses indicate that the base composition pattern of mononucleosomal DNA differs among species, and that the same sequence elements have a different impact on nucleosome positioning in different genomes despite the high level of phylogenetic conservation of histones. These studies have also shown that the DNA sequence contributes to nucleosome positioning to the point that it is possible to design synthetic DNA molecules capable of generating regular and species-specific nucleosomal patterns in vivo.

A comparison of nucleosome organization in Drosophila cell lines

Changes in the distribution of nucleosomes along the genome influence chromatin structure and impact gene expression by modulating the accessibility of DNA to transcriptional machinery. However, the role of genome-wide nucleosome positioning in gene expression and in maintaining differentiated cell states remains poorly understood. Drosophila melanogaster cell lines represent distinct tissue types and exhibit cell-type specific gene expression profiles. They thus could provide a useful tool for investigating cell-type specific nucleosome organization of an organism's genome. To evaluate this possibility, we compared genome-wide nucleosome positioning and occupancy in five different Drosophila tissue-specific cell lines, and in reconstituted chromatin, and then tested for correlations between nucleosome positioning, transcription factor binding motifs, and gene expression. Nucleosomes in all cell lines were positioned in accordance with previously known DNA-nucleosome interactions, with helically repeating A/T di-nucleotide pairs arranged within nucleosomal DNAs and AT-rich pentamers generally excluded from nucleosomal DNA. Nucleosome organization in all cell lines differed markedly from in vitro reconstituted chromatin, with highly expressed genes showing strong nucleosome organization around transcriptional start sites. Importantly, comparative analysis identified genomic regions that exhibited cell line-specific nucleosome enrichment or depletion. Further analysis of these regions identified 91 out of 16,384 possible heptamer sequences that showed differential nucleosomal occupation between cell lines, and 49 of the heptamers matched one or more known transcription factor binding sites. These results demonstrate that there is differential nucleosome positioning between these Drosophila cell lines and therefore identify a system that could be used to investigate the functional significance of differential nucleosomal positioning in cell type specification.

Spt5 Plays Vital Roles in the Control of Sense and Antisense Transcription Elongation

Spt5 is an essential and conserved factor that functions in transcription and co-transcriptional processes. However, many aspects of the requirement for Spt5 in transcription are poorly understood. We have analyzed the consequences of Spt5 depletion in Schizosaccharomyces pombe using four genome-wide approaches. Our results demonstrate that Spt5 is crucial for a normal rate of RNA synthesis and distribution of RNAPII over transcription units. In the absence of Spt5, RNAPII localization changes dramatically, with reduced levels and a relative accumulation over the first ∼500 bp, suggesting that Spt5 is required for transcription past a barrier. Spt5 depletion also results in widespread antisense transcription initiating within this barrier region. Deletions of this region alter the distribution of RNAPII on the sense strand, suggesting that the barrier observed after Spt5 depletion is normally a site at which Spt5 stimulates elongation. Our results reveal a global requirement for Spt5 in transcription elongation.

Genomes of Multicellular Organisms Have Evolved to Attract Nucleosomes to Promoter Regions

Sequences that influence nucleosome positioning in promoter regions, and their relation to gene regulation, have been the topic of much research over the last decade. In yeast, significant nucleosome-depleted regions are found, which facilitate transcription. With the arrival of nucleosome positioning maps for the human genome, it was discovered that in our genome, unlike in that of yeast, promoters encode for high nucleosome occupancy. In this work, we look at the genomes of a range of different organisms, to provide a catalog of nucleosome positioning signals in promoters across the tree of life. We utilize a computational model of the nucleosome, based on crystallographic analyses of the structure and elasticity of the nucleosome, to predict the nucleosome positioning signals in promoter regions. To be able to apply our model to large genomic datasets, we introduce an approximative scheme that makes use of the limited range of correlations in nucleosomal sequence preferences to create a computationally efficient approximation of the full biophysical model. Our predictions show that a clear distinction between unicellular and multicellular life is visible in the intrinsically encoded nucleosome affinity. Furthermore, the strength of the nucleosome positioning signals correlates with the complexity of the organism. We conclude that encoding for high nucleosome occupancy, as in the human genome, is in fact a universal feature of multicellular life.

Dissecting relative contributions of cis- and trans-determinants to nucleosome distribution by comparing Tetrahymena macronuclear and micronuclear chromatin

The ciliate protozoan Tetrahymena thermophila contains two types of structurally and functionally differentiated nuclei: the transcriptionally active somatic macronucleus (MAC) and the transcriptionally silent germ-line micronucleus (MIC). Here, we demonstrate that MAC features well-positioned nucleosomes downstream of transcription start sites and flanking splice sites. Transcription-associated trans-determinants promote nucleosome positioning in MAC. By contrast, nucleosomes in MIC are dramatically delocalized. Nucleosome occupancy in MAC and MIC are nonetheless highly correlated with each other, as well as with in vitro reconstitution and predictions based upon DNA sequence features, revealing unexpectedly strong contributions from cis-determinants. In particular, well-positioned nucleosomes are often matched with GC content oscillations. As many nucleosomes are coordinately accommodated by both cis- and trans-determinants, we propose that their distribution is shaped by the impact of these nucleosomes on the mutational and transcriptional landscape, and driven by evolutionary selection.

Spatial organization of the Schizosaccharomyces pombe genome within the nucleus

The fission yeast Schizosaccharomyces pombe is a useful experimental system for studying the organization of chromosomes within the cell nucleus. S. pombe has a small genome that is organized into three chromosomes. The small size of the genome and the small number of chromosomes are advantageous for cytological and genome-wide studies of chromosomes; however, the small size of the nucleus impedes microscopic observations owing to limits in spatial resolution during imaging. Recent advances in microscopy, such as super-resolution microscopy, have greatly expanded the use of S. pombe as a model organism in a wide range of studies. In addition, biochemical studies, such as chromatin immunoprecipitation and chromosome conformation capture, have provided complementary approaches. Here, we review the spatial organization of the S. pombe genome as determined by a combination of cytological and biochemical studies. Copyright © 2016 John Wiley & Sons, Ltd.