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

Forward genetic screen in zebrafish identifies new fungal regulators that limit host-protective Candida-innate immune interaction

Candida is one of the most frequent causes of bloodstream infections, and our first line of defense against these invasive infections is the innate immune system. The early immune response is critical in controlling Candida albicans infection, but C. albicans has several strategies to evade host immune attack. Phagocytosis of C. albicans blocks hyphal growth, limiting host damage and virulence, but how C. albicans limits early recruitment and phagocytosis in vertebrate infection is poorly understood. To study innate immune evasion by intravital imaging, we utilized the transparent larval zebrafish infection model to screen 131 C. albicans mutants for altered virulence and phagocyte response. Infections with each of the seven hypovirulent mutants led to altered phagocyte recruitment and/or phagocytosis, falling into four categories. Of particular interest among these is NMD5, a predicted β-importin and newly identified virulence factor. The nmd5∆/∆ mutant fails to limit phagocytosis, and its virulence defects are eliminated when phagocyte activity is compromised, suggesting that its role in virulence is limited to immune evasion. These quantitative intravital imaging experiments are the first to document altered Candida-phagocyte interactions for several additional mutants and clearly distinguish recruitment from phagocytic uptake, suggesting that Candida modulates both events. This initial large-scale screen of individual C. albicans mutants in a vertebrate, coupled with high-resolution imaging of Candida-phagocyte interactions, provides a more nuanced view of how diverse mutations can lead to more effective phagocytosis, a key immune process that blocks germination and drives anti-fungal immunity.

IMPORTANCE

Candida albicans is part of the human microbial community and is a dangerous opportunistic pathogen, able to prevent its elimination by the host immune system. Although Candida avoids immune attack through several strategies, we still understand little about how it regulates when immune phagocytes get recruited to the infection site and when they engulf fungal cells. We tested over 130 selected Candida mutants for their ability to cause lethal infection and found several hypovirulent mutants, which provoked altered innate immune responses, resulting in lower overall inflammation and greater host survival. Of particular interest is NMD5, which acts to limit fungal phagocytosis and is predicted to regulate the activity of stress-associated transcription factors. Our high-content screening was enabled by modeling Candida infection in transparent vertebrate zebrafish larva. Our findings help us understand how Candida survives immune attack during commensal and pathogenic growth, and may eventually inform new strategies for controlling disease.

Transcription Factor SsNdt80b Maintains Optimal Expression of SsSNF1 to Modulate Growth and Pathogenicity in Sclerotinia sclerotiorum

Microorganisms use versatile strategies to facilitate the colonisation of hosts, through remodelling transcription and metabolism to accommodate growth under harsh and hostile environments. Sclerotinia sclerotiorum is a typical necrotrophic pathogen that causes Sclerotinia stem rot in more than 700 species, resulting in serious economic losses. How S. sclerotiorum integrates mechanisms for nutrient acquisition and utilisation to maintain optimal growth and pathogenicity is still indistinct. Here, we demonstrate that Ndt80 family transcription factors (SsNdt80a,b,c) are involved in carbon source utilisation and have different roles in the growth, sclerotia formation, infection cushion development, and the virulence of S. sclerotiorum. SsNdt80b could bind the promoter of SsSNF1 and modulate the transcriptional activity of SsSNF1. Silencing SsSNF1 resulted in defects in hyphal growth and infection cushion formation, reduced cell wall-degrading enzymes, and reduced pathogenicity of S. sclerotiorum. A model is proposed in which SsNdt80b responds to carbon sources and modulates SsSnf1 to regulate the development and pathogenicity of S. sclerotiorum.

Forward genetic screen in zebrafish identifies new fungal regulators that limit host-protective Candida-innate immune interaction

Candida is one of the most frequent causes of bloodstream infections, and our first line of defense against these invasive infections is the innate immune system. The early immune response is critical in controlling C. albicans infection, but C. albicans has several strategies to evade host immune attack. Phagocytosis of C. albicans blocks hyphal growth, limiting host damage and virulence, but how C. albicans limits early recruitment and phagocytosis in vertebrate infection is poorly understood. To study innate immune evasion by intravital imaging, we utilized the transparent larval zebrafish infection model to screen 131 C. albicans mutants for altered virulence and phagocyte response. Infections with each of seven hypovirulent mutants led to altered phagocyte recruitment and/or phagocytosis, falling into four categories. Of particular interest among these is NMD5, a predicted β-importin and newly-identified virulence factor. The nmd5∆/∆ mutant fails to limit phagocytosis and its virulence defects are eliminated when phagocyte activity is compromised, suggesting that its role in virulence is limited to immune evasion. These quantitative intravital imaging experiments are the first to document altered Candida-phagocyte interactions for several additional mutants, and clearly distinguish recruitment from phagocytic uptake, suggesting that Candida modulates both events. This initial large-scale screen of individual C. albicans mutants in a vertebrate, coupled with high-resolution imaging of Candida-phagocyte interactions, provides a more nuanced view of how diverse mutations can lead to more effective phagocytosis, a key immune process which blocks germination and drives anti-fungal immunity.

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.

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

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

