How duplicated transcription regulators can diversify to govern the expression of nonoverlapping sets of genes

Adaptation of Candida albicans to specific host environments by gain-of-function mutations in transcription factors

The yeast Candida albicans is usually a harmless member of the normal microbiota in healthy persons but is also a major fungal pathogen that can colonize and infect almost every human tissue. A successful adaptation to environmental changes encountered in different host niches requires an appropriate regulation of gene expression. The zinc cluster transcription factors are the largest family of transcriptional regulators in C. albicans and are involved in the control of virtually all aspects of its biology. Under certain circumstances, mutations in these transcription factors that alter their activity and the expression of their target genes confer a selective advantage, which results in the emergence of phenotypically altered variants that are better adapted to new environmental challenges. This review describes how gain-of-function mutations in different zinc cluster transcription factors enable C. albicans to overcome antifungal therapy and to successfully establish itself in specific host niches.

Gene networks governing the response of a calcareous sponge to future ocean conditions reveal lineage-specific XBP1 regulation of the unfolded protein response

Marine sponges are predicted to be winners in the future ocean due to their exemplary adaptive capacity. However, while many sponge groups exhibit tolerance to a wide range of environmental insults, calcifying sponges may be more susceptible to thermo-acidic stress. To describe the gene regulatory networks that govern the stress response of the calcareous sponge, Leucetta chagosensis (class Calcarea, order Clathrinida), individuals were subjected to warming and acidification conditions based on the climate models for 2100. Transcriptome analysis and gene co-expression network reconstruction revealed that the unfolded protein response (UPR) was activated under thermo-acidic stress. Among the upregulated genes were two lineage-specific homologs of X-box binding protein 1 (XBP1), a transcription factor that activates the UPR. Alternative dimerization between these XBP1 gene products suggests a clathrinid-specific mechanism to reversibly sequester the transcription factor into an inactive form, enabling the rapid regulation of pathways linked to the UPR in clathrinid calcareous sponges. Our findings support the idea that transcription factor duplication events may refine evolutionarily conserved molecular pathways and contribute to ecological success.

Gene loss and cis-regulatory novelty shaped core histone gene evolution in the apiculate yeast Hanseniaspora uvarum

Core histone genes display a remarkable diversity of cis-regulatory mechanisms despite their protein sequence conservation. However, the dynamics and significance of this regulatory turnover are not well understood. Here, we describe the evolutionary history of core histone gene regulation across 400 million years in budding yeasts. We find that canonical mode of core histone regulation-mediated by the trans-regulator Spt10-is ancient, likely emerging between 320 and 380 million years ago and is fixed in the majority of extant species. Unexpectedly, we uncovered the emergence of a novel core histone regulatory mode in the Hanseniaspora genus, from its fast-evolving lineage, which coincided with the loss of 1 copy of its paralogous core histone genes. We show that the ancestral Spt10 histone regulatory mode was replaced, via cis-regulatory changes in the histone control regions, by a derived Mcm1 histone regulatory mode and that this rewiring event occurred with no changes to the trans-regulator, Mcm1, itself. Finally, we studied the growth dynamics of the cell cycle and histone synthesis in genetically modified Hanseniaspora uvarum. We find that H. uvarum divides rapidly, with most cells completing a cell cycle within 60 minutes. Interestingly, we observed that the regulatory coupling between histone and DNA synthesis was lost in H. uvarum. Our results demonstrate that core histone gene regulation was fixed anciently in budding yeasts, however it has greatly diverged in the Hanseniaspora fast-evolving lineage.

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

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

Evolution of binding preferences among whole-genome duplicated transcription factors

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

GAL Regulon in the Yeast S. cerevisiae is Highly Evolvable via Acquisition in the Coding Regions of the Regulatory Elements of the Network

GAL network in the yeast S. cerevisiae is one of the most well-characterized regulatory network. Expression of GAL genes is contingent on exposure to galactose, and an appropriate combination of the alleles of the regulatory genes GAL3, GAL1, GAL80, and GAL4. The presence of multiple regulators in the GAL network makes it unique, as compared to the many sugar utilization networks studied in bacteria. For example, utilization of lactose is controlled by a single regulator LacI, in E. coli's lac operon. Moreover, recent work has demonstrated that multiple alleles of these regulatory proteins are present in yeast isolated from ecological niches. In this work, we develop a mathematical model, and demonstrate via deterministic and stochastic runs of the model, that behavior/gene expression patterns of the cells (at a population level, and at a single-cell resolution) can be modulated by altering the binding affinities between the regulatory proteins. This adaptability is likely the key to explaining the multiple GAL regulatory alleles discovered in ecological isolates in recent years.

Evolution of a novel cell type in Dictyostelia required gene duplication of a cudA-like transcription factor

