Intersecting transcription networks constrain gene regulatory evolution

DNA sequence and chromatin differentiate sequence-specific transcription factor binding in the human malaria parasite Plasmodium falciparum

Development of the malaria parasite, Plasmodium falciparum, is regulated by a limited number of sequence-specific transcription factors (TFs). However, the mechanisms by which these TFs recognize genome-wide binding sites is largely unknown. To address TF specificity, we investigated the binding of two TF subsets that either bind CACACA or GTGCAC DNA sequence motifs and further characterized two additional ApiAP2 TFs, PfAP2-G and PfAP2-EXP, which bind unique DNA motifs (GTAC and TGCATGCA). We also interrogated the impact of DNA sequence and chromatin context on P. falciparum TF binding by integrating high-throughput in vitro and in vivo binding assays, DNA shape predictions, epigenetic post-translational modifications, and chromatin accessibility. We found that DNA sequence context minimally impacts binding site selection for paralogous CACACA-binding TFs, while chromatin accessibility, epigenetic patterns, co-factor recruitment, and dimerization correlate with differential binding. In contrast, GTGCAC-binding TFs prefer different DNA sequence context in addition to chromatin dynamics. Finally, we determined that TFs that preferentially bind divergent DNA motifs may bind overlapping genomic regions due to low-affinity binding to other sequence motifs. Our results demonstrate that TF binding site selection relies on a combination of DNA sequence and chromatin features, thereby contributing to the complexity of P. falciparum gene regulatory mechanisms.

Prevalence of and gene regulatory constraints on transcriptional adaptation in single cells

BACKGROUND

Cells and tissues have a remarkable ability to adapt to genetic perturbations via a variety of molecular mechanisms. Nonsense-induced transcriptional compensation, a form of transcriptional adaptation, has recently emerged as one such mechanism, in which nonsense mutations in a gene trigger upregulation of related genes, possibly conferring robustness at cellular and organismal levels. However, beyond a handful of developmental contexts and curated sets of genes, no comprehensive genome-wide investigation of this behavior has been undertaken for mammalian cell types and conditions. How the regulatory-level effects of inherently stochastic compensatory gene networks contribute to phenotypic penetrance in single cells remains unclear.

RESULTS

We analyze existing bulk and single-cell transcriptomic datasets to uncover the prevalence of transcriptional adaptation in mammalian systems across diverse contexts and cell types. We perform regulon gene expression analyses of transcription factor target sets in both bulk and pooled single-cell genetic perturbation datasets. Our results reveal greater robustness in expression of regulons of transcription factors exhibiting transcriptional adaptation compared to those of transcription factors that do not. Stochastic mathematical modeling of minimal compensatory gene networks qualitatively recapitulates several aspects of transcriptional adaptation, including paralog upregulation and robustness to mutation. Combined with machine learning analysis of network features of interest, our framework offers potential explanations for which regulatory steps are most important for transcriptional adaptation.

CONCLUSIONS

Our integrative approach identifies several putative hits-genes demonstrating possible transcriptional adaptation-to follow-up on experimentally and provides a formal quantitative framework to test and refine models of transcriptional adaptation.

The transcription factor STE12 influences growth on several carbon sources and production of dehydroacetic acid (DHAA) in Trichoderma reesei

The filamentous ascomycete Trichoderma reesei, known for its prolific cellulolytic enzyme production, recently also gained attention for its secondary metabolite synthesis. Both processes are intricately influenced by environmental factors like carbon source availability and light exposure. Here, we explore the role of the transcription factor STE12 in regulating metabolic pathways in T. reesei in terms of gene regulation, carbon source utilization and biosynthesis of secondary metabolites. We show that STE12 is involved in regulating cellulase gene expression and growth on carbon sources associated with iron homeostasis. STE12 impacts gene regulation in a light dependent manner on cellulose with modulation of several CAZyme encoding genes as well as genes involved in secondary metabolism. STE12 selectively influences the biosynthesis of the sorbicillinoid trichodimerol, while not affecting the biosynthesis of bisorbibutenolide, which was recently shown to be regulated by the MAPkinase pathway upstream of STE12 in the signaling cascade. We further report on the biosynthesis of dehydroacetic acid (DHAA) in T. reesei, a compound known for its antimicrobial properties, which is subject to regulation by STE12. We conclude, that STE12 exerts functions beyond development and hence contributes to balance the energy distribution between substrate consumption, reproduction and defense.

Dissecting an ancient stress resistance trait syndrome in the compost yeast Kluyveromyces marxianus

In the search to understand how evolution builds new traits, ancient events are often the hardest to dissect. Species-unique traits pose a particular challenge for geneticists-cases in which a character arose long ago and, in the modern day, is conserved within a species, distinguishing it from reproductively isolated relatives. In this work, we have developed the budding yeast genus Kluyveromyces as a model for mechanistic dissection of trait variation across species boundaries. Phenotypic profiling revealed robust heat and chemical-stress tolerance phenotypes that distinguished the compost yeast K. marxianus from the rest of the clade. We used culture-based, transcriptomic, and genetic approaches to characterize the metabolic requirements of the K. marxianus trait syndrome. We then generated a population-genomic resource for K. marxianus and harnessed it in molecular-evolution analyses, which found hundreds of housekeeping genes with evidence for adaptive protein variation unique to this species. Our data support a model in which, in the distant past, K. marxianus underwent a vastly complex remodeling of its proteome to achieve stress resistance. Such a polygenic architecture, involving nucleotide-level allelic variation on a massive scale, is consistent with theoretical models of the mechanisms of long-term adaptation, and suggests principles of broad relevance for interspecies trait genetics.

Prevalence of and gene regulatory constraints on transcriptional adaptation in single cells

Cells and tissues have a remarkable ability to adapt to genetic perturbations via a variety of molecular mechanisms. Nonsense-induced transcriptional compensation, a form of transcriptional adaptation, has recently emerged as one such mechanism, in which nonsense mutations in a gene can trigger upregulation of related genes, possibly conferring robustness at cellular and organismal levels. However, beyond a handful of developmental contexts and curated sets of genes, to date, no comprehensive genome-wide investigation of this behavior has been undertaken for mammalian cell types and contexts. Moreover, how the regulatory-level effects of inherently stochastic compensatory gene networks contribute to phenotypic penetrance in single cells remains unclear. Here we combine computational analysis of existing datasets with stochastic mathematical modeling and machine learning to uncover the widespread prevalence of transcriptional adaptation in mammalian systems and the diverse single-cell manifestations of minimal compensatory gene networks. Regulon gene expression analysis of a pooled single-cell genetic perturbation dataset recapitulates important model predictions. Our integrative approach uncovers several putative hits-genes demonstrating possible transcriptional adaptation-to follow up on experimentally, and provides a formal quantitative framework to test and refine models of transcriptional adaptation.

