From elements to modules: regulatory evolution in Ascomycota fungi

Polyphasic Characterization of the Biocontrol Potential of a Novel Strain of Trichoderma atroviride Isolated from Central Mexico

This work describes the characterization of Trichoderma atroviride strain CMU-08, isolated from Michoacán, Mexico. CMU-08 demonstrated robust growth and conidiation across a temperature range from 16 to 32 °C and a pH range from 4 to 9 on potato dextrose agar (PDA) and malt extract agar (MEA) media. The strain is an efficient antagonist of six species of phytopathogenic fungi and oomycetes in PDA, MEA, and Vogel minimal medium (VMM). Antagonist mechanisms of CMU-08 included direct mycoparasitism observed in dual-culture assays, as well as antibiosis attributed to growth inhibition via both volatile and non-volatile metabolites, with the effectiveness varying depending on the test phytopathogen and culture medium. Extracellular filtrates (ECFs) recovered from liquid cultures of CMU-08 under basal and induced conditions using Botrytis cinerea cell walls significantly inhibited their growth at a concentration of 750 µg/mL. Moreover, in detached tomato leaf assays, these ECFs reduced foliar damage caused by B. cinerea by 24-34%. The volatile organic compounds (VOCs) produced by CMU-08 also exhibited substantial efficacy, reducing foliar damage by up to 50% in similar tests. Despite showing no basal extracellular chitinase enzymatic activity, CMU-08 demonstrated significant induction of this activity in cultures supplemented with B. cinerea and Fusarium sp. cell walls. Four genes encoding extracellular chitinases (chit33, chit36, ech42, and locus 217415) showed different dynamics of transcriptional regulation during the dual-culture confrontation of strain CMU-08 with B. cinerea and Fusarium sp., varying according to the phytopathogen and the interaction stage. The CMU-08 strain shows physiological versatility and employs a variety of antagonist mechanisms toward different species of phytopathogenic microorganisms, making it a good candidate for developing a biocontrol product for field application.

Variations in Rhizospheric and Endophytic Root Fungal Communities of Scrophularia ningpoensis in Different Producing Areas

Few studies have examined the association of factors associated with soil fertility and composition with the structure of microbial communities in the rhizosphere and endosphere. Hence, this study aimed to explore the effects of geographical differences on fungal communities in the roots of Scrophularia ningpoensis and the relationship between the fungal communities and secondary metabolic components in the host plant. We found that there was greater diversity in the fungal communities of the rhizosphere compartment than in endosphere communities. Ascomycota and Basidiomycota were dominant among the endosphere fungi, whereas Mortierellomycota was distributed in the rhizosphere. The composition of bulk soil obtained from different producing areas was significantly different, and the correlation between the rhizospheric and physicochemical compartments of the soil was higher than that observed with the endophytic compartment. Redundancy analysis and canonical correspondence analysis of the rhizospheric and endophytic samples revealed that the organic matter, total organic carbon, total nitrogen, and Hg levels were adequately correlated with the composition of rhizospheric and endophytic fungal communities. Multiple linear regression analyses facilitated the identification of potentially beneficial fungi whose abundance was correlated with levels of secondary metabolites, such as harpagide and harpagoside. These fungi could potentially provide valuable information regarding the use of S. ningpoensis in the medicinal plant industry.

Rapid evolutionary repair by secondary perturbation of a primary disrupted transcriptional network

The discrete steps of transcriptional rewiring have been proposed to occur neutrally to ensure steady gene expression under stabilizing selection. A conflict-free switch of a regulon between regulators may require an immediate compensatory evolution to minimize deleterious effects. Here, we perform an evolutionary repair experiment on the Lachancea kluyveri yeast sef1Δ mutant using a suppressor development strategy. Complete loss of SEF1 forces cells to initiate a compensatory process for the pleiotropic defects arising from misexpression of TCA cycle genes. Using different selective conditions, we identify two adaptive loss-of-function mutations of IRA1 and AZF1. Subsequent analyses show that Azf1 is a weak transcriptional activator regulated by the Ras1-PKA pathway. Azf1 loss-of-function triggers extensive gene expression changes responsible for compensatory, beneficial, and trade-off phenotypes. The trade-offs can be alleviated by higher cell density. Our results not only indicate that secondary transcriptional perturbation provides rapid and adaptive mechanisms potentially stabilizing the initial stage of transcriptional rewiring but also suggest how genetic polymorphisms of pleiotropic mutations could be maintained in the population.

Species-specific effects of the introduction of Aspergillus nidulans gfdB in osmophilic aspergilli

Industrial fungi need a strong environmental stress tolerance to ensure acceptable efficiency and yields. Previous studies shed light on the important role that Aspergillus nidulans gfdB, putatively encoding a NAD⁺-dependent glycerol-3-phosphate dehydrogenase, plays in the oxidative and cell wall integrity stress tolerance of this filamentous fungus model organism. The insertion of A. nidulans gfdB into the genome of Aspergillus glaucus strengthened the environmental stress tolerance of this xerophilic/osmophilic fungus, which may facilitate the involvement of this fungus in various industrial and environmental biotechnological processes. On the other hand, the transfer of A. nidulans gfdB to Aspergillus wentii, another promising industrial xerophilic/osmophilic fungus, resulted only in minor and sporadic improvement in environmental stress tolerance and meanwhile partially reversed osmophily. Because A. glaucus and A. wentii are phylogenetically closely related species and both fungi lack a gfdB ortholog, these results warn us that any disturbance of the stress response system of the aspergilli may elicit rather complex and even unforeseeable, species-specific physiological changes. This should be taken into consideration in any future targeted industrial strain development projects aiming at the fortification of the general stress tolerance of these fungi. KEY POINTS: • A. wentii c' gfdB strains showed minor and sporadic stress tolerance phenotypes. • The osmophily of A. wentii significantly decreased in the c' gfdB strains. • Insertion of gfdB caused species-specific phenotypes in A. wentii and A. glaucus.

