Genetic diversity is a key aspect of biodiversity and has important implications for evolutionary potential and thereby the persistence of species. Improving our understanding of the factors that drive genetic structure within and between populations is, therefore, a long-standing goal in evolutionary biology. However, this is a major challenge, because of the complex interplay between genetic drift, migration, and extinction/colonization dynamics on the one hand, and the biology and ecology of species on the other hand (Romiguier et al. 2014, Ellegren and Galtier 2016, Charlesworth 2003). 

Jarne et al. (2021) studied whether environmental and demographic factors affect the population genetic structure of four species of hermaphroditic freshwater snails in a similar way, using comparative analyses of neutral genetic microsatellite markers. 

Specifically, they investigated microsatellite variability of Hygrophila in almost 280 sites in Guadeloupe, Lesser Antilles, as part of a long-term survey experiment (Lamy et al. 2013). They then modelled the influence of the mating system, local environmental characteristics and demographic factors on population genetic diversity.

Consistent with theoretical predictions (Charlesworth 2003), they detected higher genetic variation in two outcrossing species than in two selfing species, emphasizing the importance of the mating system in maintaining genetic diversity. The study further identified an important role of site connectivity, through its influences on effective population size and extinction/colonisation events. Finally, the study detects an influence of interspecific interactions caused by an ongoing invasion by one of the studied species on genetic structure, highlighting the indirect effect of changes in community composition and demography on population genetics.

Jarne et al. (2021) could address the extent to which genetic structure is determined by demographic and environmental factors in multiple species given the remarkable sampling available. Additionally, the study system is extremely suitable to address this hypothesis as species’ habitats are defined and delineated. Whilst the authors did attempt to test for across-species correlations, further investigations on this matter are required. Moreover, the effect of interactions between factors should be appropriately considered in any modelling between genetic structure and local environmental or demographic features.

The findings in this study contribute to improving our understanding of factors influencing population genetic diversity, and highlights the complexity of interacting factors, therefore also emphasizing the challenges of drawing general implications, additionally hampered by the relatively limited number of species studied. Jarne et al. (2021) provide an excellent showcase of an empirical framework to test determinants of genetic structure in natural populations. As such, this study can be an example for further attempts of comparative analysis of genetic diversity.

References

Charlesworth, D. (2003) Effects of inbreeding on the genetic diversity of populations. Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences, 358, 1051-1070. doi: https://doi.org/10.1098/rstb.2003.1296

Ellegren, H. and Galtier, N. (2016) Determinants of genetic diversity. Nature Reviews Genetics, 17, 422-433. doi: https://doi.org/10.1038/nrg.2016.58

Jarne, P., Lozano del Campo, A., Lamy, T., Chapuis, E., Dubart, M., Segard, A., Canard, E., Pointier, J.-P. and David, P. (2021) Connectivity and selfing drives population genetic structure in a patchy landscape: a comparative approach of four co-occurring freshwater snail species. HAL, hal-03295242, ver. 2 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://hal.archives-ouvertes.fr/hal-03295242

Lamy, T., Gimenez, O., Pointier, J. P., Jarne, P. and David, P. (2013). Metapopulation dynamics of species with cryptic life stages. The American Naturalist, 181, 479-491. doi: https://doi.org/10.1086/669676

Romiguier, J., Gayral, P., Ballenghien, M. et al. (2014) Comparative population genomics in animals uncovers the determinants of genetic diversity. Nature, 515, 261-263. doi: https://doi.org/10.1038/nature13685

​​Durability of plant resistance to diploid pathogen Hirohisa Kishino Based on the population genetic and epidemiologic model, Saubin et al. (2021) report that the resistant hosts generated by the breeding based on the gene-for-gene interaction is durable much longer against diploid pathogens than haploid pathogens. The avr allele of pathogen that confers the resistance is genetically recessive. The heterozygotes are not recognized by the resistant hosts and only the avr/avr homozygote is adaptive. As a result, the trajectory of avr allele frequency becomes more stochastic due to genetic drift. Although the paper focuses on the evolution of standing polymorphism, it seems obvious that the adaptive mutations in pathogen have much larger probability of being deleted from the population because the individuals own the avr allele mostly in the form of heterozygote at the initial phase after the mutation. Since only few among many models of plant resistance deployment study the case of diploid pathogen and the contribution of the pathogen life cycle, this work will add an important intellect to the literature (Rimbaud et al. 2021).

