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

Significance Of the many microbial species on earth, only a small number are able to thrive in humans and cause disease. Comparison of closely related pathogenic and nonpathogenic species can therefore be useful in identifying key features that contribute to virulence. We created interspecies hybrids between Candida albicans, a prevalent fungal pathogen of humans, and Candida dubliniensis, a close, but much less pathogenic, relative. By comparing genome-wide expression differences between the two genomes in the same cell, we surmised that since the two species diverged from a common ancestor, natural selection has acted upon the expression level of an ancient metabolic pathway, illustrating that pathogenicity traits can arise over evolutionary timescales through small expression changes in deeply conserved proteins. Candida albicans is the most common cause of systemic fungal infections in humans and is considerably more virulent than its closest known relative, Candida dubliniensis. To investigate this difference, we constructed interspecies hybrids and quantified mRNA levels produced from each genome in the hybrid. This approach systematically identified expression differences in orthologous genes arising from cis-regulatory sequence changes that accumulated since the two species last shared a common ancestor, some 10 million y ago. We documented many orthologous gene-expression differences between the two species, and we pursued one striking observation: All 15 genes coding for the enzymes of glycolysis showed higher expression from the C. albicans genome than the C. dubliniensis genome in the interspecies hybrid. This pattern requires evolutionary changes to have occurred at each gene; the fact that they all act in the same direction strongly indicates lineage-specific natural selection as the underlying cause. To test whether these expression differences contribute to virulence, we created a C. dubliniensis strain in which all 15 glycolysis genes were produced at modestly elevated levels and found that this strain had significantly increased virulence in the standard mouse model of systemic infection. These results indicate that small expression differences across a deeply conserved set of metabolism enzymes can play a significant role in the evolution of virulence in fungal pathogens.

[1]  T. Gabaldón,et al.  Genomic evidence for a hybrid origin of the yeast opportunistic pathogen Candida albicans , 2019, BMC Biology.

[2]  M. Lemaire,et al.  Hexokinase and Glucokinases Are Essential for Fitness and Virulence in the Pathogenic Yeast Candida albicans , 2018, bioRxiv.

[3]  G. Moran,et al.  Expansion of the TLO gene family enhances the virulence of Candida species , 2018, PloS one.

[4]  S. Nuzhdin,et al.  The Evolution of Gene Expression in cis and trans. , 2018, Trends in genetics : TIG.

[5]  M. Hickey,et al.  Glucose Homeostasis Is Important for Immune Cell Viability during Candida Challenge and Host Survival of Systemic Fungal Infection. , 2018, Cell metabolism.

[6]  P. Wittkopp,et al.  Evolutionary Dynamics of Regulatory Changes Underlying Gene Expression Divergence among Saccharomyces Species , 2017, Genome biology and evolution.

[7]  Justin C. Fay,et al.  Cis-Regulatory Divergence in Gene Expression between Two Thermally Divergent Yeast Species , 2017, bioRxiv.

[8]  M. Lorenz,et al.  A feast for Candida: Metabolic plasticity confers an edge for virulence , 2017, PLoS pathogens.

[9]  Marek S. Skrzypek,et al.  The Candida Genome Database (CGD): incorporation of Assembly 22, systematic identifiers and visualization of high throughput sequencing data , 2016, Nucleic Acids Res..

[10]  D. MacCallum,et al.  Amplification of TLO Mediator Subunit Genes Facilitate Filamentous Growth in Candida Spp. , 2016, PLoS genetics.

[11]  Andreas Wagner,et al.  Growthcurver: an R package for obtaining interpretable metrics from microbial growth curves , 2016, BMC Bioinformatics.

[12]  Rachel B. Brem,et al.  Polygenic evolution of a sugar specialization trade-off in yeast , 2016, Nature.

[13]  Alexander D. Johnson,et al.  Intersecting transcription networks constrain gene regulatory evolution , 2015, Nature.

[14]  Hongye Li,et al.  Evolution of Gene Regulation during Transcription and Translation , 2015, Genome biology and evolution.

[15]  Emily K. Tsang,et al.  Genetic conflict reflected in tissue-specific maps of genomic imprinting in human and mouse , 2015, Nature Genetics.

