Genome-wide signatures of population bottlenecks and diversifying selection in European wolves

Genomic resources developed for domesticated species provide powerful tools for studying the evolutionary history of their wild relatives. Here we use 61K single-nucleotide polymorphisms (SNPs) evenly spaced throughout the canine nuclear genome to analyse evolutionary relationships among the three largest European populations of grey wolves in comparison with other populations worldwide, and investigate genome-wide effects of demographic bottlenecks and signatures of selection. European wolves have a discontinuous range, with large and connected populations in Eastern Europe and relatively smaller, isolated populations in Italy and the Iberian Peninsula. Our results suggest a continuous decline in wolf numbers in Europe since the Late Pleistocene, and long-term isolation and bottlenecks in the Italian and Iberian populations following their divergence from the Eastern European population. The Italian and Iberian populations have low genetic variability and high linkage disequilibrium, but relatively few autozygous segments across the genome. This last characteristic clearly distinguishes them from populations that underwent recent drastic demographic declines or founder events, and implies long-term bottlenecks in these two populations. Although genetic drift due to spatial isolation and bottlenecks seems to be a major evolutionary force diversifying the European populations, we detected 35 loci that are putatively under diversifying selection. Two of these loci flank the canine platelet-derived growth factor gene, which affects bone growth and may influence differences in body size between wolf populations. This study demonstrates the power of population genomics for identifying genetic signals of demographic bottlenecks and detecting signatures of directional selection in bottlenecked populations, despite their low background variability.

[1]  J. Ozoliņš Status of Large Carnivore Conservation in the Baltic States , 2002 .

[2]  J. Kusak,et al.  High Genetic Variability of the Grey Wolf (Canis lupus L.) Population from Croatia as Revealed by Mitochondrial DNA Control Region Sequences , 2010 .

[3]  Huanming Yang,et al.  Whole-genome sequencing of giant pandas provides insights into demographic history and local adaptation , 2012, Nature Genetics.

[4]  I. Barja Prey and Prey-Age Preference by the Iberian Wolf Canis Lupus Signatus in a Multiple-Prey Ecosystem , 2009 .

[5]  M. Gompper Top Carnivores in the Suburbs? Ecological and Conservation Issues Raised by Colonization of North eastern North America by Coyotes , 2002 .

[6]  L. Mech,et al.  Premature Reproductive Activity in Wild Wolves , 1987 .

[7]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[8]  Marti J. Anderson,et al.  Climate and habitat barriers to dispersal in the highly mobile grey wolf , 2004, Molecular ecology.

[9]  G. Evanno,et al.  Detecting the number of clusters of individuals using the software structure: a simulation study , 2005, Molecular ecology.

[10]  Luigi Boitani,et al.  Wolves : behavior, ecology, and conservation , 2003 .

[11]  P. Visscher,et al.  Human population dispersal "Out of Africa" estimated from linkage disequilibrium and allele frequencies of SNPs. , 2011, Genome research.

[12]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[13]  Anthony D Whetton,et al.  THOC5/FMIP, an mRNA export TREX complex protein, is essential for hematopoietic primitive cell survival in vivo , 2010, BMC Biology.

[14]  N. Benecke,et al.  Late-Pleistocene and early Holocene history of the canid fauna of Europe (Canidae) , 2005 .

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

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

[17]  K. Crandall,et al.  Mitochondrial DNA phylogeography and population history of the grey wolf Canis lupus , 1999, Molecular ecology.

[18]  Theunis Piersma,et al.  The interplay between habitat availability and population differentiation , 2012 .

[19]  V. Urios,et al.  Signatures of demographic bottlenecks in European wolf populations , 2011, Conservation Genetics.

[20]  M. Hofreiter,et al.  Fossil dogs and wolves from Palaeolithic sites in Belgium, the Ukraine and Russia: osteometry, ancient DNA and stable isotopes , 2009 .

[21]  J. Pauli,et al.  Reliability of genetic bottleneck tests for detecting recent population declines , 2012, Molecular ecology.

[22]  Bertrand Servin,et al.  Genome-Wide Analysis of the World's Sheep Breeds Reveals High Levels of Historic Mixture and Strong Recent Selection , 2012, PLoS biology.