Evolutionary adaptation to juvenile malnutrition impacts adult metabolism and impairs adult fitness in Drosophila

Juvenile undernutrition has lasting effects on adult metabolism of the affected individuals, but it is unclear how adult physiology is shaped over evolutionary time by natural selection driven by juvenile undernutrition. We combined RNAseq, targeted metabolomics, and genomics to study the consequences of evolution under juvenile undernutrition for metabolism of reproductively active adult females of Drosophila melanogaster. Compared to Control populations maintained on standard diet, Selected populations maintained for over 230 generations on a nutrient-poor larval diet evolved major changes in adult gene expression and metabolite abundance, in particular affecting amino acid and purine metabolism. The evolved differences in adult gene expression and metabolite abundance between Selected and Control populations were positively correlated with the corresponding differences previously reported for Selected versus Control larvae. This implies that genetic variants affect both stages similarly. Even when well fed, the metabolic profile of Selected flies resembled that of flies subject to starvation. Finally, Selected flies had lower reproductive output than Controls even when both were raised under the conditions under which the Selected populations evolved. These results imply that evolutionary adaptation to juvenile undernutrition has large pleiotropic consequences for adult metabolism, and that they are costly rather than adaptive for adult fitness. Thus, juvenile and adult metabolism do not appear to evolve independently from each other even in a holometabolous species where the two life stages are separated by a complete metamorphosis.

Extending the reach of homology by using successive computational filters to find yeast pheromone genes

The mating of fungi depends on pheromones that mediate communication between two mating types. Most species use short peptides as pheromones, which are either unmodified (e.g., α-factor in Saccharomyces cerevisiae) or C-terminally farnesylated (e.g., a-factor in S. cerevisiae). Peptide pheromones have been found by genetics or biochemistry in a small number of fungi, but their short sequences and modest conservation make it impossible to detect homologous sequences in most species. To overcome this problem, we used a four-step computational pipeline to identify candidate a-factor genes in sequenced genomes of the Saccharomycotina, the fungal clade that contains most of the yeasts: we require that candidate genes have a C-terminal prenylation motif, are shorter than 100 amino acids long, and contain a proteolytic-processing motif upstream of the potential mature pheromone sequence and that closely related species contain highly conserved homologs of the potential mature pheromone sequence. Additional manual curation exploits the observation that many species carry more than one a-factor gene, encoding identical or nearly identical pheromones. From 332 Saccharomycotina genomes, we identified strong candidate pheromone genes in 241 genomes, covering 13 clades that are each separated from each other by at least 100 million years, the time required for evolution to remove detectable sequence homology among small pheromone genes. For one small clade, the Yarrowia, we demonstrated that our algorithm found the a-factor genes: deleting all four related genes in the a-mating type of Yarrowia lipolytica prevents mating.

Cooperative assembly confers regulatory specificity and long-term genetic circuit stability

A ubiquitous feature of eukaryotic transcriptional regulation is cooperative self-assembly between transcription factors (TFs) and DNA cis-regulatory motifs. It is thought that this strategy enables specific regulatory connections to be formed in gene networks between otherwise weakly interacting, low-specificity molecular components. Here, using synthetic gene circuits constructed in yeast, we find that high regulatory specificity can emerge from cooperative, multivalent interactions among artificial zinc-finger-based TFs. We show that circuits "wired" using the strategy of cooperative TF assembly are effectively insulated from aberrant misregulation of the host cell genome. As we demonstrate in experiments and mathematical models, this mechanism is sufficient to rescue circuit-driven fitness defects, resulting in genetic and functional stability of circuits in long-term continuous culture. Our naturally inspired approach offers a simple, generalizable means for building high-fidelity, evolutionarily robust gene circuits that can be scaled to a wide range of host organisms and applications.

Evolution of a new form of haploid-specific gene regulation appearing in a limited clade of ascomycete yeast species

Over evolutionary timescales, the logic and pattern of cell-type specific gene expression can remain constant, yet the molecular mechanisms underlying such regulation can drift between alternative forms. Here, we document a new example of this principle in the regulation of the haploid-specific genes in a small clade of fungal species. For most ascomycete fungal species, transcription of these genes is repressed in the a/α cell type by a heterodimer of two homeodomain proteins, Mata1 and Matα2. We show that in the species Lachancea kluyveri, most of the haploid-specific genes are regulated in this way, but repression of one haploid-specific gene (GPA1) requires, in addition to Mata1 and Matα2, a third regulatory protein, Mcm1. Model building, based on x-ray crystal structures of the three proteins, rationalizes the requirement for all three proteins: no single pair of the proteins is optimally arranged, and we show that no single pair can bring about repression. This case study exemplifies the idea that the energy of DNA binding can be "shared out" in different ways and can result in different DNA-binding solutions across different genes-while maintaining the same overall pattern of gene expression.

Robustness and innovation in synthetic genotype networks

Genotype networks are sets of genotypes connected by small mutational changes that share the same phenotype. They facilitate evolutionary innovation by enabling the exploration of different neighborhoods in genotype space. Genotype networks, first suggested by theoretical models, have been empirically confirmed for proteins and RNAs. Comparative studies also support their existence for gene regulatory networks (GRNs), but direct experimental evidence is lacking. Here, we report the construction of three interconnected genotype networks of synthetic GRNs producing three distinct phenotypes in Escherichia coli. Our synthetic GRNs contain three nodes regulating each other by CRISPR interference and governing the expression of fluorescent reporters. The genotype networks, composed of over twenty different synthetic GRNs, provide robustness in face of mutations while enabling transitions to innovative phenotypes. Through realistic mathematical modeling, we quantify robustness and evolvability for the complete genotype-phenotype map and link these features mechanistically to GRN motifs. Our work thereby exemplifies how GRN evolution along genotype networks might be driving evolutionary innovation.

Phase separation in fungi

Phase separation, in which macromolecules partition into a concentrated phase that is immiscible with a dilute phase, is involved with fundamental cellular processes across the tree of life. We review the principles of phase separation and highlight how it impacts diverse processes in the fungal kingdom. These include the regulation of autophagy, cell signalling pathways, transcriptional circuits and the establishment of asymmetry in fungal cells. We describe examples of stable, phase-separated assemblies including membraneless organelles such as the nucleolus as well as transient condensates that also arise through phase separation and enable cells to rapidly and reversibly respond to important environmental cues. We showcase how research into phase separation in model yeasts, such as Saccharomyces cerevisiae and Schizosaccharomyces pombe, in conjunction with that in plant and human fungal pathogens, such as Ashbya gossypii and Candida albicans, is continuing to enrich our understanding of fundamental molecular processes.

