No retraction, correction or concern found when last checked on Oct 29, 2026.

A Hidden Markov Model Approach for Simultaneously Estimating Local Ancestry and Admixture Time Using Next Generation Sequence Data in Samples of Arbitrary Ploidy

Russell Corbett-Detig, Rasmus Nielsen

PLOS Genetics · 2017 · CC BY

doi:10.1371/journal.pgen.1006529

Abstract

Admixture—the mixing of genomes from divergent populations—is increasingly appreciated as a central process in evolution. To characterize and quantify patterns of admixture across the genome, a number of methods have been developed for local ancestry inference. However, existing approaches have a number of shortcomings. First, all local ancestry inference methods require some prior assumption about the expected ancestry tract lengths. Second, existing methods generally require genotypes, which is not feasible to obtain for many next-generation sequencing projects. Third, many methods assume samples are diploid, however a wide variety of sequencing applications will fail to meet this assumption. To address these issues, we introduce a novel hidden Markov model for estimating local ancestry that models the read pileup data, rather than genotypes, is generalized to arbitrary ploidy, and can estimate the time since admixture during local ancestry inference. We demonstrate that our method can simultaneously estimate the time since admixture and local ancestry with good accuracy, and that it performs well on samples of high ploidy— i.e. 100 or more chromosomes. As this method is very general, we expect it will be useful for local ancestry inference in a wider variety of populations than what previously has been possible. We then applied our method to pooled sequencing data derived from populations of Drosophila melanogaster on an ancestry cline on the east coast of North America. We find that regions of local recombination rates are negatively correlated with the proportion of African ancestry, suggesting that selection against foreign ancestry is the least efficient in low recombination regions. Finally we show that clinal outlier loci are enriched for genes associated with gene regulatory functions, consistent with a role of regulatory evolution in ecological adaptation of admixed D. melanogaster populations. Our results illustrate the potential of local ancestry inference for elucidating fundamental evolutionary processes. Author Summary When divergent populations hybridize, their offspring obtain portions of their genomes from each parent population. Although the average ancestry proportion in each descendant is equal to the proportion of ancestors from each of the ancestral populations, the contribution of each ancestry type is variable across the genome. Estimating local ancestry within admixed individuals is a fundamental goal for evolutionary genetics, and here we develop a method for doing this that circumvents many of the problems associated with existing methods. Briefly, our method can use short read data, rather than genotypes and can be applied to samples with any number of chromosomes. Furthermore, our method simultaneously estimates local ancestry and the number of generations since admixture—the time that the two ancestral populations first encountered each other. Finally, in applying our method to data from an admixture zone between ancestral populations of Drosophila melanogaster , we find many lines of evidence consistent with natural selection operating to against the introduction of foreign ancestry into populations of one predominant ancestry type. Because of the generality of this method, we expect that it will be useful for a wide variety of existing and ongoing research projects.

References 97

