Evolutionary tinkering with conserved components of a transcriptional regulatory network

Protein moonlighting by a target gene dominates phenotypic divergence of the Sef1 transcriptional regulatory network in yeasts

Transcriptional rewiring generates phenotypic novelty, acting as an important mechanism contributing to evolutionary development, speciation, and adaptation in all organisms. The phenotypic outcomes (functions) of transcription factor (TF) activity are determined by the combined effects of all target genes in the TF's regulatory network. Plastic rewiring of target genes accumulates during species divergence and ultimately alters phenotypes, indicating a TF functional switch. We define this phenomenon as 'disruptive rewiring', where the rewiring process disrupts the link between a TF and its original target genes that determine phenotypes. Here, we investigate if 'complete' disruptive rewiring is a prerequisite for a TF functional switch by employing chromatin immunoprecipitation sequencing, RNA expression, and phenotypic assays across yeast species. In yeasts where Sef1 targets TCA (tricarboxylic acid) cycle genes, we demonstrate that Sef1 orthologs can promote and inhibit respiratory growth by modulating the moonlighting function of their conserved target, NDE1. This modulation occurs without changing the overall association of Sef1 with TCA cycle genes. We propose that phenotypic masking by NDE1 promotes 'deceptive' disruptive rewiring of the Sef1 regulatory network in Saccharomyces cerevisiae, thereby potentially constraining future evolutionary trajectories.

Alternative sulphur metabolism in the fungal pathogen Candida parapsilosis

Candida parapsilosis is an opportunistic fungal pathogen commonly isolated from the environment and associated with nosocomial infection outbreaks worldwide. We describe here the construction of a large collection of gene disruptions, greatly increasing the molecular tools available for probing gene function in C. parapsilosis. We use these to identify transcription factors associated with multiple metabolic pathways, and in particular to dissect the network regulating the assimilation of sulphur. We find that, unlike in other yeasts and filamentous fungi, the transcription factor Met4 is not the main regulator of methionine synthesis. In C. parapsilosis, assimilation of inorganic sulphur (sulphate) and synthesis of cysteine and methionine is regulated by Met28, a paralog of Met4, whereas Met4 regulates expression of a wide array of transporters and enzymes involved in the assimilation of organosulfur compounds. Analysis of transcription factor binding sites suggests that Met4 is recruited by the DNA-binding protein Met32, and Met28 is recruited by Cbf1. Despite having different target genes, Met4 and Met28 have partial functional overlap, possibly because Met4 can contribute to assimilation of inorganic sulphur in the absence of Met28.

Complexes of HMO1 with DNA: Structure and Affinity

Saccharomyces cerevisiae HMO1 is an architectural nuclear DNA-binding protein that stimulates the activity of some remodelers and regulates the transcription of ribosomal protein genes, often binding to a DNA motif called IFHL. However, the molecular mechanism dictating this sequence specificity is unclear. Our circular dichroism spectroscopy studies show that the HMO1:DNA complex forms without noticeable changes in the structure of DNA and HMO1. Molecular modeling/molecular dynamics studies of the DNA complex with HMO1 Box B reveal two extended sites at the N-termini of helices I and II of Box B that are involved in the formation of the complex and stabilize the DNA bend induced by intercalation of the F114 side chain between base pairs. A comparison of the affinities of HMO1 for 24 bp DNA fragments containing either randomized or IFHL sequences reveals a twofold increase in the stability of the complex in the latter case, which may explain the selectivity in the recognition of the IFHL-containing promoter regions.

The Adr1 transcription factor directs regulation of the ergosterol pathway and azole resistance in Candida albicans

IMPORTANCE

Research often relies on well-studied orthologs within related species, with researchers using a well-studied gene or protein to allow prediction of the function of the ortholog. In the opportunistic pathogen Candida albicans, orthologs are usually compared with Saccharomyces cerevisiae, and this approach has been very fruitful. Many transcription factors (TFs) do similar jobs in the two species, but many do not, and typically changes in function are driven not by modifications in the structures of the TFs themselves but in the connections between the transcription factors and their regulated genes. This strategy of changing TF function has been termed transcription factor rewiring. In this study, we specifically looked for rewired transcription factors, or Candida-specific TFs, that might play a role in drug resistance. We investigated 30 transcription factors that were potentially rewired or were specific to the Candida clade. We found that the Adr1 transcription factor conferred resistance to drugs like fluconazole, amphotericin B, and terbinafine when activated. Adr1 is known for fatty acid and glycerol utilization in Saccharomyces, but our study reveals that it has been rewired and is connected to ergosterol biosynthesis in Candida albicans.

Candida albicans exhibits heterogeneous and adaptive cytoprotective responses to antifungal compounds

Candida albicans, an opportunistic human pathogen, poses a significant threat to human health and is associated with significant socio-economic burden. Current antifungal treatments fail, at least in part, because C. albicans can initiate a strong drug tolerance response that allows some cells to grow at drug concentrations above their minimal inhibitory concentration. To better characterize this cytoprotective tolerance program at the molecular single-cell level, we used a nanoliter droplet-based transcriptomics platform to profile thousands of individual fungal cells and establish their subpopulation characteristics in the absence and presence of antifungal drugs. Profiles of untreated cells exhibit heterogeneous expression that correlates with cell cycle stage with distinct metabolic and stress responses. At 2 days post-fluconazole exposure (a time when tolerance is measurable), surviving cells bifurcate into two major subpopulations: one characterized by the upregulation of genes encoding ribosomal proteins, rRNA processing machinery, and mitochondrial cellular respiration capacity, termed the Ribo-dominant (Rd) state; and the other enriched for genes encoding stress responses and related processes, termed the Stress-dominant (Sd) state. This bifurcation persists at 3 and 6 days post-treatment. We provide evidence that the ribosome assembly stress response (RASTR) is activated in these subpopulations and may facilitate cell survival.

SIR telomere silencing depends on nuclear envelope lipids and modulates sensitivity to a lysolipid

