Search

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.

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.

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.

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.

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.

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.

During evolution, mutations produce new lineages that gradually diverge in sequence and regulatory properties. Related strains or species can hybridize to produce viable offspring. Hybrids often outperform their parents, producing more biomass or growing more rapidly. This superior performance, termed heterosis, contrasts the more expected clash between the genomes, and has puzzled geneticists and evolutionary biologists for many years. In this review, we describe two classes of models explaining heterosis: the prevailing view attributes heterosis to rapid repair or enhancement of growth promoting pathways. An alternative view attributes heterosis to the impairment of growth-limiting pathways. The two classes are not mutually exclusive and can result from similar types of genetic interactions. We discuss the possible implications of heterosis on tradeoffs in species evolution.

Cells grow on a wide range of carbon sources by regulating substrate flow through the metabolic network. Incoming sugar, for example, can be fermented or respired, depending on the carbon identity, cell type, or growth conditions. Despite this genetically-encoded flexibility of carbon metabolism, attempts to exogenously manipulate central carbon flux by rational design have proven difficult, suggesting a robust network structure. To examine this robustness, we characterized the ethanol yield of 411 regulatory and metabolic mutants in budding yeast. The mutants showed little variation in ethanol productivity when grown on glucose or galactose, yet diversity was revealed during growth on xylulose, a rare pentose not widely available in nature. While producing ethanol at high yield, cells grown on xylulose produced ethanol at high yields, yet induced expression of respiratory genes, and were dependent on them. Analysis of mutants that affected ethanol productivity suggested that xylulose fermentation results from metabolic overflow, whereby the flux through glycolysis is higher than the maximal flux that can enter respiration. We suggest that this overflow results from a suboptimal regulatory adjustment of the cells to this unfamiliar carbon source.

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. Electronic supplementary material The online version of this article (doi:10.1186/s12915-017-0373-7) contains supplementary material, which is available to authorized users.

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.

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.

Cell growth is driven by the synthesis of proteins, genes, and other cellular components. Defining processes that limit biosynthesis rates is fundamental for understanding the determinants of cell physiology. Here, we analyze the consequences of engineering cells to express extremely high levels of mCherry proteins, as a tool to define limiting processes that fail to adapt upon increasing biosynthetic demands. Protein-burdened cells were transcriptionally and phenotypically similar to mutants of the Mediator, a transcription coactivator complex. However, our binding data suggest that the Mediator was not depleted from endogenous promoters. Burdened cells showed an overall increase in the abundance of the majority of endogenous transcripts, except for highly expressed genes. Our results, supported by mathematical modeling, suggest that wild-type cells transcribe highly expressed genes at the maximal possible rate, as defined by the transcription machinery’s physical properties. We discuss the possible cellular benefit of maximal transcription rates to allow a coordinated optimization of cell size and cell growth.

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.

