le Long-term balancing selection on chromosomal variants associated with crypsis in a stick insect

How polymorphisms are maintained within populations over long periods of time remains debated, because genetic drift and various forms of selection are expected to reduce variation. Here, we study the genetic architecture and maintenance of phenotypic morphs that confer crypsis in Timema cristinae stick insects, combining phenotypic information and genotyping-by-sequencing data from 1,360 samples across 21 populations. We find two highly divergent chromosomal variants that span megabases of sequence and are associated with colour polymorphism. We show that these variants exhibit strongly reduced effective recombination, are geographically widespread and probably diverged millions of generations ago. We detect heterokaryotype excess and signs of balancing selection acting on these variants through the species' history. A third chromosomal variant in the same genomic region likely evolved more recently from one of the two colour variants and is associated with dorsal pattern polymorphism. Our results suggest that large-scale genetic variation associated with crypsis has been maintained for long periods of time by potentially complex processes of balancing selection.

[1]  Philipp W. Messer,et al.  Heterozygote advantage as a natural consequence of adaptation in diploids , 2011, Proceedings of the National Academy of Sciences.

[2]  G. Coop,et al.  THE SIGNATURE OF POSITIVE SELECTION ON STANDING GENETIC VARIATION , 2005, Evolution; international journal of organic evolution.

[3]  K. Mockaitis,et al.  Ecological Genomics of Anopheles gambiae Along a Latitudinal Cline: A Population-Resequencing Approach , 2012, Genetics.

[4]  G. Abecasis,et al.  A note on exact tests of Hardy-Weinberg equilibrium. , 2005, American journal of human genetics.

[5]  J. Wall,et al.  Unusual haplotype structure at the proximal breakpoint of In(2L)t in a natural population of Drosophila melanogaster. , 1999, Genetics.

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

[7]  Dolph Schluter,et al.  Genetics and ecological speciation , 2009, Proceedings of the National Academy of Sciences.

[8]  Aaron A. Comeault,et al.  Genomic consequences of multiple speciation processes in a stick insect , 2012, Proceedings of the Royal Society B: Biological Sciences.

[9]  P. Schmid-Hempel,et al.  Variation in genomic recombination rates among animal taxa and the case of social insects , 2007, Heredity.

[10]  P. Andolfatto,et al.  Inversion polymorphisms and nucleotide variability in Drosophila. , 2001, Genetical research.

[11]  B. Charlesworth,et al.  Theoretical genetics of Batesian mimicry II. Evolution of supergenes. , 1975, Journal of theoretical biology.

[12]  J. A. Allen,et al.  Frequency-dependent selection by predators. , 1988, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[13]  R. Lawson,et al.  COLOR‐PATTERN VARIATION IN LAKE ERIE WATER SNAKES: THE ROLE OF GENE FLOW , 1995, Evolution; international journal of organic evolution.

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

[15]  Timothy E. Farkas,et al.  Experimental evidence for ecological selection on genome variation in the wild , 2013, Ecology letters.

[16]  Evolution of genetic diversity and human diseases , 2016, Russian Journal of Genetics.

[17]  L. Cook A two–stage model for Cepaea polymorphism , 1998 .

[18]  M. Kirkpatrick,et al.  A sequential coalescent algorithm for chromosomal inversions , 2013, Heredity.

[19]  C. Sandoval THE EFFECTS OF THE RELATIVE GEOGRAPHIC SCALES OF GENE FLOW AND SELECTION ON MORPH FREQUENCIES IN THE WALKING‐STICK TIMEMA CRISTINAE , 1994, Evolution; international journal of organic evolution.

[20]  Yuseob Kim,et al.  Emergence of long‐term balanced polymorphism under cyclic selection of spatially variable magnitude , 2015, Evolution; international journal of organic evolution.

[21]  Xiang Zhou,et al.  Polygenic Modeling with Bayesian Sparse Linear Mixed Models , 2012, PLoS genetics.

[22]  Stuart R. Dennis,et al.  Color phenotypes are under similar genetic control in two distantly related species of Timema stick insect , 2016, Evolution; international journal of organic evolution.

[23]  M. Quail,et al.  The industrial melanism mutation in British peppered moths is a transposable element , 2016, Nature.

[24]  Lewis G. Spurgin,et al.  How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings , 2010, Proceedings of the Royal Society B: Biological Sciences.