From the study of host-parasite interaction in flax rust Melampsora lini, Flor (1942, 1955) showed that the host resistance is formed by the interaction of a host resistance gene and a corresponding pathogen gene. This gene-for-gene hypothesis has been supported by experimental evidence and has served as a basis of the methods of molecular breeding targeting the dominant R genes. However, modern agriculture provides the pathogen populations with the homogeneous environments and laid strong selection pressure on them. As a result, the newly developed resistant plants face the risk of immediate resistance breakdown (Möller and Stukenbrock 2017).

Currently, quantitative resistance is getting attention as characters as a potential target for long-life (mild) resistant breeds (Lannou, 2012). They are polygenic and controlled partly by the same genes that mediate qualitative resistance but mostly by the genes that encode defense-related outputs such as strengthening of the cell wall or defense compound biosynthesis (Corwin and Kliebenstein, 2017). Progress of molecular genetics may overcome the technical difficulty (Bakkeren and Szabo, 2020). Saubin et al. (2021) notes that the pattern of genetic inheritance of the pathogen counterparts that respond to the host traits is crucial regarding with the durability of the resistant hosts. The resistance traits for which avr alleles are predicted to be recessive may be the targets of breeding.

References

Bakkeren, G., and Szabo, L. J. (2020) Progress on molecular genetics and manipulation of rust fungi. Phytopathology, 110, 532-543. https://doi.org/10.1094/PHYTO-07-19-0228-IA

Corwin, J. A., and Kliebenstein, D. J. (2017) Quantitative resistance: more than just perception of a pathogen. The Plant Cell, 29, 655-665. https://doi.org/10.1105/tpc.16.00915

Flor, H. H. (1942) Inheritance of pathogenicity in a cross between physiological races 22 and 24 of Melampsova lini. Phytopathology, 35. Abstract.

Flor, H. H. (1955) Host-parasite interactions in flax rust-its genetics and other implications. Phytopathology, 45, 680-685.

Lannou, C. (2012) Variation and selection of quantitative traits in plant pathogens. Annual review of phytopathology, 50, 319-338. https://doi.org/10.1146/annurev-phyto-081211-173031

Möller, M. and Stukenbrock, E. H. (2017) Evolution and genome architecture in fungal plant pathogens. Nature Reviews Microbiology. 15, 756–771. https://doi.org/10.1038/nrmicro.2017.76

Rimbaud, L., Fabre, F., Papaïx, J., Moury, B., Lannou, C., Barrett, L. G., and Thrall, P. H. (2021) Models of Plant Resistance Deployment. Annual Review of Phytopathology, 59. https://doi.org/10.1146/annurev-phyto-020620-122134

Saubin, M., De Mita, S., Zhu, X., Sudret, B. and Halkett, F. (2021) Impact of ploidy and pathogen life cycle on resistance durability. bioRxiv, 2021.05.28.446112, ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.05.28.446112

In recent years many studies have documented shifts in phenology in response to climate change, be it in arrival times in migrating birds, budset in trees, adult emergence in butterflies, or flowering time in annual plants (Coen et al. 2018; Piao et al. 2019). While these changes are, in part, explained by phenotypic plasticity, more and more studies find that they involve also genetic changes, that is, they involve evolutionary change (e.g., Metz et al. 2020). Yet, evolutionary change may occur through genetic drift as well as selection. Therefore, in order to demonstrate adaptive evolutionary change in response to climate change, drift has to be excluded as an alternative explanation (Hansen et al. 2012). A new study by Gay et al. (2021) shows just how difficult this can be. 

