Adaptive Evolution and Functional Redesign of Core Metabolic Proteins in Snakes

Background Adaptive evolutionary episodes in core metabolic proteins are uncommon, and are even more rarely linked to major macroevolutionary shifts. Methodology/Principal Findings We conducted extensive molecular evolutionary analyses on snake mitochondrial proteins and discovered multiple lines of evidence suggesting that the proteins at the core of aerobic metabolism in snakes have undergone remarkably large episodic bursts of adaptive change. We show that snake mitochondrial proteins experienced unprecedented levels of positive selection, coevolution, convergence, and reversion at functionally critical residues. We examined Cytochrome C oxidase subunit I (COI) in detail, and show that it experienced extensive modification of normally conserved residues involved in proton transport and delivery of electrons and oxygen. Thus, adaptive changes likely altered the flow of protons and other aspects of function in CO, thereby influencing fundamental characteristics of aerobic metabolism. We refer to these processes as “evolutionary redesign” because of the magnitude of the episodic bursts and the degree to which they affected core functional residues. Conclusions/Significance The evolutionary redesign of snake COI coincided with adaptive bursts in other mitochondrial proteins and substantial changes in mitochondrial genome structure. It also generally coincided with or preceded major shifts in ecological niche and the evolution of extensive physiological adaptations related to lung reduction, large prey consumption, and venom evolution. The parallel timing of these major evolutionary events suggests that evolutionary redesign of metabolic and mitochondrial function may be related to, or underlie, the extreme changes in physiological and metabolic efficiency, flexibility, and innovation observed in snake evolution.

[1]  Y. Kumazawa Mitochondrial DNA sequences of five squamates: phylogenetic affiliation of snakes. , 2004, DNA research : an international journal for rapid publication of reports on genes and genomes.

[2]  S. Eddy,et al.  tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. , 1997, Nucleic acids research.

[3]  C. Gissi,et al.  Nucleotide Substitution Rate of Mammalian Mitochondrial Genomes , 1999, Journal of Molecular Evolution.

[4]  J. A. McGuire,et al.  Comparative mitochondrial genomics of snakes: extraordinary substitution rate dynamics and functionality of the duplicate control region , 2007, BMC Evolutionary Biology.

[5]  Z. Yang,et al.  Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. , 2000, Molecular biology and evolution.

[6]  Jean L. Chang,et al.  Initial sequence of the chimpanzee genome and comparison with the human genome , 2005, Nature.

[7]  E. Louis,et al.  Molecular phylogenetics of squamata: the position of snakes, amphisbaenians, and dibamids, and the root of the squamate tree. , 2004, Systematic biology.

[8]  M. Nishida,et al.  Gene rearrangements in snake mitochondrial genomes: highly concerted evolution of control-region-like sequences duplicated and inserted into a tRNA gene cluster. , 1996, Molecular biology and evolution.

[9]  Ziheng Yang,et al.  PAML: a program package for phylogenetic analysis by maximum likelihood , 1997, Comput. Appl. Biosci..

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

[11]  Bronwen L. Aken,et al.  Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences , 2007, Nature.

[12]  J. Holman Fossil snakes of North America : origin, evolution, distribution, paleoecology , 2000 .

[13]  T. Tomizaki,et al.  Structures of metal sites of oxidized bovine heart cytochrome c oxidase at 2.8 A , 1995, Science.

[14]  W R Taylor,et al.  Coevolving protein residues: maximum likelihood identification and relationship to structure. , 1999, Journal of molecular biology.

[15]  M. Suematsu,et al.  The proton pumping pathway of bovine heart cytochrome c oxidase , 2007, Proceedings of the National Academy of Sciences.

[16]  宁北芳,et al.  疟原虫var基因转换速率变化导致抗原变异[英]/Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A , 2005 .

[17]  Y. Kumazawa,et al.  Complete Mitochondrial DNA Sequences of Six Snakes: Phylogenetic Relationships and Molecular Evolution of Genomic Features , 2005, Journal of Molecular Evolution.

[18]  S. Secor,et al.  Adaptive responses to feeding in Burmese pythons: pay before pumping. , 1995, The Journal of experimental biology.

