1
|
Lucena-Perez M, Kleinman-Ruiz D, Marmesat E, Saveljev AP, Schmidt K, Godoy JA. Bottleneck-associated changes in the genomic landscape of genetic diversity in wild lynx populations. Evol Appl 2021; 14:2664-2679. [PMID: 34815746 PMCID: PMC8591332 DOI: 10.1111/eva.13302] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Revised: 08/17/2021] [Accepted: 09/08/2021] [Indexed: 01/06/2023] Open
Abstract
Demographic bottlenecks generally reduce genetic diversity through more intense genetic drift, but their net effect may vary along the genome due to the random nature of genetic drift and to local effects of recombination, mutation, and selection. Here, we analyzed the changes in genetic diversity following a bottleneck by comparing whole-genome diversity patterns in populations with and without severe recent documented declines of Iberian (Lynx pardinus, n = 31) and Eurasian lynx (Lynx lynx, n = 29). As expected, overall genomic diversity correlated negatively with bottleneck intensity and/or duration. Correlations of genetic diversity with divergence, chromosome size, gene or functional site content, GC content, or recombination were observed in nonbottlenecked populations, but were weaker in bottlenecked populations. Also, functional features under intense purifying selection and the X chromosome showed an increase in the observed density of variants, even resulting in higher θ W diversity than in nonbottlenecked populations. Increased diversity seems to be related to both a higher mutational input in those regions creating a large collection of low-frequency variants, a few of which increase in frequency during the bottleneck to the point they become detectable with our limited sample, and the reduced efficacy of purifying selection, which affects not only protein structure and function but also the regulation of gene expression. The results of this study alert to the possible reduction of fitness and adaptive potential associated with the genomic erosion in regulatory elements. Further, the detection of a gain of diversity in ultra-conserved elements can be used as a sensitive and easy-to-apply signature of genetic erosion in wild populations.
Collapse
Affiliation(s)
- Maria Lucena-Perez
- Departamento de Ecología Integrativa Estación Biológica de Doñana (CSIC) Sevilla Spain
| | - Daniel Kleinman-Ruiz
- Departamento de Ecología Integrativa Estación Biológica de Doñana (CSIC) Sevilla Spain
- Departamento de Genética Facultad de Biología Universidad Complutense Madrid Spain
| | - Elena Marmesat
- Departamento de Ecología Integrativa Estación Biológica de Doñana (CSIC) Sevilla Spain
| | - Alexander P Saveljev
- Department of Animal Ecology Russian Research Institute of Game Management and Fur Farming Kirov Russia
| | - Krzysztof Schmidt
- Mammal Research Institute Polish Academy of Sciences Białowieża Poland
| | - José A Godoy
- Departamento de Ecología Integrativa Estación Biológica de Doñana (CSIC) Sevilla Spain
| |
Collapse
|
2
|
Lucena-Perez M, Marmesat E, Kleinman-Ruiz D, Martínez-Cruz B, Węcek K, Saveljev AP, Seryodkin IV, Okhlopkov I, Dvornikov MG, Ozolins J, Galsandorj N, Paunovic M, Ratkiewicz M, Schmidt K, Godoy JA. Genomic patterns in the widespread Eurasian lynx shaped by Late Quaternary climatic fluctuations and anthropogenic impacts. Mol Ecol 2020; 29:812-828. [PMID: 31995648 PMCID: PMC7064982 DOI: 10.1111/mec.15366] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2019] [Revised: 12/27/2019] [Accepted: 01/16/2020] [Indexed: 12/28/2022]
Abstract
Disentangling the contribution of long-term evolutionary processes and recent anthropogenic impacts to current genetic patterns of wildlife species is key to assessing genetic risks and designing conservation strategies. Here, we used 80 whole nuclear genomes and 96 mitogenomes from populations of the Eurasian lynx covering a range of conservation statuses, climatic zones and subspecies across Eurasia to infer the demographic history, reconstruct genetic patterns, and discuss the influence of long-term isolation and/or more recent human-driven changes. Our results show that Eurasian lynx populations shared a common history until 100,000 years ago, when Asian and European populations started to diverge and both entered a period of continuous and widespread decline, with western populations, except Kirov, maintaining lower effective sizes than eastern populations. Population declines and increased isolation in more recent times probably drove the genetic differentiation between geographically and ecologically close westernmost European populations. By contrast, and despite the wide range of habitats covered, populations are quite homogeneous genetically across the Asian range, showing a pattern of isolation by distance and providing little genetic support for the several proposed subspecies. Mitogenomic and nuclear divergences and population declines starting during the Late Pleistocene can be mostly attributed to climatic fluctuations and early human influence, but the widespread and sustained decline since the Holocene is more probably the consequence of anthropogenic impacts which intensified in recent centuries, especially in western Europe. Genetic erosion in isolated European populations and lack of evidence for long-term isolation argue for the restoration of lost population connectivity.
Collapse
Affiliation(s)
- Maria Lucena-Perez
- Department of Integrative Ecology, Estación Biológica de Doñana (CSIC), Seville, Spain
| | - Elena Marmesat
- Department of Integrative Ecology, Estación Biológica de Doñana (CSIC), Seville, Spain
| | - Daniel Kleinman-Ruiz
- Department of Integrative Ecology, Estación Biológica de Doñana (CSIC), Seville, Spain
| | - Begoña Martínez-Cruz
- School of Biological and Environmental Sciences, Liverpool John Moores University, Liverpool, UK
| | - Karolina Węcek
- Mammal Research Institute, Polish Academy of Sciences, Białowieża, Poland
| | - Alexander P Saveljev
- Department of Animal Ecology, Russian Research Institute of Game Management and Fur Farming, Kirov, Russia.,Biological Faculty of Moscow State University, Moscow, Russia
| | - Ivan V Seryodkin
- Laboratory of Ecology and Conservation of Animals, Pacific Institute of Geography of Far East Branch of Russian Academy of Sciences, Vladivostok, Russia.,Far Eastern Federal University, Vladivostok, Russia
| | - Innokentiy Okhlopkov
- Institute for Biological Problems of Cryolithozone, Siberian Division of the Russian Academy of Sciences, Yakutsk, Russia
| | - Mikhail G Dvornikov
- Department of Hunting Resources, Russian Research Institute of Game Management and Fur Farming, Kirov, Russia
| | - Janis Ozolins
- Department of Hunting and Wildlife Management, Latvijas Valsts mežzinātnes institūts "Silava", Salaspils, Latvia
| | - Naranbaatar Galsandorj
- Institute of General and Experimental Biology, Mongolian Academy of Science, Ulaanbaatar, Mongolia
| | | | | | - Krzysztof Schmidt
- Mammal Research Institute, Polish Academy of Sciences, Białowieża, Poland
| | - José A Godoy
- Department of Integrative Ecology, Estación Biológica de Doñana (CSIC), Seville, Spain
| |
Collapse
|
3
|
Casas-Marce M, Marmesat E, Soriano L, Martínez-Cruz B, Lucena-Perez M, Nocete F, Rodríguez-Hidalgo A, Canals A, Nadal J, Detry C, Bernáldez-Sánchez E, Fernández-Rodríguez C, Pérez-Ripoll M, Stiller M, Hofreiter M, Rodríguez A, Revilla E, Delibes M, Godoy JA. Spatiotemporal Dynamics of Genetic Variation in the Iberian Lynx along Its Path to Extinction Reconstructed with Ancient DNA. Mol Biol Evol 2018; 34:2893-2907. [PMID: 28962023 PMCID: PMC5850336 DOI: 10.1093/molbev/msx222] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
There is the tendency to assume that endangered species have been both genetically and demographically healthier in the past, so that any genetic erosion observed today was caused by their recent decline. The Iberian lynx (Lynx pardinus) suffered a dramatic and continuous decline during the 20th century, and now shows extremely low genome- and species-wide genetic diversity among other signs of genomic erosion. We analyze ancient (N = 10), historical (N = 245), and contemporary (N = 172) samples with microsatellite and mitogenome data to reconstruct the species' demography and investigate patterns of genetic variation across space and time. Iberian lynx populations transitioned from low but significantly higher genetic diversity than today and shallow geographical differentiation millennia ago, through a structured metapopulation with varying levels of diversity during the last centuries, to two extremely genetically depauperate and differentiated remnant populations by 2002. The historical subpopulations show varying extents of genetic drift in relation to their recent size and time in isolation, but these do not predict whether the populations persisted or went finally extinct. In conclusion, current genetic patterns were mainly shaped by genetic drift, supporting the current admixture of the two genetic pools and calling for a comprehensive genetic management of the ongoing conservation program. This study illustrates how a retrospective analysis of demographic and genetic patterns of endangered species can shed light onto their evolutionary history and this, in turn, can inform conservation actions.
Collapse
Affiliation(s)
- Mireia Casas-Marce
- Department of Integrative Ecology, Doñana Biological Station (EBD) Spanish National Research Council (CSIC), Sevilla, Spain
| | - Elena Marmesat
- Department of Integrative Ecology, Doñana Biological Station (EBD) Spanish National Research Council (CSIC), Sevilla, Spain
| | - Laura Soriano
- Department of Integrative Ecology, Doñana Biological Station (EBD) Spanish National Research Council (CSIC), Sevilla, Spain
| | - Begoña Martínez-Cruz
- Department of Integrative Ecology, Doñana Biological Station (EBD) Spanish National Research Council (CSIC), Sevilla, Spain
| | - Maria Lucena-Perez
- Department of Integrative Ecology, Doñana Biological Station (EBD) Spanish National Research Council (CSIC), Sevilla, Spain
| | - Francisco Nocete
- Grupo de Investigación MIDAS, Departamento Historia I (Prehistoria), Universidad de Huelva, Huelva, Spain
| | - Antonio Rodríguez-Hidalgo
- Department of Prehistory, Complutense University, Madrid, Spain.,IDEA (Instituto de Evolución en África), Madrid, Spain.,Equipo de Investigación Primeros Pobladores de Extremadura, Casa de la Cultura Rodríguez Moñino, Cáceres, Spain
| | - Antoni Canals
- Equipo de Investigación Primeros Pobladores de Extremadura, Casa de la Cultura Rodríguez Moñino, Cáceres, Spain.,IPHES, Institut Català de Paleoecologia Humana i Evolució Social (Catalan Institute of Human Paleoecology and Social Evolution), Tarragona, Spain.,Area de Prehistoria (Department of Prehistory), Universitat Rovira i Virgili (URV), Tarragona, Spain
| | - Jordi Nadal
- SERP, Departament de Prehistòria, Història Antiga i Arqueologia, Universitat de Barcelona, Barcelona, Spain
| | - Cleia Detry
- UNIARQ - Centro de Arqueologia da Faculdade de Letras da Universidade de Lisboa, Alameda da Universidade, Lisboa, Portugal
| | - Eloísa Bernáldez-Sánchez
- Laboratorio de Paleobiología, Instituto Andaluz del Patrimonio Histórico, Consejería de Cultura de la Junta de Andalucía, Seville, Spain
| | | | - Manuel Pérez-Ripoll
- Prehistory and Archaeology Department, University of Valencia, Valencia, Spain
| | - Mathias Stiller
- Department of Ecology and Evolutionary Biology, University of California Santa Cruz, Santa Cruz, CA
| | - Michael Hofreiter
- Evolutionary Adaptive Genomics, Institute of Biochemistry and Biology, Faculty of Mathematics and Life Sciences, University of Potsdam, Potsdam, Germany
| | - Alejandro Rodríguez
- Department of Conservation Biology, Doñana Biological Station (EBD), Spanish National Research Council (CSIC), Sevilla, Spain
| | - Eloy Revilla
- Department of Conservation Biology, Doñana Biological Station (EBD), Spanish National Research Council (CSIC), Sevilla, Spain
| | - Miguel Delibes
- Department of Conservation Biology, Doñana Biological Station (EBD), Spanish National Research Council (CSIC), Sevilla, Spain
| | - José A Godoy
- Department of Integrative Ecology, Doñana Biological Station (EBD) Spanish National Research Council (CSIC), Sevilla, Spain
| |
Collapse
|
4
|
Lucena-Perez M, Soriano L, López-Bao JV, Marmesat E, Fernández L, Palomares F, Godoy JA. Reproductive biology and genealogy in the endangered Iberian lynx: Implications for conservation. Mamm Biol 2018. [DOI: 10.1016/j.mambio.2017.11.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
5
|
Marmesat E, Schmidt K, Saveljev AP, Seryodkin IV, Godoy JA. Retention of functional variation despite extreme genomic erosion: MHC allelic repertoires in the Lynx genus. BMC Evol Biol 2017; 17:158. [PMID: 28676046 PMCID: PMC5496644 DOI: 10.1186/s12862-017-1006-z] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2017] [Accepted: 06/23/2017] [Indexed: 01/16/2023] Open
Abstract
BACKGROUND Demographic bottlenecks erode genetic diversity and may increase endangered species' extinction risk via decreased fitness and adaptive potential. The genetic status of species is generally assessed using neutral markers, whose dynamic can differ from that of functional variation due to selection. The MHC is a multigene family described as the most important genetic component of the mammalian immune system, with broad implications in ecology and evolution. The genus Lynx includes four species differing immensely in demographic history and population size, which provides a suitable model to study the genetic consequences of demographic declines: the Iberian lynx being an extremely bottlenecked species and the three remaining ones representing common and widely distributed species. We compared variation in the most variable exon of the MHCI and MHCII-DRB loci among the four species of the Lynx genus. RESULTS The Iberian lynx was characterised by lower number of MHC alleles than its sister species (the Eurasian lynx). However, it maintained most of the functional genetic variation at MHC loci present in the remaining and genetically healthier lynx species at all nucleotide, amino acid, and supertype levels. CONCLUSIONS Species-wide functional genetic diversity can be maintained even in the face of severe population bottlenecks, which caused devastating whole genome genetic erosion. This could be the consequence of divergent alleles being retained across paralogous loci, an outcome that, in the face of frequent gene conversion, may have been favoured by balancing selection.
Collapse
Affiliation(s)
- Elena Marmesat
- Department of Integrative Ecology, Estación Biológica de Doñana (CSIC), C/Américo Vespucio, 26, 41092, Sevilla, Spain
| | - Krzysztof Schmidt
- Mammal Research Institute, Polish Academy of Sciences, 17-230, Białowieża, Poland
| | - Alexander P Saveljev
- Department of Animal Ecology, Russian Research Institute of Game Management and Fur Farming, 79 Preobrazhenskaya Str, Kirov, 610000, Russia
| | - Ivan V Seryodkin
- Laboratory of Ecology and Conservation of Animals, Pacific Institute of Geography of Far East Branch of Russian Academy of Sciences, 7 Radio Street, Vladivostok, 690041, Russia
- Far Eastern Federal University, 8 Sukhanova Street, Vladivostok, 690091, Russia
| | - José A Godoy
- Department of Integrative Ecology, Estación Biológica de Doñana (CSIC), C/Américo Vespucio, 26, 41092, Sevilla, Spain.
| |
Collapse
|
6
|
Schwensow N, Mazzoni CJ, Marmesat E, Fickel J, Peacock D, Kovaliski J, Sinclair R, Cassey P, Cooke B, Sommer S. High adaptive variability and virus-driven selection on major histocompatibility complex (MHC) genes in invasive wild rabbits in Australia. Biol Invasions 2016. [DOI: 10.1007/s10530-016-1329-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
7
|
Marmesat E, Soriano L, Mazzoni CJ, Sommer S, Godoy JA. PCR Strategies for Complete Allele Calling in Multigene Families Using High-Throughput Sequencing Approaches. PLoS One 2016; 11:e0157402. [PMID: 27294261 PMCID: PMC4905633 DOI: 10.1371/journal.pone.0157402] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2016] [Accepted: 05/27/2016] [Indexed: 11/19/2022] Open
Abstract
The characterization of multigene families with high copy number variation is often approached through PCR amplification with highly degenerate primers to account for all expected variants flanking the region of interest. Such an approach often introduces PCR biases that result in an unbalanced representation of targets in high-throughput sequencing libraries that eventually results in incomplete detection of the targeted alleles. Here we confirm this result and propose two different amplification strategies to alleviate this problem. The first strategy (called pooled-PCRs) targets different subsets of alleles in multiple independent PCRs using different moderately degenerate primer pairs, whereas the second approach (called pooled-primers) uses a custom-made pool of non-degenerate primers in a single PCR. We compare their performance to the common use of a single PCR with highly degenerate primers using the MHC class I of the Iberian lynx as a model. We found both novel approaches to work similarly well and better than the conventional approach. They significantly scored more alleles per individual (11.33 ± 1.38 and 11.72 ± 0.89 vs 7.94 ± 1.95), yielded more complete allelic profiles (96.28 ± 8.46 and 99.50 ± 2.12 vs 63.76 ± 15.43), and revealed more alleles at a population level (13 vs 12). Finally, we could link each allele's amplification efficiency with the primer-mismatches in its flanking sequences and show that ultra-deep coverage offered by high-throughput technologies does not fully compensate for such biases, especially as real alleles may reach lower coverage than artefacts. Adopting either of the proposed amplification methods provides the opportunity to attain more complete allelic profiles at lower coverages, improving confidence over the downstream analyses and subsequent applications.
Collapse
Affiliation(s)
- Elena Marmesat
- Department of Integrative Ecology, Estación Biológica de Doñana (CSIC), Sevilla, Spain
| | - Laura Soriano
- Department of Integrative Ecology, Estación Biológica de Doñana (CSIC), Sevilla, Spain
| | - Camila J. Mazzoni
- Berlin Center for Genomics in Biodiversity Research (BeGenDiv), Berlin, Germany
- Leibniz Institute for Zoo and Wildlife Research (IZW), Berlin, Germany
| | - Simone Sommer
- Institute of Evolutionary Ecology and Conservation Genomics, University of Ulm, Ulm, Germany
| | - José A. Godoy
- Department of Integrative Ecology, Estación Biológica de Doñana (CSIC), Sevilla, Spain
- * E-mail:
| |
Collapse
|