Genomic Basis of Convergent Island Phenotypes in Boa Constrictors

Abstract Convergent evolution is often documented in organisms inhabiting isolated environments with distinct ecological conditions and similar selective regimes. Several Central America islands harbor dwarf Boa populations that are characterized by distinct differences in growth, mass, and craniofacial morphology, which are linked to the shared arboreal and feast-famine ecology of these island populations. Using high-density RADseq data, we inferred three dwarf island populations with independent origins and demonstrate that selection, along with genetic drift, has produced both divergent and convergent molecular evolution across island populations. Leveraging whole-genome resequencing data for 20 individuals and a newly annotated Boa genome, we identify four genes with evidence of phenotypically relevant protein-coding variation that differentiate island and mainland populations. The known roles of these genes involved in body growth (PTPRS, DMGDH, and ARSB), circulating fat and cholesterol levels (MYLIP), and craniofacial development (DMGDH and ARSB) in mammals link patterns of molecular evolution with the unique phenotypes of these island forms. Our results provide an important genome-wide example for quantifying expectations of selection and convergence in closely related populations. We also find evidence at several genomic loci that selection may be a prominent force of evolutionary change—even for small island populations for which drift is predicted to dominate. Overall, while phenotypically convergent island populations show relatively few loci under strong selection, infrequent patterns of molecular convergence are still apparent and implicate genes with strong connections to convergent phenotypes.

[1]  B. Jayne Comparative morphology of the semispinalis‐spinalis muscle of snakes and correlations with locomotion and constriction , 1982, Journal of morphology.

[2]  R. Shine Arboreality in Snakes: Ecology of the Australian Elapid Genus Hoplocephalus , 1983 .

[3]  B. Weir,et al.  ESTIMATING F‐STATISTICS FOR THE ANALYSIS OF POPULATION STRUCTURE , 1984, Evolution; international journal of organic evolution.

[4]  D. Bramble,et al.  Functional vertebrate morphology , 1985 .

[5]  D. Shulman,et al.  The interrelationships of thyroid and growth hormones: effect of growth hormone releasing hormone in hypo- and hyperthyroid male rats. , 1986, Acta endocrinologica. Supplementum.

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

[7]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[8]  Z. Hochberg,et al.  Effects of hypo or hyper‐thyroidism on growth hormone‐binding protein , 1991, Clinical endocrinology.

[9]  Elizabeth J. Robertson,et al.  Role of insulin-like growth factors in embryonic and postnatal growth , 1993, Cell.

[10]  J. T. Collins,et al.  Snakes: Ecology and Behavior , 1993 .

[11]  R. Levins,et al.  The Island Syndrome in Rodent Populations , 1994, The Quarterly Review of Biology.

[12]  B. McNab Resource Use and the Survival of Land and Freshwater Vertebrates on Oceanic Islands , 1994, The American Naturalist.

[13]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[14]  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.

[15]  Timothy E. Kennedy,et al.  Neuroendocrine dysplasia in mice lacking protein tyrosine phosphatase σ , 1999, Nature Genetics.

[16]  International Human Genome Sequencing Consortium Initial sequencing and analysis of the human genome , 2001, Nature.

[17]  B. McNab,et al.  Minimizing energy expenditure facilitates vertebrate persistence on oceanic islands , 2002 .

[18]  S. Eddy,et al.  Automated de novo identification of repeat sequence families in sequenced genomes. , 2002, Genome research.

[19]  Harold H. Zakon,et al.  Convergent Evolution on the Molecular Level , 2002, Brain, Behavior and Evolution.

[20]  S. Asa,et al.  Pituitary, pancreatic and gut neuroendocrine defects in protein tyrosine phosphatase-sigma-deficient mice. , 2002, Molecular endocrinology.

[21]  Mario Stanke,et al.  Gene prediction with a hidden Markov model and a new intron submodel , 2003, ECCB.

[22]  S. Boback Body Size Evolution in Snakes: Evidence from Island Populations , 2003, Copeia.

[23]  M. Nachman,et al.  The genetic basis of adaptive melanism in pocket mice , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Ian Korf,et al.  Gene finding in novel genomes , 2004, BMC Bioinformatics.

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

