Evolution of nucleosome occupancy: conservation of global properties and divergence of gene-specific patterns

Prediction of G4 formation in live cells with epigenetic data: a deep learning approach

G-quadruplexes (G4s) are secondary structures abundant in DNA that may play regulatory roles in cells. Despite the ubiquity of the putative G-quadruplex-forming sequences (PQS) in the human genome, only a small fraction forms G4 structures in cells. Folded G4, histone methylation and chromatin accessibility are all parts of the complex cis regulatory landscape. We propose an approach for prediction of G4 formation in cells that incorporates epigenetic and chromatin accessibility data. The novel approach termed epiG4NN efficiently predicts cell-specific G4 formation in live cells based on a local epigenomic snapshot. Our results confirm the close relationship between H3K4me3 histone methylation, chromatin accessibility and G4 structure formation. Trained on A549 cell data, epiG4NN was then able to predict G4 formation in HEK293T and K562 cell lines. We observe the dependency of model performance with different epigenetic features on the underlying experimental condition of G4 detection. We expect that this approach will contribute to the systematic understanding of correlations between structural and epigenomic feature landscape.

A strong structural correlation between short inverted repeat sequences and the polyadenylation signal in yeast and nucleosome exclusion by these inverted repeats

DNA sequences that read the same from 5′ to 3′ in either strand are called inverted repeat sequences or simply IRs. They are found throughout a wide variety of genomes, from prokaryotes to eukaryotes. Despite extensive research, their in vivo functions, if any, remain unclear. Using Saccharomyces cerevisiae, we performed genome-wide analyses for the distribution, occurrence frequency, sequence characteristics and relevance to chromatin structure, for the IRs that reportedly have a cruciform-forming potential. Here, we provide the first comprehensive map of these IRs in the S. cerevisiae genome. The statistically significant enrichment of the IRs was found in the close vicinity of the DNA positions corresponding to polyadenylation [poly(A)] sites and ~ 30 to ~ 60 bp downstream of start codon-coding sites (referred to as ‘start codons’). In the former, ApT- or TpA-rich IRs and A-tract- or T-tract-rich IRs are enriched, while in the latter, different IRs are enriched. Furthermore, we found a strong structural correlation between the former IRs and the poly(A) signal. In the chromatin formed on the gene end regions, the majority of the IRs causes low nucleosome occupancy. The IRs in the region ~ 30 to ~ 60 bp downstream of start codons are located in the + 1 nucleosomes. In contrast, fewer IRs are present in the adjacent region downstream of start codons. The current study suggests that the IRs play similar roles in Escherichia coli and S. cerevisiae to regulate or complete transcription at the RNA level.

DNA methylation and chromatin accessibility profiling of mouse and human fetal germ cells

Chromatin remodeling is important for the epigenetic reprogramming of human primordial germ cells. However, the comprehensive chromatin state has not yet been analyzed for human fetal germ cells (FGCs). Here we use nucleosome occupancy and methylation sequencing method to analyze both the genome-wide chromatin accessibility and DNA methylome at a series of crucial time points during fetal germ cell development in both human and mouse. We find 116 887 and 137 557 nucleosome-depleted regions (NDRs) in human and mouse FGCs, covering a large set of germline-specific and highly dynamic regulatory genomic elements, such as enhancers. Moreover, we find that the distal NDRs are enriched specifically for binding motifs of the pluripotency and germ cell master regulators such as NANOG, SOX17, AP2γ and OCT4 in human FGCs, indicating the existence of a delicate regulatory balance between pluripotency-related genes and germ cell-specific genes in human FGCs, and the functional significance of these genes for germ cell development in vivo. Our work offers a comprehensive and high-resolution roadmap for dissecting chromatin state transition dynamics during the epigenomic reprogramming of human and mouse FGCs.

Genome-wide mapping of nucleosome positions in Saccharomyces cerevisiae in response to different nitrogen conditions