[23]  A. Kitchener,et al.  Wolf Canis lupus , 2008 .

[24]  Lucie M. Gattepaille,et al.  Inferring population size changes with sequence and SNP data: lessons from human bottlenecks , 2013, Heredity.

[25]  R. Nowak,et al.  The systematic status of the Talian wolf Canis lupus , 2002 .

[26]  L. Steinstraesser,et al.  Regulation of human adipose-derived stromal cell osteogenic differentiation by insulin-like growth factor-1 and platelet-derived growth factor-α. , 2011, Plastic and reconstructive surgery.

[27]  S. Mwaiko,et al.  Population genomic signatures of divergent adaptation, gene flow and hybrid speciation in the rapid radiation of Lake Victoria cichlid fishes , 2012, Molecular ecology.

[28]  K. Lindblad-Toh,et al.  The genomic signature of dog domestication reveals adaptation to a starch-rich diet , 2013, Nature.

[29]  T. C. Marshall,et al.  Statistical confidence for likelihood‐based paternity inference in natural populations , 1998, Molecular ecology.

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

[31]  D. Balding,et al.  Identifying adaptive genetic divergence among populations from genome scans , 2004, Molecular ecology.

[32]  Jeremiah D. Degenhardt,et al.  A Simple Genetic Architecture Underlies Morphological Variation in Dogs , 2010, PLoS biology.

[33]  H. Kantarjian,et al.  Biology of platelet-derived growth factor and its involvement in disease. , 2006, Mayo Clinic proceedings.

[34]  Jeremiah D. Degenhardt,et al.  Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication , 2010, Nature.

[35]  D. Reich,et al.  Population Structure and Eigenanalysis , 2006, PLoS genetics.

[36]  M. Apollonio,et al.  Prey selection and dietary response by wolves in a high-density multi-species ungulate community , 2011, European Journal of Wildlife Research.

[37]  R. Wayne,et al.  Megafaunal Extinctions and the Disappearance of a Specialized Wolf Ecomorph , 2007, Current Biology.

[38]  G. Luikart,et al.  Genomics and the future of conservation genetics , 2010, Nature Reviews Genetics.

[39]  S. Piertney,et al.  Selection maintains MHC diversity through a natural population bottleneck. , 2012, Molecular biology and evolution.

[40]  K. Schmidt,et al.  Genetic diversity and relatedness within packs in an intensely hunted population of wolvesCanis lupus , 2005, Acta Theriologica.

[41]  Jean-Pierre Bocquet-Appel,et al.  When the World’s Population Took Off: The Springboard of the Neolithic Demographic Transition , 2011, Science.

[42]  Top Carnivores in the Suburbs ? , 2002 .

[43]  M. Gautier,et al.  Inferring population histories using genome-wide allele frequency data. , 2013, Molecular biology and evolution.

[44]  F. Nociti,et al.  Platelet-derived growth factor/insulin-like growth factor-1 combination and bone regeneration around implants placed into extraction sockets: a histometric study in dogs. , 2000, Implant dentistry.

[45]  O. Gaggiotti,et al.  A Genome-Scan Method to Identify Selected Loci Appropriate for Both Dominant and Codominant Markers: A Bayesian Perspective , 2008, Genetics.

[46]  N. Lowe The world’s population , 1997 .

[47]  M. Nei,et al.  MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. , 2011, Molecular biology and evolution.

[48]  L. Boitani,et al.  The world of wolves : new perspectives on ecology, behaviour and management , 2010 .

[49]  J. Wall,et al.  Gene conversion and different population histories may explain the contrast between polymorphism and linkage disequilibrium levels. , 2001, American journal of human genetics.

[50]  W. Miller,et al.  Polar and brown bear genomes reveal ancient admixture and demographic footprints of past climate change , 2012, Proceedings of the National Academy of Sciences.

[51]  S. Funk,et al.  Ecological factors influence population genetic structure of European grey wolves , 2006, Molecular ecology.

[52]  Dmitry A. Konovalov,et al.  kingroup: a program for pedigree relationship reconstruction and kin group assignments using genetic markers , 2004 .

