The evolutionary rate variation among genes of HOG-signaling pathway in yeast genomes

BackgroundResponses to extracellular stress are required for microbes to survive in changing environments. Although the stress response mechanisms have been characterized extensively, the evolution of stress response pathway remains poorly understood. Here, we studied the evolution of High Osmolarity Glycerol (HOG) pathway, one of the important osmotic stress response pathways, across 10 yeast species and underpinned the evolutionary forces acting on the pathway evolution.ResultsAlthough the HOG pathway is well conserved across the surveyed yeast species, the evolutionary rate of the genes in this pathway varied substantially among or within different lineages. The fast divergence of MSB2 gene indicates that this gene is subjected to positive selection. Moreover, transcription factors in HOG pathway tend to evolve more rapidly, but the genes in conserved MAPK cascade underwent stronger functional selection. Remarkably, the dN/dS values are negatively correlated with pathway position along HOG pathway from Sln1 (Sho1) to Hog1 for transmitting external signal into nuclear. The increased gradient of selective constraints from upstream to downstream genes suggested that the downstream genes are more pleiotropic, being required for a wider range of pathways. In addition, protein length, codon usage, gene expression, and protein interaction appear to be important factors to determine the evolution of genes in HOG pathway.ConclusionsTaken together, our results suggest that functional constraints play a large role in the evolutionary rate variation in HOG pathway, but the genetic variation was influenced by quite complicated factors, such as pathway position, protein length and so on. These findings provide some insights into how HOG pathway genes evolved rapidly for responding to environmental osmotic stress changes.ReviewersThis article was reviewed by Han Liang (nominated by Laura Landweber), Georgy Bazykin (nominated by Mikhail Gelfand) and Zhenguo Lin (nominated by John Logsdon).

[1]  S. Ge,et al.  Evolutionary rate patterns of the Gibberellin pathway genes , 2009, BMC Evolutionary Biology.

[2]  A. Oudenaarden,et al.  A Systems-Level Analysis of Perfect Adaptation in Yeast Osmoregulation , 2009, Cell.

[3]  A. Gasch.  Comparative genomics of the environmental stress response in ascomycete fungi , 2007, Yeast.

[4]  Rodrigo Lopez,et al.  Clustal W and Clustal X version 2.0 , 2007, Bioinform..

[5]  Nick Goldman,et al.  Accuracy and Power of Statistical Methods for Detecting Adaptive Evolution in Protein Coding Sequences and for Identifying Positively Selected Sites , 2004, Genetics.

[6]  Zhan Zhou,et al.  The rapid evolution of signal peptides is mainly caused by relaxed selection on non-synonymous and synonymous sites. , 2009, Gene.

[7]  Ziheng Yang.  PAML 4: phylogenetic analysis by maximum likelihood. , 2007, Molecular biology and evolution.

[8]  J. Thorner,et al.  Function and regulation in MAPK signaling pathways: lessons learned from the yeast Saccharomyces cerevisiae. , 2007, Biochimica et biophysica acta.

[9]  Timothy B Sackton,et al.  Mutations in smooth muscle α-actin (ACTA2) lead to thoracic aortic aneurysms and dissections , 2007, Nature Genetics.

[10]  B. Birren,et al.  Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae , 2004, Nature.

[11]  Wen-Hsiung Li,et al.  Fundamentals of molecular evolution , 1990 .

[12]  N. Goldman,et al.  A codon-based model of nucleotide substitution for protein-coding DNA sequences. , 1994, Molecular biology and evolution.

[13]  Ira Herskowitz,et al.  Yeast go the whole HOG for the hyperosmotic response. , 2002, Trends in genetics : TIG.

[14]  J. Thorner,et al.  When the Stress of Your Environment Makes You Go HOG Wild , 2004, Science.

[15]  M. Aguadé,et al.  Network-level molecular evolutionary analysis of the insulin/TOR signal transduction pathway across 12 Drosophila genomes. , 2008, Genome research.

[16]  Marek S. Skrzypek,et al.  The Candida Genome Database (CGD), a community resource for Candida albicans gene and protein information , 2004, Nucleic Acids Res..

[17]  P. Phillips,et al.  High Nucleotide Divergence in Developmental Regulatory Genes Contrasts With the Structural Elements of Olfactory Pathways in Caenorhabditis , 2009, Genetics.

[18]  Marcus Krantz,et al.  Comparative genomics of the HOG-signalling system in fungi , 2006, Current Genetics.

[19]  T. A. Hall,et al.  BIOEDIT: A USER-FRIENDLY BIOLOGICAL SEQUENCE ALIGNMENT EDITOR AND ANALYSIS PROGRAM FOR WINDOWS 95/98/ NT , 1999 .

[20]  I. Herskowitz,et al.  Unique and redundant roles for HOG MAPK pathway components as revealed by whole-genome expression analysis. , 2003, Molecular biology of the cell.

[21]  H. Liang,et al.  Detecting positive selection in the budding yeast genome , 2009, Journal of evolutionary biology.

[22]  Kevin P. Byrne,et al.  The Yeast Gene Order Browser: combining curated homology and syntenic context reveals gene fate in polyploid species. , 2005, Genome research.

[23]  Eugene V Koonin,et al.  Comparable contributions of structural-functional constraints and expression level to the rate of protein sequence evolution , 2008, Biology Direct.

[24]  Patrick C Phillips,et al.  Evolutionary rates and centrality in the yeast gene regulatory network , 2009, Genome Biology.

[25]  Joseph Heitman,et al.  Sensing the environment: lessons from fungi , 2007, Nature Reviews Microbiology.

[26]  David Botstein,et al.  SGD: Saccharomyces Genome Database , 1998, Nucleic Acids Res..

[27]  E. Nevo,et al.  A MAPK gene from Dead Sea fungus confers stress tolerance to lithium salt and freezing-thawing: Prospects for saline agriculture. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[28]  A. Beskow,et al.  Comparative analysis of regulatory transcription factors in Schizosaccharomyces pombe and budding yeasts , 2006, Yeast.

[29]  R. E. Miller,et al.  Patterns of evolutionary rate variation among genes of the anthocyanin biosynthetic pathway. , 1999, Molecular biology and evolution.

[30]  S. Carroll,et al.  Genome-scale approaches to resolving incongruence in molecular phylogenies , 2003, Nature.

[31]  John P. Huelsenbeck,et al.  MrBayes 3: Bayesian phylogenetic inference under mixed models , 2003, Bioinform..

[32]  Megan N. McClean,et al.  Signal processing by the HOG MAP kinase pathway , 2008, Proceedings of the National Academy of Sciences.

[33]  N. Hannett,et al.  Activated Signal Transduction Kinases Frequently Occupy Target Genes , 2006, Science.

[34]  L. Knowles,et al.  Recombination rate and protein evolution in yeast , 2007, BMC Evolutionary Biology.

[35]  Adam P Arkin,et al.  Modularity of stress response evolution , 2008, Proceedings of the National Academy of Sciences.

[36]  Pär K Ingvarsson,et al.  Gene expression and protein length influence codon usage and rates of sequence evolution in Populus tremula. , 2006, Molecular biology and evolution.

[37]  B. Birren,et al.  Sequencing and comparison of yeast species to identify genes and regulatory elements , 2003, Nature.

[38]  I. Stansfield,et al.  Phylogenetic diversity of stress signalling pathways in fungi , 2009, BMC Evolutionary Biology.

[39]  L. Rieseberg,et al.  The correlation of evolutionary rate with pathway position in plant terpenoid biosynthesis. , 2009, Molecular biology and evolution.