The evolution of novel cell types has been proposed to result from duplication of gene regulatory networks, but proven examples are rare. In addition to stalk cells and spores that make up the fruiting bodies of three major groups of Dictyostelia, those in group 4 additionally evolved basal disc and cup cells that respectively anchor the stalk to the substratum and the spore mass to the stalk. We noted a putative group-4-specific duplication of a cudA-like transcription factor (TF) in a comparative analysis of group-representative genomes. Using increased taxon sampling, we here confirmed that this TF, cdl1, duplicated into cdl1a and cdl1b in the common ancestor to group 4. cdl1a, but not cdl1b, showed signatures of positive selection, indicative of functional innovation. Deletion of cdl1a in Dictyostelium discoideum resulted in fruiting bodies with sagging spore heads that lacked the supporting cup cells and expression of cup-specific genes. Deletion of cdl1b resulted in thinner fruiting body stalks, while a cdl1b-cdl1a- double knockout showed more severe stalk defects, suggesting an ancestral role of cdl1 in stalk formation. This was confirmed in a closely related non-group 4 species, Polysphondylium violaceum, where cdl1 knockout caused defective stalk formation. These data indicate that the group-specific duplication of cdl1 and subsequent diversification of cdl1a played a pivotal role in the evolution of a novel somatic cell type in group 4 Dictyostelia.

Three topological features of regulatory networks control life-essential and specialized subsystems

Gene regulatory networks (GRNs) play key roles in development, phenotype plasticity, and evolution. Although graph theory has been used to explore GRNs, associations amongst topological features, transcription factors (TFs), and systems essentiality are poorly understood. Here we sought the relationship amongst the main GRN topological features that influence the control of essential and specific subsystems. We found that the K_nn, page rank, and degree are the most relevant GRN features: the ones are conserved along the evolution and are also relevant in pluripotent cells. Interestingly, life-essential subsystems are governed mainly by TFs with intermediary K_nn and high page rank or degree, whereas specialized subsystems are mainly regulated by TFs with low K_nn. Hence, we suggest that the high probability of TFs be toured by a random signal, and the high probability of the signal propagation to target genes ensures the life-essential subsystems' robustness. Gene/genome duplication is the main evolutionary process to rise K_nn as the most relevant feature. Herein, we shed light on unexplored topological GRN features to assess how they are related to subsystems and how the duplications shaped the regulatory systems along the evolution. The classification model generated can be found here: https://github.com/ivanrwolf/NoC/ .

Non-canonical LexA proteins regulate the SOS response in the Bacteroidetes

Lesions to DNA compromise chromosome integrity, posing a direct threat to cell survival. The bacterial SOS response is a widespread transcriptional regulatory mechanism to address DNA damage. This response is coordinated by the LexA transcriptional repressor, which controls genes involved in DNA repair, mutagenesis and cell-cycle control. To date, the SOS response has been characterized in most major bacterial groups, with the notable exception of the Bacteroidetes. No LexA homologs had been identified in this large, diverse and ecologically important phylum, suggesting that it lacked an inducible mechanism to address DNA damage. Here, we report the identification of a novel family of transcriptional repressors in the Bacteroidetes that orchestrate a canonical response to DNA damage in this phylum. These proteins belong to the S24 peptidase family, but are structurally different from LexA. Their N-terminal domain is most closely related to CI-type bacteriophage repressors, suggesting that they may have originated from phage lytic phase repressors. Given their role as SOS regulators, however, we propose to designate them as non-canonical LexA proteins. The identification of a new class of repressors orchestrating the SOS response illuminates long-standing questions regarding the origin and plasticity of this transcriptional network.

An insight on the impact of teleost whole genome duplication on the regulation of the molecular networks controlling skeletal muscle growth

Fish muscle growth is a complex process regulated by multiple pathways, resulting on the net accumulation of proteins and the activation of myogenic progenitor cells. Around 350–320 million years ago, teleost fish went through a specific whole genome duplication (WGD) that expanded the existent gene repertoire. Duplicated genes can be retained by different molecular mechanisms such as subfunctionalization, neofunctionalization or redundancy, each one with different functional implications. While the great majority of ohnolog genes have been identified in the teleost genomes, the effect of gene duplication in the fish physiology is still not well characterized. In the present study we studied the effect of WGD on the transcription of the duplicated components controlling muscle growth. We compared the expression of lineage-specific ohnologs related to myogenesis and protein balance in the fast-skeletal muscle of pacus (Piaractus mesopotamicus—Ostariophysi) and Nile tilapias (Oreochromis niloticus—Acanthopterygii) fasted for 4 days and refed for 3 days. We studied the expression of 20 ohnologs and found that in the great majority of cases, duplicated genes had similar expression profiles in response to fasting and refeeding, indicating that their functions during growth have been conserved during the period after the WGD. Our results suggest that redundancy might play a more important role in the retention of ohnologs of regulatory pathways than initially thought. Also, comparison to non-duplicated orthologs showed that it might not be uncommon for the duplicated genes to gain or loss new regulatory elements simultaneously. Overall, several of duplicated ohnologs have similar transcription profiles in response to pro-growth signals suggesting that evolution tends to conserve ohnolog regulation during muscle development and that in the majority of ohnologs related to muscle growth their functions might be very similar.

