A Single Interacting Species Leads to Widespread Parallel Evolution of the Stickleback Genome

Biotic interactions are potent, widespread causes of natural selection and divergent phenotypic evolution and can lead to genetic differentiation with gene flow among wild populations ("isolation by ecology") [1-4]. Biotic selection has been predicted to act on more genes than abiotic selection thereby driving greater adaptation [5]. However, difficulties in isolating the genome-wide effect of single biotic agents of selection have limited our ability to identify and quantify the number and type of genetic regions responding to biotic selection [6-9]. We identified geographically interspersed lakes in which threespine stickleback fish (Gasterosteus aculeatus) have repeatedly adapted to the presence or absence of a single member of the ecological community, prickly sculpin (Cottus asper), a fish that is both a competitor and a predator of the stickleback [10]. Whole-genome sequencing revealed that sculpin presence or absence accounted for the majority of genetic divergence among stickleback populations, more so than geography. The major axis of genomic variation within and between the two lake types was correlated with multiple traits, indicating parallel natural selection across a gradient of biotic environments. A large proportion of the genome-about 1.8%, encompassing more than 600 genes-differentiated stickleback from the two biotic environments. Divergence occurred in 141 discrete genomic clumps located mainly in regions of low recombination, suggesting that genes brought to lakes by the colonizing ancestral population often evolved together in linked blocks. Strong selection and a wealth of standing genetic variation explain how a single member of the biotic community can have such a rapid and profound evolutionary impact.

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

[2]  M. Arnegard,et al.  Genetics of ecological divergence during speciation , 2014, Nature.

[3]  B. Charlesworth WHY WE ARE NOT DEAD ONE HUNDRED TIMES OVER , 2013, Evolution; international journal of organic evolution.

[4]  Nathan J B Kraft,et al.  INTRAGUILD PREDATION DRIVES EVOLUTIONARY NICHE SHIFT IN THREESPINE STICKLEBACK , 2012, Evolution; international journal of organic evolution.

[5]  Richard Durbin,et al.  Fast and accurate long-read alignment with Burrows–Wheeler transform , 2010, Bioinform..

[6]  D. Rokhsar,et al.  Genome Assembly Improvement and Mapping Convergently Evolved Skeletal Traits in Sticklebacks with Genotyping-by-Sequencing , 2015, G3: Genes, Genomes, Genetics.

[7]  A. Meyer,et al.  Recurrent origin of a sexually selected trait in Xiphophorus fishes inferred from a molecular phylogeny , 1994, Nature.

[8]  Jeremy Schmutz,et al.  Widespread Parallel Evolution in Sticklebacks by Repeated Fixation of Ectodysplasin Alleles , 2005, Science.

[9]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[10]  W. Stephan Signatures of positive selection: from selective sweeps at individual loci to subtle allele frequency changes in polygenic adaptation , 2016, Molecular ecology.

[11]  J. Wolf,et al.  Widespread evidence for incipient ecological speciation: a meta-analysis of isolation-by-ecology. , 2013, Ecology letters.

[12]  C. Peichel,et al.  The genetic and molecular architecture of phenotypic diversity in sticklebacks , 2017, Philosophical Transactions of the Royal Society B: Biological Sciences.

[13]  Korbinian Strimmer,et al.  APE: Analyses of Phylogenetics and Evolution in R language , 2004, Bioinform..

[14]  Rebecca E. Irwin,et al.  Ecological and Evolutionary Consequences of Multispecies Plant-Animal Interactions , 2004 .

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

[16]  M. Grabherr,et al.  Evolution of Darwin’s finches and their beaks revealed by genome sequencing , 2015, Nature.

[17]  J. Endler Natural selection in the wild , 1987 .

[18]  Combining Experimental Evolution and Genomics to Understand How Seed Beetles Adapt to a Marginal Host Plant , 2020, Genes.

[19]  D. Absher,et al.  A Genome-wide SNP Genotyping Array Reveals Patterns of Global and Repeated Species-Pair Divergence in Sticklebacks , 2012, Current Biology.

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

[21]  Sharon I. Greenblum,et al.  Microbiome composition shapes rapid genomic adaptation of Drosophila melanogaster , 2019, Proceedings of the National Academy of Sciences.

[22]  D. Schluter,et al.  The genetic architecture of divergence between threespine stickleback species , 2001, Nature.