The authors investigated a recent evolutionary shift in flowering time by in a population an annual plant that reproduces predominantly by self-fertilization. The population has recently been subjected to increased temperatures and reduced rainfalls both of which are believed to select for earlier flowering times. They used a “resurrection” approach (Orsini et al. 2013; Weider et al. 2018): Genotypes from the past (resurrected from seeds) were compared alongside more recent genotypes (from more recently collected seeds) under identical conditions in the greenhouse. Using an experimental design that replicated genotypes, eliminated maternal effects, and controlled for microenvironmental variation, they found said genetic change in flowering times: Genotypes obtained from recently collected seeds flowered significantly (about 2 days) earlier than those obtained 22 generations before. However, neutral markers (microsatellites) also showed strong changes in allele frequencies across the 22 generations, suggesting that effective population size, Ne, was low (i.e., genetic drift was strong), which is typical for highly self-fertilizing populations. In addition, several multilocus genotypes were present at high frequencies and persisted over the 22 generations, almost as in clonal populations (e.g., Schaffner et al. 2019). The challenge was thus to evaluate whether the observed evolutionary change was the result of an adaptive response to selection or may be explained by drift alone. 

Here, Gay et al. (2021) took a particularly careful and thorough approach. First, they carried out a selection gradient analysis, finding that earlier-flowering plants produced more seeds than later-flowering plants. This suggests that, under greenhouse conditions, there was indeed selection for earlier flowering times. Second, investigating other populations from the same region (all populations are located on the Mediterranean island of Corsica, France), they found that a concurrent shift to earlier flowering times occurred also in these populations. Under the hypothesis that the populations can be regarded as independent replicates of the evolutionary process, the observation of concurrent shifts rules out genetic drift (under drift, the direction of change is expected to be random). 

The study may well have stopped here, concluding that there is good evidence for an adaptive response to selection for earlier flowering times in these self-fertilizing plants, at least under the hypothesis that selection gradients estimated in the greenhouse are relevant to field conditions. However, the authors went one step further. They used the change in the frequencies of the multilocus genotypes across the 22 generations as an estimate of realized fitness in the field and compared them to the phenotypic assays from the greenhouse. The results showed a tendency for high-fitness genotypes (positive frequency changes) to flower earlier and to produce more seeds than low-fitness genotypes. However, a simulation model showed that the observed correlations could be explained by drift alone, as long as Ne is lower than ca. 150 individuals. The findings were thus consistent with an adaptive evolutionary change in response to selection, but drift could only be excluded as the sole explanation if the effective population size was large enough. 

The study did provide two estimates of Ne (19 and 136 individuals, based on individual microsatellite loci or multilocus genotypes, respectively), but both are problematic. First, frequency changes over time may be influenced by the presence of a seed bank or by immigration from a genetically dissimilar population, which may lead to an underestimation of Ne (Wang and Whitlock 2003). Indeed, the low effective size inferred from the allele frequency changes at microsatellite loci appears to be inconsistent with levels of genetic diversity found in the population. Moreover, high self-fertilization reduces effective recombination and therefore leads to non-independence among loci. This lowers the precision of the Ne estimates (due to a higher sampling variance) and may also violate the assumption of neutrality due to the possibility of selection (e.g., due to inbreeding depression) at linked loci, which may be anywhere in the genome in case of high degrees of self-fertilization. 

There is thus no definite answer to the question of whether or not the observed changes in flowering time in this population were driven by selection. The study sets high standards for other, similar ones, in terms of thoroughness of the analyses and care in interpreting the findings. It also serves as a very instructive reminder to carefully check the assumptions when estimating neutral expectations, especially when working on species with complicated demographies or non-standard life cycles. Indeed the issues encountered here, in particular the difficulty of establishing neutral expectations in species with low effective recombination, may apply to many other species, including partially or fully asexual ones (Hartfield 2016). Furthermore, they may not be limited to estimating Ne but may also apply, for instance, to the establishment of neutral baselines for outlier analyses in genome scans (see e.g, Orsini et al. 2012). 