[26]  M. Nachman,et al.  ADAPTIVE REPTILE COLOR VARIATION AND THE EVOLUTION OF THE MC1R GENE , 2004, Evolution; international journal of organic evolution.

[27]  Ana Martins,et al.  Clinical and biochemical study of 28 patients with mucopolysaccharidosis type VI , 2004, Clinical genetics.

[28]  B. Grant,et al.  CONVERGENT EVOLUTION OF DARWIN'S FINCHES CAUSED BY INTROGRESSIVE HYBRIDIZATION AND SELECTION , 2004, Evolution; international journal of organic evolution.

[29]  Cathy H. Wu,et al.  UniProt: the Universal Protein knowledgebase , 2004, Nucleic Acids Res..

[30]  Burkhard Morgenstern,et al.  AUGUSTUS: a web server for gene finding in eukaryotes , 2004, Nucleic Acids Res..

[31]  K. Patel,et al.  Wnts and the neural crest , 2005, Anatomy and Embryology.

[32]  J. Keogh,et al.  RAPID AND REPEATED ORIGIN OF INSULAR GIGANTISM AND DWARFISM IN AUSTRALIAN TIGER SNAKES , 2005, Evolution; international journal of organic evolution.

[33]  J. Jurka,et al.  Repbase Update, a database of eukaryotic repetitive elements , 2005, Cytogenetic and Genome Research.

[34]  A. Isogai,et al.  Self-incompatibility in plants. , 2005, Annual review of plant biology.

[35]  S. Boback Natural History and Conservation of Island Boas (Boa Constrictor) in Belize , 2005, Copeia.

[36]  Pavel A. Pevzner,et al.  De novo identification of repeat families in large genomes , 2005, ISMB.

[37]  A. Janke,et al.  A mitogenomic study on the phylogenetic position of snakes , 2006 .

[38]  S. Boback A Morphometric Comparison of Island and Mainland Boas (Boa Constrictor) in Belize , 2006, Copeia.

[39]  J Waxman,et al.  Regulation of prostate cell growth and morphogenesis by Dickkopf-3 , 2006, Oncogene.

[40]  Jerzy Jurka,et al.  Annotation, submission and screening of repetitive elements in Repbase: RepbaseSubmitter and Censor , 2006, BMC Bioinformatics.

[41]  P. Insel,et al.  A Single Amino Acid Mutation Contributes to Adaptive Beach Mouse Color Pattern , 2006, Science.

[42]  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.

[43]  C. Bustamante,et al.  A Single IGF1 Allele Is a Major Determinant of Small Size in Dogs , 2007, Science.

[44]  R. Nusse,et al.  Wnt signaling mediates regional specification in the vertebrate face , 2007, Development.

[45]  R. Reed,et al.  Ecology and Conservation of an Exploited Insular Population of Boa constrictor (Squamata: Biodae) on the Cayos Cochinos, Honduras. , 2007 .

[46]  A. Rambaut,et al.  BEAST: Bayesian evolutionary analysis by sampling trees , 2007, BMC Evolutionary Biology.

[47]  S. Boback,et al.  Body Size and Head Shape of Island Boa constrictor in Belize: Environmental Versus Genetic Contributions , 2007 .

[48]  G. Wray The evolutionary significance of cis-regulatory mutations , 2007, Nature Reviews Genetics.

[49]  D. Hillis,et al.  Toxin-resistant sodium channels: parallel adaptive evolution across a complete gene family. , 2008, Molecular biology and evolution.

[50]  L. Rugiero,et al.  Biology of the Boas and Pythons , 2008 .

[51]  Mario Cáceres,et al.  The Chromosomal Polymorphism Linked to Variation in Social Behavior in the White-Throated Sparrow (Zonotrichia albicollis) Is a Complex Rearrangement and Suppressor of Recombination , 2008, Genetics.

[52]  S. Carroll Evo-Devo and an Expanding Evolutionary Synthesis: A Genetic Theory of Morphological Evolution , 2008, Cell.

[53]  D. Frynta,et al.  Mitochondrial DNA Variation Reveals Recent Evolutionary History of Main Boa constrictor Clades , 2009, Zoological science.

[54]  E. Brodie,et al.  The evolutionary origins of beneficial alleles during the repeated adaptation of garter snakes to deadly prey , 2009, Proceedings of the National Academy of Sciences.