The nuclear envelope (NE) is important in maintaining genome organization. The role of lipids in communication between the NE and telomere regulation was investigated, including how changes in lipid composition impact gene expression and overall nuclear architecture. Yeast was treated with the non-metabolizable lysophosphatidylcholine analog edelfosine, known to accumulate at the perinuclear ER. Edelfosine induced NE deformation and disrupted telomere clustering but not anchoring. Additionally, the association of Sir4 at telomeres decreased. RNA-seq analysis showed altered expression of Sir-dependent genes located at sub-telomeric (0-10 kb) regions, consistent with Sir4 dispersion. Transcriptomic analysis revealed that two lipid metabolic circuits were activated in response to edelfosine, one mediated by the membrane sensing transcription factors, Spt23/Mga2, and the other by a transcriptional repressor, Opi1. Activation of these transcriptional programs resulted in higher levels of unsaturated fatty acids and the formation of nuclear lipid droplets. Interestingly, cells lacking Sir proteins displayed resistance to unsaturated-fatty acids and edelfosine, and this phenotype was connected to Rap1.

Erg6 Acts as a Downstream Effector of the Transcription Factor Flo8 To Regulate Biofilm Formation in Candida albicans

The yeast-to-hyphal morphotype transition and subsequent biofilm formation are important virulence factors of Candida albicans and are closely associated with ergosterol biosynthesis. Flo8 is an important transcription factor that determines filamentous growth and biofilm formation in C. albicans. However, the relationship between Flo8 and regulation of the ergosterol biosynthesis pathway remains elusive. Here, we analyzed the sterol composition of a flo8-deficient C. albicans strain by gas chromatography-mass spectrometry and observed the accumulation of the sterol intermediate zymosterol, the substrate of Erg6 (C-24 sterol methyltransferase). Accordingly, the transcription level of ERG6 was reduced in the flo8-deficient strain. Yeast one-hybrid experiments revealed that Flo8 physically interacted with the ERG6 promoter. Ectopic overexpression of ERG6 in the flo8-deficient strain partially restored biofilm formation and in vivo virulence in a Galleria mellonella infection model. These findings suggest that Erg6 is a downstream effector of the transcription factor Flo8 that mediates the cross talk between sterol synthesis and virulence factors in C. albicans. IMPORTANCE Biofilm formation by C. albicans hinders its eradication by immune cells and antifungal drugs. Flo8 is an important morphogenetic transcription factor that regulates the biofilm formation and in vivo virulence of C. albicans. However, little is known about how Flo8 regulates biofilm formation and fungal pathogenicity. Here, we determined that Flo8 directly binds to the promoter of ERG6 to positively regulate its transcriptional expression. Consistently, loss of flo8 results in the accumulation of the substrate of Erg6. Moreover, ectopic overexpression of ERG6 at least partially restores the biofilm formation and virulence of the flo8-deficient strain both in vitro and in vivo. This work provides a new perspective on the metabolic link between transcription factors and morphotypes in C. albicans.

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.

Differential Hsp90-dependent gene expression is strain-specific and common among yeast strains

Enhanced phenotypic diversity increases a population's likelihood of surviving catastrophic conditions. Hsp90, an essential molecular chaperone and a central network hub in eukaryotes, has been observed to suppress or enhance the effects of genetic variation on phenotypic diversity in response to environmental cues. Because many Hsp90-interacting genes are involved in signaling transduction pathways and transcriptional regulation, we tested how common Hsp90-dependent differential gene expression is in natural populations. Many genes exhibited Hsp90-dependent strain-specific differential expression in five diverse yeast strains. We further identified transcription factors (TFs) potentially contributing to variable expression. We found that on Hsp90 inhibition or environmental stress, activities or abundances of Hsp90-dependent TFs varied among strains, resulting in differential strain-specific expression of their target genes, which consequently led to phenotypic diversity. We provide evidence that individual strains can readily display specific Hsp90-dependent gene expression, suggesting that the evolutionary impacts of Hsp90 are widespread in nature.

Evolvability and Macroevolution: Overview and Synthesis

Evolvability is best addressed from a multi-level, macroevolutionary perspective through a comparative approach that tests for among-clade differences in phenotypic diversification in response to an opportunity, such as encountered after a mass extinction, entering a new adaptive zone, or entering a new geographic area. Analyzing the dynamics of clades under similar environmental conditions can (partially) factor out shared external drivers to recognize intrinsic differences in evolvability, aiming for a macroevolutionary analog of a common-garden experiment. Analyses will be most powerful when integrating neontological and paleontological data: determining differences among extant populations that can be hypothesized to generate large-scale, long-term contrasts in evolvability among clades; or observing large-scale differences among clade histories that can by hypothesized to reflect contrasts in genetics and development observed directly in extant populations. However, many comparative analyses can be informative on their own, as explored in this overview. Differences in clade-level evolvability can be visualized in diversity-disparity plots, which can quantify positive and negative departures of phenotypic productivity from stochastic expectations scaled to taxonomic diversification. Factors that evidently can promote evolvability include modularity—when selection aligns with modular structure or with morphological integration patterns; pronounced ontogenetic changes in morphology, as in allometry or multiphase life cycles; genome size; and a variety of evolutionary novelties, which can also be evaluated using macroevolutionary lags between the acquisition of a trait and phenotypic diversification, and dead-clade-walking patterns that may signal a loss of evolvability when extrinsic factors can be excluded. High speciation rates may indirectly foster phenotypic evolvability, and vice versa. Mechanisms are controversial, but clade evolvability may be higher in the Cambrian, and possibly early in the history of clades at other times; in the tropics; and, for marine organisms, in shallow-water disturbed habitats.

SAGA Complex Subunits in Candida albicans Differentially Regulate Filamentation, Invasiveness, and Biofilm Formation