SummaryThe extent to which carbon flux is directed toward fermentation versus respiration differs between cell types and environmental conditions. Understanding the basic cellular processes governing carbon flux is challenged by the complexity of the metabolic and regulatory networks. To reveal the genetic basis for natural diversity in channeling carbon flux, we applied quantitative trait loci analysis by phenotyping and genotyping hundreds of individual F2 segregants of budding yeast that differ in their capacity to ferment the pentose sugar xylulose. Causal alleles were mapped to the RXT3 and PHO23 genes, two components of the large Rpd3 histone deacetylation complex. We show that these allelic variants modulate the expression of SNF1/AMPK-dependent respiratory genes. Our results suggest that over close evolutionary distances, diversification of carbon flow is driven by changes in global regulators, rather than adaptation of specific metabolic nodes. Such regulators may improve the ability to direct metabolic fluxes for biotechnological applications.

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.

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., Author summary Cells must have sufficient nutrients in order to perform their different functions. This becomes particularly critical when the nutrient composition in their environment changes. To accommodate such changes, cells have evolved specialized adaptation pathways that sense the nutrients available in their surrounding and adapt their gene expression and protein functions accordingly. It is typically assumed that these different starvation pathways have been optimized to enable growth in conditions in which nutrients are limiting. We reasoned, however, that since starvation periods are necessarily followed by nutrient retrieval, the starvation response may further function to optimize the recovery from starvation once nutrients are replenished. In the present study, we show that this is indeed the case for the well-studied phosphate starvation pathway in budding yeast. Using a combination of transcription profiling and high-resolution competition assays, we have characterized, at high temporal resolution, the growth-fitness advantage cells gain from activation of the phosphate starvation program when they enter into starvation and also when they exit starvation when nutrient is replenished. We find that the cells activate the starvation program in 2 temporal waves. The first wave is activated at intermediate depletion that is not yet limiting for growth, while the second wave is activated when phosphate becomes limiting. Of note, we provide evidence that the early response promoted recovery from starvation by increasing phosphate influx. Our results corroborate the notion that cells have evolved not only to optimize their instantaneous growth but also to prepare for future, predicted challenges.

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. We describe a budding yeast hybrid that grows faster than 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 increasing signs of DNA damage. A systematic genetic screen identified hundreds of alleles affecting hybrid growth whose identity vastly differed between the hybrid and its parent and between growth conditions. This large-scale rewiring of allele effects suggests that despite showing clear heterosis, the hybrid is perturbed in multiple regulatory processes. We discuss the possibility that incompatibilities contribute to hybrid vigor by perturbing safeguard mechanisms that limit growth in the parental background.

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.

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.

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.

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.

Doubling DNA but not expression As the genome replicates, and before the cell divides, the copy number of the replicated portions of the genome doubles. In bacteria and archaea, gene expression tracks with gene dosage, both of which increase after DNA replication. Voichek et al. , however, show that an increase in DNA dosage after replication does not increase gene expression in budding yeast. This expression buffering is mediated by the acetylation of newly synthesized histone H3 deposited on the replicated DNA. This acetylation helps suppress transcription from the excess DNA. Science , this issue p. 1087

Transcription factors (TFs) bind to specific DNA sequences to induce or repress gene expression. Expression levels can be tuned by changing TF concentrations, but the precision of such tuning is limited, since the fraction of time a TF occupies its binding site is subject to stochastic fluctuations. Bicoid (Bcd) is a TF that patterns the early Drosophila embryo by establishing an anterior-to-posterior concentration gradient and activating specific gene targets ("gap genes") in a concentration-dependent manner. Recently, the Bcd gradient and its in-vivo diffusion were quantified in live embryos, raising a quandary: the precision by which the Bcd target genes are defined (one-cell resolution) appeared to exceed the physical limits set by the stochastic binding of Bcd to DNA. We hypothesize that early readout of Bcd could account for the observed precision. Specifically, we consider the possibility that gap genes begin to be expressed earlier than typically measured experimentally, at a time when the distance between the nuclei is large. At this time, the difference in Bcd concentration between adjacent nuclei is large, enabling better tolerance for measurement imprecision. We show that such early decoding can indeed increase the accuracy of gap-gene expression, and that the initial pattern can be stabilized during subsequent divisions.

Yeast cells approach a mating partner by polarizing along a gradient of mating pheromones that are secreted by cells of the opposite mating type. The Bar1 protease is secreted by a-cells and, paradoxically, degrades the α-factor pheromones which are produced by cells of the opposite mating type and trigger mating in a-cells. This degradation may assist in the recovery from pheromone signaling but has also been shown to play a positive role in mating. Previous studies suggested that widely diffusing protease can bias the pheromone gradient towards the closest secreting cell. Here, we show that restricting the Bar1 protease to the secreting cell itself, preventing its wide diffusion, facilitates discrimination between equivalent mating partners. This may be mostly relevant during spore germination, where most mating events occur in nature.

We describe the DNA replication timing programs of 14 yeast mutants with an extended S phase identified by a novel genome-wide screen. These mutants are associated with the DNA replication machinery, cell-cycle control, and dNTP synthesis and affect different parts of S phase. In 13 of the mutants, origin activation time scales with the duration of S phase. A limited number of origins become inactive in these strains, with inactive origins characterized by small replicons and distributed throughout S phase. In sharp contrast, cells deleted of MRC1, a gene implicated in replication fork stabilization and in the replication checkpoint pathway, maintained wild-type firing times despite over twofold lengthening of S phase. Numerous dormant origins were activated in this mutant. Our data suggest that most perturbations that lengthen S phase affect the entire program of replication timing, rather than a specific subset of origins, maintaining the relative order of origin firing time and delaying firing with relative proportions. Mrc1 emerges as a regulator of this robustness of the replication program.