[55]  Ryan D. Hernandez,et al.  Inferring the Joint Demographic History of Multiple Populations from Multidimensional SNP Frequency Data , 2009, PLoS genetics.

[56]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[57]  S. Moyà-Solà,et al.  Physiological and life history strategies of a fossil large mammal in a resource-limited environment , 2009, Proceedings of the National Academy of Sciences.

[58]  T. Schöneberg,et al.  Molecular and functional basis of phenotypic convergence in white lizards at White Sands , 2009, Proceedings of the National Academy of Sciences.

[59]  H. Hoekstra,et al.  On the Origin and Spread of an Adaptive Allele in Deer Mice , 2009, Science.

[60]  Todd A. Castoe,et al.  Evidence for an ancient adaptive episode of convergent molecular evolution , 2009, Proceedings of the National Academy of Sciences.

[61]  M. Lovett,et al.  Comparative gene expression analysis of avian embryonic facial structures reveals new candidates for human craniofacial disorders. , 2010, Human molecular genetics.

[62]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[63]  Emmanuel Paradis,et al.  pegas: an R package for population genetics with an integrated-modular approach , 2010, Bioinform..

[64]  Nicholas Stiffler,et al.  Population Genomics of Parallel Adaptation in Threespine Stickleback using Sequenced RAD Tags , 2010, PLoS genetics.

[65]  Günter P. Wagner,et al.  Evolutionary novelties , 2010, Current Biology.

[66]  F. Hoffmann,et al.  Gene cooption and convergent evolution of oxygen transport hemoglobins in jawed and jawless vertebrates , 2010, Proceedings of the National Academy of Sciences.

[67]  M. DePristo,et al.  A framework for variation discovery and genotyping using next-generation DNA sequencing data , 2011, Nature Genetics.

[68]  Mark Yandell,et al.  MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects , 2011, BMC Bioinformatics.

[69]  Heng Li,et al.  A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data , 2011, Bioinform..

[70]  J. Sinsheimer,et al.  The N342S MYLIP polymorphism is associated with high total cholesterol and increased LDL receptor degradation in humans. , 2011, The Journal of clinical investigation.

[71]  Camilo Salazar,et al.  Chromosomal rearrangements maintain a polymorphic supergene controlling butterfly mimicry , 2011, Nature.

[72]  J. Losos,et al.  CONVERGENCE, ADAPTATION, AND CONSTRAINT , 2011, Evolution; international journal of organic evolution.

[73]  Samuel E. Fox,et al.  Discovery of Highly Divergent Repeat Landscapes in Snake Genomes Using High-Throughput Sequencing , 2011, Genome biology and evolution.

[74]  A. Amores,et al.  Stacks: Building and Genotyping Loci De Novo From Short-Read Sequences , 2011, G3: Genes | Genomes | Genetics.

[75]  Jacob D. Jaffe,et al.  The genome of the green anole lizard and a comparative analysis with birds and mammals , 2011, Nature.

[76]  H. Zakon Adaptive evolution of voltage-gated sodium channels: The first 800 million years , 2012, Proceedings of the National Academy of Sciences.

[77]  J. Miller,et al.  Predicting the Functional Effect of Amino Acid Substitutions and Indels , 2012, PloS one.

[78]  P. Andolfatto,et al.  Parallel Molecular Evolution in an Herbivore Community , 2012, Science.

[79]  E. Brodie,et al.  Constraint shapes convergence in tetrodotoxin-resistant sodium channels of snakes , 2012, Proceedings of the National Academy of Sciences.

[80]  Alex A. Pollen,et al.  The genomic basis of adaptive evolution in threespine sticklebacks , 2012, Nature.

[81]  Yongwook Choi,et al.  A fast computation of pairwise sequence alignment scores between a protein and a set of single-locus variants of another protein , 2012, BCB.

[82]  E. Dahl,et al.  Targeting the Wnt pathway in cancer: the emerging role of Dickkopf-3. , 2012, Biochimica et biophysica acta.

[83]  H. Hoekstra,et al.  Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species , 2012, PloS one.

[84]  David Bryant,et al.  Next-generation sequencing reveals phylogeographic structure and a species tree for recent bird divergences. , 2009, Molecular phylogenetics and evolution.

