A Rarefaction Approach to Identify Local Introgression in a Three Population Tree
- 1The Pennsylvania State University
Abstract
The D statistic, also known as the $ABBA-BABA$ statistic, is widely used to detect the presence of archaic genome-wide introgression between two non-sister taxa. D counts the imbalance between the number of biallelic sites where either the second and third taxa (ABBA site) share the derived allele or the first and third taxa (BABA site) share the derived allele in a four taxa tree. Here, the fourth taxon acts as an outgroup to determine the ancestral allele. When there is no introgression, these counts are expected to be equal, and a discordance between counts suggests introgression from the third taxon into either the first or second. D is limited to the detection of genome-wide introgression and exhibits a high false-positive rate when applied to smaller genomic segments. Here, we present a new method, D STatistic with Allelic Rarefaction (dstar), to address these limitations. dstar uses multiple lineages and does not require an outgroup to calculate the imbalance between the number of alleles found exclusively in the second and third taxa and the number of alleles found exclusively in the first and third taxa. dstar employs a rarefaction technique to correct for unequal sample-size and allows multiallelic sites. We use simulations to show that dstar has better precision and recall for detecting introgressed segments of DNA when compared to other methods. We conclude by recovering Denisovan DNA related to immune function in modern day Papuans. Precompiled executables, the manual, source code, and simulation and analysis scripts used in this study can be found at \url{https://github.com/TQ-Smith/DSTAR}
Related
- Long-distance dispersal suppresses introgression of local alleles during range expansions
- Adaptive introgression in subdivided populations leaves distinct genomic footprints under global and local selection: a case study in Ciona intestinalis
- Precision-weighted mixed models identify high-yielding and temporally stable cowpea germplasm across three growing seasons
- Label-free flow operators and stability operators from single-cell time courses: what local dynamical operators can and cannot identify
- Machine learning-assisted genomic characterization of population differentiation between Chinese Kirghiz and Kyrgyzstan local sheep
BibTeX
@misc{smith2026rarefaction,
title = {{A Rarefaction Approach to Identify Local Introgression in a Three Population Tree}},
author = {T. Q. Smith and Z. A. Szpiech},
howpublished = {bioRxiv},
year = {2026},
month = may,
doi = {10.64898/2026.05.13.724952},
url = {https://www.biorxiv.org/content/10.64898/2026.05.13.724952v2},
note = {CC BY},
}
RIS
TY - UNPB
TI - A Rarefaction Approach to Identify Local Introgression in a Three Population Tree
AU - T. Q. Smith
AU - Z. A. Szpiech
PY - 2026
DA - 2026/05/13
PB - bioRxiv
DO - 10.64898/2026.05.13.724952
UR - https://www.biorxiv.org/content/10.64898/2026.05.13.724952v2
AB - The D statistic, also known as the $ABBA-BABA$ statistic, is widely used to detect the presence of archaic genome-wide introgression between two non-sister taxa. D counts the imbalance between the number of biallelic sites where either the second and third taxa (ABBA site) share the derived allele or the first and third taxa (BABA site) share the derived allele in a four taxa tree. Here, the fourth taxon acts as an outgroup to determine the ancestral allele. When there is no introgression, these counts are expected to be equal, and a discordance between counts suggests introgression from the third taxon into either the first or second. D is limited to the detection of genome-wide introgression and exhibits a high false-positive rate when applied to smaller genomic segments. Here, we present a new method, D STatistic with Allelic Rarefaction (dstar), to address these limitations. dstar uses multiple lineages and does not require an outgroup to calculate the imbalance between the number of alleles found exclusively in the second and third taxa and the number of alleles found exclusively in the first and third taxa. dstar employs a rarefaction technique to correct for unequal sample-size and allows multiallelic sites. We use simulations to show that dstar has better precision and recall for detecting introgressed segments of DNA when compared to other methods. We conclude by recovering Denisovan DNA related to immune function in modern day Papuans. Precompiled executables, the manual, source code, and simulation and analysis scripts used in this study can be found at \url{https://github.com/TQ-Smith/DSTAR}
ER -