SAGA (Spt-Ada-Gcn5-acetyltransferase) is a highly conserved, multiprotein co-activator complex that consists of five distinct modules. It has two enzymatic functions, a histone acetyltransferase (HAT) and a deubiquitinase (DUB) and plays a central role in processes such as transcription initiation, elongation, protein stability, and telomere maintenance. We analyzed conditional and null mutants of the SAGA complex module components in the fungal pathogen Candida albicans; Ngg1, (the HAT module); Ubp8, (the DUB module); Tra1, (the recruitment module), Spt7, (the architecture module) and Spt8, (the TBP interaction unit), and assessed their roles in a variety of cellular processes. We observed that spt7Δ/Δ and spt8Δ/Δ strains have a filamentous phenotype, and both are highly invasive in yeast growing conditions as compared to the wild type, while ngg1Δ/Δ and ubp8Δ/Δ are in yeast-locked state and non-invasive in both YPD media and filamentous induced conditions compared to wild type. RNA-sequencing-based transcriptional profiling of SAGA mutants reveals upregulation of hyphal specific genes in spt7Δ/Δ and spt8Δ/Δ strains and downregulation of ergosterol metabolism pathway. As well, spt7Δ/Δ and spt8Δ/Δ confer susceptibility to antifungal drugs, to acidic and alkaline pH, to high temperature, and to osmotic, oxidative, cell wall, and DNA damage stresses, indicating that these proteins are important for genotoxic and cellular stress responses. Despite having similar morphological phenotypes (constitutively filamentous and invasive) spt7 and spt8 mutants displayed variation in nuclear distribution where spt7Δ/Δ cells were frequently binucleate and spt8Δ/Δ cells were consistently mononucleate. We also observed that spt7Δ/Δ and spt8Δ/Δ mutants were quickly engulfed by macrophages compared to ngg1Δ/Δ and ubp8Δ/Δ strains. All these findings suggest that the SAGA complex modules can have contrasting functions where loss of Spt7 or Spt8 enhances filamentation and invasiveness while loss of Ngg1 or Ubp8 blocks these processes.

Modulation of the complex regulatory network for methionine biosynthesis in fungi

The assimilation of inorganic sulfate and the synthesis of the sulfur-containing amino acids methionine and cysteine is mediated by a multibranched biosynthetic pathway. We have investigated this circuitry in the fungal pathogen Candida albicans, which is phylogenetically intermediate between the filamentous fungi and Saccharomyces cerevisiae. In S. cerevisiae, this pathway is regulated by a collection of five transcription factors (Met4, Cbf1, Met28, and Met31/Met32), while in the filamentous fungi the pathway is controlled by a single Met4-like factor. We found that in C. albicans, the Met4 ortholog is also a core regulator of methionine biosynthesis, where it functions together with Cbf1. While C. albicans encodes this Met4 protein, a Met4 paralog designated Met28 (Orf19.7046), and a Met31 protein, deletion, and activation constructs suggest that of these proteins only Met4 is actually involved in the regulation of methionine biosynthesis. Both Met28 and Met31 are linked to other functions; Met28 appears essential, and Met32 appears implicated in the regulation of genes of central metabolism. Therefore, while S. cerevisiae and C. albicans share Cbf1 and Met4 as central elements of the methionine biosynthesis control, the other proteins that make up the circuit in S. cerevisiae are not members of the C. albicans control network, and so the S. cerevisiae circuit likely represents a recently evolved arrangement.

Plastic rewiring of Sef1 transcriptional networks and the potential of non-functional transcription factor binding in facilitating adaptive evolution

Prior and extensive plastic rewiring of a transcriptional network, followed by a functional switch of the conserved transcriptional regulator, can shape the evolution of a new network with diverged functions. The presence of three distinct iron regulatory systems in fungi that use orthologous transcriptional regulators suggests that these systems evolved in that manner. Orthologs of the transcriptional activator Sef1 are believed to be central to how iron regulatory systems developed in fungi, involving gene gain, plastic network rewiring, and switches in regulatory function. We show that, in the protoploid yeast Lachancea kluyveri, plastic rewiring of the L. kluyveri Sef1 (Lk-Sef1) network, together with a functional switch, enabled Lk-Sef1 to regulate TCA cycle genes, unlike Candida albicans Sef1 that mainly regulates iron-uptake genes. Moreover, we observed pervasive nonfunctional binding of Sef1 to its target genes. Enhancing Lk-Sef1 activity resuscitated the corresponding transcriptional network, providing immediate adaptive benefits in changing environments. Our study not only sheds light on the evolution of Sef1-centered transcriptional networks but also shows the adaptive potential of nonfunctional transcription factor binding for evolving phenotypic novelty and diversity.

Network motif-based analysis of regulatory patterns in paralogous gene pairs

Current high-throughput experimental techniques make it feasible to infer gene regulatory interactions at the whole-genome level with reasonably good accuracy. Such experimentally inferred regulatory networks have become available for a number of simpler model organisms such as S. cerevisiae, and others. The availability of such networks provides an opportunity to compare gene regulatory processes at the whole genome level, and in particular, to assess similarity of regulatory interactions for homologous gene pairs either from the same or from different species. We present here a new technique for analyzing the regulatory interaction neighborhoods of paralogous gene pairs. Our central focus is the analysis of S. cerevisiae gene interaction graphs, which are of particular interest due to the ancestral whole-genome duplication (WGD) that allows to distinguish between paralogous transcription factors that are traceable to this duplication event and other paralogues. Similar analysis is also applied to E. coli and C. elegans networks. We compare paralogous gene pairs according to the presence and size of bi-fan arrays, classically associated in the literature with gene duplication, within other network motifs. We further extend this framework beyond transcription factor comparison to obtain topology-based similarity metrics based on the overlap of interaction neighborhoods applicable to most genes in a given organism. We observe that our network divergence metrics show considerably larger similarity between paralogues, especially those traceable to WGD. This is the case for both yeast and C. elegans, but not for E. coli regulatory network. While there is no obvious cross-species link between metrics, different classes of paralogues show notable differences in interaction overlap, with traceable duplications tending toward higher overlap compared to genes with shared protein families. Our findings indicate that divergence in paralogous interaction networks reflects a shared genetic origin, and that our approach may be useful for investigating structural similarity in the interaction networks of paralogous genes.