The physical basis and practical consequences of biological promiscuity

Proteins interact with metabolites, nucleic acids, and other proteins to orchestrate the myriad catalytic, structural and regulatory functions that support life from the simplest microbes to the most complex multicellular organisms. These molecular interactions are often exquisitely specific, but never perfectly so. Adventitious "promiscuous" interactions are ubiquitous due to the thousands of macromolecules and small molecules crowded together in cells. Such interactions may perturb protein function at the molecular level, but as long as they do not compromise organismal fitness, they will not be removed by natural selection. Although promiscuous interactions are physiologically irrelevant, they are important because they can provide a vast reservoir of potential functions that can provide the starting point for evolution of new functions, both in nature and in the laboratory.

Resolving noise-control conflict by gene duplication

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

Diversification of DNA binding specificities enabled SREBP transcription regulators to expand the repertoire of cellular functions that they govern in fungi

The Sterol Regulatory Element Binding Proteins (SREBPs) are basic-helix-loop-helix transcription regulators that control the expression of sterol biosynthesis genes in higher eukaryotes and some fungi. Surprisingly, SREBPs do not regulate sterol biosynthesis in the ascomycete yeasts (Saccharomycotina) as this role was handed off to an unrelated transcription regulator in this clade. The SREBPs, nonetheless, expanded in fungi such as the ascomycete yeasts Candida spp., raising questions about their role and evolution in these organisms. Here we report that the fungal SREBPs diversified their DNA binding preferences concomitantly with an expansion in function. We establish that several branches of fungal SREBPs preferentially bind non-palindromic DNA sequences, in contrast to the palindromic DNA motifs recognized by most basic-helix-loop-helix proteins (including SREBPs) in higher eukaryotes. Reconstruction and biochemical characterization of the likely ancestor protein suggest that an intrinsic DNA binding promiscuity in the family was resolved by alternative mechanisms in different branches of fungal SREBPs. Furthermore, we show that two SREBPs in the human commensal yeast Candida albicans drive a transcriptional cascade that inhibits a morphological switch under anaerobic conditions. Preventing this morphological transition enhances C. albicans colonization of the mammalian intestine, the fungus' natural niche. Thus, our results illustrate how diversification in DNA binding preferences enabled the functional expansion of a family of eukaryotic transcription regulators.

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.

Asymmetry in indegree and outdegree distributions of gene regulatory networks arising from dynamical robustness

Although outdegree distributions of gene regulatory networks have scale-free characteristics similar to other biological networks, indegree distributions have single-scale characteristics with significantly lower variance than that of outdegree distributions. In this study, we mathematically explain that such asymmetric characteristics arise from dynamical robustness, which is the property of maintaining an equilibrium state of gene expressions against inevitable perturbations to the networks, such as gene dysfunction and mutation of promoters. We reveal that the expression of a single gene is robust to a perturbation for a large number of inputs and a small number of outputs. Applying these results to the networks, we also show that an equilibrium state of the networks is robust if the variance of the indegree distribution is low (i.e., single-scale characteristics) and that of the outdegree distribution is high (i.e., scale-free characteristics). These asymmetric characteristics are conserved across a wide range of species, from bacteria to humans.

Comparative Transcriptomics Highlights New Features of the Iron Starvation Response in the Human Pathogen Candida glabrata

In this work, we used comparative transcriptomics to identify regulatory outliers (ROs) in the human pathogen Candida glabrata. ROs are genes that have very different expression patterns compared to their orthologs in other species. From comparative transcriptome analyses of the response of eight yeast species to toxic doses of selenite, a pleiotropic stress inducer, we identified 38 ROs in C. glabrata. Using transcriptome analyses of C. glabrata response to five different stresses, we pointed out five ROs which were more particularly responsive to iron starvation, a process which is very important for C. glabrata virulence. Global chromatin Immunoprecipitation and gene profiling analyses showed that four of these genes are actually new targets of the iron starvation responsive Aft2 transcription factor in C. glabrata. Two of them (HBS1 and DOM34b) are required for C. glabrata optimal growth in iron limited conditions. In S. cerevisiae, the orthologs of these two genes are involved in ribosome rescue by the NO GO decay (NGD) pathway. Hence, our results suggest a specific contribution of NGD co-factors to the C. glabrata adaptation to iron starvation.