[23]  E. Birney,et al.  Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt , 2009, Nature Protocols.

[24]  Marius Roesti,et al.  Recombination in the threespine stickleback genome—patterns and consequences , 2013, Molecular ecology.

[25]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[26]  Tina T. Hu,et al.  Population resequencing reveals local adaptation of Arabidopsis lyrata to serpentine soils , 2010, Nature Genetics.

[27]  A. Kondrashov Contamination of the genome by very slightly deleterious mutations: why have we not died 100 times over? , 1995, Journal of theoretical biology.

[28]  Sara E. Miller,et al.  Intraguild predation leads to genetically based character shifts in the threespine stickleback , 2015, Evolution; international journal of organic evolution.

[29]  Jeffrey A. Walker,et al.  Ecological morphology of lacustrine threespine stickleback Gasterosteus aculeatus L. (Gasterosteidae) body shape , 1997 .

[30]  Chaochun Wei,et al.  Whole-genome sequencing of six dog breeds from continuous altitudes reveals adaptation to high-altitude hypoxia , 2014, Genome research.

[31]  M. Speekenbrink,et al.  depmixS4: An R Package for Hidden Markov Models , 2010 .

[32]  D. Schluter,et al.  GENETIC SIGNATURE OF ADAPTIVE PEAK SHIFT IN THREESPINE STICKLEBACK , 2012, Evolution; international journal of organic evolution.

[33]  J. Slate,et al.  Adaptation genomics: the next generation. , 2010, Trends in ecology & evolution.

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

[35]  L. Excoffier,et al.  Surfing during population expansions promotes genetic revolutions and structuration. , 2008, Trends in ecology & evolution.

[36]  Matteo Fumagalli,et al.  Inference of natural selection from ancient DNA , 2020, Evolution letters.

[37]  Paul Shannon,et al.  VariantAnnotation: a Bioconductor package for exploration and annotation of genetic variants , 2014, Bioinform..

[38]  Alexandros Stamatakis,et al.  RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies , 2014, Bioinform..

[39]  A. Kondrashov,et al.  Fast Evolution from Precast Bricks: Genomics of Young Freshwater Populations of Threespine Stickleback Gasterosteus aculeatus , 2014, PLoS genetics.

[40]  L. Rieseberg,et al.  Directional selection is the primary cause of phenotypic diversification , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[41]  C. Adams,et al.  Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids , 2019, BMC Genomics.

[42]  Michael Matschiner,et al.  Fitchi: haplotype genealogy graphs based on the Fitch algorithm , 2016, Bioinform..

[43]  J. M. Hoekstra,et al.  The Strength of Phenotypic Selection in Natural Populations , 2001, The American Naturalist.

[44]  S. Foster,et al.  Effects of predation risk on egg steroid profiles across multiple populations of threespine stickleback , 2020, Scientific Reports.

[45]  W. Cresko,et al.  Repeated Selection of Alternatively Adapted Haplotypes Creates Sweeping Genomic Remodeling in Stickleback , 2018, Genetics.

[46]  E. Taylor,et al.  Evolutionary history of an adaptive radiation in species pairs of threespine sticklebacks (Gasterosteus): insights from mitochondrial DNA , 1999 .

[47]  D. Hartl,et al.  Principles of population genetics , 1981 .

[48]  M. Lascoux,et al.  Ecological genomics of local adaptation , 2013, Nature Reviews Genetics.

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

[50]  S. Gavrilets,et al.  The genomic signature of parallel adaptation from shared genetic variation. , 2014, Molecular ecology.

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

[52]  T. Kocher,et al.  Structure and evolution of teleost mitochondrial control regions , 1995, Journal of Molecular Evolution.

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

[54]  D. Ebert,et al.  The genetic basis of resistance and matching-allele interactions of a host-parasite system: The Daphnia magna-Pasteuria ramosa model , 2017, PLoS genetics.

[55]  M. Fortin,et al.  Should the Mantel test be used in spatial analysis? , 2015 .

[56]  Joachim Selbig,et al.  pcaMethods - a bioconductor package providing PCA methods for incomplete data , 2007, Bioinform..

[57]  M. J. Clark,et al.  The effects of heavy metal contamination on the aquatic biota of Buttle lake and the Campbell river drainage (Canada) , 1985, Archives of environmental contamination and toxicology.