Ordered and disordered regions of the Origin Recognition Complex direct differential in vivo binding at distinct motif sequences

The Origin Recognition Complex (ORC) seeds replication-fork formation by binding to DNA replication origins, which in budding yeast contain a 17bp DNA motif. High resolution structure of the ORC-DNA complex revealed two base-interacting elements: a disordered basic patch (Orc1-BP4) and an insertion helix (Orc4-IH). To define the ORC elements guiding its DNA binding in vivo, we mapped genomic locations of 38 designed ORC mutants, revealing that different ORC elements guide binding at different sites. At silencing-associated sites lacking the motif, ORC binding and activity were fully explained by a BAH domain. Within replication origins, we reveal two dominating motif variants showing differential binding modes and symmetry: a non-repetitive motif whose binding requires Orc1-BP4 and Orc4-IH, and a repetitive one where another basic patch, Orc1-BP3, can replace Orc4-IH. Disordered basic patches are therefore key for ORC-motif binding in vivo, and we discuss how these conserved, minor-groove interacting elements can guide specific ORC-DNA recognition.

Intrinsically disordered regions of the Msn2 transcription factor encode multiple functions using interwoven sequence grammars

Intrinsically disordered regions (IDRs) are abundant in eukaryotic proteins, but their sequence-function relationship remains poorly understood. IDRs of transcription factors (TFs) can direct promoter selection and recruit coactivators, as shown for the budding yeast TF Msn2. To examine how IDRs encode both these functions, we compared genomic binding specificity, coactivator recruitment, and gene induction amongst a large set of designed Msn2-IDR mutants. We find that both functions depend on multiple regions across the > 600AA IDR. Yet, transcription activity was readily disrupted by mutations that showed no effect on the Msn2 binding specificity. Our data attribute this differential sensitivity to the integration of a relaxed, composition-based code directing binding specificity with a more stringent, motif-based code controlling the recruitment of coactivators and transcription activity. Therefore, Msn2 utilizes interwoven sequence grammars for encoding multiple functions, suggesting a new IDR design paradigm of potentially general use.

The architecture of binding cooperativity between densely bound transcription factors

The binding of transcription factors (TFs) along genomes is restricted to a subset of sites containing their preferred motifs. TF-binding specificity is often attributed to the co-binding of interacting TFs; however, apart from specific examples, this model remains untested. Here, we define dependencies among budding yeast TFs that localize to overlapping promoters by profiling the genome-wide consequences of co-depleting multiple TFs. We describe unidirectional interactions, revealing Msn2 as a central factor allowing TF binding at its target promoters. By contrast, no case of mutual cooperation was observed. Particularly, Msn2 retained binding at its preferred promoters upon co-depletion of fourteen similarly bound TFs. Overall, the consequences of TF co-depletions were moderate, limited to a subset of promoters, and failed to explain the role of regions outside the DNA-binding domain in directing TF-binding preferences. Our results call for re-evaluating the role of cooperative interactions in directing TF-binding preferences.

Nucleosome retention by histone chaperones and remodelers occludes pervasive DNA-protein binding

DNA packaging within chromatin depends on histone chaperones and remodelers that form and position nucleosomes. Cells express multiple such chromatin regulators with overlapping in-vitro activities. Defining specific in-vivo activities requires monitoring histone dynamics during regulator depletion, which has been technically challenging. We have recently generated histone-exchange sensors in Saccharomyces cerevisiae, which we now use to define the contributions of 15 regulators to histone dynamics genome-wide. While replication-independent exchange in unperturbed cells maps to promoters, regulator depletions primarily affected gene bodies. Depletion of Spt6, Spt16 or Chd1 sharply increased nucleosome replacement sequentially at the beginning, middle or end of highly expressed gene bodies. They further triggered re-localization of chaperones to affected gene body regions, which compensated for nucleosome loss during transcription complex passage, but concurred with extensive TF binding in gene bodies. We provide a unified quantitative screen highlighting regulator roles in retaining nucleosome binding during transcription and preserving genomic packaging.

Histone exchange sensors reveal variant specific dynamics in mouse embryonic stem cells

Eviction of histones from nucleosomes and their exchange with newly synthesized or alternative variants is a central epigenetic determinant. Here, we define the genome-wide occupancy and exchange pattern of canonical and non-canonical histone variants in mouse embryonic stem cells by genetically encoded exchange sensors. While exchange of all measured variants scales with transcription, we describe variant-specific associations with transcription elongation and Polycomb binding. We found considerable exchange of H3.1 and H2B variants in heterochromatin and repeat elements, contrasting the occupancy and little exchange of H3.3 in these regions. This unexpected association between H3.3 occupancy and exchange of canonical variants is also evident in active promoters and enhancers, and further validated by reduced H3.1 dynamics following depletion of H3.3-specific chaperone, HIRA. Finally, analyzing transgenic mice harboring H3.1 or H3.3 sensors demonstrates the vast potential of this system for studying histone exchange and its impact on gene expression regulation in vivo.

The molecular grammar of protein disorder guiding genome-binding locations

Intrinsically disordered regions (IDRs) direct transcription factors (TFs) towards selected genomic occurrences of their binding motif, as exemplified by budding yeast's Msn2. However, the sequence basis of IDR-directed TF binding selectivity remains unknown. To reveal this sequence grammar, we analyze the genomic localizations of >100 designed IDR mutants, each carrying up to 122 mutations within this 567-AA region. Our data points at multivalent interactions, carried by hydrophobic-mostly aliphatic-residues dispersed within a disordered environment and independent of linear sequence motifs, as the key determinants of Msn2 genomic localization. The implications of our results for the mechanistic basis of IDR-based TF binding preferences are discussed.