What makes Komagataella phaffii non-conventional?

The important industrial protein production host Komagataella phaffii (syn Pichia pastoris) is classified as a non-conventional yeast. But what exactly makes K. phaffii non-conventional? In this review, we set out to address the main differences to the 'conventional' yeast Saccharomyces cerevisiae, but also pinpoint differences to other non-conventional yeasts used in biotechnology. Apart from its methylotrophic lifestyle, K. phaffii is a Crabtree-negative yeast species. But even within the methylotrophs, K. phaffii possesses distinct regulatory features such as glycerol-repression of the methanol-utilization pathway or the lack of nitrate assimilation. Rewiring of the transcriptional networks regulating carbon (and nitrogen) source utilization clearly contributes to our understanding of genetic events occurring during evolution of yeast species. The mechanisms of mating-type switching and the triggers of morphogenic phenotypes represent further examples for how K. phaffii is distinguished from the model yeast S. cerevisiae. With respect to heterologous protein production, K. phaffii features high secretory capacity but secretes only low amounts of endogenous proteins. Different to S. cerevisiae, the Golgi apparatus of K. phaffii is stacked like in mammals. While it is tempting to speculate that Golgi architecture is correlated to the high secretion levels or the different N-glycan structures observed in K. phaffii, there is recent evidence against this. We conclude that K. phaffii is a yeast with unique features that has a lot of potential to explore both fundamental research questions and industrial applications.

Insights on life cycle and cell identity regulatory circuits for unlocking genetic improvement in Zygosaccharomyces and Kluyveromyces yeasts

Evolution has provided a vast diversity of yeasts that play fundamental roles in nature and society. This diversity is not limited to genotypically homogeneous species with natural interspecies hybrids and allodiploids that blur species boundaries frequently isolated. Thus, life cycle and the nature of breeding systems have profound effects on genome variation, shaping heterozygosity, genotype diversity and ploidy level. The apparent enrichment of hybrids in industry-related environments suggests that hybridization provides an adaptive route against stressors and creates interest in developing new hybrids for biotechnological uses. For example, in the Saccharomyces genus where regulatory circuits controlling cell identity, mating competence and meiosis commitment have been extensively studied, this body of knowledge is being used to combine interesting traits into synthetic F1 hybrids, to bypass F1 hybrid sterility and to dissect complex phenotypes by bulk segregant analysis. Although these aspects are less known in other industrially promising yeasts, advances in whole-genome sequencing and analysis are changing this and new insights are being gained, especially in the food-associated genera Zygosaccharomyces and Kluyveromyces. We discuss this new knowledge and highlight how deciphering cell identity circuits in these lineages will contribute significantly to identify the genetic determinants underpinning complex phenotypes and open new avenues for breeding programmes.

Differences in evolutionary accessibility determine which equally effective regulatory motif evolves to generate pulses

Transcriptional regulatory networks (TRNs) are enriched for certain "motifs." Motif usage is commonly interpreted in adaptationist terms, i.e., that the optimal motif evolves. But certain motifs can also evolve more easily than others. Here, we computationally evolved TRNs to produce a pulse of an effector protein. Two well-known motifs, type 1 incoherent feed-forward loops (I1FFLs) and negative feedback loops (NFBLs), evolved as the primary solutions. The relative rates at which these two motifs evolve depend on selection conditions, but under all conditions, either motif achieves similar performance. I1FFLs generally evolve more often than NFBLs. Selection for a tall pulse favors NFBLs, while selection for a fast response favors I1FFLs. I1FFLs are more evolutionarily accessible early on, before the effector protein evolves high expression; when NFBLs subsequently evolve, they tend to do so from a conjugated I1FFL-NFBL genotype. In the empirical S. cerevisiae TRN, output genes of NFBLs had higher expression levels than those of I1FFLs. These results suggest that evolutionary accessibility, and not relative functionality, shapes which motifs evolve in TRNs, and does so as a function of the expression levels of particular genes.

Regulation of mating and mating-type-specific genes in Zygosaccharomyces sp. yeast

Zygosaccharomyces sp. is an industrially important yeast for the production traditional fermented foods in Japan. At present, however, there is no easy method for mating Zygosaccharomyces sp. strains in the laboratory; furthermore, little is known about the expression of mating-type-specific genes in this yeast. Here, mating was observed when Zygosaccharomyces sp. was subjected to nitrogen-starvation conditions. The expression of mating-type-specific genes, Zygo STE6 and Zygo MFα1, was induced under nitrogen-starvation conditions, as confirmed by lacZ reporter assay. This expression was mating-type-specific: Zygo STE6 was expressed specifically for mating-type a, whereas and Zygo MFα1 was expressed specifically for mating-type α. Yeast strains Zygosaccharomyces rouxii DL2 and DA2, derived from type strain Z. rouxii CBS732, did not show mating even under nitrogen-starvation conditions. Gene sequencing revealed that the Zygo STE12 in Z. rouxii CBS732 has a frameshift mutation. Under nitrogen starvation, mating was observed in both DL2 and DA2 transformed with the wild-type Zygo STE12. The expression of Zygo STE6 in Z. rouxii DL2 transformed with wild-type Zygo STE12 under nitrogen-starvation conditions was confirmed by lacZ reporter assay. Collectively, these results revealed that, under nitrogen-starvation conditions, Zygosaccharomyces sp. can mate and mating-type-specific genes are expressed. Furthermore, Zygo Ste12 is essential for both mating and the expression of mating-type-specific genes in Zygosaccharomyces sp.

The Arabidopsis zinc finger proteins SRG2 and SRG3 are positive regulators of plant immunity and are differentially regulated by nitric oxide