We show that, in yeast, the divergence rate of gene expression is not correlated with that of its associated coding sequence. Gene essentiality influences both modes of evolution, but other properties related to protein structure or promoter composition are only correlated with coding-sequence divergence or gene expression divergence, respectively. Based on these findings, we discuss the possibilities of neutral evolution of gene expression and of different modes of evolution in unicellular versus multicellular organisms.

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.

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 s-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 doublecortin-like (DCX) domains serve as protein-interaction platforms. DCX tandem domains appear in the product of the X-linked doublecortin (DCX) gene, in retinitis pigmentosa-1 (RP1), as well as in other gene products. Mutations in the human DCX gene are associated with abnormal neuronal migration, epilepsy, and mental retardation; mutations in RP1 are associated with a form of inherited blindness, while DCDC2 has been associated with dyslectic reading disabilities. Motivated by the possible importance of this gene family, a thorough analysis to detect all family members in the mouse was conducted. The DCX-repeat gene superfamily is composed of eleven paralogs, and we cloned the DCX domains from nine different genes. Our study questioned which functions attributed to the DCX domain, are conserved among the different members. Our results suggest that the proteins with the DCX-domain have conserved and unique roles in microtubule regulation and signal transduction. All the tested proteins stimulated microtubule assembly in vitro. Proteins with tandem repeats stabilized the microtubule cytoskeleton in transfected cells, while those with single repeats localized to actin-rich subcellular structures, or the nucleus. All tested proteins interacted with components of the JNK/MAP-kinase pathway, while only a subset interacted with Neurabin 2, and a nonoverlapping group demonstrated actin association. The sub-specialization of some members due to confined intracellular localization, and protein interactions may explain the success of this superfamily.

Generating clones of mutated cells within a wild-type tissue is a powerful experimental paradigm for elucidating gene function. Recently, this approach was employed for identifying genes that shape morphogen profiles in the Drosophila wing-imaginal disc. Interpreting such experiments poses a theoretical challenge. We present a general framework that links specific features of the morphogen profile in the clone vicinity to three basic morphogen properties: diffusion, degradation, and binding to immobile elements. Our results provide rigorous criteria to examine existing data and can facilitate the design and interpretation of future clone experiments.

Background In budding yeast, perturbations that prolong S phase lead to a proportionate delay in the activation times of most origins. The DNA replication checkpoint was implicated in this scaling phenotype, as an intact checkpoint was shown to be required for the delayed activation of late origins in response to hydroxyurea treatment. In support of that, scaling is lost in cells deleted of mrc1, a mediator of the replication checkpoint signal. Mrc1p, however, also plays a role in normal replication. Results To examine whether the replication checkpoint is required for scaling the replication profile with S phase duration we measured the genome-wide replication profile of different MRC1 alleles that separate its checkpoint function from its role in normal replication, and further analyzed the replication profiles of S phase mutants that are checkpoint deficient. We found that the checkpoint is not required for scaling; rather the unique replication phenotype of mrc1 deleted cells is attributed to the role of Mrc1 in normal replication. This is further supported by the replication profiles of tof1Δ which functions together with Mrc1p in normal replication, and by the distinct replication profiles of specific POL2 alleles which differ in their interaction with Mrc1p. Conclusions We suggest that the slow fork progression in mrc1 deleted cells reduces the likelihood of passive replication leading to the activation of origins that remain mostly dormant in wild-type cells. Electronic supplementary material The online version of this article (doi:10.1186/s12915-014-0079-z) contains supplementary material, which is available to authorized users.