Complementary strategies for directing in vivo transcription factor binding through DNA binding domains and intrinsically disordered regions

DNA binding domains (DBDs) of transcription factors (TFs) recognize DNA sequence motifs that are highly abundant in genomes. Within cells, TFs bind a subset of motif-containing sites as directed by either their DBDs or DBD-external (nonDBD) sequences. To define the relative roles of DBDs and nonDBDs in directing binding preferences, we compared the genome-wide binding of 48 (∼30%) budding yeast TFs with their DBD-only, nonDBD-truncated, and nonDBD-only mutants. With a few exceptions, binding locations differed between DBDs and TFs, resulting from the cumulative action of multiple determinants mapped mostly to disordered nonDBD regions. Furthermore, TFs' preferences for promoters of the fuzzy nucleosome architecture were lost in DBD-only mutants, whose binding spread across promoters, implicating nonDBDs' preferences in this hallmark of budding yeast regulatory design. We conclude that DBDs and nonDBDs employ complementary DNA-targeting strategies, whose balance defines TF binding specificity along genomes.

Evolution of transcription factor binding through sequence variations and turnover of binding sites

Variations in noncoding regulatory sequences play a central role in evolution. Interpreting such variations, however, remains difficult even in the context of defined attributes such as transcription factor (TF) binding sites. Here, we systematically link variations in cis-regulatory sequences to TF binding by profiling the allele-specific binding of 27 TFs expressed in a yeast hybrid, in which two related genomes are present within the same nucleus. TFs localize preferentially to sites containing their known consensus motifs but occupy only a small fraction of the motif-containing sites available within the genomes. Differential binding of TFs to the orthologous alleles was well explained by variations that alter motif sequence, whereas differences in chromatin accessibility between alleles were of little apparent effect. Motif variations that abolished binding when present in only one allele were still bound when present in both alleles, suggesting evolutionary compensation, with a potential role for sequence conservation at the motif's vicinity. At the level of the full promoter, we identify cases of binding-site turnover, in which binding sites are reciprocally gained and lost, yet most interspecific differences remained uncompensated. Our results show the flexibility of TFs to bind imprecise motifs and the fast evolution of TF binding sites between related species.

Rtt109 promotes nucleosome replacement ahead of the replication fork

DNA replication perturbs chromatin by triggering the eviction, replacement, and incorporation of nucleosomes. How this dynamic is orchestrated in time and space is poorly understood. Here, we apply a genetically encoded sensor for histone exchange to follow the time-resolved histone H3 exchange profile in budding yeast cells undergoing slow synchronous replication in nucleotide-limiting conditions. We find that new histones are incorporated not only behind, but also ahead of the replication fork. We provide evidence that Rtt109, the S-phase-induced acetyltransferase, stabilizes nucleosomes behind the fork but promotes H3 replacement ahead of the fork. Increased replacement ahead of the fork is independent of the primary Rtt109 acetylation target H3K56 and rather results from Vps75-dependent Rtt109 activity toward the H3 N terminus. Our results suggest that, at least under nucleotide-limiting conditions, selective incorporation of differentially modified H3s behind and ahead of the replication fork results in opposing effects on histone exchange, likely reflecting the distinct challenges for genome stability at these different regions.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge was studied in a few cases only. To provide a genome-scale view, we considered the set of budding yeast TFs classified as whole-genome duplication (WGD)-retained paralogs (~35% of all specific TFs). Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA-binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed uneven splitting of ancestral preferences between duplicates, and the preferential acquiring of new targets by the least conserved paralog (biased neo/sub-functionalization). Interactions between paralogs were rare, and, when present, occurred through weak competition for DNA-binding or dependency between dimer-forming paralogs. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

Generation and timing of graded responses to morphogen gradients

Morphogen gradients are known to subdivide a naive cell field into distinct zones of gene expression. Here, we examine whether morphogens can also induce a graded response within such domains. To this end, we explore the role of the Dorsal protein nuclear gradient along the dorsoventral axis in defining the graded pattern of actomyosin constriction that initiates gastrulation in early Drosophila embryos. Two complementary mechanisms for graded accumulation of mRNAs of crucial zygotic Dorsal target genes were identified. First, activation of target-gene expression expands over time from the ventral-most region of high nuclear Dorsal to lateral regions, where the levels are lower, as a result of a Dorsal-dependent activation probability of transcription sites. Thus, sites that are activated earlier will exhibit more mRNA accumulation. Second, once the sites are activated, the rate of RNA Polymerase II loading is also dependent on Dorsal levels. Morphological restrictions require that translation of the graded mRNA be delayed until completion of embryonic cell formation. Such timing is achieved by large introns, which provide a delay in production of the mature mRNAs. Spatio-temporal regulation of key zygotic genes therefore shapes the pattern of gastrulation.

A regulatory phosphorylation site on Mec1 controls chromatin occupancy of RNA polymerases during replication stress

Upon replication stress, budding yeast checkpoint kinase Mec1^ATR triggers the downregulation of transcription, thereby reducing the level of RNA polymerase (RNAP) on chromatin to facilitate replication fork progression. Here, we identify a hydroxyurea-induced phosphorylation site on Mec1, Mec1-S1991, that contributes to the eviction of RNAPII and RNAPIII during replication stress. The expression of the non-phosphorylatable mec1-S1991A mutant reduces replication fork progression genome-wide and compromises survival on hydroxyurea. This defect can be suppressed by destabilizing chromatin-bound RNAPII through a TAP fusion to its Rpb3 subunit, suggesting that lethality in mec1-S1991A mutants arises from replication-transcription conflicts. Coincident with a failure to repress gene expression on hydroxyurea in mec1-S1991A cells, highly transcribed genes such as GAL1 remain bound at nuclear pores. Consistently, we find that nuclear pore proteins and factors controlling RNAPII and RNAPIII are phosphorylated in a Mec1-dependent manner on hydroxyurea. Moreover, we show that Mec1 kinase also contributes to reduced RNAPII occupancy on chromatin during an unperturbed S phase by promoting degradation of the Rpb1 subunit.