Nitric oxide (NO) regulates the deployment of a phalanx of immune responses, chief among which is the activation of a constellation of defence-related genes. However, the underlying molecular mechanisms remain largely unknown. The Arabidopsis thaliana zinc finger transcription factor (ZF-TF), S-nitrosothiol (SNO) Regulated 1 (SRG1), is a central target of NO bioactivity during plant immunity. Here we characterize the remaining members of the SRG gene family. Both SRG2 and, especially, SRG3 were positive regulators of salicylic acid-dependent plant immunity. Analysis of SRG single, double and triple mutants implied that SRG family members have additive functions in plant immunity and, surprisingly, are under reciprocal regulation. SRG2 and SRG3 localized to the nucleus and functioned as ethylene-responsive element binding factor-associated amphiphilic repression (EAR) domain-dependent transcriptional repressors: NO abolished this activity for SRG3 but not for SRG2. Consistently, loss of GSNOR function, resulting in increased (S)NO concentrations, fully suppressed the disease resistance phenotype established from SRG3 but not SRG2 overexpression. Remarkably, SRG3 but not SRG2 was S-nitrosylated in vitro and in vivo. Our findings suggest that the SRG family has separable functions in plant immunity, and, surprisingly, these ZF-TFs exhibit reciprocal regulation. It is remarkable that, through neofunctionalization, the SRG family has evolved to become differentially regulated by the key immune-related redox cue, NO.

Lineage-specific selection and the evolution of virulence in the Candida clade

Of the many microbial species on earth, only a small number are able to thrive in humans and cause disease. Comparison of closely related pathogenic and nonpathogenic species can therefore be useful in identifying key features that contribute to virulence. We created interspecies hybrids between Candida albicans, a prevalent fungal pathogen of humans, and Candida dubliniensis, a close, but much less pathogenic, relative. By comparing genome-wide expression differences between the two genomes in the same cell, we surmised that since the two species diverged from a common ancestor, natural selection has acted upon the expression level of an ancient metabolic pathway, illustrating that pathogenicity traits can arise over evolutionary timescales through small expression changes in deeply conserved proteins.

Candida albicans is the most common cause of systemic fungal infections in humans and is considerably more virulent than its closest known relative, Candida dubliniensis. To investigate this difference, we constructed interspecies hybrids and quantified mRNA levels produced from each genome in the hybrid. This approach systematically identified expression differences in orthologous genes arising from cis-regulatory sequence changes that accumulated since the two species last shared a common ancestor, some 10 million y ago. We documented many orthologous gene-expression differences between the two species, and we pursued one striking observation: All 15 genes coding for the enzymes of glycolysis showed higher expression from the C. albicans genome than the C. dubliniensis genome in the interspecies hybrid. This pattern requires evolutionary changes to have occurred at each gene; the fact that they all act in the same direction strongly indicates lineage-specific natural selection as the underlying cause. To test whether these expression differences contribute to virulence, we created a C. dubliniensis strain in which all 15 glycolysis genes were produced at modestly elevated levels and found that this strain had significantly increased virulence in the standard mouse model of systemic infection. These results indicate that small expression differences across a deeply conserved set of metabolism enzymes can play a significant role in the evolution of virulence in fungal pathogens.

Crossbreeding of Yeasts Domesticated for Fermentation: Infertility Challenges

Sexual reproduction is almost a universal feature of eukaryotic organisms, which allows the reproduction of new organisms by combining the genetic information from two individuals of different sexes. Based on the mechanism of sexual reproduction, crossbreeding provides an attractive opportunity to improve the traits of animals, plants, and fungi. The budding yeast Saccharomyces cerevisiae has been widely utilized in fermentative production since ancient times. Currently it is still used for many essential biotechnological processes including the production of beer, wine, and biofuels. It is surprising that many yeast strains used in the industry exhibit low rates of sporulation resulting in limited crossbreeding efficiency. Here, I provide an overview of the recent findings about infertility challenges of yeasts domesticated for fermentation along with the progress in crossbreeding technologies. The aim of this review is to create an opportunity for future crossbreeding of yeasts used for fermentation.

A conceptual framework of evolutionary novelty and innovation

Since 1990 the recognition of deep homologies among metazoan developmental processes and the spread of more mechanistic approaches to developmental biology have led to a resurgence of interest in evolutionary novelty and innovation. Other evolutionary biologists have proposed central roles for behaviour and phenotypic plasticity in generating the conditions for the construction of novel morphologies, or invoked the accessibility of new regions of vast sequence spaces. These approaches contrast with more traditional emphasis on the exploitation of ecological opportunities as the primary source of novelty. This definitional cornucopia reflects differing stress placed on three attributes of novelties: their radical nature, the generation of new taxa, and ecological and evolutionary impact. Such different emphasis has led to conflating four distinct issues: the origin of novel attributes (genes, developmental processes, phenotypic characters), new functions, higher clades and the ecological impact of new structures and functions. Here I distinguish novelty (the origin of new characters, deep character transformations, or new combinations) from innovation, the ecological and evolutionary success of clades. Evidence from the fossil record of macroevolutionary lags between the origin of a novelty and its ecological success demonstrates that novelty may be decoupled from innovation, and only definitions of novelty based on radicality (rather than generativity or consequentiality) can be assessed without reference to the subsequent history of the clade to which a novelty belongs. These considerations suggest a conceptual framework for novelty and innovation, involving: (i) generation of the potential for novelty; (ii) the formation of novel attributes; (iii) refinement of novelties through adaptation; (iv) exploitation of novelties by a clade, which may coincide with a new round of ecological or environmental potentiation; followed by (v) the establishment of innovations through ecological processes. This framework recognizes that there is little empirical support for either the dominance of ecological opportunity, nor abrupt discontinuities (often caricatured as 'hopeful monsters'). This general framework may be extended to aspects of cultural and social innovation.

Predicting Evolutionary Constraints by Identifying Conflicting Demands in Regulatory Networks

Gene regulation networks allow organisms to adapt to diverse environmental niches. However, the constraints underlying the evolution of gene regulation remain ill defined. Here, we show that partial order-a concept that ranks network output levels as a function of different input signals-identifies such constraints. We tested our predictions by experimentally evolving an engineered signal-integrating network in multiple environments. We find that populations: (1) expand in fitness space along the Pareto-optimal front associated with conflicts in regulatory demands, by fine-tuning binding affinities within the network, and (2) expand beyond the Pareto-optimal front through changes in the network structure. Our constraint predictions are based only on partial order and do not require information on the network architecture or underlying genetics. Overall, our findings show that limited knowledge of current regulatory phenotypes can provide predictions on future evolutionary constraints.

Predicting Evolution Using Regulatory Architecture

The limits of evolution have long fascinated biologists. However, the causes of evolutionary constraint have remained elusive due to a poor mechanistic understanding of studied phenotypes. Recently, a range of innovative approaches have leveraged mechanistic information on regulatory networks and cellular biology. These methods combine systems biology models with population and single-cell quantification and with new genetic tools, and they have been applied to a range of complex cellular functions and engineered networks. In this article, we review these developments, which are revealing the mechanistic causes of epistasis at different levels of biological organization-in molecular recognition, within a single regulatory network, and between different networks-providing first indications of predictable features of evolutionary constraint.

