Evolutionary divergence of intrinsic and trans-regulated nucleosome positioning sequences reveals plastic rules for chromatin organization

Species-aware DNA language models capture regulatory elements and their evolution

BACKGROUND

The rise of large-scale multi-species genome sequencing projects promises to shed new light on how genomes encode gene regulatory instructions. To this end, new algorithms are needed that can leverage conservation to capture regulatory elements while accounting for their evolution.

RESULTS

Here, we introduce species-aware DNA language models, which we trained on more than 800 species spanning over 500 million years of evolution. Investigating their ability to predict masked nucleotides from context, we show that DNA language models distinguish transcription factor and RNA-binding protein motifs from background non-coding sequence. Owing to their flexibility, DNA language models capture conserved regulatory elements over much further evolutionary distances than sequence alignment would allow. Remarkably, DNA language models reconstruct motif instances bound in vivo better than unbound ones and account for the evolution of motif sequences and their positional constraints, showing that these models capture functional high-order sequence and evolutionary context. We further show that species-aware training yields improved sequence representations for endogenous and MPRA-based gene expression prediction, as well as motif discovery.

CONCLUSIONS

Collectively, these results demonstrate that species-aware DNA language models are a powerful, flexible, and scalable tool to integrate information from large compendia of highly diverged genomes.

Establishment and function of chromatin organization at replication origins

The origin recognition complex (ORC) is essential for initiation of eukaryotic chromosome replication as it loads the replicative helicase-the minichromosome maintenance (MCM) complex-at replication origins1. Replication origins display a stereotypic nucleosome organization with nucleosome depletion at ORC-binding sites and flanking arrays of regularly spaced nucleosomes2-4. However, how this nucleosome organization is established and whether this organization is required for replication remain unknown. Here, using genome-scale biochemical reconstitution with approximately 300 replication origins, we screened 17 purified chromatin factors from budding yeast and found that the ORC established nucleosome depletion over replication origins and flanking nucleosome arrays by orchestrating the chromatin remodellers INO80, ISW1a, ISW2 and Chd1. The functional importance of the nucleosome-organizing activity of the ORC was demonstrated by orc1 mutations that maintained classical MCM-loader activity but abrogated the array-generation activity of ORC. These mutations impaired replication through chromatin in vitro and were lethal in vivo. Our results establish that ORC, in addition to its canonical role as the MCM loader, has a second crucial function as a master regulator of nucleosome organization at the replication origin, a crucial prerequisite for efficient chromosome replication.

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.

Ruler elements in chromatin remodelers set nucleosome array spacing and phasing

Arrays of regularly spaced nucleosomes dominate chromatin and are often phased by alignment to reference sites like active promoters. How the distances between nucleosomes (spacing), and between phasing sites and nucleosomes are determined remains unclear, and specifically, how ATP-dependent chromatin remodelers impact these features. Here, we used genome-wide reconstitution to probe how Saccharomyces cerevisiae ATP-dependent remodelers generate phased arrays of regularly spaced nucleosomes. We find that remodelers bear a functional element named the 'ruler' that determines spacing and phasing in a remodeler-specific way. We use structure-based mutagenesis to identify and tune the ruler element residing in the Nhp10 and Arp8 modules of the INO80 remodeler complex. Generally, we propose that a remodeler ruler regulates nucleosome sliding direction bias in response to (epi)genetic information. This finally conceptualizes how remodeler-mediated nucleosome dynamics determine stable steady-state nucleosome positioning relative to other nucleosomes, DNA bound factors, DNA ends and DNA sequence elements.

Approaching Protein Barriers: Emerging Mechanisms of Replication Pausing in Eukaryotes

During nuclear DNA replication multiprotein replisome machines have to jointly traverse and duplicate the total length of each chromosome during each cell cycle. At certain genomic locations replisomes encounter tight DNA-protein complexes and slow down. This fork pausing is an active process involving recognition of a protein barrier by the approaching replisome via an evolutionarily conserved Fork Pausing/Protection Complex (FPC). Action of the FPC protects forks from collapse at both programmed and accidental protein barriers, thus promoting genome integrity. In addition, FPC stimulates the DNA replication checkpoint and regulates topological transitions near the replication fork. Eukaryotic cells have been proposed to employ physiological programmed fork pausing for various purposes, such as maintaining copy number at repetitive loci, precluding replication-transcription encounters, regulating kinetochore assembly, or controlling gene conversion events during mating-type switching. Here we review the growing number of approaches used to study replication pausing in vivo and in vitro as well as the characterization of additional factors recently reported to modulate fork pausing in different systems. Specifically, we focus on the positive role of topoisomerases in fork pausing. We describe a model where replisome progression is inherently cautious, which ensures general preservation of fork stability and genome integrity but can also carry out specialized functions at certain loci. Furthermore, we highlight classical and novel outstanding questions in the field and propose venues for addressing them. Given how little is known about replisome pausing at protein barriers in human cells more studies are required to address how conserved these mechanisms are.

Genome information processing by the INO80 chromatin remodeler positions nucleosomes

The fundamental molecular determinants by which ATP-dependent chromatin remodelers organize nucleosomes across eukaryotic genomes remain largely elusive. Here, chromatin reconstitutions on physiological, whole-genome templates reveal how remodelers read and translate genomic information into nucleosome positions. Using the yeast genome and the multi-subunit INO80 remodeler as a paradigm, we identify DNA shape/mechanics encoded signature motifs as sufficient for nucleosome positioning and distinct from known DNA sequence preferences of histones. INO80 processes such information through an allosteric interplay between its core- and Arp8-modules that probes mechanical properties of nucleosomal and linker DNA. At promoters, INO80 integrates this readout of DNA shape/mechanics with a readout of co-evolved sequence motifs via interaction with general regulatory factors bound to these motifs. Our findings establish a molecular mechanism for robust and yet adjustable +1 nucleosome positioning and, more generally, remodelers as information processing hubs that enable active organization and allosteric regulation of the first level of chromatin.

DNA sequence preferences or statistical positioning of histones has not explained genomic patterns of nucleosome organisation in vivo. Here, the authors establish DNA shape/mechanics as key elements that have evolved together with binding sites of DNA sequence-specific barriers so that such information directs nucleosome positioning by chromatin remodelers.

Delving into Eukaryotic Origins of Replication Using DNA Structural Features