Evolution of gene regulatory network of C4 photosynthesis in the genus Flaveria reveals the evolutionary status of C3-C4 intermediate species

C₄ photosynthesis evolved from ancestral C₃ photosynthesis by recruiting pre-existing genes to fulfill new functions. The enzymes and transporters required for the C₄ metabolic pathway have been intensively studied and well documented; however, the transcription factors (TFs) that regulate these C₄ metabolic genes are not yet well understood. In particular, how the TF regulatory network of C₄ metabolic genes was rewired during the evolutionary process is unclear. Here, we constructed gene regulatory networks (GRNs) for four closely evolutionarily related species from the genus Flaveria, which represent four different evolutionary stages of C₄ photosynthesis: C₃ (F. robusta), type I C₃-C₄ (F. sonorensis), type II C₃-C₄ (F. ramosissima), and C₄ (F. trinervia). Our results show that more than half of the co-regulatory relationships between TFs and core C₄ metabolic genes are species specific. The counterparts of the C₄ genes in C₃ species were already co-regulated with photosynthesis-related genes, whereas the required TFs for C₄ photosynthesis were recruited later. The TFs involved in C₄ photosynthesis were widely recruited in the type I C₃-C₄ species; nevertheless, type II C₃-C₄ species showed a divergent GRN from C₄ species. In line with these findings, a ¹³CO₂ pulse-labeling experiment showed that the CO₂ initially fixed into C₄ acid was not directly released to the Calvin-Benson-Bassham cycle in the type II C₃-C₄ species. Therefore, our study uncovered dynamic changes in C₄ genes and TF co-regulation during the evolutionary process; furthermore, we showed that the metabolic pathway of the type II C₃-C₄ species F. ramosissima represents an alternative evolutionary solution to the ammonia imbalance in C₃-C₄ intermediate species.

The transcriptional response to oxidative stress is independent of stress-granule formation

Cells respond to stress with translational arrest, robust transcriptional changes, and transcription-independent formation of mRNP assemblies termed stress granules (SGs). Despite considerable interest in the role of SGs in oxidative, unfolded protein and viral stress responses, whether and how SGs contribute to stress-induced transcription have not been rigorously examined. To address this, we characterized transcriptional changes in Drosophila S2 cells induced by acute oxidative-stress and assessed how these were altered under conditions that disrupted SG assembly. Oxidative stress for 3 h predominantly resulted in induction or up-regulation of stress-responsive mRNAs whose levels peaked during recovery after stress cessation. The stress transcriptome is enriched in mRNAs coding for chaperones including HSP70s, small heat shock proteins, glutathione transferases, and several noncoding RNAs. Oxidative stress also induced cytoplasmic SGs that disassembled 3 h after stress cessation. As expected, RNAi-mediated knockdown of the conserved G3BP1/Rasputin protein inhibited SG assembly. However, this disruption had no significant effect on the stress-induced transcriptional response or stress-induced translational arrest. Thus SG assembly and stress-induced gene expression alterations appear to be driven by distinctive signaling processes. We suggest that while SG assembly represents a fast, transient mechanism, the transcriptional response enables a slower, longer-lasting mechanism for adaptation to and recovery from cell stress.

Identification of Taxonomically Restricted Transcripts from Illumina RNA Sequencing Data

In order to perform a well-balanced comparative transcriptomic analysis, the reference genome and annotations for all species included in the comparison must be of a similar quality and completeness. Frequently, comparative transcriptomic analyses include non-model organisms whose annotations are not as well curated; this inequality can lead to biases.To avoid potential biases stemming from incomplete annotations, a comparative transcriptomic analysis can incorporate de novo transcriptome assemblies for each species, which reduces this disparity. This chapter covers all of the steps which are necessary to run a comparative transcriptomic analysis with de novo transcriptome assemblies, from the first step of the experimental design to the sequencing, and ultimately the bioinformatic analysis.

Homeobox Transcription Factors Are Required for Fungal Development and the Suppression of Host Defense Mechanisms in the Colletotrichum scovillei-Pepper Pathosystem

Colletotrichum scovillei, an ascomycete phytopathogenic fungus, is the main causal agent of serious yield losses of economic crops worldwide. The fungus causes anthracnose disease on several fruits, including peppers. However, little is known regarding the underlying molecular mechanisms involved in the development of anthracnose caused by this fungus. In an initial step toward understanding the development of anthracnose on pepper fruits, we retrieved 624 transcription factors (TFs) from the whole genome of C. scovillei and comparatively analyzed the entire repertoire of TFs among phytopathogenic fungi. Evolution and proliferation of members of the homeobox-like superfamily, including homeobox (HOX) TFs that regulate the development of eukaryotic organisms, were demonstrated in the genus Colletotrichum. C. scovillei was found to contain 10 HOX TF genes (CsHOX1 to CsHOX10), which were functionally characterized using deletion mutants of each CsHOX gene. Notably, CsHOX1 was identified as a pathogenicity factor required for the suppression of host defense mechanisms, which represents a new role for HOX TFs in pathogenic fungi. CsHOX2 and CsHOX7 were found to play essential roles in conidiation and appressorium development, respectively, in a stage-specific manner in C. scovillei. Our study provides a molecular basis for understanding the mechanisms associated with the development of anthracnose on fruits caused by C. scovillei, which will aid in the development of novel approaches for disease management. IMPORTANCE The ascomycete phytopathogenic fungus, Colletotrichum scovillei, causes serious yield loss on peppers. However, little is known about molecular mechanisms involved in the development of anthracnose caused by this fungus. We analyzed whole-genome sequences of C. scovillei and isolated 624 putative TFs, revealing the existence of 10 homeobox (HOX) transcription factor (TF) genes. We found that CsHOX1 is a pathogenicity factor required for the suppression of host defense mechanism, which represents a new role for HOX TFs in pathogenic fungi. We also found that CsHOX2 and CsHOX7 play essential roles in conidiation and appressorium development, respectively, in a stage-specific manner in C. scovillei. Our study contributes to understanding the mechanisms associated with the development of anthracnose on fruits caused by C. scovillei, which will aid for initiating novel approaches for disease management.