Protein-coding changes preceded cis-regulatory gains in a newly evolved transcription circuit

Changes in both the coding sequence of transcriptional regulators and in the cis-regulatory sequences recognized by these regulators have been implicated in the evolution of transcriptional circuits. However, little is known about how they evolved in concert. We describe an evolutionary pathway in fungi where a new transcriptional circuit (a-specific gene repression by the homeodomain protein Matα2) evolved by coding changes in this ancient regulator, followed millions of years later by cis-regulatory sequence changes in the genes of its future regulon. By analyzing a group of species that has acquired the coding changes but not the cis-regulatory sites, we show that the coding changes became necessary for the regulator's deeply conserved function, thereby poising the regulator to jump-start formation of the new circuit.

Identifying the drivers of computationally detected correlated evolution among sites under antibiotic selection

The ultimate causes of correlated evolution among sites in a genome remain difficult to tease apart. To address this problem directly, we performed a high-throughput search for correlated evolution among sites associated with resistance to a fluoroquinolone antibiotic using whole-genome data from clinical strains of Pseudomonas aeruginosa, before validating our computational predictions experimentally. We show that for at least two sites, this correlation is underlain by epistasis. Our analysis also revealed eight additional pairs of synonymous substitutions displaying correlated evolution underlain by physical linkage, rather than selection associated with antibiotic resistance. Our results provide direct evidence that both epistasis and physical linkage among sites can drive the correlated evolution identified by high-throughput computational tools. In other words, the observation of correlated evolution is not by itself sufficient evidence to guarantee that the sites in question are epistatic; such a claim requires additional evidence, ideally coming from direct estimates of epistasis, based on experimental evidence.

G-protein-coupled Receptors in Fungi

G-protein-coupled receptors (GPCRs) are the largest family of transmembrane receptors in fungi. These receptors have an important role in the transduction of extracellular signals into intracellular sites in response to diverse stimuli. They enable fungi to coordinate cell function and metabolism, thereby promoting their survival and propagation, and sense certain fundamentally conserved elements, such as nutrients, pheromones, and stress, for adaptation to their niches, environmental stresses, and host environment, causing disease and pathogen virulence. This chapter highlights the role of GPCRs in fungi in coordinating cell function and metabolism. Fungal cells sense the molecular interactions between extracellular signals. Their respective sensory systems are described here in detail.

Phyllotaxis: is the golden angle optimal for light capture?

Phyllotactic patterns are some of the most conspicuous in nature. To create these patterns plants must control the divergence angle between the appearance of successive organs, sometimes to within a fraction of a degree. The most common angle is the Fibonacci or golden angle, and its prevalence has led to the hypothesis that it has been selected by evolution as optimal for plants with respect to some fitness benefits, such as light capture. We explore arguments for and against this idea with computer models. We have used both idealized and scanned leaves from Arabidopsis thaliana and Cardamine hirsuta to measure the overlapping leaf area of simulated plants after varying parameters such as leaf shape, incident light angles, and other leaf traits. We find that other angles generated by Fibonacci-like series found in nature are equally optimal for light capture, and therefore should be under similar evolutionary pressure. Our findings suggest that the iterative mechanism for organ positioning itself is a more likely target for evolutionary pressure, rather than a specific divergence angle, and our model demonstrates that the heteroblastic progression of leaf shape in A. thaliana can provide a potential fitness benefit via light capture.

The Causes and Consequences of Genetic Interactions (Epistasis)

The same mutation can have different effects in different individuals. One important reason for this is that the outcome of a mutation can depend on the genetic context in which it occurs. This dependency is known as epistasis. In recent years, there has been a concerted effort to quantify the extent of pairwise and higher-order genetic interactions between mutations through deep mutagenesis of proteins and RNAs. This research has revealed two major components of epistasis: nonspecific genetic interactions caused by nonlinearities in genotype-to-phenotype maps, and specific interactions between particular mutations. Here, we provide an overview of our current understanding of the mechanisms causing epistasis at the molecular level, the consequences of genetic interactions for evolution and genetic prediction, and the applications of epistasis for understanding biology and determining macromolecular structures.

Evolution of resilience in protein interactomes across the tree of life

Phenotype robustness to environmental fluctuations is a common biological phenomenon. Although most phenotypes involve multiple proteins that interact with each other, the basic principles of how such interactome networks respond to environmental unpredictability and change during evolution are largely unknown. Here we study interactomes of 1,840 species across the tree of life involving a total of 8,762,166 protein-protein interactions. Our study focuses on the resilience of interactomes to network failures and finds that interactomes become more resilient during evolution, meaning that interactomes become more robust to network failures over time. In bacteria, we find that a more resilient interactome is in turn associated with the greater ability of the organism to survive in a more complex, variable, and competitive environment. We find that at the protein family level proteins exhibit a coordinated rewiring of interactions over time and that a resilient interactome arises through gradual change of the network topology. Our findings have implications for understanding molecular network structure in the context of both evolution and environment.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

Intrinsic cooperativity potentiates parallel cis-regulatory evolution

Convergent evolutionary events in independent lineages provide an opportunity to understand why evolution favors certain outcomes over others. We studied such a case where a large set of genes-those coding for the ribosomal proteins-gained cis-regulatory sequences for a particular transcription regulator (Mcm1) in independent fungal lineages. We present evidence that these gains occurred because Mcm1 shares a mechanism of transcriptional activation with an ancestral regulator of the ribosomal protein genes, Rap1. Specifically, we show that Mcm1 and Rap1 have the inherent ability to cooperatively activate transcription through contacts with the general transcription factor TFIID. Because the two regulatory proteins share a common interaction partner, the presence of one ancestral cis-regulatory sequence can 'channel' random mutations into functional sites for the second regulator. At a genomic scale, this type of intrinsic cooperativity can account for a pattern of parallel evolution involving the fixation of hundreds of substitutions.

Synthetic circuits reveal how mechanisms of gene regulatory networks constrain evolution

Phenotypic variation is the raw material of adaptive Darwinian evolution. The phenotypic variation found in organismal development is biased towards certain phenotypes, but the molecular mechanisms behind such biases are still poorly understood. Gene regulatory networks have been proposed as one cause of constrained phenotypic variation. However, most pertinent evidence is theoretical rather than experimental. Here, we study evolutionary biases in two synthetic gene regulatory circuits expressed in Escherichia coli that produce a gene expression stripe-a pivotal pattern in embryonic development. The two parental circuits produce the same phenotype, but create it through different regulatory mechanisms. We show that mutations cause distinct novel phenotypes in the two networks and use a combination of experimental measurements, mathematical modelling and DNA sequencing to understand why mutations bring forth only some but not other novel gene expression phenotypes. Our results reveal that the regulatory mechanisms of networks restrict the possible phenotypic variation upon mutation. Consequently, seemingly equivalent networks can indeed be distinct in how they constrain the outcome of further evolution.