DNA replication in eukaryotes is an intricate process, which is precisely synchronized by a set of regulatory proteins, and the replication fork emanates from discrete sites on chromatin called origins of replication (Oris). These spots are considered as the gateway to chromosomal replication and are stereotyped by sequence motifs. The cognate sequences are noticeable in a small group of entire origin regions or totally absent across different metazoans. Alternatively, the use of DNA secondary structural features can provide additional information compared to the primary sequence. In this article, we report the trends in DNA sequence-based structural properties of origin sequences in nine eukaryotic systems representing different families of life. Biologically relevant DNA secondary structural properties, namely, stability, propeller twist, flexibility, and minor groove shape were studied in the sequences flanking replication start sites. Results indicate that Oris in yeasts show lower stability, more rigidity, and narrow minor groove preferences compared to genomic sequences surrounding them. Yeast Oris also show preference for A-tracts and the promoter element TATA box in the vicinity of replication start sites. On the contrary, Drosophila melanogaster, humans, and Arabidopsis thaliana do not have such features in their Oris, and instead, they show high preponderance of G-rich sequence motifs such as putative G-quadruplexes or i-motifs and CpG islands. Our extensive study applies the DNA structural feature computation to delve into origins of replication across organisms ranging from yeasts to mammals and including a plant. Insights from this study would be significant in understanding origin architecture and help in designing new algorithms for predicting DNA trans-acting factor recognition events.

Beads on a string-nucleosome array arrangements and folding of the chromatin fiber

Understanding how the genome is structurally organized as chromatin is essential for understanding its function. Here, we review recent developments that allowed the readdressing of old questions regarding the primary level of chromatin structure, the arrangement of nucleosomes along the DNA and the folding of the nucleosome fiber in nuclear space. In contrast to earlier views of nucleosome arrays as uniformly regular and folded, recent findings reveal heterogeneous array organization and diverse modes of folding. Local structure variations reflect a continuum of functional states characterized by differences in post-translational histone modifications, associated chromatin-interacting proteins and nucleosome-remodeling enzymes.

At the Beginning of the End and in the Middle of the Beginning: Structure and Maintenance of Telomeric DNA Repeats and Interstitial Telomeric Sequences

Tandem DNA repeats derived from the ancestral (TTAGGG)n run were first detected at chromosome ends of the majority of living organisms, hence the name telomeric DNA repeats. Subsequently, it has become clear that telomeric motifs are also present within chromosomes, and they were suitably called interstitial telomeric sequences (ITSs). It is well known that telomeric DNA repeats play a key role in chromosome stability, preventing end-to-end fusions and precluding the recurrent DNA loss during replication. Recent data suggest that ITSs are also important genomic elements as they confer its karyotype plasticity. In fact, ITSs appeared to be among the most unstable microsatellite sequences as they are highly length polymorphic and can trigger chromosomal fragility and gross chromosomal rearrangements. Importantly, mechanisms responsible for their instability appear to be similar to the mechanisms that maintain the length of genuine telomeres. This review compares the mechanisms of maintenance and dynamic properties of telomeric repeats and ITSs and discusses the implications of these dynamics on genome stability.

Systematic Study of Nucleosome-Displacing Factors in Budding Yeast

Nucleosomes present a barrier for the binding of most transcription factors (TFs). However, special TFs known as nucleosome-displacing factors (NDFs) can access embedded sites and cause the depletion of the local nucleosomes as well as repositioning of the neighboring nucleosomes. Here, we developed a novel high-throughput method in yeast to identify NDFs among 104 TFs and systematically characterized the impact of orientation, affinity, location, and copy number of their binding motifs on the nucleosome occupancy. Using this assay, we identified 29 NDF motifs and divided the nuclear TFs into three groups with strong, weak, and no nucleosome-displacing activities. Further studies revealed that tight DNA binding is the key property that underlies NDF activity, and the NDFs may partially rely on the DNA replication to compete with nucleosome. Overall, our study presents a framework to functionally characterize NDFs and elucidate the mechanism of nucleosome invasion.

DMS-Seq for In Vivo Genome-wide Mapping of Protein-DNA Interactions and Nucleosome Centers

Protein-DNA interactions provide the basis for chromatin structure and gene regulation. Comprehensive identification of protein-occupied sites is thus vital to an in-depth understanding of genome function. Dimethyl sulfate (DMS) is a chemical probe that has long been used to detect footprints of DNA-bound proteins in vitro and in vivo. Here, we describe a genomic footprinting method, dimethyl sulfate sequencing (DMS-seq), which exploits the cell-permeable nature of DMS to obviate the need for nuclear isolation. This feature makes DMS-seq simple in practice and removes the potential risk of protein re-localization during nuclear isolation. DMS-seq successfully detects transcription factors bound to cis-regulatory elements and non-canonical chromatin particles in nucleosome-free regions. Furthermore, an unexpected preference of DMS confers on DMS-seq a unique potential to directly detect nucleosome centers without using genetic manipulation. We expect that DMS-seq will serve as a characteristic method for genome-wide interrogation of in vivo protein-DNA interactions.

Conserved noncoding transcription and core promoter regulatory code in early Drosophila development

Multicellular development is driven by regulatory programs that orchestrate the transcription of protein-coding and noncoding genes. To decipher this genomic regulatory code, and to investigate the developmental relevance of noncoding transcription, we compared genome-wide promoter activity throughout embryogenesis in 5 Drosophila species. Core promoters, generally not thought to play a significant regulatory role, in fact impart restrictions on the developmental timing of gene expression on a global scale. We propose a hierarchical regulatory model in which core promoters define broad windows of opportunity for expression, by defining a range of transcription factors from which they can receive regulatory inputs. This two-tiered mechanism globally orchestrates developmental gene expression, including extremely widespread noncoding transcription. The sequence and expression specificity of noncoding RNA promoters are evolutionarily conserved, implying biological relevance. Overall, this work introduces a hierarchical model for developmental gene regulation, and reveals a major role for noncoding transcription in animal development.

Parallel mapping with site-directed hydroxyl radicals and micrococcal nuclease reveals structural features of positioned nucleosomes in vivo

Micrococcal nuclease (MNase) has been widely used for analyses of nucleosome locations in many organisms. However, due to its sequence preference, the interpretations of the positions and occupancies of nucleosomes using MNase have remained controversial. Next-generation sequencing (NGS) has also been utilized for analyses of MNase-digests, but some technical biases are commonly present in the NGS experiments. Here, we established a gel-based method to map nucleosome positions in Saccharomyces cerevisiae, using isolated nuclei as the substrate for the histone H4 S47C-site-directed chemical cleavage in parallel with MNase digestion. The parallel mapping allowed us to compare the chemically and enzymatically cleaved sites by indirect end-labeling and primer extension mapping, and thus we could determine the nucleosome positions and the sizes of the nucleosome-free regions (or nucleosome-depleted regions) more accurately, as compared to nucleosome mapping by MNase alone. The analysis also revealed that the structural features of the nucleosomes flanked by the nucleosome-free region were different from those within regularly arrayed nucleosomes, showing that the structures and dynamics of individual nucleosomes strongly depend on their locations. Moreover, we demonstrated that the parallel mapping results were generally consistent with the previous genome-wide chemical mapping and MNase-Seq results. Thus, the gel-based parallel mapping will be useful for the analysis of a specific locus under various conditions.

Shelterin components mediate genome reorganization in response to replication stress