GhLUX1 and GhELF3 Are Two Components of the Circadian Clock That Regulate Flowering Time of Gossypium hirsutum

Photoperiod is an important external factor that regulates flowering time, the core mechanism of which lies in the circadian clock-controlled expression of FLOWERING LOCUS T (FT) and its upstream regulators. However, the roles of the circadian clock in regulating cotton flowering time are largely unknown. In this study, we cloned two circadian clock genes in cotton, GhLUX1 and GhELF3. The physicochemical and structural properties of their putative proteins could satisfy the prerequisites for the interaction between them, which was proved by yeast two-hybrid (Y2H) and Bimolecular Fluorescent Complimentary (BiFC) assays. Phylogenetic analysis of LUXs and ELF3s indicated that the origin of LUXs was earlier than that of ELF3s, but ELF3s were more divergent and might perform more diverse functions. GhLUX1, GhELF3, GhCOL1, and GhFT exhibited rhythmic expression and were differentially expressed in the early flowering and late-flowering cotton varieties under different photoperiod conditions. Both overexpression of GhLUX1 and overexpression of GhELF3 in Arabidopsis delayed flowering probably by changing the oscillation phases and amplitudes of the key genes in the photoperiodic flowering pathway. Both silencing of GhLUX1 and silencing of GhELF3 in cotton increased the expression of GhCOL1 and GhFT and resulted in early flowering. In summary, the circadian clock genes were involved in regulating cotton flowering time and could be the candidate targets for breeding early maturing cotton varieties.

A network-based comparative framework to study conservation and divergence of proteomes in plant phylogenies

Comparative functional genomics offers a powerful approach to study species evolution. To date, the majority of these studies have focused on the transcriptome in mammalian and yeast phylogenies. Here, we present a novel multi-species proteomic dataset and a computational pipeline to systematically compare the protein levels across multiple plant species. Globally we find that protein levels diverge according to phylogenetic distance but is more constrained than the mRNA level. Module-level comparative analysis of groups of proteins shows that proteins that are more highly expressed tend to be more conserved. To interpret the evolutionary patterns of conservation and divergence, we develop a novel network-based integrative analysis pipeline that combines publicly available transcriptomic datasets to define co-expression modules. Our analysis pipeline can be used to relate the changes in protein levels to different species-specific phenotypic traits. We present a case study with the rhizobia-legume symbiosis process that supports the role of autophagy in this symbiotic association.

Invariable stoichiometry of ribosomal proteins in mouse brain tissues with aging

Across phyla, the ribosomes-the central molecular machines for translation of genetic information-exhibit an overall preserved architecture and a conserved functional core. The natural heterogeneity of the ribosome periodically phases a debate on their functional specialization and the tissue-specific variations of the ribosomal protein (RP) pool. Using sensitive differential proteomics, we performed a thorough quantitative inventory of the protein composition of ribosomes from 3 different mouse brain tissues, i.e., hippocampus, cortex, and cerebellum, across various ages, i.e., juvenile, adult, and middle-aged mouse groups. In all 3 brain tissues, in both monosomal and polysomal ribosome fractions, we detected an invariant set of 72 of 79 core RPs, RACK1 and 2 of the 8 RP paralogs, the stoichiometry of which remained constant across different ages. The amount of a few RPs punctually varied in either one tissue or one age group, but these fluctuations were within the tight bounds of the measurement noise. Further comparison with the ribosomes from a high-metabolic-rate organ, e.g., the liver, revealed protein composition identical to that of the ribosomes from the 3 brain tissues. Together, our data show an invariant protein composition of ribosomes from 4 tissues across different ages of mice and support the idea that functional heterogeneity may arise from factors other than simply ribosomal protein stoichiometry.

Transcriptional Regulation in Archaea: From Individual Genes to Global Regulatory Networks

Archaea are major contributors to biogeochemical cycles, possess unique metabolic capabilities, and resist extreme stress. To regulate the expression of genes encoding these unique programs, archaeal cells use gene regulatory networks (GRNs) composed of transcription factor proteins and their target genes. Recent developments in genetics, genomics, and computational methods used with archaeal model organisms have enabled the mapping and prediction of global GRN structures. Experimental tests of these predictions have revealed the dynamical function of GRNs in response to environmental variation. Here, we review recent progress made in this area, from investigating the mechanisms of transcriptional regulation of individual genes to small-scale subnetworks and genome-wide global networks. At each level, archaeal GRNs consist of a hybrid of bacterial, eukaryotic, and uniquely archaeal mechanisms. We discuss this theme from the perspective of the role of individual transcription factors in genome-wide regulation, how these proteins interact to compile GRN topological structures, and how these topologies lead to emergent, high-level GRN functions. We conclude by discussing how systems biology approaches are a fruitful avenue for addressing remaining challenges, such as discovering gene function and the evolution of GRNs.