Genomewide mechanisms of chronological longevity by dietary restriction in budding yeast

Dietary restriction is arguably the most promising nonpharmacological intervention to extend human life and health span. Yet, only few genetic regulators mediating the cellular response to dietary restriction are known, and the question remains which other regulatory factors are involved. Here, we measured at the genomewide level the chronological lifespan of Saccharomyces cerevisiae gene deletion strains under two nitrogen source regimens, glutamine (nonrestricted) and γ‐aminobutyric acid (restricted). We identified 473 mutants with diminished or enhanced extension of lifespan. Functional analysis of such dietary restriction genes revealed novel processes underlying longevity by the nitrogen source quality, which also allowed us to generate a prioritized catalogue of transcription factors orchestrating the dietary restriction response. Importantly, deletions of transcription factors Msn2, Msn4, Snf6, Tec1, and Ste12 resulted in diminished lifespan extension and defects in cell cycle arrest upon nutrient starvation, suggesting that regulation of the cell cycle is a major mechanism of chronological longevity. We further show that STE12 overexpression is enough to extend lifespan, linking the pheromone/invasive growth pathway with cell survivorship. Our global picture of the genetic players of longevity by dietary restriction highlights intricate regulatory cross‐talks in aging cells.

Sign epistasis caused by hierarchy within signalling cascades

Sign epistasis is a central evolutionary constraint, but its causal factors remain difficult to predict. Here we use the notion of parameterised optima to explain epistasis within a signalling cascade, and test these predictions in Escherichia coli. We show that sign epistasis arises from the benefit of tuning phenotypic parameters of cascade genes with respect to each other, rather than from their complex and incompletely known genetic bases. Specifically, sign epistasis requires only that the optimal phenotypic parameters of one gene depend on the phenotypic parameters of another, independent of other details, such as activating or repressing nature, position within the cascade, intra-genic pleiotropy or genotype. Mutational effects change sign more readily in downstream genes, indicating that optimising downstream genes is more constrained. The findings show that sign epistasis results from the inherent upstream-downstream hierarchy between signalling cascade genes, and can be addressed without exhaustive genotypic mapping.

Sign epistasis clearly constrains evolution, but its causes are difficult to decipher. Here, the authors study epistasis in a signalling cascade, and arrive at a general criterion and understanding of sign epistasis as arising from the inherent hierarchy between signalling cascade components.

Flip/flop mating-type switching in the methylotrophic yeast Ogataea polymorpha is regulated by an Efg1-Rme1-Ste12 pathway

In haploid cells of Ogataea (Hansenula) polymorpha an environmental signal, nitrogen starvation, induces a reversible change in the structure of a chromosome. This process, mating-type switching, inverts a 19-kb DNA region to place either MATa or MATα genes under centromeric repression of transcription, depending on the orientation of the region. Here, we investigated the genetic pathway that controls switching. We characterized the transcriptomes of haploid and diploid O. polymorpha by RNAseq in rich and nitrogen-deficient media, and found that there are no constitutively a-specific or α-specific genes other than the MAT genes themselves. We mapped a switching defect in a sibling species (O. parapolymorpha strain DL-1) by interspecies bulk segregant analysis to a frameshift in the transcription factor EFG1, which in Candida albicans regulates filamentous growth and white-opaque switching. Gene knockout, overexpression and ChIPseq experiments show that EFG1 regulates RME1, which in turn regulates STE12, to achieve mating-type switching. All three genes are necessary both for switching and for mating. Overexpression of RME1 or STE12 is sufficient to induce switching without a nitrogen depletion signal. The homologous recombination genes RAD51 and RAD17 are also necessary for switching. The pathway controlling switching in O. polymorpha shares no components with the regulation of HO in S. cerevisiae, which does not involve any environmental signal, but it shares some components with mating-type switching in Kluyveromyces lactis and with white-opaque phenotypic switching in C. albicans.

The molecular mechanisms of self-fertility (homothallism) vary enormously among fungal species. We previously found that in the yeast Ogataea polymorpha, homothallism is achieved by a novel mating-type switching mechanism that exchanges the locations of MATa and MATα genes between expression and repression contexts. Switching in this species is induced by nitrogen depletion, unlike the analogous process in Saccharomyces cerevisiae. Here, we show that the upstream parts of the genetic pathway controlling the environmental induction of switching in O. polymorpha are the same as the environmental pathway that induces competence for mating in this species.

The rewiring of transcription circuits in evolution

The binding of transcription regulators to cis-regulatory sequences is a key step through which all cells regulate expression of their genes. Due to gains and losses of cis-regulatory sequences and changes in the transcription regulators themselves, the binding connections between regulators and their target genes rapidly change over evolutionary time and constitute a major source of biological novelty. This review covers recent work, carried out in a wide range of species, that addresses the overall extent of these evolutionary changes, their consequences, and some of the molecular mechanisms that lie behind them.

Conservation and Divergence of p53 Oscillation Dynamics across Species

The tumor-suppressing transcription factor p53 is highly conserved at the protein level and plays a key role in the DNA damage response. One important aspect of p53 regulation is its dynamics in response to DNA damage, which include oscillations. Here, we observe that, while the qualitative oscillatory nature of p53 dynamics is conserved across cell lines derived from human, monkey, dog, mouse, and rat, the oscillation period is variable. Specifically, rodent cells exhibit rapid p53 oscillations, whereas dog, monkey, and human cells show slower oscillations. Computational modeling and experiments identify stronger negative feedback between p53 and MDM2 as the driver of faster oscillations in rodents, suggesting that the period of oscillation is a network-level property. In total, our study shows that despite highly conserved signaling, the quantitative features of p53 oscillations can diverge across evolution. We caution that strong amino acid conservation of proteins and transcriptional network similarity do not necessarily imply conservation of time dynamics.

An Evolutionary Perspective on Yeast Mating-Type Switching

