Cited by 28

The DOI matches its Crossref record; 15 of the 15 references it cites were found in a registry or an index.

Detecting Signatures of Positive Selection along Defined Branches of a Population Tree Using LSD

Librado, Pablo, Orlando, Ludovic·2018·Molecular Biology and Evolution

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1093/molbev/msy053

Cites

15 of 15 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1. Overview of methods for identifying lineage-specific adaptations. (A) PBS calculations for population A, applying the Gromov product on pairs of FST indexes. (B) The time spanned from the MRCA to the split of populations A and B determines the total branch length on which mutations could have been accumulated in the ancestral population, and thus its diversity. (C) Five population tree, with popA-popE representing external populations, AB and CD internal balanced branches, and ABCD an internal unbalanced branch. (D) Collection of gene trees, each summarizing local evolutionary forces operating on particular genomic regions, potentially including lineage-specific selection. (E) Local gene tree highlighting the footprint of selection on the individuals belonging to population A. The right panel provides a zoom into the population A undergoing selection at a particular locus, and its sister population B, altogether with the equations underlying the calculation of LSD(A). Note this local gene tree exemplifies strong selection signatures, characterized by low intrapopulation diversity and high interpopulation divergence. LSD can also be calculated from sister clades that are not reciprocally monophyletic.
Fig. 1. Overview of methods for identifying lineage-specific adaptations. (A) PBS calculations for population A, applying the Gromov product on pairs of FST indexes. (B) The time spanned from the MRCA to the split of populations A and B determines the total branch length on which mutations could have been accumulated in the ancestral population, and thus its diversity. (C) Five population tree, with popA-popE representing external populations, AB and CD internal balanced branches, and ABCD an internal unbalanced branch. (D) Collection of gene trees, each summarizing local evolutionary forces operating on particular genomic regions, potentially including lineage-specific selection. (E) Local gene tree highlighting the footprint of selection on the individuals belonging to population A. The right panel provides a zoom into the population A undergoing selection at a particular locus, and its sister population B, altogether with the equations underlying the calculation of LSD(A). Note this local gene tree exemplifies strong selection signatures, characterized by low intrapopulation diversity and high interpopulation divergence. LSD can also be calculated from sister clades that are not reciprocally monophyletic.CC BY-NC · PubMed Central
Fig. 2. Performance of LSDnorm in trifurcating population trees. Fifty chromosomes (25 diploid individuals) were sampled for populations A and B. The initial frequency and dominance coefficient of beneficial mutations were assumed to be 0.05 and 0.5, respectively. (A) The left panel shows the simulated scenario, where beneficial alleles were introduced in population A. The x–y scatterplot shows that only LSDnorm(A) scales with increasing selection coefficients. (B) The same adaptive mutations appeared, independently, in populations A and B. LSDnorm(A) reached lower values than in (A). (C) Convergent evolution from variants that were already segregating in the ancestral population, and became beneficial after the split of populations A and B.
Fig. 2. Performance of LSDnorm in trifurcating population trees. Fifty chromosomes (25 diploid individuals) were sampled for populations A and B. The initial frequency and dominance coefficient of beneficial mutations were assumed to be 0.05 and 0.5, respectively. (A) The left panel shows the simulated scenario, where beneficial alleles were introduced in population A. The x–y scatterplot shows that only LSDnorm(A) scales with increasing selection coefficients. (B) The same adaptive mutations appeared, independently, in populations A and B. LSDnorm(A) reached lower values than in (A). (C) Convergent evolution from variants that were already segregating in the ancestral population, and became beneficial after the split of populations A and B.CC BY-NC · PubMed Central
Fig. 3. Comparing LSDnorm.phased, LSDnorm.unphased, and PBS, across a wide range of initial allele frequencies (iAFs) and sample sizes. Intermediate dominance was assumed. The y-axis measures the difference between the scores estimated at s = 0.03 and s = 0.00, as proxies for the discriminatory power of LSDnorm.phased, LSDnorm.unphased, and PBS. LSDnorm.phased was generally associated with increased accuracy, compared with LSDnorm.unphased, and this, in turn, to PBS. Wider confidence intervals for LSDnorm than for PBS were partly due to different sample sizes from the outgroup population.
Fig. 3. Comparing LSDnorm.phased, LSDnorm.unphased, and PBS, across a wide range of initial allele frequencies (iAFs) and sample sizes. Intermediate dominance was assumed. The y-axis measures the difference between the scores estimated at s = 0.03 and s = 0.00, as proxies for the discriminatory power of LSDnorm.phased, LSDnorm.unphased, and PBS. LSDnorm.phased was generally associated with increased accuracy, compared with LSDnorm.unphased, and this, in turn, to PBS. Wider confidence intervals for LSDnorm than for PBS were partly due to different sample sizes from the outgroup population.CC BY-NC · PubMed Central
Fig. 4. Performance of LSDnorm in five-population trees. Fifty chromosomes (25 diploid individuals) were sampled for all populations. The initial frequency and dominance coefficient of beneficial mutations were assumed to be 0.05 and 0.5, respectively. (A) The left panel shows the simulated scenario, where beneficial alleles are introduced in the unbalanced lineage ABCD. The x–y scatterplot shows that only LSDnorm(ABCD) scales with increasing selection coefficients. (B) Beneficial alleles, associated with a range of selection coefficients, are introduced in population A, increasing LSDnorm(A) scores. (C) The unbalanced lineage ABCD experiences introgression from its sister lineage, leading to population E. As the proportion of migrants increases towards one, LSDnorm(ABCD) approaches to zero, which is indicative of the absence of genetic differentiation from its sister lineage. (D) Population A experiences introgression from the distant population population E. LSDnorm(ABCD) values can eventually become negative.
Fig. 4. Performance of LSDnorm in five-population trees. Fifty chromosomes (25 diploid individuals) were sampled for all populations. The initial frequency and dominance coefficient of beneficial mutations were assumed to be 0.05 and 0.5, respectively. (A) The left panel shows the simulated scenario, where beneficial alleles are introduced in the unbalanced lineage ABCD. The x–y scatterplot shows that only LSDnorm(ABCD) scales with increasing selection coefficients. (B) Beneficial alleles, associated with a range of selection coefficients, are introduced in population A, increasing LSDnorm(A) scores. (C) The unbalanced lineage ABCD experiences introgression from its sister lineage, leading to population E. As the proportion of migrants increases towards one, LSDnorm(ABCD) approaches to zero, which is indicative of the absence of genetic differentiation from its sister lineage. (D) Population A experiences introgression from the distant population population E. LSDnorm(ABCD) values can eventually become negative.CC BY-NC · PubMed Central
Fig. 5. LSDnorm applied to the 1000 human genome project (phase III). (A) Human population tree assumed to investigate selection along the lineage ancestral to Europeans (CEU) and Asians (PJB and CHB) (red arrow), and the European population (green arrow). (B) Histogram showing the empirical LSDnorm distribution for the CEU population. The black arrow highlights the excess of 10-kb windows with a maximum score of one. (C) LSDnorm scores along the SH2B3-ATXN2 locus, altogether with the local gene tree at chr12: 111,440,001–111,450,000, which shows widely spread derived alleles amongst CEU (green), PJL (dark yellow) and CHB (dark brown) populations.
Fig. 5. LSDnorm applied to the 1000 human genome project (phase III). (A) Human population tree assumed to investigate selection along the lineage ancestral to Europeans (CEU) and Asians (PJB and CHB) (red arrow), and the European population (green arrow). (B) Histogram showing the empirical LSDnorm distribution for the CEU population. The black arrow highlights the excess of 10-kb windows with a maximum score of one. (C) LSDnorm scores along the SH2B3-ATXN2 locus, altogether with the local gene tree at chr12: 111,440,001–111,450,000, which shows widely spread derived alleles amongst CEU (green), PJL (dark yellow) and CHB (dark brown) populations.CC BY-NC · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 15 references

