The study generates the first ancient DNA from the northern Indian Himalayas: 8 usable ancient genomes (~3444-21 cal yr B.P.) from Himachal Pradesh, Ladakh, and Mustang, plus 10 present-day Himachal individuals sequenced with community engagement. Three well-covered ancient individuals (Sanglung 16.6Γ; Kibber 26.5Γ; Tashigang) modeled via qpAdm as ~80.6-85.8% Tibetan-related and ~14.2-19.4% Afanasievo/Steppe-related ancestry, with DATES dating the admixture to ~2700-3800 yr B.P. Present-day HIM-A (Lahaul/Spiti) shows genetic continuity with Sanglung over ~1300 years, while two Kinnaur individuals (HIM-B, HIM-C) carry ~27.6-37.9% and ~68.9-83.4% ANI-ASI South Asian ancestry respectively, dated ~2400-3000 and ~700 yr B.P. Denisovan-derived EPAS1 introgression and EGLN1/TIMP3 selection signals confirm Tibetan-ancestry-mediated high-altitude adaptation despite admixture .
Point estimates from reported qpAdm ranges; HIM-A modeled as 100% Sanglung-related (continuity); HIM-B/HIM-C values are range midpoints. Steppe/South Asian split for HIM-B reflects different tested two-source models.
Strengths include rigorous authenticity controls (damage patterns, trimming sensitivity analyses, PMD filtering, <5% mitochondrial contamination with outlier exclusion), complementary methods (PCA, f3, D-statistics, qpAdm/qpWave, DATES, FEEMS, MSMC2, ROH), exemplary community engagement with Bhoti-language consent and results dissemination, and transparent data deposition (SRA PRJNA1250143; eGA EGAS50000001342). The authors candidly flag key limitations: very small sample sizes (n=8 ancient, n=10 present-day), low coverage for Kanam/Lippa/Choskar precluding qpAdm modeling, and large standard errors on DATES admixture estimates .
Remaining blind spots: (1) the Tibetan source group Laga is itself low-coverage, weakening source attribution; (2) Afanasievo is a distant proxyβthe true Steppe vector could be any related Central Asian group, and the authors cannot distinguish direct Steppe contact versus mediation through already-admixed Tibetan groups; (3) Tashigang's admixture date is nonsignificant and single-source models cannot be excluded; (4) the D-statistic depends on Lubrak as comparator, which shares ancestry with targets, limiting resolution; (5) two ancient individuals' contamination estimates were too imprecise to compute. The TIMP3 selection signal at rs5998622 is intriguing but rests on eQTL annotation from GTEx without functional validationβappropriately framed as hypothesis-generating .
The authors specify falsification: future genomes showing exclusively Tibetan ancestry with no Steppe/South Asian signal across ancient and modern samples would refute the central claim. High-coverage genomes from the Zhangzhung-period region or earlier northern Himalayan sites, and better source populations than Laga, could substantially revise the admixture timing and route. One notable conflict-of-interest disclosure: co-author M. Aldenderfer is a deputy editor at Science Advances (the publishing journal)βa transparency point, not evidence of misconduct .
Know what changed, what holds up, and what remains uncertain. Every Friday. No ads.