Molecular Determinants Involved in Candida albicans Biofilm Formation and Regulation

Candida albicans is known for its pathogenicity, although it lives within the human body as a commensal member. The commensal nature of C. albicans is well controlled and regulated by the host's immune system as they live in the harmonized microenvironment. However, the development of certain unusual microhabitat conditions (change in pH, co-inhabiting microorganisms' population ratio, debilitated host-immune system) pokes this commensal fungus to transform into a pathogen in such a way that it starts to propagate very rapidly and tries to breach the epithelial barrier to enter the host's systemic circulations. In addition, Candida is infamous as a major nosocomial (hospital-acquired infection) agent because it enters the human body through venous catheters or medical prostheses. The hysterical mode of C. albicans growth builds its microcolony or biofilm, which is pathogenic for the host. Biofilms propose additional resistance mechanisms from host immunity or extracellular chemicals to aid their survival. Differential gene expressions and regulations within the biofilms cause altered morphology and metabolism. The genes associated with adhesiveness, hyphal/pseudo-hyphal growth, persister cell transformation, and biofilm formation by C. albicans are controlled by myriads of cell-signaling regulators. These genes' transcription is controlled by different molecular determinants like transcription factors and regulators. Therefore, this review has focused discussion on host-immune-sensing molecular determinants of Candida during biofilm formation, regulatory descriptors (secondary messengers, regulatory RNAs, transcription factors) of Candida involved in biofilm formation that could enable small-molecule drug discovery against these molecular determinants, and lead to disrupt the well-structured Candida biofilms effectively.

RFX transcription factor in the human-associated yeast Candida albicans regulates adhesion to oral epithelium

Adhesion to mucosal surfaces is a critical step in many bacterial and fungal infections. Here, using a mouse model of oral infection by the human fungal pathobiont Candida albicans, we report the identification of a novel regulator of C. albicans adhesion to the oral mucosa. The regulator is a member of the regulatory factor X (RFX) family of transcription factors, which control cellular processes ranging from genome integrity in model yeasts to tissue differentiation in vertebrates. Mice infected with the C. albicans rfx1 deletion mutant displayed increased fungal burden in tongues compared to animals infected with the reference strain. High-resolution imaging revealed RFX1 transcripts being expressed by C. albicans cells during infection. Concomitant with the increase in fungal burden, the rfx1 mutant elicited an enhanced innate immune response. Transcriptome analyses uncovered HWP1, a gene encoding an adhesion protein that mediates covalent attachment to buccal cells, as a major RFX1-regulated locus. Consistent with this result, we establish that C. albicans adhesion to oral cells is modulated by RFX1 in an HWP1-dependent manner. Our findings expand the repertoire of biological processes controlled by the RFX family and illustrate a mechanism whereby C. albicans can adjust adhesion to the oral epithelium.

Inferring gene regulatory networks using transcriptional profiles as dynamical attractors

Genetic regulatory networks (GRNs) regulate the flow of genetic information from the genome to expressed messenger RNAs (mRNAs) and thus are critical to controlling the phenotypic characteristics of cells. Numerous methods exist for profiling mRNA transcript levels and identifying protein-DNA binding interactions at the genome-wide scale. These enable researchers to determine the structure and output of transcriptional regulatory networks, but uncovering the complete structure and regulatory logic of GRNs remains a challenge. The field of GRN inference aims to meet this challenge using computational modeling to derive the structure and logic of GRNs from experimental data and to encode this knowledge in Boolean networks, Bayesian networks, ordinary differential equation (ODE) models, or other modeling frameworks. However, most existing models do not incorporate dynamic transcriptional data since it has historically been less widely available in comparison to "static" transcriptional data. We report the development of an evolutionary algorithm-based ODE modeling approach (named EA) that integrates kinetic transcription data and the theory of attractor matching to infer GRN architecture and regulatory logic. Our method outperformed six leading GRN inference methods, none of which incorporate kinetic transcriptional data, in predicting regulatory connections among TFs when applied to a small-scale engineered synthetic GRN in Saccharomyces cerevisiae. Moreover, we demonstrate the potential of our method to predict unknown transcriptional profiles that would be produced upon genetic perturbation of the GRN governing a two-state cellular phenotypic switch in Candida albicans. We established an iterative refinement strategy to facilitate candidate selection for experimentation; the experimental results in turn provide validation or improvement for the model. In this way, our GRN inference approach can expedite the development of a sophisticated mathematical model that can accurately describe the structure and dynamics of the in vivo GRN.

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.

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

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