The ancestral levels of transcription and the evolution of sexual phenotypes in filamentous fungi

Changes in gene expression have been hypothesized to play an important role in the evolution of divergent morphologies. To test this hypothesis in a model system, we examined differences in fruiting body morphology of five filamentous fungi in the Sordariomycetes, culturing them in a common garden environment and profiling genome-wide gene expression at five developmental stages. We reconstructed ancestral gene expression phenotypes, identifying genes with the largest evolved increases in gene expression across development. Conducting knockouts and performing phenotypic analysis in two divergent species typically demonstrated altered fruiting body development in the species that had evolved increased expression. Our evolutionary approach to finding relevant genes proved far more efficient than other gene deletion studies targeting whole genomes or gene families. Combining gene expression measurements with knockout phenotypes facilitated the refinement of Bayesian networks of the genes underlying fruiting body development, regulation of which is one of the least understood processes of multicellular development.

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.

Recurrent rewiring and emergence of RNA regulatory networks

Alterations in regulatory networks contribute to evolutionary change. Transcriptional networks are reconfigured by changes in the binding specificity of transcription factors and their cognate sites. The evolution of RNA-protein regulatory networks is far less understood. The PUF (Pumilio and FBF) family of RNA regulatory proteins controls the translation, stability, and movements of hundreds of mRNAs in a single species. We probe the evolution of PUF-RNA networks by direct identification of the mRNAs bound to PUF proteins in budding and filamentous fungi and by computational analyses of orthologous RNAs from 62 fungal species. Our findings reveal that PUF proteins gain and lose mRNAs with related and emergent biological functions during evolution. We demonstrate at least two independent rewiring events for PUF3 orthologs, independent but convergent evolution of PUF4/5 binding specificity and the rewiring of the PUF4/5 regulons in different fungal lineages. These findings demonstrate plasticity in RNA regulatory networks and suggest ways in which their rewiring occurs.

Natural gene expression variation studies in yeast

The rise of sequence information across different yeast species and strains is driving an increasing number of studies in the emerging field of genomics to associate polymorphic variants, mRNA abundance and phenotypic differences between individuals. Here, we gathered evidence from recent studies covering several layers that define the genotype-phenotype gap, such as mRNA abundance, allele-specific expression and translation efficiency to demonstrate how genetic variants co-evolve and define an individual's genome. Moreover, we exposed several antecedents where inter- and intra-specific studies led to opposite conclusions, probably owing to genetic divergence. Future studies in this area will benefit from the access to a massive array of well-annotated genomes and new sequencing technologies, which will allow the fine breakdown of the complex layers that delineate the genotype-phenotype map. Copyright © 2016 John Wiley & Sons, Ltd.

Transcriptional rewiring over evolutionary timescales changes quantitative and qualitative properties of gene expression

Evolutionary changes in transcription networks are an important source of diversity across species, yet the quantitative consequences of network evolution have rarely been studied. Here we consider the transcriptional 'rewiring' of the three GAL genes that encode the enzymes needed for cells to convert galactose to glucose. In Saccharomyces cerevisiae, the transcriptional regulator Gal4 binds and activates these genes. In the human pathogen Candida albicans (which last shared a common ancestor with S. cerevisiae some 300 million years ago), we show that different regulators, Rtg1 and Rtg3, activate the three GAL genes. Using single-cell dynamics and RNA-sequencing, we demonstrate that although the overall logic of regulation is the same in both species-the GAL genes are induced by galactose-there are major differences in both the quantitative response of these genes to galactose and in the position of these genes in the overall transcription network structure of the two species.

Reconstruction and Analysis of the Evolution of Modular Transcriptional Regulatory Programs Using Arboretum

Comparative functional genomics aims to measure and compare genome-wide functional data such as transcriptomes, proteomes, and epigenomes across multiple species to study the conservation and divergence patterns of such quantitative measurements. However, computational methods to systematically compare these quantitative genomic profiles across multiple species are in their infancy. We developed Arboretum, a novel algorithm to identify modules of co-expressed genes and trace their evolutionary history across multiple species from a complex phylogeny. To interpret the results from Arboretum we developed several measures to examine the extent of conservation and divergence in modules and their relationship to species lifestyle, cis-regulatory elements, and gene duplication. We applied Arboretum to study the evolution of modular transcriptional regulatory programs controlling transcriptional response to different environmental stresses in the yeast Ascomycota phylogeny. We found that modules of similar patterns of expression captured the transcriptional responses to different stresses across species; however, the genes exhibiting these patterns were not the same. Divergence in module membership was associated with changes in lifestyle and specific clades and that gene duplication was a major factor contributing to the divergence of module membership.

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.

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.

Comparative analysis of gene regulatory networks: from network reconstruction to evolution

Regulation of gene expression is central to many biological processes. Although reconstruction of regulatory circuits from genomic data alone is therefore desirable, this remains a major computational challenge. Comparative approaches that examine the conservation and divergence of circuits and their components across strains and species can help reconstruct circuits as well as provide insights into the evolution of gene regulatory processes and their adaptive contribution. In recent years, advances in genomic and computational tools have led to a wealth of methods for such analysis at the sequence, expression, pathway, module, and entire network level. Here, we review computational methods developed to study transcriptional regulatory networks using comparative genomics, from sequence to functional data. We highlight how these methods use evolutionary conservation and divergence to reliably detect regulatory components as well as estimate the extent and rate of divergence. Finally, we discuss the promise and open challenges in linking regulatory divergence to phenotypic divergence and adaptation.