Cell differentiation in yeast species is controlled by a reversible, programmed DNA-rearrangement process called mating-type switching. Switching is achieved by two functionally similar but structurally distinct processes in the budding yeast Saccharomyces cerevisiae and the fission yeast Schizosaccharomyces pombe. In both species, haploid cells possess one active and two silent copies of the mating-type locus (a three-cassette structure), the active locus is cleaved, and synthesis-dependent strand annealing is used to replace it with a copy of a silent locus encoding the opposite mating-type information. Each species has its own set of components responsible for regulating these processes. In this review, we summarize knowledge about the function and evolution of mating-type switching components in these species, including mechanisms of heterochromatin formation, MAT locus cleavage, donor bias, lineage tracking, and environmental regulation of switching. We compare switching in these well-studied species to others such as Kluyveromyces lactis and the methylotrophic yeasts Ogataea polymorpha and Komagataella phaffii. We focus on some key questions: Which cells switch mating type? What molecular apparatus is required for switching? Where did it come from? And what is the evolutionary purpose of switching?

An Essential Regulatory System Originating from Polygenic Transcriptional Rewiring of PhoP-PhoQ of Xanthomonas campestris

How essential, regulatory genes originate and evolve is intriguing because mutations of these genes not only lead to lethality in organisms, but also have pleiotropic effects since they control the expression of multiple downstream genes. Therefore, the evolution of essential, regulatory genes is not only determined by genetic variations of their own sequences, but also by the biological function of downstream genes and molecular mechanisms of regulation. To understand the origin of essential, regulatory genes, experimental dissection of the complete regulatory cascade is needed. Here, we provide genetic evidences to reveal that PhoP-PhoQ is an essential two-component signal transduction system in the gram-negative bacterium Xanthomonas campestris, but that its orthologs in other bacteria belonging to Proteobacteria are nonessential. Mutational, biochemical, and chromatin immunoprecipitation together with high-throughput sequencing analyses revealed that phoP and phoQ of X. campestris and its close relative Pseudomonas aeruginosa are replaceable, and that the consensus binding motifs of the transcription factor PhoP are also highly conserved. PhoP _Xcc in X. campestris regulates the transcription of a number of essential, structural genes by directly binding to cis-regulatory elements (CREs); however, these CREs are lacking in the orthologous essential, structural genes in P. aeruginosa, and thus the regulatory relationships between PhoP _Pae and these downstream essential genes are disassociated. Our findings suggested that the recruitment of regulatory proteins by critical structural genes via transcription factor-CRE rewiring is a driving force in the origin and functional divergence of essential, regulatory genes.

Polymorphisms in the yeast galactose sensor underlie a natural continuum of nutrient-decision phenotypes

In nature, microbes often need to "decide" which of several available nutrients to utilize, a choice that depends on a cell's inherent preference and external nutrient levels. While natural environments can have mixtures of different nutrients, phenotypic variation in microbes' decisions of which nutrient to utilize is poorly studied. Here, we quantified differences in the concentration of glucose and galactose required to induce galactose-responsive (GAL) genes across 36 wild S. cerevisiae strains. Using bulk segregant analysis, we found that a locus containing the galactose sensor GAL3 was associated with differences in GAL signaling in eight different crosses. Using allele replacements, we confirmed that GAL3 is the major driver of GAL induction variation, and that GAL3 allelic variation alone can explain as much as 90% of the variation in GAL induction in a cross. The GAL3 variants we found modulate the diauxic lag, a selectable trait. These results suggest that ecological constraints on the galactose pathway may have led to variation in a single protein, allowing cells to quantitatively tune their response to nutrient changes in the environment.

A spectrum of modularity in multi-functional gene circuits

A major challenge in systems biology is to understand the relationship between a circuit's structure and its function, but how is this relationship affected if the circuit must perform multiple distinct functions within the same organism? In particular, to what extent do multi‐functional circuits contain modules which reflect the different functions? Here, we computationally survey a range of bi‐functional circuits which show no simple structural modularity: They can switch between two qualitatively distinct functions, while both functions depend on all genes of the circuit. Our analysis reveals two distinct classes: hybrid circuits which overlay two simpler mono‐functional sub‐circuits within their circuitry, and emergent circuits, which do not. In this second class, the bi‐functionality emerges from more complex designs which are not fully decomposable into distinct modules and are consequently less intuitive to predict or understand. These non‐intuitive emergent circuits are just as robust as their hybrid counterparts, and we therefore suggest that the common bias toward studying modular systems may hinder our understanding of real biological circuits.

Gene regulatory network plasticity predates a switch in function of a conserved transcription regulator

The rewiring of gene regulatory networks can generate phenotypic novelty. It remains an open question, however, how the large number of connections needed to form a novel network arise over evolutionary time. Here, we address this question using the network controlled by the fungal transcription regulator Ndt80. This conserved protein has undergone a dramatic switch in function-from an ancestral role regulating sporulation to a derived role regulating biofilm formation. This switch in function corresponded to a large-scale rewiring of the genes regulated by Ndt80. However, we demonstrate that the Ndt80-target gene connections were undergoing extensive rewiring prior to the switch in Ndt80's regulatory function. We propose that extensive drift in the Ndt80 regulon allowed for the exploration of alternative network structures without a loss of ancestral function, thereby facilitating the formation of a network with a new function.

Widespread Historical Contingency in Influenza Viruses

In systems biology and genomics, epistasis characterizes the impact that a substitution at a particular location in a genome can have on a substitution at another location. This phenomenon is often implicated in the evolution of drug resistance or to explain why particular "disease-causing" mutations do not have the same outcome in all individuals. Hence, uncovering these mutations and their locations in a genome is a central question in biology. However, epistasis is notoriously difficult to uncover, especially in fast-evolving organisms. Here, we present a novel statistical approach that replies on a model developed in ecology and that we adapt to analyze genetic data in fast-evolving systems such as the influenza A virus. We validate the approach using a two-pronged strategy: extensive simulations demonstrate a low-to-moderate sensitivity with excellent specificity and precision, while analyses of experimentally validated data recover known interactions, including in a eukaryotic system. We further evaluate the ability of our approach to detect correlated evolution during antigenic shifts or at the emergence of drug resistance. We show that in all cases, correlated evolution is prevalent in influenza A viruses, involving many pairs of sites linked together in chains; a hallmark of historical contingency. Strikingly, interacting sites are separated by large physical distances, which entails either long-range conformational changes or functional tradeoffs, for which we find support with the emergence of drug resistance. Our work paves a new way for the unbiased detection of epistasis in a wide range of organisms by performing whole-genome scans.

Global Identification of Biofilm-Specific Proteolysis in Candida albicans