The dynamic nature of genome organization impacts critical nuclear functions including the regulation of gene expression, replication, and DNA damage repair. Despite significant progress, the mechanisms responsible for reorganization of the genome in response to cellular stress, such as aberrant DNA replication, are poorly understood. Here, we show that fission yeast cells carrying a mutation in the DNA-binding protein Sap1 show defects in DNA replication progression and genome stability and display extensive changes in genome organization. Chromosomal regions such as subtelomeres that show defects in replication progression associate with the nuclear envelope in sap1 mutant cells. Moreover, high-resolution, genome-wide chromosome conformation capture (Hi-C) analysis revealed prominent contacts between telomeres and chromosomal arm regions containing replication origins proximal to binding sites for Taz1, a component of the Shelterin telomere protection complex. Strikingly, we find that Shelterin components are required for interactions between Taz1-associated chromosomal arm regions and telomeres. These analyses reveal an unexpected role for Shelterin components in genome reorganization in cells experiencing replication stress, with important implications for understanding the mechanisms governing replication and genome stability.

Cross-Regulation between Transposable Elements and Host DNA Replication

Transposable elements subvert host cellular functions to ensure their survival. Their interaction with the host DNA replication machinery indicates that selective pressures lead them to develop ancestral and convergent evolutionary adaptations aimed at conserved features of this fundamental process. These interactions can shape the co-evolution of the transposons and their hosts.

Regulation of DNA Replication through Natural Impediments in the Eukaryotic Genome

All living organisms need to duplicate their genetic information while protecting it from unwanted mutations, which can lead to genetic disorders and cancer development. Inaccuracies during DNA replication are the major cause of genomic instability, as replication forks are prone to stalling and collapse, resulting in DNA damage. The presence of exogenous DNA damaging agents as well as endogenous difficult-to-replicate DNA regions containing DNA–protein complexes, repetitive DNA, secondary DNA structures, or transcribing RNA polymerases, increases the risk of genomic instability and thus threatens cell survival. Therefore, understanding the cellular mechanisms required to preserve the genetic information during S phase is of paramount importance. In this review, we will discuss our current understanding of how cells cope with these natural impediments in order to prevent DNA damage and genomic instability during DNA replication.

Nucleosomal signatures impose nucleosome positioning in coding and noncoding sequences in the genome

In the yeast genome, a large proportion of nucleosomes occupy well-defined and stable positions. While the contribution of chromatin remodelers and DNA binding proteins to maintain this organization is well established, the relevance of the DNA sequence to nucleosome positioning in the genome remains controversial. Through quantitative analysis of nucleosome positioning, we show that sequence changes distort the nucleosomal pattern at the level of individual nucleosomes in three species of Schizosaccharomyces and in Saccharomyces cerevisiae. This effect is equally detected in transcribed and nontranscribed regions, suggesting the existence of sequence elements that contribute to positioning. To identify such elements, we incorporated information from nucleosomal signatures into artificial synthetic DNA molecules and found that they generated regular nucleosomal arrays indistinguishable from those of endogenous sequences. Strikingly, this information is species-specific and can be combined with coding information through the use of synonymous codons such that genes from one species can be engineered to adopt the nucleosomal organization of another. These findings open the possibility of designing coding and noncoding DNA molecules capable of directing their own nucleosomal organization.

Transcription factor Reb1 is required for proper transcriptional start site usage at the divergently transcribed TFC6-ESC2 locus in Saccharomyces cerevisiae

Eukaryotic promoters generally contain nucleosome depleted regions near their transcription start sites. In the model organism Saccharomyces cerevisiae, these regions are adjacent to binding sites for general regulatory transcription factors, and the Reb1 protein is commonly bound to promoter DNA near such regions. The yeast TFC6 promoter is a unique RNA polymerase II promoter in that it is autoregulated by its own gene product Tfc6p, which is part of the RNA polymerase III transcription factor complex TFIIIC. We previously demonstrated that mutation of a potential Reb1 binding site adjacent to the TFIIIC binding site in the TFC6 promoter modestly reduces transcript levels, but leads to a severe decrease in Tfc6 protein levels due to an upstream shift in the TFC6 transcription start site. Here we confirm that Reb1p indeed binds to the TFC6 promoter, and is important for proper transcription start site selection and protein expression. Interestingly, loss of Reb1p association at this site has a similar effect on the adjacent divergently transcribed ESC2 promoter, resulting in a significant increase of 5'-extended ESC2 transcripts and reduction of Esc2 protein levels. This altered divergent transcription may be the result of changes in nucleosome positioning at this locus in the absence of Reb1p binding. We speculate that an important function of general regulatory factors such as Reb1p is to establish and maintain proper transcription start sites at promoters, and that when binding of such factors is compromised, resulting effects on mRNA translation may be an underappreciated aspect of gene regulation studies.

Evidence of selection for an accessible nucleosomal array in human

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Nucleosome Stability Distinguishes Two Different Promoter Types at All Protein-Coding Genes in Yeast

Previous studies indicate that eukaryotic promoters display a stereotypical chromatin landscape characterized by a well-positioned +1 nucleosome near the transcription start site and an upstream -1 nucleosome that together demarcate a nucleosome-free (or -depleted) region. Here we present evidence that there are two distinct types of promoters distinguished by the resistance of the -1 nucleosome to micrococcal nuclease digestion. These different architectures are characterized by two sequence motifs that are broadly deployed at one set of promoters where a nuclease-sensitive ("fragile") nucleosome forms, but concentrated in a narrower, nucleosome-free region at all other promoters. The RSC nucleosome remodeler acts through the motifs to establish stable +1 and -1 nucleosome positions, while binding of a small set of general regulatory (pioneer) factors at fragile nucleosome promoters plays a key role in their destabilization. We propose that the fragile nucleosome promoter architecture is adapted for regulation of highly expressed, growth-related genes.

A species-specific nucleosomal signature defines a periodic distribution of amino acids in proteins

Nucleosomes are the basic structural units of chromatin. Most of the yeast genome is organized in a pattern of positioned nucleosomes that is stably maintained under a wide range of physiological conditions. In this work, we have searched for sequence determinants associated with positioned nucleosomes in four species of fission and budding yeasts. We show that mononucleosomal DNA follows a highly structured base composition pattern, which differs among species despite the high degree of histone conservation. These nucleosomal signatures are present in transcribed and non-transcribed regions across the genome. In the case of open reading frames, they correctly predict the relative distribution of codons on mononucleosomal DNA, and they also determine a periodicity in the average distribution of amino acids along the proteins. These results establish a direct and species-specific connection between the position of each codon around the histone octamer and protein composition.

Crop epigenetics and the molecular hardware of genotype × environment interactions