[25]  M. Wellenreuther,et al.  Sexual selection and genetic colour polymorphisms in animals , 2014, Molecular ecology.

[26]  S. Ho,et al.  Relaxed Phylogenetics and Dating with Confidence , 2006, PLoS biology.

[27]  N. Besansky,et al.  The Population Genomics of Trans-Specific Inversion Polymorphisms in Anopheles gambiae , 2009, Genetics.

[28]  B. Crespi,et al.  Host-plant adaptation drives the parallel evolution of reproductive isolation , 2002, Nature.

[29]  E. Betrán,et al.  Recombination and gene flux caused by gene conversion and crossing over in inversion heterokaryotypes. , 1997, Genetics.

[30]  J. Wakeley Coalescent Theory: An Introduction , 2008 .

[31]  J K Kelly,et al.  A test of neutrality based on interlocus associations. , 1997, Genetics.

[32]  Remco Bouckaert,et al.  Evolutionary Rates and Hbv: Issues of Rate Estimation with Bayesian Molecular Methods , 2013, Antiviral therapy.

[33]  A. Whibley,et al.  Genetic architecture and balancing selection: the life and death of differentiated variants , 2017, Molecular ecology.

[34]  M. Blaxter,et al.  RAD-Seq derived markers flank the shell colour and banding loci of the Cepaea nemoralis supergene , 2013, Molecular ecology.

[35]  P. Hedrick What is the evidence for heterozygote advantage selection? , 2012, Trends in ecology & evolution.

[36]  B. Crespi,et al.  Molecular phylogenetics of sexual and parthenogenetic Timema walking–sticks , 1998, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[37]  B. Charlesworth,et al.  Elements of Evolutionary Genetics , 2010 .

[38]  J. Wakeley,et al.  Estimating ancestral population parameters. , 1997, Genetics.

[39]  J. Pritchard,et al.  A Map of Recent Positive Selection in the Human Genome , 2006, PLoS biology.

[40]  Montgomery Slatkin,et al.  Linkage disequilibrium — understanding the evolutionary past and mapping the medical future , 2008, Nature Reviews Genetics.

[41]  Philip W. Hedrick,et al.  Genetic Polymorphism in Heterogeneous Environments: The Age of Genomics , 2006 .

[42]  P. Hedrick,et al.  Genetic Polymorphism in Heterogeneous Environments , 1976 .

[43]  M. Stephens,et al.  Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. , 2003, Genetics.

[44]  A. Drummond,et al.  Bayesian inference of population size history from multiple loci , 2008, BMC Evolutionary Biology.

[45]  A. Bond The Evolution of Color Polymorphism: Crypticity, Searching Images, and Apostatic Selection , 2007 .

[46]  D. Charlesworth,et al.  The genetics of inbreeding depression , 2009, Nature Reviews Genetics.

[47]  V. Gavrilovic,et al.  Allopatric genetic origins for sympatric host-plant shifts and race formation in Rhagoletis , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[48]  M. Nei,et al.  Mathematical model for studying genetic variation in terms of restriction endonucleases. , 1979, Proceedings of the National Academy of Sciences of the United States of America.

[49]  T. Clutton‐Brock,et al.  Life history trade-offs at a single locus maintain sexually selected genetic variation , 2013, Nature.

[50]  B. Charlesworth,et al.  The effects of local selection, balanced polymorphism and background selection on equilibrium patterns of genetic diversity in subdivided populations. , 1997, Genetical research.

[51]  José A. M. Borghans,et al.  Heterozygote advantage fails to explain the high degree of polymorphism of the MHC , 2004, Immunogenetics.

[52]  Dong Xie,et al.  BEAST 2: A Software Platform for Bayesian Evolutionary Analysis , 2014, PLoS Comput. Biol..

[53]  M. Kirkpatrick,et al.  Coalescent patterns for chromosomal inversions in divergent populations , 2012, Philosophical Transactions of the Royal Society B: Biological Sciences.

[54]  Aaron A. Comeault,et al.  Selection on a Genetic Polymorphism Counteracts Ecological Speciation in a Stick Insect , 2015, Current Biology.

[55]  M. Kirkpatrick,et al.  Chromosome Inversions, Local Adaptation and Speciation , 2017, Genetics.