References

Cohen JM, Lajeunesse MJ, Rohr JR (2018) A global synthesis of animal phenological responses to climate change. Nature Climate Change, 8, 224–228. https://doi.org/10.1038/s41558-018-0067-3

Gay L, Dhinaut J, Jullien M, Vitalis R, Navascués M, Ranwez V, Ronfort J (2021) Evolution of flowering time in a selfing annual plant: Roles of adaptation and genetic drift. bioRxiv, 2020.08.21.261230, ver. 4 recommended and peer-reviewed by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2020.08.21.261230

Hansen MM, Olivieri I, Waller DM, Nielsen EE (2012) Monitoring adaptive genetic responses to environmental change. Molecular Ecology, 21, 1311–1329. https://doi.org/10.1111/j.1365-294X.2011.05463.x

Hartfield M (2016) Evolutionary genetic consequences of facultative sex and outcrossing. Journal of Evolutionary Biology, 29, 5–22. https://doi.org/10.1111/jeb.12770

Metz J, Lampei C, Bäumler L, Bocherens H, Dittberner H, Henneberg L, Meaux J de, Tielbörger K (2020) Rapid adaptive evolution to drought in a subset of plant traits in a large-scale climate change experiment. Ecology Letters, 23, 1643–1653. https://doi.org/10.1111/ele.13596

Orsini L, Schwenk K, De Meester L, Colbourne JK, Pfrender ME, Weider LJ (2013) The evolutionary time machine: using dormant propagules to forecast how populations can adapt to changing environments. Trends in Ecology & Evolution, 28, 274–282. https://doi.org/10.1016/j.tree.2013.01.009

Orsini L, Spanier KI, Meester LD (2012) Genomic signature of natural and anthropogenic stress in wild populations of the waterflea Daphnia magna: validation in space, time and experimental evolution. Molecular Ecology, 21, 2160–2175. https://doi.org/10.1111/j.1365-294X.2011.05429.x

Piao S, Liu Q, Chen A, Janssens IA, Fu Y, Dai J, Liu L, Lian X, Shen M, Zhu X (2019) Plant phenology and global climate change: Current progresses and challenges. Global Change Biology, 25, 1922–1940. https://doi.org/10.1111/gcb.14619

Schaffner LR, Govaert L, De Meester L, Ellner SP, Fairchild E, Miner BE, Rudstam LG, Spaak P, Hairston NG (2019) Consumer-resource dynamics is an eco-evolutionary process in a natural plankton community. Nature Ecology & Evolution, 3, 1351–1358. https://doi.org/10.1038/s41559-019-0960-9

Wang J, Whitlock MC (2003) Estimating Effective Population Size and Migration Rates From Genetic Samples Over Space and Time. Genetics, 163, 429–446. PMID: 12586728

Weider LJ, Jeyasingh PD, Frisch D (2018) Evolutionary aspects of resurrection ecology: Progress, scope, and applications—An overview. Evolutionary Applications, 11, 3–10. https://doi.org/10.1111/eva.12563

Within-host parasite density is tightly linked to parasite fitness often determining both transmission success and virulence (parasite-induced harm to the host) (Alizon et al., 2009, Anderson & May, 1982). Parasite density may thus be controlled by selection balancing these conflicting pressures. Actual within-host density regulation may be under host or parasite control, or due to other environmental factors (Wale et al., 2019, Vale et al., 2011, Chrostek et al., 2013). Vertically transmitted parasites may also be more vulnerable to drift associated with bottlenecks between generations, which may also determine within-host population size (Mathe-Hubert et al., 2019, Mira & Moran, 2002).

Bénard et al. (2021) use 3 experiments to disentangle the role of drift versus host factors in the control of within-host Wolbachia growth in Drosophila melanogaster. They use the wMelPop Wolbachia strain in which virulence (fly longevity) and within-host growth correlate positively with copy number in the genomic region Octomom (Chrostek et al., 2013, Chrostek & Teixeira, 2015). Octomom copy number can be used as a marker for different genetic lineages within the wMelPop strain.