Crop plants encounter thermal environments which fluctuate on a diurnal and seasonal basis. Future climate resilient cultivars will need to respond to thermal profiles reflecting more variable conditions, and harness plasticity that involves regulation of epigenetic processes and complex genomic regulatory networks. Compartmentalization within plant cells insulates the genomic central processing unit within the interphase nucleus. This review addresses the properties of the chromatin hardware in which the genome is embedded, focusing on the biophysical and thermodynamic properties of DNA, histones and nucleosomes. It explores the consequences of thermal and ionic variation on the biophysical behavior of epigenetic marks such as DNA cytosine methylation (5mC), and histone variants such as H2A.Z, and how these contribute to maintenance of chromatin integrity in the nucleus, while enabling specific subsets of genes to be regulated. Information is drawn from theoretical molecular in vitro studies as well as model and crop plants and incorporates recent insights into the role epigenetic processes play in mediating between environmental signals and genomic regulation. A preliminary speculative framework is outlined, based on the evidence of what appears to be a cohesive set of interactions at molecular, biophysical and electrostatic level between the various components contributing to chromatin conformation and dynamics. It proposes that within plant nuclei, general and localized ionic homeostasis plays an important role in maintaining chromatin conformation, whilst maintaining complex genomic regulation that involves specific patterns of epigenetic marks. More generally, reversible changes in DNA methylation appear to be consistent with the ability of nuclear chromatin to manage variation in external ionic and temperature environment. Whilst tentative, this framework provides scope to develop experimental approaches to understand in greater detail the internal environment of plant nuclei. It is hoped that this will generate a deeper understanding of the molecular mechanisms underlying genotype × environment interactions that may be beneficial for long-term improvement of crop performance in less predictable climates.

Comparative Transcriptomics in Yeasts

Comparative functional genomics approaches have already shed an important light on the evolution of gene expression that underlies phenotypic diversity. However, comparison across many species in a phylogeny presents several major challenges. Here, we describe our experimental framework for comparative transcriptomics in a complex phylogeny.

Arrested replication forks guide retrotransposon integration

Long terminal repeat (LTR) retrotransposons are an abundant class of genomic parasites that replicate by insertion of new copies into the host genome. Fungal LTR retrotransposons prevent mutagenic insertions through diverse targeting mechanisms that avoid coding sequences, but conserved principles guiding their target site selection have not been established. Here, we show that insertion of the fission yeast LTR retrotransposon Tf1 is guided by the DNA binding protein Sap1 and that the efficiency and location of the targeting depend on the activity of Sap1 as a replication fork barrier. We propose that Sap1 and the fork arrest it causes guide insertion of Tf1 by tethering the integration complex to target sites.

Nucleosome positioning: resources and tools online

Nucleosome positioning is an important process required for proper genome packing and its accessibility to execute the genetic program in a cell-specific, timely manner. In the recent years hundreds of papers have been devoted to the bioinformatics, physics and biology of nucleosome positioning. The purpose of this review is to cover a practical aspect of this field, namely, to provide a guide to the multitude of nucleosome positioning resources available online. These include almost 300 experimental datasets of genome-wide nucleosome occupancy profiles determined in different cell types and more than 40 computational tools for the analysis of experimental nucleosome positioning data and prediction of intrinsic nucleosome formation probabilities from the DNA sequence. A manually curated, up to date list of these resources will be maintained at http://generegulation.info.

High-Throughput Identification of Cis-Regulatory Rewiring Events in Yeast

A coregulated module of genes ("regulon") can have evolutionarily conserved expression patterns and yet have diverged upstream regulators across species. For instance, the ribosomal genes regulon is regulated by the transcription factor (TF) TBF1 in Candida albicans, while in Saccharomyces cerevisiae it is regulated by RAP1. Only a handful of such rewiring events have been established, and the prevalence or conditions conducive to such events are not well known. Here, we develop a novel probabilistic scoring method to comprehensively screen for regulatory rewiring within regulons across 23 yeast species. Investigation of 1,713 regulons and 176 TFs yielded 5,353 significant rewiring events at 5% false discovery rate (FDR). Besides successfully recapitulating known rewiring events, our analyses also suggest TF candidates for certain processes reported to be under distinct regulatory controls in S. cerevisiae and C. albicans, for which the implied regulators are not known: 1) Oxidative stress response (Sc-MSN2 to Ca-FKH2) and 2) nutrient modulation (Sc-RTG1 to Ca-GCN4/Ca-UME6). Furthermore, a stringent screen to detect TF rewiring at individual genes identified 1,446 events at 10% FDR. Overall, these events are supported by strong coexpression between the predicted regulator and its target gene(s) in a species-specific fashion (>50-fold). Independent functional analyses of rewiring TF pairs revealed greater functional interactions and shared biological processes between them (P = 1 × 10(-3)).Our study represents the first comprehensive assessment of regulatory rewiring; with a novel approach that has generated a unique high-confidence resource of several specific events, suggesting that evolutionary rewiring is relatively frequent and may be a significant mechanism of regulatory innovation.

Single-Nucleotide-Specific Targeting of the Tf1 Retrotransposon Promoted by the DNA-Binding Protein Sap1 of Schizosaccharomyces pombe

Transposable elements (TEs) constitute a substantial fraction of the eukaryotic genome and, as a result, have a complex relationship with their host that is both adversarial and dependent. To minimize damage to cellular genes, TEs possess mechanisms that target integration to sequences of low importance. However, the retrotransposon Tf1 of Schizosaccharomyces pombe integrates with a surprising bias for promoter sequences of stress-response genes. The clustering of integration in specific promoters suggests that Tf1 possesses a targeting mechanism that is important for evolutionary adaptation to changes in environment. We report here that Sap1, an essential DNA-binding protein, plays an important role in Tf1 integration. A mutation in Sap1 resulted in a 10-fold drop in Tf1 transposition, and measures of transposon intermediates support the argument that the defect occurred in the process of integration. Published ChIP-Seq data on Sap1 binding combined with high-density maps of Tf1 integration that measure independent insertions at single-nucleotide positions show that 73.4% of all integration occurs at genomic sequences bound by Sap1. This represents high selectivity because Sap1 binds just 6.8% of the genome. A genome-wide analysis of promoter sequences revealed that Sap1 binding and amounts of integration correlate strongly. More important, an alignment of the DNA-binding motif of Sap1 revealed integration clustered on both sides of the motif and showed high levels specifically at positions +19 and -9. These data indicate that Sap1 contributes to the efficiency and position of Tf1 integration.

Comparative analysis of gene regulatory networks: from network reconstruction to evolution

Regulation of gene expression is central to many biological processes. Although reconstruction of regulatory circuits from genomic data alone is therefore desirable, this remains a major computational challenge. Comparative approaches that examine the conservation and divergence of circuits and their components across strains and species can help reconstruct circuits as well as provide insights into the evolution of gene regulatory processes and their adaptive contribution. In recent years, advances in genomic and computational tools have led to a wealth of methods for such analysis at the sequence, expression, pathway, module, and entire network level. Here, we review computational methods developed to study transcriptional regulatory networks using comparative genomics, from sequence to functional data. We highlight how these methods use evolutionary conservation and divergence to reliably detect regulatory components as well as estimate the extent and rate of divergence. Finally, we discuss the promise and open challenges in linking regulatory divergence to phenotypic divergence and adaptation.