Environment-dependent epistasis increases phenotypic diversity in gene regulatory networks

Mutations to gene regulatory networks can be maladaptive or a source of evolutionary novelty. Epistasis confounds our understanding of how mutations affect the expression patterns of gene regulatory networks, a challenge exacerbated by the dependence of epistasis on the environment. We used the toolkit of synthetic biology to systematically assay the effects of pairwise and triplet combinations of mutant genotypes on the expression pattern of a gene regulatory network expressed in Escherichia coli that interprets an inducer gradient across a spatial domain. We uncovered a preponderance of epistasis that can switch in magnitude and sign across the inducer gradient to produce a greater diversity of expression pattern phenotypes than would be possible in the absence of such environment-dependent epistasis. We discuss our findings in the context of the evolution of hybrid incompatibilities and evolutionary novelties.

Robustness and innovation in synthetic genotype networks

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

Primate protein-ligand interfaces exhibit significant conservation and unveil human-specific evolutionary drivers

Despite the vast phenotypic differences observed across primates, their protein products are largely similar to each other at the sequence level. We hypothesized that, since proteins accomplish all their functions via interactions with other molecules, alterations in the sites that participate in these interactions may be of critical importance. To uncover the extent to which these sites evolve across primates, we built a structurally-derived dataset of ~4,200 one-to-one orthologous sequence groups across 18 primate species, consisting of ~68,000 ligand-binding sites that interact with DNA, RNA, small molecules, ions, or peptides. Using this dataset, we identify functionally important patterns of conservation and variation within the amino acid residues that facilitate protein-ligand interactions across the primate phylogeny. We uncover that interaction sites are significantly more conserved than other sites, and that sites binding DNA and RNA further exhibit the lowest levels of variation. We also show that the subset of ligand-binding sites that do vary are enriched in components of gene regulatory pathways and uncover several instances of human-specific ligand-binding site changes within transcription factors. Altogether, our results suggest that ligand-binding sites have experienced selective pressure in primates and propose that variation in these sites may have an outsized effect on phenotypic variation in primates through pleiotropic effects on gene regulation.

The GRN concept as a guide for evolutionary developmental biology

Organismal phenotypes result largely from inherited developmental programs, usually executed during embryonic and juvenile life stages. These programs are not blank slates onto which natural selection can draw arbitrary forms. Rather, the mechanisms of development play an integral role in shaping phenotypic diversity and help determine the evolutionary trajectories of species. Modern evolutionary biology must, therefore, account for these mechanisms in both theory and in practice. The gene regulatory network (GRN) concept represents a potent tool for achieving this goal whose utility has grown in tandem with advances in "omic" technologies and experimental techniques. However, while the GRN concept is widely utilized, it is often less clear what practical implications it has for conducting research in evolutionary developmental biology. In this Perspective, we attempt to provide clarity by discussing how experiments and projects can be designed in light of the GRN concept. We first map familiar biological notions onto the more abstract components of GRN models. We then review how diverse functional genomic approaches can be directed toward the goal of constructing such models and discuss current methods for functionally testing evolutionary hypotheses that arise from them. Finally, we show how the major steps of GRN model construction and experimental validation suggest generalizable workflows that can serve as a scaffold for project design. Taken together, the practical implications that we draw from the GRN concept provide a set of guideposts for studies aiming at unraveling the molecular basis of phenotypic diversity.

Recognition of galactose by a scaffold protein recruits a transcriptional activator for the GAL regulon induction in Candida albicans