Metabolic regulation in model ascomycetes--adjusting similar genomes to different lifestyles

The related yeasts Saccharomyces cerevisiae and Candida albicans have similar genomes but very different lifestyles. These fungi have modified transcriptional and post-translational regulatory processes to adapt their similar genomes to the distinct biological requirements of the two yeasts. We review recent findings comparing the differences between these species, highlighting how they have achieved specialized metabolic capacities tailored to their lifestyles despite sharing similar genomes. Studying this transcriptional and post-transcriptional rewiring may improve our ability to interpret phenotype from genotype.

Insights into molecular evolution from yeast genomics

Enabled by comparative genomics, yeasts have increasingly developed into a powerful model system for molecular evolution. Here we survey several areas in which yeast studies have made important contributions, including regulatory evolution, gene duplication and divergence, evolution of gene order and evolution of complexity. In each area we highlight key studies and findings based on techniques ranging from statistical analysis of large datasets to direct laboratory measurements of fitness. Future work will combine traditional evolutionary genetics analysis and experimental evolution with tools from systems biology to yield mechanistic insight into complex phenotypes.

The transcription factors Tec1 and Ste12 interact with coregulators Msa1 and Msa2 to activate adhesion and multicellular development

In Saccharomyces cerevisiae and related yeast species, the TEA transcription factor Tec1, together with a second transcription factor, Ste12, controls development, including cell adhesion and filament formation. Tec1-Ste12 complexes control target genes through Tec1 binding sites (TEA consensus sequences [TCSs]) that can be further combined with Ste12 binding sites (pheromone response elements [PREs]) for cooperative DNA binding. The activity of Tec1-Ste12 complexes is known to be under negative control of the Dig1 and Dig2 (Dig1/2) transcriptional corepressors that confer regulation by upstream signaling pathways. Here, we found that Tec1 and Ste12 can associate with the transcriptional coregulators Msa1 and Msa2 (Msa1/2), which were previously found to associate with the cell cycle transcription factor complexes SBF (Swi4/Swi6 cell cycle box binding factor) and MBF (Mbp1/Swi6 cell cycle box binding factor) to control G1-specific transcription. We further show that Tec1-Ste12-Msa1/2 complexes (i) do not contain Swi4 or Mbp1, (ii) assemble at single TCSs or combined TCS-PREs in vitro, and (iii) coregulate genes involved in adhesive and filamentous growth by direct promoter binding in vivo. Finally, we found that, in contrast to Dig proteins, Msa1/2 seem to act as coactivators that enhance the transcriptional activity of Tec1-Ste12. Taken together, our findings add an additional layer of complexity to our understanding of the control mechanisms exerted by the evolutionarily conserved TEA domain and Ste12-like transcription factors.

Integrated module and gene-specific regulatory inference implicates upstream signaling networks

Regulatory networks that control gene expression are important in diverse biological contexts including stress response and development. Each gene's regulatory program is determined by module-level regulation (e.g. co-regulation via the same signaling system), as well as gene-specific determinants that can fine-tune expression. We present a novel approach, Modular regulatory network learning with per gene information (MERLIN), that infers regulatory programs for individual genes while probabilistically constraining these programs to reveal module-level organization of regulatory networks. Using edge-, regulator- and module-based comparisons of simulated networks of known ground truth, we find MERLIN reconstructs regulatory programs of individual genes as well or better than existing approaches of network reconstruction, while additionally identifying modular organization of the regulatory networks. We use MERLIN to dissect global transcriptional behavior in two biological contexts: yeast stress response and human embryonic stem cell differentiation. Regulatory modules inferred by MERLIN capture co-regulatory relationships between signaling proteins and downstream transcription factors thereby revealing the upstream signaling systems controlling transcriptional responses. The inferred networks are enriched for regulators with genetic or physical interactions, supporting the inference, and identify modules of functionally related genes bound by the same transcriptional regulators. Our method combines the strengths of per-gene and per-module methods to reveal new insights into transcriptional regulation in stress and development.

Evolutionary principles of modular gene regulation in yeasts

Divergence in gene regulation can play a major role in evolution. Here, we used a phylogenetic framework to measure mRNA profiles in 15 yeast species from the phylum Ascomycota and reconstruct the evolution of their modular regulatory programs along a time course of growth on glucose over 300 million years [corrected]. We found that modules have diverged proportionally to phylogenetic distance, with prominent changes in gene regulation accompanying changes in lifestyle and ploidy, especially in carbon metabolism. Paralogs have significantly contributed to regulatory divergence, typically within a very short window from their duplication. Paralogs from a whole genome duplication (WGD) event have a uniquely substantial contribution that extends over a longer span. Similar patterns occur when considering the evolution of the heat shock regulatory program measured in eight of the species, suggesting that these are general evolutionary principles. DOI:http://dx.doi.org/10.7554/eLife.00603.001.

Arboretum: reconstruction and analysis of the evolutionary history of condition-specific transcriptional modules