Structured nucleosome fingerprints enable high-resolution mapping of chromatin architecture within regulatory regions

Transcription factors canonically bind nucleosome-free DNA, making the positioning of nucleosomes within regulatory regions crucial to the regulation of gene expression. Using the assay of transposase accessible chromatin (ATAC-seq), we observe a highly structured pattern of DNA fragment lengths and positions around nucleosomes in Saccharomyces cerevisiae, and use this distinctive two-dimensional nucleosomal “fingerprint” as the basis for a new nucleosome-positioning algorithm called NucleoATAC. We show that NucleoATAC can identify the rotational and translational positions of nucleosomes with up to base-pair resolution and provide quantitative measures of nucleosome occupancy in S. cerevisiae, Schizosaccharomyces pombe, and human cells. We demonstrate the application of NucleoATAC to a number of outstanding problems in chromatin biology, including analysis of sequence features underlying nucleosome positioning, promoter chromatin architecture across species, identification of transient changes in nucleosome occupancy and positioning during a dynamic cellular response, and integrated analysis of nucleosome occupancy and transcription factor binding.

Regulation of Ace2-dependent genes requires components of the PBF complex in Schizosaccharomyces pombe

The division cycle of unicellular yeasts is completed with the activation of a cell separation program that results in the dissolution of the septum assembled during cytokinesis between the 2 daughter cells, allowing them to become independent entities. Expression of the eng1(+) and agn1(+) genes, encoding the hydrolytic enzymes responsible for septum degradation, is activated at the end of each cell cycle by the transcription factor Ace2. Periodic ace2(+) expression is regulated by the transcriptional complex PBF (PCB Binding Factor), composed of the forkhead-like proteins Sep1 and Fkh2 and the MADS box-like protein Mbx1. In this report, we show that Ace2-dependent genes contain several combinations of motifs for Ace2 and PBF binding in their promoters. Thus, Ace2, Fkh2 and Sep1 were found to bind in vivo to the eng1(+) promoter. Ace2 binding was coincident with maximum level of eng1(+) expression, whereas Fkh2 binding was maximal when mRNA levels were low, supporting the notion that they play opposing roles. In addition, we found that the expression of eng1(+) and agn1(+) was differentially affected by mutations in PBF components. Interestingly, agn1(+) was a major target of Mbx1, since its ectopic expression resulted in the suppression of Mbx1 deletion phenotypes. Our results reveal a complex regulation system through which the transcription factors Ace2, Fkh2, Sep1 and Mbx1 in combination control the expression of the genes involved in separation at the end of the cell division cycle.

Comparative Genomics Reveals Chd1 as a Determinant of Nucleosome Spacing in Vivo

Packaging of genomic DNA into nucleosomes is nearly universally conserved in eukaryotes, and many features of the nucleosome landscape are quite conserved. Nonetheless, quantitative aspects of nucleosome packaging differ between species because, for example, the average length of linker DNA between nucleosomes can differ significantly even between closely related species. We recently showed that the difference in nucleosome spacing between two Hemiascomycete species—Saccharomyces cerevisiae and Kluyveromyces lactis—is established by trans-acting factors rather than being encoded in cis in the DNA sequence. Here, we generated several S. cerevisiae strains in which endogenous copies of candidate nucleosome spacing factors are deleted and replaced with the orthologous factors from K. lactis. We find no change in nucleosome spacing in such strains in which H1 or Isw1 complexes are swapped. In contrast, the K. lactis gene encoding the ATP-dependent remodeler Chd1 was found to direct longer internucleosomal spacing in S. cerevisiae, establishing that this remodeler is partially responsible for the relatively long internucleosomal spacing observed in K. lactis. By analyzing several chimeric proteins, we find that sequence differences that contribute to the spacing activity of this remodeler are dispersed throughout the coding sequence, but that the strongest spacing effect is linked to the understudied N-terminal end of Chd1. Taken together, our data find a role for sequence evolution of a chromatin remodeler in establishing quantitative aspects of the chromatin landscape in a species-specific manner.

Control of relative timing and stoichiometry by a master regulator

Developmental processes in cells require a series of complex steps. Often only a single master regulator activates genes in these different steps. This poses several challenges: some targets need to be ordered temporally, while co-functional targets may need to be synchronized in both time and expression level. Here we study in single cells the dynamic activation patterns of early meiosis genes in budding yeast, targets of the meiosis master regulator Ime1. We quantify the individual roles of the promoter and protein levels in expression pattern control, as well as the roles of individual promoter elements. We find a consistent expression pattern difference between a non-cofunctional pair of genes, and a highly synchronized activation of a co-functional pair. We show that dynamic control leading to these patterns is distributed between promoter, gene and external regions. Through specific reciprocal changes to the promoters of pairs of genes, we show that different genes can use different promoter elements to reach near identical activation patterns.

The chromatin remodelers RSC and ISW1 display functional and chromatin-based promoter antagonism

ISWI family chromatin remodelers typically organize nucleosome arrays, while SWI/SNF family remodelers (RSC) typically disorganize and eject nucleosomes, implying an antagonism that is largely unexplored in vivo. Here, we describe two independent genetic screens for rsc suppressors that yielded mutations in the promoter-focused ISW1a complex or mutations in the ‘basic patch’ of histone H4 (an epitope that regulates ISWI activity), strongly supporting RSC-ISW1a antagonism in vivo. RSC and ISW1a largely co-localize, and genomic nucleosome studies using rsc isw1 mutant combinations revealed opposing functions: promoters classified with a nucleosome-deficient region (NDR) gain nucleosome occupancy in rsc mutants, but this gain is attenuated in rsc isw1 double mutants. Furthermore, promoters lacking NDRs have the highest occupancy of both remodelers, consistent with regulation by nucleosome occupancy, and decreased transcription in rsc mutants. Taken together, we provide the first genetic and genomic evidence for RSC-ISW1a antagonism and reveal different mechanisms at two different promoter architectures.

DOI:http://dx.doi.org/10.7554/eLife.06073.001

The genome of an organism can contain hundreds to thousands of genes. However, these genes are not all active at the same time. Instead, mechanisms exist that control when genes are switched off or on. One such mechanism alters how tightly DNA is packaged into a structure called chromatin. To form chromatin, DNA is wrapped around histone proteins at different points along its length; these structures are known as nucleosomes. Once formed, chromatin can either adopt a tightly packed form that represses gene activity or a less compact form associated with gene activation.

The proteins that control how chromatin is packed are called ‘chromatin remodelers’. These proteins work in complexes that either disassemble chromatin—for example, by repositioning nucleosomes or removing histones—or organize chromatin by replacing nucleosomes and making it more compact.