Evolution of new regulatory functions on biophysically realistic fitness landscapes

Gene expression is controlled by networks of regulatory proteins that interact specifically with external signals and DNA regulatory sequences. These interactions force the network components to co-evolve so as to continually maintain function. Yet, existing models of evolution mostly focus on isolated genetic elements. In contrast, we study the essential process by which regulatory networks grow: the duplication and subsequent specialization of network components. We synthesize a biophysical model of molecular interactions with the evolutionary framework to find the conditions and pathways by which new regulatory functions emerge. We show that specialization of new network components is usually slow, but can be drastically accelerated in the presence of regulatory crosstalk and mutations that promote promiscuous interactions between network components.Gene networks evolve by transcription factor (TF) duplication and divergence of their binding site specificities, but little is known about the global constraints at play. Here, the authors study the coevolution of TFs and binding sites using a biophysical-evolutionary approach, and show that the emerging complex fitness landscapes strongly influence regulatory evolution with a role for crosstalk.

Modelling the evolution of transcription factor binding preferences in complex eukaryotes

Transcription factors (TFs) exert their regulatory action by binding to DNA with specific sequence preferences. However, different TFs can partially share their binding sequences due to their common evolutionary origin. This "redundancy" of binding defines a way of organizing TFs in "motif families" by grouping TFs with similar binding preferences. Since these ultimately define the TF target genes, the motif family organization entails information about the structure of transcriptional regulation as it has been shaped by evolution. Focusing on the human TF repertoire, we show that a one-parameter evolutionary model of the Birth-Death-Innovation type can explain the TF empirical repartition in motif families, and allows to highlight the relevant evolutionary forces at the origin of this organization. Moreover, the model allows to pinpoint few deviations from the neutral scenario it assumes: three over-expanded families (including HOX and FOX genes), a set of "singleton" TFs for which duplication seems to be selected against, and a higher-than-average rate of diversification of the binding preferences of TFs with a Zinc Finger DNA binding domain. Finally, a comparison of the TF motif family organization in different eukaryotic species suggests an increase of redundancy of binding with organism complexity.

Inference and Evolutionary Analysis of Genome-Scale Regulatory Networks in Large Phylogenies

Changes in transcriptional regulatory networks can significantly contribute to species evolution and adaptation. However, identification of genome-scale regulatory networks is an open challenge, especially in non-model organisms. Here, we introduce multi-species regulatory network learning (MRTLE), a computational approach that uses phylogenetic structure, sequence-specific motifs, and transcriptomic data, to infer the regulatory networks in different species. Using simulated data from known networks and transcriptomic data from six divergent yeasts, we demonstrate that MRTLE predicts networks with greater accuracy than existing methods because it incorporates phylogenetic information. We used MRTLE to infer the structure of the transcriptional networks that control the osmotic stress responses of divergent, non-model yeast species and then validated our predictions experimentally. Interrogating these networks reveals that gene duplication promotes network divergence across evolution. Taken together, our approach facilitates study of regulatory network evolutionary dynamics across multiple poorly studied species.

Ongoing resolution of duplicate gene functions shapes the diversification of a metabolic network

The evolutionary mechanisms leading to duplicate gene retention are well understood, but the long-term impacts of paralog differentiation on the regulation of metabolism remain underappreciated. Here we experimentally dissect the functions of two pairs of ancient paralogs of the GALactose sugar utilization network in two yeast species. We show that the Saccharomyces uvarum network is more active, even as over-induction is prevented by a second co-repressor that the model yeast Saccharomyces cerevisiae lacks. Surprisingly, removal of this repression system leads to a strong growth arrest, likely due to overly rapid galactose catabolism and metabolic overload. Alternative sugars, such as fructose, circumvent metabolic control systems and exacerbate this phenotype. We further show that S. cerevisiae experiences homologous metabolic constraints that are subtler due to how the paralogs have diversified. These results show how the functional differentiation of paralogs continues to shape regulatory network architectures and metabolic strategies long after initial preservation.

Systematic Genetic Screen for Transcriptional Regulators of the Candida albicans White-Opaque Switch