In a first experiment, they introgressed and backcrossed this Wolbachia strain into 6 different host genetic backgrounds and show striking differences in within-host symbiont densities which correlate positively with Octomom copy number. This is consistent with host genotype selecting different Wolbachia strains, but also with bottlenecks and drift between generations. To distinguish between these possibilities, they perform 2 further experiments. 

A second experiment repeated experiment 1, but this time introgression was into 3 independent lines of the Bolivia and USA Drosophila populations; those that, respectively, exhibited the lowest and highest Wolbachia density and Octomom copy number. In this experiment, growth and Octomom copy number were measured across the 3 lines, for each population, after 1, 13 and 25 generations. Although there were little differences between replicates at generation 1, there were differences at generations 13 and 25 among the replicates of both the Bolivia and USA lines. These results are indicative of parasite control, or drift being responsible for within-host growth rather than host factors. 

A third experiment tested whether Wolbachia density and copy number were under host or parasite control. This was done, again using the USA and Bolivia lines, but this time those from the first experiment, several generations following the initial introgression and backcrossing. The newly introgressed lines were again followed for 25 generations. At generation 1, Wolbachia phenotypes resembled those of the donor parasite population and not the recipient host population indicating a possible maternal effect, but a lack of host control over the parasite. Furthermore, Wolbachia densities and Octomom number differed among replicate lines through time for Bolivia populations and from the donor parasite lines for both populations. These differences among replicate lines that share both host and parasite origins suggest that drift and/or maternal effects are responsible for within-host Wolbachia density and Octomom number. 

These findings indicate that drift appears to play a role in shaping Wolbachia evolution in this system. Nevertheless, completely ruling out the role of the host or parasite in controlling densities will require further study. The findings of Bénard and coworkers (2021) should stimulate future work on the contribution of drift to the evolution of vertically transmitted parasites.

References

Alizon S, Hurford A, Mideo N, Baalen MV (2009) Virulence evolution and the trade-off hypothesis: history, current state of affairs and the future. Journal of Evolutionary Biology, 22, 245–259. https://doi.org/10.1111/j.1420-9101.2008.01658.x

Anderson RM, May RM (1982) Coevolution of hosts and parasites. Parasitology, 85, 411–426. https://doi.org/10.1017/S0031182000055360

Bénard A, Henri H, Noûs C, Vavre F, Kremer N (2021) Wolbachia load variation in Drosophila is more likely caused by drift than by host genetic factors. bioRxiv, 2020.11.29.402545, ver. 4  recommended and peer-reviewed by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2020.11.29.402545

Chrostek E, Marialva MSP, Esteves SS, Weinert LA, Martinez J, Jiggins FM, Teixeira L (2013) Wolbachia Variants Induce Differential Protection to Viruses in Drosophila melanogaster: A Phenotypic and Phylogenomic Analysis. PLOS Genetics, 9, e1003896. https://doi.org/10.1371/journal.pgen.1003896

Chrostek E, Teixeira L (2015) Mutualism Breakdown by Amplification of Wolbachia Genes. PLOS Biology, 13, e1002065. https://doi.org/10.1371/journal.pbio.1002065

Mathé‐Hubert H, Kaech H, Hertaeg C, Jaenike J, Vorburger C (2019) Nonrandom associations of maternally transmitted symbionts in insects: The roles of drift versus biased cotransmission and selection. Molecular Ecology, 28, 5330–5346. https://doi.org/10.1111/mec.15206

Mira A, Moran NA (2002) Estimating Population Size and Transmission Bottlenecks in Maternally Transmitted Endosymbiotic Bacteria. Microbial Ecology, 44, 137–143. https://doi.org/10.1007/s00248-002-0012-9

Vale PF, Wilson AJ, Best A, Boots M, Little TJ (2011) Epidemiological, Evolutionary, and Coevolutionary Implications of Context-Dependent Parasitism. The American Naturalist, 177, 510–521. https://doi.org/10.1086/659002