Studies in many organisms have identified two key families of chromatin remodelers. In yeast, the ISWI family of complexes assembles chromatin and a protein complex called RSC disassembles chromatin. Parnell, Schlichter et al. used a range of genetic techniques to investigate whether these two chromatin-remodeling complexes work against each other in a species of yeast called Saccharomyces cerevisiae. The results suggest that this is indeed the case.

Both the ISWI complex and the RSC complex bind to regions of DNA called promoters, which are found at the start of a gene and help to regulate its activity. Parnell, Schlichter et al. found that the RSC complex helps to activate genes by establishing or maintaining regions of nucleosome-poor chromatin at a promoter. The chromatin is relaxed in these regions, which allows the proteins that activate genes to access the DNA. This effect is partially counteracted by the ISWI complex, which repositions nucleosomes across the promoters to fill the gaps created by the RSC complex.

In comparison, Parnell, Schlichter et al. found that promoters that do not have regions of nucleosome-poor chromatin contain a larger number of both of the remodeling complexes and have a high turnover of histone proteins. This suggests that at these sites, the RSC proteins are needed to overcome the assembly of nucleosomes by the ISWI complex in order to activate the gene. Thus, these two chromatin remodelers, ISWI and RSC, compete at promoters to determine whether they contain or lack nucleosomes, which helps determine whether the gene is silent or active, respectively. Future work will look further at how the ‘winner’ is determined: how transcription factors or signaling systems within the cell bias the recruitment or activity of RSC or ISWI at particular promoters, to determine gene activity.

DOI:http://dx.doi.org/10.7554/eLife.06073.002

Ubiquitous human 'master' origins of replication are encoded in the DNA sequence via a local enrichment in nucleosome excluding energy barriers

As the elementary building block of eukaryotic chromatin, the nucleosome is at the heart of the compromise between the necessity of compacting DNA in the cell nucleus and the required accessibility to regulatory proteins. The recent availability of genome-wide experimental maps of nucleosome positions for many different organisms and cell types has provided an unprecedented opportunity to elucidate to what extent the DNA sequence conditions the primary structure of chromatin and in turn participates in the chromatin-mediated regulation of nuclear functions, such as gene expression and DNA replication. In this study, we use in vivo and in vitro genome-wide nucleosome occupancy data together with the set of nucleosome-free regions (NFRs) predicted by a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix, to investigate the role of intrinsic nucleosome occupancy in the regulation of the replication spatio-temporal programme in human. We focus our analysis on the so-called replication U/N-domains that were shown to cover about half of the human genome in the germline (skew-N domains) as well as in embryonic stem cells, somatic and HeLa cells (mean replication timing U-domains). The 'master' origins of replication (MaOris) that border these megabase-sized U/N-domains were found to be specified by a few hundred kb wide regions that are hyper-sensitive to DNase I cleavage, hypomethylated, and enriched in epigenetic marks involved in transcription regulation, the hallmarks of localized open chromatin structures. Here we show that replication U/N-domain borders that are conserved in all considered cell lines have an environment highly enriched in nucleosome-excluding-energy barriers, suggesting that these ubiquitous MaOris have been selected during evolution. In contrast, MaOris that are cell-type-specific are mainly regulated epigenetically and are no longer favoured by a local abundance of intrinsic NFRs encoded in the DNA sequence. At the smaller few hundred bp scale of gene promoters, CpG-rich promoters of housekeeping genes found nearby ubiquitous MaOris as well as CpG-poor promoters of tissue-specific genes found nearby cell-type-specific MaOris, both correspond to in vivo NFRs that are not coded as nucleosome-excluding-energy barriers. Whereas the former promoters are likely to correspond to high occupancy transcription factor binding regions, the latter are an illustration that gene regulation in human is typically cell-type-specific.

Nucleosome positioning in yeasts: methods, maps, and mechanisms

Eukaryotic nuclear DNA is packaged into nucleosomes. During the past decade, genome-wide nucleosome mapping across species revealed the high degree of order in nucleosome positioning. There is a conserved stereotypical nucleosome organization around transcription start sites (TSSs) with a nucleosome-depleted region (NDR) upstream of the TSS and a TSS-aligned regular array of evenly spaced nucleosomes downstream over the gene body. As nucleosomes largely impede access to DNA and thereby provide an important level of genome regulation, it is of general interest to understand the mechanisms generating nucleosome positioning and especially the stereotypical NDR-array pattern. We focus here on the most advanced models, unicellular yeasts, and review the progress in mapping nucleosomes and which nucleosome positioning mechanisms are discussed. There are four mechanistic aspects: How are NDRs generated? How are individual nucleosomes positioned, especially those flanking the NDRs? How are nucleosomes evenly spaced leading to regular arrays? How are regular arrays aligned at TSSs? The main candidates for nucleosome positioning determinants are intrinsic DNA binding preferences of the histone octamer, specific DNA binding factors, nucleosome remodeling enzymes, transcription, and statistical positioning. We summarize the state of the art in an integrative model where nucleosomes are positioned by a combination of all these candidate determinants. We highlight the predominance of active mechanisms involving nucleosome remodeling enzymes which may be recruited by DNA binding factors and the transcription machinery. While this mechanistic framework emerged clearly during recent years, the involved factors and their mechanisms are still poorly understood and require future efforts combining in vivo and in vitro approaches.

Transcription mediated insulation and interference direct gene cluster expression switches

In yeast, many tandemly arranged genes show peak expression in different phases of the metabolic cycle (YMC) or in different carbon sources, indicative of regulation by a bi-modal switch, but it is not clear how these switches are controlled. Using native elongating transcript analysis (NET-seq), we show that transcription itself is a component of bi-modal switches, facilitating reciprocal expression in gene clusters. HMS2, encoding a growth-regulated transcription factor, switches between sense- or antisense-dominant states that also coordinate up- and down-regulation of transcription at neighbouring genes. Engineering HMS2 reveals alternative mono-, di- or tri-cistronic and antisense transcription units (TUs), using different promoter and terminator combinations, that underlie state-switching. Promoters or terminators are excluded from functional TUs by read-through transcriptional interference, while antisense TUs insulate downstream genes from interference. We propose that the balance of transcriptional insulation and interference at gene clusters facilitates gene expression switches during intracellular and extracellular environmental change.

DOI:http://dx.doi.org/10.7554/eLife.03635.001

A DNA double helix is made up of two DNA strands, which in turn are made of molecules that are each known by a single letter—A, T, C, or G. The sequence of these ‘letters’ in each DNA strand contains biological information.

Genes are sections of DNA that can be ‘expressed’ to produce proteins and RNA molecules. To express a gene, the DNA strands in the double helix must first be partially separated so that one of them can be used as a template to build an RNA molecule in a process called transcription. Either of the DNA strands in a helix can be used as an RNA template, but contain different genes and are read in opposite directions. One of the two strands is called the ‘sense’ strand, the other the ‘antisense’ strand.