The GAL pathway of yeasts has long served as a model system for understanding of how regulatory mode of eukaryotic metabolic pathways evolves. While Gal4 mode has been well-characterized in Saccharomycetaceae clade, little is known about the regulation of the GAL pathway in other yeasts. Here, we find that Rep1, a Ndt80-like family transcription factor, serves as a galactose sensor in the commensal-pathogenic fungus Candida albicans. It is presented at the GAL gene promoters independent of the presence of galactose. Rep1 recognizes galactose via a direct physical interaction. The net result of this interaction is the recruitment of a transcriptional activator Cga1 (Candida galactose gene activator, orf19.4959) and transcription of the GAL genes proceeds. Rep1 and Cga1 are conserved across the CTG species. Rep1 itself does not possess transcriptional activity. Instead, it provides a scaffold to recruit different factors for transcriptional regulation. Rep1-Cga1 mode of regulation represents a new example of network rewiring in fungi, which provides insight into how C. albicans evolves transcriptional programs to colonize diverse host niches.

Evolution of yeast hybrids by aborted meiosis

Sterile hybrids are broadly considered evolutionary dead-ends because of their faulty sexual reproduction. While sterility in obligate sexual organisms is a clear constraint in perpetuating the species, some facultative sexual microbes such as yeasts can propagate asexually and maintain genome plasticity. Moreover, incomplete meiotic pathways in yeasts represent alternative routes to the standard meiosis that generates genetic combinations in the population and fuel adaptation. Here, we review how aborting meiosis promotes genome-wide allele shuffling in sterile Saccharomyces hybrids and describe approaches to identify evolved clones in a cell population. We further discuss possible implications of this process in generating phenotypic novelty and report cases of abortive meiosis across yeast species.

Transcription Factor Mavib-1 Negatively Regulates Conidiation by Affecting Utilization of Carbon and Nitrogen Source in Metarhizium acridum

Conidium is the main infection unit and reproductive unit of pathogenic fungi. Exploring the mechanism of conidiation and its regulation contributes to understanding the pathogenicity of pathogenic fungi. Vib-1, a transcription factor, was reported to participate in the conidiation process. However, the regulation mechanism of Vib-1 in conidiation is still unclear. In this study, we analyzed the function of Vib-1 and its regulation mechanism in conidiation through knocking out and overexpression of Vib-1 in entomopathogenic fungus Metarhizium acridum. Results showed that the colonial growth of Mavib-1 disruption mutant (ΔMavib-1) was significantly decreased, and conidiation was earlier compared to wild type (WT), while overexpression of Mavib-1 led to a delayed conidiation especially when carbon or nitrogen sources were insufficient. Overexpression of Mavib-1 resulted in a conidiation pattern shift from microcycle conidiation to normal conidiation on nutrient-limited medium. These results indicated that Mavib-1 acted as a positive regulator in vegetative growth and a negative regulator in conidiation by affecting utilization of carbon and nitrogen sources in M. acridum. Transcription profile analysis demonstrated that many genes related to carbon and nitrogen source metabolisms were differentially expressed in ΔMavib-1 and OE strains compared to WT. Moreover, Mavib-1 affects the conidial germination, tolerance to UV-B and heat stresses, cell wall integrity, conidial surface morphology and conidial hydrophobicity in M. acridum. These findings unravel the regulatory mechanism of Mavib-1 in fungal growth and conidiation, and enrich the knowledge to conidiation pattern shift of filamentous fungi.

Intrinsic transcriptional heterogeneity in neuroblastoma guides mechanistic and therapeutic insights

Cell state is controlled by master transcription factors (mTFs) that determine the cellular gene expression program. Cancer cells acquire dysregulated gene expression programs by mutational and non-mutational processes. Intratumoral heterogeneity can result from cells displaying distinct mTF-regulated cell states, which co-exist within the tumor. One archetypal tumor associated with transcriptionally regulated heterogeneity is high-risk neuroblastoma (NB). Patients with NB have poor overall survival despite intensive therapies, and relapsed patients are commonly refractory to treatment. The cellular populations that comprise NB are marked by different cohorts of mTFs and differential sensitivity to conventional therapies. Recent studies have highlighted mechanisms by which NB cells dynamically shift the cell state with treatment, revealing new opportunities to control the cellular response to treatment by manipulating cell-state-defining transcriptional programs. Here, we review recent advances in understanding transcriptionally defined cancer heterogeneity. We offer challenges to the field to encourage translation of basic science into clinical benefit.

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.

A Fungal Transcription Regulator of Vacuolar Function Modulates Candida albicans Interactions with Host Epithelial Cells