Comparative functional genomics studies the evolution of biological processes by analyzing functional data, such as gene expression profiles, across species. A major challenge is to compare profiles collected in a complex phylogeny. Here, we present Arboretum, a novel scalable computational algorithm that integrates expression data from multiple species with species and gene phylogenies to infer modules of coexpressed genes in extant species and their evolutionary histories. We also develop new, generally applicable measures of conservation and divergence in gene regulatory modules to assess the impact of changes in gene content and expression on module evolution. We used Arboretum to study the evolution of the transcriptional response to heat shock in eight species of Ascomycota fungi and to reconstruct modules of the ancestral environmental stress response (ESR). We found substantial conservation in the stress response across species and in the reconstructed components of the ancestral ESR modules. The greatest divergence was in the most induced stress, primarily through module expansion. The divergence of the heat stress response exceeds that observed in the response to glucose depletion in the same species. Arboretum and its associated analyses provide a comprehensive framework to systematically study regulatory evolution of condition-specific responses.

Beyond Arabidopsis: the circadian clock in non-model plant species

Circadian clocks allow plants to temporally coordinate many aspects of their biology with the diurnal cycle derived from the rotation of Earth on its axis. Although there is a rich history of the study of clocks in many plant species, in recent years much progress in elucidating the architecture and function of the plant clock has emerged from studies of the model plant, Arabidopsis thaliana. There is considerable interest in extending this knowledge of the circadian clock into diverse plant species in order to address its role in topics as varied as agricultural productivity and the responses of individual species and plant communities to global climate change and environmental degradation. The analysis of circadian clocks in the green lineage provides insight into evolutionary processes in plants and throughout the eukaryotes.

Saccharomyces diversity and evolution: a budding model genus

Saccharomyces cerevisiae is one of the best-understood and most powerful genetic model systems. Several disciplines are now converging to turn Saccharomyces into an exciting model genus for evolutionary genetics and genomics. Yeast taxonomists and ecologists have dramatically expanded and clarified Saccharomyces diversity, more than doubling the number of bona fide species since 2000. High-quality genome sequences are available (or soon will be) for all seven known species. Haploid laboratory strains are enabling a deep integration of classic genetic approaches with modern genomic tools. Population genomic surveys and quantitative trait mapping of variation within species are underway across the genus. Finally, several case studies have illuminated general and novel genetic mechanisms of evolution. Expanding strain collections, low-cost genome sequencing, and tools for precise genetic manipulation promise to usher in a golden era for this surprisingly diverse genus as an evolutionary model.

Protein modularity, cooperative binding, and hybrid regulatory states underlie transcriptional network diversification

We examine how different transcriptional network structures can evolve from an ancestral network. By characterizing how the ancestral mode of gene regulation for genes specific to a-type cells in yeast species evolved from an activating paradigm to a repressing one, we show that regulatory protein modularity, conversion of one cis-regulatory sequence to another, distribution of binding energy among protein-protein and protein-DNA interactions, and exploitation of ancestral network features all contribute to the evolution of a novel regulatory mode. The formation of this derived mode of regulation did not disrupt the ancestral mode and thereby created a hybrid regulatory state where both means of transcription regulation (ancestral and derived) contribute to the conserved expression pattern of the network. Finally, we show how this hybrid regulatory state has resolved in different ways in different lineages to generate the diversity of regulatory network structures observed in modern species.

Coevolution within and between regulatory loci can preserve promoter function despite evolutionary rate acceleration

Phenotypes that appear to be conserved could be maintained not only by strong purifying selection on the underlying genetic systems, but also by stabilizing selection acting via compensatory mutations with balanced effects. Such coevolution has been invoked to explain experimental results, but has rarely been the focus of study. Conserved expression driven by the unc-47 promoters of Caenorhabditis elegans and C. briggsae persists despite divergence within a cis-regulatory element and between this element and the trans-regulatory environment. Compensatory changes in cis and trans are revealed when these promoters are used to drive expression in the other species. Functional changes in the C. briggsae promoter, which has experienced accelerated sequence evolution, did not lead to alteration of gene expression in its endogenous environment. Coevolution among promoter elements suggests that complex epistatic interactions within cis-regulatory elements may facilitate their divergence. Our results offer a detailed picture of regulatory evolution in which subtle, lineage-specific, and compensatory modifications of interacting cis and trans regulators together maintain conserved gene expression patterns.

Some phenotypes, including gene expression patterns, are conserved between distantly related species. However, the molecular bases of those phenotypes are not necessarily conserved. Instead, regulatory DNA sequences and the proteins with which they interact can change over time with balanced effects, preserving expression patterns and concealing regulatory divergence. Coevolution between interacting molecules makes gene regulation highly species-specific, and it can be detected when the cis-regulatory DNA of one species is used to drive expression in another species. In this way, we identified regions of the C. elegans and C. briggsae unc-47 promoters that have coevolved with the lineage-specific trans-regulatory environments of these organisms. The C. briggsae promoter experienced accelerated sequence change relative to related species. All of this evolution occurred without changing the expression pattern driven by the promoter in its endogenous environment.

A functional selection model explains evolutionary robustness despite plasticity in regulatory networks

Evolutionary rewiring of regulatory networks is an important source of diversity among species. Previous evidence suggested substantial divergence of regulatory networks across species. However, systematically assessing the extent of this plasticity and its functional implications has been challenging due to limited experimental data and the noisy nature of computational predictions. Here, we introduce a novel approach to study cis-regulatory evolution, and use it to trace the regulatory history of 88 DNA motifs of transcription factors across 23 Ascomycota fungi. While motifs are conserved, we find a pervasive gain and loss in the regulation of their target genes. Despite this turnover, the biological processes associated with a motif are generally conserved. We explain these trends using a model with a strong selection to conserve the overall function of a transcription factor, and a much weaker selection over the specific genes it targets. The model also accounts for the turnover of bound targets measured experimentally across species in yeasts and mammals. Thus, selective pressures on regulatory networks mostly tolerate local rewiring, and may allow for subtle fine-tuning of gene regulation during evolution.