Order through disorder: The role of intrinsically disordered regions in transcription factor binding specificity

Transcription factors (TFs) must bind at specific genomic locations to accurately regulate gene expression. The ability of TFs to recognize specific DNA sequence motifs arises from the inherent preferences of their globular DNA-binding domains (DBDs). Yet, these preferences are insufficient to explain the in vivo TF binding site selection. TFs are enriched with intrinsically disordered regions (IDRs), most of which are poorly characterized. While not generally considered as determinants of TF binding specificity, IDRs guide protein-protein interactions within transcriptional condensates, and multiple examples exist in which short IDRs flanking the DBD contribute to binding specificity via direct contact with the DNA. We recently reported that long IDRs, present away from the DBD, act as major specificity determinants at the genomic scale. Here, we discuss mechanisms through which IDRs contribute to DNA binding specificity, highlighting the role of long IDRs in dictating the in vivo binding site selection.

Transcription Factor Binding to Replicated DNA

Genome replication perturbs the DNA regulatory environment by displacing DNA-bound proteins, replacing nucleosomes, and introducing dosage imbalance between regions replicating at different S-phase stages. Recently, we showed that these effects are integrated to maintain transcription homeostasis: replicated genes increase in dosage, but their expression remains stable due to replication-dependent epigenetic changes that suppress transcription. Here, we examine whether reduced transcription from replicated DNA results from limited accessibility to regulatory factors by measuring the time-resolved binding of RNA polymerase II (Pol II) and specific transcription factors (TFs) to DNA during S phase in budding yeast. We show that the Pol II binding pattern is largely insensitive to DNA dosage, indicating limited binding to replicated DNA. In contrast, binding of three TFs (Reb1, Abf1, and Rap1) to DNA increases with the increasing DNA dosage. We conclude that the replication-specific chromatin environment remains accessible to regulatory factors but suppresses RNA polymerase recruitment.

Rtt109 slows replication speed by histone N-terminal acetylation

The wrapping of DNA around histone octamers challenges processes that use DNA as their template. In vitro, DNA replication through chromatin depends on histone modifiers, raising the possibility that cells modify histones to optimize fork progression. Rtt109 is an acetyl transferase that acetylates histone H3 before its DNA incorporation on the K56 and N-terminal residues. We previously reported that, in budding yeast, a wave of histone H3 K9 acetylation progresses ∼3–5 kb ahead of the replication fork. Whether this wave contributes to replication dynamics remained unknown. Here, we show that the replication fork velocity increases following deletion of RTT109, the gene encoding the enzyme required for the prereplication H3 acetylation wave. By using histone H3 mutants, we find that Rtt109-dependent N-terminal acetylation regulates fork velocity, whereas K56 acetylation contributes to replication dynamics only when N-terminal acetylation is compromised. We propose that acetylation of newly synthesized histones slows replication by promoting replacement of nucleosomes evicted by the incoming fork, thereby protecting genome integrity.

Accumulation of cis- and trans-regulatory variations is associated with phenotypic divergence of a complex trait between yeast species

Gene regulatory variations accumulate during evolution and alter gene expression. While the importance of expression variation in phenotypic evolution is well established, the molecular basis remains largely unknown. Here, we examine two closely related yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus, which show phenotypical differences in morphology and cell cycle progression when grown in the same environment. By profiling the cell cycle transcriptome and binding of key transcription factors (TFs) in the two species and their hybrid, we show that changes in expression levels and dynamics of oscillating genes are dominated by upstream trans-variations. We find that multiple cell cycle regulators show both cis- and trans-regulatory variations, which alters their expression in favor of the different cell cycle phenotypes. Moreover, we show that variations in the cell cycle TFs, Fkh1, and Fkh2 affect both the expression of target genes, and the binding specificity of an interacting TF, Ace2. Our study reveals how multiple variations accumulate and propagate through the gene regulatory network, alter TFs binding, contributing to phenotypic changes in cell cycle progression.

Speed-Specificity Trade-Offs in the Transcription Factors Search for Their Genomic Binding Sites

Transcription factors (TFs) regulate gene expression by binding DNA sequences recognized by their DNA-binding domains (DBDs). DBD-recognized motifs are short and highly abundant in genomes. The ability of TFs to bind a specific subset of motif-containing sites, and to do so rapidly upon activation, is fundamental for gene expression in all eukaryotes. Despite extensive interest, our understanding of the TF-target search process is fragmented; although binding specificity and detection speed are two facets of this same process, trade-offs between them are rarely addressed. In this opinion article, we discuss potential speed-specificity trade-offs in the context of existing models. We further discuss the recently described 'distributed specificity' paradigm, suggesting that intrinsically disordered regions (IDRs) promote specificity while reducing the TF-target search time.

Independent evolution of transcript abundance and gene regulatory dynamics