Microorganisms typically maintain cellular homeostasis despite facing large fluctuations in their surroundings. Microbes that reside on human mucosal surfaces may experience significant variations in nutrient and ion availability as well as pH. Whether the mechanisms employed by these microbial cells to sustain homeostasis directly impact on the interplay with the host’s mucosae remains unclear. Here, we report that the previously uncharacterized transcription regulator ZCF8 in the human-associated yeast Candida albicans maintains vacuole homeostasis when the fungus faces fluctuations in nitrogen. Genome-wide identification of genes directly regulated by Zcf8p followed by fluorescence microscopy to define their subcellular localization uncovered the fungal vacuole as a top target of Zcf8p regulation. Deletion and overexpression of ZCF8 resulted in alterations in vacuolar morphology and luminal pH and rendered the fungus resistant or susceptible to nigericin and brefeldin A, two drugs that impair vacuole and associated functions. Furthermore, we establish that the regulator modulates C. albicans attachment to epithelial cells in a manner that depends on the status of the fungal vacuole. Our findings, therefore, suggest that fungal vacuole physiology regulation is intrinsically linked to, and shapes to a significant extent, the physical interactions that Candida cells establish with mammalian mucosal surfaces.

Recurrent Loss of abaA, a Master Regulator of Asexual Development in Filamentous Fungi, Correlates with Changes in Genomic and Morphological Traits

Gene regulatory networks (GRNs) drive developmental and cellular differentiation, and variation in their architectures gives rise to morphological diversity. Pioneering studies in Aspergillus fungi, coupled with subsequent work in other filamentous fungi, have shown that the GRN governed by the BrlA, AbaA, and WetA proteins controls the development of the asexual fruiting body or conidiophore. A specific aspect of conidiophore development is the production of phialides, conidiophore structures that are under the developmental control of AbaA and function to repetitively generate spores. Fungal genome sequencing has revealed that some filamentous fungi lack abaA, and also produce asexual structures that lack phialides, raising the hypothesis that abaA loss is functionally linked to diversity in asexual fruiting body morphology. To examine this hypothesis, we carried out an extensive search for the abaA gene across 241 genomes of species from the fungal subphylum Pezizomycotina. We found that abaA was independently lost in four lineages of Eurotiomycetes, including from all sequenced species within the order Onygenales, and that all four lineages that have lost abaA also lack the ability to form phialides. Genetic restoration of abaA from Aspergillus nidulans into Histoplasma capsulatum, a pathogenic species from the order Onygenales that lacks an endogenous copy of abaA, did not alter Histoplasma conidiation morphology but resulted in a marked increase in spore viability. We also discovered that species lacking abaA contain fewer AbaA binding motifs in the regulatory regions of orthologs of some AbaA target genes, suggesting that the asexual fruiting body GRN of organisms that have lost abaA has likely been rewired. Our results provide an illustration of how repeated losses of a key regulatory transcription factor have contributed to the diversity of an iconic fungal morphological trait.

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.

Evolution of the complex transcription network controlling biofilm formation in Candida species

We examine how a complex transcription network composed of seven 'master' regulators and hundreds of target genes evolved over a span of approximately 70 million years. The network controls biofilm formation in several Candida species, a group of fungi that are present in humans both as constituents of the microbiota and as opportunistic pathogens. Using a variety of approaches, we observed two major types of changes that have occurred in the biofilm network since the four extant species we examined last shared a common ancestor. Master regulator 'substitutions' occurred over relatively long evolutionary times, resulting in different species having overlapping but different sets of master regulators of biofilm formation. Second, massive changes in the connections between the master regulators and their target genes occurred over much shorter timescales. We believe this analysis is the first detailed, empirical description of how a complex transcription network has evolved.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

High-resolution analysis of cell-state transitions in yeast suggests widespread transcriptional tuning by alternative starts

BACKGROUND

The start and end sites of messenger RNAs (TSSs and TESs) are highly regulated, often in a cell-type-specific manner. Yet the contribution of transcript diversity in regulating gene expression remains largely elusive. We perform an integrative analysis of multiple highly synchronized cell-fate transitions and quantitative genomic techniques in Saccharomyces cerevisiae to identify regulatory functions associated with transcribing alternative isoforms.

RESULTS

Cell-fate transitions feature widespread elevated expression of alternative TSS and, to a lesser degree, TES usage. These dynamically regulated alternative TSSs are located mostly upstream of canonical TSSs, but also within gene bodies possibly encoding for protein isoforms. Increased upstream alternative TSS usage is linked to various effects on canonical TSS levels, which range from co-activation to repression. We identified two key features linked to these outcomes: an interplay between alternative and canonical promoter strengths, and distance between alternative and canonical TSSs. These two regulatory properties give a plausible explanation of how locally transcribed alternative TSSs control gene transcription. Additionally, we find that specific chromatin modifiers Set2, Set3, and FACT play an important role in mediating gene repression via alternative TSSs, further supporting that the act of upstream transcription drives the local changes in gene transcription.