[19]  R. Adkins,et al.  Evolution of eutherian cytochrome c oxidase subunit II: heterogeneous rates of protein evolution and altered interaction with cytochrome c. , 1996, Molecular biology and evolution.

[20]  S. Yoshikawa,et al.  The low-spin heme of cytochrome c oxidase as the driving element of the proton-pumping process , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[21]  R. Christensen,et al.  Physicochemical evolution and molecular adaptation of the cetacean and artiodactyl cytochrome b proteins. , 2005, Molecular biology and evolution.

[22]  Ruta,et al.  Nice snake, shame about the legs. , 2000, Trends in ecology & evolution.

[23]  A. F. Bennett,et al.  Physiology: Postprandial cardiac hypertrophy in pythons , 2005, Nature.

[24]  Vincent J. Lynch,et al.  Inventing an arsenal: adaptive evolution and neofunctionalization of snake venom phospholipase A2 genes , 2007, BMC Evolutionary Biology.

[25]  T. Andrews,et al.  Accelerated Evolution of Cytochrome b in Simian Primates: Adaptive Evolution in Concert with Other Mitochondrial Proteins? , 1998, Journal of Molecular Evolution.

[26]  Allan C. Wilson,et al.  Adaptive evolution in the stomach lysozymes of foregut fermenters , 1987, Nature.

[27]  T. Tomizaki,et al.  The Whole Structure of the 13-Subunit Oxidized Cytochrome c Oxidase at 2.8 Å , 1996, Science.

[28]  S. Hedges,et al.  The phylogeny of squamate reptiles (lizards, snakes, and amphisbaenians) inferred from nine nuclear protein-coding genes. , 2005, Comptes rendus biologies.

[29]  D. Wildman,et al.  Molecular evolution of aerobic energy metabolism in primates. , 2001, Molecular phylogenetics and evolution.

[30]  Zhengyuan O. Wang,et al.  Context dependence and coevolution among amino acid residues in proteins. , 2005, Methods in enzymology.

[31]  R. Gennis,et al.  Possible proton relay pathways in cytochrome c oxidase. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[32]  S. Muse,et al.  A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome. , 1994, Molecular biology and evolution.

[33]  D. Penny,et al.  Two new avian mitochondrial genomes (penguin and goose) and a summary of bird and reptile mitogenomic features. , 2003, Gene.

[34]  S. Secor,et al.  A vertebrate model of extreme physiological regulation , 1998, Nature.

[35]  S. Hedges,et al.  Early evolution of the venom system in lizards and snakes , 2006, Nature.

[36]  B. Schmidt,et al.  Role of the conserved arginine pair in proton and electron transfer in cytochrome C oxidase. , 2004, Biochemistry.

[37]  J. Faith,et al.  Evolution of base-substitution gradients in primate mitochondrial genomes. , 2005, Genome research.

[38]  B. Fry From genome to "venome": molecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins. , 2005, Genome research.

[39]  Harry W. Greene,et al.  Snakes: The Evolution of Mystery in Nature , 1997 .

[40]  S. Muse,et al.  Site-to-site variation of synonymous substitution rates. , 2005, Molecular biology and evolution.

[41]  J. Thompson,et al.  The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. , 1997, Nucleic acids research.

[42]  Sergei L. Kosakovsky Pond,et al.  HyPhy: hypothesis testing using phylogenies , 2005, Bioinform..

[43]  R. Gennis,et al.  Substitution of asparagine for aspartate-135 in subunit I of the cytochrome bo ubiquinol oxidase of Escherichia coli eliminates proton-pumping activity. , 1993, Biochemistry.

[44]  M. McCue Snakes survive starvation by employing supply- and demand-side economic strategies. , 2007, Zoology.

[45]  A. Dalziel,et al.  Molecular Evolution of Cytochrome c Oxidase in High-Performance Fish (Teleostei: Scombroidei) , 2006, Journal of Molecular Evolution.

[46]  R. Nielsen,et al.  Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. , 2002, Molecular biology and evolution.