The human fungal pathogen Candida albicans can reversibly switch between two cell types named "white" and "opaque," each of which is stable through many cell divisions. These two cell types differ in their ability to mate, their metabolic preferences and their interactions with the mammalian innate immune system. A highly interconnected network of eight transcriptional regulators has been shown to control switching between these two cell types. To identify additional regulators of the switch, we systematically and quantitatively measured white-opaque switching rates of 196 strains, each deleted for a specific transcriptional regulator. We identified 19 new regulators with at least a 10-fold effect on switching rates and an additional 14 new regulators with more subtle effects. To investigate how these regulators affect switching rates, we examined several criteria, including the binding of the eight known regulators of switching to the control region of each new regulatory gene, differential expression of the newly found genes between cell types, and the growth rate of each mutant strain. This study highlights the complexity of the transcriptional network that regulates the white-opaque switch and the extent to which switching is linked to a variety of metabolic processes, including respiration and carbon utilization. In addition to revealing specific insights, the information reported here provides a foundation to understand the highly complex coupling of white-opaque switching to cellular physiology.

Reconfiguration of Transcriptional Control of Lysine Biosynthesis in Candida albicans Involves a Central Role for the Gcn4 Transcriptional Activator

Microbes evolve rapidly so as to reconfigure their gene expression to adapt to the metabolic demands in diverse environmental niches. Here, we explored how conditions of nutrient deprivation regulate lysine biosynthesis in the human fungal pathogen Candida albicans. We show that although both Saccharomyces cerevisiae and C. albicans respond to lysine deprivation by transcriptional upregulation of lysine biosynthesis, the regulatory factors required for this control have been reconfigured in these species. We found that Gcn4 is an essential and direct transcriptional regulator of the expression of lysine biosynthetic genes under lysine starvation conditions in C. albicans. Our results therefore suggest that the regulation of the lysine biosynthetic pathway in Candida clade genomes involves gain of function by the master transcriptional regulator Gcn4, coincident with the neofunctionalization of the S. cerevisiae pathway-specific regulator Lys14.

Evolution of transcriptional control is essential for organisms to cope with diversification into a spectrum of environments, including environments with limited nutrients. Lysine biosynthesis in fungi occurs in eight enzymatic steps. In Saccharomyces cerevisiae, amino acid starvation elicits the induction of LYS gene expression, mediated by the master regulator Gcn4 and the pathway-specific transcriptional regulator Lys14. Here, we have shown that the activation of LYS gene expression in the human fungal pathogen Candida albicans is predominantly controlled by Gcn4 under amino acid starvation conditions. Multiple lines of study showed that the four C. albicans LYS14-like genes have no role in the regulation of lysine biosynthesis. Whereas Gcn4 is dispensable for the growth of S. cerevisiae under lysine deprivation conditions, it is an essential regulator required for the growth of C. albicans under these conditions, as gcn4 deletion caused lysine auxotrophy. Gcn4 is required for the induction of increased LYS2 and LYS9 mRNA but not for the induction of increased LYS4 mRNA. Under lysine or isoleucine-valine deprivation conditions, Gcn4 recruitment to LYS2 and LYS9 promoters was induced in C. albicans. Indeed, in contrast to the S. cerevisiae LYS gene promoters, all LYS gene promoters in C. albicans harbored a Gcn4 binding site but not all harbored the S. cerevisiae Lys14 binding site, indicating the evolutionary divergence of cis-regulatory motifs. Thus, the transcriptional rewiring of the lysine biosynthetic pathway in C. albicans involves not only neofunctionalization of the four LYS14-like genes but the attendant strengthening of control by Gcn4, indicating a coordinated response with a much broader scope for control of amino acid biosynthesis in this human pathogen.

IMPORTANCE Microbes evolve rapidly so as to reconfigure their gene expression to adapt to the metabolic demands in diverse environmental niches. Here, we explored how conditions of nutrient deprivation regulate lysine biosynthesis in the human fungal pathogen Candida albicans. We show that although both Saccharomyces cerevisiae and C. albicans respond to lysine deprivation by transcriptional upregulation of lysine biosynthesis, the regulatory factors required for this control have been reconfigured in these species. We found that Gcn4 is an essential and direct transcriptional regulator of the expression of lysine biosynthetic genes under lysine starvation conditions in C. albicans. Our results therefore suggest that the regulation of the lysine biosynthetic pathway in Candida clade genomes involves gain of function by the master transcriptional regulator Gcn4, coincident with the neofunctionalization of the S. cerevisiae pathway-specific regulator Lys14.

How Transcription Networks Evolve and Produce Biological Novelty

The rewiring of gene regulatory networks over evolutionary timescales produces changes in the patterns of gene expression and is a major source of diversity among species. Yet the molecular mechanisms underlying evolutionary rewiring are only beginning to be understood. Here, we discuss recent analyses in ascomycete yeasts that have revealed several general principles of network rewiring. Specifically, we discuss how transcription networks can maintain a functional output despite changes in mechanism, how specific types of constraints alter available evolutionary trajectories, and how regulatory rewiring can ultimately lead to phenotypic novelty. We also argue that the structure and "logic" of extant gene regulatory networks can largely be accounted for by constraints that shape their evolutionary trajectories.