Genic introgression from an invasive exotic fungal forest pathogen increases the establishment potential of a sibling native pathogen

Significant hybridization between the invasive North American fungal plant pathogen Heterobasidion irregulare and its Eurasian sister species H. annosum is ongoing in Italy. Whole genomes of nine natural hybrids were sequenced, assembled and compared with those of three genotypes each of the two parental species. Genetic relationships among hybrids and their level of admixture were determined. A multi-approach pipeline was used to assign introgressed genomic blocks to each of the two species. Alleles that introgressed from H. irregulare to H. annosum were associated with pathways putatively related to saprobic processes, while alleles that introgressed from the native to the invasive species were mainly linked to gene regulation. There was no overlap of allele categories introgressed in the two directions. Phenotypic experiments documented a fitness increase in H. annosum genotypes characterized by introgression of alleles from the invasive species, supporting the hypothesis that hybridization results in putatively adaptive introgression. Conversely, introgression from the native into the exotic species appeared to be driven by selection on genes favoring genome stability. Since the introgression of specific alleles from the exotic H. irregulare into the native H. annosum increased the invasiveness of the latter species, we propose that two invasions may be co-occurring: the first one by genotypes of the exotic species, and the second one by alleles belonging to the exotic species. Given that H. irregulare represents a threat to European forests, monitoring programs need to track not only exotic genotypes in native forest stands, but also exotic alleles introgressed in native genotypes.

Physiological function of FKBP12, a primary target of rapamycin/FK506: a newly identified role in transcription of ribosomal protein genes in yeast

In this review, we have summarized the information from a study on FKBP12 (FK506 binding protein 12 kDa) with a view to understand its drug-free, physiological roles in transcription of ribosomal protein gene in Saccharomyces cerevisiae. FKBP12 with peptidyl-prolylisomerase (PPIase) activity is widely conserved among many eukaryotes. FKBP12 is a primary target for the two structurally related drugs, FK506 and rapamycin. FKBP12 bound with FK506 or rapamycin inhibits calcineurin and target of rapamycin complex 1 (TORC1), respectively. The molecular mechanisms of the effect of FKBP12 in the presence of these drugs have been elucidated. Conversely, the physiological role of FKBP12 has been unclear, especially in yeast. Our study revealed that the deletion of FPR1 (FK506-sensitive prolinerotamase 1 gene), which encodes yeast FKBP12, induced severe growth defect synthetically with deletion of HMO1 (high mobility group family 1). HMO1 encodes an HMGB family protein involved in transcription of ribosomal component genes. Fpr1 was shown to bind specifically to the promoters of ribosomal protein genes (RPGs) dependent on Rap1 (repressor/activator binding protein 1). Importantly, Fpr1 and Hmo1 promote the binding of Fhl1/Ifh1 (forkhead-like 1/interacts with forkhead 1), key regulators of RPG transcription, to certain RPG promoters independently and/or cooperatively with each other. Taken together, we conclude that Fpr1 physiologically functions as transcription factor of RPGs in S. cerevisiae. To our knowledge, this is the first study to demonstrate that FKBP12 participates in ribosome synthesis independently of drugs, and it may also provide a clue to the unidentified function of other PPIase proteins.

The RSC (Remodels the Structure of Chromatin) complex of Candida albicans shows compositional divergence with distinct roles in regulating pathogenic traits

Regulation of gene expression programs is crucial for the survival of microbial pathogens in host environments and for their ability to cause disease. Here we investigated the epigenetic regulator RSC (Remodels the Structure of Chromatin) in the most prevalent human fungal pathogen Candida albicans. Biochemical analysis showed that CaRSC comprises 13 subunits and contains two novel non-essential members, which we named Nri1 and Nri2 (Novel RSC Interactors) that are exclusive to the CTG clade of Saccharomycotina. Genetic analysis showed distinct essentiality of C. albicans RSC subunits compared to model fungal species suggesting functional and structural divergence of RSC functions in this fungal pathogen. Transcriptomic and proteomic profiling of a conditional mutant of the essential catalytic subunit gene STH1 demonstrated global roles of RSC in C. albicans biology, with the majority of growth-related processes affected, as well as mis-regulation of genes involved in morphotype switching, host-pathogen interaction and adaptive fitness. We further assessed the functions of non-essential CaRSC subunits, showing that the novel subunit Nri1 and the bromodomain subunit Rsc4 play roles in filamentation and stress responses; and also interacted at the genetic level to regulate cell viability. Consistent with these roles, Rsc4 is required for full virulence of C. albicans in the murine model of systemic infection. Taken together, our data builds the first comprehensive study of the composition and roles of RSC in C. albicans, showing both conserved and distinct features compared to model fungal systems. The study illuminates how C. albicans uses RSC-dependent transcriptional regulation to respond to environmental signals and drive survival fitness and virulence in mammals.

Regulation of ribosomal protein genes: An ordered anarchy

Ribosomal protein genes are among the most highly expressed genes in most cell types. Their products are generally essential for ribosome synthesis, which is the cornerstone for cell growth and proliferation. Many cellular resources are dedicated to producing ribosomal proteins and thus this process needs to be regulated in ways that carefully balance the supply of nascent ribosomal proteins with the demand for new ribosomes. Ribosomal protein genes have classically been viewed as a uniform interconnected regulon regulated in eukaryotic cells by target of rapamycin and protein kinase A pathway in response to changes in growth conditions and/or cellular status. However, recent literature depicts a more complex picture in which the amount of ribosomal proteins produced varies between genes in response to two overlapping regulatory circuits. The first includes the classical general ribosome‐producing program and the second is a gene‐specific feature responsible for fine‐tuning the amount of ribosomal proteins produced from each individual ribosomal gene. Unlike the general pathway that is mainly controlled at the level of transcription and translation, this specific regulation of ribosomal protein genes is largely achieved through changes in pre‐mRNA splicing efficiency and mRNA stability. By combining general and specific regulation, the cell can coordinate ribosome production, while allowing functional specialization and diversity. Here we review the many ways ribosomal protein genes are regulated, with special focus on the emerging role of posttranscriptional regulatory events in fine‐tuning the expression of ribosomal protein genes and its role in controlling the potential variation in ribosome functions.

