Estimating dispersal rates and locating genetic ancestors with genome-wide genealogies
- University of Toronto, Canada
- University of California, Davis, United States
Abstract
Spatial patterns in genetic diversity are shaped by individuals dispersing from their parents and larger-scale population movements. It has long been appreciated that these patterns of movement shape the underlying genealogies along the genome leading to geographic patterns of isolation by distance in contemporary population genetic data. However, extracting the enormous amount of information contained in genealogies along recombining sequences has, until recently, not been computationally feasible. Here we capitalize on important recent advances in genome-wide gene-genealogy reconstruction and develop methods to use thousands of trees to estimate per-generation dispersal rates and to locate the genetic ancestors of a sample back through time. We take a likelihood approach in continuous space using a simple approximate model (branching Brownian motion) as our prior distribution of spatial genealogies. After testing our method with simulations we apply it to Arabidopsis thaliana. We estimate a dispersal rate of roughly 60km2 per generation, slightly higher across latitude than across longitude, potentially reflecting a northward post-glacial expansion. Locating ancestors allows us to visualize major geographic movements, alternative geographic histories, and admixture. Our method highlights the huge amount of information about past dispersal events and population movements contained in genome-wide genealogies.
Data availability
All code used to perform the analyses in this study can be found at \url{https://github.com/mmosmond/spacetrees-ms}. More information on how to run our method, \texttt{spacetrees}, is available at \url{https://github.com/osmond-lab/spacetrees}.
-
Relate-inferred genealogies for 66 longread Arabidopsis thaliana genomesZenodo, doi:10.5281/zenodo.5099656.
-
Towards telomere-to-telomere assemblies of Arabidopsis thaliana genomes: CLR vs HiFi readsNational Center for Biotechnology Information.
-
ddAraThal4 (thale cress)National Center for Biotechnology Information.
Article and author information
Author details
Funding
Natural Sciences and Engineering Research Council of Canada (RGPIN-2021-03207)
- Matthew Osmond
Natural Sciences and Engineering Research Council of Canada (DGECR-2021-00114)
- Matthew Osmond
Banting Research Foundation
- Matthew Osmond
National Institute of General Medical Sciences (R01 GM108779)
- Graham Coop
National Institute of General Medical Sciences (R35 GM136290)
- Graham Coop
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2024, Osmond & Coop
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 668
- views
-
- 101
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Estimating dispersal rates and locating genetic ancestors with genome-wide genealogies
eLife 13:e72177.
https://doi.org/10.7554/eLife.72177
Further reading
-
- Evolutionary Biology
- Neuroscience
The first complete 3D reconstruction of the compound eye of a minute wasp species sheds light on the nuts and bolts of size reduction.
-
- Evolutionary Biology
- Genetics and Genomics
Copy number variants (CNVs) are an important source of genetic variation underlying rapid adaptation and genome evolution. Whereas point mutation rates vary with genomic location and local DNA features, the role of genome architecture in the formation and evolutionary dynamics of CNVs is poorly understood. Previously, we found the GAP1 gene in Saccharomyces cerevisiae undergoes frequent amplification and selection in glutamine-limitation. The gene is flanked by two long terminal repeats (LTRs) and proximate to an origin of DNA replication (autonomously replicating sequence, ARS), which likely promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we evolved engineered strains lacking either the adjacent LTRs, ARS, or all elements in glutamine-limited chemostats. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. Removal of local DNA elements significantly impacts the fitness effect of GAP1 CNVs and the rate of adaptation. In 177 CNV lineages, across all four strains, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, distal ones mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination can mediate gene amplification following de novo retrotransposon events. Our study reveals that template switching during DNA replication is a prevalent source of adaptive CNVs.