Cells adapt to environmental changes through genetic mutations that stabilize novel phenotypes. Often, this adaptation involves regulatory changes which modulate gene expression. In the budding yeast, ribosomal-related gene expression correlates with cell growth rate across different environments. To examine whether the same relationship between gene expression and growth rate is observed also across natural populations, we measured gene expression, growth rate and ethanol production of twenty-four wild type yeast strains originating from diverse habitats, grown on the pentose sugar xylulose. We found that expression of ribosome-related genes did not correlate with growth rate. Rather, growth rate was correlated with the expression of amino acid biosynthesis genes. Searching other databases, we observed a similar correlation between growth rate and amino-acid biosyntehsis genes in a library of gene deletions. We discuss the implications of our results for understanding how cells coordinate their translation capacity with available nutrient resources.

Various histone modifications are widely associated with gene expression, but their functional selectivity at individual genes remains to be characterized. Here, we identify widespread differences between genome-wide patterns of two prominent marks, H3K9ac and H3K4me3, in budding yeasts. As well as characteristic gene profiles, relative modification levels vary significantly amongst genes, irrespective of expression. Interestingly, we show that these differences couple to contrasting features: higher methylation to essential, periodically expressed, 'DPN' (Depleted Proximal Nucleosome) genes, and higher acetylation to non-essential, responsive, 'OPN' (Occupied Proximal Nucleosome) genes. Thus, H3K4me3 may generally associate with expression stability, and H3K9ac, with variability. To evaluate this notion, we examine their association with expression divergence between the closely related species, S. cerevisiae and S. paradoxus. Although individually well conserved at orthologous genes, changes between modifications are mostly uncorrelated, indicating largely non-overlapping regulatory mechanisms. Notably, we find that inter-species differences in methylation, but not acetylation, are well correlated with expression changes, thereby proposing H3K4me3 as a candidate regulator of expression divergence. Taken together, our results suggest distinct evolutionary roles for expression-linked modifications, wherein H3K4me3 may contribute to stabilize average expression, whilst H3K9ac associates with more indirect aspects such as responsiveness.

Morphogen gradients are used to pattern a field of cells according to variations in the concentration of a signaling molecule. Typically, the morphogen emanates from a confined group of cells. During early embryogenesis, however, the ability to define a restricted source for morphogen production is limited. Thus, various early patterning systems rely on a broadly expressed morphogen that generates an activation gradient within its expression domain. Computational and experimental work has shed light on how a sharp and robust gradient can be established under those situations, leading to a mechanism termed 'morphogen shuttling'. This mechanism relies on an extracellular shuttling molecule that forms an inert, highly diffusible complex with the morphogen. Morphogen release from the complex following cleavage of the shuttling molecule by an extracellular protease leads to the accumulation of free ligand at the center of its expression domain and a graded activation of the developmental pathway that decreases significantly even within the morphogen-expression domain.

When stressed by metal depletion, budding yeast adopt an asymmetric division pattern whereby vacuoles are maintained within dividing mother cells while the vacuoles-deprived daughter cells arrest division. This linear growth mode represents a bet-hedging strategy beneficial at the population level., Budding yeast restricts division to a subpopulation of mother cells when metal is depleted. This population splitting into dividing mothers and arrested daughters implements a bet-hedging strategy beneficial for population long-term survival. Proliferation is limited by the availability of vacuoles, which are asymmetrically segregated to mother cells in a WHI5-dependent manner. Asymmetric resource distribution increases population growth under limiting conditions, defining a novel stress-response strategy., We report that when budding yeast are transferred to low-metal environment, they adopt a proliferation pattern in which division is restricted to the subpopulation of mother cells which were born in rich conditions, before the shift. Mother cells continue to divide multiple times following the shift, generating at each division a single daughter cell, which arrests in G1. The transition to a mother-restricted proliferation pattern is characterized by asymmetric segregation of the vacuole to the mother cell and requires the transcription repressor Whi5. Notably, while deletion of WHI5 alleviates daughter cell division arrest in low-zinc conditions, it results in a lower final population size, as cell division rate becomes progressively slower. Our data suggest a new stress-response strategy, in which the dilution of a limiting cellular resource is prevented by maintaining it within a subset of dividing cells, thereby increasing population growth.