Changes in gene expression drive novel phenotypes, raising interest in how gene expression evolves. In contrast to the static genome, cells modulate gene expression in response to changing environments. Previous comparative studies focused on specific conditions, describing interspecies variation in expression levels, but providing limited information about variation across different conditions. To close this gap, we profiled mRNA levels of two related yeast species in hundreds of conditions and used coexpression analysis to distinguish variation in the dynamic pattern of gene expression from variation in expression levels. The majority of genes whose expression varied between the species maintained a conserved dynamic pattern. Cases of diverged dynamic pattern correspond to genes that were induced under distinct subsets of conditions in the two species. Profiling the interspecific hybrid allowed us to distinguish between genes with predominantly cis- or trans-regulatory variation. We find that trans-varying alleles are dominantly inherited, and that cis-variations are often complemented by variations in trans Based on these results, we suggest that gene expression diverges primarily through changes in expression levels, but does not alter the pattern by which these levels are dynamically regulated.

A Visual Framework for Classifying Determinants of Cell Size

Cells control their size by coordinating cell cycle progression with volume growth. Size control is typically studied at specific cell cycle transitions that are delayed or accelerated depending on size. This focus is well suited for revealing mechanisms acting at these transitions, but neglects the dynamics in other cell cycle phases, and is therefore inherently limited for studying how the characteristic cell size is determined. We address this limitation through a formalism that intuitively visualizes the characteristic size emerging from integrated cell cycle dynamics of individual cells. Applying this formalism to budding yeast, we describe the contributions of the un-budded (G1) and budded (S-G2-M) phase to size adjustments following environmental or genetic perturbations. We show that although the budded phase can be perturbed with little consequences for G1 dynamics, perturbations in G1 propagate to the budded phase. Our study provides an integrated view on cell size determinants in budding yeast.

•

An intuitive visualization framework for cell size control is described

•

Cell size control in different environments or mutant backgrounds can be compared

•

Mutual dependencies between size control at different cell cycle phases are described

Resolving noise-control conflict by gene duplication

Gene duplication promotes adaptive evolution in two main ways: allowing one duplicate to evolve a new function and splitting ancestral functions between the duplicates. The second scenario may resolve adaptive conflicts that can rise when one gene performs different functions. In an apparent departure from both scenarios, low-expressing transcription factor (TF) duplicates commonly bind to the same DNA motifs and act in overlapping conditions. To examine for possible benefits of this apparent redundancy, we examined the Msn2 and Msn4 duplicates in budding yeast. We show that Msn2,4 function as one unit by inducing the same set of target genes in overlapping conditions. Yet, the two-factor composition allows this unit’s expression to be both environmentally responsive and with low noise, resolving an adaptive conflict that limits expression of single genes. We propose that duplication can provide adaptive benefit through cooperation rather than functional divergence, allowing two-factor dynamics with beneficial properties that cannot be achieved by a single gene.

Gene expression tuning is limited by a conflict between expression plasticity and expression noise. This study shows that duplicated genes can resolve this conflict by functioning as one unit with expression control that is both environmentally responsive and low-noise, revealing a new adaptive benefit of duplication.

Dynamics of Spaetzle morphogen shuttling in the Drosophila embryo shapes gastrulation patterning

Establishment of morphogen gradients in the early Drosophila embryo is challenged by a diffusible sextracellular milieu, and by rapid nuclear divisions that occur at the same time. To understand how a sharp gradient is formed within this dynamic environment, we followed the generation of graded nuclear Dorsal protein, the hallmark of pattern formation along the dorso-ventral axis, in live embryos. The dynamics indicate that a sharp extracellular gradient is formed through diffusion-based shuttling of the Spaetzle (Spz) morphogen that progresses through several nuclear divisions. Perturbed shuttling in wntD mutant embryos results in a flat activation peak and aberrant gastrulation. Re-entry of Dorsal into the nuclei at the final division cycle plays an instructive role, as the residence time of Dorsal in each nucleus is translated to the amount of zygotic transcript that will be produced, thereby guiding graded accumulation of specific zygotic transcripts that drive patterned gastrulation. We conclude that diffusion-based ligand shuttling, coupled with dynamic readout, establishes a refined pattern within the diffusible environment of early embryos.

A repressor-decay timer for robust temporal patterning in embryonic Drosophila neuroblast lineages

Biological timers synchronize patterning processes during embryonic development. In the Drosophila embryo, neural progenitors (neuroblasts; NBs) produce a sequence of unique neurons whose identities depend on the sequential expression of temporal transcription factors (TTFs). The stereotypy and precision of NB lineages indicate reproducible TTF timer progression. We combine theory and experiments to define the timer mechanism. The TTF timer is commonly described as a relay of activators, but its regulatory circuit is also consistent with a repressor-decay timer, where TTF expression begins when its repressor decays. Theory shows that repressor-decay timers are more robust to parameter variations than activator-relay timers. This motivated us to experimentally compare the relative importance of the relay and decay interactions in vivo. Comparing WT and mutant NBs at high temporal resolution, we show that the TTF sequence progresses primarily by repressor-decay. We suggest that need for robust performance shapes the evolutionary-selected designs of biological circuits.

Epigenetic Control of Expression Homeostasis during Replication Is Stabilized by the Replication Checkpoint

DNA replication introduces a dosage imbalance between early and late replicating genes. In budding yeast, buffering gene expression against this imbalance depends on marking replicated DNA by H3K56 acetylation (H3K56ac). Whether additional processes are required for suppressing transcription from H3K56ac-labeled DNA remains unknown. Here, using a database-guided candidate screen, we find that COMPASS, the H3K4 methyltransferase, and its upstream effector, PAF1C, act downstream of H3K56ac to buffer expression. Replicated genes show reduced abundance of the transcription activating mark H3K4me3 and accumulate the transcription inhibitory mark H3K4me2 near transcription start sites. Notably, in hydroxyurea-exposed cells, the S phase checkpoint stabilizes H3K56ac and becomes essential for buffering. We suggest that H3K56ac suppresses transcription of replicated genes by interfering with post-replication recovery of epigenetic marks and assign a new function for the S phase checkpoint in stabilizing this mechanism during persistent dosage imbalance.

