1
|
Novo I, Ordás P, Moraga N, Santiago E, Quesada H, Caballero A. Impact of population structure in the estimation of recent historical effective population size by the software GONE. Genet Sel Evol 2023; 55:86. [PMID: 38049712 PMCID: PMC10694967 DOI: 10.1186/s12711-023-00859-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Accepted: 11/20/2023] [Indexed: 12/06/2023] Open
Abstract
BACKGROUND Effective population size (Ne) is a crucial parameter in conservation genetics and animal breeding. A recent method, implemented by the software GONE, has been shown to be rather accurate in estimating recent historical changes in Ne from a single sample of individuals. However, GONE estimations assume that the population being studied has remained isolated for a period of time, that is, without migration or confluence of other populations. If this occurs, the estimates of Ne can be heavily biased. In this paper, we evaluate the impact of migration and admixture on the estimates of historical Ne provided by GONE through a series of computer simulations considering several scenarios: (a) the mixture of two or more ancestral populations; (b) subpopulations that continuously exchange individuals through migration; (c) populations receiving migrants from a large source; and (d) populations with balanced systems of chromosomal inversions, which also generate genetic structure. RESULTS Our results indicate that the estimates of historical Ne provided by GONE may be substantially biased when there has been a recent mixture of populations that were previously separated for a long period of time. Similarly, biases may occur when the rate of continued migration between populations is low, or when chromosomal inversions are present at high frequencies. However, some biases due to population structuring can be eliminated by conducting population structure analyses and restricting the estimation to the differentiated groups. In addition, disregarding the genomic regions that are involved in inversions can also remove biases in the estimates of Ne. CONCLUSIONS Different kinds of deviations from isolation and panmixia of the populations can generate biases in the recent historical estimates of Ne. Therefore, estimation of past demography could benefit from performing population structure analyses beforehand, by mitigating the impact of these biases on historical Ne estimates.
Collapse
Affiliation(s)
- Irene Novo
- Centro de Investigación Mariña, Universidade de Vigo, Facultade de Bioloxía, 36310, Vigo, Spain.
| | - Pilar Ordás
- Centro de Investigación Mariña, Universidade de Vigo, Facultade de Bioloxía, 36310, Vigo, Spain
| | - Natalia Moraga
- Centro de Investigación Mariña, Universidade de Vigo, Facultade de Bioloxía, 36310, Vigo, Spain
| | - Enrique Santiago
- Departamento de Biología Funcional, Facultad de Biología, Universidad de Oviedo, 33006, Oviedo, Spain
| | - Humberto Quesada
- Centro de Investigación Mariña, Universidade de Vigo, Facultade de Bioloxía, 36310, Vigo, Spain
| | - Armando Caballero
- Centro de Investigación Mariña, Universidade de Vigo, Facultade de Bioloxía, 36310, Vigo, Spain
| |
Collapse
|
2
|
Duckett DJ, Calder K, Sullivan J, Tank DC, Carstens BC. Reduced representation approaches produce similar results to whole genome sequencing for some common phylogeographic analyses. PLoS One 2023; 18:e0291941. [PMID: 38032899 PMCID: PMC10688678 DOI: 10.1371/journal.pone.0291941] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2023] [Accepted: 09/09/2023] [Indexed: 12/02/2023] Open
Abstract
When designing phylogeographic investigations researchers can choose to collect many different types of molecular markers, including mitochondrial genes or genomes, SNPs from reduced representation protocols, large sequence capture data sets, and even whole genomes. Given that the statistical power and accuracy of various analyses are expected to differ depending on both the type of marker and the amount of data collected, an exploration of the variance across methodological results as a function of marker type should provide valuable information to researchers. Here we collect mitochondrial Cytochrome b sequences, whole mitochondrial genomes, single nucleotide polymorphisms (SNP)s isolated using a genotype by sequencing (GBS) protocol, sequences from ultraconserved elements, and low-coverage nuclear genomes from the North American water vole (Microtus richardsoni). We estimate genetic distances, population genetic structure, and historical demography using data from each of these datasets and compare the results across markers. As anticipated, the results exhibit differences across marker types, particularly in terms of the resolution offered by different analyses. A cost-benefit analysis indicates that SNPs collected using a GBS protocol are the most cost-effective molecular marker, with inferences that mirror those collected from the whole genome data at a fraction of the cost per sample.
Collapse
Affiliation(s)
- Drew J. Duckett
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University, Columbus, OH, United States of America
| | - Kailee Calder
- College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, United States of America
| | - Jack Sullivan
- Department of Biological Sciences, University of Idaho, Moscow, ID, United States of America
| | - David C. Tank
- Department of Botany, University of Wyoming, Laramie, WY, United States of America
| | - Bryan C. Carstens
- Department of Evolution, Ecology, and Organismal Biology, The Ohio State University, Columbus, OH, United States of America
| |
Collapse
|
3
|
Scerri EML. One species, many roots? Nat Ecol Evol 2023; 7:975-976. [PMID: 37198291 DOI: 10.1038/s41559-023-02080-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/19/2023]
Affiliation(s)
- Eleanor M L Scerri
- Pan-African Evolution Research Group, Max Planck Institute of Geoanthropology, Jena, Germany.
- Department of Prehistory, University of Cologne, Cologne, Germany.
- Department of Classics and Archaeology, University of Malta, Msida, Malta.
| |
Collapse
|
4
|
Scerri EML, Will M. The revolution that still isn't: The origins of behavioral complexity in Homo sapiens. J Hum Evol 2023; 179:103358. [PMID: 37058868 DOI: 10.1016/j.jhevol.2023.103358] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Revised: 03/13/2023] [Accepted: 03/14/2023] [Indexed: 04/16/2023]
Abstract
The behavioral origins of Homo sapiens can be traced back to the first material culture produced by our species in Africa, the Middle Stone Age (MSA). Beyond this broad consensus, the origins, patterns, and causes of behavioral complexity in modern humans remain debated. Here, we consider whether recent findings continue to support popular scenarios of: (1) a modern human 'package,' (2) a gradual and 'pan-African' emergence of behavioral complexity, and (3) a direct connection to changes in the human brain. Our geographically structured review shows that decades of scientific research have continuously failed to find a discrete threshold for a complete 'modernity package' and that the concept is theoretically obsolete. Instead of a continent-wide, gradual accumulation of complex material culture, the record exhibits a predominantly asynchronous presence and duration of many innovations across different regions of Africa. The emerging pattern of behavioral complexity from the MSA conforms to an intricate mosaic characterized by spatially discrete, temporally variable, and historically contingent trajectories. This archaeological record bears no direct relation to a simplistic shift in the human brain but rather reflects similar cognitive capacities that are variably manifested. The interaction of multiple causal factors constitutes the most parsimonious explanation driving the variable expression of complex behaviors, with demographic processes such as population structure, size, and connectivity playing a key role. While much emphasis has been given to innovation and variability in the MSA record, long periods of stasis and a lack of cumulative developments argue further against a strictly gradualistic nature in the record. Instead, we are confronted with humanity's deep, variegated roots in Africa, and a dynamic metapopulation that took many millennia to reach the critical mass capable of producing the ratchet effect commonly used to define contemporary human culture. Finally, we note a weakening link between 'modern' human biology and behavior from around 300 ka ago.
Collapse
Affiliation(s)
- Eleanor M L Scerri
- Pan-African Evolution Research Group, Max Planck Institute for Geoanthropology, Kahlaische Str. 10, 07749, Jena, Germany; Department of Classics and Archaeology, University of Malta, Msida, MSD 2080, Malta; Department of Prehistory, University of Cologne, 50931, Cologne, Germany.
| | - Manuel Will
- Department of Early Prehistory and Quaternary Ecology, University of Tübingen, Schloss Hohentübingen, Burgsteige 11, 72070, Tübingen, Germany
| |
Collapse
|
5
|
Lesturgie P, Braun CD, Clua E, Mourier J, Thorrold SR, Vignaud T, Planes S, Mona S. Like a rolling stone: Colonization and migration dynamics of the gray reef shark ( Carcharhinus amblyrhynchos). Ecol Evol 2023; 13:e9746. [PMID: 36644707 PMCID: PMC9831972 DOI: 10.1002/ece3.9746] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Revised: 12/18/2022] [Accepted: 12/27/2022] [Indexed: 01/13/2023] Open
Abstract
Designing appropriate management plans requires knowledge of both the dispersal ability and what has shaped the current distribution of the species under consideration. Here, we investigated the evolutionary history of the endangered gray reef shark (Carcharhinus amblyrhynchos) across its range by sequencing thousands of RADseq loci in 173 individuals in the Indo-Pacific (IP). We first bring evidence of the occurrence of a range expansion (RE) originating close to the Indo-Australian Archipelago (IAA) where two stepping-stone waves (east and westward) colonized almost the entire IP. Coalescent modeling additionally highlighted a homogenous connectivity (Nm ~ 10 per generation) throughout the range, and isolation by distance model suggested the absence of barriers to dispersal despite the affinity of C. amblyrhynchos to coral reefs. This coincides with long-distance swims previously recorded, suggesting that the strong genetic structure at the IP scale (F ST ~ 0.56 between its ends) is the consequence of its broad current distribution and organization in a large number of demes. Our results strongly suggest that management plans for the gray reef shark should be designed on a range-wide rather than a local scale due to its continuous genetic structure. We further contrasted these results with those obtained previously for the sympatric but strictly lagoon-associated Carcharhinus melanopterus, known for its restricted dispersal ability. Carcharhinus melanopterus exhibits a similar RE dynamic but is characterized by a stronger genetic structure and a nonhomogeneous connectivity largely dependent on local coral reefs availability. This sheds new light on shark evolution, emphasizing the roles of IAA as source of biodiversity and of life-history traits in shaping the extent of genetic structure and diversity.
Collapse
Affiliation(s)
- Pierre Lesturgie
- Institut de Systématique, Evolution, Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, EPHE‐PSLUniversité PSL, CNRS, SU, UAParisFrance
| | - Camrin D. Braun
- Biology DepartmentWoods Hole Oceanographic InstitutionWoods HoleMassachusettsUSA
| | - Eric Clua
- Laboratoire d'Excellence CORAILPapetoaiFrench Polynesia
- EPHE, PSL Research UniversityParisFrance
| | - Johann Mourier
- Laboratoire d'Excellence CORAILPapetoaiFrench Polynesia
- Université de Corse Pasquale Paoli, UMS 3514 Plateforme Marine Stella MareBigugliaFrance
| | - Simon R. Thorrold
- Biology DepartmentWoods Hole Oceanographic InstitutionWoods HoleMassachusettsUSA
| | | | - Serge Planes
- Laboratoire d'Excellence CORAILPapetoaiFrench Polynesia
- EPHE, PSL Research UniversityParisFrance
| | - Stefano Mona
- Institut de Systématique, Evolution, Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, EPHE‐PSLUniversité PSL, CNRS, SU, UAParisFrance
- EPHE, PSL Research UniversityParisFrance
| |
Collapse
|
6
|
Legried B, Terhorst J. Rates of convergence in the two-island and isolation-with-migration models. Theor Popul Biol 2022; 147:16-27. [PMID: 36007782 DOI: 10.1016/j.tpb.2022.08.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2021] [Revised: 08/10/2022] [Accepted: 08/11/2022] [Indexed: 11/25/2022]
Abstract
A number of powerful demographic inference methods have been developed in recent years, with the goal of fitting rich evolutionary models to genetic data obtained from many populations. In this paper we investigate the statistical performance of these methods in the specific case where there is continuous migration between populations. Compared with earlier work, migration significantly complicates the theoretical analysis and requires new techniques. We employ the theories of phase-type distributions and concentration of measure in order to study the two-island and isolation-with-migration models, resulting in both upper and lower bounds on rates of convergence for parametric estimators in migration models. For the upper bounds, we consider inferring rates of coalescent and migration on the basis of directly observing pairwise coalescent times, and, more realistically, when (conditionally) Poisson-distributed mutations dropped on latent trees are observed. We complement these upper bounds with information-theoretic lower bounds which establish a limit, in terms of sample size, below which inference is effectively impossible.
Collapse
Affiliation(s)
- Brandon Legried
- Department of Statistics, University of Michigan, United States of America
| | - Jonathan Terhorst
- Department of Statistics, University of Michigan, United States of America.
| |
Collapse
|
7
|
Boitard S, Arredondo A, Chikhi L, Mazet O. Heterogeneity in effective size across the genome: effects on the inverse instantaneous coalescence rate (IICR) and implications for demographic inference under linked selection. Genetics 2022; 220:6512058. [PMID: 35100421 PMCID: PMC8893248 DOI: 10.1093/genetics/iyac008] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Accepted: 01/01/2022] [Indexed: 01/22/2023] Open
Abstract
The relative contribution of selection and neutrality in shaping species genetic diversity is one of the most central and controversial questions in evolutionary theory. Genomic data provide growing evidence that linked selection, i.e. the modification of genetic diversity at neutral sites through linkage with selected sites, might be pervasive over the genome. Several studies proposed that linked selection could be modeled as first approximation by a local reduction (e.g. purifying selection, selective sweeps) or increase (e.g. balancing selection) of effective population size (Ne). At the genome-wide scale, this leads to variations of Ne from one region to another, reflecting the heterogeneity of selective constraints and recombination rates between regions. We investigate here the consequences of such genomic variations of Ne on the genome-wide distribution of coalescence times. The underlying motivation concerns the impact of linked selection on demographic inference, because the distribution of coalescence times is at the heart of several important demographic inference approaches. Using the concept of inverse instantaneous coalescence rate, we demonstrate that in a panmictic population, linked selection always results in a spurious apparent decrease of Ne along time. Balancing selection has a particularly large effect, even when it concerns a very small part of the genome. We also study more general models including genuine population size changes, population structure or transient selection and find that the effect of linked selection can be significantly reduced by that of population structure. The models and conclusions presented here are also relevant to the study of other biological processes generating apparent variations of Ne along the genome.
Collapse
Affiliation(s)
- Simon Boitard
- CBGP, Université de Montpellier, CIRAD, INRAE, Institut Agro, IRD, Montferrier-sur-Lez 34988, France
- Corresponding author: Université de Montpellier, CIRAD, INRAE, Institut Agro, IRD, 755 Avenue du Campus Agropolis, CS 30016, Montferrier-sur-Lez 34988, France.
| | - Armando Arredondo
- Institut National des Sciences Appliquées, Institut de Mathématiques de Toulouse, Université de Toulouse,Toulouse 31062, France
| | - Lounès Chikhi
- Instituto Gulbenkian de Ciência, Oeiras P-2780-156, Portugal
- Laboratoire Évolution & Diversité Biologique (EDB UMR 5174), CNRS, IRD, UPS, Université de Toulouse Midi-Pyrénées, Toulouse 31062, France
| | - Olivier Mazet
- Institut National des Sciences Appliquées, Institut de Mathématiques de Toulouse, Université de Toulouse,Toulouse 31062, France
| |
Collapse
|
8
|
Nadachowska‐Brzyska K, Konczal M, Babik W. Navigating the temporal continuum of effective population size. Methods Ecol Evol 2021. [DOI: 10.1111/2041-210x.13740] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
| | | | - Wieslaw Babik
- Jagiellonian University in Kraków Faculty of Biology Institute of Environmental Sciences Kraków Poland
| |
Collapse
|
9
|
Lesturgie P, Planes S, Mona S. Coalescence times, life history traits and conservation concerns: An example from four coastal shark species from the Indo-Pacific. Mol Ecol Resour 2021; 22:554-566. [PMID: 34407294 DOI: 10.1111/1755-0998.13487] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Revised: 07/27/2021] [Accepted: 08/12/2021] [Indexed: 11/30/2022]
Abstract
Dispersal abilities play a crucial role in shaping the extent of population genetic structure, with more mobile species being panmictic over large geographical ranges and less mobile ones organized in metapopulations exchanging migrants to different degrees. In turn, population structure directly influences the coalescence pattern of the sampled lineages, but the consequences on the estimated variation of the effective population size (Ne ) over time obtained by means of unstructured demographic models remain poorly understood. However, this knowledge is crucial for biologically interpreting the observed Ne trajectory and further devising conservation strategies in endangered species. Here we investigated the demographic history of four shark species (Carharhinus melanopterus, Carharhinus limbatus, Carharhinus amblyrhynchos, Galeocerdo cuvier) with different degrees of endangered status and life history traits related to dispersal distributed in the Indo-Pacific and sampled off New Caledonia. We compared several evolutionary scenarios representing both structured (metapopulation) and unstructured models and then inferred the Ne variation through time. By performing extensive coalescent simulations, we provided a general framework relating the underlying population structure and the observed Ne dynamics. On this basis, we concluded that the recent decline observed in three out of the four considered species when assuming unstructured demographic models can be explained by the presence of population structure. Furthermore, we also demonstrated the limits of the inferences based on the sole site frequency spectrum and warn that statistics based on linkage disequilibrium will be needed to exclude recent demographic events affecting meta-populations.
Collapse
Affiliation(s)
- Pierre Lesturgie
- Institut de Systématique, Evolution, Biodiversité, ISYEB (UMR 7205), Muséum National d'Histoire Naturelle, CNRS, Sorbonne Université, EPHE, Université des Antilles, Paris, France
| | - Serge Planes
- PSL Research University: EPHE-UPVD-CNRS, USR 3278 CRIOBE, Université de Perpignan, Perpignan, France.,Laboratoire d'Excellence CORAIL, Papetoai, French Polynesia
| | - Stefano Mona
- Institut de Systématique, Evolution, Biodiversité, ISYEB (UMR 7205), Muséum National d'Histoire Naturelle, CNRS, Sorbonne Université, EPHE, Université des Antilles, Paris, France.,Laboratoire d'Excellence CORAIL, Papetoai, French Polynesia.,EPHE, PSL Research University, Paris, France
| |
Collapse
|
10
|
Bourgeois YXC, Warren BH. An overview of current population genomics methods for the analysis of whole-genome resequencing data in eukaryotes. Mol Ecol 2021; 30:6036-6071. [PMID: 34009688 DOI: 10.1111/mec.15989] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2020] [Revised: 04/26/2021] [Accepted: 05/11/2021] [Indexed: 01/01/2023]
Abstract
Characterizing the population history of a species and identifying loci underlying local adaptation is crucial in functional ecology, evolutionary biology, conservation and agronomy. The constant improvement of high-throughput sequencing techniques has facilitated the production of whole genome data in a wide range of species. Population genomics now provides tools to better integrate selection into a historical framework, and take into account selection when reconstructing demographic history. However, this improvement has come with a profusion of analytical tools that can confuse and discourage users. Such confusion limits the amount of information effectively retrieved from complex genomic data sets, and impairs the diffusion of the most recent analytical tools into fields such as conservation biology. It may also lead to redundancy among methods. To address these isssues, we propose an overview of more than 100 state-of-the-art methods that can deal with whole genome data. We summarize the strategies they use to infer demographic history and selection, and discuss some of their limitations. A website listing these methods is available at www.methodspopgen.com.
Collapse
Affiliation(s)
| | - Ben H Warren
- Institut de Systématique, Evolution, Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, CNRS, Sorbonne Université, EPHE, UA, CP 51, Paris, France
| |
Collapse
|
11
|
Sellinger TPP, Abu-Awad D, Tellier A. Limits and convergence properties of the sequentially Markovian coalescent. Mol Ecol Resour 2021; 21:2231-2248. [PMID: 33978324 DOI: 10.1111/1755-0998.13416] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Revised: 04/19/2021] [Accepted: 04/29/2021] [Indexed: 02/07/2023]
Abstract
Several methods based on the sequentially Markovian coalescent (SMC) make use of full genome sequence data from samples to infer population demographic history including past changes in population size, admixture, migration events and population structure. More recently, the original theoretical framework has been extended to allow the simultaneous estimation of population size changes along with other life history traits such as selfing or seed banking. The latter developments enhance the applicability of SMC methods to nonmodel species. Although convergence proofs have been given using simulated data in a few specific cases, an in-depth investigation of the limitations of SMC methods is lacking. In order to explore such limits, we first develop a tool inferring the best case convergence of SMC methods assuming the true underlying coalescent genealogies are known. This tool can be used to quantify the amount and type of information that can be confidently retrieved from given data sets prior to the analysis of the real data. Second, we assess the inference accuracy when the assumptions of SMC approaches are violated due to departures from the model, namely the presence of transposable elements, variable recombination and mutation rates along the sequence, and SNP calling errors. Third, we deliver a new interpretation of SMC methods by highlighting the importance of the transition matrix, which we argue can be used as a set of summary statistics in other statistical inference methods, uncoupling the SMC from hidden Markov models (HMMs). We finally offer recommendations to better apply SMC methods and build adequate data sets under budget constraints.
Collapse
Affiliation(s)
| | - Diala Abu-Awad
- Department of Life Science Systems, Technical University of Munich, Munchen, Germany
| | - Aurélien Tellier
- Department of Life Science Systems, Technical University of Munich, Munchen, Germany
| |
Collapse
|
12
|
Arredondo A, Mourato B, Nguyen K, Boitard S, Rodríguez W, Noûs C, Mazet O, Chikhi L. Inferring number of populations and changes in connectivity under the n-island model. Heredity (Edinb) 2021; 126:896-912. [PMID: 33846579 PMCID: PMC8178352 DOI: 10.1038/s41437-021-00426-9] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2020] [Revised: 03/11/2021] [Accepted: 03/12/2021] [Indexed: 12/11/2022] Open
Abstract
Inferring the demographic history of species is one of the greatest challenges in populations genetics. This history is often represented as a history of size changes, ignoring population structure. Alternatively, when structure is assumed, it is defined a priori as a population tree and not inferred. Here we propose a framework based on the IICR (Inverse Instantaneous Coalescence Rate). The IICR can be estimated for a single diploid individual using the PSMC method of Li and Durbin (2011). For an isolated panmictic population, the IICR matches the population size history, and this is how the PSMC outputs are generally interpreted. However, it is increasingly acknowledged that the IICR is a function of the demographic model and sampling scheme with limited connection to population size changes. Our method fits observed IICR curves of diploid individuals with IICR curves obtained under piecewise stationary symmetrical island models. In our models we assume a fixed number of time periods during which gene flow is constant, but gene flow is allowed to change between time periods. We infer the number of islands, their sizes, the periods at which connectivity changes and the corresponding rates of connectivity. Validation with simulated data showed that the method can accurately recover most of the scenario parameters. Our application to a set of five human PSMCs yielded demographic histories that are in agreement with previous studies using similar methods and with recent research suggesting ancient human structure. They are in contrast with the view of human evolution consisting of one ancestral population branching into three large continental and panmictic populations with varying degrees of connectivity and no population structure within each continent.
Collapse
Affiliation(s)
- Armando Arredondo
- Université de Toulouse, Institut National des Sciences Appliquées, Institut de Mathématiques de Toulouse, Toulouse, France. .,Institut de Mathématiques de Toulouse; UMR5219. Université de Toulouse, Toulouse, France.
| | - Beatriz Mourato
- Institut de Mathématiques de Toulouse; UMR5219. Université de Toulouse, Toulouse, France.,Instituto Gulbenkian de Ciência, Oeiras, Portugal
| | - Khoa Nguyen
- Université de Toulouse, Institut National des Sciences Appliquées, Institut de Mathématiques de Toulouse, Toulouse, France
| | - Simon Boitard
- CBGP, Université de Montpellier, CIRAD, INRAE, Institut Agro, IRD, Montpellier, France
| | - Willy Rodríguez
- Institut de Mathématiques de Toulouse; UMR5219. Université de Toulouse, Toulouse, France.,ENAC - Ecole Nationale de l'Aviation Civile, Université de Toulouse, Toulouse, France
| | | | - Olivier Mazet
- Université de Toulouse, Institut National des Sciences Appliquées, Institut de Mathématiques de Toulouse, Toulouse, France.,Institut de Mathématiques de Toulouse; UMR5219. Université de Toulouse, Toulouse, France
| | - Lounès Chikhi
- Instituto Gulbenkian de Ciência, Oeiras, Portugal. .,Laboratoire Évolution & Diversité Biologique (EDB UMR 5174), CNRS, IRD, UPS, Université de Toulouse Midi-Pyrénées, Toulouse, France.
| |
Collapse
|
13
|
Kerdoncuff E, Lambert A, Achaz G. Testing for population decline using maximal linkage disequilibrium blocks. Theor Popul Biol 2020; 134:171-181. [DOI: 10.1016/j.tpb.2020.03.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2019] [Revised: 03/26/2020] [Accepted: 03/29/2020] [Indexed: 02/02/2023]
|
14
|
Feutry P, Devloo-Delva F, Tran Lu Y A, Mona S, Gunasekera RM, Johnson G, Pillans RD, Jaccoud D, Kilian A, Morgan DL, Saunders T, Bax NJ, Kyne PM. One panel to rule them all: DArTcap genotyping for population structure, historical demography, and kinship analyses, and its application to a threatened shark. Mol Ecol Resour 2020; 20:1470-1485. [PMID: 32492756 DOI: 10.1111/1755-0998.13204] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Revised: 05/26/2020] [Accepted: 05/28/2020] [Indexed: 01/25/2023]
Abstract
With recent advances in sequencing technology, genomic data are changing how important conservation management decisions are made. Applications such as Close-Kin Mark-Recapture demand large amounts of data to estimate population size and structure, and their full potential can only be realised through ongoing improvements in genotyping strategies. Here we introduce DArTcap, a cost-efficient method that combines DArTseq and sequence capture, and illustrate its use in a high resolution population analysis of Glyphis garricki, a rare, poorly known and threatened euryhaline shark. Clustering analyses and spatial distribution of kin pairs from four different regions across northern Australia and one in Papua New Guinea, representing its entire known range, revealed that each region hosts at least one distinct population. Further structuring is likely within Van Diemen Gulf, the region that included the most rivers sampled, suggesting additional population structuring would be found if other rivers were sampled. Coalescent analyses and spatially explicit modelling suggest that G. garricki experienced a recent range expansion during the opening of the Gulf of Carpentaria following the conclusion of the Last Glacial Maximum. The low migration rates between neighbouring populations of a species that is found only in restricted coastal and riverine habitats show the importance of managing each population separately, including careful monitoring of local and remote anthropogenic activities that may affect their environments. Overall we demonstrated how a carefully chosen SNP panel combined with DArTcap can provide highly accurate kinship inference and also support population structure and historical demography analyses, therefore maximising cost-effectiveness.
Collapse
Affiliation(s)
| | - Floriaan Devloo-Delva
- CSIRO Oceans and Atmosphere, Hobart, TAS, Australia.,School of Natural Sciences - Quantitative Marine Science, University of Tasmania, Hobart, TAS, Australia
| | - Adrien Tran Lu Y
- Institut de Systématique, Évolution, Biodiversité (ISYEB), UMR 7205, MNHN, CNRS, EPHE, Sorbonne Université, Paris, France.,EPHE, PSL University, Paris, France
| | - Stefano Mona
- Institut de Systématique, Évolution, Biodiversité (ISYEB), UMR 7205, MNHN, CNRS, EPHE, Sorbonne Université, Paris, France.,EPHE, PSL University, Paris, France.,Laboratoire d'Excellence CORAIL, Papetoai, French Polynesia
| | | | - Grant Johnson
- Department of Primary Industry and Fisheries, Aquatic Resource Research Unit, Darwin, NT, Australia
| | | | - Damian Jaccoud
- Diversity Arrays Technology Pty Ltd, University of Canberra, Bruce, ACT, Australia
| | - Andrzej Kilian
- Diversity Arrays Technology Pty Ltd, University of Canberra, Bruce, ACT, Australia
| | - David L Morgan
- Centre for Sustainable Aquatic Ecosystems, Harry Butler Institute, Murdoch University, Murdoch, WA, Australia
| | - Thor Saunders
- Department of Primary Industry and Fisheries, Aquatic Resource Research Unit, Darwin, NT, Australia
| | - Nicholas J Bax
- CSIRO Oceans and Atmosphere, Hobart, TAS, Australia.,Institute for Marine and Antarctic Science, University of Tasmania, Hobart, TAS, Australia
| | - Peter M Kyne
- Research Institute for the Environment and Livelihoods, Charles Darwin University, Darwin, NT, Australia
| |
Collapse
|
15
|
Poelstra JW, Salmona J, Tiley GP, Schüßler D, Blanco MB, Andriambeloson JB, Bouchez O, Campbell CR, Etter PD, Hohenlohe PA, Hunnicutt KE, Iribar A, Johnson EA, Kappeler PM, Larsen PA, Manzi S, Ralison JM, Randrianambinina B, Rasoloarison RM, Rasolofoson DW, Stahlke AR, Weisrock DW, Williams RC, Chikhi L, Louis EE, Radespiel U, Yoder AD. Cryptic Patterns of Speciation in Cryptic Primates: Microendemic Mouse Lemurs and the Multispecies Coalescent. Syst Biol 2020; 70:203-218. [PMID: 32642760 DOI: 10.1093/sysbio/syaa053] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2019] [Revised: 06/13/2020] [Accepted: 06/23/2020] [Indexed: 12/21/2022] Open
Abstract
Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic revision has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we investigate mouse lemur diversity in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for two pairs of sister lineages that include three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among the two pairs of lineages. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species, M. mittermeieri and M. lehilahytsara is not supported-a result that is particularly striking when using the genealogical discordance index (gdi). Nonsister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated relaxed-clock estimates. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case. [Cryptic species; effective population size; microendemism; multispecies coalescent; speciation; species delimitation.].
Collapse
Affiliation(s)
| | - Jordi Salmona
- CNRS, Université Paul Sabatier, IRD; UMR5174 EDB (Laboratoire Évolution & Diversité Biologique), 118 route de Narbonne, 31062 Toulouse, France
| | - George P Tiley
- Department of Biology, Duke University, Durham, NC 27708, USA
| | - Dominik Schüßler
- Research Group Ecology and Environmental Education, Department of Biology, University of Hildesheim, Universitaetsplatz 1, 31141 Hildesheim, Germany
| | - Marina B Blanco
- Department of Biology, Duke University, Durham, NC 27708, USA.,Duke Lemur Center, Duke University, Durham, NC 27705, USA
| | - Jean B Andriambeloson
- Department of Zoology and Animal Biodiversity, University of Antananarivo, Antananarivo 101, Madagascar
| | - Olivier Bouchez
- INRA, US 1426, GeT-PlaGe, Genotoul, Castanet-Tolosan, France
| | - C Ryan Campbell
- Department of Biology, Duke University, Durham, NC 27708, USA.,Department of Evolutionary Anthropology, Duke University, Durham, NC 27708, USA
| | - Paul D Etter
- Institute of Molecular Biology, University of Oregon, Eugene, OR, USA
| | - Paul A Hohenlohe
- Department of Biological Sciences, Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID 83844, USA
| | - Kelsie E Hunnicutt
- Department of Biology, Duke University, Durham, NC 27708, USA.,Department of Biological Sciences, University of Denver, Denver, CO 80208, USA
| | - Amaia Iribar
- CNRS, Université Paul Sabatier, IRD; UMR5174 EDB (Laboratoire Évolution & Diversité Biologique), 118 route de Narbonne, 31062 Toulouse, France
| | - Eric A Johnson
- Institute of Molecular Biology, University of Oregon, Eugene, OR, USA
| | - Peter M Kappeler
- Behavioral Ecology and Sociobiology Unit, German Primate Center, Kellnerweg 6, 37077 Göttingen, Germany
| | - Peter A Larsen
- Department of Biology, Duke University, Durham, NC 27708, USA.,Department of Veterinary and Biomedical Sciences, University of Minnesota, Saint Paul, MN 55108, USA
| | - Sophie Manzi
- CNRS, Université Paul Sabatier, IRD; UMR5174 EDB (Laboratoire Évolution & Diversité Biologique), 118 route de Narbonne, 31062 Toulouse, France
| | - JosÉ M Ralison
- Department of Zoology and Animal Biodiversity, University of Antananarivo, Antananarivo 101, Madagascar
| | - Blanchard Randrianambinina
- Groupe d'Etude et de Recherche sur les Primates de Madagascar (GERP), BP 779, Antananarivo 101, Madagascar.,Faculté des Sciences, University of Mahajanga, Mahajanga, Madagascar
| | - Rodin M Rasoloarison
- Behavioral Ecology and Sociobiology Unit, German Primate Center, Kellnerweg 6, 37077 Göttingen, Germany
| | - David W Rasolofoson
- Groupe d'Etude et de Recherche sur les Primates de Madagascar (GERP), BP 779, Antananarivo 101, Madagascar
| | - Amanda R Stahlke
- Department of Biological Sciences, Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID 83844, USA
| | - David W Weisrock
- Department of Biology, University of Kentucky, Lexington, KY, 40506, USA
| | - Rachel C Williams
- Department of Biology, Duke University, Durham, NC 27708, USA.,Duke Lemur Center, Duke University, Durham, NC 27705, USA
| | - LounÈs Chikhi
- CNRS, Université Paul Sabatier, IRD; UMR5174 EDB (Laboratoire Évolution & Diversité Biologique), 118 route de Narbonne, 31062 Toulouse, France.,Instituto Gulbenkian de Ciência, Rua da Quinta Grande, 6, 2780-156 Oeiras, Portugal
| | - Edward E Louis
- Grewcock Center for Conservation and Research, Omaha's Henry Doorly Zoo and Aquarium, Omaha, NE, USA
| | - Ute Radespiel
- Institute of Zoology, University of Veterinary Medicine Hannover, Buenteweg 17, 30559 Hannover, Germany Jelmer Poelstra, Jordi Salmona, George P. Tiley are the joint first authors. Ute Radespiel and Anne D. Yoder are the joint senior authors
| | - Anne D Yoder
- Department of Biology, Duke University, Durham, NC 27708, USA
| |
Collapse
|
16
|
Sellinger TPP, Abu Awad D, Moest M, Tellier A. Inference of past demography, dormancy and self-fertilization rates from whole genome sequence data. PLoS Genet 2020; 16:e1008698. [PMID: 32251472 PMCID: PMC7173940 DOI: 10.1371/journal.pgen.1008698] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2019] [Revised: 04/21/2020] [Accepted: 02/24/2020] [Indexed: 02/04/2023] Open
Abstract
Several methods based on the Sequential Markovian coalescence (SMC) have been developed that make use of genome sequence data to uncover population demographic history, which is of interest in its own right and is a key requirement to generate a null model for selection tests. While these methods can be applied to all possible kind of species, the underlying assumptions are sexual reproduction in each generation and non-overlapping generations. However, in many plants, invertebrates, fungi and other taxa, those assumptions are often violated due to different ecological and life history traits, such as self-fertilization or long term dormant structures (seed or egg-banking). We develop a novel SMC-based method to infer 1) the rates/parameters of dormancy and of self-fertilization, and 2) the populations' past demographic history. Using simulated data sets, we demonstrate the accuracy of our method for a wide range of demographic scenarios and for sequence lengths from one to 30 Mb using four sampled genomes. Finally, we apply our method to a Swedish and a German population of Arabidopsis thaliana demonstrating a selfing rate of ca. 0.87 and the absence of any detectable seed-bank. In contrast, we show that the water flea Daphnia pulex exhibits a long lived egg-bank of three to 18 generations. In conclusion, we here present a novel method to infer accurate demographies and life-history traits for species with selfing and/or seed/egg-banks. Finally, we provide recommendations for the use of SMC-based methods for non-model organisms, highlighting the importance of the per site and the effective ratios of recombination over mutation.
Collapse
Affiliation(s)
| | - Diala Abu Awad
- Department of Population Genetics, Technische Universitaet Muenchen, Freising, Germany
| | - Markus Moest
- Department of Ecology, University of Innsbruck, Innsbruck, Austria
| | - Aurélien Tellier
- Department of Population Genetics, Technische Universitaet Muenchen, Freising, Germany
| |
Collapse
|
17
|
Khimoun A, Doums C, Molet M, Kaufmann B, Peronnet R, Eyer PA, Mona S. Urbanization without isolation: the absence of genetic structure among cities and forests in the tiny acorn ant Temnothorax nylanderi. Biol Lett 2020; 16:20190741. [PMID: 31992150 DOI: 10.1098/rsbl.2019.0741] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Urban alteration of neutral and adaptive evolutionary processes is still underexplored. Using a genome-wide SNP dataset, we investigated (i) urban-induced modifications of population demography, genetic diversity and population structure and (ii) signature of divergent selection between urban and forest populations in the ant species, Temnothorax nylanderi. Our results did not reveal an impact of urbanization on neutral processes since we observed: (i) analogous genetic diversity among paired urban/forest sites and two control populations; (ii) weak population genetic structure explained neither by habitat (urban versus forest) nor by geography; (iii) a remarkably similar demographic history across populations with an ancestral growth followed by a recent decline, regardless of their current habitat or geographical location. The micro-geographical home range of ants may explain their resilience to urbanization. Finally, we detected 19 candidate loci discriminating urban/forest populations and associated with core cellular components, molecular function or biological process. Two of these loci were associated with a gene ontology term that was previously found to belong to a module of co-expressed genes related to caste phenotype. These results call for transcriptomics analyses to identify genes associated with ant social traits and to infer their potential role in urban adaptation.
Collapse
Affiliation(s)
- A Khimoun
- Biogéosciences, UMR 6282 CNRS, Université Bourgogne Franche-Comté, 6 Boulevard Gabriel, 21000 Dijon, France
| | - C Doums
- Institut de Systématique, Évolution, Biodiversité (ISYEB), UMR 7205, MNHN, CNRS, EPHE, Sorbonne Université, 75005 Paris, France.,EPHE, PSL University, 75005 Paris, France
| | - M Molet
- Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), UMR 7618, Sorbonne Université, Université Paris Est Créteil, Université Paris Diderot, CNRS, INRA, IRD, 75005 Paris, France
| | - B Kaufmann
- Université de Lyon, UMR5023 Ecologie des Hydrosystèmes Naturels et Anthropisés, Université Lyon 1, ENTPE, CNRS, Villeurbanne 69622, France
| | - R Peronnet
- Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), UMR 7618, Sorbonne Université, Université Paris Est Créteil, Université Paris Diderot, CNRS, INRA, IRD, 75005 Paris, France
| | - P A Eyer
- Department of Entomology, Texas A&M University, 2143 TAMU, College Station, TX 77843-2143, USA
| | - S Mona
- Institut de Systématique, Évolution, Biodiversité (ISYEB), UMR 7205, MNHN, CNRS, EPHE, Sorbonne Université, 75005 Paris, France.,EPHE, PSL University, 75005 Paris, France
| |
Collapse
|
18
|
Blath J, Buzzoni E, Koskela J, Wilke Berenguer M. Statistical tools for seed bank detection. Theor Popul Biol 2020; 132:1-15. [PMID: 31945384 DOI: 10.1016/j.tpb.2020.01.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2019] [Revised: 12/24/2019] [Accepted: 01/02/2020] [Indexed: 10/25/2022]
Abstract
We derive statistical tools to analyze the patterns of genetic variability produced by models related to seed banks; in particular the Kingman coalescent, its time-changed counterpart describing so-called weak seed banks, the strong seed bank coalescent, and the two-island structured coalescent. As (strong) seed banks stratify a population, we expect them to produce a signal comparable to population structure. We present tractable formulas for Wright's FST and the expected site frequency spectrum for these models, and show that they can distinguish between some models for certain ranges of parameters. We then use pseudo-marginal MCMC to show that the full likelihood can reliably distinguish between all models in the presence of parameter uncertainty under moderate stratification, and point out statistical pitfalls arising from stratification that is either too strong or too weak. We further show that it is possible to infer parameters, and in particular determine whether mutation is taking place in the (strong) seed bank.
Collapse
Affiliation(s)
- Jochen Blath
- Institut für Mathematik, Technische Universität Berlin, Straße des 17. Juni 136, 10623 Berlin, Germany.
| | - Eugenio Buzzoni
- Institut für Mathematik, Technische Universität Berlin, Straße des 17. Juni 136, 10623 Berlin, Germany.
| | - Jere Koskela
- Department of Statistics, University of Warwick, Coventry CV4 7AL, UK.
| | - Maite Wilke Berenguer
- Fakultät für Mathematik, Ruhr-Universität Bochum, Universitätsstraße 150, 44801 Bochum, Germany.
| |
Collapse
|
19
|
|