Biological processes are often executed by a small number of molecules per individual cell, leading to significant cell-to-cell variability (“noise”) in gene expression (Paulsson 2004; Maheshri and O'Shea 2007; Raj and van Oudenaarden 2008; Tawfik 2010). When analyzing gene expression noise, it is convenient to distinguish between intrinsic variations, resulting from stochastic production, and extrinsic variations propagating from global (e.g., ribosomes, polymerases, metabolites, etc.) or pathway-specific factors (Elowitz et al. 2002). Intrinsic noise is of particular interest, as it reflects on the transcription process itself (Paulsson 2004, 2005; Raser and O'Shea 2004; Raj and van Oudenaarden 2008; Rinott et al. 2011; So et al. 2011). The prevailing model of gene expression noise assumes that genes transit stochastically between states that are permissive or nonpermissive to transcription (Paulsson 2004, 2005; Raser and O'Shea 2004; Friedman et al. 2006; Raj et al. 2006; Zenklusen et al. 2008; So et al. 2011). This two-state model predicts a scaling relationship between mean expression m and the coefficient of variation (noise, SD/mean) η: η2 = b/m + ηext2 where b is the typical number of protein molecules made during a single “on” state (“burst size”) and ηext denotes the extrinsic noise (Paulsson 2004, 2005; Raser and O'Shea 2004; Bar-Even et al. 2006; Raj et al. 2006; Pedraza and Paulsson 2008; Tan and van Oudenaarden 2010; Taniguchi et al. 2010). Note that burst size, b, accounts for all the transcription-translation processes following the main stochastic event (burst initiation), integrating the number of mRNA molecules produced per burst and the number of protein molecules made per each mRNA molecule. Upon a perturbation, the noise–mean relationship may change, depending on whether burst size or burst frequency were modulated (Pedraza and Paulsson 2008; Zenklusen et al. 2008; Tan and van Oudenaarden 2010). Genome-wide analysis of the noise–mean relationship in yeast (Bar-Even et al. 2006; Newman et al. 2006) or Escherichia coli (Taniguchi et al. 2010; So et al. 2011) genes reported a general dependency that was well defined by the scaling relation η2 = b/m + ηext2, suggesting a similar burst size for many genes. The expression of genes deviating from the scaling curve, displaying higher-than-expected noise (Bar-Even et al. 2006; Newman et al. 2006), was more responsive to changing conditions and also diverged more between related species (Tirosh and Barkai 2008; Choi and Kim 2009; Lehner 2010). Notably, high noise, responsiveness, and divergence were all correlated with the organization of gene promoters: All three measures were low in promoters lacking a TATA box and containing a nucleosome free region (NFR) proximal to the transcription start site (referred to as DPN promoters—depleted proximal nucleosome), and were high in TATA-containing promoters that lack NFR (OPN, occupied proximal nucleosome) (Field et al. 2008; Tirosh and Barkai 2008; Choi and Kim 2009). TATA box was also shown to increase noise in Pho5 expression (Raser and O'Shea 2004) and in synthetic promoters (Blake et al. 2006; Murphy et al. 2010). The observation that genes with a characteristic promoter structure had a high noise (relative to that expected given their mean expression) is consistent with the idea that promoter sequence influences not only burst frequency but also burst size. Still, the principles by which promoter sequences regulate those two processes are not understood, primarily because most studies analyzing the interplay between promoter sequence variations and gene expression consider mean expression only (e.g., Yun et al. 2012). To distinguish between the effects of promoter sequence on burst size and burst frequency, we generated large libraries of sequence-mutated promoters. Specifically, we chose 22 yeast promoters that span a range of expression and noise levels. Using mutagenic PCR, we generated hundreds of sequence variants of each promoter. Each variant was fused to a fluorescent reporter, and the associated mean expression and noise (coefficient of variation) in a population of identical cells were measured. We found that sequence variants in each of the libraries defined a scaling curve η2 = b/m + ηext2, with a constant estimated burst size b, that was largely promoter-specific and was particularly large for OPN promoters containing a TATA box. A small fraction of sequence mutations leading to a large change in burst size was identified in the OPN-type promoters containing a TATA box. These changes were biased toward reducing burst size, and were almost fully explained by elimination of a TATA box or insertion of a new out-of-frame translation start site. Interestingly, mutations that deleted a TATA box in low-noise DPN-type promoters did not reduce burst size. Our results suggest that burst size is a promoter-specific property that is insensitive to most sequence mutations but is largely influenced by the interaction between TATA box and promoter nucleosomes.

Accurate models of the cross-talk between signaling pathways and transcriptional regulatory networks within cells are essential to understand complex response programs. We present a new computational method that combines condition-specific time-series expression data with general protein interaction data to reconstruct dynamic and causal stress response networks. These networks characterize the pathways involved in the response, their time of activation, and the affected genes. The signaling and regulatory components of our networks are linked via a set of common transcription factors that serve as targets in the signaling network and as regulators of the transcriptional response network. Detailed case studies of stress responses in budding yeast demonstrate the predictive power of our method. Our method correctly identifies the core signaling proteins and transcription factors of the response programs. It further predicts the involvement of additional transcription factors and other proteins not previously implicated in the response pathways. We experimentally verify several of these predictions for the osmotic stress response network. Our approach requires little condition-specific data: only a partial set of upstream initiators and time-series gene expression data, which are readily available for many conditions and species. Consequently, our method is widely applicable and can be used to derive accurate, dynamic response models in several species.

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 twenty-two 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.

Closely related species show a high degree of differences in gene expression, but the functional significance of these differences remains unclear. Similarly, stress responses in yeast typically involve differential expression of numerous genes, and it is unclear how many of these are functionally significant. To address these issues, we compared the expression programs of four yeast species under different growth conditions, and found that the response of these species to stress has diverged extensively. On an individual gene basis, most transcriptional responses are not conserved in any pair of species, and there are very limited common responses among all four species. We present evidence that many evolutionary changes in stress responses are compensated either ( i ) by the response of related genes or ( ii ) by changes in the basal expression levels of the genes whose responses have diverged. Thus, stress-related genes are often induced upon stress in some species but maintain high levels even in the absence of stress at other species, indicating a transition between induced and constitutive activation. In addition, ∼15% of the stress responses are specific to only one of the four species, with no evidence for compensating effects or stress-related annotations, and these may reflect fortuitous regulation that is unimportant for the stress response (i.e., biological noise). Frequent compensatory changes and biological noise may explain how diverged expression responses support similar physiological responses.

To examine the role of nucleosome occupancy in the evolution of gene expression, we measured the genome-wide nucleosome profiles of four yeast species, three belonging to the Saccharomyces sensu stricto lineage and the more distantly related Candida glabrata. Nucleosomes and associated promoter elements at C. glabrata genes are typically shifted upstream by ∼20 bp, compared to their orthologs from sensu stricto species. Nonetheless, all species display the same global organization features first described for Saccharomyces cerevisiae: a stereotypical nucleosome organization along genes and a division of promoters into those that contain or lack a pronounced nucleosome-depleted region (NDR), with the latter displaying a more dynamic pattern of gene expression. Despite this global similarity, however, nucleosome occupancy at specific genes diverged extensively between sensu stricto and C. glabrata orthologs (∼50 million years). Orthologs with dynamic expression patterns tend to maintain their lack of NDR, but apart from that, sensu stricto and C. glabrata orthologs are nearly as similar in nucleosome occupancy patterns as nonorthologous genes. This extensive divergence in nucleosome occupancy contrasts with a conserved pattern of gene expression. Thus, while some evolutionary changes in nucleosome occupancy contribute to gene expression divergence, nucleosome occupancy often diverges extensively with apparently little impact on 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., Author Summary 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.

Individuals of the same or closely related species can vary substantially in size. Still, the proportions within and between tissues are precisely kept. This adaptation of pattern with size termed scaling, is receiving a growing attention. We review experimental evidence for scaling, and describe theoretical models for mechanisms that scale morphogen gradients. We particularly note the Expansion-Repression mechanism, in which a diffusible molecule that positively regulates the morphogen gradient width is repressed by morphogen signaling. The Expansion-Repression circuit provides scaling in a robust manner and is readily implemented by a host of molecular mechanisms. We suggest means for identifying such a circuit in a system of interest.