[85]  David L. Stern,et al.  The genetic causes of convergent evolution , 2013, Nature Reviews Genetics.

[86]  Mauricio O. Carneiro,et al.  From FastQ Data to High‐Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline , 2013, Current protocols in bioinformatics.

[87]  P. Stadler,et al.  MITOS: improved de novo metazoan mitochondrial genome annotation. , 2013, Molecular phylogenetics and evolution.

[88]  Arndt von Haeseler,et al.  NextGenMap: fast and accurate read mapping in highly polymorphic genomes , 2013, Bioinform..

[89]  J. J. Emerson,et al.  Comparative Sex Chromosome Genomics in Snakes: Differentiation, Evolutionary Strata, and Lack of Global Dosage Compensation , 2013, PLoS biology.

[90]  M. Bouchard,et al.  Inactivation of LAR family phosphatase genes Ptprs and Ptprf causes craniofacial malformations resembling Pierre-Robin sequence , 2013, Development.

[91]  R. Dudley,et al.  Repeated elevational transitions in hemoglobin function during the evolution of Andean hummingbirds , 2013, Proceedings of the National Academy of Sciences.

[92]  Mark V. Lomolino,et al.  Of mice and mammoths: generality and antiquity of the island rule , 2013 .

[93]  Jeffrey D. Jensen,et al.  Adaptive Evolution of Multiple Traits Through Multiple Mutations at a Single Gene , 2013, Science.

[94]  Drew R. Schield,et al.  The Burmese python genome reveals the molecular basis for extreme adaptation in snakes , 2013, Proceedings of the National Academy of Sciences.

[95]  Colin N. Dewey,et al.  De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis , 2013, Nature Protocols.

[96]  E. Heyer,et al.  The role of GHR and IGF1 genes in the genetic determination of African pygmies’ short stature , 2012, European Journal of Human Genetics.

[97]  Angel Amores,et al.  Stacks: an analysis tool set for population genomics , 2013, Molecular ecology.

[98]  Inanç Birol,et al.  Assemblathon 2: evaluating de novo methods of genome assembly in three vertebrate species , 2013, GigaScience.

[99]  Tanya M. Teslovich,et al.  Discovery and refinement of loci associated with lipid levels , 2013, Nature Genetics.

[100]  R. Dahiya,et al.  microRNA-183 is an oncogene targeting Dkk-3 and SMAD4 in prostate cancer , 2013, British Journal of Cancer.

[101]  Michael R. Sussman,et al.  Genomic basis for the convergent evolution of electric organs , 2014, Science.

[102]  P. Trainor,et al.  Disrupting hedgehog and WNT signaling interactions promotes cleft lip pathogenesis. , 2014, The Journal of clinical investigation.

[103]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[104]  E. Rosenblum,et al.  The Molecular Basis of Phenotypic Convergence , 2014 .

[105]  L. Revell,et al.  Toward a Tree-of-Life for the boas and pythons: multilocus species-level phylogeny with unprecedented taxon sampling. , 2014, Molecular phylogenetics and evolution.

[106]  S. Edwards,et al.  Parallel Evolution of Tetrodotoxin Resistance in Three Voltage-Gated Sodium Channel Genes in the Garter Snake Thamnophis sirtalis , 2014, Molecular biology and evolution.

[107]  Wei Zhang,et al.  doublesex is a mimicry supergene , 2014, Nature.

[108]  Matthew Fraser,et al.  InterProScan 5: genome-scale protein function classification , 2014, Bioinform..

[109]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[110]  M. J. Thompson,et al.  Supergenes and their role in evolution , 2014, Heredity.

[111]  Zhihai Ma,et al.  Widespread contribution of transposable elements to the innovation of gene regulatory networks , 2014, Genome research.

[112]  F. Burbrink,et al.  Phylogeographical structure within Boa constrictor imperator across the lowlands and mountains of Central America and Mexico , 2014 .

[113]  K. Parsons,et al.  Wnt signalling underlies the evolution of new phenotypes and craniofacial variability in Lake Malawi cichlids , 2014, Nature Communications.

[114]  Sara M. Willems,et al.  The impact of low-frequency and rare variants on lipid levels , 2015, Nature Genetics.

[115]  Huaiyu Mi,et al.  The InterPro protein families database: the classification resource after 15 years , 2014, Nucleic Acids Res..