Research works printed but not linked: 57

  1. Abbas, 1983, Persistence of high intestinal lactase activity in Pakistan, Hum Genet, 64, 277, 10.1007/BF00279410 · DOI
  2. Agyemang, 2003, Is the blood pressure of people from African origin adults in the UK higher or lower than that in European origin white people? A review of cross-sectional data, J Hum Hypertens, 17, 523, 10.1038/sj.jhh.1001586 · DOI
  3. Auburger, 2014, 12q24 locus association with type 1 diabetes: SH2B3 or ATXN2?, World J Diabetes, 5, 316, 10.4239/wjd.v5.i3.316 · DOI
  4. Berg, 2015, A coalescent model for a sweep of a unique standing variant, Genetics, 201, 707, 10.1534/genetics.115.178962 · DOI
  5. Bersaglieri, 2004, Genetic signatures of strong recent positive selection at the lactase gene, Am J Hum Genet, 74, 1111, 10.1086/421051 · DOI
  6. Blanton, 1991, Linkage mapping of autosomal dominant retinitis pigmentosa (RP1) to the pericentric region of human chromosome 8, Genomics, 11, 857, 10.1016/0888-7543(91)90008-3 · DOI
  7. Bonhomme, 2010, Detecting selection in population trees: The Lewontin and Krakauer test extended, Genetics, 186, 241, 10.1534/genetics.110.117275 · DOI
  8. Brinkworth, 2014, The contribution of natural selection to present-day susceptibility to chronic inflammatory and autoimmune disease, Curr Opin Immunol, 31, 66, 10.1016/j.coi.2014.09.008 · DOI
  9. Chen, 2010, Population differentiation as a test for selective sweeps, Genome Res, 20, 393, 10.1101/gr.100545.109 · DOI
  10. Clément, 2009, In obese and non-obese adults, the cis-regulatory rs361072 promoter variant of PIK3CB is associated with insulin resistance not with type 2 diabetes, Mol Genet Metab, 96, 129, 10.1016/j.ymgme.2008.11.160 · DOI
  11. Colagiuri, 2002, The “carnivore connection”–evolutionary aspects of insulin resistance, Eur J Clin Nutr, 56, S30, 10.1038/sj.ejcn.1601351 · DOI
  12. DeGiorgio, 2016, SweepFinder2: Increased sensitivity, robustness and flexibility, Bioinformatics, 32, 1895, 10.1093/bioinformatics/btw051 · DOI
  13. Fan, 2016, Going global by adapting local: A review of recent human adaptation, Science, 354, 54, 10.1126/science.aaf5098 · DOI
  14. Fariello, 2013, Detecting signatures of selection through haplotype differentiation among hierarchically structured populations, Genetics, 193, 929, 10.1534/genetics.112.147231 · DOI
  15. Ferrer-Admetlla, 2014, On detecting incomplete soft or hard selective sweeps using haplotype structure, Mol Biol Evol, 31, 1275, 10.1093/molbev/msu077 · DOI
  16. Field, 2016, Detection of human adaptation during the past 2000 years, Science, 354, 760, 10.1126/science.aag0776 · DOI · Paper page
  17. Garud, 2015, Recent selective sweeps in North American Drosophila melanogaster show signatures of soft sweeps, PLoS Genet., 11, e1005004., 10.1371/journal.pgen.1005004 · DOI
  18. Green, 2010, A draft sequence of the Neandertal genome, Science, 328, 710, 10.1126/science.1188021 · DOI · Paper page
  19. Grün, 2005, U-series and ESR analyses of bones and teeth relating to the human burials from Skhul, J Hum Evol, 49, 316, 10.1016/j.jhevol.2005.04.006 · DOI · Paper page
  20. Hermisson, 2005, Soft sweeps: Molecular population genetics of adaptation from standing genetic variation, Genetics, 169, 2335, 10.1534/genetics.104.036947 · DOI
  21. Hershkovitz, 2018, The earliest modern humans outside Africa, Science, 359, 456, 10.1126/science.aap8369 · DOI
  22. Heyer, 2011, Lactase persistence in central Asia: Phenotype, genotype, and evolution, Hum Biol, 83, 379, 10.3378/027.083.0304 · DOI
  23. Hofer, 2009, Large allele frequency differences between human continental groups are more likely to have occurred by drift during range expansions than by selection, Ann Hum Genet, 73, 95, 10.1111/j.1469-1809.2008.00489.x · DOI
  24. Hudson, 1987, A test of neutral molecular evolution based on nucleotide data, Genetics, 116, 153, 10.1093/genetics/116.1.153 · DOI
  25. Hudson, 1992, Estimation of levels of gene flow from DNA sequence data, Genetics, 132, 583, 10.1093/genetics/132.2.583 · DOI
  26. Ihara, 2007, Sept4, a component of presynaptic scaffold and Lewy bodies, is required for the suppression of alpha-synuclein neurotoxicity, Neuron, 53, 519, 10.1016/j.neuron.2007.01.019 · DOI
  27. Kaplan, 1988, The coalescent process in models with selection, Genetics, 120, 819, 10.1093/genetics/120.3.819 · DOI
  28. Kaplan, 1989, The “hitchhiking effect” revisited, Genetics, 123, 887, 10.1093/genetics/123.4.887 · DOI · Paper page
  29. Kim, 2002, Detecting a local signature of genetic hitchhiking along a recombining chromosome, Genetics, 160, 765, 10.1093/genetics/160.2.765 · DOI
  30. Kullo, 2014, The ATXN2-SH2B3 locus is associated with peripheral arterial disease: An electronic medical record-based genome-wide association study, Front Genet, 5, 166., 10.3389/fgene.2014.00166 · DOI
  31. Lefort, 2015, FastME 2.0: A comprehensive, accurate, and fast distance-based phylogeny inference program, Mol Biol Evol, 32, 2798, 10.1093/molbev/msv150 · DOI
  32. Lewontin, 1973, Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms, Genetics, 74, 175, 10.1093/genetics/74.1.175 · DOI
  33. Marciniak, 2017, Harnessing ancient genomes to study the history of human adaptation, Nat. Rev. Genet, 18, 659, 10.1038/nrg.2017.65 · DOI
  34. McDonald, 1991, Adaptive protein evolution at the Adh locus in Drosophila, Nature, 351, 652, 10.1038/351652a0 · DOI
  35. Nei, 1975, Lewontin-Krakauer test for neutral genes, Genetics, 80, 395., 10.1093/genetics/80.2.395 · DOI
  36. Nielsen, 2017, Tracing the peopling of the world through genomics, Nature, 541, 302, 10.1038/nature21347 · DOI · Paper page
  37. Orrù, 2013, Genetic variants regulating immune cell levels in health and disease, Cell, 155, 242, 10.1016/j.cell.2013.08.041 · DOI
  38. Pace, 2014, Leishmaniasis, J. Infect, 69(1 Suppl), S10, 10.1016/j.jinf.2014.07.016 · DOI
  39. Patterson, 2012, Ancient admixture in human history, Genetics, 192, 1065, 10.1534/genetics.112.145037 · DOI · Paper page
  40. Pavlidis, 2013, SweeD: Likelihood-based detection of selective sweeps in thousands of genomes, Mol Biol Evol, 30, 2224, 10.1093/molbev/mst112 · DOI
  41. Pierce, 1999, Mutations in a gene encoding a new oxygen-regulated photoreceptor protein cause dominant retinitis pigmentosa, Nat Genet, 22, 248, 10.1038/10305 · DOI
  42. Price, 2010, FastTree 2—Approximately maximum-likelihood trees for large alignments, PLoS One, 5, e9490., 10.1371/journal.pone.0009490 · DOI · Paper page
  43. Racimo, 2016, Testing for ancient selection using cross-population allele frequency differentiation, Genetics, 202, 733, 10.1534/genetics.115.178095 · DOI · Paper page
  44. Ronquist, 2003, MrBayes 3: Bayesian phylogenetic inference under mixed models, Bioinforma Oxf Engl, 19, 1572, 10.1093/bioinformatics/btg180 · DOI
  45. Shriver, 2004, The genomic distribution of population substructure in four populations using 8,525 autosomal SNPs, Hum Genomics, 1, 274., 10.1186/1479-7364-1-4-274 · DOI
  46. Sikora, 2017, Ancient genomes show social and reproductive behavior of early Upper Paleolithic foragers, Science, 358, 659, 10.1126/science.aao1807 · DOI · Paper page
  47. Smith, 1974, The hitch-hiking effect of a favourable gene, Genet Res, 23, 23, 10.1017/S0016672300014634 · DOI · Paper page
  48. Stamatakis, 2014, RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies, Bioinformatics, 30, 1312, 10.1093/bioinformatics/btu033 · DOI · Paper page
  49. Sun, 2017, Adaptation to deep-sea chemosynthetic environments as revealed by mussel genomes, Nat Ecol Evol, 1, 121, 10.1038/s41559-017-0121 · DOI
  50. Tataru, 2017, Statistical Inference in the Wright–Fisher Model Using Allele Frequency Data, e30
  51. The 1000 Genomes Project Consortium, 2015, A global reference for human genetic variation, Nature, 526, 68, 10.1038/nature15393 · DOI · Paper page
  52. Vatsiou, 2016, Changes in selective pressures associated with human population expansion may explain metabolic and immune related pathways enriched for signatures of positive selection, BMC Genomics, 17, 504., 10.1186/s12864-016-2783-2 · DOI
  53. Vitti, 2013, Detecting natural selection in genomic data, Annu Rev Genet, 47, 97, 10.1146/annurev-genet-111212-133526 · DOI · Paper page
  54. Voight, 2006, A map of recent positive selection in the human genome, PLoS Biol, 4, e72., 10.1371/journal.pbio.0040072 · DOI
  55. Yang, 1997, Bayesian phylogenetic inference using DNA sequences: A Markov chain Monte Carlo method, Mol Biol Evol, 14, 717, 10.1093/oxfordjournals.molbev.a025811 · DOI · Paper page
  56. Zhang, 2005, WebGestalt: An integrated system for exploring gene sets in various biological contexts, Nucleic Acids Res., 33, W741, 10.1093/nar/gki475 · DOI
  57. Zhang, 2013, Signatures of natural selection on genetic variants affecting complex human traits, Appl. Transl. Genomics, 2, 78, 10.1016/j.atg.2013.10.002 · DOI

Cited by

Sorted by citations

  • Ancient DNA analysis

    Orlando, L., Allaby, R., Skoglund, P., Der Sarkissian, C., Stockhammer, P. W., Ávila-Arcos, M. C., Fu, Q., Krause, J., Willerslev, E., Stone, A. C., et al. (2021) · Nature Reviews Methods Primers · doi:10.1038/s43586-020-00011-0 · Cited by 411

    Self-citation

Metadata

Authors
Librado, P., & Orlando, L.
Affiliations
Centre for GeoGenetics, Natural History Museum of Denmark, Copenhagen, Denmark
Venue
Molecular Biology and Evolution
Published
2018-03-29
Volume
35
Issue
6
Pages
1520-1535
ISSN
  • 0737-4038
  • 1537-1719
Cited by
28
References
15
Language
en

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review