CONCLUSIONS

The integrative analysis of multiple cell-fate transitions suggests the presence of a regulatory control system of alternative TSSs that is important for dynamic tuning of gene expression. Our work provides a framework for understanding how TSS heterogeneity governs eukaryotic gene expression, particularly during cell-fate changes.

Meiotic Cells Counteract Programmed Retrotransposon Activation via RNA-Binding Translational Repressor Assemblies

Retrotransposon proliferation poses a threat to germline integrity. While retrotransposons must be activated in developing germ cells in order to survive and propagate, how they are selectively activated in the context of meiosis is unclear. We demonstrate that the transcriptional activation of Ty3/Gypsy retrotransposons and host defense are controlled by master meiotic regulators. We show that budding yeast Ty3/Gypsy co-opts binding sites of the essential meiotic transcription factor Ndt80 upstream of the integration site, thereby tightly linking its transcriptional activation to meiotic progression. We also elucidate how yeast cells thwart Ty3/Gypsy proliferation by blocking translation of the retrotransposon mRNA using amyloid-like assemblies of the RNA-binding protein Rim4. In mammals, several inactive Ty3/Gypsy elements are undergoing domestication. We show that mammals utilize equivalent master meiotic regulators (Stra8, Mybl1, Dazl) to regulate Ty3/Gypsy-derived genes in developing gametes. Our findings inform how genes that are evolving from retrotransposons can build upon existing regulatory networks during domestication.

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

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

Which came first, the transcriptional regulator or its target genes? An evolutionary perspective into the construction of eukaryotic regulons

Eukaryotic regulons are regulatory units formed by a set of genes under the control of the same transcription factor (TF). Despite the functional plasticity, TFs are highly conserved and recognize the same DNA sequences in different organisms. One of the main factors that confer regulatory specificity is the distribution of the binding sites of the TFs along the genome, allowing the configuration of different transcriptional regulatory networks (TRNs) from the same regulator. A similar scenario occurs between tissues of the same organism, where a TRN can be rewired by epigenetic factors, modulating the accessibility of the TF to its binding sites. In this article we discuss concepts that can help to formulate testable hypotheses about the construction of regulons, exploring the presence and absence of the elements that form a TRN throughout the evolution of an ancestral lineage. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

A 'parameiosis' drives depolyploidization and homologous recombination in Candida albicans

Meiosis is a conserved tenet of sexual reproduction in eukaryotes, yet this program is seemingly absent from many extant species. In the human fungal pathogen Candida albicans, mating of diploid cells generates tetraploid products that return to the diploid state via a non-meiotic process of depolyploidization known as concerted chromosome loss (CCL). Here, we report that recombination rates are more than three orders of magnitude higher during CCL than during normal mitotic growth. Furthermore, two conserved ‘meiosis-specific’ factors play central roles in CCL as SPO11 mediates DNA double-strand break formation while both SPO11 and REC8 regulate chromosome stability and promote inter-homolog recombination. Unexpectedly, SPO11 also promotes DNA repair and recombination during normal mitotic divisions. These results indicate that C. albicans CCL represents a ‘parameiosis’ that blurs the conventional boundaries between mitosis and meiosis. They also reveal parallels with depolyploidization in mammalian cells and provide potential insights into the evolution of meiosis.

Mating of Candida albicans produces tetraploid products that return to the diploid state via a non-meiotic process known as concerted chromosome loss (CCL). Here, Anderson et al. show high recombination rates during CCL and identify factors that are essential for chromosome stability and recombination during CCL.

The causes of evolvability and their evolution

Evolvability is the ability of a biological system to produce phenotypic variation that is both heritable and adaptive. It has long been the subject of anecdotal observations and theoretical work. In recent years, however, the molecular causes of evolvability have been an increasing focus of experimental work. Here, we review recent experimental progress in areas as different as the evolution of drug resistance in cancer cells and the rewiring of transcriptional regulation circuits in vertebrates. This research reveals the importance of three major themes: multiple genetic and non-genetic mechanisms to generate phenotypic diversity, robustness in genetic systems, and adaptive landscape topography. We also discuss the mounting evidence that evolvability can evolve and the question of whether it evolves adaptively.

