Reply to Barton et al: signatures of natural selection during the Black Death
Vilgalys, T. P.; Klunk, J.; Demeure, C. E.; Cheng, X.; Shiratori, M.; Madej, J.; Beau, R.; Elli, D.; Patino, M. I.; Redfern, R.; DeWitte, S. N.; Gamble, J.; Boldsen, J. L.; Carmichael, A.; Varlik, N.; Eaton, K.; Grenier, J.-C.; Golding, G. B.; Devault, A.; Rouillard, J.-M.; Yotova, V.; Sindeaux, R.; Ye, C. J.; Bikaran, M.; Dumaine, A.; Brinkworth, J. F.; Missiakas, D.; Rouleau, G. A.; Steinrücken, M.; Pizarro-Cerda, J.; Poinar, H.; Barreiro, L. B.
Show abstract
Barton et al.1 raise several statistical concerns regarding our original analyses2 that highlight the challenge of inferring natural selection using ancient genomic data. We show here that these concerns have limited impact on our original conclusions. Specifically, we recover the same signature of enrichment for high FST values at the immune loci relative to putatively neutral sites after switching the allele frequency estimation method to a maximum likelihood approach, filtering to only consider known human variants, and down-sampling our data to the same mean coverage across sites. Furthermore, using permutations, we show that the rs2549794 variant near ERAP2 continues to emerge as the strongest candidate for selection (p = 1.2x10-5), falling below the Bonferroni-corrected significance threshold recommended by Barton et al. Importantly, the evidence for selection on ERAP2 is further supported by functional data demonstrating the impact of the ERAP2 genotype on the immune response to Y. pestis and by epidemiological data from an independent group showing that the putatively selected allele during the Black Death protects against severe respiratory infection in contemporary populations.
Matching journals
The top 5 journals account for 50% of the predicted probability mass.
Similar papers in this journal
Similar papers in this journal
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.