Jointly representing long-range genetic similarity and spatially heterogeneous isolation-by-distance
- 1. University of Chicago
- 2. University of Bologna
Description
Isolation-by-distance patterns in genetic variation are a widespread feature of the geographic structure of genetic variation in many species, and many methods have been developed to illuminate such patterns in genetic data. However, long-range genetic similarities also exist, often as a result of rare or episodic long-range gene flow. Jointly characterizing patterns of isolation-by-distance and long-range genetic similarity in genetic data is an open data analysis challenge that, if resolved, could help produce more complete representations of the geographic structure of genetic data in any given species. Here, we present a computationally tractable method that identifies long-range genetic similarities in a background of spatially heterogeneous isolation-by-distance variation. The method uses a coalescent-based framework, and models long-range genetic similarity in terms of directional events with source fractions describing the fraction of ancestry at a location tracing back to a remote source. The method produces geographic maps annotated with inferred long-range edges, as well as maps of uncertainty in the geographic location of each source of long-range gene flow. We have implemented the method in a package called FEEMSmix (an extension to FEEMS), and validated its implementation using simulations representative of typical data applications. We also apply this method to two empirical data sets. In a data set of over 4,000 humans (Homo sapiens) across Afro-Eurasia, we recover many known signals of long-distance dispersal from recent centuries. Similarly, in a data set of over 100 gray wolves (Canis lupus) across North America, we identify several previously unknown long-range connections, some of which were attributable to recording errors in sampling locations. Therefore, beyond identifying genuine long-range dispersals, our approach also serves as a useful tool for quality control in spatial genetic studies.
Data availability
The wolves data set is provided as part of the FEEMS package in (https://doi.org/10.7554/eLife.61927) (and is also publicly available from the original publication, https://doi.org/10.1111/mec.13364). This data set can be found at https://doi.org/10.5061/dryad.c9b25. The corrected wolves data set and the human data set used in this study can be found at https://doi.org/10.5061/dryad.p8cz8wb18 and https://zenodo.org/records/15007585. All simulated data can be reproduced using code in https://github.com/VivaswatS/feems/tree/admixture_edge. Finally, FEEMSmix is readily available as a complete python package from https://github.com/NovembreLab/feems.
Additional details
Identifiers
- DOI
- 10.1371/journal.pgen.1011612
- Other
- oai:uchicago.tind.io:16377
Funding
- National Institute of General Medical Sciences
- R35 GM149521
- National Institute of General Medical Sciences
- R01 GM132383
- NextGenerationEU, National Recovery and Resilience Plan
- National Biodiversity Future Center -NBFC