[16]  M. Netea,et al.  Metabolism impacts upon Candida immunogenicity and pathogenicity at multiple levels , 2014, Trends in microbiology.

[17]  David W Williams,et al.  Comparative adherence of Candida albicans and Candida dubliniensis to human buccal epithelial cells and extracellular matrix proteins. , 2014, Medical mycology.

[18]  Hunter B. Fraser,et al.  Evolution at two levels of gene expression in yeast , 2013, Genome research.

[19]  Hunter B. Fraser,et al.  The Molecular Mechanism of a Cis-Regulatory Adaptation in Yeast , 2013, PLoS genetics.

[20]  B. Böttcher,et al.  Global Transcriptome Sequencing Identifies Chlamydospore Specific Markers in Candida albicans and Candida dubliniensis , 2013, PloS one.

[21]  S. Lorenz,et al.  Species and condition specific adaptation of the transcriptional landscapes in Candida albicans and Candida dubliniensis , 2013, BMC Genomics.

[22]  N. Gow,et al.  The Evolutionary Rewiring of Ubiquitination Targets Has Reprogrammed the Regulation of Carbon Assimilation in the Pathogenic Yeast Candida albicans , 2012, mBio.

[23]  M. Netea,et al.  Growth of Candida albicans Cells on the Physiologically Relevant Carbon Source Lactate Affects Their Recognition and Phagocytosis by Immune Cells , 2012, Infection and Immunity.

[24]  A. Peña,et al.  Characterization of glycolytic metabolism and ion transport of Candida albicans , 2012, Yeast.

[25]  N. Gow,et al.  Host carbon sources modulate cell wall architecture, drug resistance and virulence in a fungal pathogen , 2012, Cellular microbiology.

[26]  Sasha F. Levy,et al.  Polygenic cis-regulatory adaptation in the evolution of yeast pathogenicity , 2012, Genome research.

[27]  Rachel B. Brem,et al.  Evolution of a Membrane Protein Regulon in Saccharomyces , 2012, Molecular biology and evolution.

[28]  Giovanni Parmigiani,et al.  Integrating diverse genomic data using gene sets , 2011, Genome Biology.

[29]  G. Moran,et al.  Candida albicans versus Candida dubliniensis: Why Is C. albicans More Pathogenic? , 2011, International journal of microbiology.

[30]  Hunter B. Fraser,et al.  Genome-wide approaches to the study of adaptive gene expression evolution: systematic studies of evolutionary adaptations involving gene expression will allow many fundamental questions in evolutionary biology to be addressed. , 2011, BioEssays : news and reviews in molecular, cellular and developmental biology.

[31]  C. d’Enfert,et al.  Contribution of the glycolytic flux and hypoxia adaptation to efficient biofilm formation by Candida albicans , 2011, Molecular microbiology.

[32]  J. Heitman,et al.  Calcineurin Controls Drug Tolerance, Hyphal Growth, and Virulence in Candida dubliniensis , 2011, Eukaryotic Cell.

[33]  C Joel McManus,et al.  Regulatory divergence in Drosophila revealed by mRNA-seq. , 2010, Genome research.

[34]  Victoria Chen,et al.  Systematic screens of a Candida albicans homozygous deletion library decouple morphogenetic switching and pathogenicity , 2010, Nature Genetics.

[35]  M. Robinson,et al.  A scaling normalization method for differential expression analysis of RNA-seq data , 2010, Genome Biology.

[36]  Sandrine Dudoit,et al.  Polygenic and directional regulatory evolution across pathways in Saccharomyces , 2010, Proceedings of the National Academy of Sciences.

[37]  K. Hokamp,et al.  Comparative Transcript Profiling of Candida albicans and Candida dubliniensis Identifies SFL2, a C. albicans Gene Required for Virulence in a Reconstituted Epithelial Infection Model , 2009, Eukaryotic Cell.

[38]  Brice Enjalbert,et al.  Glucose promotes stress resistance in the fungal pathogen Candida albicans. , 2009, Molecular biology of the cell.

[39]  Davis J. McCarthy,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[40]  A. Sellam,et al.  Transcriptional Regulation of Carbohydrate Metabolism in the Human Pathogen Candida albicans , 2009, PLoS pathogens.