The RNA molecule does not transcribe a whole DNA strand; instead, it transcribes a section of DNA, known as a transcription unit, which contains at least one gene. The end of a transcription unit is marked by certain signals that stop transcription. However, some transcription units in a DNA strand overlap, so there must be some way that the transcription machinery can sometimes ignore these stop signals.

The activity of some genes is linked to the activity of their immediate neighbours. Furthermore, some genes are expressed in different amounts in response to changes in environmental conditions. Researchers have previously suggested that there must be some form of switch that controls when these genes are expressed.

Nguyen et al. now engineer start and stop signals at a neighbouring pair of genes, called HMS2 and BAT2, in yeast. When one gene is switched on, the other is switched off and which gene is active depends on the diet of the yeast cells.

On the antisense DNA strand opposite to HMS2 is another gene, SUT650. Nguyen et al. show that when this gene is transcribed, the transcription of HMS2 on the other DNA strand is blocked. This has the knock-on effect of turning on BAT2. Conversely, transcribing HMS2 switches off SUT650 and BAT2 because the end of HMS2 overlaps with the beginning of both SUT650 and BAT2. Switching between different genes relies on loops that physically link the start and stop signals of the gene to be transcribed while ignoring the start and stop signals for neighbouring genes.

Proteins called transcription factors can bind to DNA and affect whether a gene is transcribed. Nguyen et al. found that a transcription factor that binds near the start of the HMS2 gene helps to control which DNA strand is transcribed. When transcription factors do not bind to the start of HMS2, antisense transcription—and the expression of SUT650—occurs instead.

Overall, Nguyen et al. show that the transcription process itself makes up part of a switch that can control the expression of several genes on both the sense and antisense strands of a DNA double helix. This may also explain how many other, more complex, gene networks are activated in response to changes in the environment.

DOI:http://dx.doi.org/10.7554/eLife.03635.002

Comparative analysis of methods for genome-wide nucleosome cartography

Nucleosomes contribute to compacting the genome into the nucleus and regulate the physical access of regulatory proteins to DNA either directly or through the epigenetic modifications of the histone tails. Precise mapping of nucleosome positioning across the genome is, therefore, essential to understanding the genome regulation. In recent years, several experimental protocols have been developed for this purpose that include the enzymatic digestion, chemical cleavage or immunoprecipitation of chromatin followed by next-generation sequencing of the resulting DNA fragments. Here, we compare the performance and resolution of these methods from the initial biochemical steps through the alignment of the millions of short-sequence reads to a reference genome to the final computational analysis to generate genome-wide maps of nucleosome occupancy. Because of the lack of a unified protocol to process data sets obtained through the different approaches, we have developed a new computational tool (NUCwave), which facilitates their analysis, comparison and assessment and will enable researchers to choose the most suitable method for any particular purpose. NUCwave is freely available at http://nucleosome.usal.es/nucwave along with a step-by-step protocol for its use.

A snapshot of Snf2 enzymes in fission yeast

Eukaryotic chromatin is remodelled by the evolutionarily conserved Snf2 family of enzymes in an ATP-dependent manner. Several Snf2 enzymes are part of CRCs (chromatin remodelling complexes). In the present review we focus our attention on the functions of Snf2 enzymes and CRCs in fission yeast. We discuss their molecular mechanisms and roles and in regulating gene expression, DNA recombination, euchromatin and heterochromatin structure.

Purification of G1 daughter cells from different Saccharomycetes species through an optimized centrifugal elutriation procedure

Centrifugal elutriation discriminates cells according to their sedimentation coefficients, generating homogeneous samples well suited for genomic comparative approaches. It can, for instance, isolate G1 daughter cells from a Saccharomyces cerevisiae unsynchronized population, alleviating ageing and cell-cycle biases when conducting genome-wide/single-cell studies. The present report describes a straightforward and robust procedure to determine whether a cell population of virtually any yeast species can be efficiently elutriated, while offering solutions to optimize success. This approach was used to characterize elutriation parameters and S-phase progression of four yeast species (S. cerevisiae, Candida glabrata, Lachancea kluyveri and Pichia sorbitophila) and could theoretically be applied to any culture of single, individual cells.

Mapping yeast transcriptional networks

The term “transcriptional network” refers to the mechanism(s) that underlies coordinated expression of genes, typically involving transcription factors (TFs) binding to the promoters of multiple genes, and individual genes controlled by multiple TFs. A multitude of studies in the last two decades have aimed to map and characterize transcriptional networks in the yeast Saccharomyces cerevisiae. We review the methodologies and accomplishments of these studies, as well as challenges we now face. For most yeast TFs, data have been collected on their sequence preferences, in vivo promoter occupancy, and gene expression profiles in deletion mutants. These systematic studies have led to the identification of new regulators of numerous cellular functions and shed light on the overall organization of yeast gene regulation. However, many yeast TFs appear to be inactive under standard laboratory growth conditions, and many of the available data were collected using techniques that have since been improved. Perhaps as a consequence, comprehensive and accurate mapping among TF sequence preferences, promoter binding, and gene expression remains an open challenge. We propose that the time is ripe for renewed systematic efforts toward a complete mapping of yeast transcriptional regulatory mechanisms.

Fuzziness and noise in nucleosomal architecture

Nucleosome organization plays a key role in the regulation of gene expression. However, despite the striking advances in the accuracy of nucleosome maps, there are still severe discrepancies on individual nucleosome positioning and how this influences gene regulation. The variability among nucleosome maps, which precludes the fine analysis of nucleosome positioning, might emerge from diverse sources. We have carefully inspected the extrinsic factors that may induce diversity by the comparison of microccocal nuclease (MNase)-Seq derived nucleosome maps generated under distinct conditions. Furthermore, we have also explored the variation originated from intrinsic nucleosome dynamics by generating additional maps derived from cell cycle synchronized and asynchronous yeast cultures. Taken together, our study has enabled us to measure the effect of noise in nucleosome occupancy and positioning and provides insights into the underlying determinants. Furthermore, we present a systematic approach that may guide the standardization of MNase-Seq experiments in order to generate reproducible genome-wide nucleosome patterns.

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

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.

Clustered regulatory elements at nucleosome-depleted regions punctuate a constant nucleosomal landscape in Schizosaccharomyces pombe

BACKGROUND

Nucleosomes facilitate the packaging of the eukaryotic genome and modulate the access of regulators to DNA. A detailed description of the nucleosomal organization under different transcriptional programmes is essential to understand their contribution to genomic regulation.

RESULTS

To visualize the dynamics of individual nucleosomes under different transcriptional programmes we have generated high-resolution nucleosomal maps in Schizosaccharomyces pombe. We show that 98.5% of the genome remains almost invariable during mitosis and meiosis while remodelling is limited to approximately 1100 nucleosomes in the promoters of a subset of meiotic genes. These inducible nucleosome-depleted regions (NDR) and also those constitutively present in the genome overlap precisely with clusters of binding sites for transcription factors (TF) specific for meiosis and for different functional classes of genes, respectively. Deletion of two TFs affects only a small fraction of all the NDRs to which they bind in vivo, indicating that TFs collectively contribute to NDR maintenance.