This article is categorized under:

Translation > Ribosome Biogenesis

Translation > Ribosome Structure/Function

Translation > Translation Regulation

A Theoretical Framework for Evolutionary Cell Biology

One of the last uncharted territories in evolutionary biology concerns the link with cell biology. Because all phenotypes ultimately derive from events at the cellular level, this connection is essential to building a mechanism-based theory of evolution. Given the impressive developments in cell biological methodologies at the structural and functional levels, the potential for rapid progress is great. The primary challenge for theory development is the establishment of a quantitative framework that transcends species boundaries. Two approaches to the problem are presented here: establishing the long-term steady-state distribution of mean phenotypes under specific regimes of mutation, selection, and drift and evaluating the energetic costs of cellular structures and functions. Although not meant to be the final word, these theoretical platforms harbor potential for generating insight into a diversity of unsolved problems, ranging from genome structure to cellular architecture to aspects of motility in organisms across the Tree of Life.

Protein-Protein Interactions in Candida albicans

Despite being one of the most important human fungal pathogens, Candida albicans has not been studied extensively at the level of protein-protein interactions (PPIs) and data on PPIs are not readily available in online databases. In January 2018, the database called "Biological General Repository for Interaction Datasets (BioGRID)" that contains the most PPIs for C. albicans, only documented 188 physical or direct PPIs (release 3.4.156) while several more can be found in the literature. Other databases such as the String database, the Molecular INTeraction Database (MINT), and the Database for Interacting Proteins (DIP) database contain even fewer interactions or do not even include C. albicans as a searchable term. Because of the non-canonical codon usage of C. albicans where CUG is translated as serine rather than leucine, it is often problematic to use the yeast two-hybrid system in Saccharomyces cerevisiae to study C. albicans PPIs. However, studying PPIs is crucial to gain a thorough understanding of the function of proteins, biological processes and pathways. PPIs can also be potential drug targets. To aid in creating PPI networks and updating the BioGRID, we performed an exhaustive literature search in order to provide, in an accessible format, a more extensive list of known PPIs in C. albicans.

Functional duplication of Rap1 in methylotrophic yeasts

The telomere regulator and transcription factor Rap1 is the only telomere protein conserved in yeasts and mammals. Its functional repertoire in budding yeasts is a particularly interesting field for investigation, given the high evolutionary diversity of this group of unicellular organisms. In the methylotrophic thermotolerant species Hansenula polymorpha DL-1 the RAP1 gene is duplicated (HpRAP1A and HpRAP1B). Here, we report the functional characterization of the two paralogues from H. polymorpha DL-1. We uncover distinct (but overlapping) DNA binding preferences of HpRap1A and HpRap1B proteins. We show that only HpRap1B is able to recognize telomeric DNA directly and to protect it from excessive recombination, whereas HpRap1A is associated with subtelomere regions. Furthermore, we identify specific binding sites for both HpRap1A and HpRap1B within promoters of a large number of ribosomal protein genes (RPGs), implicating Rap1 in the control of the RP regulon in H. polymorpha. Our bioinformatic analysis suggests that RAP1 was duplicated early in the evolution of the “methylotrophs” clade, and the two genes evolved independently. Therefore, our characterization of Rap1 paralogues in H. polymorpha may be relevant to other “methylotrophs”, yielding valuable insights into the evolution of budding yeasts.

The p38/HOG stress-activated protein kinase network couples growth to division in Candida albicans

Cell size is a complex trait that responds to developmental and environmental cues. Quantitative size analysis of mutant strain collections disrupted for protein kinases and transcriptional regulators in the pathogenic yeast Candida albicans uncovered 66 genes that altered cell size, few of which overlapped with known size genes in the budding yeast Saccharomyces cerevisiae. A potent size regulator specific to C. albicans was the conserved p38/HOG MAPK module that mediates the osmostress response. Basal HOG activity inhibited the SBF G1/S transcription factor complex in a stress-independent fashion to delay the G1/S transition. The HOG network also governed ribosome biogenesis through the master transcriptional regulator Sfp1. Hog1 bound to the promoters and cognate transcription factors for ribosome biogenesis regulons and interacted genetically with the SBF G1/S machinery, and thereby directly linked cell growth and division. These results illuminate the evolutionary plasticity of size control and identify the HOG module as a nexus of cell cycle and growth regulation.

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.

Integration of Growth and Cell Size via the TOR Pathway and the Dot6 Transcription Factor in Candida albicans

In most species, size homeostasis appears to be exerted in late G1 phase as cells commit to division, called Start in yeast and the Restriction Point in metazoans. This size threshold couples cell growth to division, and, thereby, establishes long-term size homeostasis. Our former investigations have shown that hundreds of genes markedly altered cell size under homeostatic growth conditions in the opportunistic yeast Candida albicans, but surprisingly only few of these overlapped with size control genes in the budding yeast Saccharomyces cerevisiae Here, we investigated one of the divergent potent size regulators in C. albicans, the Myb-like HTH transcription factor Dot6. Our data demonstrated that Dot6 is a negative regulator of Start, and also acts as a transcriptional activator of ribosome biogenesis (Ribi) genes. Genetic epistasis uncovered that Dot6 interacted with the master transcriptional regulator of the G1 machinery, SBF complex, but not with the Ribi and cell size regulators Sch9, Sfp1, and p38/Hog1. Dot6 was required for carbon-source modulation of cell size, and it is regulated at the level of nuclear localization by the TOR pathway. Our findings support a model where Dot6 acts as a hub that integrates growth cues directly via the TOR pathway to control the commitment to mitotic division at G1.