Dual role of starvation signaling in promoting growth and recovery

Growing cells are subject to cycles of nutrient depletion and repletion. A shortage of nutrients activates a starvation program that promotes growth in limiting conditions. To examine whether nutrient-deprived cells prepare also for their subsequent recovery, we followed the transcription program activated in budding yeast transferred to low-phosphate media and defined its contribution to cell growth during phosphate limitation and upon recovery. An initial transcription wave was induced by moderate phosphate depletion that did not affect cell growth. A second transcription wave followed when phosphate became growth limiting. The starvation program contributed to growth only in the second, growth-limiting phase. Notably, the early response, activated at moderate depletion, promoted recovery from starvation by increasing phosphate influx upon transfer to rich medium. Our results suggest that cells subject to nutrient depletion prepare not only for growth in the limiting conditions but also for their predicted recovery once nutrients are replenished.

Principles of cellular resource allocation revealed by condition-dependent proteome profiling

Growing cells coordinate protein translation with metabolic rates. Central to this coordination is ribosome production. Ribosomes drive cell growth, but translation of ribosomal proteins competes with production of non-ribosomal proteins. Theory shows that cell growth is maximized when all expressed ribosomes are constantly translating. To examine whether budding yeast function at this limit of full ribosomal usage, we profiled the proteomes of cells growing in different environments. We find that cells produce excess ribosomal proteins, amounting to a constant ≈8% of the proteome. Accordingly, ≈25% of ribosomal proteins expressed in rapidly growing cells does not contribute to translation. Further, this fraction increases as growth rate decreases and these excess ribosomal proteins are employed when translation demands unexpectedly increase. We suggest that steadily growing cells prepare for conditions that demand increased translation by producing excess ribosomes, at the expense of lower steady-state growth rate.

Lighting Up ERK Activity

Activation of extracellular signal regulated kinase (ERK) is used by many signaling pathways to control tissue patterning in a broad range of multicellular organisms. In this issue of Developmental Cell, Johnson et al. (2017) provide an optogenetic approach to manipulate this pathway with high precision and explore its signaling code.

A role for Rtt109 in buffering gene-dosage imbalance during DNA replication

Chromatin can function as an integrator of DNA-related processes, allowing communication, for example, between DNA replication and gene transcription. Such communication is needed to overcome the gene-dosage imbalance introduced during DNA replication, when certain genes are replicated prior to others. Increased transcription of early replicating genes could alter regulatory balances. This does not occur, suggesting a mechanism that suppresses expression from newly replicated DNA. Critical to this buffering is Rtt109, which acetylates the internal K56 residue of newly synthesized histone H3 prior to incorporation onto DNA. H3K56ac distinguishes replicated from non-replicated DNA, communicating this information to the transcription machinery to ensure expression homeostasis during S phase.

Heterosis as a consequence of regulatory incompatibility

BACKGROUND

The merging of genomes in inter-specific hybrids can result in novel phenotypes, including increased growth rate and biomass yield, a phenomenon known as heterosis. Heterosis is typically viewed as the opposite of hybrid incompatibility. In this view, the superior performance of the hybrid is attributed to heterozygote combinations that compensate for deleterious mutations accumulating in each individual genome, or lead to new, over-dominating interactions with improved performance. Still, only fragmented knowledge is available on genes and processes contributing to heterosis.

RESULTS

We describe a budding yeast hybrid that grows faster than both its parents under different environments. Phenotypically, the hybrid progresses more rapidly through cell cycle checkpoints, relieves the repression of respiration in fast growing conditions, does not slow down its growth when presented with ethanol stress, and shows increased signs of DNA damage. A systematic genetic screen identified hundreds of S. cerevisiae alleles whose deletion reduced growth of the hybrid. These growth-affecting alleles were condition-dependent, and differed greatly from alleles that reduced the growth of the S. cerevisiae parent.

CONCLUSIONS

Our results define a budding yeast hybrid that is perturbed in multiple regulatory processes but still shows a clear growth heterosis. We propose that heterosis results from incompatibilities that perturb regulatory mechanisms, which evolved to protect cells against damage or prepare them for future challenges by limiting cell growth.

Coupling phenotypic persistence to DNA damage increases genetic diversity in severe stress

Mutation rate balances the need to protect genome integrity with the advantage of evolutionary innovations. Microorganisms increase their mutation rate when stressed, perhaps addressing the growing need for evolutionary innovation. Such a strategy, however, is only beneficial under moderate stresses that allow cells to divide and realize their mutagenic potential. In contrast, severe stresses rapidly kill the majority of the population with the exception of a small minority of cells that are in a phenotypically distinct state termed persistence. Although persisters were discovered many decades ago, the stochastic event triggering persistence is poorly understood. We report that spontaneous DNA damage triggers persistence in Saccharomyces cerevisiae by activating the general stress response, providing protection against a range of harsh stress and drug environments. We further show that the persister subpopulation carries an increased load of genetic variants in the form of insertions, deletions or large structural variations, which are unrelated to their stress survival. This coupling of DNA damage to phenotypic persistence may increase genetic diversity specifically in severe stress conditions, where diversity is beneficial but the ability to generate de novo mutations is limited.