[56]  A. Whibley,et al.  Evolution of dominance mechanisms at a butterfly mimicry supergene , 2014, Nature Communications.

[57]  M. Nei,et al.  THE BOTTLENECK EFFECT AND GENETIC VARIABILITY IN POPULATIONS , 1975, Evolution; international journal of organic evolution.

[58]  M. Wade,et al.  Molecular evolution near a two-locus balanced polymorphism. , 2000, Journal of Theoretical Biology.

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

[60]  S. Kerje,et al.  Structural genomic changes underlie alternative reproductive strategies in the ruff (Philomachus pugnax) , 2015, Nature Genetics.

[61]  M. Stephens,et al.  Inference of population structure using multilocus genotype data: dominant markers and null alleles , 2007, Molecular ecology notes.

[62]  B. Crespi,et al.  The evolution of geographic parthenogenesis in Timema walking‐sticks , 2002, Molecular ecology.

[63]  D. Reich,et al.  Principal components analysis corrects for stratification in genome-wide association studies , 2006, Nature Genetics.

[64]  J. Endler An overview of the relationships between mimicry and crypsis , 1981 .

[65]  J. Endler,et al.  The complex business of survival by aposematism. , 2005, Trends in ecology & evolution.

[66]  Shamil Sunyaev,et al.  A limited role for balancing selection. , 2005, Trends in genetics : TIG.

[67]  M. Przeworski,et al.  A new approach to estimate parameters of speciation models with application to apes. , 2007, Genome research.

[68]  H. A. Orr,et al.  The evolutionary genetics of speciation. , 1998, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[69]  Katalin Csill'ery,et al.  abc: an R package for approximate Bayesian computation (ABC) , 2011, 1106.2793.

[70]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[71]  N. Patterson,et al.  Estimating and interpreting FST: The impact of rare variants , 2013, Genome research.

[72]  M. Zody,et al.  The Rose-comb Mutation in Chickens Constitutes a Structural Rearrangement Causing Both Altered Comb Morphology and Defective Sperm Motility , 2012, PLoS genetics.

[73]  J. Sese,et al.  Sequencing of the genus Arabidopsis identifies a complex history of nonbifurcating speciation and abundant trans-specific polymorphism , 2016, Nature Genetics.

[74]  C. Clarke,et al.  The genetics of the mimetic butterfly Papilio polytes L. , 1972, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[75]  J. Slate,et al.  Heterozygote Advantage for Fecundity , 2006, PloS one.

[76]  P. Nosil Reproductive isolation caused by visual predation on migrants between divergent environments , 2004, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[77]  B. Clarke,et al.  Evidence for Apostatic Selection by Wild Passerines , 1968, Nature.

[78]  B. Crespi,et al.  MULTIPLE DIRECT TRANSITIONS FROM SEXUAL REPRODUCTION TO APOMICTIC PARTHENOGENESIS IN Timema STICK INSECTS , 2009, Evolution; international journal of organic evolution.

[79]  H. Levene,et al.  Genetic Equilibrium When More Than One Ecological Niche is Available , 1953, The American Naturalist.

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

[81]  C. Landry,et al.  Chromosomal variation segregates within incipient species and correlates with reproductive isolation , 2014, Molecular ecology.

[82]  A. Ruíz,et al.  Effect of inversion polymorphism on the neutral nucleotide variability of linked chromosomal regions in Drosophila. , 2000, Genetics.

[83]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.

[84]  G. Churchill,et al.  Properties of statistical tests of neutrality for DNA polymorphism data. , 1995, Genetics.

[85]  Peter L. Ralph,et al.  Patterns of Neutral Diversity Under General Models of Selective Sweeps , 2011, Genetics.

[86]  A. Fijarczyk,et al.  Detecting balancing selection in genomes: limits and prospects , 2015, Molecular ecology.

[87]  B. Weir Inferences about linkage disequilibrium. , 1979, Biometrics.

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

[89]  P. Pamilo,et al.  Associative overdominance, heterozygosity and fitness , 1998, Heredity.

[90]  C. Sandoval Persistence of a walking-stick population (Phasmatoptera: Timematodea) after a wildfire. , 2000 .

[91]  Ziheng Yang,et al.  Exploring uncertainty in the calibration of the molecular clock , 2012, Biology Letters.