The Genetic Makeup and Expression of the Glycolytic and Fermentative Pathways Are Highly Conserved Within the Saccharomyces Genus

The ability of the yeast Saccharomyces cerevisiae to convert glucose, even in the presence of oxygen, via glycolysis and the fermentative pathway to ethanol has played an important role in its domestication. Despite the extensive knowledge on these pathways in S. cerevisiae, relatively little is known about their genetic makeup in other industrially relevant Saccharomyces yeast species. In this study we explore the diversity of the glycolytic and fermentative pathways within the Saccharomyces genus using S. cerevisiae, Saccharomyces kudriavzevii, and Saccharomyces eubayanus as paradigms. Sequencing data revealed a highly conserved genetic makeup of the glycolytic and fermentative pathways in the three species in terms of number of paralogous genes. Although promoter regions were less conserved between the three species as compared to coding sequences, binding sites for Rap1, Gcr1 and Abf1, main transcriptional regulators of glycolytic and fermentative genes, were highly conserved. Transcriptome profiling of these three strains grown in aerobic batch cultivation in chemically defined medium with glucose as carbon source, revealed a remarkably similar expression of the glycolytic and fermentative genes across species, and the conserved classification of genes into major and minor paralogs. Furthermore, transplantation of the promoters of major paralogs of S. kudriavzevii and S. eubayanus into S. cerevisiae demonstrated not only the transferability of these promoters, but also the similarity of their strength and response to various environmental stimuli. The relatively low homology of S. kudriavzevii and S. eubayanus promoters to their S. cerevisiae relatives makes them very attractive alternatives for strain construction in S. cerevisiae, thereby expanding the S. cerevisiae molecular toolbox.

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.

Host-Pathogen Interactions Mediated by MDR Transporters in Fungi: As Pleiotropic as it Gets!

Fungal infections caused by Candida, Aspergillus, and Cryptococcus species are an increasing problem worldwide, associated with very high mortality rates. The successful prevalence of these human pathogens is due to their ability to thrive in stressful host niche colonization sites, to tolerate host immune system-induced stress, and to resist antifungal drugs. This review focuses on the key role played by multidrug resistance (MDR) transporters, belonging to the ATP-binding cassette (ABC), and the major facilitator superfamilies (MFS), in mediating fungal resistance to pathogenesis-related stresses. These clearly include the extrusion of antifungal drugs, with C. albicans CDR1 and MDR1 genes, and corresponding homologs in other fungal pathogens, playing a key role in this phenomenon. More recently, however, clues on the transcriptional regulation and physiological roles of MDR transporters, including the transport of lipids, ions, and small metabolites, have emerged, linking these transporters to important pathogenesis features, such as resistance to host niche environments, biofilm formation, immune system evasion, and virulence. The wider view of the activity of MDR transporters provided in this review highlights their relevance beyond drug resistance and the need to develop therapeutic strategies that successfully face the challenges posed by the pleiotropic nature of these transporters.

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.

How transcription circuits explore alternative architectures while maintaining overall circuit output

This review by Dalal and Johnson focuses on the evolutionary rewiring of transcription regulators and the conservation of patterns of gene expression. They describe how preservation of gene expression patterns in the wake of extensive rewiring is a general feature of transcription circuit evolution.

Transcription regulators bind to cis-regulatory sequences and thereby control the expression of target genes. While transcription regulators and the target genes that they regulate are often deeply conserved across species, the connections between the two change extensively over evolutionary timescales. In this review, we discuss case studies where, despite this extensive evolutionary rewiring, the resulting patterns of gene expression are preserved. We also discuss in silico models that reach the same general conclusions and provide additional insights into how this process occurs. Together, these approaches make a strong case that the preservation of gene expression patterns in the wake of extensive rewiring is a general feature of transcription circuit evolution.

Abf1 and other general regulatory factors control ribosome biogenesis gene expression in budding yeast

Ribosome biogenesis in Saccharomyces cerevisiae involves a regulon of >200 genes (Ribi genes) coordinately regulated in response to nutrient availability and cellular growth rate. Two cis-acting elements called PAC and RRPE are known to mediate Ribi gene repression in response to nutritional downshift. Here, we show that most Ribi gene promoters also contain binding sites for one or more General Regulatory Factors (GRFs), most frequently Abf1 and Reb1, and that these factors are enriched in vivo at Ribi promoters. Abf1/Reb1/Tbf1 promoter association was required for full Ribi gene expression in rich medium and for its modulation in response to glucose starvation, characterized by a rapid drop followed by slow recovery. Such a response did not entail changes in Abf1 occupancy, but it was paralleled by a quick increase, followed by slow decrease, in Rpd3L histone deacetylase occupancy. Remarkably, Abf1 site disruption also abolished Rpd3L complex recruitment in response to starvation. Extensive mutational analysis of the DBP7 promoter revealed a complex interplay of Tbf1 sites, PAC and RRPE in the transcriptional regulation of this Ribi gene. Our observations point to GRFs as new multifaceted players in Ribi gene regulation both during exponential growth and under repressive conditions.

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.

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

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

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.

Functional Characterization of Populus PsnSHN2 in Coordinated Regulation of Secondary Wall Components in Tobacco

Wood formation is a biological process during which the most abundant lignocellulosic biomass on earth is produced. Although a number of transcription factors have been linked to the regulation of wood formation process, none of them has been demonstrated to be a higher hierarchical regulator that coordinately regulates secondary wall biosynthesis genes. Here, we identified a Populus gene, PsnSHN2, a counterpart of the Arabidopsis AP2/ERF type transcription factor, SHINE2. PsnSHN2 is predominantly expressed in xylem tissues and acted evidently as a high hierarchical transcriptional activator. Overexpression of PsnSHN2 in tobacco significantly altered the expression of both transcription factors and biosynthesis genes involved in secondary wall formation, leading to the thickened secondary walls and the changed cell wall composition. The most significant changes occurred in the contents of cellulose and hemicellulose that increased 37% and 28%, respectively, whereas the content of lignin that decreased 34%. Furthermore, PsnSHN2 activated or repressed the promoter activities of transcription factors involved in secondary wall biosynthesis and bound to five cis-acting elements enriched in the promoter regions of these transcription factors. Taken together, our results suggest PsnSHN2 coordinately regulate secondary wall formation through selective up/down-regulation of its downstream transcription factors that control secondary wall formation.