Genomics and the making of yeast biodiversity

Yeasts are unicellular fungi that do not form fruiting bodies. Although the yeast lifestyle has evolved multiple times, most known species belong to the subphylum Saccharomycotina (syn. Hemiascomycota, hereafter yeasts). This diverse group includes the premier eukaryotic model system, Saccharomyces cerevisiae; the common human commensal and opportunistic pathogen, Candida albicans; and over 1000 other known species (with more continuing to be discovered). Yeasts are found in every biome and continent and are more genetically diverse than angiosperms or chordates. Ease of culture, simple life cycles, and small genomes (∼10-20Mbp) have made yeasts exceptional models for molecular genetics, biotechnology, and evolutionary genomics. Here we discuss recent developments in understanding the genomic underpinnings of the making of yeast biodiversity, comparing and contrasting natural and human-associated evolutionary processes. Only a tiny fraction of yeast biodiversity and metabolic capabilities has been tapped by industry and science. Expanding the taxonomic breadth of deep genomic investigations will further illuminate how genome function evolves to encode their diverse metabolisms and ecologies.

Evolutionary Conservation and Diversification of Puf RNA Binding Proteins and Their mRNA Targets

Reprogramming of a gene’s expression pattern by acquisition and loss of sequences recognized by specific regulatory RNA binding proteins may be a major mechanism in the evolution of biological regulatory programs. We identified that RNA targets of Puf3 orthologs have been conserved over 100–500 million years of evolution in five eukaryotic lineages. Focusing on Puf proteins and their targets across 80 fungi, we constructed a parsimonious model for their evolutionary history. This model entails extensive and coordinated changes in the Puf targets as well as changes in the number of Puf genes and alterations of RNA binding specificity including that: 1) Binding of Puf3 to more than 200 RNAs whose protein products are predominantly involved in the production and organization of mitochondrial complexes predates the origin of budding yeasts and filamentous fungi and was maintained for 500 million years, throughout the evolution of budding yeast. 2) In filamentous fungi, remarkably, more than 150 of the ancestral Puf3 targets were gained by Puf4, with one lineage maintaining both Puf3 and Puf4 as regulators and a sister lineage losing Puf3 as a regulator of these RNAs. The decrease in gene expression of these mRNAs upon deletion of Puf4 in filamentous fungi (N. crassa) in contrast to the increase upon Puf3 deletion in budding yeast (S. cerevisiae) suggests that the output of the RNA regulatory network is different with Puf4 in filamentous fungi than with Puf3 in budding yeast. 3) The coregulated Puf4 target set in filamentous fungi expanded to include mitochondrial genes involved in the tricarboxylic acid (TCA) cycle and other nuclear-encoded RNAs with mitochondrial function not bound by Puf3 in budding yeast, observations that provide additional evidence for substantial rewiring of post-transcriptional regulation. 4) Puf3 also expanded and diversified its targets in filamentous fungi, gaining interactions with the mRNAs encoding the mitochondrial electron transport chain (ETC) complex I as well as hundreds of other mRNAs with nonmitochondrial functions. The many concerted and conserved changes in the RNA targets of Puf proteins strongly support an extensive role of RNA binding proteins in coordinating gene expression, as originally proposed by Keene. Rewiring of Puf-coordinated mRNA targets and transcriptional control of the same genes occurred at different points in evolution, suggesting that there have been distinct adaptations via RNA binding proteins and transcription factors. The changes in Puf targets and in the Puf proteins indicate an integral involvement of RNA binding proteins and their RNA targets in the adaptation, reprogramming, and function of gene expression.

A map of the evolutionary history of Puf proteins and their RNA targets shows that reprogramming of global gene expression programs via adaptive mutations that affect protein-RNA interactions is an important source of biological diversity.

We set out to trace the evolutionary history of an RNA binding protein and how its interactions with targets change over evolution. Identifying this natural history is a step toward understanding the critical differences between organisms and how gene expression programs are rewired during evolution. Using bioinformatics and experimental approaches, we broadly surveyed the evolution of binding targets of a particular family of RNA binding proteins—the Puf proteins, whose protein sequences and target RNA sequences are relatively well-characterized—across 99 eukaryotic species. We found five groups of species in which targets have been conserved for at least 100 million years and then took advantage of genome sequences from a large number of fungal species to deeply investigate the conservation and changes in Puf proteins and their RNA targets. Our analyses identified multiple and extensive reconfigurations during the natural history of fungi and suggest that RNA binding proteins and their RNA targets are profoundly involved in evolutionary reprogramming of gene expression and help define distinct programs unique to each organism. Continuing to uncover the natural history of RNA binding proteins and their interactions will provide a unique window into the gene expression programs of present day species and point to new ways to engineer gene expression programs.