Well-organized chromatin is involved in a number of various transcriptional regulation and gene expression. We used genome-wide mapping of nucleosomes in response to different nitrogen conditions to determine both nucleosome profiles and gene expression events in Saccharomyces cerevisiae. Nitrogen conditions influence general nucleosome profiles and the expression of nitrogen catabolite repression (NCR) sensitive genes. The nucleosome occupancy of TATA-containing genes was higher compared to TATA-less genes. TATA-less genes in high or low nucleosome occupancy, showed a significant change in gene coding regions when shifting cells from glutamine to proline as the sole nitrogen resource. Furthermore, a correlation between the expression of nucleosome occupancy induced NCR sensitive genes or TATA containing genes in NCR sensitive genes, and nucleosome prediction were found when cells were cultured in proline or shifting from glutamine to proline as the sole nitrogen source compared to glutamine. These results also showed that variation of nucleosome occupancy accompany with chromatin-dependent transcription factor could influence the expression of a series of genes involved in the specific regulation of nitrogen utilization.

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.

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.

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.

Chromatin mediation of a transcriptional memory effect in yeast

Previous studies have described a transcriptional "memory effect," whereby transcript levels of many Abf1-regulated genes in the budding yeast Saccharomyces cerevisiae are undiminished even after Abf1 has dissociated from its regulatory sites. Here we provide additional support for this effect and investigate its molecular basis. We show that the effect is observed in a distinct abf1 ts mutant from that used in earlier studies, demonstrating that it is robust, and use chromatin immunoprecipitation to show that Abf1 association is decreased similarly from memory effect and transcriptionally responsive genes at the restrictive temperature. We also demonstrate that the association of TATA-binding protein and Pol II decreases after the loss of Abf1 binding for transcriptionally responsive genes but not for memory effect genes. Examination of genome-wide nucleosome occupancy data reveals that although transcriptionally responsive genes exhibit increased nucleosome occupancy in abf1 ts yeast, the promoter regions of memory effect targets show no change in abf1 ts mutants, maintaining an open chromatin conformation even after Abf1 eviction. This contrasting behavior reflects different inherent propensity for nucleosome formation between the two classes, driven by the presence of A/T-rich sequences upstream of the Abf1 site in memory effect gene promoters. These sequence-based differences show conservation in closely related fungi and also correlate with different gene expression noise, suggesting a physiological basis for greater access to "memory effect" promoter regions. Thus, our results establish a conserved mechanism underlying a transcriptional memory effect whereby sequences surrounding Abf1 binding sequences affect local nucleosome occupancy following loss of Abf1 binding. Furthermore, these findings demonstrate that sequence-based differences in the propensity for nucleosome occupancy can influence the transcriptional response of genes to an altered regulatory signal.

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.

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.

Genome-wide survey of transcriptional initiation in the pathogenic fungus, Candida glabrata