Genome-Wide Screen for Haploinsufficient Cell Size Genes in the Opportunistic Yeast Candida albicans

One of the most critical but still poorly understood aspects of eukaryotic cell proliferation is the basis for commitment to cell division in late G1 phase, called Start in yeast and the Restriction Point in metazoans. In all species, a critical cell size threshold coordinates cell growth with cell division and thereby establishes a homeostatic cell size. While a comprehensive survey of cell size genetic determinism has been performed in the saprophytic yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe, very little is known in pathogenic fungi. As a number of critical Start regulators are haploinsufficient for cell size, we applied a quantitative analysis of the size phenome, using elutriation-barcode sequencing methodology, to 5639 barcoded heterozygous deletion strains of the opportunistic yeast Candida albicans. Our screen identified conserved known regulators and biological processes required to maintain size homeostasis in the opportunistic yeast C. albicans. We also identified novel C. albicans-specific size genes and provided a conceptual framework for future mechanistic studies. Interestingly, some of the size genes identified were required for fungal pathogenicity suggesting that cell size homeostasis may be elemental to C. albicans fitness or virulence inside the host.

Metabolic adaptation via regulated enzyme degradation in the pathogenic yeast Candida albicans

The virulence of Candida albicans is dependent upon fitness attributes as well as virulence factors. These attributes include robust stress responses and metabolic flexibility. The assimilation of carbon sources is important for growth and essential for the establishment of infections by C. albicans. Previous studies showed that the C. albicans ICL1 genes, which encode the glyoxylate cycle enzymes isocitratelyase are required for growth on non-fermentable carbon sources such as lactate and oleic acid and were repressed by 2% glucose. In contrast to S. cerevsiae, the enzyme CaIcl1 was not destabilised by glucose, resulting with its metabolite remaining at high levels. Further glucose addition has caused CaIcl1 to lose its signal and mechanisms that trigger destabilization in response to glucose. Another purpose of this study was to test the stability of the Icl1 enzyme in response to the dietary sugars, fructose, and galactose. In the present study, the ICL1 mRNAs expression was quantified using Quantitative Real Time PCR, whereby the stability of protein was measured and quantified using Western blot and phosphoimager, and the replacing and cloning of ICL1 ORF by gene recombination and ubiquitin binding was conducted via co-immuno-precipitation. Following an analogous experimental approach, the analysis was repeated using S. cerevisiaeas a control. Both galactose and fructose were found to trigger the degradation of the ICL1 transcript in C. albicans. The Icl1 enzyme was stable following galactose addition but was degraded in response to fructose. C. albicans Icl1 (CaIcl1) was also subjected to fructose-accelerated degradation when expressed in S. cerevisiae, indicating that, although it lacks a ubiquitination site, CaIcl1 is sensitive to fructose-accelerated protein degradation. The addition of an ubiquitination site to CaIcl1 resulted in this enzyme becoming sensitive to galactose-accelerated degradation and increases its rate of degradation in the presence of fructose. It can be concluded that ubiquitin-independent pathways of fructose-accelerated enzyme degradation exist in C. albicans.

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.

Beauvericin Potentiates Azole Activity via Inhibition of Multidrug Efflux, Blocks Candida albicans Morphogenesis, and Is Effluxed via Yor1 and Circuitry Controlled by Zcf29

Invasive fungal infections are a leading cause of human mortality. Effective treatment is hindered by the rapid emergence of resistance to the limited number of antifungal drugs, demanding new strategies to treat life-threatening fungal infections. Here, we explore a powerful strategy to enhance antifungal efficacy against leading human fungal pathogens by using the natural product beauvericin. We found that beauvericin potentiates the activity of azole antifungals against azole-resistant Candida isolates via inhibition of multidrug efflux and that beauvericin itself is effluxed via Yor1. As observed in Saccharomyces cerevisiae, we determined that beauvericin inhibits TOR signaling in Candida albicans To further characterize beauvericin activity in C. albicans, we leveraged genome sequencing of beauvericin-resistant mutants. Resistance was conferred by mutations in transcription factor genes TAC1, a key regulator of multidrug efflux, and ZCF29, which was uncharacterized. Transcriptional profiling and chromatin immunoprecipitation coupled to microarray analyses revealed that Zcf29 binds to and regulates the expression of multidrug transporter genes. Beyond drug resistance, we also discovered that beauvericin blocks the C. albicans morphogenetic transition from yeast to filamentous growth in response to diverse cues. We found that beauvericin represses the expression of many filament-specific genes, including the transcription factor BRG1 Thus, we illuminate novel circuitry regulating multidrug efflux and establish that simultaneously targeting drug resistance and morphogenesis provides a promising strategy to combat life-threatening fungal infections.

Dual-species transcriptional profiling during systemic candidiasis reveals organ-specific host-pathogen interactions

Candida albicans is a common cause of life-threatening fungal bloodstream infections. In the murine model of systemic candidiasis, the kidney is the primary target organ while the fungal load declines over time in liver and spleen. To better understand these organ-specific differences in host-pathogen interaction, we performed gene expression profiling of murine kidney, liver and spleen and determined the fungal transcriptome in liver and kidney. We observed a delayed transcriptional immune response accompanied by late induction of fungal stress response genes in the kidneys. In contrast, early upregulation of the proinflammatory response in the liver was associated with a fungal transcriptome resembling response to phagocytosis, suggesting that phagocytes contribute significantly to fungal control in the liver. Notably, C. albicans hypha-associated genes were upregulated in the absence of visible filamentation in the liver, indicating an uncoupling of gene expression and morphology and a morphology-independent effect by hypha-associated genes in this organ. Consistently, integration of host and pathogen transcriptional data in an inter-species gene regulatory network indicated connections of C. albicans cell wall remodelling and metabolism to the organ-specific immune responses.