Environmental and Evolutionary Drivers of the Modular Gene Regulatory Network Underlying Phenotypic Plasticity for Stress Resistance in the Nematode Caenorhabditis remanei

Organisms can cope with stressful environments via a combination of phenotypic plasticity at the individual level and adaptation at the population level. Changes in gene expression can play an important role in both. Significant advances in our understanding of gene regulatory plasticity and evolution have come from comparative studies in the field and laboratory. Experimental evolution provides another powerful path by which to learn about how differential regulation of genes and pathways contributes to both acclimation and adaptation. Here we present results from one such study using the nematode Caenorhabditis remanei We selected one set of lines to withstand heat stress and another oxidative stress. We then compared transcriptional responses to acute heat stress of both and an unselected control to the ancestral population using a weighted gene coexpression network analysis, finding that the transcriptional response is primarily dominated by a plastic response that is shared in the ancestor and all evolved populations. In addition, we identified several modules that respond to artificial selection by (1) changing the baseline level of expression, (2) altering the magnitude of the plastic response, or (3) a combination of the two. Our findings therefore reveal that while patterns of transcriptional response can be perturbed with short bouts of intense selection, the overall ancestral structure of transcriptional plasticity is largely maintained over time.

Candida albicans: An Emerging Yeast Model to Study Eukaryotic Genome Plasticity

Saccharomyces cerevisiae and Schizosaccharomyces pombe have served as uncontested unicellular model organisms, as major discoveries made in the field of genome biology using yeast genetics have proved to be relevant from yeast to humans. The yeast Candida albicans has attracted much attention because of its ability to switch between a harmless commensal and a dreaded human pathogen. C. albicans bears unique features regarding its life cycle, genome structure, and dynamics, and their links to cell biology and adaptation to environmental challenges. Examples include a unique reproduction cycle with haploid, diploid, and tetraploid forms; a distinctive organisation of chromosome hallmarks; a highly dynamic genome, with extensive karyotypic variations, including aneuploidies, isochromosome formation, and loss-of-heterozygosity; and distinctive links between the response to DNA alterations and cell morphology. These features have made C. albicans emerge as a new and attractive unicellular model to study genome biology and dynamics in eukaryotes.

Developmental regulation of an organelle tether coordinates mitochondrial remodeling in meiosis

Cellular differentiation involves remodeling cellular architecture to transform one cell type to another. By investigating mitochondrial dynamics during meiotic differentiation in budding yeast, we sought to understand how organelle morphogenesis is developmentally controlled in a system where regulators of differentiation and organelle architecture are known, but the interface between them remains unexplored. We analyzed the regulation of mitochondrial detachment from the cell cortex, a known meiotic alteration to mitochondrial morphology. We found that mitochondrial detachment is enabled by the programmed destruction of the mitochondria-endoplasmic reticulum-cortex anchor (MECA), an organelle tether that bridges mitochondria and the plasma membrane. MECA regulation is governed by a meiotic transcription factor, Ndt80, which promotes the activation of a conserved kinase, Ime2. We further present evidence for Ime2-dependent phosphorylation and degradation of MECA in a temporally controlled manner. Our study defines a key mechanism that coordinates mitochondrial morphogenesis with the landmark events of meiosis and demonstrates that cells can developmentally regulate tethering to induce organelle remodeling.

Genetic Analysis of NDT80 Family Transcription Factors in Candida albicans Using New CRISPR-Cas9 Approaches