[53]  J. Kusak,et al.  Genetic structure of expanding wolf (Canis lupus) populations in Italy and Croatia, and the early steps of the recolonization of the Eastern Alps , 2014 .

[54]  D. Demontis,et al.  North-South Differentiation and a Region of High Diversity in European Wolves (Canis lupus) , 2013, PloS one.

[55]  Roland Kays,et al.  A genome-wide perspective on the evolutionary history of enigmatic wolf-like canids. , 2011, Genome research.

[56]  Ryan D. Hernandez,et al.  Genome-Wide Patterns of Nucleotide Polymorphism in Domesticated Rice , 2007, PLoS genetics.

[57]  M. Longaker,et al.  Regulation of Human Adipose-Derived Stromal Cell Osteogenic Differentiation by Insulin-Like Growth Factor-1 and Platelet-Derived Growth Factor-&agr; , 2010, Plastic and reconstructive surgery.

[58]  David H. Alexander,et al.  Fast model-based estimation of ancestry in unrelated individuals. , 2009, Genome research.

[59]  P. Visscher,et al.  Novel multilocus measure of linkage disequilibrium to estimate past effective population size. , 2003, Genome research.

[60]  J. Brigham‐Grette,et al.  Rapid sea-level rise and Holocene climate in the Chukchi Sea , 2006 .

[61]  F. Gugerli,et al.  Are adaptive loci transferable across genomes of related species? Outlier and environmental association analyses in Alpine Brassicaceae species , 2013, Molecular ecology.

[62]  Peter M Visscher,et al.  Recent human effective population size estimated from linkage disequilibrium. , 2007, Genome research.

[63]  Bridgett M. vonHoldt,et al.  STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method , 2011, Conservation Genetics Resources.

[64]  James A. Cuff,et al.  Genome sequence, comparative analysis and haplotype structure of the domestic dog , 2005, Nature.

[65]  E. Randi,et al.  Concordant mitochondrial and microsatellite DNA structuring between Polish lowland and Carpathian Mountain wolves , 2013, Conservation Genetics.

[66]  N. Cunningham,et al.  Bone induction by osteogenin and bone morphogenetic proteins. , 1990, Biomaterials.

[67]  V. Lucchini,et al.  Evidence of genetic distinction and long‐term population decline in wolves (Canis lupus) in the Italian Apennines , 2004, Molecular ecology.

[68]  S. Lynch,et al.  A review of the effects of insulin-like growth factor and platelet derived growth factor on in vivo cartilage healing and repair. , 2006, Osteoarthritis and cartilage.

[69]  Mukesh A. Patel,et al.  A Review on Effects , 2012 .

[70]  S. Kalinowski,et al.  The computer program STRUCTURE does not reliably identify the main genetic clusters within species: simulations and implications for human population structure , 2011, Heredity.

[71]  E. Ostrander,et al.  The IGF1 small dog haplotype is derived from Middle Eastern gray wolves , 2010, BMC Biology.

[72]  M. Hofreiter,et al.  Diversity lost: are all Holarctic large mammal species just relict populations? , 2010, BMC Biology.

[73]  Ilan Gronau,et al.  Genome Sequencing Highlights the Dynamic Early History of Dogs , 2014, PLoS genetics.

[74]  B. Jędrzejewska,et al.  Wolf Canis lupus numbers, diet and damage to livestock in relation to hunting and ungulate abundance in northeastern Belarus during 1990–2000 , 2003, Wildlife Biology.

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

[76]  H. Okarma,et al.  Population dynamics of Wolves Canis lupus in Bialowieża Primeval Forest (Poland and Belarus) in relation to hunting by humans, 1847–1993 , 1996 .

[77]  Baltic States Status of Large Carnivore Conservation in the , 2001 .

[78]  Lynn E Eberly,et al.  Body size and predatory performance in wolves: is bigger better? , 2009, The Journal of animal ecology.

[79]  E. Randi Genetics and conservation of wolves Canis lupus in Europe , 2011 .

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

[81]  W. Branicki,et al.  articlePhylogeographic history of grey wolves in Europe , 2010 .