[116]  D. Hillis,et al.  Convergence of ion channel genome content in early animal evolution , 2015, Proceedings of the National Academy of Sciences.

[117]  O. Kohany,et al.  Repbase Update, a database of repetitive elements in eukaryotic genomes , 2015, Mobile DNA.

[118]  Cathy H. Wu,et al.  UniProt: the universal protein knowledgebase , 2016, Nucleic Acids Research.

[119]  Drew R. Schield,et al.  Phylogeographic and population genetic analyses reveal multiple species of Boa and independent origins of insular dwarfism. , 2016, Molecular phylogenetics and evolution.

[120]  F. Cunningham,et al.  The Ensembl Variant Effect Predictor , 2016, Genome Biology.

[121]  E. Tuttle,et al.  Divergence and Functional Degradation of a Sex Chromosome-like Supergene , 2016, Current Biology.

[122]  C. Witt,et al.  Predictable convergence in hemoglobin function has unpredictable molecular underpinnings , 2016, Science.

[123]  D. Absher,et al.  Genome-wide association study and targeted metabolomics identifies sex-specific association of CPS1 with coronary artery disease , 2016, Nature Communications.

[124]  M. Julião,et al.  Climbing , 2016, Palliative and Supportive Care.

[125]  S. Musani,et al.  Dietary choline and betaine; associations with subclinical markers of cardiovascular disease risk and incidence of CVD, coronary heart disease and stroke: the Jackson Heart Study , 2018, European Journal of Nutrition.

[126]  A. Whitehead,et al.  The genomic landscape of rapid repeated evolutionary adaptation to toxic pollution in wild fish , 2016, Science.

[127]  Meng-Pin Weng,et al.  modPhEA: model organism Phenotype Enrichment Analysis of eukaryotic gene sets , 2017, Bioinform..

[128]  A. Leaché,et al.  Evaluating mechanisms of diversification in a Guineo‐Congolian tropical forest frog using demographic model selection , 2017, Molecular ecology.

[129]  M. Pellegrini,et al.  Altered methylation of specific DNA loci in the liver of Bhmt‐null mice results in repression of Iqgap2 and F2rl2 and is associated with development of preneoplastic foci , 2017, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[130]  Drew R. Schield,et al.  Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes , 2017, Ecology and evolution.

[131]  S. Bhattacharyya,et al.  Chondroitin sulfatases differentially regulate Wnt signaling in prostate stem cells through effects on SHP2, phospho-ERK1/2, and Dickkopf Wnt signaling pathway inhibitor (DKK3) , 2017, Oncotarget.

[132]  Sara M. Schaal,et al.  Composite measures of selection can improve the signal‐to‐noise ratio in genome scans , 2017 .

[133]  C. Feschotte,et al.  Regulatory activities of transposable elements: from conflicts to benefits , 2016, Nature Reviews Genetics.

[134]  Sara M. Schaal,et al.  MINOTAUR: A platform for the analysis and visualization of multivariate results from genome scans with R Shiny , 2016, bioRxiv.

[135]  Richard H. Adams,et al.  GppFst: genomic posterior predictive simulations of FST and dXY for identifying outlier loci from population genomic data , 2016, Bioinform..

[136]  Drew R. Schield,et al.  Squamate reptiles challenge paradigms of genomic repeat element evolution set by birds and mammals , 2018, Nature Communications.

[137]  Judith A. Blake,et al.  Mouse Genome Database (MGD)-2018: knowledgebase for the laboratory mouse , 2017, Nucleic Acids Res..

[138]  U. Giger,et al.  Mucopolysaccharidosis Type VI in a Great Dane Caused by a Nonsense Mutation in the ARSB Gene , 2018, Veterinary pathology.

[139]  Jeffery P. Demuth,et al.  Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome , 2018, Genome biology and evolution.

[140]  Hyeonhae Choi,et al.  Effect of radioactive iodine-induced hypothyroidism on longitudinal bone growth during puberty in immature female rats , 2018, Experimental animals.

[141]  Richard H. Adams,et al.  ThetaMater: Bayesian estimation of population size parameter θ from genomic data , 2018, Bioinform..

[142]  J. Mallatt,et al.  Predictable convergence in hemoglobin function has unpredictable molecular underpinnings , 2022 .