Early commitment and robust differentiation in colonic crypts

Tissue stem cells produce a constant flux of differentiated cells with distinct proportions. Here, we show that stem cells in colonic crypts differentiate early to form precisely 1:3 ratio of secretory to absorptive cells. This precision is surprising, as there are only eight stem cells making irreversible fate decisions, and so large stochastic effects of this small pool should have yielded much larger noise in cell proportions. We use single molecule FISH, lineage‐tracing mice and simulations to identify the homeostatic mechanisms facilitating robust proportions. We find that Delta‐Notch lateral inhibition operates in a restricted spatial zone to reduce initial noise in cell proportions. Increased dwell time and dispersive migration of secretory cells further averages additional variability added during progenitor divisions and breaks up continuous patches of same‐fate cells. These noise‐reducing mechanisms resolve the trade‐off between early commitment and robust differentiation and ensure spatially uniform spread of secretory cells. Our findings may apply to other cases where small progenitor pools expand to give rise to precise tissue cell proportions.

Model-based analysis of DNA replication profiles: predicting replication fork velocity and initiation rate by profiling free-cycling cells

Eukaryotic cells initiate DNA synthesis by sequential firing of hundreds of origins. This ordered replication is described by replication profiles, which measure the DNA content within a cell population. Here, we show that replication dynamics can be deduced from replication profiles of free-cycling cells. While such profiles lack explicit temporal information, they are sensitive to fork velocity and initiation capacity through the passive replication pattern, namely the replication of origins by forks emanating elsewhere. We apply our model-based approach to a compendium of profiles that include most viable budding yeast mutants implicated in replication. Predicted changes in fork velocity or initiation capacity are verified by profiling synchronously replicating cells. Notably, most mutants implicated in late (or early) origin effects are explained by global modulation of fork velocity or initiation capacity. Our approach provides a rigorous framework for analyzing DNA replication profiles of free-cycling cells.

Rethinking cell growth models

The minimal description of a growing cell consists of self-replicating ribosomes translating the cellular proteome. While neglecting all other cellular components, this model provides key insights into the control and limitations of growth rate. It shows, for example, that growth rate is maximized when ribosomes work at full capacity, explains the linear relation between growth rate and the ribosome fraction of the proteome and defines the maximal possible growth rate. This ribosome-centered model also highlights the challenge of coordinating cell growth with related processes such as cell division or nutrient production. Coordination is promoted when ribosomes don't translate at maximal capacity, as it allows escaping strict exponential growth. Recent data support the notion that multiple cellular processes limit growth. In particular, increasing transcriptional demand may be as deleterious as increasing translational demand, depending on growth conditions. Consistent with the idea of trade-off, cells may forgo maximal growth to enable more efficient interprocess coordination and faster adaptation to changing conditions.

Dealing with Gene-Dosage Imbalance during S Phase

DNA replication perturbs the dosage balance between genes that replicate early during S phase and those that replicate late. If propagated to influence protein content, this dosage imbalance could influence cellular functions. In bacteria, mechanisms have evolved to use this imbalance to tune certain processes with the rate of cell growth. By contrast, eukaryotes buffer this dosage imbalance to ensure gene expression homeostasis also during S phase. Here, we outline classical and more recent studies describing how different organisms deal with this replication-dependent dosage imbalance, and describe recent results linking the eukaryotic buffering mechanism to replication-dependent histone acetylation. Finally, we discuss the possible implications of this buffering mechanism and speculate why it is specific to eukaryote cells.

Rapid evolutionary adaptation to growth on an 'unfamiliar' carbon source

Background

Cells constantly adapt to changes in their environment. When environment shifts between conditions that were previously encountered during the course of evolution, evolutionary-programmed responses are possible. Cells, however, may also encounter a new environment to which a novel response is required. To characterize the first steps in adaptation to a novel condition, we studied budding yeast growth on xylulose, a sugar that is very rarely found in the wild.

Results

We previously reported that growth on xylulose induces the expression of amino acid biosynthesis genes in multiple natural yeast isolates. This induction occurs despite the presence of amino acids in the growth medium and is a unique response to xylulose, not triggered by naturally available carbon sources. Propagating these strains for ~300 generations on xylulose significantly improved their growth rate. Notably, the most significant change in gene expression was the loss of amino acid biosynthesis gene induction. Furthermore, the reduction in amino-acid biosynthesis gene expression on xylulose was tightly correlated with the improvement in growth rate, suggesting that internal depletion of amino-acids presented a major bottleneck limiting growth in xylulose.

Conclusions

We discuss the possible implications of our results for explaining how cells maintain the balance between supply and demand of amino acids during growth in evolutionary ‘familiar’ vs. ‘novel’ conditions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3010-x) contains supplementary material, which is available to authorized users.

A WntD-Dependent Integral Feedback Loop Attenuates Variability in Drosophila Toll Signaling

Patterning by morphogen gradients relies on the capacity to generate reproducible distribution profiles. Morphogen spread depends on kinetic parameters, including diffusion and degradation rates, which vary between embryos, raising the question of how variability is controlled. We examined this in the context of Toll-dependent dorsoventral (DV) patterning of the Drosophila embryo. We find that low embryo-to-embryo variability in DV patterning relies on wntD, a Toll-target gene expressed initially at the posterior pole. WntD protein is secreted and disperses in the extracellular milieu, associates with its receptor Frizzled4, and inhibits the Toll pathway by blocking the Toll extracellular domain. Mathematical modeling predicts that WntD accumulates until the Toll gradient narrows to its desired spread, and we support this feedback experimentally. This circuit exemplifies a broadly applicable induction-contraction mechanism, which reduces patterning variability through a restricted morphogen-dependent expression of a secreted diffusible inhibitor.