Gal3 Binds Gal80 Tighter than Gal1 Indicating Adaptive Protein Changes Following Duplication

Derived from the yeast whole-genome duplication, Saccharomyces cerevisiae GAL1 and GAL3 encode the catabolic enzyme galactokinase (Gal1) and its transcriptional coinducer (Gal3), whereas the ancestral, preduplicated GAL1 gene performed both functions. Previous studies indicated that divergence was primarily driven by changes in upstream promoter elements, and changes in GAL3's coding region are assumed to be the result of drift. We show that replacement of GAL3's open-reading-frame with GAL1's results in an extended lag phase upon switching to growth on galactose with up to 2.5-fold differences in the initial cell masses. Accordingly, the binding affinity of Gal3 to Gal80 was found to be greater than 10-folds higher than that of Gal1, with both a higher association rate (ka) and lower dissociation (kd) rate. Thus, while changes in the noncoding, regulatory regions were the initial driving force for GAL3's subfunctionalization as a coinducer, adaptive changes in the protein sequence seem to have followed.

Mechanisms of mutational robustness in transcriptional regulation

Robustness is the invariance of a phenotype in the face of environmental or genetic change. The phenotypes produced by transcriptional regulatory circuits are gene expression patterns that are to some extent robust to mutations. Here we review several causes of this robustness. They include robustness of individual transcription factor binding sites, homotypic clusters of such sites, redundant enhancers, transcription factors, redundant transcription factors, and the wiring of transcriptional regulatory circuits. Such robustness can either be an adaptation by itself, a byproduct of other adaptations, or the result of biophysical principles and non-adaptive forces of genome evolution. The potential consequences of such robustness include complex regulatory network topologies that arise through neutral evolution, as well as cryptic variation, i.e., genotypic divergence without phenotypic divergence. On the longest evolutionary timescales, the robustness of transcriptional regulation has helped shape life as we know it, by facilitating evolutionary innovations that helped organisms such as flowering plants and vertebrates diversify.

Evolution of regulatory networks in Candida glabrata: learning to live with the human host

The opportunistic human fungal pathogen Candida glabrata is second only to C. albicans as the cause of Candida infections and yet is more closely related to Saccharomyces cerevisiae. Recent advances in functional genomics technologies and computational approaches to decipher regulatory networks, and the comparison of these networks among these and other Ascomycete species, have revealed both unique and shared strategies in adaptation to a human commensal/opportunistic pathogen lifestyle and antifungal drug resistance in C. glabrata. Recently, several C. glabrata sister species in the Nakeseomyces clade representing both human associated (commensal) and environmental isolates have had their genomes sequenced and analyzed. This has paved the way for comparative functional genomics studies to characterize the regulatory networks in these species to identify informative patterns of conservation and divergence linked to phenotypic evolution in the Nakaseomyces lineage.

Budding off: bringing functional genomics to Candida albicans

Candida species are the most prevalent human fungal pathogens, with Candida albicans being the most clinically relevant species. Candida albicans resides as a commensal of the human gastrointestinal tract but is a frequent cause of opportunistic mucosal and systemic infections. Investigation of C. albicans virulence has traditionally relied on candidate gene approaches, but recent advances in functional genomics have now facilitated global, unbiased studies of gene function. Such studies include comparative genomics (both between and within Candida species), analysis of total RNA expression, and regulation and delineation of protein-DNA interactions. Additionally, large collections of mutant strains have begun to aid systematic screening of clinically relevant phenotypes. Here, we will highlight the development of functional genomics in C. albicans and discuss the use of these approaches to addressing both commensalism and pathogenesis in this species.

Intermolecular epistasis shaped the function and evolution of an ancient transcription factor and its DNA binding sites

Complexes of specifically interacting molecules, such as transcription factor proteins (TFs) and the DNA response elements (REs) they recognize, control most biological processes, but little is known concerning the functional and evolutionary effects of epistatic interactions across molecular interfaces. We experimentally characterized all combinations of genotypes in the joint protein-DNA sequence space defined by an historical transition in TF-RE specificity that occurred some 500 million years ago in the DNA-binding domain of an ancient steroid hormone receptor. We found that rampant epistasis within and between the two molecules was essential to specific TF-RE recognition and to the evolution of a novel TF-RE complex with unique derived specificity. Permissive and restrictive epistatic mutations across the TF-RE interface opened and closed potential evolutionary paths accessible by the other, making the evolution of each molecule contingent on its partner's history and allowing a molecular complex with novel specificity to evolve.

Mechanically Induced Trapping of Molecular Interactions and Its Applications