Yeast: an experimental organism for 21st Century biology

In this essay, we revisit the status of yeast as a model system for biology. We first summarize important contributions of yeast to eukaryotic biology that we anticipated in 1988 in our first article on the subject. We then describe transformative developments that we did not anticipate, most of which followed the publication of the complete genomic sequence of Saccharomyces cerevisiae in 1996. In the intervening 23 years it appears to us that yeast has graduated from a position as the premier model for eukaryotic cell biology to become the pioneer organism that has facilitated the establishment of the entirely new fields of study called "functional genomics" and "systems biology." These new fields look beyond the functions of individual genes and proteins, focusing on how these interact and work together to determine the properties of living cells and organisms.

A conserved developmental patterning network produces quantitatively different output in multiple species of Drosophila

Differences in the level, timing, or location of gene expression can contribute to alternative phenotypes at the molecular and organismal level. Understanding the origins of expression differences is complicated by the fact that organismal morphology and gene regulatory networks could potentially vary even between closely related species. To assess the scope of such changes, we used high-resolution imaging methods to measure mRNA expression in blastoderm embryos of Drosophila yakuba and Drosophila pseudoobscura and assembled these data into cellular resolution atlases, where expression levels for 13 genes in the segmentation network are averaged into species-specific, cellular resolution morphological frameworks. We demonstrate that the blastoderm embryos of these species differ in their morphology in terms of size, shape, and number of nuclei. We present an approach to compare cellular gene expression patterns between species, while accounting for varying embryo morphology, and apply it to our data and an equivalent dataset for Drosophila melanogaster. Our analysis reveals that all individual genes differ quantitatively in their spatio-temporal expression patterns between these species, primarily in terms of their relative position and dynamics. Despite many small quantitative differences, cellular gene expression profiles for the whole set of genes examined are largely similar. This suggests that cell types at this stage of development are conserved, though they can differ in their relative position by up to 3–4 cell widths and in their relative proportion between species by as much as 5-fold. Quantitative differences in the dynamics and relative level of a subset of genes between corresponding cell types may reflect altered regulatory functions between species. Our results emphasize that transcriptional networks can diverge over short evolutionary timescales and that even small changes can lead to distinct output in terms of the placement and number of equivalent cells.

For a gene to function properly, it must be active in the right place, at the right time, and in the right amount. Changes in any of these features can lead to observable differences between individuals and species and in some cases can lead to disease. We do not currently understand how the position, timing, and amount of gene expression is encoded in DNA sequence. One approach to this problem is to compare how gene expression differs between species and to try to relate changes in DNA sequence to changes in gene expression. Here, we take the first step by comparing gene expression patterns at high spatial and temporal resolution between embryos of three species of fruit flies. We develop methods for comparing gene expression in individual cells, which allow us to control for variation in the size, shape, and number of nuclei between embryos. We find measurable quantitative differences in the patterns for all individual genes that we have examined. However, by considering all genes in our dataset at once, we show that many genes are changing together, leading to largely equivalent types of cells in these three species.

The evolution of gene expression regulatory networks in yeasts

Gene regulation is a major source of phenotypic diversity between and within species. This aspect of evolution has long been addressed from the sole point of view of the genome sequence. The incredible development of transcriptomics approaches now allows one to actually study the topology and the properties of regulatory networks on an evolutionary perspective. This new discipline is called comparative functional genomics or comparative transcriptomics. This article reviews some of the main advances made in this field, using yeast species, and especially the species sequenced in the frame of the Genolevures program, as a model.

Gene acquisition, duplication and metabolic specification: the evolution of fungal methylisocitrate lyases

Gene duplication represents an evolutionary mechanism for expanding metabolic potential. Here we analysed the evolutionary relatedness of isocitrate and methylisocitrate lyases, which are key enzymes of the glyoxylate and methylcitrate cycle respectively. Phylogenetic analyses imply that ancient eukaryotes acquired an isocitrate lyase gene from a prokaryotic source, but it was lost in some eukaryotic lineages. However, protists, oomycetes and most fungi maintained this gene and successfully integrated the corresponding enzyme into the glyoxylate cycle. A second gene, encoding a highly related enzyme, is present in fungi, but absent from other eukaryotes. This methylisocitrate lyase is specifically involved in propionyl-CoA degradation via the methylcitrate cycle. Although bacteria possess methylisocitrate lyases with a structural fold similar to that of isocitrate lyases, their sequence identity to fungal methylisocitrate lyases is low. Phylogenetic analyses imply that fungal methylisocitrate lyases arose from gene duplication of an ancient isocitrate lyase gene from the basidiomycete lineage. Mutagenesis of active-site residues of a bacterial and fungal isocitrate lyase, which have been predicted to direct the substrate specificity of iso- and methylisocitrate lyases, experimentally confirmed the possibility of direct evolution of methylisocitrate lyases from isocitrate lyases. Thus, gene duplication has increased the metabolic capacity of fungi.

Pyrenophora bromi, causal agent of brownspot of bromegrass, expresses a gene encoding a protein with homology and similar activity to Ptr ToxB, a host-selective toxin of wheat