Rewiring of the Ppr1 Zinc Cluster Transcription Factor from Purine Catabolism to Pyrimidine Biogenesis in the Saccharomycetaceae

Metabolic pathways are largely conserved in eukaryotes, but the transcriptional regulation of these pathways can sometimes vary between species; this has been termed "rewiring." Recently, it has been established that in the Saccharomyces lineage starting from Naumovozyma castellii, genes involved in allantoin breakdown have been genomically relocated to form the DAL cluster. The formation of the DAL cluster occurred along with the loss of urate permease (UAP) and urate oxidase (UOX), reducing the requirement for oxygen and bypassing the candidate Ppr1 inducer, uric acid. In Saccharomyces cerevisiae, this allantoin catabolism cluster is regulated by the transcription factor Dal82, which is not present in many of the pre-rearrangement fungal species. We have used ChIP-chip analysis, transcriptional profiling of an activated Ppr1 protein, bioinformatics, and nitrogen utilization studies to establish that in Candida albicans the zinc cluster transcription factor Ppr1 controls this allantoin catabolism regulon. Intriguingly, in S. cerevisiae, the Ppr1 ortholog binds the same DNA motif (CGG(N6)CCG) as in C. albicans but serves as a regulator of pyrimidine biosynthesis. This transcription factor rewiring appears to have taken place at the same phylogenetic step as the formation of the rearranged DAL cluster. This transfer of the control of allantoin degradation from Ppr1 to Dal82, together with the repositioning of Ppr1 to the regulation of pyrimidine biosynthesis, may have resulted from a switch to a metabolism that could exploit hypoxic conditions in the lineage leading to N. castellii and S. cerevisiae.

A Network of Paralogous Stress Response Transcription Factors in the Human Pathogen Candida glabrata

The yeast Candida glabrata has become the second cause of systemic candidemia in humans. However, relatively few genome-wide studies have been conducted in this organism and our knowledge of its transcriptional regulatory network is quite limited. In the present work, we combined genome-wide chromatin immunoprecipitation (ChIP-seq), transcriptome analyses, and DNA binding motif predictions to describe the regulatory interactions of the seven Yap (Yeast AP1) transcription factors of C. glabrata. We described a transcriptional network containing 255 regulatory interactions and 309 potential target genes. We predicted with high confidence the preferred DNA binding sites for 5 of the 7 CgYaps and showed a strong conservation of the Yap DNA binding properties between S. cerevisiae and C. glabrata. We provided reliable functional annotation for 3 of the 7 Yaps and identified for Yap1 and Yap5 a core regulon which is conserved in S. cerevisiae, C. glabrata, and C. albicans. We uncovered new roles for CgYap7 in the regulation of iron-sulfur cluster biogenesis, for CgYap1 in the regulation of heme biosynthesis and for CgYap5 in the repression of GRX4 in response to iron starvation. These transcription factors define an interconnected transcriptional network at the cross-roads between redox homeostasis, oxygen consumption, and iron metabolism.

Making sense of transcription networks

When transcription regulatory networks are compared among distantly related eukaryotes, a number of striking similarities are observed: a larger-than-expected number of genes, extensive overlapping connections, and an apparently high degree of functional redundancy. It is often assumed that the complexity of these networks represents optimized solutions, precisely sculpted by natural selection; their common features are often asserted to be adaptive. Here, we discuss support for an alternative hypothesis: the common structural features of transcription networks arise from evolutionary trajectories of "least resistance"--that is, the relative ease with which certain types of network structures are formed during their evolution.

Promoter architecture and transcriptional regulation of Abf1-dependent ribosomal protein genes in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, ribosomal protein gene (RPG) promoters display binding sites for either Rap1 or Abf1 transcription factors. Unlike Rap1-associated promoters, the small cohort of Abf1-dependent RPGs (Abf1-RPGs) has not been extensively investigated. We show that RPL3, RPL4B, RPP1A, RPS22B and RPS28A/B share a common promoter architecture, with an Abf1 site upstream of a conserved element matching the sequence recognized by Fhl1, a transcription factor which together with Ifh1 orchestrates Rap1-associated RPG regulation. Abf1 and Fhl1 promoter association was confirmed by ChIP and/or gel retardation assays. Mutational analysis revealed a more severe requirement of Abf1 than Fhl1 binding sites for RPG transcription. In the case of RPS22B an unusual Tbf1 binding site promoted both RPS22B and intron-hosted SNR44 expression. Abf1-RPG down-regulation upon TOR pathway inhibition was much attenuated at defective mutant promoters unable to bind Abf1. TORC1 inactivation caused the expected reduction of Ifh1 occupancy at RPS22B and RPL3 promoters, but unexpectedly it entailed largely increased Abf1 association with Abf1-RPG promoters. We present evidence that Abf1 recruitment upon nutritional stress, also observed for representative ribosome biogenesis genes, favours RPG transcriptional rescue upon nutrient replenishment, thus pointing to nutrient-regulated Abf1 dynamics at promoters as a novel mechanism in ribosome biogenesis control.

The Power of Natural Variation for Model Organism Biology

Genetic background effects have long been recognized and, in some cases studied, but they are often viewed as a nuisance by molecular biologists. We suggest that genetic variation currently represents a critical frontier for molecular studies. Human genetics has seen a surge of interest in genetic variation and its contributions to disease, but insights into disease mechanisms are difficult since information about gene function is lacking. By contrast, model organism genetics has excelled at revealing molecular mechanisms of cellular processes, but often de-emphasizes genetic variation and its functional consequences. We argue that model organism biology would benefit from incorporating natural variation, both to capture how well laboratory lines exemplify the species they represent and to inform on molecular processes and their variability. Such a synthesis would also greatly expand the relevance of model systems for studies of complex trait variation, including disease.

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.

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.

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.