Measuring binding affinities and association/dissociation rates of molecular interactions is important for a quantitative understanding of cellular mechanisms. Many low-throughput methods have been developed throughout the years to obtain these parameters. Acquiring data with higher accuracy and throughput is, however, necessary to characterize complex biological networks. Here, we provide an overview of a high-throughput microfluidic method based on mechanically induced trapping of molecular interactions (MITOMI). MITOMI can be used to obtain affinity constants and kinetic rates of hundreds of protein-ligand interactions in parallel. It has been used in dozens of studies to measure binding affinities of transcription factors, map protein interaction networks, identify pharmacological inhibitors, and perform high-throughput, low-cost molecular diagnostics. This article covers the technological aspects of MITOMI and its applications.

How do regulatory networks evolve and expand throughout evolution?

Throughout evolution, regulatory networks need to expand and adapt to accommodate novel genes and gene functions. However, the molecular details explaining how gene networks evolve remain largely unknown. Recent studies demonstrate that changes in transcription factors contribute to the evolution of regulatory networks. In particular, duplication of transcription factors followed by specific mutations in their DNA-binding or interaction domains propels the divergence and emergence of new networks. The innate promiscuity and modularity of regulatory networks contributes to their evolvability: duplicated promiscuous regulators and their target promoters can acquire mutations that lead to gradual increases in specificity, allowing neofunctionalization or subfunctionalization.

Evolution of transcription factor function as a mechanism for changing metazoan developmental gene regulatory networks

The form that an animal takes during development is directed by gene regulatory networks (GRNs). Developmental GRNs interpret maternally deposited molecules and externally supplied signals to direct cell-fate decisions, which ultimately leads to the arrangements of organs and tissues in the organism. Genetically encoded modifications to these networks have generated the wide range of metazoan diversity that exists today. Most studies of GRN evolution focus on changes to cis-regulatory DNA, and it was historically theorized that changes to the transcription factors that bind to these cis-regulatory modules (CRMs) contribute to this process only rarely. A growing body of evidence suggests that changes to the coding regions of transcription factors play a much larger role in the evolution of developmental gene regulatory networks than originally imagined. Just as cis-regulatory changes make use of modular binding site composition and tissue-specific modules to avoid pleiotropy, transcription factor coding regions also predominantly evolve in ways that limit the context of functional effects. Here, we review the recent works that have led to this unexpected change in the field of Evolution and Development (Evo-Devo) and consider the implications these studies have had on our understanding of the evolution of developmental processes.

Reshuffling transcriptional circuits: how microorganisms adapt to colonize the human body

Several hundred taxa of microorganisms-including bacteria, archaea and eukaryotes-inhabit the human body. What did it take for these species to become stable residents of humans? Recent reports illustrate how evolutionary changes in transcriptional circuits played a pivotal role in the adaptation of single-celled eukaryotes to colonize mammals.

Gene duplication of type-B ARR transcription factors systematically extends transcriptional regulatory structures in Arabidopsis

Many of duplicated genes are enriched in signaling pathways. Recently, gene duplication of kinases has been shown to provide genetic buffering and functional diversification in cellular signaling. Transcription factors (TFs) are also often duplicated. However, how duplication of TFs affects their regulatory structures and functions of target genes has not been explored at the systems level. Here, we examined regulatory and functional roles of duplication of three major ARR TFs (ARR1, 10, and 12) in Arabidopsis cytokinin signaling using wild-type and single, double, and triple deletion mutants of the TFs. Comparative analysis of gene expression profiles obtained from Arabidopsis roots in wild-type and these mutants showed that duplication of ARR TFs systematically extended their transcriptional regulatory structures, leading to enhanced robustness and diversification in functions of target genes, as well as in regulation of cellular networks of target genes. Therefore, our results suggest that duplication of TFs contributes to robustness and diversification in functions of target genes by extending transcriptional regulatory structures.

Duplication of a promiscuous transcription factor drives the emergence of a new regulatory network

The emergence of new genes throughout evolution requires rewiring and extension of regulatory networks. However, the molecular details of how the transcriptional regulation of new gene copies evolves remain largely unexplored. Here we show how duplication of a transcription factor gene allowed the emergence of two independent regulatory circuits. Interestingly, the ancestral transcription factor was promiscuous and could bind different motifs in its target promoters. After duplication, one paralogue evolved increased binding specificity so that it only binds one type of motif, whereas the other copy evolved a decreased activity so that it only activates promoters that contain multiple binding sites. Interestingly, only a few mutations in both the DNA-binding domains and in the promoter binding sites were required to gradually disentangle the two networks. These results reveal how duplication of a promiscuous transcription factor followed by concerted cis and trans mutations allows expansion of a regulatory network.