DNA sequencing of the 5'-flanking region of the transcriptome effectively identifies transcription initiation sites and also aids in identifying unknown genes. This study describes a comprehensive polling of transcription start sites and an analysis of full-length complementary DNAs derived from the genome of the pathogenic fungus Candida glabrata. A comparison of the sequence reads derived from a cDNA library prepared from cells grown under different culture conditions against the reference genomic sequence of the Candida Genome Database (CGD: http://www.candidagenome.org/) revealed the expression of 4316 genes and their acknowledged transcription start sites (TSSs). In addition this analysis also predicted 59 new genes including 22 that showed no homology to the genome of Saccharomyces cerevisiae, a genetically close relative of C. glabrata. Furthermore, comparison of the 5'-untranslated regions (5'-UTRs) and core promoters of C. glabrata to those of S. cerevisiae showed various global similarities and differences among orthologous genes. Thus, the C. glabrata transcriptome can complement the annotation of the genome database and should provide new insights into the organization, regulation, and function of genes of this important human pathogen.

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.

Conserved substitution patterns around nucleosome footprints in eukaryotes and Archaea derive from frequent nucleosome repositioning through evolution

Nucleosomes, the basic repeat units of eukaryotic chromatin, have been suggested to influence the evolution of eukaryotic genomes, both by altering the propensity of DNA to mutate and by selection acting to maintain or exclude nucleosomes in particular locations. Contrary to the popular idea that nucleosomes are unique to eukaryotes, histone proteins have also been discovered in some archaeal genomes. Archaeal nucleosomes, however, are quite unlike their eukaryotic counterparts in many respects, including their assembly into tetramers (rather than octamers) from histone proteins that lack N- and C-terminal tails. Here, we show that despite these fundamental differences the association between nucleosome footprints and sequence evolution is strikingly conserved between humans and the model archaeon Haloferax volcanii. In light of this finding we examine whether selection or mutation can explain concordant substitution patterns in the two kingdoms. Unexpectedly, we find that neither the mutation nor the selection model are sufficient to explain the observed association between nucleosomes and sequence divergence. Instead, we demonstrate that nucleosome-associated substitution patterns are more consistent with a third model where sequence divergence results in frequent repositioning of nucleosomes during evolution. Indeed, we show that nucleosome repositioning is both necessary and largely sufficient to explain the association between current nucleosome positions and biased substitution patterns. This finding highlights the importance of considering the direction of causality between genetic and epigenetic change.

Genome sequences as well as epigenetic states, such as DNA methylation or nucleosome binding patterns, change during evolution. But what is the causal relationship between the two? We already know that nucleotide variation within and between species is distributed unevenly around nucleosome footprints, but does this mean that sequence evolution follows a biased course because the presence of nucleosomes affects mutation and DNA repair dynamics? Or is it, in fact, the other way around, i.e. changes happen at the DNA level and prompt shifts in nucleosome positioning? To investigate the direction of causality in genetic versus epigenetic evolution, we analyze substitutions patterns in eukaryotes as well as the archaeon Haloferax volcanii in the context of genome-wide nucleosome binding maps. We demonstrate that the relationship between nucleosome positions and between-species divergence patterns, strikingly similar in eukaryotes and archaea, can be explained in large parts by nucleosomes shifting positions in response to substitution, although both mutation and selection biases might still exist. Our results illustrate that it is important to consider the direction of causality between epigenetic and genetic change when analyzing patterns of sequence divergence and using sequence conservation to infer selection on epigenetic states.

NCS: incorporating positioning data to quantify nucleosome stability in yeast

MOTIVATION

With the spreading technique of mass sequencing, nucleosome positions and scores for their intensity have become available through several previous studies in yeast, but relatively few studies have specifically aimed to determine the score of nucleosome stability. Based on mass sequencing data, we proposed a nucleosome center score (NCS) for quantifying nucleosome stability by measuring shifts of the nucleosome center, and then mapping NCS scores to nucleosome positions in Brogaard et al.'s study.

RESULTS

We demonstrated the efficiency of NCS by known preference of A/T-based tracts for nucleosome formation, and showed that central nucleosomal DNA is more sensitive to A/T-based tracts than outer regions, which corresponds to the central histone tetramer-dominated region. We also found significant flanking preference around nucleosomal DNA for A/T-based dinucleotides, suggesting that neighboring sequences could affect nucleosome stability. Finally, the difference between results of NCS and Brogaard et al.'s scores was addressed and discussed.

CONTACTS

jchiang@mail.ncku.edu.tw

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Abundant gene-by-environment interactions in gene expression reaction norms to copper within Saccharomyces cerevisiae

Genetic variation for plastic phenotypes potentially contributes phenotypic variation to populations that can be selected during adaptation to novel ecological contexts. However, the basis and extent of plastic variation that manifests in diverse environments remains elusive. Here, we characterize copper reaction norms for mRNA abundance among five Saccharomyces cerevisiae strains to 1) describe population variation across the full range of ecologically relevant copper concentrations, from starvation to toxicity, and 2) to test the hypothesis that plastic networks exhibit increased population variation for gene expression. We find that although the vast majority of the variation is small in magnitude (considerably <2-fold), not just some, but most genes demonstrate variable expression across environments, across genetic backgrounds, or both. Plastically expressed genes included both genes regulated directly by copper-binding transcription factors Mac1 and Ace1 and genes indirectly responding to the downstream metabolic consequences of the copper gradient, particularly genes involved in copper, iron, and sulfur homeostasis. Copper-regulated gene networks exhibited more similar behavior within the population in environments where those networks have a large impact on fitness. Nevertheless, expression variation in genes like Cup1, important to surviving copper stress, was linked with variation in mitotic fitness and in the breadth of differential expression across the genome. By revealing a broader and deeper range of population variation, our results provide further evidence for the interconnectedness of genome-wide mRNA levels, their dependence on environmental context and genetic background, and the abundance of variation in gene expression that can contribute to future evolution.