CONCLUSIONS

Our results show that the nucleosomal profile in S. pombe is largely maintained under different physiological conditions and patterns of gene expression. This relatively constant landscape favours the concentration of regulators in constitutive and inducible NDRs. The combinatorial analysis of binding motifs in this discrete fraction of the genome will facilitate the definition of the transcriptional regulatory networks.

DNA replication origins

The onset of genomic DNA synthesis requires precise interactions of specialized initiator proteins with DNA at sites where the replication machinery can be loaded. These sites, defined as replication origins, are found at a few unique locations in all of the prokaryotic chromosomes examined so far. However, replication origins are dispersed among tens of thousands of loci in metazoan chromosomes, thereby raising questions regarding the role of specific nucleotide sequences and chromatin environment in origin selection and the mechanisms used by initiators to recognize replication origins. Close examination of bacterial and archaeal replication origins reveals an array of DNA sequence motifs that position individual initiator protein molecules and promote initiator oligomerization on origin DNA. Conversely, the need for specific recognition sequences in eukaryotic replication origins is relaxed. In fact, the primary rule for origin selection appears to be flexibility, a feature that is modulated either by structural elements or by epigenetic mechanisms at least partly linked to the organization of the genome for gene expression.

Characteristics of nucleosomes and linker DNA regions on the genome of the basidiomycete Mixia osmundae revealed by mono- and dinucleosome mapping

We present findings on the nucleosomal arrangement in the genome of the basidiomycete Mixia osmundae, focusing on nucleosomal linker DNA regions. We have assembled the genomic sequences of M. osmundae, annotated genes and transcription start sites (TSSs) on the genome, and created a detailed nucleosome map based on sequencing mono- and dinucleosomal DNA fragments. The nucleosomal DNA length distribution of M. osmundae is similar to that of the filamentous ascomycete Aspergillus fumigatus, but differs from that of ascomycetous yeasts, strongly suggesting that nucleosome positioning has evolved primarily through neutral drift in fungal species. We found clear association between dinucleotide frequencies and linker DNA regions mapped as the midpoints of dinucleosomes. We also describe a unique pattern found in the nucleosome-depleted region upstream of the TSS observed in the dinucleosome map and the precursor status of dinucleosomes prior to the digestion into mononucleosomes by comparing the mono- and dinucleosome maps. We demonstrate that observation of dinucleosomes as well as of mononucleosomes is valuable in investigating nucleosomal organization of the genome.

Yeast rises to the occasion

Genetic analyses of 15 species of yeast have shed new light on the divergence of gene regulation during evolution, with significant changes occurring after an event in which a whole genome was duplicated.

Evolutionary principles of modular gene regulation in yeasts

Divergence in gene regulation can play a major role in evolution. Here, we used a phylogenetic framework to measure mRNA profiles in 15 yeast species from the phylum Ascomycota and reconstruct the evolution of their modular regulatory programs along a time course of growth on glucose over 300 million years. We found that modules have diverged proportionally to phylogenetic distance, with prominent changes in gene regulation accompanying changes in lifestyle and ploidy, especially in carbon metabolism. Paralogs have significantly contributed to regulatory divergence, typically within a very short window from their duplication. Paralogs from a whole genome duplication (WGD) event have a uniquely substantial contribution that extends over a longer span. Similar patterns occur when considering the evolution of the heat shock regulatory program measured in eight of the species, suggesting that these are general evolutionary principles.

DOI:http://dx.doi.org/10.7554/eLife.00603.001

The incredible diversity of living creatures belies the fact that their genes are quite similar. In the 1970s Mary-Claire King and Allan Wilson proposed that a process called gene regulation—which determines when, where and how genes are expressed as proteins—is responsible for this diversity. Four decades later, the central role of gene regulation in evolution has been confirmed in a wide range of species including bacteria, fungi, flies and mammals, although the details remain poorly understood. In recent years it has been suggested that the duplication of genes—and sometimes the duplication of whole genomes—has had a crucial influence on the part played by gene regulation in the evolution of many different species.

Ascomycota fungi are uniquely suited to the study of genetics and evolution because of their diversity—they include C. albicans, a fungus that is found in the human mouth and gut, and various species of yeast—and because many of their genomes have already been sequenced. Moreover, their genomes are relatively small, which simplifies the task of working out how it has changed over the course of evolution. It is also known that species in this branch of the tree of life diverged before and after an event in which a whole genome was duplicated.

Ascomycota fungi use glucose as a source of carbon in different ways during aerobic growth. Most, including C. albicans, are respiratory and rely on oxidative phosphorylation processes to produce energy. However, a small number—including S. cerevisiae and S. pombe, two types of yeast that are widely used as model organisms—prefer to ferment glucose, even when oxygen is available. Species that favor the latter respiro-fermentative lifestyle have evolved independently at least twice: once after the whole genome duplication event that lead to S. cerevisiae, and once when S. pombe and the other fission yeasts evolved.

Thompson et al. have measured mRNA profiles in 15 different species of yeast and reconstructed how the regulation of groups of genes (modules) have evolved over a period of more than 300 million years. They found that modules have diverged proportionally to evolutionary time, with prominent changes in gene regulation being associated with changes in lifestyle (especially changes in carbon metabolism) and a whole genome duplication event.

Gene duplication events result in gene paralogs—identical genes at different places in the genome—and these have made significant contributions to the evolution of different forms of gene regulation, especially just after the duplication event. Moreover, the paralogs produced in whole genome duplication events have resulted in bigger changes over longer periods of time. Similar patterns were observed in the regulation of the genes involved in the response to heat shock in eight of the species, which suggests that these are general evolutionary principles.

The changes in gene expression associated with the respiro-fermentative lifestyle may also have implications for our understanding of cancer: healthy cells rely on oxidative phosphorylation to produce energy whereas, similar to yeast cells, most cancerous cells rely on respiro-fermentation. Furthermore, yeast cells and cancer cells both support their rapid growth and proliferation by using glucose for biosynthesis to support cell division, although this process is not fully understood. Normal cells, on the other hand, use glucose primarily for energy and tend not to divide rapidly.

Thompson et al. found that the genes encoding enzymes in two biosynthetic pathways—one that produces the nucleotides necessary for DNA replication, and one that synthesizes glycine—are induced in respiro-fermentative yeasts but repressed in respiratory yeast cells. The fact that similar changes are observed in the same two pathways when normal cells become cancer cells suggests that these pathways have an important role in the development of cancer. The framework developed by Thompson et al. could also be used to explore the evolution of gene regulation in other species and biological processes.

DOI:http://dx.doi.org/10.7554/eLife.00603.002