Chromatin dynamics during DNA replication

Chromatin is composed of DNA and histones, which provide a unified platform for regulating DNA-related processes, mostly through their post-translational modification. During DNA replication, histone arrangement is perturbed, first to allow progression of DNA polymerase and then during repackaging of the replicated DNA. To study how DNA replication influences the pattern of histone modification, we followed the cell-cycle dynamics of 10 histone marks in budding yeast. We find that histones deposited on newly replicated DNA are modified at different rates: While some marks appear immediately upon replication (e.g., H4K16ac, H3K4me1), others increase with transcription-dependent delays (e.g., H3K4me3, H3K36me3). Notably, H3K9ac was deposited as a wave preceding the replication fork by ∼5–6 kb. This replication-guided H3K9ac was fully dependent on the acetyltransferase Rtt109, while expression-guided H3K9ac was deposited by Gcn5. Further, topoisomerase depletion intensified H3K9ac in front of the replication fork and in sites where RNA polymerase II was trapped, suggesting supercoiling stresses trigger H3K9 acetylation. Our results assign complementary roles for DNA replication and gene expression in defining the pattern of histone modification.

Expression homeostasis during DNA replication

Genome replication introduces a stepwise increase in the DNA template available for transcription. Genes replicated early in S phase experience this increase before late-replicating genes, raising the question of how expression levels are affected by DNA replication. We show that in budding yeast, messenger RNA (mRNA) synthesis rate is buffered against changes in gene dosage during S phase. This expression homeostasis depends on acetylation of H3 on its internal K56 site by Rtt109/Asf1. Deleting these factors, mutating H3K56 or up-regulating its deacetylation, increases gene expression in S phase in proportion to gene replication timing. Therefore, H3K56 acetylation on newly deposited histones reduces transcription efficiency from replicated DNA, complementing its role in guarding genome stability. Our study provides molecular insight into the mechanism maintaining expression homeostasis during DNA replication.

A two-step patterning process increases the robustness of periodic patterning in the fly eye

Complex periodic patterns can self-organize through dynamic interactions between diffusible activators and inhibitors. In the biological context, self-organized patterning is challenged by spatial heterogeneities ('noise') inherent to biological systems. How spatial variability impacts the periodic patterning mechanism and how it can be buffered to ensure precise patterning is not well understood. We examine the effect of spatial heterogeneity on the periodic patterning of the fruit fly eye, an organ composed of ∼800 miniature eye units (ommatidia) whose periodic arrangement along a hexagonal lattice self-organizes during early stages of fly development. The patterning follows a two-step process, with an initial formation of evenly spaced clusters of ∼10 cells followed by a subsequent refinement of each cluster into a single selected cell. Using a probabilistic approach, we calculate the rate of patterning errors resulting from spatial heterogeneities in cell size, position and biosynthetic capacity. Notably, error rates were largely independent of the desired cluster size but followed the distributions of signaling speeds. Pre-formation of large clusters therefore greatly increases the reproducibility of the overall periodic arrangement, suggesting that the two-stage patterning process functions to guard the pattern against errors caused by spatial heterogeneities. Our results emphasize the constraints imposed on self-organized patterning mechanisms by the need to buffer stochastic effects. Author summary Complex periodic patterns are common in nature and are observed in physical, chemical and biological systems. Understanding how these patterns are generated in a precise manner is a key challenge. Biological patterns are especially intriguing, as they are generated in a noisy environment; cell position and cell size, for example, are subject to stochastic variations, as are the strengths of the chemical signals mediating cell-to-cell communication. The need to generate a precise and robust pattern in this 'noisy' environment restricts the space of patterning mechanisms that can function in the biological setting. Mathematical modeling is useful in comparing the sensitivity of different mechanisms to such variations, thereby highlighting key aspects of their design.We use mathematical modeling to study the periodic patterning of the fruit fly eye. In this system, a highly ordered lattice of differentiated cells is generated in a two-dimensional cell epithelium. The pattern is first observed by the appearance of evenly spaced clusters of ∼10 cells that express specific genes. Each cluster is subsequently refined into a single cell, which initiates the formation and differentiation of a miniature eye unit, the ommatidium. We formulate a mathematical model based on the known molecular properties of the patterning mechanism, and use a probabilistic approach to calculate the errors in cluster formation and refinement resulting from stochastic cell-to-cell variations ('noise') in different quantitative parameters. This enables us to define the parameters most influencing noise sensitivity. Notably, we find that this error is roughly independent of the desired cluster size, suggesting that large clusters are beneficial for ensuring the overall reproducibility of the periodic cluster arrangement. For the stage of cluster refinement, we find that rapid communication between cells is critical for reducing error. Our work provides new insights into the constraints imposed on mechanisms generating periodic patterning in a realistic, noisy environment, and in particular, discusses the different considerations in achieving optimal design of the patterning network.

Periodic patterning of the Drosophila eye is stabilized by the diffusible activator Scabrous

Generation of periodic patterns is fundamental to the differentiation of multiple tissues during development. How such patterns form robustly is still unclear. The Drosophila eye comprises ∼750 units, whose crystalline order is set during differentiation of the eye imaginal disc: an activation wave sweeping across the disc is coupled to lateral inhibition, sequentially selecting pro-neural cells. Using mathematical modelling, here we show that this template-based lateral inhibition is highly sensitive to spatial variations in biochemical parameters and cell sizes. We reveal the basis of this sensitivity, and suggest that it can be overcome by assuming a short-range diffusible activator. Clonal experiments identify Scabrous, a previously implicated inhibitor, as the predicted activator. Our results reveal the mechanism by which periodic patterning in the fly eye is stabilized against spatial variations, highlighting how the need to maintain robustness shapes the design of patterning circuits.