Nucleosome Positioning

Nucleosome positioning is not only related to genomic DNA compaction but also to other biological functions. After the chromatin is digested by micrococcal nuclease, nucleosomal (nucleosome-bound) DNA fragments can be sequenced and mapped on the genomic DNA sequence. Due to the development of modern DNA sequencing technology, genome-wide nucleosome mapping has been performed in a wide range of eukaryotic species. Comparative analyses of the nucleosome positions have revealed that the nucleosome is more frequently formed in exonic than intronic regions, and that most of transcription start and translation (or transcription) end sites are located in nucleosome linker DNA regions, indicating that nucleosome positioning influences transcription initiation, transcription termination, and gene splicing. In addition, nucleosomal DNA contains guanine and cytosine (G + C)-rich sequences and a high level of cytosine methylation. Thus, the nucleosome positioning system has been conserved during eukaryotic evolution.

Noise-mean relationship in mutated promoters

Gene expression depends on the frequency of transcription events (burst frequency) and on the number of mRNA molecules made per event (burst size). Both processes are encoded in promoter sequence, yet their dependence on mutations is poorly understood. Theory suggests that burst size and frequency can be distinguished by monitoring the stochastic variation (noise) in gene expression: Increasing burst size will increase mean expression without changing noise, while increasing burst frequency will increase mean expression and decrease noise. To reveal principles by which promoter sequence regulates burst size and frequency, we randomly mutated 22 yeast promoters chosen to span a range of expression and noise levels, generating libraries of hundreds of sequence variants. In each library, mean expression (m) and noise (coefficient of variation, η) varied together, defining a scaling curve: η(2) = b/m + η(ext)(2). This relation is expected if sequence mutations modulate burst frequency primarily. The estimated burst size (b) differed between promoters, being higher in promoter containing a TATA box and lacking a nucleosome-free region. The rare variants that significantly decreased b were explained by mutations in TATA, or by an insertion of an out-of-frame translation start site. The decrease in burst size due to mutations in TATA was promoter-dependent, but independent of other mutations. These TATA box mutations also modulated the responsiveness of gene expression to changing conditions. Our results suggest that burst size is a promoter-specific property that is relatively robust to sequence mutations but is strongly dependent on the interaction between the TATA box and promoter nucleosomes.

Comparative studies of gene expression and the evolution of gene regulation

The hypothesis that differences in gene regulation have an important role in speciation and adaptation is more than 40 years old. With the advent of new sequencing technologies, we are able to characterize and study gene expression levels and associated regulatory mechanisms in a large number of individuals and species at an unprecedented resolution and scale. We have thus gained new insights into the evolutionary pressures that shape gene expression levels and have developed an appreciation for the relative importance of evolutionary changes in different regulatory genetic and epigenetic mechanisms. The current challenge is to link gene regulatory changes to adaptive evolution of complex phenotypes. Here we mainly focus on comparative studies in primates and how they are complemented by studies in model organisms.

Comparative analyses of homocitrate synthase genes of ascomycetous yeasts

Most ascomycetous yeasts have 2 homocitrate synthases (HCSs). Among the fungal lysine biosynthesis-related genes, only the HCS gene was duplicated in the course of evolution. It was recently reported that HCS of Saccharomyces cerevisiae has an additional function in nuclear activities involving chromatin regulation related to DNA damage repair, which is not related to lysine biosynthesis. Thus, it is possible that the bifunctionality is associated with HCS gene duplication. Phylogenetic analysis showed that duplication has occurred multiple times during evolution of the ascomycetous yeasts. It is likely that the HCS gene duplication in S. cerevisiae occurred in the course of Saccharomyces evolution. Although the nucleosome position profiles of the two S. cerevisiae HCS genes were similar in the coding regions, they were different in the promoter regions, suggesting that they are subject to different regulatory controls. S. cerevisiae has maintained HCS activity for lysine biosynthesis and has obtained bifunctionality.