[41]  Marie-Adèle Rajandream,et al.  Comparative genomics of the fungal pathogens Candida dubliniensis and Candida albicans. , 2009, Genome research.

[42]  S. Brunke,et al.  Identifying infection-associated genes of Candida albicans in the postgenomic era. , 2009, FEMS yeast research.

[43]  K. Kristiansen,et al.  SOAP2: an improved ultrafast tool for short read alignment , 2009, Bioinform..

[44]  Naama Barkai,et al.  A Yeast Hybrid Provides Insight into the Evolution of Gene Expression Regulation , 2009, Science.

[45]  Clifford A. Meyer,et al.  Model-based Analysis of ChIP-Seq (MACS) , 2008, Genome Biology.

[46]  D. MacCallum,et al.  Differential regulation of the transcriptional repressor NRG1 accounts for altered host‐cell interactions in Candida albicans and Candida dubliniensis , 2007, Molecular microbiology.

[47]  Alistair J. P. Brown,et al.  Effects of Depleting the Essential Central Metabolic Enzyme Fructose-1,6-Bisphosphate Aldolase on the Growth and Viability of Candida albicans: Implications for Antifungal Drug Target Discovery , 2006, Eukaryotic Cell.

[48]  G. Moran,et al.  Candida dubliniensis: ten years on. , 2005, FEMS microbiology letters.

[49]  S. Filler,et al.  Mice with disseminated candidiasis die of progressive sepsis. , 2005, The Journal of infectious diseases.

[50]  S. Noble,et al.  Strains and Strategies for Large-Scale Gene Deletion Studies of the Diploid Human Fungal Pathogen Candida albicans , 2005, Eukaryotic Cell.

[51]  Gary Moran,et al.  Comparative genomics using Candida albicans DNA microarrays reveals absence and divergence of virulence-associated genes in Candida dubliniensis. , 2004, Microbiology.

[52]  Gerald R. Fink,et al.  Transcriptional Response of Candida albicans upon Internalization by Macrophages , 2004, Eukaryotic Cell.

[53]  D. Soll,et al.  The Closely Related Species Candida albicans and Candida dubliniensis Can Mate , 2004, Eukaryotic Cell.

[54]  Andrew G. Clark,et al.  Evolutionary changes in cis and trans gene regulation , 2004, Nature.

[55]  John D. Storey,et al.  Statistical significance for genomewide studies , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[56]  Tom H. Pringle,et al.  The human genome browser at UCSC. , 2002, Genome research.

[57]  G. Gilfillan,et al.  Candida dubliniensis: phylogeny and putative virulence factors. , 1998, Microbiology.

[58]  Y Hochberg,et al.  Mathematical formulae for the prediction of the residual beta cell function during the first two years of disease in children and adolescents with insulin-dependent diabetes mellitus. , 1995, Medical hypotheses.

[59]  D. Coleman,et al.  Candida dubliniensis sp. nov.: phenotypic and molecular characterization of a novel species associated with oral candidosis in HIV-infected individuals. , 1995, Microbiology.

[60]  A. Brown,et al.  Fluctuations in glycolytic mRNA levels during morphogenesis in Candida albicans reflect underlying changes in growth and are not a response to cellular dimorphism , 1994, Molecular microbiology.

[61]  D. Fell.  Metabolic control analysis: a survey of its theoretical and experimental development. , 1992, The Biochemical journal.

[62]  Q. Mitrovich,et al.  Genetics and molecular biology in Candida albicans. , 2010, Methods in enzymology.

[63]  A. Gillum,et al.  Isolation of the Candida albicans gene for orotidine-5′-phosphate decarboxylase by complementation of S. cerevisiae ura3 and E. coli pyrF mutations , 2004, Molecular and General Genetics MGG.

[64]  G. Moran,et al.  Comparison of the epidemiology, drug resistance mechanisms, and virulence of Candida dubliniensis and Candida albicans. , 2004, FEMS yeast research.

[65]  M. Miyaji,et al.  Pathogenicity and virulence of Candida dubliniensis: comparison with C. albicans. , 2002, Medical mycology.

[66]  RESEARCH ARTICLE Open Access Research article Analysis of gene evolution and metabolic pathways using the Candida Gene Order Browser , 2022 .