[92]  Kevin R. Thornton,et al.  A New Approach for Using Genome Scans to Detect Recent Positive Selection in the Human Genome , 2007, PLoS biology.

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

[94]  P. Nosil,et al.  COUNTERACTING SELECTIVE REGIMES AND HOST PREFERENCE EVOLUTION IN ECOTYPES OF TWO SPECIES OF WALKING‐STICKS , 2005, Evolution; international journal of organic evolution.

[95]  M. Nordborg,et al.  The genealogy of sequences containing multiple sites subject to strong selection in a subdivided population. , 2003, Genetics.

[96]  A. Forsman,et al.  Inter-individual variation promotes ecological success of populations and species: evidence from experimental and comparative studies , 2016 .

[97]  Paul Scheet,et al.  Linkage Disequilibrium-Based Quality Control for Large-Scale Genetic Studies , 2008, PLoS genetics.

[98]  N. Barton,et al.  The effects of multilocus balancing selection on neutral variability. , 2002, Genetics.

[99]  J. Wingfield,et al.  A supergene determines highly divergent male reproductive morphs in the ruff , 2015, Nature Genetics.

[100]  Thomas Nagylaki,et al.  Polymorphism in multiallelic migration-selection models with dominance. , 2009, Theoretical population biology.

[101]  C. Alex Buerkle,et al.  Stick Insect Genomes Reveal Natural Selection’s Role in Parallel Speciation , 2014, Science.

[102]  M. Kirkpatrick,et al.  Chromosome inversions, adaptive cassettes and the evolution of species' ranges , 2015, Molecular ecology.

[103]  J. Hermisson,et al.  A general condition for adaptive genetic polymorphism in temporally and spatially heterogeneous environments. , 2014, Theoretical population biology.

[104]  C. Wiuf,et al.  The Probability and Chromosomal Extent of trans-specific Polymorphism , 2004, Genetics.

[105]  Candy Rowe,et al.  Cognition and the evolution of camouflage , 2016, Proceedings of the Royal Society B: Biological Sciences.

[106]  L H. Rieseberg,et al.  Chromosomal rearrangements and speciation. , 2001, Trends in ecology & evolution.

[107]  M. Kirkpatrick How and Why Chromosome Inversions Evolve , 2010, PLoS biology.

[108]  J. Mckinnon,et al.  Linking color polymorphism maintenance and speciation. , 2007, Trends in ecology & evolution.

[109]  R. Nielsen,et al.  A Model-Based Approach for Identifying Signatures of Ancient Balancing Selection in Genetic Data , 2014, PLoS genetics.

[110]  S. Merilaita,et al.  Animal camouflage: current issues and new perspectives , 2009, Philosophical Transactions of the Royal Society B: Biological Sciences.

[111]  Sharon R Grossman,et al.  Detecting natural selection in genomic data. , 2013, Annual review of genetics.

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

[113]  B. Cassone,et al.  Localization of Candidate Regions Maintaining a Common Polymorphic Inversion (2La) in Anopheles gambiae , 2007, PLoS genetics.

[114]  A. Bond,et al.  Apostatic selection by blue jays produces balanced polymorphism in virtual prey , 1998, Nature.

[115]  P. Hedrick,et al.  PERSPECTIVE: DETECTING ADAPTIVE MOLECULAR POLYMORPHISM: LESSONS FROM THE MHC , 2003, Evolution; international journal of organic evolution.

[116]  S. Schaeffer Selection in Heterogeneous Environments Maintains the Gene Arrangement Polymorphism of Drosophila pseudoobscura , 2008, Evolution; international journal of organic evolution.

[117]  S. Schaeffer,et al.  Molecular Population Genetics of Inversion Breakpoint Regions in Drosophila pseudoobscura , 2013, G3: Genes, Genomes, Genetics.

[118]  D. Charlesworth Balancing Selection and Its Effects on Sequences in Nearby Genome Regions , 2006, PLoS genetics.

[119]  Russell B. Corbett-Detig,et al.  Population Genomics of Inversion Polymorphisms in Drosophila melanogaster , 2012, PLoS genetics.

[120]  R. Faria,et al.  Chromosomal speciation revisited: rearranging theory with pieces of evidence. , 2010, Trends in ecology & evolution.

[121]  D. Stuart-Fox,et al.  Accelerated speciation in colour-polymorphic birds , 2012, Nature.