Wale N, Jones MJ, Sim DG, Read AF, King AA (2019) The contribution of host cell-directed vs. parasite-directed immunity to the disease and dynamics of malaria infections. Proceedings of the National Academy of Sciences, 116, 22386–22392. https://doi.org/10.1073/pnas.1908147116

The advent of genome-wide association studies (GWAS) has been a great promise for our understanding of the connection between genotype and phenotype. Today, the NHGRI-EBI GWAS catalog contains 251,401 associations from 4,961 studies (1). This wealth of studies has also generated interest to use the summary statistics beyond the few top hits in order to make predictions for individuals without known phenotype, e.g. to predict polygenic risk scores or to study polygenic selection by comparing different groups. For instance, polygenic selection acting on the most studied polygenic trait, height, has been subject to multiple studies during the past decade (e.g. 2–6). They detected north-south gradients in Europe which were consistent with expectations. However, their GWAS summary statistics were based on the GIANT consortium data set, a meta-analysis of GWAS conducted in different European cohorts (7,8). The availability of large data sets with less stratification such as the UK Biobank (9) has led to a re-evaluation of those results. The nature of the GIANT consortium data set was realized to represent a potential problem for studies of polygenic adaptation which led several of the authors of the original articles to caution against the interpretations of polygenic selection on height (10,11). This was a great example on how the scientific community assessed their own earlier results in a critical way as more data became available. At the same time it left the question whether there is detectable polygenic selection separating populations more open than ever.

Generally, recent years have seen several articles critically assessing the portability of GWAS results and risk score predictions to other populations (12–14). Refoyo-Martínez et al. (15) are now presenting a systematic assessment on the robustness of cross-population signatures of polygenic adaptation in humans. They compiled GWAS results for complex traits which have been studied in more than one cohort and then use allele frequencies from the 1000 Genomes Project data (16) set to detect signals of polygenic score overdispersion. As the source for the allele frequencies is kept the same across all tests, differences between the signals must be caused by the underlying GWAS. The results are concerning as the level of overdispersion largely depends on the choice of GWAS cohort. Cohorts with homogenous ancestries show little to no overdispersion compared to cohorts of mixed ancestries such as meta-analyses. It appears that the meta-analyses fail to fully account for stratification in their data sets.

The authors based most of their analyses on the heavily studied trait height. Additionally, they use educational attainment (measured as the number of school years of an individual) as an example. This choice was due to the potential over- or misinterpretation of results by the media, the general public and by far right hate groups. Such traits are potentially confounded by unaccounted cultural and socio-economic factors. Showing that previous results about polygenic selection on educational attainment are not robust is an important result that needs to be communicated well. This forms a great example for everyone working in human genomics. We need to be aware that our results can sometimes be misinterpreted. And we need to make an effort to write our papers and communicate our results in a way that is honest about the limitations of our research and that prevents the misuse of our results by hate groups.

This article represents an important contribution to the field. It is cruicial to be aware of potential methodological biases and technical artifacts. Future studies of polygenic adaptation need to be cautious with their interpretations of polygenic score overdispersion. A recommendation would be to use GWAS results obtained in homogenous cohorts. But even if different biobank-scale cohorts of homogeneous ancestry are employed, there will always be some remaining risk of unaccounted stratification. These conclusions may seem sobering but they are part of the scientific process. We need additional controls and new, different methods than polygenic score overdispersion for assessing polygenic selection. Last year also saw the presentation of a novel approach using sequence data and GWAS summary statistics to detect directional selection on a polygenic trait (17). This new method appears to be robust to bias stemming from stratification in the GWAS cohort as well as other confounding factors. Such new developments show light at the end of the tunnel for the use of GWAS summary statistics in the study of polygenic adaptation.

References

1. Buniello A, MacArthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Research. 2019 Jan 8;47(D1):D1005–12. doi: https://doi.org/10.1093/nar/gky1120