Patterning in the Drosophila eye is achieved by a series of signalling cascades over several cell distances. Here Gavish et al. model lateral inhibition in the developing eye to understand how developmental noise refines such patterning, identifying a novel activator required for buffering spatial variability.

DNA Crossover Motifs Associated with Epigenetic Modifications Delineate Open Chromatin Regions in Arabidopsis

The rate of crossover, the reciprocal exchanges of homologous chromosomal segments, is not uniform along chromosomes differing between male and female meiocytes. To better understand the factors regulating this variable landscape, we performed a detailed genetic and epigenetic analysis of 737 crossover events in Arabidopsis thaliana. Crossovers were more frequent than expected in promoters. Three DNA motifs enriched in crossover regions and less abundant in crossover-poor pericentric regions were identified. One of these motifs, the CCN repeat, was previously unknown in plants. The A-rich motif was preferentially associated with promoters, while the CCN repeat and the CTT repeat motifs were preferentially associated with genes. Analysis of epigenetic modifications around the motifs showed, in most cases, a specific epigenetic architecture. For example, we show that there is a peak of nucleosome occupancy and of H3K4me3 around the CCN and CTT repeat motifs while nucleosome occupancy was lowest around the A-rich motif. Cytosine methylation levels showed a gradual decrease within ∼2 kb of the three motifs, being lowest at sites where crossover occurred. This landscape was conserved in the decreased DNA methylation1 mutant. In summary, the crossover motifs are associated with epigenetic landscapes corresponding to open chromatin and contributing to the nonuniformity of crossovers in Arabidopsis.

The Evolutionary Potential of Phenotypic Mutations

Errors in protein synthesis, so-called phenotypic mutations, are orders-of-magnitude more frequent than genetic mutations. Here, we provide direct evidence that alternative protein forms and phenotypic variability derived from translational errors paved the path to genetic, evolutionary adaptations via gene duplication. We explored the evolutionary origins of Saccharomyces cerevisiae IDP3 - an NADP-dependent isocitrate dehydrogenase mediating fatty acids ß-oxidation in the peroxisome. Following the yeast whole genome duplication, IDP3 diverged from a cytosolic ancestral gene by acquisition of a C-terminal peroxisomal targeting signal. We discovered that the pre-duplicated cytosolic IDPs are partially localized to the peroxisome owing to +1 translational frameshifts that bypass the stop codon and unveil cryptic peroxisomal targeting signals within the 3’-UTR. Exploring putative cryptic signals in all 3’-UTRs of yeast genomes, we found that other enzymes related to NADPH production such as pyruvate carboxylase 1 (PYC1) might be prone to peroxisomal localization via cryptic signals. Using laboratory evolution we found that these translational frameshifts are rapidly imprinted via genetic single base deletions occurring within the very same gene location. Further, as exemplified here, the sequences that promote translational frameshifts are also more prone to genetic deletions. Thus, genotypes conferring higher phenotypic variability not only meet immediate challenges by unveiling cryptic 3’-UTR sequences, but also boost the potential for future genetic adaptations.

The rarity of genetic mutations limits the likelihood of adaptation. However, transcriptional and translational errors, so-called phenotypic mutations, are >10⁵-fold more frequent, thus generating protein mutants from unmodified genes. We provide the first evidence that phenotypic mutations paved the path to what later, after gene duplication, became newly compartmentalized enzymes. Thus, gene duplication followed rather than initiated the divergence of this new trait. Our findings also show that translational infidelity and phenotypic variability comprise the origins of evolutionary innovations, and how selection for enhanced phenotypic variability also promotes the appearance of genetic mutations that lead to the very same outcome.

Loss of growth homeostasis by genetic decoupling of cell division from biomass growth: implication for size control mechanisms

Growing cells adjust their division time with biomass accumulation to maintain growth homeostasis. Size control mechanisms, such as the size checkpoint, provide an inherent coupling of growth and division by gating certain cell cycle transitions based on cell size. We describe genetic manipulations that decouple cell division from cell size, leading to the loss of growth homeostasis, with cells becoming progressively smaller or progressively larger until arresting. This was achieved by modulating glucose influx independently of external glucose. Division rate followed glucose influx, while volume growth was largely defined by external glucose. Therefore, the coordination of size and division observed in wild-type cells reflects tuning of two parallel processes, which is only refined by an inherent feedback-dependent coupling. We present a class of size control models explaining the observed breakdowns of growth homeostasis.

Systematic identification of cell size regulators in budding yeast

Cell size is determined by a complex interplay between growth and division, involving multiple cellular pathways. To identify systematically processes affecting size control in G1 in budding yeast, we imaged and analyzed the cell cycle of millions of individual cells representing 591 mutants implicated in size control. Quantitative metric distinguished mutants affecting the mechanism of size control from the majority of mutants that have a perturbed size due to indirect effects modulating cell growth. Overall, we identified 17 negative and dozens positive size control regulators, with the negative regulators forming a small network centered on elements of mitotic exit network. Some elements of the translation machinery affected size control with a notable distinction between the deletions of parts of small and large ribosomal subunit: parts of small ribosomal subunit tended to regulate size control, while parts of the large subunit affected cell growth. Analysis of small cells revealed additional size control mechanism that functions in G2/M, complementing the primary size control in G1. Our study provides new insights about size control mechanisms in budding yeast.