Nucleosome organization affects the sensitivity of gene expression to promoter mutations

Gene expression diverges rapidly between related species, playing a key role in the evolution of new phenotypes. The extent of divergence differs greatly between genes and is correlated to promoter nucleosome organization. We hypothesized that this may be partially explained by differential sensitivity of expression to mutations in the promoter region. We measured the sensitivity of 22 yeast promoters with varying nucleosome patterns to random mutations in sequence. Mutation sensitivity differed by up to 10-fold between promoters. This difference could not be explained by the abundance of transcription factor binding sites. Rather, mutation sensitivity positively correlated with the relative occupancy of nucleosomes at the proximal promoter region. Furthermore, mutation sensitivity was reduced upon introduction of a binding site for Reb1, a factor that blocks nucleosome formation, suggesting that nucleosome organization directly regulates mutation sensitivity. Our study suggests an important role for chromatin structure in the evolution of gene expression.

Promoter nucleosome organization shapes the evolution of gene expression

Understanding why genes evolve at different rates is fundamental to evolutionary thinking. In species of the budding yeast, the rate at which genes diverge in expression correlates with the organization of their promoter nucleosomes: genes lacking a nucleosome-free region (denoted OPN for “Occupied Proximal Nucleosomes”) vary widely between the species, while the expression of those containing NFR (denoted DPN for “Depleted Proximal Nucleosomes”) remains largely conserved. To examine if early evolutionary dynamics contributes to this difference in divergence, we artificially selected for high expression of GFP–fused proteins. Surprisingly, selection was equally successful for OPN and DPN genes, with ∼80% of genes in each group stably increasing in expression by a similar amount. Notably, the two groups adapted by distinct mechanisms: DPN–selected strains duplicated large genomic regions, while OPN–selected strains favored trans mutations not involving duplications. When selection was removed, DPN (but not OPN) genes reverted rapidly to wild-type expression levels, consistent with their lower diversity between species. Our results suggest that promoter organization constrains the early evolutionary dynamics and in this way biases the path of long-term evolution.

Species diverge by mutations that change protein structure or protein regulation. While the evolution of protein sequence was studied extensively, much less is known about the divergence of gene expression. To better understand the process of gene expression evolution, we characterized the early genomic response of yeast cells to selection for high gene expression. Notably, the response to selection was strongly dependent on the organization of the gene promoter: genes whose promoters had a pronounced nucleosome-free region (NFR) primarily duplicated the chromosome containing the gene of interest, while genes whose promoters lacked a pronounced NFR adapted by trans mutations not involving duplications. Further, when selection was removed, the former (but not the later) evolved strains reverted rapidly to wild-type expression levels, consistent with their lower diversity between species. Together, our study provides strong support to the idea that physiological regulation impacts the evolutionary path and suggests that, by regulating promoter nucleosomes, cells can regulate the response to selection and control the long-term stability of the selected changes.

Nucleosome positioning: bringing order to the eukaryotic genome

Nucleosomes are an essential component of eukaryotic chromosomes. The impact of nucleosomes is seen not just on processes that directly access the genome, such as transcription, but also on an evolutionary timescale. Recent studies in various organisms have provided high-resolution maps of nucleosomes throughout the genome. Computational analysis, in conjunction with many other kinds of data, has shed light on several aspects of nucleosome biology. Nucleosomes are positioned by several means, including intrinsic sequence biases, by stacking against a fixed barrier, by DNA-binding proteins and by chromatin remodelers. These studies underscore the important organizational role of nucleosomes in all eukaryotic genomes. This paper reviews recent genomic studies that have shed light on the determinants of nucleosome positioning and their impact on the genome.