2. Turchin MC, Chiang CW, Palmer CD, Sankararaman S, Reich D, Hirschhorn JN. Evidence of widespread selection on standing variation in Europe at height-associated SNPs. Nature Genetics. 2012 Sep;44(9):1015–9. doi: https://doi.org/10.1038/ng.2368

3. Berg JJ, Coop G. A Population Genetic Signal of Polygenic Adaptation. PLOS Genetics. 2014 Aug 7;10(8):e1004412. doi: https://doi.org/10.1371/journal.pgen.1004412

4. Robinson MR, Hemani G, Medina-Gomez C, Mezzavilla M, Esko T, Shakhbazov K, et al. Population genetic differentiation of height and body mass index across Europe. Nature Genetics. 2015 Nov;47(11):1357–62. doi: https://doi.org/10.1038/ng.3401

5. Mathieson I, Lazaridis I, Rohland N, Mallick S, Patterson N, Roodenberg SA, et al. Genome-wide patterns of selection in 230 ancient Eurasians. Nature. 2015 Dec;528(7583):499–503. doi: https://doi.org/10.1038/nature16152

6. Racimo F, Berg JJ, Pickrell JK. Detecting polygenic adaptation in admixture graphs. Genetics. 2018. Arp;208(4):1565–1584. doi: https://doi.org/10.1534/genetics.117.300489

7. Lango Allen H, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, et al. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature. 2010 Oct;467(7317):832–8. doi: https://doi.org/10.1038/nature09410

8. Wood AR, Esko T, Yang J, Vedantam S, Pers TH, Gustafsson S, et al. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014 Nov;46(11):1173–86. doi: https://doi.org/10.1038/ng.3097

9. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018 Oct;562(7726):203–9. doi: https://doi.org/10.1038/s41586-018-0579-z

10. Berg JJ, Harpak A, Sinnott-Armstrong N, Joergensen AM, Mostafavi H, Field Y, et al. Reduced signal for polygenic adaptation of height in UK Biobank. eLife. 2019 Mar 21;8:e39725. doi: https://doi.org/10.7554/eLife.39725

11. Sohail M, Maier RM, Ganna A, Bloemendal A, Martin AR, Turchin MC, et al. Polygenic adaptation on height is overestimated due to uncorrected stratification in genome-wide association studies. eLife. 2019 Mar 21;8:e39702. doi: https://doi.org/10.7554/eLife.39702

12. Martin AR, Kanai M, Kamatani Y, Okada Y, Neale BM, Daly MJ. Clinical use of current polygenic risk scores may exacerbate health disparities. Nature Genetics. 2019 Apr;51(4):584–91. doi: https://doi.org/10.1038/s41588-019-0379-x

13. Bitarello BD, Mathieson I. Polygenic Scores for Height in Admixed Populations. G3: Genes, Genomes, Genetics. 2020 Nov 1;10(11):4027–36. doi: https://doi.org/10.1534/g3.120.401658

14. Uricchio LH, Kitano HC, Gusev A, Zaitlen NA. An evolutionary compass for detecting signals of polygenic selection and mutational bias. Evolution Letters. 2019;3(1):69–79. doi: https://doi.org/10.1002/evl3.97

15. Refoyo-Martínez A, Liu S, Jørgensen AM, Jin X, Albrechtsen A, Martin AR, Racimo F. How robust are cross-population signatures of polygenic adaptation in humans? bioRxiv, 2021, 2020.07.13.200030, version 5 peer-reviewed and recommended by Peer community in Evolutionary Biology. doi: https://doi.org/10.1101/2020.07.13.200030

16. Auton A, Abecasis GR, Altshuler DM, Durbin RM, Abecasis GR, Bentley DR, et al. A global reference for human genetic variation. Nature. 2015 Sep 30;526(7571):68–74. doi: https://doi.org/10.1038/nature15393

17. Stern AJ, Speidel L, Zaitlen NA, Nielsen R. Disentangling selection on genetically correlated polygenic traits using whole-genome genealogies. bioRxiv. 2020 May 8;2020.05.07.083402. doi: https://doi.org/10.1101/2020.05.07.083402