[122]  Paul Scheet,et al.  A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. , 2006, American journal of human genetics.

[123]  Pardis C Sabeti,et al.  Detecting recent positive selection in the human genome from haplotype structure , 2002, Nature.

[124]  Aaron A. Comeault,et al.  Transitions between phases of genomic differentiation during stick-insect speciation , 2017, Nature Ecology &Evolution.

[125]  M. Lynch,et al.  Mutation Accumulation and the Extinction of Small Populations , 1995, The American Naturalist.

[126]  S. Otto,et al.  Evolution of sex: Resolving the paradox of sex and recombination , 2002, Nature Reviews Genetics.

[127]  Renaud Vitalis,et al.  rehh: an R package to detect footprints of selection in genome-wide SNP data from haplotype structure , 2012, Bioinform..

[128]  Peter Donnelly,et al.  Multiple Instances of Ancient Balancing Selection Shared Between Humans and Chimpanzees , 2013, Science.

[129]  Crespi,et al.  Phylogenetic evidence for the evolution of ecological specialization in Timema walking‐sticks , 2000 .

[130]  J. Krzywinski,et al.  Semipermeable species boundaries between Anopheles gambiae and Anopheles arabiensis: Evidence from multilocus DNA sequence variation , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[131]  J. Fordyce,et al.  GENOMIC REGIONS WITH A HISTORY OF DIVERGENT SELECTION AFFECT FITNESS OF HYBRIDS BETWEEN TWO BUTTERFLY SPECIES , 2012, Evolution; international journal of organic evolution.

[132]  F. Tajima Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. , 1989, Genetics.

[133]  A. Savitzky,et al.  Smoothing and Differentiation of Data by Simplified Least Squares Procedures. , 1964 .

[134]  M. E. Galassi,et al.  GNU SCIENTI C LIBRARY REFERENCE MANUAL , 2005 .

[135]  P. Hedrick,et al.  Negative‐assortative mating for color in wolves , 2016, Evolution; international journal of organic evolution.

[136]  N. Barton,et al.  ACCUMULATING POSTZYGOTIC ISOLATION GENES IN PARAPATRY: A NEW TWIST ON CHROMOSOMAL SPECIATION , 2003, Evolution; international journal of organic evolution.

[137]  R. Bouckaert,et al.  Looking for trees in the forest: summary tree from posterior samples , 2013, BMC Evolutionary Biology.

[138]  B. Crespi,et al.  Experimental evidence that predation promotes divergence in adaptive radiation. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[139]  R. Butlin,et al.  What explains rare and conspicuous colours in a snail? A test of time-series data against models of drift, migration or selection , 2016, Heredity.

[140]  D. Charlesworth The status of supergenes in the 21st century: recombination suppression in Batesian mimicry and sex chromosomes and other complex adaptations , 2015, Evolutionary applications.

[141]  R. Butlin Recombination and speciation , 2005, Molecular ecology.

[142]  A. D. de Roos,et al.  The Effect of Population Size and Recombination on Delayed Evolution of Polymorphism and Speciation in Sexual Populations , 2008, The American Naturalist.

[143]  C. Sandoval Differential visual predation on morphs of Timema cristinae (Phasmatodeae:Timemidae) and its consequences for host range , 1994 .

[144]  Jason E Stajich,et al.  Likelihoods From Summary Statistics: Recent Divergence Between Species , 2005, Genetics.

[145]  J. Kelly,et al.  Effects of Spatially Varying Selection on Nucleotide Diversity and Linkage Disequilibrium: Insights From Deer Mouse Globin Genes , 2008, Genetics.

[146]  S. Yeaman Genomic rearrangements and the evolution of clusters of locally adaptive loci , 2013, Proceedings of the National Academy of Sciences.

[147]  P. Hedrick Population genetics of malaria resistance in humans , 2011, Heredity.

[148]  Joachim Hermisson,et al.  Soft Sweeps III: The Signature of Positive Selection from Recurrent Mutation , 2006, PLoS genetics.

[149]  N L Kaplan,et al.  The coalescent process in models with selection and recombination. , 1988, Genetics.

[150]  M. Stephens,et al.  Interpreting principal component analyses of spatial population genetic variation , 2008, Nature Genetics.

[151]  L. Keller,et al.  A Y-like social chromosome causes alternative colony organization in fire ants , 2013, Nature.