1.
Kronforst MR, Young LG, Blume LM, Gilbert LE. Multilocus analyses of admixture and introgression among hybridizing Heliconius butterflies. Evolution. 2006;60: 1254–16.
Not in our library
2.
Population Genomics of Sub-Saharan Drosophila melanogaster: African Diversity and Non-African Admixture · 2012
No notice
3.
The Genomic Signature of Crop-Wild Introgression in Maize · 2013
Corrected
4.
Speciation and Introgression between Mimulus nasutus and Mimulus guttatus · 2014
No notice
5.
Rieseberg LH, Raymond O, Rosenthal DM, Lai Z, Livingstone K, Nakazato T, et al Major ecological transitions in wild sunflowers facilitated by hybridization. Science. 2003;301: 1211–1216.
Not in our library10.1126/science.1086949
6.
Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, et al A draft sequence of the neandertal genome. Science. 2010;328: 710–722.
Not in our library10.1126/science.1188021
7.
Sankararaman S, Mallick S, Dannemann M, Prüfer K, Kelso J, Pääbo S, et al The genomic landscape of Neanderthal ancestry in present-day humans. Nature. 2014;507: 354–357.
Not in our library10.1038/nature12961
8.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38: 904–909.
Not in our library10.1038/ng1847
9.
Lohmueller KE, Pearce CL, Pike M, Lander ES, Hirschhorn JN. Meta-analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat Genet. 2003;33: 177–182.
Not in our library10.1038/ng1071
10.
Freedman ML, Reich D, Penney KL, McDonald GJ, Mignault AA, Patterson N, et al Assessing the impact of population stratification on genetic association studies. Nat Genet. 2004;36: 388–393.
Not in our library10.1038/ng1333
11.
Caracristi G, Schlotterer C. Genetic differentiation between American and European Drosophila melanogaster populations could be attributed to admixture of African alleles. Mol Biol Evol. 2003;20: 792–799.
Not in our library10.1093/molbev/msg091
12.
Kolbe JJ, Glor RE, Schettino LR, Lara AC, Larson A. Genetic variation increases during biological invasion by a Cuban lizard. Nature. 2004;431: 177–181.
Not in our library10.1038/nature02807
13.
Consortium THG, Consortium G. Butterfly genome reveals promiscuous exchange of mimicry adaptations among species. Nature. 2012;487: 94–98.
Not in our library10.1038/nature11041
14.
Racimo F, Sankararaman S, Nielsen R, Huerta-Sanchez E. Evidence for archaic adaptive introgression in humans. Nat Genet. 2015;16: 359–371.
Not in our library10.1038/nrg3936
15.
The Strength of Selection against Neanderthal Introgression · 2016
No notice
16.
Harris K, Nielsen R. The genetic cost of neanderthal introgression. Genetics. 2016;203: 881–891.
Not in our library10.1534/genetics.116.186890
17.
Pool JE. The mosaic ancestry of the Drosophila Genetic Reference Panel and the D. melanogaster reference genome reveals a network of epistatic fitness interactions. Mol Biol Evol. 2015;: msv194–16.
Not in our library10.1093/molbev/msv194
18.
High-resolution mapping reveals hundreds of genetic incompatibilities in hybridizing fish species · 2014
No notice
19.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155: 945–959.
Not in our library10.1093/genetics/155.2.945
20.
Skotte L, Korneliussen TS, Albrechtsen A. Estimating individual admixture proportions from next generation sequencing data. Genetics. 2013;195: 693–702.
Not in our library10.1534/genetics.113.154138
21.
Bergland AO, Tobler R, González J, Schmidt P, Petrov D. Secondary contact and local adaptation contribute to genome-wide patterns of clinal variation in Drosophila melanogaster. Mol Ecol. 2016;25: 1157–1174.
Not in our library10.1111/mec.13455
22.
Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003;164: 1567–1587.
Not in our library10.1093/genetics/164.4.1567
23.
Sankararaman S, Sridhar S, Kimmel G, Halperin E. Estimating local ancestry in admixed populations. Am J Hum Genet. 2008;82: 290–303.
Not in our library10.1016/j.ajhg.2007.09.022
24.
Maples BK, Gravel S, Kenny EE, Bustamante CD. RFMix: a discriminative modeling approach for rapid and robust local-ancestry inference. Am J Hum Genet. 2013;93: 278–288.
Not in our library10.1016/j.ajhg.2013.06.020
25.
Baran Y, Pasaniuc B, Sankararaman S, Torgerson DG, Gignoux C, Eng C, et al Fast and accurate inference of local ancestry in Latino populations. Bioinformatics. 2012;28: 1359–1367.
Not in our library10.1093/bioinformatics/bts144
26.
Pool JE, Nielsen R. Inference of historical changes in migration rate from the lengths of migrant tracts. Genetics. 2009;181: 711–719.
Not in our library10.1534/genetics.108.098095
27.
Koopman W, Li Y, Coart E. Linked vs. unlinked markers: multilocus microsatellite haplotype‐sharing as a tool to estimate gene flow and introgression. Mol Ecol. 2007;16: 243–256.
Not in our library10.1111/j.1365-294x.2006.03137.x
28.
Patterson N, Hattangadi N, Lane B. Methods for high-density admixture mapping of disease genes. Am J Hum Genet. 2004;74: 979–1000.
Not in our library10.1086/420871
29.
Gravel S. Population Genetics Models of Local Ancestry. Genetics. 2012;191: 607–619.
Not in our library10.1534/genetics.112.139808
30.
The History of African Gene Flow into Southern Europeans, Levantines, and Jews · 2011
No notice
31.
Loh PR, Lipson M, Patterson N, Moorjani P. Inferring admixture histories of human populations using linkage disequilibrium. Genetics. 2013;193: 1233–1254.
Not in our library10.1534/genetics.112.147330
32.
Hellenthal G, Busby GBJ, Band G, Wilson JF, Capelli C, Falush D, et al A Genetic Atlas of Human Admixture History. Science. 2014;343: 747–751.
Not in our library10.1126/science.1243518
33.
Consortium 1GP, BGI-Shenzhen, European Bioinformatics Institute, Illumina, Hospital BAW, College B, et al An integrated map of genetic variation from 1,092 human genomes. Nature. 2012.
Not in our library10.1038/nature11632
34.
Pasaniuc B, Sankararaman S, Kimmel G, Halperin E. Inference of locus-specific ancestry in closely related populations. Bioinformatics. 2009;25: i213–i221.
Not in our library10.1093/bioinformatics/btp197
35.
Haddrill PR, Thornton KR, Charlesworth B, Andolfatto P. Multilocus patterns of nucleotide variability and the demographic and selection history of Drosophila melanogaster populations. Genome Res. 2005;15: 790–799.
Not in our library10.1101/gr.3541005
36.
Thornton K, Andolfatto P. Approximate Bayesian inference reveals evidence for a recent, severe bottleneck in a Netherlands population of Drosophila melanogaster. Genetics. 2006;172: 1607–1619.
Not in our library10.1534/genetics.105.048223
37.
Inferring the Demographic History and Rate of Adaptive Substitution in Drosophila · 2006
No notice
38.
Duchen P, Živković D, Hutter S, Stephan W. Demographic inference reveals African and European admixture in the North American Drosophila melanogaster population. Genetics. 2013;193: 291–301.
Not in our library10.1534/genetics.112.145912
39.
Langley CH, Stevens K, Cardeno C, Lee YCG, Schrider DR, Pool JE, et al Genomic variation in natural populations of Drosophila melanogaster. Genetics. 2012;192: 533–598.
Not in our library10.1534/genetics.112.142018
40.
Genomic Evidence of Rapid and Stable Adaptive Oscillations over Seasonal Time Scales in Drosophila · 2014
No notice
41.
Schrider DR, Hahn MW, Begun DJ. Parallel evolution of copy-number variation across continents in Drosophila melanogaster. Mol Biol Evol. 2016;33: 1308–1316.
Not in our library10.1093/molbev/msw014
42.
Reinhardt JA, Kolaczkowski B, Jones CD, Begun DJ, Kern AD. Parallel geographic variation in Drosophila melanogaster. Genetics. 2014;197: 361–373.
Not in our library10.1534/genetics.114.161463
43.
Genome‐wide patterns of latitudinal differentiation among populations of D rosophila melanogaster from N orth A merica · 2012
No notice
44.
Oakeshott JG, Gibson JB, Anderson PR, Knibb WR. Alcohol dehydrogenase and glycerol-3-phosphate dehydrogenase clines in Drosophila melanogaster on different continents. Evolution. 1982;36: 86–96.
Not in our library10.1111/j.1558-5646.1982.tb05013.x
45.
Berry A, Kreitman M. Molecular analysis of an allozyme cline: alcohol dehydrogenase in Drosophila melanogaster on the east coast of North America. Genetics. 1993;134: 869–893.
Not in our library10.1093/genetics/134.3.869
46.
Yukilevich R, True JR. African morphology, behavior, and phermones underlie incipient sexual isolation between US and Caribbean Drosophila melanogaster. Evolution. 2008;62: 2807–2828.
Not in our library10.1111/j.1558-5646.2008.00488.x
47.
Lachance J, True JR. X-autosome incompatibilities in Drosophila melanogaster: test of Haldane's rue and geographic patterns within species. Evolution. 2010;64: 3035–3046.
Not in our library10.1111/j.1558-5646.2010.01028.x
48.
Corbett-Detig RB, Zhou J, Clark AG, Hartl DL, Ayroles JF. Genetic incompatibilities are widespread within species. Nature. 2013;504: 135–137.
Not in our library10.1038/nature12678
49.
Langley CH, Crepeau M, Cardeno C, Corbett-Detig R, Stevens K. Circumventing heterozygosity: sequencing the amplified genome of a single haploid Drosophila melanogaster embryo. Genetics. 2011;188: 239–246.
Not in our library10.1534/genetics.111.127530
50.
Sequencing of Pooled DNA Samples (Pool-Seq) Uncovers Complex Dynamics of Transposable Element Insertions in Drosophila melanogaster · 2012
No notice
51.
Adaptation of Drosophila to a novel laboratory environment reveals temporally heterogeneous trajectories of selected alleles · 2012
No notice
52.
Kolaczkowski B, Kern AD, Holloway AK, Begun DJ. Genomic differentiation between temperate and tropical Australian populations of Drosophila melanogaster. Genetics. 2011;187: 245–260.
Not in our library10.1534/genetics.110.123059
53.
Kapun M, Schalkwyk H, McAllister B, Flatt T, Schlotterer C. Inference of chromosomal inversion dynamics from Pool‐Seq data in natural and laboratory populations of Drosophila melanogaster. Mol Ecol. 2014;23: 1813–1827.
Not in our library10.1111/mec.12594
54.
Zhu Y, Bergland AO, González J, Petrov DA. Empirical validation of pooled whole genome population re-sequencing in Drosophila melanogaster. PLoS ONE. 2012;7: e41901
Not in our library10.1371/journal.pone.0041901
55.
Chen GK, Marjoram P, Wall JD. Fast and flexible simulation of DNA sequence data. Genome Res. 2008;19: 136–142.
Not in our library10.1101/gr.083634.108
56.
Corbett-Detig R, Jones M. SELAM: simulation of epistasis and local adaptation during admixture with mate choice. Bioinformatics. 2016;btw365.
Not in our library10.1093/bioinformatics/btw365
57.
Liang M, Nielsen R. The lengths of admixture tracts. Genetics. 2014;197: 953–967.
Not in our library10.1534/genetics.114.162362
58.
The Date of Interbreeding between Neandertals and Modern Humans · 2012
No notice
59.
The Many Landscapes of Recombination in Drosophila melanogaster · 2012
No notice
60.
Moltke I, Grarup N, Jørgensen ME, Bjerregaard P, Treebak JT, Fumagalli M, et al A common Greenlandic TBC1D4 variant confers muscle insulin resistance and type 2 diabetes. Nature. 2014;512: 190–193.
Not in our library10.1038/nature13425
61.
Moltke I, Fumagalli M, Korneliussen TS. Uncovering the genetic history of the present-day Greenlandic population. Am J Hum Genet. 2015;96: 54–69.
Not in our library10.1016/j.ajhg.2014.11.012
62.
The Metabochip, a Custom Genotyping Array for Genetic Studies of Metabolic, Cardiovascular, and Anthropometric Traits · 2012
Corrected
63.
Population Genomics of Inversion Polymorphisms in Drosophila melanogaster · 2012
Corrected
64.
Kapun M, Fabian DK, Goudet J, Flatt T. Genomic evidence for adaptive inversion clines in Drosophila melanogaster. Mol Biol Evol. 2016;33: 1317–1336.
Not in our library10.1093/molbev/msw016
65.
Mackay TFC, Richards S, Stone EA, Barbadilla A, Ayroles JF, Zhu D, et al The Drosophila melanogaster Genetic Reference Panel. Nature. 2012;482: 173–178.
Not in our library10.1038/nature10811
66.
Natural variation in genome architecture among 205 Drosophila melanogaster Genetic Reference Panel lines · 2014
No notice
67.
Corbett-Detig RB, Cardeno C, Langley CH. Sequence-based detection and breakpoint assembly of polymorphic inversions. Genetics. 2012;192: 131–137.
Not in our library10.1534/genetics.112.141622
68.
The Genomics of Speciation in Drosophila: Diversity, Divergence, and Introgression Estimated Using Low-Coverage Genome Sequencing · 2009
No notice
69.
Krimbas CB, Powell JR. Drosophila Inversion Polymorphism. CRC Press; 1992.
Not in our library
70.
Aulard S, David JR, Lemenieur F. Chromosomal inversion polymorphism in Afrotropical populations of Drosophila melanogaster. Genet Res. 2002;79: 49–63.
Not in our library10.1017/s0016672301005407
71.
Pasaniuc B, Sankararaman S, Torgerson DG, Gignoux C, Zaitlen N, Eng C, et al Analysis of Latino populations from GALA and MEC studies reveals genomic loci with biased local ancestry estimation. Bioinformatics. 2013;29: 1407–1415.
Not in our library10.1093/bioinformatics/btt166
72.
Ting CT, Takahashi A, Wu CI. Incipient speciation by sexual isolation in Drosophila: concurrent evolution at multiple loci. PNAS. 2001;98: 6709–6713.
Not in our library10.1073/pnas.121418898
73.
Hollocher H, Ting CT, Wu ML, Wu CI. Incipient speciation by sexual isolation in Drosophila melanogaster: extensive genetic divergence without reinforcement. Genetics. 1997;147: 1191–1201.
Not in our library10.1093/genetics/147.3.1191
74.
Hollocher H, Ting CT, Pollack F, Wu CI. Incipient speciation by sexual isolation in Drosophila melanogaster: variation in mating preference and correlation between sexes. Evolution. 1997;51: 1175–1181.
Not in our library10.1111/j.1558-5646.1997.tb03965.x
75.
Coyne JA. Genetics and speciation. Nature. 1992;335: 511–515.
Not in our library10.1038/355511a0
76.
Coyne JA, Orr HA. Patterns of speciation in Drosophila. Evolution. 1989;43: 362–381.
Not in our library10.1111/j.1558-5646.1989.tb04233.x
77.
Charlesworth B, Coyne JA, Barton NH. The relative rates of evolution of sex chromosomes and autosomes. American Naturalist. 1987;130: 113–146.
Not in our library
78.
Levine MT, Begun DJ. Evidence of Spatially Varying Selection Acting on Four Chromatin-Remodeling Loci in Drosophila melanogaster. Genetics. 2008;179: 475–485.
Not in our library10.1534/genetics.107.085423
79.
Ellis LL, Carney GE. Socially-Responsive Gene Expression in Male Drosophila melanogaster Is Influenced by the Sex of the Interacting Partner. Genetics. 2011;187: 157–169.
Not in our library10.1534/genetics.110.122754
80.
Edwards AC, Zwarts L, Yamamoto A, Callaerts P, Mackay TF. Mutations in many genes affect aggressive behavior in Drosophila melanogaster. BMC Biol. 2009;7: 29–13.
Not in our library10.1186/1741-7007-7-29
81.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. 1995;57: 289–300.
Not in our library
82.
SNP Calling, Genotype Calling, and Sample Allele Frequency Estimation from New-Generation Sequencing Data · 2012
No notice
83.
Fumagalli M, Vieira FG, Linderoth T, Nielsen R. ngsTools: methods for population genetics analyses from next-generation sequencing data. Bioinformatics. 2014;30: 1486–1487.
Not in our library10.1093/bioinformatics/btu041
84.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43: 491–498.
Not in our library10.1038/ng.806
85.
Baird SJE, Barton NH, Etheridge AM. The distribution of surviving blocks of an ancestral genome. Theoretical Population Biology. 2003;64: 451–471.
Not in our library10.1016/s0040-5809(03)00098-4
86.
Genomic Ancestry of North Africans Supports Back-to-Africa Migrations · 2012
No notice
87.
Fast "coalescent" simulation · 2006
No notice
88.
Kiefer J. Sequential minimax search for a maximum. Proceedings of the American Mathematical Society. 1953;4: 502–506.
Not in our library
89.
Inferring the Joint Demographic History of Multiple Populations from Multidimensional SNP Frequency Data · 2009
No notice
90.
Pool JE, Aquadro CF. History and Structure of Sub-Saharan Populations of Drosophila melanogaster. Genetics. 2006;174: 915–929.
Not in our library10.1534/genetics.106.058693
91.
Lack JB, Cardeno CM, Crepeau MW, Taylor W, Corbett-Detig RB, Stevens KA, et al The Drosophila Genome Nexus: a population genomic resource of 623 Drosophila melanogaster genomes, including 197 from a single ancestral range population. Genetics. 2015;199: 1229–1241.
Not in our library10.1534/genetics.115.174664
92.
Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv preprint arXiv:13033997. 2013.
Not in our library
93.
Adams MD, Celniker SE, Holt RA, Evans CA, Gocayne JD, Amanatides PG. The genome sequence of Drosophila melanogaster. Science. 2000;287: 2185–2195.
Not in our library10.1126/science.287.5461.2185
94.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25: 2078–2079.
Not in our library10.1093/bioinformatics/btp352
95.
Best DJ, Roberts DE. Algorithm AS 89: the upper tail probabilities of Spearman's rho. Journal of the Royal Statistical Society Series C. 1975;24: 377–379.
Not in our library
96.
Kao JY, Lymer S, Hwang SH, Sung A, Nuzhdin SV. Postmating reproductive barriers contribute to the incipient sexual isolation of the United States and Caribbean Drosophila melanogaster. Ecol Evol. 2015;5: 3171–3182.
Not in our library10.1002/ece3.1596
97.
Kofler R, Schlötterer C. Gowinda: unbiased analysis of gene set enrichment for genome-wide association studies. Bioinformatics. 2012;28: 2084–2085.
Not in our library10.1093/bioinformatics/bts315

Cited by papers we hold 3

The genomic consequences of hybridization · 2021
No notice
Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass · 2021
No notice
The Unreasonable Effectiveness of Convolutional Neural Networks in Population Genetic Inference · 2018
No notice

Reads

Read by people on this page

None yet.

Checked or read by AI through naulon

None yet.

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.

Are you an author?

Is this your paper?

Claim it to add your note and be paid when agents cite it. Free.