Ptr ToxB, encoded by ToxB, is one of multiple host-selective toxins (HST) produced by the wheat pathogen Pyrenophora tritici-repentis. Homologs of ToxB are found in several ascomycetes, including sister species Pyrenophora bromi, causal agent of brownspot of bromegrass. Due to the close evolutionary relatedness of P. tritici-repentis and P. bromi and that of their grass hosts, we hypothesized that homologs of ToxB in P. bromi may act as HST in the disease interaction between P. bromi and bromegrass. A representative set of transcriptionally active P. bromi ToxB genes were heterologously expressed in Pichia pastoris and the resultant proteins tested for their ability to act as HST on bromegrass. The tested Pyrenophora bromi ToxB (Pb ToxB) proteins were not toxic to bromegrass; thus, Pb ToxB does not appear to function as an HST in the P. bromi-bromegrass interaction. Instead, we revealed that the Pb ToxB proteins can be toxic to Ptr ToxB-sensitive wheat, at levels similar to Ptr ToxB, and the corresponding P. bromi ToxB genes are expressed in P. bromi-inoculated wheat. Our data suggest that P. bromi possesses the potential to become a wheat pathogen and highlights the importance of investigating the interaction between P. bromi and wheat.

Evolution of transcription networks--lessons from yeasts

That regulatory evolution is important in generating phenotypic diversity was suggested soon after the discovery of gene regulation. In the past few decades, studies in animals have provided a number of examples in which phenotypic changes can be traced back to specific alterations in transcriptional regulation. Recent advances in DNA sequencing technology and functional genomics have stimulated a new wave of investigation in simple model organisms. In particular, several genome-wide comparative analyses of transcriptional circuits across different yeast species have been performed. These studies have revealed that transcription networks are remarkably plastic: large scale rewiring in which target genes move in and out of regulons through changes in cis-regulatory sequences appears to be a general phenomenon. Transcription factor substitution and the formation of new combinatorial interactions are also important contributors to the rewiring. In several cases, a transition through intermediates with redundant regulatory programs has been suggested as a mechanism through which rewiring can occur without a loss in fitness. Because the basic features of transcriptional regulation are deeply conserved, we speculate that large scale rewiring may underlie the evolution of complex phenotypes in multi-cellular organisms; if so, such rewiring may leave traceable changes in the genome from which the genetic basis of functional innovation can be detected.

The role of nucleosome positioning in the evolution of gene regulation

Chromatin organization plays a major role in gene regulation and can affect the function and evolution of new transcriptional programs. However, it can be difficult to decipher the basis of changes in chromatin organization and their functional effect on gene expression. Here, we present a large-scale comparative genomic analysis of the relationship between chromatin organization and gene expression, by measuring mRNA abundance and nucleosome positions genome-wide in 12 Hemiascomycota yeast species. We found substantial conservation of global and functional chromatin organization in all species, including prominent nucleosome-free regions (NFRs) at gene promoters, and distinct chromatin architecture in growth and stress genes. Chromatin organization has also substantially diverged in both global quantitative features, such as spacing between adjacent nucleosomes, and in functional groups of genes. Expression levels, intrinsic anti-nucleosomal sequences, and trans-acting chromatin modifiers all play important, complementary, and evolvable roles in determining NFRs. We identify five mechanisms that couple chromatin organization to evolution of gene regulation and have contributed to the evolution of respiro-fermentation and other key systems, including (1) compensatory evolution of alternative modifiers associated with conserved chromatin organization, (2) a gradual transition from constitutive to trans-regulated NFRs, (3) a loss of intrinsic anti-nucleosomal sequences accompanying changes in chromatin organization and gene expression, (4) re-positioning of motifs from NFRs to nucleosome-occluded regions, and (5) the expanded use of NFRs by paralogous activator-repressor pairs. Our study sheds light on the molecular basis of chromatin organization, and on the role of chromatin organization in the evolution of gene regulation.

Gene duplication and the evolution of ribosomal protein gene regulation in yeast

Coexpression of genes within a functional module can be conserved at great evolutionary distances, whereas the associated regulatory mechanisms can substantially diverge. For example, ribosomal protein (RP) genes are tightly coexpressed in Saccharomyces cerevisiae, but the cis and trans factors associated with them are surprisingly diverged across Ascomycota fungi. Little is known, however, about the functional impact of such changes on actual expression levels or about the selective pressures that affect them. Here, we address this question in the context of the evolution of the regulation of RP gene expression by using a comparative genomics approach together with cross-species functional assays. We show that an activator (Ifh1) and a repressor (Crf1) that control RP gene regulation in normal and stress conditions in S. cerevisiae are derived from the duplication and subsequent specialization of a single ancestral protein. We provide evidence that this regulatory innovation coincides with the duplication of RP genes in a whole-genome duplication (WGD) event and may have been important for tighter control of higher levels of RP transcripts. We find that subsequent loss of the derived repressor led to the loss of a stress-dependent repression of RPs in the fungal pathogen Candida glabrata. Our comparative computational and experimental approach shows how gene duplication can constrain and drive regulatory evolution and provides a general strategy for reconstructing the evolutionary trajectory of gene regulation across species.

Fungal regulatory evolution: cis and trans in the balance

Regulatory divergence is likely a major driving force in evolution. Comparative genomics is being increasingly used to infer the evolution of gene regulation. Ascomycota fungi are uniquely suited among eukaryotes for regulatory evolution studies, due to broad phylogenetic scope, many sequenced genomes, and tractability of genomic analysis. Here we review recent advances in the identification of the contribution of cis- and trans-factors to expression divergence. Whereas current strategies have led to the discovery of surprising signatures and mechanisms, we still understand very little about the adaptive role of regulatory evolution. Empirical studies including experimental evolution, comparative functional genomics and hybrid and engineered strains are showing early promise toward deciphering the contribution of regulatory divergence to adaptation.