Candida albicans is a fungal species that is part of the normal human microbiota and also an opportunistic pathogen capable of causing mucosal and systemic infections. C. albicans cells proliferate in a planktonic (suspension) state, but they also form biofilms, organized and tightly packed communities of cells attached to a solid surface. Biofilms colonize many niches of the human body and persist on implanted medical devices, where they are a major source of new C. albicans infections. Here, we used an unbiased and global substrate-profiling approach to discover proteolytic activities produced specifically by C. albicans biofilms, compared to planktonic cells, with the goal of identifying potential biofilm-specific diagnostic markers and targets for therapeutic intervention. This activity-based profiling approach, coupled with proteomics, identified Sap5 (Candidapepsin-5) and Sap6 (Candidapepsin-6) as major biofilm-specific proteases secreted by C. albicans. Fluorogenic peptide substrates with selectivity for Sap5 or Sap6 confirmed that their activities are highly upregulated in C. albicans biofilms; we also show that these activities are upregulated in other Candida clade pathogens. Deletion of the SAP5 and SAP6 genes in C. albicans compromised biofilm development in vitro in standard biofilm assays and in vivo in a rat central venous catheter biofilm model. This work establishes secreted proteolysis as a promising enzymatic marker and potential therapeutic target for Candida biofilm formation.

Biofilm formation by the opportunistic fungal pathogen C. albicans is a major cause of life-threatening infections. This work provides a global characterization of secreted proteolytic activity produced specifically by C. albicans biofilms. We identify activity from the proteases Sap5 and Sap6 as highly upregulated during C. albicans biofilm formation and develop Sap-cleavable fluorogenic substrates that enable the detection of biofilms from C. albicans and also from additional pathogenic Candida species. Furthermore, SAP5 and SAP6 deletions confirm that both proteases are required for proper biofilm development in vitro and in vivo. We propose that secreted proteolysis is a promising marker for the diagnosis and potential therapeutic targeting of Candida biofilm-associated infections.

Novelty and Innovation in the History of Life

The history of life as documented by the fossil record encompasses evolutionary diversifications at scales ranging from the Ediacaran-Cambrian explosion of animal life and the invasion of land by vascular plants, insects and vertebrates to the diversification of flowering plants over the past 100 million years and the radiation of horses. Morphological novelty and innovation has been a recurrent theme. The architects of the modern synthesis of evolutionary theory made three claims about evolutionary novelty and innovation: first, that all diversifications in the history of life represent adaptive radiations; second, that adaptive radiations are driven principally by ecological opportunity rather than by the supply of new morphological novelties, thus the primary questions about novelty and innovation focus on their ecological and evolutionary success; and third, that the rate of morphological divergence between taxa was more rapid early in the history of a clade but slowed over time as ecological opportunities declined. These claims have strongly influenced subsequent generations of evolutionary biologists, yet over the past two decades each has been challenged by data from the fossil record, by the results of comparative phylogenetic analyses and through insights from evolutionary developmental biology. Consequently a broader view of novelty and innovation is required. An outstanding issue for future work is identifying the circumstances associated with different styles of diversification and whether their frequency has changed through the history of life.

Disruption of promoter memory by synthesis of a long noncoding RNA

The yeast HO endonuclease is expressed in late G1 in haploid mother cells to initiate mating-type interconversion. Cells can be arrested in G1 by nutrient deprivation or by pheromone exposure, but cells that resume cycling after nutrient deprivation or cyclin-dependent kinase (CDK) inactivation express HO in the first cell cycle, whereas HO is not expressed until the second cycle after release from pheromone arrest. Here, we show that transcription of a long noncoding RNA (lncRNA) mediates this differential response. The SBF and Mediator factors remain bound to the inactive promoter during arrest due to CDK inactivation, and these bound factors allow the cell to remember a transcriptional decision made before arrest. If the presence of mating pheromone indicates that this decision is no longer appropriate, a lncRNA originating at -2700 upstream of the HO gene is induced, and the transcription machinery displaces promoter-bound SBF, preventing HO transcription in the subsequent cell cycle. Further, we find that the displaced SBF is blocked from rebinding due to incorporation of its recognition sites within nucleosomes. Expressing the pHO-lncRNA in trans is ineffective, indicating that transcription in cis is required. Factor displacement during lncRNA transcription could be a general mechanism for regulating memory of previous events at promoters.

Epistasis in protein evolution

The structure, function, and evolution of proteins depend on physical and genetic interactions among amino acids. Recent studies have used new strategies to explore the prevalence, biochemical mechanisms, and evolutionary implications of these interactions-called epistasis-within proteins. Here we describe an emerging picture of pervasive epistasis in which the physical and biological effects of mutations change over the course of evolution in a lineage-specific fashion. Epistasis can restrict the trajectories available to an evolving protein or open new paths to sequences and functions that would otherwise have been inaccessible. We describe two broad classes of epistatic interactions, which arise from different physical mechanisms and have different effects on evolutionary processes. Specific epistasis-in which one mutation influences the phenotypic effect of few other mutations-is caused by direct and indirect physical interactions between mutations, which nonadditively change the protein's physical properties, such as conformation, stability, or affinity for ligands. In contrast, nonspecific epistasis describes mutations that modify the effect of many others; these typically behave additively with respect to the physical properties of a protein but exhibit epistasis because of a nonlinear relationship between the physical properties and their biological effects, such as function or fitness. Both types of interaction are rampant, but specific epistasis has stronger effects on the rate and outcomes of evolution, because it imposes stricter constraints and modulates evolutionary potential more dramatically; it therefore makes evolution more contingent on low-probability historical events and leaves stronger marks on the sequences, structures, and functions of protein families.

Differences in environmental stress response among yeasts is consistent with species-specific lifestyles

Defining how organisms respond to environmental change has always been an important step toward understanding their adaptive capacity and physiology. Variation in transcription during stress has been widely described in model species, especially in the yeast Saccharomyces cerevisiae, which helped to shape general rules regarding how cells cope with environmental constraints, as well as to decipher the functions of many genes. Comparison of the environmental stress response (ESR) across species is essential to obtaining better insight into the common and species-specific features of stress defense. In this context, we explored the transcriptional landscape of the yeast Lachancea kluyveri (formerly Saccharomyces kluyveri) in response to diverse stresses, using RNA sequencing. We investigated variation in gene expression and observed a link between genetic plasticity and environmental sensitivity. We identified the ESR genes in this species and compared them to those already found in S. cerevisiae We observed common features between the two species, as well as divergence in the regulatory networks involved. Of interest, some changes were related to differences in species lifestyle. Thus we were able to decipher how adaptation to stress has evolved among different yeast species. Finally, by analyzing patterns of coexpression, we were able to propose potential biological functions for 42% of genes and also annotate 301 genes for which no function could be assigned by homology. This large data set allowed for the characterization of the evolution of gene regulation and provides an efficient tool for assessing gene function.