Ndt80 family transcription factors are highly conserved in fungi, where they regulate diverse processes. The human fungal pathogen Candida albicans contains three genes (NDT80, REP1, and RON1) that encode proteins with similarity to Saccharomyces cerevisiae Ndt80, although the homology is restricted to the DNA binding domain. To better understand their role in virulence functions, we used clustered regularly interspaced short palindromic repeat/CRISPR-associated gene 9 (CRISPR/Cas9) to delete the three NDT80-family genes. An ndt80Δ mutant showed strong defects in forming hyphae in response to serum or N-acetylglucosamine (GlcNAc), which was linked to the ability of Ndt80 to regulate the expression of RAS1, an upstream regulator of hyphal signaling. Conversely, the ndt80Δ mutant formed hyphal cells on glycerol medium, indicating that Ndt80 is not required for hyphal growth under all conditions. In contrast to our previously published data, a ron1Δ single mutant could grow and form hyphae in response to GlcNAc. However, deleting RON1 partially restored the ability of an ndt80Δ mutant to form hyphae in response to GlcNAc, indicating a link to GlcNAc signaling. REP1 was required for growth on GlcNAc, as expected, but not for GlcNAc or serum to induce hyphae. The ndt80Δ mutant was defective in growing under stressful conditions, such as elevated temperature, but not the ron1Δ mutant or rep1Δ mutant. Quantitative assays did not reveal any significant differences in the fluconazole susceptibility of the NDT80-family mutants. Interestingly, double and triple mutant analysis did not identify significant genetic interactions for these NDT80 family genes, indicating that they mainly function independently, in spite of their conserved DNA binding domain.IMPORTANCE Transcription factors play key roles in regulating virulence of the human fungal pathogen C. albicans In addition to regulating the expression of virulence factors, they also control the ability of C. albicans to switch to filamentous hyphal growth, which facilitates biofilm formation on medical devices and invasion into tissues. We therefore used new CRISPR/Cas9 methods to examine the effects of deleting three C. albicans genes (NDT80, REP1, and RON1) that encode transcription factors with similar DNA binding domains. Interestingly, double and triple mutant strains mostly showed the combined properties of the single mutants; there was only very limited evidence of synergistic interactions in regulating morphogenesis, stress resistance, and ability to metabolize different sugars. These results demonstrate that NDT80, REP1, and RON1 have distinct functions in regulating C. albicans virulence functions.

Global analysis of genetic circuitry and adaptive mechanisms enabling resistance to the azole antifungal drugs

Invasive fungal infections caused by the pathogen Candida albicans have transitioned from a rare curiosity to a major cause of human mortality. This is in part due to the emergence of resistance to the limited number of antifungals available to treat fungal infections. Azoles function by targeting the biosynthesis of ergosterol, a key component of the fungal cell membrane. Loss-of-function mutations in the ergosterol biosynthetic gene ERG3 mitigate azole toxicity and enable resistance that depends upon fungal stress responses. Here, we performed a genome-wide synthetic genetic array screen in Saccharomyces cerevisiae to map ERG3 genetic interactors and uncover novel circuitry important for azole resistance. We identified nine genes that enabled erg3-mediated azole resistance in the model yeast and found that only two of these genes had a conserved impact on resistance in C. albicans. Further, we screened a C. albicans homozygous deletion mutant library and identified 13 genes for which deletion enhances azole susceptibility. Two of the genes, RGD1 and PEP8, were also important for azole resistance acquired by diverse mechanisms. We discovered that loss of function of retrograde transport protein Pep8 overwhelms the functional capacity of the stress response regulator calcineurin, thereby abrogating azole resistance. To identify the mechanism through which the GTPase activator protein Rgd1 enables azole resistance, we selected for mutations that restore resistance in strains lacking Rgd1. Whole genome sequencing uncovered parallel adaptive mechanisms involving amplification of both chromosome 7 and a large segment of chromosome 3. Overexpression of a transporter gene on the right portion of chromosome 3, NPR2, was sufficient to enable azole resistance in the absence of Rgd1. Thus, we establish a novel mechanism of adaptation to drug-induced stress, define genetic circuitry underpinning azole resistance, and illustrate divergence in resistance circuitry over evolutionary time.

Fungal infections caused by the pathogen Candida albicans pose a serious threat to human health. Treating these infections relies heavily on the azole antifungals, however, resistance to these drugs develops readily demanding novel therapeutic strategies. We performed large-scale systematic screens in both C. albicans and the model yeast Saccharomyces cerevisiae to identify genes that enable azole resistance. Our genome-wide screen in S. cerevisiae identified nine determinants of azole resistance, only two of which were important for resistance in C. albicans. Our screen of C. albicans mutants identified 13 genes for which deletion enhances susceptibility to azoles, including RGD1 and PEP8. We found that loss of Pep8 overwhelms the functional capacity of a key stress response regulator, calcineurin. In contrast, amplification of chromosome 7 and the right portion of chromosome 3 can restore resistance in strains lacking Rgd1, suggesting that Rgd1 may enable azole resistance by inducing genes in these amplified regions. Specifically, overexpression of a gene involved in transport on chromosome 3, NPR2, was sufficient to restore azole resistance in the absence of Rgd1. Thus, we establish novel circuitry important for antifungal drug resistance, and uncover adaptive mechanisms involving genomic plasticity that occur in response to drug induced stress.

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.

The rewiring of transcription circuits in evolution

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