Research Article

Genetics and Genomics

Variable prediction accuracy of polygenic scores within an ancestry group

Department of Biological Sciences, Columbia University, United States
Department of Sociology, Princeton University, United States
Office of Population Research, Princeton University, United States
Department of Genetics, Stanford University, United States
Department of Biology, Stanford University, United States
Howard Hughes Medical Institute, Stanford University, United States
Department of Systems Biology, Columbia University, United States

Jan 30, 2020

Open access
Copyright information

Abstract
eLife digest
Introduction
Results
Discussion
Materials and methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

Fields as diverse as human genetics and sociology are increasingly using polygenic scores based on genome-wide association studies (GWAS) for phenotypic prediction. However, recent work has shown that polygenic scores have limited portability across groups of different genetic ancestries, restricting the contexts in which they can be used reliably and potentially creating serious inequities in future clinical applications. Using the UK Biobank data, we demonstrate that even within a single ancestry group (i.e., when there are negligible differences in linkage disequilibrium or in causal alleles frequencies), the prediction accuracy of polygenic scores can depend on characteristics such as the socio-economic status, age or sex of the individuals in which the GWAS and the prediction were conducted, as well as on the GWAS design. Our findings highlight both the complexities of interpreting polygenic scores and underappreciated obstacles to their broad use.

eLife digest

Complex diseases like cancer and heart disease are caused by the interplay of many factors: the variants of genes we inherit, the lifestyles we lead and the environments we inhabit, plus the interaction of all these factors. In fact, almost every trait, even how many years we will spend studying, is influenced both by our environment and our genes.

To identify some of the genetic factors at play, scientists perform analyses known as genome-wide association studies, or GWAS for short. In these studies, the genomes from many different people are scanned to look for genetic differences associated with differences in traits. By summing up all the small genetic differences, so-called “polygenic scores” can be calculated. When there is a large genetic component to a trait, polygenic scores can be useful predictive tools.

But there is a catch: polygenic scores make less accurate predictions for individuals of a different ancestry than those involved in the GWAS, which limits the use of these tools around the world. Mostafavi, Harpak et al. set out to understand if there are other factors in addition to ancestry that could influence the performance of polygenic scores.

Using data from the UK Biobank, an international health resource that pairs genomic data and clinical information, Mostafavi, Harpak et al. examined polygenic scores among individuals that share a single, common ancestry. These polygenic scores were used to predict three traits (blood pressure, body mass index and educational attainment) in individuals and the predictions were then compared to the actual trait values to see how accurate they were. The analysis revealed that even within a group of people with similar ancestry, the accuracy of polygenic scores can vary, depending on characteristics such as the sex, age or socioeconomic status of the individuals.

This analysis emphasises how variable GWAS and their predictive value can be even within seemingly similar population groups. It further highlights both the complexities of interpreting polygenic scores and underappreciated obstacles to their broad use in medical and social sciences.

Introduction

Genome-wide association studies (GWAS) have now been conducted for thousands of human complex traits, revealing that the genetic architecture is almost always highly polygenic, that is that the bulk of the heritable variation is due to thousands of genetic variants, each with tiny marginal effects (Boyle et al., 2017; Bulik-Sullivan et al., 2015). These findings make it difficult to interpret the molecular basis for variation in a trait, but they lend themselves more immediately to another use: phenotypic prediction. Under the assumption that alleles act additively, a 'polygenic score' (PGS) can be created by summing the effects of the alleles carried by an individual; this score can then be used to predict that individual’s phenotype (Henderson, 1984; Meuwissen et al., 2001; Kathiresan et al., 2008; Lynch and Walsh, 1998). For highly heritable traits, such scores already provide informative predictions in some contexts: for example, prediction accuracies are 24.4% for height (using $R^{2}$ as a measure) (Yengo et al., 2018) and up to 13% for educational attainment (using incremental $R^{2}$ ) (Lee et al., 2018).

This genomic approach to phenotypic prediction has been rapidly adopted in three distinct fields. In human genetics, PGS have been shown to help identify individuals that are more likely to be at risk of diseases such as breast cancer and cardiovascular disease (Khera et al., 2018; Inouye et al., 2018; Mavaddat et al., 2019; Khera et al., 2019). Based on these findings, a number of papers have advocated that PGS be adopted in designing clinical studies, and by clinicians as additional risk factors to consider in treating patients (Torkamani et al., 2018; Khera et al., 2018). In human evolutionary genetics, several lines of evidence suggest that adaptation may often take the form of shifts in the optimum of a polygenic phenotype and hence act jointly on the many variants that influence the phenotype (Pritchard and Di Rienzo, 2010; Berg and Coop, 2014; Höllinger et al., 2019; Sella and Barton, 2019). In this context, the goal is to test whether the set of variants that influence a trait are rapidly evolving across populations or over time (Field et al., 2016; Berg et al., 2019; Uricchio et al., 2019; Edge and Coop, 2019; Racimo et al., 2018; Berg and Coop, 2014). Finally, in various disciplines of the social sciences, PGS are increasingly used to distinguish environmental from genetic sources of variability (Conley, 2016), as well as to understand how genetic variation among individuals may cause heterogeneous treatment effects when studying how an environmental influence (e.g., a schooling reform) affects an outcome (such as BMI) (Barcellos et al., 2018; Davies et al., 2018). In all these applications, the premise is that PGS will ‘port’ well across groups—that is that they remain predictive not only in samples very similar to the ones in which the GWAS was conducted, but also in other sets of individuals (henceforth ‘prediction sets’).

As recent papers have highlighted, however, PGS are not as predictive in individuals whose genetic ancestry differs substantially from the ancestry of individuals in the original GWAS (reviewed in Martin et al., 2019). As one illustration, PGS calculated in the UK Biobank predict phenotypes of individuals sampled in the UK Biobank better than those of individuals sampled in the BioBank Japan Project: for instance, the incremental $R^{2}$ for height is approximately 11% in the UK versus 3% in Japan (Martin et al., 2019). Similarly, using PGS based on Europeans and European-Americans, the largest educational attainment GWAS to date ('EA3') reported an incremental $R^{2}$ of 10.6% for European-Americans but only 1.6% for African-Americans (Lee et al., 2018).

To date, such observations have been discussed mainly in terms of population genetic factors that reduce portability (Martin et al., 2017; Kim et al., 2018; Duncan et al., 2018; De La Vega and Bustamante, 2018; Sirugo et al., 2019; Martin et al., 2019). Notably, GWAS does not pinpoint causal variants, but instead implicates a set of possible causal variants that lie in close physical proximity in the genome. The estimated effect of a given SNP depends on the extent of linkage disequilibrium (LD) with the causal sites (Pritchard and Przeworski, 2001; Bulik-Sullivan et al., 2015). LD differences between populations that arose from their distinct demographic and recombination histories will lead to variation in the estimated effect sizes and hence to variable phenotypic prediction accuracies (Rosenberg et al., 2019). Populations will also differ in the allele frequencies of causal variants. This problem is particularly acute for alleles that are rare in the population in which the GWAS was conducted but common in the population in which the trait is being predicted. Such variants are likely to have noisy effect size estimates in the estimation sample or may not be included in the PGS at all, and yet they contribute substantially to heritability in the target population. Furthermore, causal loci or effect sizes may differ among populations, for instance if the effect of an allele depends on the genetic background on which it arises (e.g., Adhikari et al., 2019). For all these reasons, we should expect PGS to be less predictive across ancestries.

In practice, given that most individuals (about 80%) included in current GWAS are of European ancestry (Popejoy and Fullerton, 2016; Martin et al., 2019), PGS are systematically more predictive in European-ancestry individuals than among other people. As a consequence, the clinical applications and scientific understanding to be gained from PGS will predominantly and unfairly benefit a small subset of humanity. A number of papers have therefore highlighted the importance of expanding GWAS efforts to include more diverse ancestries (Martin et al., 2018; Bien et al., 2019; Wojcik et al., 2019; Martin et al., 2019; Sirugo et al., 2019).

Importantly, factors other than ancestry could also impact the accuracy and portability of PGS. For example, the educational attainment of an individual depends not only on their own genotype, but on the genotypes of their parents, due to nurturing effects (Kong et al., 2018), and of their peers, due to social genetic effects (Domingue et al., 2018), and of course on non-genetic factors. Also, traits such as height and educational attainment show strong patterns of assortative mating, which can distort effect size estimates in GWAS (Domingue et al., 2014; Robinson et al., 2017; Ruby et al., 2018). To what extent these effects remain the same across cultures and environments is unknown, but if they differ, so will the prediction accuracy. More generally, while we still know little about genotype-environment interactions (GxE) in humans, they are well-documented in other species—notably in experimental settings—and would further reduce the portability of PGS across environments (Gibson, 2008; Tropf et al., 2017; Mills and Rahal, 2019; Lynch and Walsh, 1998). In addition, the extent of environmental variability could differ between GWAS and prediction groups, which would change the proportion of the variance in the trait explained by a PGS (i.e., the prediction accuracy). PGS for some traits may also include a component of environmental or cultural confounding with population structure (Sohail et al., 2019; Haworth et al., 2019; Lawson et al., 2020; Kerminen et al., 2018; Berg et al., 2019); this source of confounding can increase or decrease prediction accuracy, depending on the structure in the prediction samples.

Given these considerations, it is important to ask to what extent PGS are portable among groups within the same ancestry. To explore this question, we stratified the subset of UK Biobank samples designated as ‘White British’ (WB) according to some of the standard sample characteristics of GWAS studies: the ages of the individuals, their sex, and socio-economic status. We chose to focus on these particular characteristics because they vary among GWAS samples depending on sample ascertainment procedures. Furthermore, these characteristics have been shown to influence heritability for some traits in a study of a subset of the UK Biobank (Ge et al., 2017), raising the possibility that these choices also influence prediction accuracy. Indeed, for three example traits, we show that there exist major differences in the prediction accuracy of the PGS among these groups, even though they share highly similar genetic ancestries. We further demonstrate for a variety of traits that prediction accuracy differs markedly depending on whether the GWAS is conducted in unrelated individuals or in pairs of siblings, even when controlling for the precision of the estimates. This finding is again unexpected under standard GWAS assumptions; it underscores the importance of genetic effects that are included in estimates from some study designs and not others and highlights underappreciated challenges with GWAS-based phenotypic prediction.

At present, it is difficult to determine the reasons why we see such variable prediction accuracy across these strata and study designs. Contributing factors probably include indirect genetic effects from relatives, assortative mating, varying levels of genetic and environmental variance, GxE interaction effects and perhaps undetected confounding. Nonetheless, our results make clear that the prediction accuracy of PGS can be affected in unpredictable ways by known—and presumably unknown—factors in addition to genetic ancestry.

Results

Sample characteristics of the GWAS and prediction set can influence prediction accuracy even within a single ancestry

We examined how PGS for a few example traits port across samples that are of similar genetic ancestry but differ in terms of some common study characteristics, such as the male:female ratio (henceforth ‘sex ratio’), age distribution, or socio-economic status (SES). To this end, we limited our analysis to the largest subset of individuals in the UKB with a relatively homogeneous ancestry: 337,536 unrelated individuals that were characterized by the UKB, based on self-reported ethnicities as well as genetic analysis, as ‘White British’ (WB) (Bycroft et al., 2018). In all analyses, we further adjusted for the first 20 principal components of the genotype data, to account for population structure within this set of individuals (Materials and methods).

In all analyses, we randomly selected a subset of individuals to be the prediction set; we then conducted GWAS using the remaining individuals and built a PGS model by LD-based clumping of the associations (Materials and methods). To examine the reliability of the prediction, we considered the incremental $R^{2}$ , that is the $R^{2}$ increment obtained when adding the PGS to a model with other covariates (referred to as 'prediction accuracy' henceforth). Whether this measure is appropriate depends on how PGS are to be used; it is not always the most obvious choice in human genetics, where the goal is often to identify individuals at high risk of developing a particular disease (i.e., in the tail of the polygenic score distribution). Nonetheless, because it has been widely reported in discussions of portability across genetic ancestries (e.g., Lee et al., 2018; Martin et al., 2019), we also used it here; later, we also present some results on binary traits using incremental area under the receiver operator curve (AUC).

As a first case, we considered the prediction accuracy of a PGS for diastolic blood pressure in prediction sets stratified by sex, motivated by reports that variation in this trait may arise for somewhat distinct reasons in the two sexes (Reckelhoff, 2001; Zhou et al., 2017). We randomly selected males and females as prediction sets (20K individuals each), and used a subset of the rest of the individuals for GWAS, matching the numbers of females and males in the GWAS set (total sample size 122,774); we refer to this mixed set, somewhat loosely, as the 'diverse GWAS.' Adjusting for mean sex effects and medication use (see Materials and methods), the prediction accuracy is about 1.15-fold higher for females than for males (Mann-Whitney $p = 1.1 \cdot 1 0^{- 5}$ ; Figure 1A). Thus, despite equal representation of males and females in the GWAS set, the prediction accuracy varies depending on the sex ratio of prediction samples. To examine this further, we repeated the same analysis but performed the GWAS in only one sex (which we refer to as 'stratified GWAS' using the same sample size as in the diverse GWAS). [Note that the diverse GWAS sample is not a merge of the stratified GWAS samples but a mixed-sex sample of equal sample size to that used in the women-only and the men-only GWAS, to allow for direct comparison between GWASs. Results for the merged GWAS (with a much larger sample size) are presented in Appendix 1—figure 1A.] When the GWAS is conducted only in females, the prediction accuracy is about 1.35-fold higher for females than for males; in turn, when GWAS was done in only males, the prediction accuracy in both sexes is similar, as well as somewhat decreased (Figure 1A).

Figure 1

Download asset Open asset

Variable prediction accuracy of polygenic scores within an ancestry group.

Shown are incremental $R^{2}$ values (i.e., the increment in $R^{2}$ obtained by adding a polygenic score predictor to a model with covariates alone) in different prediction sets. Each box and whiskers plot is computed based on 20 iterations of resampling GWAS and prediction sets. Thick horizontal lines denote the medians. The polygenic scores were estimated in samples of unrelated WB individuals. Phenotypes were then predicted in distinct samples of unrelated WB individuals, stratified by sex (A), age (B) or Townsend deprivation index, a measure of SES (C). In red and green cases, polygenic scores are based on a GWAS in a sample limited to one sex, age or SES group (a 'stratum'). In blue, polygenic scores are based on a GWAS in a diverse sample matching the number of individuals in each stratum. GWAS samples sizes are: 122,774 for all three diastolic blood pressure GWAS samples, 72,328 for all three BMI GWAS samples, 73,280 for years of schooling GWAS in the diverse sample and 73,283 for GWAS in the low SES and high SES samples.

We then considered two other cases, evaluating prediction accuracy in groups stratified by age for BMI—since the UK Biobank participants were enrolled within about a five-year span, differences in age could in principle also be reflective of cohort effects—and by adult SES for years of schooling, using the Townsend deprivation index as a measure; our choices were motivated by prior evidence suggesting that these characteristics of the GWAS influence estimates of SNP-heritability (Branigan et al., 2013; Conley et al., 2015; Belsky et al., 2018; Elks et al., 2012; Ge et al., 2017). We withheld a random set of 10K individuals in each quartile of age and SES for prediction and performed GWAS using a subset of the remaining individuals, matching the sample sizes across quartiles in the GWAS set (total sample sizes of 72,328 and 73,280 for BMI and years of schooling GWAS, respectively). Similar to our observation for diastolic blood pressure, the prediction accuracy varies across prediction sets: it is 1.4-fold higher for BMI in the youngest quartile compared to the oldest (Mann-Whitney $p = 1.1 \cdot 1 0^{- 5}$ ; Figure 1B), and 2-fold higher for years of schooling in the lowest SES quartile compared to the highest (Mann-Whitney $p = 2.9 \cdot 1 0^{- 6}$ ; Figure 1C). Furthermore, the differences across groups are again sensitive to the choice of the GWAS set: the differences are marked when GWAS is restricted to the youngest quartile for BMI and the lowest SES quartile for years of schooling, but diminished when the GWAS is performed in the oldest and the highest SES quartiles for BMI and years of schooling, respectively (Figure 1B, C). These results remained qualitatively unchanged when we used $R^{2}$ instead of incremental $R^{2}$ to measure prediction accuracy (Appendix 1—figure 2).

In these analyses, we used a p-value threshold of 10^-4 for inclusion of a SNP in the PGS. The choice of how stringent to make the GWAS p-value threshold is important but somewhat arbitrary, with approaches ranging from requiring genome-wide significance to including all SNPs (Weedon et al., 2008; Pharoah et al., 2008; Euesden et al., 2015; Vilhjálmsson et al., 2015; Ware et al., 2017; Mostafavi et al., 2017; Speidel et al., 2019). Often, this threshold is chosen to maximize prediction accuracy in an independent validation set. When the goal is to compare prediction performance across different groups, there is no obvious optimal choice of the p-value threshold. [The optimal p-value in this context will differ across studies, as it depends not only on the genetic architecture and heritability of the trait, but also on the GWAS sample size, that is power (Dudbridge, 2013).] As we show, however, the qualitative trends reported in Figure 1 do not depend on the p-value threshold choice (Appendix 1—figure 3); moreover, the qualitative trends remain when LDpred is used (with a prior probability of 1 on loci being causal; Vilhjálmsson et al., 2015) instead of pruning approaches (Appendix 1—figure 3).

These results pertain to three exemplar traits and do not speak to the prevalence of this phenomenon. Nonetheless, they demonstrate that the prediction accuracy of a polygenic score can vary markedly depending on sample characteristics of both the original GWAS and the prediction set, even within a single ancestry, and that this variation in prediction accuracy can be substantial—on the same order as reported for different continental ancestries within the UK Biobank (Martin et al., 2019). As one example, the prediction accuracy in East Asian samples, averaged across a number of traits, is about half of that in European samples when GWAS was European-based; when the GWAS is done in the lowest SES group for years of schooling, prediction accuracy in the highest SES group is less than half of that in the lowest SES (Figure 1C). Moreover, whereas for these traits, we had prior information about which characteristics may be relevant, other aspects that vary across sets of individuals are undoubtedly important as well (e.g., smoking behavior and diet may modify genetic effects on lipid traits; Bentley et al., 2019; Telkar et al., 2019), and for other traits of interest, much less may be known a priori.

Possible explanations for the variable prediction accuracy

Our goal in this paper is to highlight that prediction accuracies can vary across groups of highly similar ancestry, rather than to investigate the likely causes for any particular phenotype. Nonetheless, we provide some observations that may cast light on these results. We first note that in these three examples, the prediction accuracies track SNP heritability differences across strata (Figure 2A,B,C). This relationship should be expected, given that the estimation noise decreases with heritability (Appendix 1), and potentially underlies the observation that prediction accuracies using the diverse GWAS sample are often intermediate between those obtained from stratified GWAS samples of equal sample size (Figure 1).

Figure 2

Download asset Open asset

Differences in environmental variance alone do not explain the variable prediction accuracy.

(**A,B,C**) The x-axes show heritability estimates (± SE) based on LD score regression in each set. The y-axes show incremental $R^{2}$ values obtained using the procedure described in Figure 1, with GWAS performed in a pooled sample of all strata and testing in stratified prediction sets (see Materials and methods); points and bars show mean and central 80% range computed based on 20 iterations of resampling GWAS and prediction sets. ‘Q’ denotes quartile of age and SES in (**B,E**) and (**C,F**), respectively. (**D,E,F**) The x-axes show phenotypic variance estimates (± SE) across strata after adjusting for covariates (sex, age and 20 PCs). If the heritability differences across strata are due to differences in environmental variance alone, with genetic variance constant, then heritability should be inversely proportional to phenotypic variance. The best-fitting model for this inverse proportionality (dashed line, simple linear regression) provides a poor fit to the observations.

Perhaps the simplest explanation for these findings would be that heritabilities, and hence prediction accuracies, vary only because of differences in the extent of environmental variance across strata, while the genetic variance is the same. We can test this hypothesis by examining whether the heritability decreases with increasing phenotypic variance (more precisely whether it is inversely proportional to it), as expected if the genetic variance is fixed across strata. What we find instead is that the estimated SNP heritabilities for all three traits increase or remain the same with increasing phenotypic variance (Figure 2D,E,F). Thus, for these traits at least, the variable prediction accuracy is not simply the result of differences in the extent of environmental heterogeneity across strata.

Another possibility is that there is an interaction between genetic effects and sample characteristics, for instance that different sets of genetic variants contribute to blood pressure levels in males and females or to BMI across different stages of life. [Although such interactions could in some contexts be thought of as reflecting GxE, we use the term ‘sample characteristic’ rather than ‘environment’, as environment has different meaning across disciplines, referring in some contexts only to factors that are exogenous to genetics. Viewed in this lens, SES in adulthood cannot be interpreted as exogenous, because it is in part determined by educational achievement, which is itself influenced by genetic factors, and similarly it is questionable whether age or sex are environments.] This explanation is not supported by bivariate LD score regression, which indicates that the genetic correlations across strata are close to 1 (Appendix 1—table 2; Materials and methods). Yet when we re-estimate individual SNP effects in the prediction sets for SNPs ascertained in the original GWAS, the estimated effects of trait-increasing alleles are larger in the groups with higher prediction accuracy (Appendix 1—figure 4; Materials and methods).

One simple model that could reconcile these findings is if effect sizes are highly correlated across the groups, but systematically larger in those groups with higher prediction accuracy. This explanation is reminiscent of the ‘amplification’ model of genetic influences on cognition during development (Briley and Tucker-Drob, 2013).

Other factors complicate interpretation, however, and may also contribute to our observations. In particular, for the case of years of schooling, conditioning on adult SES induces a form of range restriction, which could contribute to variable prediction accuracy across strata. We note, however, that we see highly variable prediction accuracies across SES strata even when the GWAS is conducted in a diverse sample (i.e., including individuals from all strata) (Figure 1C); in that regard, our approach mimics what happens in practice when polygenic scores are used to predict phenotypes in a sample with a smaller range of SES (e.g., Rimfeld et al., 2018). More generally, although this type of range restriction is artificially amplified in our example, SES differences may often be a problem for GWAS in which the sample is not representative of the population; for instance, the most recent major GWAS of educational attainment (Lee et al., 2018) included numerous medical data sets and the 23andMe data set, which are not representative of the national population.

Another potentially important factor is that the adjustment for PCs may not be a sufficient control for the different ways in which population structure can confound GWAS results (Vilhjálmsson and Nordborg, 2013), leading to variable prediction accuracy across strata if they differ in their population structure. To examine this possibility, we repeated the analysis in Figure 1 but using a linear mixed model (LMM) approach (including PCs among other covariates; see Materials and methods), and obtained qualitatively similar results (Appendix 1—figure 5). Although not a perfect fix (Listgarten et al., 2013; Mathieson and McVean, 2013), the fact that we obtain similar results using PCs and LMM suggests that confounding due to population stratification in the UK Biobank alone does not explain the variable prediction accuracies across strata.

Obstacles to portability explored through a comparison of standard and family-based GWAS

Beyond sample characteristics such as age or sex, a number of other factors may shape the portability of scores across groups of similar ancestry. Standard GWAS is done in samples of individuals that deliberately exclude close relatives; as implemented, it detects direct effects of the genetic variants, but also any indirect genetic effects of parents, siblings, or peers, effects of assortative mating among parents, and potentially environmental differences associated with fine-scale population structure (Young et al., 2018; Trejo and Benjamin, 2019; Kong et al., 2018; Lee et al., 2018; Berg et al., 2019). Given that many of these effects are likely to be culturally mediated (Stulp et al., 2017; Selzam et al., 2019), it seems plausible that they may vary within as well as across groups of individuals with different ancestries. If culturally-contingent effects contribute to GWAS estimates (and hence to PGS), they may lead to differences in the prediction accuracy in samples unlike the original GWAS.

To demonstrate that these considerations are not just hypothetical, we compared the prediction accuracy when the PGS is trained on ‘unrelated’ individuals such as those used in a standard GWAS to one obtained from a sibling-based (or ‘sib-based’) GWAS (Materials and methods). In the latter, genotype differences between sibs, a result of random Mendelian segregation in the parents, are tested for association with the phenotypic differences between them. Because the tests depend on phenotypic differences between siblings who, of course, have the same parents, these tests are conditioned on the parental genotypes and hence exclude many of the indirect effects signals that may be picked up in standard GWAS (Appendix 1). Differences between standard and sib-based GWAS are thus informative about the presence of factors other than direct genetic effects (Wood et al., 2014; Trejo and Benjamin, 2019; Lee et al., 2018; Berg et al., 2019; Selzam et al., 2019).

A challenge in this comparison is that the UKB contains only ~22K sibling pairs, ~19K of whom are labeled as ‘White British’ (WB). The siblings are similar to the unrelated individuals in terms of ages, SES distributions and genetic ancestries (Appendix 1—figures 6 and 7) but include a higher proportion of females; this difference is unlikely to influence our analyses (see below). While a large number, 19K pairs is still too few to have adequate power to discover trait-associated SNPs, when compared to a standard GWAS using the much larger sample of unrelated WB individuals (~340K).

To increase power and enable a direct comparison between the two designs, we split the SNP ascertainment and effect estimation steps as follows (Figure 3A): we identified SNPs using a standard GWAS with a large sample size (median ~270K across the traits considered) (see Materials and methods). We then estimated the effect of each significant SNP using (i) a sib-based association test and (ii) a standard association test. We chose the size of the estimation set in (ii) such that the median standard error of effect estimates in (i) and (ii) is approximately equal. We then compared the prediction accuracy of the two PGS obtained in this way (‘standard PGS’ and ‘sib-based PGS’) in an independent prediction set of unrelated individuals; as we show in Appendix 1, our approach leads to highly similar prediction accuracies of the two approaches under a model with direct effects only (see Materials and methods for details). A further advantage is that the two scores are compared for the same set of SNPs, such that LD patterns and allele frequency differences do not come into play.

Figure 3

Download asset Open asset

Comparison of prediction accuracy of standard and sib-based polygenic scores.

(A) After ascertaining SNPs in a large sample of unrelated individuals, we estimated the effects of these SNPs with a standard regression using unrelated individuals and, independently, using sib-regression. We then used the polygenic scores for prediction in a third sample of unrelated individuals. We chose the sample size of the standard PGS estimation set such that median effect estimate SEs are equal in the two designs, thereby ensuring equal prediction accuracy under a vanilla model with no indirect effects or assortative mating. Numbers in parentheses are median sample size in each set across 20 traits (see Materials and methods and Appendix 1—table 1 for the definition of each trait, and Appendix 1—table 3 for sample sizes for each trait). (B) Ratio of prediction accuracy in the two designs across 20 traits. For each trait, we performed 10 resampling iterations of unrelated individuals into three sets for discovery, estimation and prediction (small points). Large points show median values. (**C-F**) We repeated this procedure with different discovery-set p-value thresholds for including a SNP in the polygenic score. The higher the p-value threshold is, the more SNPs are included. For each p-value threshold, points show 10 iterations as described and large points show median values. Shown are a subset of traits, with traits appearing in (B) but not shown here presented in Appendix 1—figure 12.

We applied the approach to 20 traits, focusing on traits with relatively high heritability estimates as well as social and behavioral traits that have been the focus of recent attention in social sciences. For the majority of the traits, such as diastolic blood pressure, BMI, and hair color, the prediction accuracies of standard and sib-based PGS were similar (Figure 3B), as expected under standard GWAS assumptions and as observed for traits simulated under these assumptions (Appendix 1—figure 8). However, for height and for a range of social and behavioral traits, such as years of schooling, pack years of smoking and household income, the prediction accuracy of the sib-based PGS was substantially lower than that of the standard PGS (Figure 3B). [We caution that, because the first step of our study design is to identify SNPs that are associated with the trait in a large set of unrelated individuals and we subsequently match the sampling variances of sib- and standard GWAS, rather than identify distinct sets of SNPs separately in the two designs, the ratio of prediction accuracies that we obtain cannot be directly compared to those reported in other studies.]

A number of factors could contribute to the differences between prediction accuracies for PGS based on sibs versus unrelated individuals, including confounding effects of population stratification, indirect genetic effects from parents and assortative mating. The relative importance of each factor will vary across traits (Rosenberg et al., 2019; Kong et al., 2018; Haworth et al., 2019; Ruby et al., 2018; Selzam et al., 2019). For educational attainment, this gap is likely to reflect at least in part the documented contribution of indirect genetic effects to the standard PGS (Lee et al., 2018; Kong et al., 2018; Young et al., 2018). We show in Appendix 1 that in the presence of indirect genetic effects mediated through parents, standard PGS outperforms sib-based PGS unless direct and indirect effects are strongly anticorrelated (Appendix 1—figure 9), which seems unlikely to be the case for years of schooling. The difference in the performance of sib-based and standard PGS observed for other social and behavioral outcomes, such as household income and age at first sexual intercourse (Figure 3B), may reflect a similar phenomenon. An additional contribution to divergent prediction accuracies could come from indirect effects among siblings, which would also contribute differentially to standard and sibling-based PGS. For height, there may be an important contribution of assortative mating to the difference in prediction accuracies (Wood et al., 2014; Robinson et al., 2017; Lee et al., 2018). In Appendix 1, we show that under a simple model of positive assortative mating, the prediction accuracy based on a standard PGS is higher than that of a sib-based PGS (Appendix 1—figure 10). We further confirmed that the difference in the sex ratio of the siblings and unrelated individuals, mentioned earlier, has a negligible effect on these differences, though it may underlie the slightly lower prediction accuracy of the standard PGS for pulse rate (Appendix 1—figure 11).

The lower prediction accuracies for PGS based on sib-based GWAS indicate that complications such as assortative mating or indirect effects contribute to the standard GWAS estimates. In the absence of these complications, we ensure that prediction accuracies are comparable by matching the sampling errors of the two approaches (Figure 3A). In the presence of these complications, the magnitude of the ratio of prediction accuracies should reflect the strength of assortative mating, the relative contribution of indirect genetic effects compared to direct effects, and so forth. However, interpreting the magnitude of the deviation from 1 is far from straightforward: as we show in Appendix 1, the relative difference in prediction accuracies between the two approaches stems in part from the noise-to-signal ratio for the effect estimates in sib-based versus standard GWAS (Appendix 1, Appendix 1—figures 9 and 10), and as a result also depends on features of the comparison like the sample sizes used and the PGS model.

Motivated by these considerations, we examined how the prediction accuracy varies when progressively relaxing the GWAS p-value threshold for inclusion of SNPs, that is when including more weakly associated SNPs in the PGS. [In Figure 3B, results are shown for the p-value threshold that maximizes the prediction accuracy of the standard PGS, replicating the practice when comparing populations of different ancestry; Martin et al., 2019.] For hair color and diastolic blood pressure, there is little to no difference in prediction accuracy between the two estimation methods, regardless of the number of SNPs included in the score (Figure 3C,D). In contrast, for height, standard and sib-based PGS perform similarly when based on the most significantly associated SNPs, but standard PGS progressively outperforms sib-based PGS when more SNPs are included (Figure 3E). Similarly, the difference in prediction accuracy between sib-based and standard PGS changes markedly for years of schooling, household income and other social and behavioral traits (Figure 3F and Appendix 1—figure 12). The growing gap in performance with increasing p-value threshold likely reflects a combination of an increasing noise-to-signal ratio for the effect estimates in sib-based versus standard GWAS (see Appendix 1) and changes in the relative importance of direct effects versus other factors such as indirect parental effects and assortative mating.

In summary, the differences between the prediction accuracies of standard and sib-based PGS seen for a number of traits (Figure 3B), notably social and behavioral ones, demonstrate that standard GWAS estimates often include a substantial contribution of factors other than direct effects. In these cases, even if the power to detect direct effects were comparable, standard GWAS would lead to higher prediction accuracy than sib-GWAS. In some contexts that may be a sufficient reason to rely on PGS derived from standard GWAS. However, that gain stems from the inclusion of factors such as indirect effects and assortative mating that are likely to be modulated by SES, environment and culture (e.g., Selzam et al., 2019; Stulp et al., 2017). Thus, the increased prediction accuracy likely comes at a cost of not always porting well across groups, even of the same ancestry, in ways that may be difficult to anticipate.

Discussion

Although the conversation around the portability of PGS has largely focused on genetic ancestries, our results show that prediction accuracy can also differ, in some cases substantially, across groups of similar ancestry—even due to basic study design differences such as age, sex or SES composition. When due only to increased environmental variance, such decreased accuracy may not pose a problem, at least for certain applications. But as we have shown, differences in the degree of environmental variance are not the primary explanation for the patterns we report (Figure 2), and other factors, including differences in the magnitude of genetic effects among groups, indirect effects and assortative mating, also lead to differences in the prediction accuracy of PGS, in ways that may make applications of phenotypic prediction less reliable, even within a single ancestry group. For some traits, there is prior information about which factors are likely to be important, but not always, and even for well-studied traits, it may be difficult to enumerate all the influential factors. As an example, we considered the accuracy of the polygenic score for years of schooling and found that it also varies somewhat depending on whether individuals have no sibling or one sibling in the prediction sets (Materials and methods; Appendix 1—figure 13).

Following the discussion of portability across ancestries, we have focused on incremental $R^{2}$ as a measure of portability. This measure is less directly informative when the goal is to use PGS to reliably identify individuals in the tails of the distribution, that is those at elevated risk of developing a disease—the main application of PGS in human genetics, as distinct from social science or evolutionary biology. Nonetheless, the same concerns raised here are likely to apply. To illustrate that point, we considered binary outcomes of the traits considered in Figure 1, 'hypertension' (defined as diastolic blood pressure > 110 mmHG), 'obesity' (defined as BMI > 35 kg/m²), and 'college completion', and evaluated the prediction accuracy as measured by incremental AUC (Appendix 1—figure 14).The qualitative results are the same as in Figure 1. We also examined how incremental AUC varies by sex for five binary disease traits that we chose because they have relatively high heritability. For three of them, hypothyroidism and two cardiovascular outcomes, prediction accuracy varies depending on both the GWAS and prediction sets (Appendix 1—figure 15).

Thus, for both quantitative and binary traits, the question of the domain over which a PGS applies is not just about LD patterns, allele frequencies or GxG effects but also about the extent of environmental and genetic variance, GxE, as well as the contribution of direct effects versus indirect effects, assortative mating and environmental confounding. An important implication is that differences in prediction accuracies among groups with distinct ancestries cannot be interpreted exclusively or even primarily in terms of population genetic parameters when these groups differ dramatically in their SES (Chetty and Hendren, 2018; Conley, 2010; Nuru-Jeter et al., 2018; Reich, 2017) and other factors that may affect portability—especially when the relative contribution of these factors to GWAS signals remains unknown (Young et al., 2019; Mills and Rahal, 2019). Thus, efforts to conduct GWAS in groups that vary in ancestry and geographic locations will need to be accompanied by a careful examination of variation in portability along other dimensions.

While these results raise the question of how to best construct a PGS, the answer is not obvious, and likely depends on the specific trait and samples. For example, for the three cases shown in Figure 1, considering a fixed GWAS sample size, the highest prediction accuracy is attained with a GWAS sample limited to some stratum (e.g., women for diastolic blood pressure). Yet a much larger merged data set containing the union of strata generates the most predictive PGS (Appendix 1—figure 1). Together, these observations suggest a trade-off between the factors that are shared among strata and lead to increased power with sample size and those that differ across strata and underlie the variable prediction accuracy. In principle then, if influential factors were known, the composition of the GWAS sample could be optimized to yield the highest accuracy in a given prediction set, but how much each stratum should be weighted will depend on a number of factors such as the genetic and environmental variance in each stratum, genetic correlation across strata, and sample sizes. Moreover, factors such as assortative mating and indirect effects are soaked up into the GWAS estimates—and critically also into the SNP heritability estimates. Thus, the choice of a GWAS sample is about more than power; it is implicitly making a choice about all sorts of sample characteristics that may or may not hold true of the prediction set.

In that regard, it is worth noting that while classical twin studies were often constituted to be representative of a reference population (often national in nature) (Polderman et al., 2015; Branigan et al., 2013), the same is not true of most contemporary human genetic datasets, which are skewed towards medical case-control studies, biobanks that are opt-in (and thus tend to include individuals who are wealthier and better educated than the population average) or direct-to-consumer proprietary genetic databases (which are even more skewed along these dimensions) (Lee et al., 2018). For instance, individuals in UK Biobank have higher SES than the rest of the British population (Fry et al., 2017) and are presumably self-selected for a certain level of interest in biomedical research. These factors alone raise challenges as to the broad portability of PGS derived from them. More generally, it seems plausible that individuals included in a GWAS differ from those that, for myriad reasons, do not end up participating (Taylor et al., 2018), in ways that make it difficult to predict the domain over which GWAS-based estimates can be reliably generalized.

One fruitful way forward may be to study data from related individuals, in which it should be possible to decompose the components of the signals identified in GWAS into direct and indirect effects, the degree of assortative mating and the contribution of residual stratification (Zhang et al., 2015; Young et al., 2018; Kong et al., 2018). Not only will this decomposition help us to better interpret the results of GWAS and the resulting PGS, it will make it possible to examine under which circumstances, and for which phenotypes, components port more reliably to other sets of individuals, both unrelated and related. Ultimately, we envisage that in order to be broadly applicable, GWAS-based phenotypic prediction models will need to include not only a PGS but some study characteristics, other social and environmental measures and, perhaps crucially, their interactions.

Materials and methods

UK biobank

The UK Biobank (UKB) is a large study of about half a million United Kingdom residents, recruited between years 2006 to 2010 (Bycroft et al., 2018). In addition to genetic data, hundreds of phenotypes were collected through measurements and questionnaires at assessment centers, and by accessing medical records of the participants.

Inclusion criteria

Request a detailed protocol

In this study, we focused on 408,434 participants who passed quality control (QC) measures provided by UKB; specifically, for whom the reported sex (QC parameter ‘Submitted.Gender’) matched their inferred sex from genotype data (QC parameter ‘Inferred.Gender’); who were not identified as outliers based on heterozygosity and missing rate (QC parameter ‘het.missing.outliers’==0); and did not have an excessive number of relatives in the database (QC parameter ‘excess.relatives’==0). We further selected individuals identified by UKB to be of ‘White British’ (WB) ancestry (QC parameter ‘in.white.British.ancestry.subset’==1), which is a label that refers to those who, when given a set of choices, self-reported to be of ‘White’ and ‘British’ ethnic backgrounds and, in addition, were tightly clustered in a principal component analysis of the genotype data, as detailed in Bycroft et al. (2018). We excluded individuals that had withdrawn from the UK Biobank by the time of the analyses here. For a given trait, we further conditioned on individuals for whom the trait value was reported.

Phenotype data

Request a detailed protocol

We focused on 25 traits, including traits with relatively high heritability estimates as well as social and behavioral traits that have been the focus of recent attention in social sciences (see Appendix 1—table 1 for a complete list of phenotype data used in this work, and their corresponding numeric field codes in the UKB data showcase). We calculated the phenotype ‘years of schooling’ by converting the maximal educational qualification of the participants to years following Okbay et al. (2016) (Appendix 1—table 4). For diastolic blood pressure, pulse rate, and forced vital capacity, we took the average of the first two rounds of measurement taken during the same examination at UKB assessment centers. We adjusted the diastolic blood pressure levels for blood pressure lowering medication following Evangelou et al. (2018) by shifting the values upward by 10 mmHg for individuals taking medication. For hand grip strength, we took the average of the measurements for the two hands. For categorical phenotypes, we assigned integer values to each category (Appendix 1—table 1). For hair color, individuals who reported hair color variable ‘Other’ were excluded from the analyses. We considered binary traits, ‘hypertension’ defined as diastolic blood pressure >110 mmHG, ‘obesity’ defined as BMI >35 kg/m², and ‘college completion’ defined based on attainment of a college or a university degree. Disease outcomes were ascertained using self-reported information and/or using the hospital inpatient main and secondary diagnoses coded according to the International Classification of Diseases (ICD-9 and ICD-10). Hypothyroidism, type 2 diabetes, and rheumatoid arthritis were ascertained based on ICD-10 codes of E03.X, E11.X and M06.X, respectively. Myocardial infarction was ascertained based on ICD-9 codes of 410.9, 411.9, 412.9, or ICD-10 codes of I21.X, I22.X, I23.X, I24.1, I25.2 following Khera et al. (2018), or participants with myocardial infarction outcome data among the UK Biobank’s algorithmically-defined outcomes. We also considered the binary outcome of ever being diagnosed to have had a heart attack, angina or stroke. For a subset of individuals, multiple measurements of a phenotype were provided, corresponding to multiple visits to UKB assessment centers; in those cases, we used the measurements during the first visit.

Genotype data

Request a detailed protocol

UKB participants were genotyped on either of two similar genotyping arrays, UK Biobank Axiom and UK BiLEVE arrays, at a total of ~850K markers. We focused on autosomal bi-allelic SNPs shared between both arrays, and used plink v. 1.90b5 (Chang et al., 2015) to filter SNPs with calling rate >0.95, minor allele frequency >10⁻³, and Hardy-Weinberg equilibrium test p-val >10⁻¹⁰ among the WB samples, resulting in 616,323 SNPs.

GWAS and trait prediction methods

GWAS by sample characteristics

Request a detailed protocol

We focused on a set of 337,488 WB samples that were identified by the UKB to be ‘unrelated’ (sample QC parameter ‘used.in.pca.calculation’==1 as provided by UKB), defined such that no pairs of individuals are inferred to be 3^rd degree relatives or closer. We split the sample into non-overlapping sets of individuals by one of the following factors: age at recruitment (in years), sex, and Townsend deprivation index at recruitment (used as a proxy for socio-economic status or SES, specifically we take the negative of the Townsend deprivation index as a measure of SES). For SES and age, we divided the sample into four sets: Q1 [minimum value, first quartile], Q2 (first quartile, second quartile], Q3 (second quartile, third quartile], and Q4 (third quartile, maximum value]. We randomly selected 10K samples in each SES and age group, and 20K of males and 20K of females as held-out prediction sets, and performed GWAS using the remaining samples, matching sample sizes across groups in the GWAS set. We performed nine GWASs: for years of schooling in SES Q1 and SES Q4 (sample size 73,283 for each), and in a diverse sample with equal number of individuals from all four groups (sample size 73,280); for body mass index (BMI) in Q1, Q4, and in a diverse sample with equal number of individuals from all four groups (sample size 72,328 for each); and for diastolic blood pressure in males, females, and in a diverse sample with equal number of males and females (sample size 122,774 for each). We performed all GWASs using plink v. 2.0 (with the flag --linear), adjusting for sex, age (at recruitment) and first 20 PCs as covariates. PCs are principal components of the genotype data, as provided by UKB, calculated using the entire cohort (not just WB individuals). For a subset of cases (where GWAS was performed in samples restricted by characteristics described above), we additionally performed association tests using a linear mixed model (LMM) as implemented in BOLT-LMM v. 2.3.2 (Loh et al., 2015), using LD scores computed from 1000 Genomes European-ancestry samples, with sex, age and first 20 PCs as covariates. The GWAS summary statistics were used to construct PGS for the samples in the prediction sets.

To better understand the performance of PGS across the strata (see ‘Possible explanations for the variable prediction accuracy’), we estimated the mean effect sizes of significant SNPs in each of the strata. To avoid overfitting, we first performed an association test in the pooled sample of all strata excluding individuals in the prediction sets and matching the number of individuals per stratum; sample size 293,132 for years of schooling, 272,456 for BMI, and 245,548 for diastolic blood pressure. Then for significantly associated SNPs (LD pruned as described in ‘Polygenic score construction and trait prediction’), we re-estimated the effect sizes in each of the strata in the prediction sets (see Appendix 1—figure 4). We also used these pooled GWASs to explore the relationship between prediction accuracy and SNP heritability (as shown in Figure 2) and with GWAS sample size (Appendix 1—figure 1). We performed 20 iterations of all above steps.

In addition to above examples, we explored the prediction accuracy for years of schooling when GWAS and prediction sets are stratified based the participants’ number of full siblings. Specifically, we performed GWAS using individuals who had exactly one sibling (sample size 90,417), and evaluated prediction in two independent samples of individuals who reported having no siblings or having one sibling (sample size 20K for each) (see Appendix 1—figure 13).

We also considered five binary disease outcomes stratified by sex. Specifically, we performed GWAS in equally sized samples of males and females for hypothyroidism (sample size 135,526), type 2 diabetes (sample size 136,061), rheumatoid arthritis (sample size 136,039), myocardial infarction (sample size 136,061) and having been diagnosed with a heart attack or angina or stroke (sample size 135,833), leaving out 20K samples of males and females for prediction (see Appendix 1—figures 14 and 15). For these traits we used a logistic regression model for GWAS (using plink v. 2.0 with the flag --logistic). An important caveat to analyses of disease outcomes recorded during multiple follow-ups is that for ‘age’, we could only consider the age at recruitment in the GWAS; that approach is not ideal, considering that a fraction of individuals died during the course of the study (about 20K individuals in the full cohort).

Standard versus sibling-based polygenic score

Request a detailed protocol

We used the genetic relatedness information provided by UKB to infer sibling pairs among the WB samples. Following Bycroft et al. (2018), we marked pairs with $\frac{1}{2^{5 / 2}} < ϕ < \frac{1}{2^{3 / 2}}$ and IBS0 > 0.0012 as siblings, where $ϕ$ is the estimated kinship coefficient and IBS0 is the fraction of loci at which individuals share no alleles. By this approach, we identified 19,329 sibling pairs including 35,634 individuals across 17,328 families. For a given trait, we included pairs with the property that trait values for both individuals were reported. We then formed two sets of individuals: 'Siblings' set, including the sibling pairs randomly sampled to include only one pair per family, and an 'Unrelateds' set, including the unrelated individuals identified by the UKB (see section 'GWAS by sample characteristics' above), but excluding the Siblings and 6,911 individuals that were related to the Siblings (3^rd degree or closer).

We focused on 20 quantitative traits (see Figure 3B for the list of traits considered in this analysis) and a number of simulated traits (see below). For each trait, we first downsampled the Unrelateds set to a sample size $n^{*}$ such that the median standard error of effect estimates roughly matched the median standard error in the sibling-based regression (see 'Estimating $n^{*}$ ' below). We then divided the Unrelateds set into three non-overlapping sets: after sampling $n^{*}$ individuals (Unrelateds- $n^{*}$ set), we randomly split the rest of the Unrelateds set into an Unrelateds-prediction set (10% of the samples) to be used as a sample for trait prediction ('prediction set'), and an Unrelateds-discovery set (90% of the samples) to be used for the discovery of trait associated variants (see Appendix 1—figure 3 for sample sizes in each set). For each trait, we performed standard GWAS in the Unrelateds-discovery set, and ascertained SNPs by thresholding on association p-values. We then estimated the effect sizes for these ascertained SNPs in two ways: by a sibling-based association test in the Siblings set (using plink v. 1.90b5’s QFAM procedure with the flag --qfam), and by a standard association test in the Unrelateds- $n^{*}$ set (using plink v. 2.0). Subsequently, for each set of ascertained SNPs in the Unrelateds-discovery set, two PGS were constructed for the samples in the Unrelateds-prediction set (see Figure 3A for overview of the pipeline). We performed 10 iterations of the above sampling, ascertainment and estimation steps, except for simulated traits where we performed 30 iterations.

Estimating $n^{*}$

Request a detailed protocol

In order to compare the performance of sibling-based and standard GWAS designs, we wanted to match both analyses to have similar prediction accuracy under a vanilla model of no assortative mating, population structure stratification or indirect effects. In Appendix 1, we show that this could be achieved by matching median effect estimate standard errors. For each trait, we therefore calculated $n^{*}$ , the sample size of a standard GWAS that yields roughly equal standard errors in the standard and sibling-based regressions. Specifically, for each trait, we first performed sibling-based GWAS in the Siblings using plink’s QFAM procedure (with the flag --qfam mperm=100000 emp-se). We then randomly sampled a range of sample sizes from the set of Unrelateds, from 5K to 20K in 1K increments. Following Wood et al. (2014), for each sample size, we performed a standard GWAS, and investigated the linear relationship between the square root of the sample size and the inverse of the median standard error of the effect size estimates. We then used this linear relationship to estimate the sample size of a standard GWAS that corresponds to the inverse of the median standard error of the effect sizes estimate in the sibling-based GWAS.

All standard association tests were performed using plink v. 2.0 (with the flag --linear), adjusting for sex, age and first 20 PCs as covariates. For sibling-based association tests we first residualized the phenotypic values on age and sex, and then regressed the sibling differences in residuals on sibling genotypic differences using plink’s QFAM procedure as described above.

We also considered a version of the analysis described above, in which we first residualized the phenotypes on covariates in the pooled sample of all WB individuals, and then ran the pipeline on the residuals without further adjustment for covariates in the GWAS or prediction evaluation. As shown in Appendix 1—figure 16, this approach produced results that are qualitatively the same to what we present in Figure 3.

Simulated traits

Request a detailed protocol

We wanted to check that given the study design described above, sibling-based and standard PGS perform similarly with respect to trait prediction, under the vanilla model of no population stratification, assortative mating or indirect genetic effects (Figure 3). To this end, we simulated traits with heritability $h^{2} =$ 0.1 or 0.5 and either 10K or 100K causal SNPs. For each set of parameters, we simulated three replicates giving a total of 12 simulated traits.

We randomly selected the causal SNPs from a set of 10,879,183 imputed SNPs, considering that most causal variants are plausibly not directly genotyped on SNP arrays. We used a set of SNPs that passed quality control procedures by the Neale lab (http://www.nealelab.is/uk-biobank), namely autosomal SNPs, imputed using the haplotype reference consortium (HRC) panel, which have INFO score > 0.8 and have minor allele frequency > 10⁻⁴; we further limited the SNP set to ones that were bi-allelic in the WB sample. As in Martin et al. (2017), we randomly assigned effect sizes to these causal SNPs as $β ~ N (0, \frac{h^{2}}{m})$ , and zero for non-causal SNPs. We then calculated genetic component of the trait, $g$ , for all WB samples under an additive model by summing the allelic counts weighted by their effect sizes using plink (with the flag --score). Allelic counts were determined by converting imputation dosages to genotype calls with no hard calling threshold. We also assigned environmental contributions as $ε ~ N (0,1 - h^{2})$ , and then constructed the PGS for each individual,

g = \sum_{i = 1}^{m} β_{i} X_{i},

where $X_{i}$ is the number of minor alleles at SNP $i$ carried by the individual, and the trait value for the individual is calculated as the sum of genetic and environmental contributions:

y = \sqrt{h^{2}} (\frac{g - \bar{g}}{σ_{g}}) + \sqrt{1 - h^{2}} (\frac{ε - \bar{ε}}{σ_{ε}})

where bars represent averages, $σ_{g}$ is the standard deviation of PGS across individuals and $σ_{ε}$ is the standard deviation of environmental contributions across individuals. These simulated traits were then analyzed using the same pipelines as the other traits (e.g., adjusting for covariates etc.). Importantly, SNP discovery and effect size estimations in GWAS were performed without knowledge of the causal SNPs.

Polygenic score construction and trait prediction

Request a detailed protocol

For all GWAS designs described above, we used p-value thresholding followed by clumping to choose sets of roughly independent SNPs to build PGS. We considered a logarithmically-spaced range of p-values: 10⁻⁸, 10⁻⁷, 10⁻⁶, 10⁻⁵, 10⁻⁴, 10⁻³, and 10⁻² (or a subset if no SNP reached that significance level). We then used plink’s clumping procedure (with the flag --clump) with LD threshold $r^{2}$ < 0.1 (using 10,000 randomly selected unrelated WB samples as a reference for LD structure) and physical distance threshold of >1MB. The selected SNPs were then used to calculate PGS for individuals in the prediction sets, by summing the allelic counts weighted by their estimated effect sizes (log of the odds ratios in the case of binary traits) using plink (with the flag --score). In a subset of cases, we also calculated polygenic scores using LDpred assuming all loci are causal (Vilhjálmsson et al., 2015). To evaluate prediction accuracy, we calculated the incremental $R^{2}$ : we first determined $R^{2}$ in a regression of the phenotype to the covariates, and then calculated the change in $R^{2}$ when including the PGS as a predictor. For binary traits, we calculated the incremental area under the receiver operator curve (AUC).

Estimating heritability and genetic correlation

Request a detailed protocol

We calculated SNP heritability across sex, age and SES groups for diastolic blood pressure, BMI and years of schooling, respectively (as described in the section ‘GWAS by sample characteristics’) as well as genetic correlations across pairs of groups: we first performed GWAS using all unrelated WB individuals in each group. We then used the GWAS summary statistics to perform LD score regression with LD scores computed from the 1000 Genomes European-ancestry samples (Bulik-Sullivan et al., 2015).

Appendix 1

1 Prediction accuracies of polygenic scores based on standard and sib-GWAS

1.1 Overview of derived results

In the main text, we compare the prediction accuracies of polygenic scores (PGS) based on a standard GWAS of unrelated individuals and a GWAS based on sibling differences, for a number of traits. Here, we describe how this comparison is implemented, and how indirect effects and assortative mating manifest in this comparison.

Matching standard and sib-based prediction accuracies

Current standard GWAS are based on huge sample sizes, leading to less noisy estimates than are afforded by family association studies such as those based on sib-differences, which are typically much smaller. This difference in precision needs to be taken into account in making comparisons between the prediction accuracy of scores derived from the two approaches. We show that under a vanilla additive model with no assortative mating, indirect effects, population structure (or other complications), and if the standard GWAS is subsampled to a sample size

n^{*} \approx \frac{1}{1 + (1 - h^{2}) (1 - 2 ρ_{s i b})} n^{p a i r s},

where $n^{p a i r s}$ is the number of sib pairs, $h^{2}$ is the heritability and $ρ_{s i b}$ is the correlation in environmental effects experienced by siblings, the two study designs are expected to have the same (out-of-sample) prediction accuracy (see Section 1.2). This analytic result is not that useful in practice, however; in particular, it requires prior knowledge about the extent to which environmental effects correlate among siblings. Instead, we took an empirical approach to match the prediction accuracies in the two approaches: following Wood et al. (2014), we subsampled the regular GWAS to match the median standard errors of the sib-GWAS. As we show in Section 1.2.3, under our vanilla model, we then expect equal out-of-sample prediction accuracies for polygenic scores derived from the two study designs.

Indirect parental effects

In the presence of indirect parental effects, the out-of-sample prediction accuracy takes a simple form. For a polygenic score based on a standard GWAS, we obtain

E [R_{u r}^{2}] = τ^{2} \frac{1}{1 + c},

where $τ^{2}$ is the ratio of the variance in the trait due to both direct effects and indirect effects of transmitted parental alleles over the total phenotypic variance; and $c$ is a term representing the noise-to-signal ratio in a standard GWAS. For the polygenic score based on sib-GWAS, we obtain

E [R_{s i b}^{2}] = {(1 + ρ \frac{σ_{η}}{σ_{β}})}^{2} h_{β}^{2} \frac{1}{1 + c τ^{2} / h_{β}^{2}} .

where $σ_{β}^{2}$ and $σ_{η}^{2}$ are the variances of random direct and indirect effects, respectively, $ρ$ is the correlation between direct and indirect effects, and $h_{β}^{2}$ is the proportion of the phenotypic variance explained by direct effects. Our results suggest that under plausible conditions, the presence of indirect effects would lead to higher prediction accuracy in a standard GWAS. This result holds whether direct and indirect effects are positively correlated, uncorrelated or even somewhat negatively correlated (Appendix 1—figure 9).

Assortative mating

We investigated several models of assortative mating by simulation. Standard GWAS-based polygenic scores have greater prediction accuracies than those based on sib-GWAS when the parental phenotypes are positively correlated, and the reverse is true if they are negatively correlated (Appendix 1—figure 10 A,B). The relative difference in prediction accuracies of the two study designs grows with the inclusion of more SNPs in the polygenic score model (Appendix 1—figure 10 D,F).

In our analytic model, we ignored the ascertainment step of our study design, in which it is decided which SNPs to include in the polygenic score. We assumed that SNPs are pre-ascertained and that the set of ascertained SNPs includes all causal ones. In a subset of simulations, we implemented the ascertainment step based on an independent simulated GWAS (see below). In both settings, we refer (somewhat loosely therefore) to the regression on ascertained SNPs in a sample of unrelated individuals as ‘standard GWAS’ and the regression of the difference in phenotypes on the difference in sib genotypes as ‘sib-GWAS.’

1.2 Picking the sample size of the standard GWAS to match the prediction accuracy of the score based on the sib-GWAS

We look for the sample size $n^{*}$ of a standard GWAS performed on sample of unrelated individuals such that, under our vanilla model, the resulting polygenic score has the same (out-of-sample) prediction accuracy as the polygenic score obtained from a sib-GWAS with sample size $n_{p a i r s}$ . We begin by assuming that all causal sites $i$ are known; that they are unlinked; that they have only additive, direct effects on the phenotype; and that there is no population stratification or assortative mating. We first find the sampling variance of the effect size estimate for a single site obtained from each of the two study designs. We then examine (and ultimately match) the prediction accuracy of the polygenic scores obtained from effect sizes estimated in the estimation sets, ${\hat{β}}_{u r}, {\hat{β}}_{s i b}$ , on a new, independent prediction sample of unrelated individuals ${(x^{'}, y^{'})}$ .

1.2.1 Sampling error of the estimated effect size at a single site

Our model for the phenotypic value $y$ is

y = g + e

where $e$ is a Normally distributed environmental effect (which includes all sources of random noise) and

g = β_{0}^{u r} + \sum_{i} β_{i} x_{i}

where $x_{i} \in {0, 1, 2}$ are random genotypes. The genotype is coded as the the number of alleles with effect $β_{i}$ carried by the individual at site $i$ . Effect sizes $β = {β_{i}}$ are treated as fixed parameters throughout (except when noted otherwise in the very last step leading to Equation 23). We can rewrite our model to focus on the effect size at a single site $i$ :

y = β_{0} + β_{i} x_{i} + ϵ_{i},

where

ϵ_{i} = g - β_{i} x_{i} + e,

with variance

V a r [ϵ_{i}] = V a r [g - β_{i} x_{i}] + V a r [e] = V a r [y] - β_{i}^{2} V a r [x_{i}]

In an OLS regression, the standard error for the effect of an allele at site $i$ is

V a r [{\hat{β}}_{i}^{u r}] = \frac{V a r [ϵ_{i}]}{(n - 1) V a r [x_{i}]} = \frac{V a r [y] - β_{i}^{2} V a r [x_{i}]}{(n - 1) V a r [x_{i}]},

where $n$ is the sample size and ${\hat{β}}^{u r}$ denotes that the estimate was obtained using a sample of unrelated individuals. In sib-GWAS, our model for site $i$ is

Δ y = β_{0}^{s i b} + β_{i} Δ x_{i} + Δ ϵ_{i},

with variance

V a r [Δ ϵ_{i}] = V a r [Δ g - β_{i} Δ x_{i}] + V a r [Δ e] =

V a r [Δ g] + β_{i}^{2} V a r [Δ x_{i}] - 2 β_{i}^{2} V a r [Δ x_{i}] + V a r [Δ e] .

Recall that for siblings (denoted with subscripts $A$ and $B$ ), we expect

C o v [x_{i, A}, x_{i, B}] = \frac{1}{2} V a r [x_{i}],

C o v [g_{A}, g_{B}] = \frac{1}{2} V a r [g] .

Plugging these back in, we obtain

V a r [Δ ϵ_{i}] = V a r [g] - β_{i}^{2} V a r [x_{i}] + 2 V a r [e] (1 - ρ_{s i b})

where $ρ_{s i b} = C o r [e_{A}, e_{B}]$ is the correlation in environmental effects between sibs. The variance of the estimated effect size in sib-GWAS is therefore

V a r [{\hat{β}}_{i}^{s i b}] = \frac{V a r [Δ ϵ_{i}]}{(n_{p a i r s} - 1) V a r [Δ x_{i}]} = \frac{V a r [y] - β_{i}^{2} V a r [x_{i}] + V a r [e] (1 - 2 ρ_{s i b})}{(n_{p a i r s} - 1) V a r [x_{i}]} .

1.2.2 Sample size required for equal prediction accuracy

We measure prediction accuracy as the expected squared correlation between polygenic scores $\hat{g}$ and phenotypic values in an independent prediction set of unrelated individuals, denoted ${(x^{'}, y^{'})}$ ,

E_{{(x^{'}, y^{'})}} [R^{2}] = \frac{C o v^{2} [\hat{g} (x^{'}), y^{'}]}{V a r [y^{'}] V a r [\hat{g} (x^{'})]},

To incorporate randomness both in the estimation set (summarized by the Multivariate Normal distribution of $\hat{β}$ ) and the prediction set ${(x^{'}, y^{'})}$ , we will require

E_{{\hat{β}}^{u r} (n^{*})} [E_{{(x^{'}, y^{'})}} [R^{2}]] \overset{!}{=} E_{{\hat{β}}^{s i b} (n^{p a i r s})} [E_{{(x^{'}, y^{'})}} [R^{2}]]

where $\hat{β} (n)$ is a set ${{\hat{β}}_{i}}$ estimated in a GWAS with sample size $n$ . Equivalently,

E_{{\hat{β}}^{s i b}} [\frac{C o v^{2} [{\hat{g}}_{s i b} (x^{'}), y^{'}]}{V a r [{\hat{g}}_{s i b} (x^{'})]}] \overset{!}{=} E_{{\hat{β}}^{u r}} [\frac{C o v^{2} [{\hat{g}}_{u r} (x^{'}), y^{'}]}{V a r [{\hat{g}}_{u r} (x^{'})]}],

where we left out the sample sizes for brevity, and $V a r [y^{'}]$ was cancelled out. Finally, we can replace Equation 4 by its first order Taylor approximation to get the requirement

\frac{E_{\hat{β}} [C o v_{{(x^{'}, y^{'})}} [{\hat{g}}_{s i b} (x^{'}), y^{'}]]^{2}}{E_{\hat{β}} [V a r_{{(x^{'}, y^{'})}} [{\hat{g}}_{s i b} (x^{'})]]} \overset{!}{=} \frac{E_{\hat{β}} [C o v_{{(x^{'}, y^{'})}} [{\hat{g}}_{u r} (x^{'}), y^{'}]]^{2}}{E_{\hat{β}} [V a r_{{(x^{'}, y^{'})}} [{\hat{g}}_{u r} (x^{'})]]} .

We solve Equation 4 for a sample size $n^{*}$ to be used for estimation of the polygenic score in a standard GWAS that satisfies Equation 4. We note that if the vector of estimates $\hat{β}$ is given, then

\begin{array}{ll} C o v_{{(x^{'}, y^{'})}} [y^{'}, \hat{g} (x^{'}) | \hat{β}] = C o v_{{(x^{'}, y^{'})}} [g (x^{'}), g (x^{'}) + \sum_{i}^{m} x_{i}^{'} ({\hat{β}}_{i} - β_{i}) | \hat{β}] = \\ V a r_{{(x^{'}, y^{'})}} [g (x^{'}) | \hat{β}] + \sum_{i}^{m} C o v_{{(x^{'}, y^{'})}} [β_{i} x_{i}^{'}, ({\hat{β}}_{i} - β_{i}) x_{i}^{'} | \hat{β}] = \sum_{i}^{m} V a r [x_{i}^{'}] β_{i} {\hat{β}}_{i} . \end{array}

Since for every $i$ , we have

E [{\hat{β}}_{i}^{u r}] = E [{\hat{β}}_{i}^{s i b}] = β_{i},

we obtain

E_{{\hat{β}}^{s i b}} [C o v [y^{'}, {\hat{g}}_{s i b} (x^{'}) | {\hat{β}}^{s i b}]] = \sum_{i}^{m} V a r [x_{i}^{'}] β_{i}^{2} = E_{{\hat{β}}^{u r}} [C o v [y^{'}, {\hat{g}}_{u r} (x^{'}) | {\hat{β}}^{u r}]],

which turns the requirement of Equation 5 into

E_{{\hat{β}}^{s i b}} [V a r_{{(x^{'}, y^{'})}} [{\hat{g}}_{s i b} (x^{'})]] \overset{!}{=} E_{{\hat{β}}^{u r}} [V a r_{{(x^{'}, y^{'})}} [{\hat{g}}_{u r} (x^{'})]],

or simply

\sum_{i}^{m} V a r [x_{i}] V a r [{\hat{β}}_{i}^{u r}] \overset{!}{=} \sum_{i}^{m} V a r [x_{i}] V a r [{\hat{β}}_{i}^{s i b}] .

Plugging the sampling variance results from Equation 2 and Equation 3 into Equation 7 and reordering, we obtain

\frac{n^{*} - 1}{n_{p a i r s} - 1} = \frac{\sum_{i}^{m} V a r [y] - β_{i}^{2} V a r [x_{i}]}{\sum_{i}^{m} V a r [y] - β_{i}^{2} V a r [x_{i}] + V a r [e] (1 - 2 ρ_{s i b})},

or, assuming that the trait is polygenic such that $m ≫ 1$ ,

\frac{n^{*}}{n^{p a i r s}} \approx \frac{1}{1 + (1 - h^{2}) (1 - 2 ρ_{s i b})} .

Equation 8 can in principle be applied to the estimation of $ρ_{s i b}$ for a given trait, under our model assumptions, and given an independent estimate of $h^{2} .$

1.2.3 Empirical matching of standard errors

The result of Equation 8 is the same as we would obtain if we required

\forall i V a r [{\hat{β}}_{i}^{s i b} (x_{i})] \overset{!}{=} V a r [{\hat{β}}_{i}^{u r} (x_{i}^{s i b})]

without taking into account randomness in the prediction set. In practice (and in the results shown in the main text), we have no prior knowledge about $ρ_{s i b}$ and instead we find a sample size $n^{*}$ for the standard GWAS such that

m e d i a n_{{s i t e s i}} (V a r [{\hat{β}}_{i}^{s i b} (x)]) \overset{!}{=} m e d i a n_{{s i t e s i}} (V a r [{\hat{β}}_{i}^{u r} (x)])

We note that the condition in Equation 9 is approximately met because, if we assume that $y$ is a highly polygenic trait where

\forall i β_{i}^{2} V a r [x_{i}] << V a r [y],

then, if for one site $j$ , $n^{*}$ satisfies

V a r [{\hat{β}}_{j}^{s i b} (x)] = V a r [{\hat{β}}_{j}^{u r} (x)] = \frac{D (n^{*})}{V a r [x_{j}]}

such that $D (n^{*})$ is the same for sib-GWAS and standard GWAS, then for all sites $D (n^{*}) = \frac{V a r [y]}{n^{*} - 1}$ is the same, namely,

\forall i V a r [{\hat{β}}_{i}^{s i b} (x)] = V a r [{\hat{β}}_{i}^{u r} (x)] = \frac{D (n^{*})}{V a r [x_{i}]}

Equation 10 can therefore be thought of as using a weighted-median to estimate $n^{*}$ where each site $i$ is weighted by $\frac{1}{V a r [x_{i}]}$ . In conclusion, the requirement of Equation 10 leads to equal prediction accuracy of standard and sib-GWAS under the vanilla model assumptions. We note further that in the main text (Figure 3), to follow common practice, we use incremental $R^{2}$ throughout rather than $R^{2}$ . However, as we show in Appendix 1—figure 16, using $R^{2}$ instead gives highly similar qualitative results.

1.3 Indirect parental effects

1.3.1 Distribution of the effect size estimate at a single site

We consider an additive model with direct effects as well as indirect parental effects, assuming no interaction between the parents and the polygenic score of the children and ignoring possible indirect effects of siblings on each other. The other assumptions from the previous section—for example independent segregation of alleles across sites—remain. We start by considering the model

y = β_{0} + g + n + e

where $g$ is the sum of direct effects in an individual with genotype (effect-allele count) $x_{i}$ at each site $i$ ,

g = \sum_{i}^{m} β_{i} x_{i},

and

n = \sum_{i}^{m} η_{i} (x_{i} + {\tilde{x}}_{i}^{m} + {\tilde{x}}_{i}^{p})

is the sum of parental indirect effects, with overall parental effect allele count $x_{i} + {\tilde{x}}_{i}^{p} + {\tilde{x}}_{i}^{m}$ at each site, where ${\tilde{x}}_{i}^{m}$ is the untransmitted maternal effect allele count, and ${\tilde{x}}_{i}^{p}$ the untransmitted paternal effect allele count, with ${\tilde{x}}_{i}^{m}, {\tilde{x}}_{i}^{p} \in {0, 1}$ . As we show, when we choose the standard GWAS sample size $n^{*}$ such that the sampling error of the effect size estimates matches that of the sib-GWAS, the prediction accuracies of the two polygenic scores differ in an independent sample: unless there is a large, negative correlation between indirect and direct effects, the polygenic score from standard GWAS is expected to outperform the one based on sib-GWAS.

We first examine the distribution of an estimated effect size of $x_{i}$ on the phenotype. The OLS regression for a single site in a standard GWAS follows Equation 1 and can be rewritten as

y = β_{0} + (β_{i} + η_{i}) x_{i} + η_{i} ({\tilde{x}}_{i}^{p} + {\tilde{x}}_{i}^{m}) + ϵ_{i}

with

ϵ_{i} = g + n + e - (β_{i} + η_{i}) x_{i} - η_{i} ({\tilde{x}}_{i}^{p} + {\tilde{x}}_{i}^{m}) .

By the assumption of no assortative mating or other population structure,

C o v [{\tilde{x}}_{i}^{p}, {\tilde{x}}_{i}^{m}] = C o v [x_{i}, {\tilde{x}}_{i}^{m}] = C o v [x_{i}, {\tilde{x}}_{i}^{p}] = 0 .

It directly follows that under the generative model specified by Equation 11, the OLS regression of $y$ to $x_{i}$ and ${\tilde{x}}_{i}^{p} + {\tilde{x}}_{i}^{m}$ is a regression involving two independent variables. Therefore, ${\hat{β}}_{i}^{u r}$ is Normally distributed with expectation

E [{\hat{β}}_{i}^{u r}] = β_{i} + η_{i} .

We next calculate the variance of ${\hat{β}}_{i}^{u r}$ . From Equation 12 and

V a r [{\tilde{x}}_{i}^{m} + {\tilde{x}}_{i}^{p}] = V a r [x_{i}],

we obtain

V a r [ϵ_{i}] = V a r [y] + {(β_{i} + η_{i})}^{2} V a r [x_{i}] + η_{i}^{2} V a r [x_{i}] - 2 C o v [g + n, (β_{i} + η_{i}) x_{i}] - 2 C o v [n, η_{i} ({\tilde{x}}_{i}^{m} + {\tilde{x}}_{i}^{p})] =

= V a r [y] - V a r [x_{i}] (β_{i}^{2} + 2 β_{i} η_{i} + 2 η_{i}^{2}) .

Finally,

V a r [{\hat{β}}_{i}^{u r}] = \frac{V a r [ϵ_{i}]}{(n - 1) V a r [x_{i}]} = \frac{V a r [y] - V a r [x_{i}] (β_{i}^{2} + 2 β_{i} η_{i} + 2 η_{i}^{2})}{(n - 1) V a r [x_{i}]} .

In sib regression, we have

Δ y = Δ g + Δ e

since indirect parental effects cancel out when taking the difference between siblings (as siblings have the same parental effect allele count). Thus, the expected estimate is the same as it was in the absence of indirect effects. Using the same considerations as in Section 1.2 for the variance in sib differences, we obtain

{\hat{β}}_{i}^{s i b} \sim N (β_{i}, \frac{V a r [g] - β_{i}^{2} V a r [x_{i}] + V a r [e] (1 - 2 ρ_{s i b})}{(n_{p a i r s} - 1) V a r [x_{i}]}),

where $ρ_{s i b}$ is again the correlation in environmental effects between siblings.

1.3.2 Polygenic score prediction accuracy

We now examine the difference in prediction accuracies of ${\hat{g}}^{u r}$ and ${\hat{g}}^{s i b}$ after matching

V a r [{\hat{β}}_{i}^{u r}] \overset{!}{=} V a r [{\hat{β}}_{i}^{s i b}]

by choosing a standard GWAS sample size $n^{*}$ that empirically satisfies the condition, as we do in the main text (see also Section 1.2.3).

We can derive the expected prediction accuracy by averaging over both the estimation set (which we again shorthand as the distribution of $\hat{β}$ ) and the prediction set ${(x^{'}, y^{'})}$ . By the law of total expectation,

E [R^{2}] = E_{\hat{β}} [E_{{(x^{'}, y^{'})}} [R^{2}]] = E_{\hat{β}} [\frac{C o v_{{(x^{'}, y^{'})}}^{2} [\hat{g} (x^{'}), y^{'} | \hat{β}]}{V a r_{{(x^{'}, y^{'})}} [y^{'} | \hat{β}] V a r_{{(x^{'}, y^{'})}} [\hat{g} (x^{'}) | \hat{β}]}] \approx

\approx \frac{E_{\hat{β}} {[C o v_{{(x^{'}, y^{'})}} [\hat{g} (x^{'}), y^{'} | \hat{β}]]}^{2}}{V a r_{{(x^{'}, y^{'})}} [y^{'} | \hat{β}] E_{\hat{β}} [V a r_{{(x^{'}, y^{'})}} [\hat{g} (x^{'}) | \hat{β}]]},

where the last step is an approximation of the expectation of ratio by its first-order Taylor expansion, a ratio of expectations. The numerator of Equation 15 is

E_{\hat{β}} {[C o v_{{(x^{'}, y^{'})}} [\hat{g} (x^{'}), y^{'} | \hat{β}]]}^{2} = E_{\hat{β}} [\sum_{i}^{m} (β_{i} + η_{i}) {\hat{β}}_{i} C o v_{{(x^{'}, y^{'})}} {[x_{i}^{'}, x_{j}^{'} | \hat{β]}]}^{2} =

= E_{\hat{β}} {[\sum_{i}^{m} V a r [x_{i}] (β_{i} + η_{i}) {\hat{β}}_{i}]}^{2} =

= {(\sum_{i}^{m} V a r [x_{i}] (β_{i} + η_{i}) E [{\hat{β}}_{i}])}^{2} .

The terms in the denominator of Equation 15 are

V a r_{{(x^{'}, y^{'})}} [y^{'} | \hat{β}] = V a r [y]

and

E_{\hat{β}} [V a r_{{(x^{'}, y^{'})}} [\hat{g} (x^{'}) | \hat{β}] = E_{\hat{β}} [\sum_{i}^{m} V a r [x_{i}] {\hat{β}}_{i}^{2}] = \sum_{i}^{m} V a r [x_{i}] (E {[{\hat{β}}_{i}]}^{2} + V a r [{\hat{β}}_{i}]) .

Plugging Equations 16,17,18 back into Equation 15, we obtain

E [R^{2}] \approx \frac{{(\sum_{i}^{m} V a r [x_{i}] (β_{i} + η_{i}) E [{\hat{β}}_{i}])}^{2}}{V a r [y] (\sum_{i}^{m} V a r [x_{i}] V a r [{\hat{β}}_{i}] + \sum_{i}^{m} V a r [x_{i}] E {[{\hat{β}}_{i}]}^{2})} .

We note that

\tilde{C} := V a r [y] \sum_{i}^{m} V a r [x_{i}] V a r [{\hat{β}}_{i}]

is the same for sib-GWAS and standard GWAS under the requirement of Equation 14. We therefore have

E [R_{u r}^{2}] \approx \frac{{(\sum_{i}^{m} V a r [x_{i}] {(β_{i} + η_{i})}^{2})}^{2}}{\tilde{C} + V a r [y] \sum_{i}^{m} V a r [x_{i}] {(β_{i} + η_{i})}^{2}},

and

E [R_{s i b}^{2}] \approx \frac{{(\sum_{i}^{m} V a r [x_{i}] (β_{i} + η_{i}) β_{i})}^{2}}{\tilde{C} + V a r [y] \sum_{i}^{m} V a r [x_{i}] β_{i}^{2}} .

If we denote the proportion of the phenotypic variance explained by direct effects by

h_{β}^{2} := \frac{\sum_{i}^{m} V a r [x_{i}] β_{i}^{2}}{V a r [y]},

the proportion of the phenotypic variance explained by indirect effects of transmitted parental alleles by

τ_{η}^{2} := \frac{\sum_{i}^{m} V a r [x_{i}] η_{i}^{2}}{V a r [y]},

and the proportion of phenotypic variance explained by both direct and indirect effects of transmitted alleles by

τ^{2} := \frac{\sum_{i}^{m} V a r [x_{i}] {(β_{i} + η_{i})}^{2}}{V a r [y]}

then Equation 20 can be written as

E [R_{u r}^{2}] \approx τ^{2} \frac{1}{1 + c},

where we defined

c := \frac{\sum_{i}^{m} V a r [x_{i}] V a r [{\hat{β}}_{i}]}{\sum_{i}^{m} V a r [x_{i}] {(β_{i} + η_{i})}^{2}} .

Here, $c$ can be thought of as a summary of the noise-to-signal ratio, with respect to the signal coming from both direct and indirect effects of transmitted alleles. If we consider effects $β$ and $η$ as random, treating results obtained thus far as conditional on $β$ and $η$ , and further assume that effects are i.i.d. across sites (implying, in particular, that effect sizes and allele frequencies are independent),

(\begin{array}{ll} β_{i} \\ η_{i} \end{array}) \sim ((\begin{array}{ll} 0 \\ 0 \end{array}), (\begin{array}{ll} σ_{β}^{2} & ρ σ_{β} σ_{η} \\ ρ σ_{β} σ_{η} & σ_{η}^{2} \end{array})),

the expectation of the numerator of Equation 21 is

E_{β, η} [\sum_{i}^{m} V a r [x_{i}] β_{i} (β_{i} + η_{i}) | β, η] = \sum_{i}^{m} V a r [x_{i}] E_{β_{i}, η_{i}} [β_{i}^{2} + β_{i} η_{i}] = \sum_{i}^{m} V a r [x_{i}] (σ_{β}^{2} + ρ σ_{β} σ_{η})

and thus Equation 21, in expectation, is:

E [R_{s i b}^{2}] \approx E_{β, η} [E [R_{s i b}^{2} | β, η]] = {(1 + ρ \frac{σ_{η}}{σ_{β}})}^{2} h_{β}^{2} \frac{1}{1 + c / α} .

where

α := h_{β}^{2} / τ^{2} = \frac{\sum_{i}^{m} V a r [x_{i}] β_{i}^{2}}{\sum_{i}^{m} V a r [x_{i}] {(β_{i} + η_{i})}^{2}} .

We examined the fit of this prediction to simulated data. Specifically, we ran simulations to estimate effect sizes in a sib-GWAS and in a standard GWAS, after choosing $n^{*}$ to match their sampling variances. Finally, we used the polygenic scores to predict phenotypic values in a sample of unrelated individuals (see Section 1.3.3 for further detail).

Appendix 1—figure 9 A,C,D show the analytic result alongside simulation results, for different correlation coefficients between indirect and direct effect sizes. Even in the absence of a correlation between indirect and direct effect sizes, the polygenic score based on standard GWAS outperforms the polygenic score based on sib-GWAS.

To understand this behavior and dependency of the $\frac{R_{s i b}^{2}}{R_{u r}^{2}}$ ratio on other parameters, we divide Equation 23 by Equation 22 and obtain

E [\frac{R_{s i b}^{2}}{R_{u r}^{2}}] \approx \frac{E [R_{s i b}^{2}]}{E [R_{u r}^{2}]} \approx {(1 + ρ \frac{σ_{η}}{σ_{β}})}^{2} α \frac{1 + c}{1 + c / α} .

Noting further that

{(1 + ρ \frac{σ_{η}}{σ_{β}})}^{2} α = {(\frac{σ_{β} + ρ σ_{η}}{σ_{β}})}^{2} \frac{σ_{β}^{2}}{σ_{β}^{2} + 2 ρ σ_{β} σ_{η} + σ_{η}^{2}} = 1 - (1 - ρ^{2}) \frac{τ_{η}^{2}}{τ^{2}},

we obtain

E [\frac{R_{s i b}^{2}}{R_{u r}^{2}}] \approx [1 - (1 - ρ^{2}) \frac{τ_{η}^{2}}{τ^{2}}] \frac{1 + c}{1 + c \frac{τ^{2}}{h_{β}^{2}}} .

A few conclusions emerge from Equation 24 and the accompanying simulations. First, the sib-GWAS based polygenic score will outperform the standard GWAS-based polygenic score only if direct and indirect effects are strongly negatively correlated (see Appendix 1—figure 9A-D for illustration). Second, the term

\frac{1 + c}{1 + c \frac{τ^{2}}{h_{β}^{2}}} = \frac{1 + \frac{\sum_{i}^{m} V a r [{\hat{β}}_{i}] V a r [x_{i}]}{τ^{2}}}{1 + \frac{\sum_{i}^{m} V a r [{\hat{β}}_{i}] V a r [x_{i}]}{h_{β}^{2}}}

can be interpreted as the dependence on the noise-to-signal ratio (where the signals are the proportions of phenotypic variance explained by direct and indirect effects of transmitted alleles). For a given sampling variance (matched across the two study designs), the extent of the signal will differ between standard GWAS and sib-GWAS. Importantly, the sampling variance influences the ratio of prediction accuracies. If indirect effects do not exist or make negligible contributions to the trait in question, then the ratio of prediction accuracies is expected to be close to one. In the presence of indirect effects, however, the magnitude of the deviation from one depends on the relationship between direct and indirect effects (and their covariance) as well as on the (matched) sampling variance. Simulations of several parameter combinations suggest that the overall effect of this dependence on the noise-to-signal ratio is a decrease in $R_{s i b}^{2} / R_{u r}^{2}$ as noise increases; as more SNPs are included in the polygenic scores, the advantage of the standard GWAS-based polygenic score over that of the sib-GWAS grows larger (Appendix 1—figure 9 E-H). These considerations inform the interpretation of patterns observed in Figure 3C–F of the main text.

1.3.3 Simulations of indirect effects

For each set of simulated individuals (discovery, estimation and prediction sets), we first simulated mother-father pairs, assigning parental alleles from $B e r n o u l l i (p_{i})$ , where $p_{i}$ denotes the allele frequency at site $i$ . We then sampled the parental alleles at random to generate offspring (one offspring per each mother-father pair to simulate a sample of unrelated individuals and two offspring to generate sibling pairs). Phenotypes of the offspring were assigned under an additive model, sampling from a Normal distribution with mean

\sum_{i}^{m} β_{i} x_{i} + η_{i} (x_{i}^{p} + x_{i}^{m})

(where $x_{i}^{m}$ and $x_{i}^{p}$ are the maternal and paternal effect allele counts, respectively) and variance $σ_{e}^{2},$ representing the total variance of environmental effects. When there is no correlation between direct and indirect effects, $σ_{e}^{2} = 1 - h_{β}^{2} - 2 τ_{η}^{2}$ . Using this approach, we generated a set of sibling pairs and estimated SNP effect sizes from these simulated data using a sib-GWAS. We calculated $n^{*}$ as follows: we simulated sets of unrelated individuals with a range of sample sizes. In each set, we performed a simple linear regression of the phenotypic values on the genotypes. We then estimated a linear relationship between the inverse of the median standard error of effect size estimates (as a dependent variable) and the square root of the sample size. Using this linear relationship, we predicted the sample size for the unrelated set that gives a median standard error equal to the median standard error of sib-GWAS effect size estimates ( $n^{*}$ ). Finally, we simulated a set of unrelated individuals with sample size $n^{*}$ and compared the prediction accuracy ( $R^{2}$ ) of the polygenic score based on standard GWAS on this sample with the one obtained from sib-GWAS.

We additionally investigated the effect of the number of SNPs included in the polygenic scores. For this analysis, we sorted the SNPs based on the association p-value obtained in an independent simulated set of unrelated individuals.

In these simulations, we used the following parameter values:

The ratio of the phenotypic variance accounted for by direct effects versus by indirect effects ( $h_{β}^{2} / τ_{η}^{2}$ ): 5
The phenotypic variance explained by offspring and parental alleles, given no correlation between direct and indirect effects ( $h_{β}^{2} + 2 τ_{η}^{2}$ ): 0.25 or 0.5
The ratio of the variance of direct effects to the variance of indirect effects ( $σ_{β}^{2} / σ_{η}^{2}$ ): 5
Allele frequencies, $p$ , drawn from a truncated exponential distribution, truncated on the left such that the minimum allele frequency is 1%.
The number of loci, assumed independent (i.e., in linkage equilibrium): 100 (all causal), or 10,000 (all causal) or 10,000 (20% causal)
SNP effect sizes drawn as
- $(\begin{array}{ll} β_{i} \\ η_{i} \end{array}) \sim N ((\begin{array}{ll} 0 \\ 0 \end{array}), (\begin{array}{ll} σ_{β}^{2} & ρ σ_{β} σ_{η} \\ ρ σ_{β} σ_{η} & σ_{η}^{2} \end{array})),$

where $ρ$ is the correlation between direct and indirect effect sizes. Effects sizes were then re-scaled to satisfy $\sum_{i}^{m} 2 β_{i}^{2} p_{i} (1 - p_{i}) = h_{β}^{2}$ and $\sum_{i}^{m} 2 η_{i}^{2} p_{i} (1 - p_{i}) = τ_{η}^{2}$ . Effects were set to 0 for non-causal loci.

The number of sibling pairs for sib GWAS: 10,000
The number of unrelated individuals for prediction: 10,000
The number of unrelated individuals for discovery GWAS (i.e., to decide which SNPs to include): 20,000
Number of iterations used to estimate $n^{*}$ and $R^{2}$ for a given set of parameters: 10

1.4 Assortative mating

We consider assortative mating with regard to a phenotype, whereby the parents of individuals were more likely to mate if they were similar with respect to that phenotype. This process generates a correlation between genetic variants that contribute to the phenotype (i.e., linkage disequilibrium). Consequently, in a standard GWAS, the effect sizes of causal SNPs will partially capture the effect of other causal SNPs as well. Estimated effect sizes are thus expected to be inflated under positive assortative mating (mating of similar individuals) and deflated under negative assortative mating (mating of dissimilar individuals). In turn, in a sib-GWAS, the estimates are in expectation unaffected by assortative mating, because genetic differences between siblings arise from random Mendelian segregation in the parents.

1.4.1 Simulations of assortative mating

We used simulations to examine the phenotypic prediction accuracies of polygenic scores based on sib- and standard GWAS under a model with assortative mating (assuming no indirect effects or population stratification beyond assortative mating); to this end, we considered a sample of unrelated individuals, varying the degree of correlation between parental phenotypes $ρ_{a}$ . Similar to our simulations for indirect effects (Section 1.3.3), we first simulated the estimation procedure in a sibling-based and in a standard GWAS (with sample size $n^{*}$ ). We then computed the prediction accuracy $R^{2}$ in an independent sample of unrelated individuals (see ‘Further simulation details’ below).

We first considered the simple case of a single generation of assortative mating. In the presence of positive assortative mating ( $ρ_{a} > 0$ ), polygenic scores based on standard GWAS outperform those based on sib-GWAS, whereas the opposite is true in the case of negative assortative mating ( $ρ_{a} < 0$ ) (Appendix 1—figure 10 A). In simulations of two generations of assortative mating, the gap between the prediction accuracies of scores based on standard and sib-GWAS (Appendix 1—figure 10 B) widens, suggesting that our qualitative findings apply to scenarios of sustained assortative mating as well.

We further investigated prediction accuracy as a function of the number of SNPs included in the polygenic scores, by progressively increasing the p-value threshold, using p-values obtained from an independent GWAS in unrelated samples (similar to our analysis in Figure 3). We considered two genetic architectures scenarios: (i) in which all SNPs are causal and (ii) the case in which 20% of of SNPs are causal (leading polygenic scores to include non-causal SNPs). Under both scenarios, the gap in prediction accuracies between standard and sib-GWAS grows with the number of SNPs (Appendix 1—figure 10 C-F).

Further simulation details

We simulated parental and offspring alleles as described for indirect effects in Section 1.3.3. To mimic assortative mating between parents, we first simulated i.i.d. genotypes (with effect allele counts $x_{i}$ at each SNP $i$ ) and randomly assigned ‘mother’ and ‘father’ labels to each individual. We then generated corresponding parental phenotypes under an additive model as

N (\sum_{i}^{m} β_{i} x_{i}, \sqrt{1 - h^{2}})

where $β_{i}$ is the effect size of SNP $i$ , and $h^{2}$ is the heritability. The same model was used to generate offspring phenotypes.

To mimic the assortative mating process, we induced a given correlation between parental phenotypes, $ρ_{a}$ , by paring mothers and fathers as follows: we first generated a random matrix

(\begin{array}{ll} u_{m, i} \\ u_{p, i} \end{array}) \sim N ((\begin{array}{ll} {\bar{y}}_{m} \\ {\bar{y}}_{p} \end{array}), (\begin{array}{ll} σ_{y_{m}}^{2} & ρ_{a} σ_{y_{m}} σ_{y_{p}} \\ ρ_{a} σ_{y_{m}} σ_{y_{p}} & σ_{y_{p}}^{2} \end{array})),

where ${\bar{y}}_{m}$ and ${\bar{y}}_{p}$ are the average phenotypes of mothers and fathers, respectively, $σ_{y_{m}}$ and $σ_{y_{p}}$ are the standard deviation of the phenotypes of mothers and fathers, respectively. We then sorted the mothers and fathers sets such that the ranks of values in $y_{m}$ and $y_{p}$ match the ranks of values in $u_{m}$ and $u_{p}$ , respectively, to obtain $c o r (y_{m}, y_{p}) \approx c o r (u_{m}, u_{p}) = ρ_{a}$ . In the case of two generations of assortative mating, we simulated the generation of the grandparents similarly. We compared the performance of polygenic scores based on standard and sib-GWAS as described in Section 1.3.3. In the simulations, we used the following parameter values:

Heritability under random mating ( $h^{2}$ ): 0.5
The number of loci, assumed independent (i.e., in linkage equilibrium) under random mating: 10,000 (all causal) or 10,000 (20% causal)
Allele frequencies, $p$ , drawn from a truncated exponential distribution, truncated on the left such that the minimum allele frequency is 1%.
SNP effect sizes set to 0 for non-causal loci and drawn as $β_{i} \sim N (0, σ^{2})$ , choosing $σ^{2}$ to satisfy $\sum_{i}^{m} 2 β_{i}^{2} p_{i} (1 - p_{i}) = h^{2}$ for causal loci.
The number of sibling pairs for sib-GWAS: 10,000
The number of unrelated individuals for prediction: 10,000
The number of unrelated individuals for discovery GWAS (i.e., to decide which SNPs to include in the polygenic score): 20,000
The number of iterations used to estimate $n^{*}$ and $R^{2}$ for a given set of parameters: 10

Appendix 1—figure 1

Download asset Open asset

Variable prediction accuracy within an ancestry group.

This figure extends Figure 1 of the main text, showing prediction accuracies based on large-scale diverse GWAS that are the union of all strata matching the number of individuals in each stratum. The numbers in parentheses show GWAS sample sizes (see Materials and methods for details). Each box and whiskers plot was computed based on 20 iterations of resampling estimation and prediction sets. Thick horizontal lines denote the medians. The polygenic scores were estimated in samples of unrelated WB individuals. Phenotypes were then predicted in distinct samples of unrelated WB individuals, stratified by sex (A), age (B) or Townsend deprivation index, a measure of SES (C). In red and green cases, polygenic scores are based on a GWAS in a sample limited to one sex, age or SES group (a 'stratum’). In black, polygenic scores are based on a diverse GWAS in a pooled sample of all strata. In blue, polygenic scores are based on a diverse GWAS in a pooled sample of all strata but downsampled to match the size of the stratified GWAS.

Appendix 1—figure 2

Download asset Open asset

Variable prediction accuracy (measured as $R^{2}$ ) within an ancestry group.

This figure mirrors Figure 1 of the main text, except for the y-axis showing $R^{2}$ values (squared correlation between polygenic score and phenotype residualized on covariates), rather than incremental $R^{2}$ . Each box and whiskers plot was computed based on 20 iterations of resampling estimation and prediction sets. Thick horizontal lines denote the medians. The polygenic scores were estimated in samples of unrelated WB individuals. Phenotypes were then predicted in distinct samples of unrelated WB individuals, stratified by sex (A), age (B) or Townsend deprivation index, a measure of SES (C). In red and green cases, polygenic scores are based on a GWAS in a sample limited to one sex, age or SES group (a 'stratum’). In blue, polygenic scores are based on a GWAS in a diverse sample of all strata downsampled to match the size of the stratified GWAS.

Appendix 1—figure 3

Download asset Open asset

Dependence on the polygenic score model.

This figure extends Figure 1 of the main text, showing the prediction accuracies as a function of the p-value threshold for inclusion of a SNP in the polygenic score when based on a pruning and thresholding approach. The higher the p-value threshold is, the more SNPs are included. Last points on the x-axis correspond to a polygenic score model based on the LDpred approach (Vilhjálmsson et al., 2015) with a prior probability of 1 on loci being causal. Shown are incremental $R^{2}$ values in different prediction sets. Points and error bars are mean and central 80% range computed based on 20 iterations of resampling estimation and prediction sets. (**A–C**) The polygenic scores were estimated in samples of unrelated WB individuals. Phenotypes were then predicted in distinct samples of unrelated WB individuals, stratified by sex (A), age (B) or Townsend deprivation index, a measure of SES (C). (**D–I**) Same as in **A-C**, but here the polygenic scores are based on a GWAS in a sample limited to one sex, age or SES group. The trends shown in Figure 1 of the main text are for p-value threshold of 10⁻⁴, and are qualitatively similar to the trends for other choices of the polygenic score model. For each trait, sample sizes are matched across all GWAS sets.

Appendix 1—figure 4

Download asset Open asset

Estimating mean effect size across strata.

SNPs were ascertained in large samples of unrelated WB individuals. The effects of trait-increasing alleles were then re-estimated in an independent set of unrelated WB individuals (that were excluded from the original GWAS) stratified by sex for diastolic blood pressure (A), by age for BMI (B) and by Townsend deprivation index, a measure of SES for years of schooling (C). Points and error bars are mean and central 80% range computed based on 20 iterations of resampling ascertainment and estimation sets, plotted as a function of the p-value threshold (for p-values obtained in the discovery GWAS).

Appendix 1—figure 5

Download asset Open asset

Variable prediction accuracy within an ancestry also seen using a linear mixed model.

This figure mirrors the last two columns in Appendix 1—figure 3, except that here, the GWAS estimates were obtained from a linear mixed model (LMM) (Loh et al., 2015). Shown are the prediction accuracies, measured as incremental $R^{2}$ , as a function of the p-value threshold for inclusion of a SNP in the polygenic score. Points and error bars are mean and central 80% range computed based on 20 iterations of resampling estimation and prediction sets. The polygenic scores are based on a GWAS in a sample limited to one sex, age or SES group. Phenotypes are then predicted in distinct samples of unrelated individuals, stratified by sex (**A,B**), age (**C,D**) or Townsend deprivation index, as a measure of SES (**E,F**). The qualitative trends are similar to those in Appendix 1—figure 3, which uses a standard linear regression with PCs (principal components of the genotype data) as a control for population structure when testing for an association between the phenotypes and genotypes. The similarity suggests that the observed differences in prediction accuracies across strata are not driven to a large degree by population structure confounding.

Appendix 1—figure 6

Download asset Open asset

Comparison of siblings and unrelated individuals in the UK Biobank with respect to age, SES, and sex ratio.

Panels show the distribution of Townsend deprivation index, a measure of SES (A), the age distribution (B), and the proportion of males (C) for the siblings and unrelated sets used in the analysis described for Figure 3 of the main text. For each sibling pair, one sibling was randomly selected for these comparisons. The asterisk symbol marks a significant difference at the 1% level between siblings and unrelated individuals, as assessed by a Mann-Whitney test. SES and age distributions are quite similar in siblings and unrelated sets, whereas the proportion of males is significantly smaller in the siblings.

Appendix 1—figure 7

Download asset Open asset

Comparison of siblings and unrelated individuals in the UK Biobank with respect to population structure.

Panels show the distribution of PCs (principal components of the genotype data) for the siblings and unrelated sets used in the analysis described for Figure 3 of the main text. For each sibling pair, one sibling was randomly selected for these comparisons. The asterisk symbol marks a significant difference at the 1% level between siblings and unrelated individuals, as assessed by a Mann-Whitney test. Despite slight but significant differences, siblings and unrelated sets are broadly similar with respect to their genetic ancestries.

Appendix 1—figure 8

Download asset Open asset

Comparison of prediction accuracies of polygenic scores based on standard and sib-GWAS for simulated traits.

This figure mirrors Figure 3B of the main text, but here plotted for 12 simulated traits. The numbers in parentheses are the heritability, the number of causal loci considered, and the simulation replicate number, respectively. Three traits were simulated for each pair of heritability and number of causal loci parameters (see Materials and methods for simulation details). Small points show the ratio of the prediction accuracies in the two designs across 30 iterations; in each iteration, we resample sets of unrelated individuals to constitute three sets for discovery, estimation and prediction. Larger points show median values.

Appendix 1—figure 9

Download asset Open asset

Simulation results for polygenic scores based on standard GWAS and sib-GWAS in the presence of indirect effects.

(**A,B**) Simulation results as a function of the correlation between direct and indirect effects, $ρ$ . Simulations were performed with $h_{β}^{2} = 0.5$ , $τ_{η}^{2} = 0.1$ , and $σ_{β}^{2} / σ_{η}^{2} = 5$ . The size of the estimation set in the sib-GWAS is 10,000, and the size of the estimation set in the standard GWAS is chosen to match sampling variances between the two study designs. The polygenic scores is based on 10,000 causal loci; its performance was evaluated in an independent set of 10,000 unrelated individuals. As long as direct and indirect effects are not strongly negatively correlated, the out of sample prediction accuracy is higher for the polygenic scores based on standard GWAS. (C) Same as (A) but with three-fold greater environmental noise. (D) Same as (A) but with 100 causal loci. In (**A–D**) points are mean ± 2 SD in 10 simulation iterations. Solid lines are values based on analytic expressions derived in Section 1.3.2. (**E–H**) Simulation results, with the same parameters as in (A) but $ρ = 0.5$ , as a function of the number of SNPs included in the polygenic scores, with all loci being causal (**E,F**), or with 20% of loci being causal (**G,H**). SNPs are added in increasing order of their association p-value in an independent set of 20,000 unrelated individuals. In both cases, the ratio of prediction accuracies of polygenic scores based on sib- versus standard GWAS becomes smaller with the inclusion of more weakly associated SNPs, a behavior qualitatively similar to observations in Figure 3 in the main text. Points are mean ± 2 SD in 10 simulations. See Section 1.3.3 for simulation details.

Appendix 1—figure 10

Download asset Open asset

Simulation results for polygenic scores based on standard GWAS and sib-GWAS in the presence of assortative mating.

(A) Simulation results as a function of the approximate correlation between parental phenotypes, $ρ_{a}$ . Simulations were performed with $h^{2} = 0.5$ under random mating. The size of the estimation set in the sib-GWAS is 10,000, and the size of the estimation set in the standard GWAS is chosen to match sampling variances between the two study designs. The polygenic score is based on 10,000 causal loci; its performance was evaluated in an independent set of 10,000 unrelated individuals. Standard-GWAS based polygenic scores outperforms (underperforms) sib-GWAS based polygenic scores under positive (negative) assortative mating. (B) Ratio of prediction accuracies of the polygenic scores based on sib- versus standard GWAS, as a function of $ρ_{a}$ , for two sets of simulations with one or two generations of assortative mating, with same parameters as in (A). (**C–F**) Simulation results, with the same parameters as in (A) but $ρ_{a} = 0.5$ , as a function of the number of SNPs included in the polygenic score, with all loci being causal (**C,D**), or with 20% of loci being causal (**E,F**). SNPs are added in the order of their association p-value in an independent set of 20,000 unrelated individuals. In both cases, the ratio of prediction accuracies for scores based on sib-GWAS versus standard GWAS becomes smaller with the inclusion of more weakly associated SNPs, a behavior that is qualitatively similar to observations in Figure 3 in the main text. Points are mean ± 2 SD in 10 simulation iterations. See Section 1.4.1 for simulation details.

Appendix 1—figure 11

Download asset Open asset

Comparison of prediction accuracies of polygenic scores based on standard and sib-GWAS matched for sex ratio.

This figure mirrors Figure 3B of the main text, but here the samples of siblings and unrelated individuals used in the analysis are matched for their sex ratios. Results are shown for diastolic blood pressure, as the prediction accuracy differed between sexes (Figure 1); the related phenotype of pulse rate; and a subset of the traits for which the prediction accuracy varied by GWAS design (Figure 3B). Small points show the ratio of the prediction accuracies in the two designs across 10 iterations; in each iteration, we resample sets of unrelated individuals to constitute three sets for discovery, estimation and prediction. Larger points show median values. We note that pulse rate is now similarly predicted by the two GWAS approaches, suggesting that perhaps the slightly higher prediction accuracy of the sib-GWAS shown in the main text Figure 3B are due to the sex ratio difference; for other traits, results are qualitatively unchanged.

Appendix 1—figure 12

Download asset Open asset

Prediction accuracy of polygenic scores based on sib-and standard GWAS, for a range of traits.

This figure complements Figure 3C–F of the main text, showing the results of the study design depicted in Figure 3A for all traits presented in Figure 3. As described for Figure 3, we randomly divided unrelated individuals to constitute three non-overlapping sets for discovery, estimation and prediction. Small points correspond to 10 iterations of resmapling these three sets. The prediction accuracy is plotted as a function of the p-value threshold, where p-values come from the discovery GWAS. Lines show median values.

Appendix 1—figure 13

Download asset Open asset

Prediction accuracy for years of schooling, for individuals with 0 or 1 full sibling.

(A) The y-axis shows the prediction accuracy, measured as incremental $R^{2}$ , in prediction sets stratified by participants’ number of siblings, using a polygenic score for years of schooling based on a GWAS performed using individuals who reported to have exactly 1 sibling. The x-axis shows the p-value threshold for inclusion of a SNP in the polygenic score when based on a pruning and thresholding approach. Last points on the x-axis correspond to a polygenic score model based on the LDpred approach (Vilhjálmsson et al., 2015) with a prior probability of 1 on loci being causal. Points are values based on 10 iterations of resampling estimation and prediction sets. Thick horizontal lines denote the mean values. (**B–E**) Comparison of the distribution of Townsend deprivation index (B) the age distribution (C), the proportion of males (D), and mean years of schooling (± 2 SD) between individuals who reported having no sibling and those who reported having 1 sibling. The two sets have somewhat different distributions of ages (or possibly come from somewhat different birth cohorts), a feature that could contribute to the patterns seen in panel A, but are otherwise similar with respect to the other features considered.

Appendix 1—figure 14

Download asset Open asset

Variable prediction accuracy for binary traits, when measured as incremental AUC.

This figure is analogous to the one shown in Figure 1 of the main text, but considering dichotomized versions of the traits presented in Figure 1 in the prediction sets, and with the y-axis showing incremental AUC values rather than incremental $R^{2}$ . The polygenic scores are based on GWAS using the quantitative trait values as in Figure 1. The traits are (A) diastolic blood pressure of over 110 mmHg, (B) BMI of over 35 Kg/m², and (C) completing a college or a university degree. Each box and whiskers plot was computed based on 20 iterations of resampling estimation and prediction sets. Thick horizontal lines denote the medians.

Appendix 1—figure 15

Download asset Open asset

Variable prediction accuracy for binary disease phenotypes, measured as incremental AUC, in men versus women.

This figure is analogous to the one shown in Figure 1 of the main text, but looking at disease traits, and with the y-axis showing incremental AUC rather than incremental $R^{2}$ . Each box and whiskers plot was computed based on 20 iterations of resampling estimation and prediction sets. Thick horizontal lines denote the medians. The variable prediction accuracy of PGS based on GWAS in men only versus women only could be driven in part by the differences in ratios of cases to controls (and hence by differences in the precision of the effect size estimates). However, we also observe that the prediction accuracy can vary depending on the sex composition of the prediction set (e.g., for cardiovascular outcomes), an observation that cannot be attributed to differences in case:control ratios of the GWAS.

Appendix 1—figure 16

Download asset Open asset

Comparison of prediction accuracies of polygenic scores (measured as $R^{2}$ ) based on standard and sib-GWAS.

This figure mirrors Figure 3B of the main text, but here we first residualized the phenotypes on covariates, and then ran the same pipeline described as that used to generate Figure 3B on the residuals without further adjustment for covariates in the GWAS or prediction evaluation. Thus, this figure relates more directly to the analytical derivation in Section 1.2. However, the results in Figure 3B and here are qualitatively similar. Small points show the ratio of the prediction accuracies in the two designs across 10 iterations; in each iteration, we resample sets of unrelated individuals to constitute three sets for discovery, estimation and prediction. Larger points show median values.

Appendix 1—table 1

UK Biobank phenotype data used in this study and their corresponding data fields.

In parentheses are the units of measurements.

Trait	Description	UKB data field
Age	Age when attended assessment center (years)	21003
Age at first sex	Self-reported age at first sexual intercourse (years)	2139
Alcohol intake frequency	Self-reported category, encoded as an integer: 1, 'Daily or almost daily’; 2, 'Three or four times a week’; 3, 'Once or twice a week’; 4, 'One to three times a month’; 5, 'Special occasions only’; 6, 'Never’	1558
Basal metabolic rate	Estimated from body composition impedance measurements (KJ)	23105
Birth weight	Self-reported birth weight (Kg)	20022
Body mass index	Constructed from height and weight measurements (Kg/m²)	21001
Diastolic blood pressure	Measured using automated devices (mmHg); values are adjusted for medicine use (see Materials and methods)	4079, 6153, 6177
Fluid intelligence	Unweighted sum of the number of correct answers given to 13 fluid intelligence questions	20016
Forced vital capacity	Calculated from breath spirometry (liters)	3062
Hair color	Self-reported category, encoded as an integer: 1, 'Blonde’; 2, 'Red’; 3, 'Light brown’; 4, 'Dark brown’; 5, 'Black’; none, 'Other’	1747
Hand grip strength	Measured right and left hand isometric grip strength (Kg)	46, 47
Height	Measured standing height (cm)	50
Hip circumference	Measured hip circumference (cm)	49
Hospital inpatient diagnoses	Diagnoses made during hospital inpatient admissions, coded according to the International Classification of Diseases (ICD-9 and ICD-10)	41202, 41203, 41204, 41205, 41270, 41271
Household income	Self-reported average total annual household income before tax category, encoded as an integer: 1, 'Less than £18,000'; 2, '£18,000 to £30,999'; 3, '£31,000 to £51,999'; 4, '£52,000 to £100,000'; 5, 'Greater than £100,000'	738
Myocardial infarction outcomes	Algorithmically-defined myocardial infarction outcomes obtained through combinations of UK Biobank's assessment data collection (e.g., self-reported conditions), and data from hospital admissions	42000, 42001
Neuroticism score	Derived summary score, based on participants’ responses to 12 neurotic behaviour-related questions	20127
Number of full siblings	Sum of self-reported number of full brothers and full sisters	1873, 1883
Overall health rating	Self-reported category, encoded as an integer: 1, 'Excellent'; 2, 'Good'; 3, 'Fair'; 4, 'Poor’	2178
Pack years of smoking	Calculated for individuals who have ever smoked as the number of cigarettes smoked per day, divided by twenty, multiplied by the number of years of smoking (years)	20161
Pulse rate	Measured during the automated blood pressure readings (bpm)	102
Qualifications	Self-reported educational or professional qualifications, selected from: 'College or University degree', 'NVQ or HND or HNC or equivalent', 'Other professional qualifications eg: nursing, teaching', 'A levels/AS levels or equivalent', 'O levels/GCSEs or equivalent', 'CSEs or equivalent', or 'None of the above'	6138
Sex	Self-reported sex and as determined from genotyping analysis	31, 22001
Skin color	Self-reported category, encoded as an integer: 1, 'Very fair'; 2, 'Fair'; 3, 'Light olive'; 4, 'Dark olive'; 5, 'Brown'; 6, 'Black'	1717
Townsend deprivation index	Townsend deprivation index at recruitment	189
Vascular/heart problems	Self-reported vascular/heart problems diagnosed by doctor selected from the categories: 'Heart attack’, 'Angina', 'Stroke', 'High blood pressure', and 'None of the above'	6150
Waist circumference	Measured waist circumference (cm)	48

Appendix 1—table 2

Genetic correlations across samples that vary by a study characteristic.

Numbers are genetic correlations estimated using LD score regression for BMI, years of schooling and diastolic blood pressure, across samples stratified by age, Townsend deprivation index (a measure of socioeconomic status, SES), and sex, respectively. ’Q’ denotes quartile of age or SES.

Trait/characteristic	Pair of strata	Genetic correlation (s.e.)
BMI/Age	(Q1,Q2)	0.93 (0.036)
	(Q1,Q3)	0.95 (0.035)
	(Q1,Q4)	0.95 (0.038)
	(Q2,Q3)	0.89 (0.032)
	(Q2,Q4)	0.91 (0.036)
	(Q3,Q4)	1.00 (0.040)
Years of schooling/SES	(Q1,Q2)	0.98 (0.054)
	(Q1,Q3)	1.00 (0.067)
	(Q1,Q4)	0.93 (0.068)
	(Q2,Q3)	0.97 (0.064)
	(Q2,Q4)	1.09 (0.074)
	(Q3,Q4)	1.04 (0.074)
Diastolic blood pressure/Sex	(male,female)	0.93 (0.031)

Appendix 1—table 3

Sample sizes used for siblings and unrelated sets.

As described in Figure 3A, for the comparison of prediction accuracies of polygenic scores based on standard and sib-GWAS, we first ascertain SNPs in a large sample of unrelated individuals (‘Unrelated-discovery’) and then estimate the effect of these SNPs with a standard regression using unrelated individuals (‘Unrelated-n*') and, independently, using sib-regression (in the ‘Siblings’ set). Finally, we used the polygenic scores for prediction in a third sample of unrelated individuals (‘Unrelated-prediction’). This table shows sample sizes used for each set across the traits analyzed. For simulated traits, the numbers in parentheses are heritability, number of causal loci, and simulation replicate, respectively (three traits were simulated for each pair of heritability and number of causal loci parameters, see Materials and methods for simulation details).

Trait	Siblings (pairs)	Unrelated-n*	Unrelated- discovery	Unrelated- prediction
Age at first sex	13675	8746	244988	27220
Alcohol intake frequency	17282	10923	276885	30764
Basal metabolic rate	16802	13467	269750	29972
Birth weight	6750	5766	159074	17674
BMI	17217	12359	274868	30540
Diastolic blood pressure	14791	9514	253227	28136
Fluid intelligence	3889	2979	101016	11223
Forced vital capacity	14605	10009	252576	28064
Hair color	16859	11763	272209	30245
Hand grip strength	17070	10832	275117	30568
Height	17242	18147	269973	29997
Hip circumference	17254	11648	275930	30658
Household income	13240	8704	239326	26591
Neuroticism score	11756	6909	227010	25223
Overall health rating	17189	10378	276581	30731
Pack years of smoking	2307	1682	85544	9504
Pulse rate	14791	8812	253859	28206
Skin color	16903	10334	274159	30462
Waist circumference	17257	11749	275873	30652
Years of schooling	17037	11885	273553	30394
Simulated trait (0.5,10K,1)	17299	11685	276404	30711
Simulated trait (0.5,10K,2)	17299	11505	276566	30729
Simulated trait (0.5,10K,3)	17299	11422	276641	30737
Simulated trait (0.5,100K,1)	17299	11814	276288	30698
Simulated trait (0.5,100K,2)	17299	11833	276271	30696
Simulated trait (0.5,100K,3)	17299	11490	276579	30731
Simulated trait (0.1,10K,1)	17299	9072	278756	30972
Simulated trait (0.1,10K,2)	17299	9158	278678	30964
Simulated trait (0.1,10K,3)	17299	9111	278721	30968
Simulated trait (0.1,100K,1)	17299	9133	278701	30966
Simulated trait (0.1,100K,2)	17299	9069	278758	30973
Simulated trait (0.1,100K,3)	17299	9108	278723	30969

Appendix 1—table 4

Qualifications to years of schooling conversion table.

Educational or professional qualifications were converted to years of schooling following Okbay et al. (2016).

Qualifications (UKB data field 6138)	Years of schooling
College or University degree	20
NVQ or HND or HNC or equivalent	19
Other professional qualifications eg: nursing, teaching	15
A levels/AS levels or equivalent	13
O levels/GCSEs or equivalent	10
CSEs or equivalent	10
None of the above	7

Data availability

The GWAS summary statistics generated in this study have been uploaded to Dryad.

The following data sets were generated

(2019) Dryad Digital Repository
Variable prediction accuracy of polygenic scores within an ancestry group.

https://doi.org/10.5061/dryad.66t1g1jxs

References

(2019) A GWAS in latin americans highlights the convergent evolution of lighter skin pigmentation in eurasia
Nature Communications 10:385.

https://doi.org/10.1038/s41467-018-08147-0
- Google Scholar
Preprint
(2018) Education can reduce health disparities related to genetic risk of obesity: evidence from a british reform
bioRxiv.

https://doi.org/10.1101/260463
- Google Scholar
1. Belsky DW
2. Domingue BW
3. Wedow R
4. Arseneault L
5. Boardman JD
6. Caspi A
7. Conley D
8. Fletcher JM
9. Freese J
10. Herd P
11. Moffitt TE
12. Poulton R
13. Sicinski K
14. Wertz J
15. Harris KM
(2018) Genetic analysis of social-class mobility in five longitudinal studies
PNAS 115:E7275–E7284.

https://doi.org/10.1073/pnas.1801238115
- PubMed
- Google Scholar
1. Bentley AR
2. Sung YJ
3. Brown MR
4. Winkler TW
5. Kraja AT
6. Ntalla I
7. Schwander K
8. Chasman DI
9. Lim E
10. Deng X
11. Guo X
12. Liu J
13. Lu Y
14. Cheng CY
15. Sim X
16. Vojinovic D
17. Huffman JE
18. Musani SK
19. Li C
20. Feitosa MF
21. Richard MA
22. Noordam R
23. Baker J
24. Chen G
25. Aschard H
26. Bartz TM
27. Ding J
28. Dorajoo R
29. Manning AK
30. Rankinen T
31. Smith AV
32. Tajuddin SM
33. Zhao W
34. Graff M
35. Alver M
36. Boissel M
37. Chai JF
38. Chen X
39. Divers J
40. Evangelou E
41. Gao C
42. Goel A
43. Hagemeijer Y
44. Harris SE
45. Hartwig FP
46. He M
47. Horimoto A
48. Hsu FC
49. Hung YJ
50. Jackson AU
51. Kasturiratne A
52. Komulainen P
53. Kühnel B
54. Leander K
55. Lin KH
56. Luan J
57. Lyytikäinen LP
58. Matoba N
59. Nolte IM
60. Pietzner M
61. Prins B
62. Riaz M
63. Robino A
64. Said MA
65. Schupf N
66. Scott RA
67. Sofer T
68. Stancáková A
69. Takeuchi F
70. Tayo BO
71. van der Most PJ
72. Varga TV
73. Wang TD
74. Wang Y
75. Ware EB
76. Wen W
77. Xiang YB
78. Yanek LR
79. Zhang W
80. Zhao JH
81. Adeyemo A
82. Afaq S
83. Amin N
84. Amini M
85. Arking DE
86. Arzumanyan Z
87. Aung T
88. Ballantyne C
89. Barr RG
90. Bielak LF
91. Boerwinkle E
92. Bottinger EP
93. Broeckel U
94. Brown M
95. Cade BE
96. Campbell A
97. Canouil M
98. Charumathi S
99. Chen YI
100. Christensen K
101. Concas MP
102. Connell JM
103. de Las Fuentes L
104. de Silva HJ
105. de Vries PS
106. Doumatey A
107. Duan Q
108. Eaton CB
109. Eppinga RN
110. Faul JD
111. Floyd JS
112. Forouhi NG
113. Forrester T
114. Friedlander Y
115. Gandin I
116. Gao H
117. Ghanbari M
118. Gharib SA
119. Gigante B
120. Giulianini F
121. Grabe HJ
122. Gu CC
123. Harris TB
124. Heikkinen S
125. Heng CK
126. Hirata M
127. Hixson JE
128. Ikram MA
129. Jia Y
130. Joehanes R
131. Johnson C
132. Jonas JB
133. Justice AE
134. Katsuya T
135. Khor CC
136. Kilpeläinen TO
137. Koh WP
138. Kolcic I
139. Kooperberg C
140. Krieger JE
141. Kritchevsky SB
142. Kubo M
143. Kuusisto J
144. Lakka TA
145. Langefeld CD
146. Langenberg C
147. Launer LJ
148. Lehne B
149. Lewis CE
150. Li Y
151. Liang J
152. Lin S
153. Liu CT
154. Liu J
155. Liu K
156. Loh M
157. Lohman KK
158. Louie T
159. Luzzi A
160. Mägi R
161. Mahajan A
162. Manichaikul AW
163. McKenzie CA
164. Meitinger T
165. Metspalu A
166. Milaneschi Y
167. Milani L
168. Mohlke KL
169. Momozawa Y
170. Morris AP
171. Murray AD
172. Nalls MA
173. Nauck M
174. Nelson CP
175. North KE
176. O'Connell JR
177. Palmer ND
178. Papanicolau GJ
179. Pedersen NL
180. Peters A
181. Peyser PA
182. Polasek O
183. Poulter N
184. Raitakari OT
185. Reiner AP
186. Renström F
187. Rice TK
188. Rich SS
189. Robinson JG
190. Rose LM
191. Rosendaal FR
192. Rudan I
193. Schmidt CO
194. Schreiner PJ
195. Scott WR
196. Sever P
197. Shi Y
198. Sidney S
199. Sims M
200. Smith JA
201. Snieder H
202. Starr JM
203. Strauch K
204. Stringham HM
205. Tan NYQ
206. Tang H
207. Taylor KD
208. Teo YY
209. Tham YC
210. Tiemeier H
211. Turner ST
212. Uitterlinden AG
213. van Heemst D
214. Waldenberger M
215. Wang H
216. Wang L
217. Wang L
218. Wei WB
219. Williams CA
220. Wilson G
221. Wojczynski MK
222. Yao J
223. Young K
224. Yu C
225. Yuan JM
226. Zhou J
227. Zonderman AB
228. Becker DM
229. Boehnke M
230. Bowden DW
231. Chambers JC
232. Cooper RS
233. de Faire U
234. Deary IJ
235. Elliott P
236. Esko T
237. Farrall M
238. Franks PW
239. Freedman BI
240. Froguel P
241. Gasparini P
242. Gieger C
243. Horta BL
244. Juang JJ
245. Kamatani Y
246. Kammerer CM
247. Kato N
248. Kooner JS
249. Laakso M
250. Laurie CC
251. Lee IT
252. Lehtimäki T
253. Magnusson PKE
254. Oldehinkel AJ
255. Penninx B
256. Pereira AC
257. Rauramaa R
258. Redline S
259. Samani NJ
260. Scott J
261. Shu XO
262. van der Harst P
263. Wagenknecht LE
264. Wang JS
265. Wang YX
266. Wareham NJ
267. Watkins H
268. Weir DR
269. Wickremasinghe AR
270. Wu T
271. Zeggini E
272. Zheng W
273. Bouchard C
274. Evans MK
275. Gudnason V
276. Kardia SLR
277. Liu Y
278. Psaty BM
279. Ridker PM
280. van Dam RM
281. Mook-Kanamori DO
282. Fornage M
283. Province MA
284. Kelly TN
285. Fox ER
286. Hayward C
287. van Duijn CM
288. Tai ES
289. Wong TY
290. Loos RJF
291. Franceschini N
292. Rotter JI
293. Zhu X
294. Bierut LJ
295. Gauderman WJ
296. Rice K
297. Munroe PB
298. Morrison AC
299. Rao DC
300. Rotimi CN
301. Cupples LA
302. COGENT-Kidney Consortium, EPIC-InterAct Consortium, Understanding Society Scientific Group, Lifelines Cohort
(2019) Multi-ancestry genome-wide gene-smoking interaction study of 387,272 individuals identifies new loci associated with serum lipids
Nature Genetics 51:636–648.

https://doi.org/10.1038/s41588-019-0378-y
- PubMed
- Google Scholar
1. Berg JJ
2. Harpak A
3. Sinnott-Armstrong N
4. Joergensen AM
5. Mostafavi H
6. Field Y
7. Boyle EA
8. Zhang X
9. Racimo F
10. Pritchard JK
11. Coop G
(2019) Reduced signal for polygenic adaptation of height in UK biobank
eLife 8:e39725.

https://doi.org/10.7554/eLife.39725
- PubMed
- Google Scholar
1. Berg JJ
2. Coop G
(2014) A population genetic signal of polygenic adaptation
PLOS Genetics 10:e1004412.

https://doi.org/10.1371/journal.pgen.1004412
- PubMed
- Google Scholar
1. Bien SA
2. Wojcik GL
3. Hodonsky CJ
4. Gignoux CR
5. Cheng I
6. Matise TC
7. Peters U
8. Kenny EE
9. North KE
(2019) The future of genomic studies must be globally representative: perspectives from PAGE
Annual Review of Genomics and Human Genetics 20:181–200.

https://doi.org/10.1146/annurev-genom-091416-035517
- PubMed
- Google Scholar
(2017) An expanded view of complex traits: from polygenic to omnigenic
Cell 169:1177–1186.

https://doi.org/10.1016/j.cell.2017.05.038
- PubMed
- Google Scholar
(2013) Variation in the heritability of educational attainment: an international Meta-Analysis
Social Forces 92:109–140.

https://doi.org/10.1093/sf/sot076
- Google Scholar
1. Briley DA
2. Tucker-Drob EM
(2013) Explaining the increasing heritability of cognitive ability across development: a meta-analysis of longitudinal twin and adoption studies
Psychological Science 24:1704–1713.

https://doi.org/10.1177/0956797613478618
- PubMed
- Google Scholar
(2015) LD score regression distinguishes confounding from polygenicity in genome-wide association studies
Nature Genetics 47:291–295.

https://doi.org/10.1038/ng.3211
- PubMed
- Google Scholar
1. Bycroft C
2. Freeman C
3. Petkova D
4. Band G
5. Elliott LT
6. Sharp K
7. Motyer A
8. Vukcevic D
9. Delaneau O
10. O'Connell J
11. Cortes A
12. Welsh S
13. Young A
14. Effingham M
15. McVean G
16. Leslie S
17. Allen N
18. Donnelly P
19. Marchini J
(2018) The UK biobank resource with deep phenotyping and genomic data
Nature 562:203–209.

https://doi.org/10.1038/s41586-018-0579-z
- PubMed
- Google Scholar
1. Chang CC
2. Chow CC
3. Tellier LC
4. Vattikuti S
5. Purcell SM
6. Lee JJ
(2015) Second-generation PLINK: rising to the challenge of larger and richer datasets
GigaScience 4:7.

https://doi.org/10.1186/s13742-015-0047-8
- PubMed
- Google Scholar
Report
1. Chetty R
2. Hendren N
(2018)
Race and Economic Opportunity in the United States: An Intergenerational Perspective

National Bureau of Economic Research.
- Google Scholar
Book
1. Conley D
(2010)
Being Black, Living in the Red: Race, Wealth, and Social Policy in America

University of California Press.
- Google Scholar
(2015) Is the effect of parental education on offspring biased or moderated by genotype?
Sociological Science 2:82–105.

https://doi.org/10.15195/v2.a6
- PubMed
- Google Scholar
1. Conley D
(2016) Socio-Genomic research using Genome-Wide molecular data
Annual Review of Sociology 42:275–299.

https://doi.org/10.1146/annurev-soc-081715-074316
- Google Scholar
(2018) The causal effects of education on health outcomes in the UK biobank
Nature Human Behaviour 2:117–125.

https://doi.org/10.1038/s41562-017-0279-y
- PubMed
- Google Scholar
1. De La Vega FM
2. Bustamante CD
(2018) Polygenic risk scores: a biased prediction?
Genome Medicine 10:100.

https://doi.org/10.1186/s13073-018-0610-x
- PubMed
- Google Scholar
(2014) Genetic and educational assortative mating among US adults
PNAS 111:7996–8000.

https://doi.org/10.1073/pnas.1321426111
- PubMed
- Google Scholar
(2018) The social genome of friends and schoolmates in the national longitudinal study of adolescent to adult health
PNAS 115:702–707.

https://doi.org/10.1073/pnas.1711803115
- Google Scholar
1. Dudbridge F
(2013) Power and predictive accuracy of polygenic risk scores
PLOS Genetics 9:e1003348.

https://doi.org/10.1371/journal.pgen.1003348
- PubMed
- Google Scholar
Preprint
1. Duncan L
2. Shen H
3. Gelaye B
4. Ressler K
5. Feldman M
6. Peterson R
7. Domingue B
(2018) Analysis of polygenic score usage and performance across diverse human populations
bioRxiv.

https://doi.org/10.1101/398396
- Google Scholar
1. Edge MD
2. Coop G
(2019) Reconstructing the history of polygenic scores using coalescent trees
Genetics 211:235–262.

https://doi.org/10.1534/genetics.118.301687
- PubMed
- Google Scholar
1. Elks CE
2. den Hoed M
3. Zhao JH
4. Sharp SJ
5. Wareham NJ
6. Loos RJ
7. Ong KK
(2012) Variability in the heritability of body mass index: a systematic review and meta-regression
Frontiers in Endocrinology 3:29.

https://doi.org/10.3389/fendo.2012.00029
- PubMed
- Google Scholar
(2015) PRSice: polygenic risk score software
Bioinformatics 31:1466–1468.

https://doi.org/10.1093/bioinformatics/btu848
- PubMed
- Google Scholar
1. Evangelou E
2. Warren HR
3. Mosen-Ansorena D
4. Mifsud B
5. Pazoki R
6. Gao H
7. Ntritsos G
8. Dimou N
9. Cabrera CP
10. Karaman I
11. Ng FL
12. Evangelou M
13. Witkowska K
14. Tzanis E
15. Hellwege JN
16. Giri A
17. Velez Edwards DR
18. Sun YV
19. Cho K
20. Gaziano JM
21. Wilson PWF
22. Tsao PS
23. Kovesdy CP
24. Esko T
25. Mägi R
26. Milani L
27. Almgren P
28. Boutin T
29. Debette S
30. Ding J
31. Giulianini F
32. Holliday EG
33. Jackson AU
34. Li-Gao R
35. Lin WY
36. Luan J
37. Mangino M
38. Oldmeadow C
39. Prins BP
40. Qian Y
41. Sargurupremraj M
42. Shah N
43. Surendran P
44. Thériault S
45. Verweij N
46. Willems SM
47. Zhao JH
48. Amouyel P
49. Connell J
50. de Mutsert R
51. Doney ASF
52. Farrall M
53. Menni C
54. Morris AD
55. Noordam R
56. Paré G
57. Poulter NR
58. Shields DC
59. Stanton A
60. Thom S
61. Abecasis G
62. Amin N
63. Arking DE
64. Ayers KL
65. Barbieri CM
66. Batini C
67. Bis JC
68. Blake T
69. Bochud M
70. Boehnke M
71. Boerwinkle E
72. Boomsma DI
73. Bottinger EP
74. Braund PS
75. Brumat M
76. Campbell A
77. Campbell H
78. Chakravarti A
79. Chambers JC
80. Chauhan G
81. Ciullo M
82. Cocca M
83. Collins F
84. Cordell HJ
85. Davies G
86. de Borst MH
87. de Geus EJ
88. Deary IJ
89. Deelen J
90. Del Greco M F
91. Demirkale CY
92. Dörr M
93. Ehret GB
94. Elosua R
95. Enroth S
96. Erzurumluoglu AM
97. Ferreira T
98. Frånberg M
99. Franco OH
100. Gandin I
101. Gasparini P
102. Giedraitis V
103. Gieger C
104. Girotto G
105. Goel A
106. Gow AJ
107. Gudnason V
108. Guo X
109. Gyllensten U
110. Hamsten A
111. Harris TB
112. Harris SE
113. Hartman CA
114. Havulinna AS
115. Hicks AA
116. Hofer E
117. Hofman A
118. Hottenga JJ
119. Huffman JE
120. Hwang SJ
121. Ingelsson E
122. James A
123. Jansen R
124. Jarvelin MR
125. Joehanes R
126. Johansson Å
127. Johnson AD
128. Joshi PK
129. Jousilahti P
130. Jukema JW
131. Jula A
132. Kähönen M
133. Kathiresan S
134. Keavney BD
135. Khaw KT
136. Knekt P
137. Knight J
138. Kolcic I
139. Kooner JS
140. Koskinen S
141. Kristiansson K
142. Kutalik Z
143. Laan M
144. Larson M
145. Launer LJ
146. Lehne B
147. Lehtimäki T
148. Liewald DCM
149. Lin L
150. Lind L
151. Lindgren CM
152. Liu Y
153. Loos RJF
154. Lopez LM
155. Lu Y
156. Lyytikäinen LP
157. Mahajan A
158. Mamasoula C
159. Marrugat J
160. Marten J
161. Milaneschi Y
162. Morgan A
163. Morris AP
164. Morrison AC
165. Munson PJ
166. Nalls MA
167. Nandakumar P
168. Nelson CP
169. Niiranen T
170. Nolte IM
171. Nutile T
172. Oldehinkel AJ
173. Oostra BA
174. O'Reilly PF
175. Org E
176. Padmanabhan S
177. Palmas W
178. Palotie A
179. Pattie A
180. Penninx B
181. Perola M
182. Peters A
183. Polasek O
184. Pramstaller PP
185. Nguyen QT
186. Raitakari OT
187. Ren M
188. Rettig R
189. Rice K
190. Ridker PM
191. Ried JS
192. Riese H
193. Ripatti S
194. Robino A
195. Rose LM
196. Rotter JI
197. Rudan I
198. Ruggiero D
199. Saba Y
200. Sala CF
201. Salomaa V
202. Samani NJ
203. Sarin AP
204. Schmidt R
205. Schmidt H
206. Shrine N
207. Siscovick D
208. Smith AV
209. Snieder H
210. Sõber S
211. Sorice R
212. Starr JM
213. Stott DJ
214. Strachan DP
215. Strawbridge RJ
216. Sundström J
217. Swertz MA
218. Taylor KD
219. Teumer A
220. Tobin MD
221. Tomaszewski M
222. Toniolo D
223. Traglia M
224. Trompet S
225. Tuomilehto J
226. Tzourio C
227. Uitterlinden AG
228. Vaez A
229. van der Most PJ
230. van Duijn CM
231. Vergnaud AC
232. Verwoert GC
233. Vitart V
234. Völker U
235. Vollenweider P
236. Vuckovic D
237. Watkins H
238. Wild SH
239. Willemsen G
240. Wilson JF
241. Wright AF
242. Yao J
243. Zemunik T
244. Zhang W
245. Attia JR
246. Butterworth AS
247. Chasman DI
248. Conen D
249. Cucca F
250. Danesh J
251. Hayward C
252. Howson JMM
253. Laakso M
254. Lakatta EG
255. Langenberg C
256. Melander O
257. Mook-Kanamori DO
258. Palmer CNA
259. Risch L
260. Scott RA
261. Scott RJ
262. Sever P
263. Spector TD
264. van der Harst P
265. Wareham NJ
266. Zeggini E
267. Levy D
268. Munroe PB
269. Newton-Cheh C
270. Brown MJ
271. Metspalu A
272. Hung AM
273. O'Donnell CJ
274. Edwards TL
275. Psaty BM
276. Tzoulaki I
277. Barnes MR
278. Wain LV
279. Elliott P
280. Caulfield MJ
281. Million Veteran Program
(2018) Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits
Nature Genetics 50:1412–1425.

https://doi.org/10.1038/s41588-018-0205-x
- PubMed
- Google Scholar
1. Field Y
2. Boyle EA
3. Telis N
4. Gao Z
5. Gaulton KJ
6. Golan D
7. Yengo L
8. Rocheleau G
9. Froguel P
10. McCarthy MI
11. Pritchard JK
(2016) Detection of human adaptation during the past 2000 years
Science 354:760–764.

https://doi.org/10.1126/science.aag0776
- PubMed
- Google Scholar
1. Fry A
2. Littlejohns TJ
3. Sudlow C
4. Doherty N
5. Adamska L
6. Sprosen T
7. Collins R
8. Allen NE
(2017) Comparison of Sociodemographic and Health-Related characteristics of UK biobank participants with those of the general population
American Journal of Epidemiology 186:1026–1034.

https://doi.org/10.1093/aje/kwx246
- PubMed
- Google Scholar
1. Ge T
2. Chen CY
3. Neale BM
4. Sabuncu MR
5. Smoller JW
(2017) Phenome-wide heritability analysis of the UK biobank
PLOS Genetics 13:e1006711.

https://doi.org/10.1371/journal.pgen.1006711
- PubMed
- Google Scholar
1. Gibson G
(2008) The environmental contribution to gene expression profiles
Nature Reviews Genetics 9:575–581.

https://doi.org/10.1038/nrg2383
- PubMed
- Google Scholar
1. Haworth S
2. Mitchell R
3. Corbin L
4. Wade KH
5. Dudding T
6. Budu-Aggrey A
7. Carslake D
8. Hemani G
9. Paternoster L
10. Smith GD
11. Davies N
12. Lawson DJ
13. J Timpson N
(2019) Apparent latent structure within the UK biobank sample has implications for epidemiological analysis
Nature Communications 10:333.

https://doi.org/10.1038/s41467-018-08219-1
- PubMed
- Google Scholar
Book
1. Henderson CR
(1984)
Applications of Linear Models in Animal Breeding, 462

University of Guelph Guelph.
- Google Scholar
(2019) Polygenic adaptation: From sweeps to subtle frequency shifts
PLOS Genetics 15:e1008035.

https://doi.org/10.1371/journal.pgen.1008035
- Google Scholar
(2018) Genomic risk prediction of coronary artery disease in 480,000 adults: implications for primary prevention
Journal of the American College of Cardiology 72:1883–1893.

https://doi.org/10.1016/j.jacc.2018.07.079
- PubMed
- Google Scholar
1. Kathiresan S
2. Melander O
3. Anevski D
4. Guiducci C
5. Burtt NP
6. Roos C
7. Hirschhorn JN
8. Berglund G
9. Hedblad B
10. Groop L
11. Altshuler DM
12. Newton-Cheh C
13. Orho-Melander M
(2008) Polymorphisms associated with cholesterol and risk of cardiovascular events
New England Journal of Medicine 358:1240–1249.

https://doi.org/10.1056/NEJMoa0706728
- PubMed
- Google Scholar
Preprint
(2018) Geographic variation and Bias in Polygenic scores of complex diseases and traits in Finland
bioRxiv.

https://doi.org/10.1101/485441
- Google Scholar
1. Khera AV
2. Chaffin M
3. Aragam KG
4. Haas ME
5. Roselli C
6. Choi SH
7. Natarajan P
8. Lander ES
9. Lubitz SA
10. Ellinor PT
11. Kathiresan S
(2018) Genome-wide polygenic scores for common diseases identify individuals with risk equivalent to monogenic mutations
Nature Genetics 50:1219–1224.

https://doi.org/10.1038/s41588-018-0183-z
- Google Scholar
1. Khera AV
2. Chaffin M
3. Wade KH
4. Zahid S
5. Brancale J
6. Xia R
7. Distefano M
8. Senol-Cosar O
9. Haas ME
10. Bick A
11. Aragam KG
12. Lander ES
13. Smith GD
14. Mason-Suares H
15. Fornage M
16. Lebo M
17. Timpson NJ
18. Kaplan LM
19. Kathiresan S
(2019) Polygenic prediction of weight and obesity trajectories from birth to adulthood
Cell 177:587–596.

https://doi.org/10.1016/j.cell.2019.03.028
- PubMed
- Google Scholar
1. Kim MS
2. Patel KP
3. Teng AK
4. Berens AJ
5. Lachance J
(2018) Genetic disease risks can be misestimated across global populations
Genome Biology 19:179.

https://doi.org/10.1186/s13059-018-1561-7
- PubMed
- Google Scholar
(2018) The nature of nurture: effects of parental genotypes
Science 359:424–428.

https://doi.org/10.1126/science.aan6877
- PubMed
- Google Scholar
1. Lawson DJ
2. Davies NM
3. Haworth S
4. Ashraf B
5. Howe L
6. Crawford A
7. Hemani G
8. Davey Smith G
9. Timpson NJ
(2020) Is population structure in the genetic biobank era irrelevant, a challenge, or an opportunity?
Human Genetics 139:1–2.

https://doi.org/10.1007/s00439-019-02014-8
- PubMed
- Google Scholar
1. Lee JJ
2. Wedow R
3. Okbay A
4. Kong E
5. Maghzian O
6. Zacher M
7. Nguyen-Viet TA
8. Bowers P
9. Sidorenko J
10. Karlsson Linnér R
11. Fontana MA
12. Kundu T
13. Lee C
14. Li H
15. Li R
16. Royer R
17. Timshel PN
18. Walters RK
19. Willoughby EA
20. Yengo L
21. Alver M
22. Bao Y
23. Clark DW
24. Day FR
25. Furlotte NA
26. Joshi PK
27. Kemper KE
28. Kleinman A
29. Langenberg C
30. Mägi R
31. Trampush JW
32. Verma SS
33. Wu Y
34. Lam M
35. Zhao JH
36. Zheng Z
37. Boardman JD
38. Campbell H
39. Freese J
40. Harris KM
41. Hayward C
42. Herd P
43. Kumari M
44. Lencz T
45. Luan J
46. Malhotra AK
47. Metspalu A
48. Milani L
49. Ong KK
50. Perry JRB
51. Porteous DJ
52. Ritchie MD
53. Smart MC
54. Smith BH
55. Tung JY
56. Wareham NJ
57. Wilson JF
58. Beauchamp JP
59. Conley DC
60. Esko T
61. Lehrer SF
62. Magnusson PKE
63. Oskarsson S
64. Pers TH
65. Robinson MR
66. Thom K
67. Watson C
68. Chabris CF
69. Meyer MN
70. Laibson DI
71. Yang J
72. Johannesson M
73. Koellinger PD
74. Turley P
75. Visscher PM
76. Benjamin DJ
77. Cesarini D
78. 23andMe Research Team, COGENT (Cognitive Genomics Consortium), Social Science Genetic Association Consortium
(2018) Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals
Nature Genetics 50:1112–1121.

https://doi.org/10.1038/s41588-018-0147-3
- PubMed
- Google Scholar
(2013) FaST-LMM-Select for addressing confounding from spatial structure and rare variants
Nature Genetics 45:470–471.

https://doi.org/10.1038/ng.2620
- PubMed
- Google Scholar
1. Loh PR
2. Tucker G
3. Bulik-Sullivan BK
4. Vilhjálmsson BJ
5. Finucane HK
6. Salem RM
7. Chasman DI
8. Ridker PM
9. Neale BM
10. Berger B
11. Patterson N
12. Price AL
(2015) Efficient bayesian mixed-model analysis increases association power in large cohorts
Nature Genetics 47:284–290.

https://doi.org/10.1038/ng.3190
- PubMed
- Google Scholar
Book
1. Lynch M
2. Walsh B
(1998)
Genetics and Analysis of Quantitative Traits, 1

Sunderland, MA: Sinauer.
- Google Scholar
1. Martin AR
2. Gignoux CR
3. Walters RK
4. Wojcik GL
5. Neale BM
6. Gravel S
7. Daly MJ
8. Bustamante CD
9. Kenny EE
(2017) Human demographic history impacts genetic risk prediction across diverse populations
The American Journal of Human Genetics 100:635–649.

https://doi.org/10.1016/j.ajhg.2017.03.004
- PubMed
- Google Scholar
1. Martin AR
2. Teferra S
3. Möller M
4. Hoal EG
5. Daly MJ
(2018) The critical needs and challenges for genetic architecture studies in Africa
Current Opinion in Genetics & Development 53:113–120.

https://doi.org/10.1016/j.gde.2018.08.005
- PubMed
- Google Scholar
1. Martin AR
2. Kanai M
3. Kamatani Y
4. Okada Y
5. Neale BM
6. Daly MJ
(2019) Clinical use of current polygenic risk scores may exacerbate health disparities
Nature Genetics 51:584–591.

https://doi.org/10.1038/s41588-019-0379-x
- PubMed
- Google Scholar
1. Mathieson I
2. McVean G
(2013) Reply to: "FaST-LMM-Select for addressing confounding from spatial structure and rare variants"
Nature Genetics 45:471.

https://doi.org/10.1038/ng.2619
- PubMed
- Google Scholar
1. Mavaddat N
2. Michailidou K
3. Dennis J
4. Lush M
5. Fachal L
6. Lee A
7. Tyrer JP
8. Chen TH
9. Wang Q
10. Bolla MK
11. Yang X
12. Adank MA
13. Ahearn T
14. Aittomäki K
15. Allen J
16. Andrulis IL
17. Anton-Culver H
18. Antonenkova NN
19. Arndt V
20. Aronson KJ
21. Auer PL
22. Auvinen P
23. Barrdahl M
24. Beane Freeman LE
25. Beckmann MW
26. Behrens S
27. Benitez J
28. Bermisheva M
29. Bernstein L
30. Blomqvist C
31. Bogdanova NV
32. Bojesen SE
33. Bonanni B
34. Børresen-Dale AL
35. Brauch H
36. Bremer M
37. Brenner H
38. Brentnall A
39. Brock IW
40. Brooks-Wilson A
41. Brucker SY
42. Brüning T
43. Burwinkel B
44. Campa D
45. Carter BD
46. Castelao JE
47. Chanock SJ
48. Chlebowski R
49. Christiansen H
50. Clarke CL
51. Collée JM
52. Cordina-Duverger E
53. Cornelissen S
54. Couch FJ
55. Cox A
56. Cross SS
57. Czene K
58. Daly MB
59. Devilee P
60. Dörk T
61. Dos-Santos-Silva I
62. Dumont M
63. Durcan L
64. Dwek M
65. Eccles DM
66. Ekici AB
67. Eliassen AH
68. Ellberg C
69. Engel C
70. Eriksson M
71. Evans DG
72. Fasching PA
73. Figueroa J
74. Fletcher O
75. Flyger H
76. Försti A
77. Fritschi L
78. Gabrielson M
79. Gago-Dominguez M
80. Gapstur SM
81. García-Sáenz JA
82. Gaudet MM
83. Georgoulias V
84. Giles GG
85. Gilyazova IR
86. Glendon G
87. Goldberg MS
88. Goldgar DE
89. González-Neira A
90. Grenaker Alnæs GI
91. Grip M
92. Gronwald J
93. Grundy A
94. Guénel P
95. Haeberle L
96. Hahnen E
97. Haiman CA
98. Håkansson N
99. Hamann U
100. Hankinson SE
101. Harkness EF
102. Hart SN
103. He W
104. Hein A
105. Heyworth J
106. Hillemanns P
107. Hollestelle A
108. Hooning MJ
109. Hoover RN
110. Hopper JL
111. Howell A
112. Huang G
113. Humphreys K
114. Hunter DJ
115. Jakimovska M
116. Jakubowska A
117. Janni W
118. John EM
119. Johnson N
120. Jones ME
121. Jukkola-Vuorinen A
122. Jung A
123. Kaaks R
124. Kaczmarek K
125. Kataja V
126. Keeman R
127. Kerin MJ
128. Khusnutdinova E
129. Kiiski JI
130. Knight JA
131. Ko YD
132. Kosma VM
133. Koutros S
134. Kristensen VN
135. Krüger U
136. Kühl T
137. Lambrechts D
138. Le Marchand L
139. Lee E
140. Lejbkowicz F
141. Lilyquist J
142. Lindblom A
143. Lindström S
144. Lissowska J
145. Lo WY
146. Loibl S
147. Long J
148. Lubiński J
149. Lux MP
150. MacInnis RJ
151. Maishman T
152. Makalic E
153. Maleva Kostovska I
154. Mannermaa A
155. Manoukian S
156. Margolin S
157. Martens JWM
158. Martinez ME
159. Mavroudis D
160. McLean C
161. Meindl A
162. Menon U
163. Middha P
164. Miller N
165. Moreno F
166. Mulligan AM
167. Mulot C
168. Muñoz-Garzon VM
169. Neuhausen SL
170. Nevanlinna H
171. Neven P
172. Newman WG
173. Nielsen SF
174. Nordestgaard BG
175. Norman A
176. Offit K
177. Olson JE
178. Olsson H
179. Orr N
180. Pankratz VS
181. Park-Simon TW
182. Perez JIA
183. Pérez-Barrios C
184. Peterlongo P
185. Peto J
186. Pinchev M
187. Plaseska-Karanfilska D
188. Polley EC
189. Prentice R
190. Presneau N
191. Prokofyeva D
192. Purrington K
193. Pylkäs K
194. Rack B
195. Radice P
196. Rau-Murthy R
197. Rennert G
198. Rennert HS
199. Rhenius V
200. Robson M
201. Romero A
202. Ruddy KJ
203. Ruebner M
204. Saloustros E
205. Sandler DP
206. Sawyer EJ
207. Schmidt DF
208. Schmutzler RK
209. Schneeweiss A
210. Schoemaker MJ
211. Schumacher F
212. Schürmann P
213. Schwentner L
214. Scott C
215. Scott RJ
216. Seynaeve C
217. Shah M
218. Sherman ME
219. Shrubsole MJ
220. Shu XO
221. Slager S
222. Smeets A
223. Sohn C
224. Soucy P
225. Southey MC
226. Spinelli JJ
227. Stegmaier C
228. Stone J
229. Swerdlow AJ
230. Tamimi RM
231. Tapper WJ
232. Taylor JA
233. Terry MB
234. Thöne K
235. Tollenaar R
236. Tomlinson I
237. Truong T
238. Tzardi M
239. Ulmer HU
240. Untch M
241. Vachon CM
242. van Veen EM
243. Vijai J
244. Weinberg CR
245. Wendt C
246. Whittemore AS
247. Wildiers H
248. Willett W
249. Winqvist R
250. Wolk A
251. Yang XR
252. Yannoukakos D
253. Zhang Y
254. Zheng W
255. Ziogas A
256. Dunning AM
257. Thompson DJ
258. Chenevix-Trench G
259. Chang-Claude J
260. Schmidt MK
261. Hall P
262. Milne RL
263. Pharoah PDP
264. Antoniou AC
265. Chatterjee N
266. Kraft P
267. García-Closas M
268. Simard J
269. Easton DF
270. ABCTB Investigators, kConFab/AOCS Investigators, NBCS Collaborators
(2019) Polygenic risk scores for prediction of breast Cancer and breast Cancer subtypes
The American Journal of Human Genetics 104:21–34.

https://doi.org/10.1016/j.ajhg.2018.11.002
- PubMed
- Google Scholar
(2001)
Prediction of total genetic value using Genome-Wide dense marker maps

Genetics 157:1819–1829.
- PubMed
- Google Scholar
1. Mills MC
2. Rahal C
(2019) A scientometric review of genome-wide association studies
Communications Biology 2:9.

https://doi.org/10.1038/s42003-018-0261-x
- PubMed
- Google Scholar
(2017) Identifying genetic variants that affect viability in large cohorts
PLOS Biology 15:e2002458.

https://doi.org/10.1371/journal.pbio.2002458
- PubMed
- Google Scholar
(2018) Relative roles of race versus socioeconomic position in studies of health inequalities: a matter of interpretation
Annual Review of Public Health 39:169–188.

https://doi.org/10.1146/annurev-publhealth-040617-014230
- PubMed
- Google Scholar
1. Okbay A
2. Beauchamp JP
3. Fontana MA
4. Lee JJ
5. Pers TH
6. Rietveld CA
7. Turley P
8. Chen GB
9. Emilsson V
10. Meddens SF
11. Oskarsson S
12. Pickrell JK
13. Thom K
14. Timshel P
15. de Vlaming R
16. Abdellaoui A
17. Ahluwalia TS
18. Bacelis J
19. Baumbach C
20. Bjornsdottir G
21. Brandsma JH
22. Pina Concas M
23. Derringer J
24. Furlotte NA
25. Galesloot TE
26. Girotto G
27. Gupta R
28. Hall LM
29. Harris SE
30. Hofer E
31. Horikoshi M
32. Huffman JE
33. Kaasik K
34. Kalafati IP
35. Karlsson R
36. Kong A
37. Lahti J
38. van der Lee SJ
39. deLeeuw C
40. Lind PA
41. Lindgren KO
42. Liu T
43. Mangino M
44. Marten J
45. Mihailov E
46. Miller MB
47. van der Most PJ
48. Oldmeadow C
49. Payton A
50. Pervjakova N
51. Peyrot WJ
52. Qian Y
53. Raitakari O
54. Rueedi R
55. Salvi E
56. Schmidt B
57. Schraut KE
58. Shi J
59. Smith AV
60. Poot RA
61. St Pourcain B
62. Teumer A
63. Thorleifsson G
64. Verweij N
65. Vuckovic D
66. Wellmann J
67. Westra HJ
68. Yang J
69. Zhao W
70. Zhu Z
71. Alizadeh BZ
72. Amin N
73. Bakshi A
74. Baumeister SE
75. Biino G
76. Bønnelykke K
77. Boyle PA
78. Campbell H
79. Cappuccio FP
80. Davies G
81. De Neve JE
82. Deloukas P
83. Demuth I
84. Ding J
85. Eibich P
86. Eisele L
87. Eklund N
88. Evans DM
89. Faul JD
90. Feitosa MF
91. Forstner AJ
92. Gandin I
93. Gunnarsson B
94. Halldórsson BV
95. Harris TB
96. Heath AC
97. Hocking LJ
98. Holliday EG
99. Homuth G
100. Horan MA
101. Hottenga JJ
102. de Jager PL
103. Joshi PK
104. Jugessur A
105. Kaakinen MA
106. Kähönen M
107. Kanoni S
108. Keltigangas-Järvinen L
109. Kiemeney LA
110. Kolcic I
111. Koskinen S
112. Kraja AT
113. Kroh M
114. Kutalik Z
115. Latvala A
116. Launer LJ
117. Lebreton MP
118. Levinson DF
119. Lichtenstein P
120. Lichtner P
121. Liewald DC
122. Loukola A
123. Madden PA
124. Mägi R
125. Mäki-Opas T
126. Marioni RE
127. Marques-Vidal P
128. Meddens GA
129. McMahon G
130. Meisinger C
131. Meitinger T
132. Milaneschi Y
133. Milani L
134. Montgomery GW
135. Myhre R
136. Nelson CP
137. Nyholt DR
138. Ollier WE
139. Palotie A
140. Paternoster L
141. Pedersen NL
142. Petrovic KE
143. Porteous DJ
144. Räikkönen K
145. Ring SM
146. Robino A
147. Rostapshova O
148. Rudan I
149. Rustichini A
150. Salomaa V
151. Sanders AR
152. Sarin AP
153. Schmidt H
154. Scott RJ
155. Smith BH
156. Smith JA
157. Staessen JA
158. Steinhagen-Thiessen E
159. Strauch K
160. Terracciano A
161. Tobin MD
162. Ulivi S
163. Vaccargiu S
164. Quaye L
165. van Rooij FJ
166. Venturini C
167. Vinkhuyzen AA
168. Völker U
169. Völzke H
170. Vonk JM
171. Vozzi D
172. Waage J
173. Ware EB
174. Willemsen G
175. Attia JR
176. Bennett DA
177. Berger K
178. Bertram L
179. Bisgaard H
180. Boomsma DI
181. Borecki IB
182. Bültmann U
183. Chabris CF
184. Cucca F
185. Cusi D
186. Deary IJ
187. Dedoussis GV
188. van Duijn CM
189. Eriksson JG
190. Franke B
191. Franke L
192. Gasparini P
193. Gejman PV
194. Gieger C
195. Grabe HJ
196. Gratten J
197. Groenen PJ
198. Gudnason V
199. van der Harst P
200. Hayward C
201. Hinds DA
202. Hoffmann W
203. Hyppönen E
204. Iacono WG
205. Jacobsson B
206. Järvelin MR
207. Jöckel KH
208. Kaprio J
209. Kardia SL
210. Lehtimäki T
211. Lehrer SF
212. Magnusson PK
213. Martin NG
214. McGue M
215. Metspalu A
216. Pendleton N
217. Penninx BW
218. Perola M
219. Pirastu N
220. Pirastu M
221. Polasek O
222. Posthuma D
223. Power C
224. Province MA
225. Samani NJ
226. Schlessinger D
227. Schmidt R
228. Sørensen TI
229. Spector TD
230. Stefansson K
231. Thorsteinsdottir U
232. Thurik AR
233. Timpson NJ
234. Tiemeier H
235. Tung JY
236. Uitterlinden AG
237. Vitart V
238. Vollenweider P
239. Weir DR
240. Wilson JF
241. Wright AF
242. Conley DC
243. Krueger RF
244. Davey Smith G
245. Hofman A
246. Laibson DI
247. Medland SE
248. Meyer MN
249. Yang J
250. Johannesson M
251. Visscher PM
252. Esko T
253. Koellinger PD
254. Cesarini D
255. Benjamin DJ
256. LifeLines Cohort Study
(2016) Genome-wide association study identifies 74 loci associated with educational attainment
Nature 533:539–542.

https://doi.org/10.1038/nature17671
- PubMed
- Google Scholar
(2008) Polygenes, risk prediction, and targeted prevention of breast Cancer
New England Journal of Medicine 358:2796–2803.

https://doi.org/10.1056/NEJMsa0708739
- PubMed
- Google Scholar
(2015) Meta-analysis of the heritability of human traits based on fifty years of twin studies
Nature Genetics 47:702–709.

https://doi.org/10.1038/ng.3285
- PubMed
- Google Scholar
1. Popejoy AB
2. Fullerton SM
(2016) Genomics is failing on diversity
Nature 538:161–164.

https://doi.org/10.1038/538161a
- PubMed
- Google Scholar
1. Pritchard JK
2. Di Rienzo A
(2010) Adaptation - not by sweeps alone
Nature Reviews Genetics 11:665–667.

https://doi.org/10.1038/nrg2880
- PubMed
- Google Scholar
1. Pritchard JK
2. Przeworski M
(2001) Linkage disequilibrium in humans: models and data
The American Journal of Human Genetics 69:1–14.

https://doi.org/10.1086/321275
- PubMed
- Google Scholar
(2018) Detecting polygenic adaptation in admixture graphs
Genetics 208:1565–1584.

https://doi.org/10.1534/genetics.117.300489
- PubMed
- Google Scholar
1. Reckelhoff JF
(2001) Gender differences in the regulation of blood pressure
Hypertension 37:1199–1208.

https://doi.org/10.1161/01.HYP.37.5.1199
- PubMed
- Google Scholar
Book
1. Reich M
(2017)
Racial Inequality: A Political-Economic Analysis

Princeton University Press.
- Google Scholar
1. Rimfeld K
2. Krapohl E
3. Trzaskowski M
4. Coleman JRI
5. Selzam S
6. Dale PS
7. Esko T
8. Metspalu A
9. Plomin R
(2018) Genetic influence on social outcomes during and after the soviet era in Estonia
Nature Human Behaviour 2:269–275.

https://doi.org/10.1038/s41562-018-0332-5
- PubMed
- Google Scholar
1. Robinson MR
2. Kleinman A
3. Graff M
4. Vinkhuyzen AAE
5. Couper D
6. Miller MB
7. Peyrot WJ
8. Abdellaoui A
9. Zietsch BP
10. Nolte IM
11. van Vliet-Ostaptchouk JV
12. Snieder H
13. Medland SE
14. Martin NG
15. Magnusson PKE
16. Iacono WG
17. McGue M
18. North KE
19. Yang J
20. Visscher PM
(2017) Genetic evidence of assortative mating in humans
Nature Human Behaviour 1:16.

https://doi.org/10.1038/s41562-016-0016
- Google Scholar
(2019) Interpreting polygenic scores, polygenic adaptation, and human phenotypic differences
Evolution, Medicine, and Public Health 2019:26–34.

https://doi.org/10.1093/emph/eoy036
- PubMed
- Google Scholar
1. Ruby JG
2. Wright KM
3. Rand KA
4. Kermany A
5. Noto K
6. Curtis D
7. Varner N
8. Garrigan D
9. Slinkov D
10. Dorfman I
11. Granka JM
12. Byrnes J
13. Myres N
14. Ball C
(2018) Estimates of the heritability of human longevity are substantially inflated due to assortative mating
Genetics 210:1109–1124.

https://doi.org/10.1534/genetics.118.301613
- PubMed
- Google Scholar
1. Sella G
2. Barton NH
(2019) Thinking about the evolution of complex traits in the era of Genome-Wide association studies
Annual Review of Genomics and Human Genetics 20:461–493.

https://doi.org/10.1146/annurev-genom-083115-022316
- PubMed
- Google Scholar
(2019) Comparing within- and Between-Family polygenic score prediction
The American Journal of Human Genetics 105:351–363.

https://doi.org/10.1016/j.ajhg.2019.06.006
- PubMed
- Google Scholar
(2019) The missing diversity in human genetic studies
Cell 177:26–31.

https://doi.org/10.1016/j.cell.2019.02.048
- PubMed
- Google Scholar
1. Sohail M
2. Maier RM
3. Ganna A
4. Bloemendal A
5. Martin AR
6. Turchin MC
7. Chiang CW
8. Hirschhorn J
9. Daly MJ
10. Patterson N
11. Neale B
12. Mathieson I
13. Reich D
14. Sunyaev SR
(2019) Polygenic adaptation on height is overestimated due to uncorrected stratification in genome-wide association studies
eLife 8:e39702.

https://doi.org/10.7554/eLife.39702
- PubMed
- Google Scholar
Preprint
1. Speidel L
2. Forest M
3. Shi S
4. Myers S
(2019) A method for Genome-Wide genealogy estimation for thousands of samples
bioRxiv.

https://doi.org/10.1101/550558
- Google Scholar
(2017) Assortative mating for human height: a meta-analysis
American Journal of Human Biology 29:e22917.

https://doi.org/10.1002/ajhb.22917
- Google Scholar
1. Taylor AE
2. Jones HJ
3. Sallis H
4. Euesden J
5. Stergiakouli E
6. Davies NM
7. Zammit S
8. Lawlor DA
9. Munafò MR
10. Davey Smith G
11. Tilling K
(2018) Exploring the association of genetic factors with participation in the avon longitudinal study of parents and children
International Journal of Epidemiology 47:1207–1216.

https://doi.org/10.1093/ije/dyy060
- PubMed
- Google Scholar
Preprint
1. Telkar N
2. Reiker T
3. Walters RG
4. Lin K
5. Eriksson A
6. Gurdasani D
7. Gilly A
(2019) The transferability of lipid loci across african, asian and european cohorts
bioRxiv.

https://doi.org/10.1101/525170
- Google Scholar
(2018) The personal and clinical utility of polygenic risk scores
Nature Reviews Genetics 19:581–590.

https://doi.org/10.1038/s41576-018-0018-x
- PubMed
- Google Scholar
Preprint
1. Trejo S
2. Benjamin WD
(2019) Genetic nature or genetic nurture? quantifying Bias in analyses using polygenic scores
bioRxiv.

https://doi.org/10.1101/524850
- Google Scholar
1. Tropf FC
2. Lee SH
3. Verweij RM
4. Stulp G
5. van der Most PJ
6. de Vlaming R
7. Bakshi A
8. Briley DA
9. Rahal C
10. Hellpap R
11. Iliadou AN
12. Esko T
13. Metspalu A
14. Medland SE
15. Martin NG
16. Barban N
17. Snieder H
18. Robinson MR
19. Mills MC
(2017) Hidden heritability due to heterogeneity across seven populations
Nature Human Behaviour 1:757–765.

https://doi.org/10.1038/s41562-017-0195-1
- PubMed
- Google Scholar
(2019) An evolutionary compass for detecting signals of polygenic selection and mutational Bias
Evolution Letters 3:69–79.

https://doi.org/10.1002/evl3.97
- PubMed
- Google Scholar
(2015) Modeling linkage disequilibrium increases accuracy of polygenic risk scores
The American Journal of Human Genetics 97:576–592.

https://doi.org/10.1016/j.ajhg.2015.09.001
- PubMed
- Google Scholar
1. Vilhjálmsson BJ
2. Nordborg M
(2013) The nature of confounding in genome-wide association studies
Nature Reviews Genetics 14:1–2.

https://doi.org/10.1038/nrg3382
- PubMed
- Google Scholar
Preprint
1. Ware EB
2. Schmitz LL
3. Faul JD
4. Gard A
5. Mitchell C
6. Smith JA
7. Zhao W
8. Weir D
9. Kardia SLR
(2017) Heterogeneity in Polygenic scores for common human traits
bioRxiv.

https://doi.org/10.1101/106062
- Google Scholar
(2008) Genome-wide association analysis identifies 20 loci that influence adult height
Nature Genetics 40:575–583.

https://doi.org/10.1038/ng.121
- PubMed
- Google Scholar
1. Wojcik GL
2. Graff M
3. Nishimura KK
4. Tao R
5. Haessler J
6. Gignoux CR
7. Highland HM
8. Patel YM
9. Sorokin EP
10. Avery CL
11. Belbin GM
12. Bien SA
13. Cheng I
14. Cullina S
15. Hodonsky CJ
16. Hu Y
17. Huckins LM
18. Jeff J
19. Justice AE
20. Kocarnik JM
21. Lim U
22. Lin BM
23. Lu Y
24. Nelson SC
25. Park SL
26. Poisner H
27. Preuss MH
28. Richard MA
29. Schurmann C
30. Setiawan VW
31. Sockell A
32. Vahi K
33. Verbanck M
34. Vishnu A
35. Walker RW
36. Young KL
37. Zubair N
38. Acuña-Alonso V
39. Ambite JL
40. Barnes KC
41. Boerwinkle E
42. Bottinger EP
43. Bustamante CD
44. Caberto C
45. Canizales-Quinteros S
46. Conomos MP
47. Deelman E
48. Do R
49. Doheny K
50. Fernández-Rhodes L
51. Fornage M
52. Hailu B
53. Heiss G
54. Henn BM
55. Hindorff LA
56. Jackson RD
57. Laurie CA
58. Laurie CC
59. Li Y
60. Lin DY
61. Moreno-Estrada A
62. Nadkarni G
63. Norman PJ
64. Pooler LC
65. Reiner AP
66. Romm J
67. Sabatti C
68. Sandoval K
69. Sheng X
70. Stahl EA
71. Stram DO
72. Thornton TA
73. Wassel CL
74. Wilkens LR
75. Winkler CA
76. Yoneyama S
77. Buyske S
78. Haiman CA
79. Kooperberg C
80. Le Marchand L
81. Loos RJF
82. Matise TC
83. North KE
84. Peters U
85. Kenny EE
86. Carlson CS
(2019) Genetic analyses of diverse populations improves discovery for complex traits
Nature 570:514–518.

https://doi.org/10.1038/s41586-019-1310-4
- PubMed
- Google Scholar
1. Wood AR
2. Esko T
3. Yang J
4. Vedantam S
5. Pers TH
6. Gustafsson S
7. Chu AY
8. Estrada K
9. Luan J
10. Kutalik Z
11. Amin N
12. Buchkovich ML
13. Croteau-Chonka DC
14. Day FR
15. Duan Y
16. Fall T
17. Fehrmann R
18. Ferreira T
19. Jackson AU
20. Karjalainen J
21. Lo KS
22. Locke AE
23. Mägi R
24. Mihailov E
25. Porcu E
26. Randall JC
27. Scherag A
28. Vinkhuyzen AA
29. Westra HJ
30. Winkler TW
31. Workalemahu T
32. Zhao JH
33. Absher D
34. Albrecht E
35. Anderson D
36. Baron J
37. Beekman M
38. Demirkan A
39. Ehret GB
40. Feenstra B
41. Feitosa MF
42. Fischer K
43. Fraser RM
44. Goel A
45. Gong J
46. Justice AE
47. Kanoni S
48. Kleber ME
49. Kristiansson K
50. Lim U
51. Lotay V
52. Lui JC
53. Mangino M
54. Mateo Leach I
55. Medina-Gomez C
56. Nalls MA
57. Nyholt DR
58. Palmer CD
59. Pasko D
60. Pechlivanis S
61. Prokopenko I
62. Ried JS
63. Ripke S
64. Shungin D
65. Stancáková A
66. Strawbridge RJ
67. Sung YJ
68. Tanaka T
69. Teumer A
70. Trompet S
71. van der Laan SW
72. van Setten J
73. Van Vliet-Ostaptchouk JV
74. Wang Z
75. Yengo L
76. Zhang W
77. Afzal U
78. Arnlöv J
79. Arscott GM
80. Bandinelli S
81. Barrett A
82. Bellis C
83. Bennett AJ
84. Berne C
85. Blüher M
86. Bolton JL
87. Böttcher Y
88. Boyd HA
89. Bruinenberg M
90. Buckley BM
91. Buyske S
92. Caspersen IH
93. Chines PS
94. Clarke R
95. Claudi-Boehm S
96. Cooper M
97. Daw EW
98. De Jong PA
99. Deelen J
100. Delgado G
101. Denny JC
102. Dhonukshe-Rutten R
103. Dimitriou M
104. Doney AS
105. Dörr M
106. Eklund N
107. Eury E
108. Folkersen L
109. Garcia ME
110. Geller F
111. Giedraitis V
112. Go AS
113. Grallert H
114. Grammer TB
115. Gräßler J
116. Grönberg H
117. de Groot LC
118. Groves CJ
119. Haessler J
120. Hall P
121. Haller T
122. Hallmans G
123. Hannemann A
124. Hartman CA
125. Hassinen M
126. Hayward C
127. Heard-Costa NL
128. Helmer Q
129. Hemani G
130. Henders AK
131. Hillege HL
132. Hlatky MA
133. Hoffmann W
134. Hoffmann P
135. Holmen O
136. Houwing-Duistermaat JJ
137. Illig T
138. Isaacs A
139. James AL
140. Jeff J
141. Johansen B
142. Johansson Å
143. Jolley J
144. Juliusdottir T
145. Junttila J
146. Kho AN
147. Kinnunen L
148. Klopp N
149. Kocher T
150. Kratzer W
151. Lichtner P
152. Lind L
153. Lindström J
154. Lobbens S
155. Lorentzon M
156. Lu Y
157. Lyssenko V
158. Magnusson PK
159. Mahajan A
160. Maillard M
161. McArdle WL
162. McKenzie CA
163. McLachlan S
164. McLaren PJ
165. Menni C
166. Merger S
167. Milani L
168. Moayyeri A
169. Monda KL
170. Morken MA
171. Müller G
172. Müller-Nurasyid M
173. Musk AW
174. Narisu N
175. Nauck M
176. Nolte IM
177. Nöthen MM
178. Oozageer L
179. Pilz S
180. Rayner NW
181. Renstrom F
182. Robertson NR
183. Rose LM
184. Roussel R
185. Sanna S
186. Scharnagl H
187. Scholtens S
188. Schumacher FR
189. Schunkert H
190. Scott RA
191. Sehmi J
192. Seufferlein T
193. Shi J
194. Silventoinen K
195. Smit JH
196. Smith AV
197. Smolonska J
198. Stanton AV
199. Stirrups K
200. Stott DJ
201. Stringham HM
202. Sundström J
203. Swertz MA
204. Syvänen AC
205. Tayo BO
206. Thorleifsson G
207. Tyrer JP
208. van Dijk S
209. van Schoor NM
210. van der Velde N
211. van Heemst D
212. van Oort FV
213. Vermeulen SH
214. Verweij N
215. Vonk JM
216. Waite LL
217. Waldenberger M
218. Wennauer R
219. Wilkens LR
220. Willenborg C
221. Wilsgaard T
222. Wojczynski MK
223. Wong A
224. Wright AF
225. Zhang Q
226. Arveiler D
227. Bakker SJ
228. Beilby J
229. Bergman RN
230. Bergmann S
231. Biffar R
232. Blangero J
233. Boomsma DI
234. Bornstein SR
235. Bovet P
236. Brambilla P
237. Brown MJ
238. Campbell H
239. Caulfield MJ
240. Chakravarti A
241. Collins R
242. Collins FS
243. Crawford DC
244. Cupples LA
245. Danesh J
246. de Faire U
247. den Ruijter HM
248. Erbel R
249. Erdmann J
250. Eriksson JG
251. Farrall M
252. Ferrannini E
253. Ferrières J
254. Ford I
255. Forouhi NG
256. Forrester T
257. Gansevoort RT
258. Gejman PV
259. Gieger C
260. Golay A
261. Gottesman O
262. Gudnason V
263. Gyllensten U
264. Haas DW
265. Hall AS
266. Harris TB
267. Hattersley AT
268. Heath AC
269. Hengstenberg C
270. Hicks AA
271. Hindorff LA
272. Hingorani AD
273. Hofman A
274. Hovingh GK
275. Humphries SE
276. Hunt SC
277. Hypponen E
278. Jacobs KB
279. Jarvelin MR
280. Jousilahti P
281. Jula AM
282. Kaprio J
283. Kastelein JJ
284. Kayser M
285. Kee F
286. Keinanen-Kiukaanniemi SM
287. Kiemeney LA
288. Kooner JS
289. Kooperberg C
290. Koskinen S
291. Kovacs P
292. Kraja AT
293. Kumari M
294. Kuusisto J
295. Lakka TA
296. Langenberg C
297. Le Marchand L
298. Lehtimäki T
299. Lupoli S
300. Madden PA
301. Männistö S
302. Manunta P
303. Marette A
304. Matise TC
305. McKnight B
306. Meitinger T
307. Moll FL
308. Montgomery GW
309. Morris AD
310. Morris AP
311. Murray JC
312. Nelis M
313. Ohlsson C
314. Oldehinkel AJ
315. Ong KK
316. Ouwehand WH
317. Pasterkamp G
318. Peters A
319. Pramstaller PP
320. Price JF
321. Qi L
322. Raitakari OT
323. Rankinen T
324. Rao DC
325. Rice TK
326. Ritchie M
327. Rudan I
328. Salomaa V
329. Samani NJ
330. Saramies J
331. Sarzynski MA
332. Schwarz PE
333. Sebert S
334. Sever P
335. Shuldiner AR
336. Sinisalo J
337. Steinthorsdottir V
338. Stolk RP
339. Tardif JC
340. Tönjes A
341. Tremblay A
342. Tremoli E
343. Virtamo J
344. Vohl MC
345. Amouyel P
346. Asselbergs FW
347. Assimes TL
348. Bochud M
349. Boehm BO
350. Boerwinkle E
351. Bottinger EP
352. Bouchard C
353. Cauchi S
354. Chambers JC
355. Chanock SJ
356. Cooper RS
357. de Bakker PI
358. Dedoussis G
359. Ferrucci L
360. Franks PW
361. Froguel P
362. Groop LC
363. Haiman CA
364. Hamsten A
365. Hayes MG
366. Hui J
367. Hunter DJ
368. Hveem K
369. Jukema JW
370. Kaplan RC
371. Kivimaki M
372. Kuh D
373. Laakso M
374. Liu Y
375. Martin NG
376. März W
377. Melbye M
378. Moebus S
379. Munroe PB
380. Njølstad I
381. Oostra BA
382. Palmer CN
383. Pedersen NL
384. Perola M
385. Pérusse L
386. Peters U
387. Powell JE
388. Power C
389. Quertermous T
390. Rauramaa R
391. Reinmaa E
392. Ridker PM
393. Rivadeneira F
394. Rotter JI
395. Saaristo TE
396. Saleheen D
397. Schlessinger D
398. Slagboom PE
399. Snieder H
400. Spector TD
401. Strauch K
402. Stumvoll M
403. Tuomilehto J
404. Uusitupa M
405. van der Harst P
406. Völzke H
407. Walker M
408. Wareham NJ
409. Watkins H
410. Wichmann HE
411. Wilson JF
412. Zanen P
413. Deloukas P
414. Heid IM
415. Lindgren CM
416. Mohlke KL
417. Speliotes EK
418. Thorsteinsdottir U
419. Barroso I
420. Fox CS
421. North KE
422. Strachan DP
423. Beckmann JS
424. Berndt SI
425. Boehnke M
426. Borecki IB
427. McCarthy MI
428. Metspalu A
429. Stefansson K
430. Uitterlinden AG
431. van Duijn CM
432. Franke L
433. Willer CJ
434. Price AL
435. Lettre G
436. Loos RJ
437. Weedon MN
438. Ingelsson E
439. O'Connell JR
440. Abecasis GR
441. Chasman DI
442. Goddard ME
443. Visscher PM
444. Hirschhorn JN
445. Frayling TM
446. Electronic Medical Records and Genomics (eMEMERGEGE) Consortium, MIGen Consortium, PAGEGE Consortium, LifeLines Cohort Study.
(2014) Defining the role of common variation in the genomic and biological architecture of adult human height
Nature Genetics 46:1173–1186.

https://doi.org/10.1038/ng.3097
- PubMed
- Google Scholar
1. Yengo L
2. Sidorenko J
3. Kemper KE
4. Zheng Z
5. Wood AR
6. Weedon MN
7. Frayling TM
8. Hirschhorn J
9. Yang J
10. Visscher PM
11. GIANT Consortium
(2018) Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of european ancestry
Human Molecular Genetics 27:3641–3649.

https://doi.org/10.1093/hmg/ddy271
- PubMed
- Google Scholar
(2018) Relatedness disequilibrium regression estimates heritability without environmental Bias
Nature Genetics 50:1304–1310.

https://doi.org/10.1038/s41588-018-0178-9
- PubMed
- Google Scholar
(2019) Deconstructing the sources of genotype-phenotype associations in humans
Science 365:1396–1400.

https://doi.org/10.1126/science.aax3710
- PubMed
- Google Scholar
1. Zhang G
2. Bacelis J
3. Lengyel C
4. Teramo K
5. Hallman M
6. Helgeland Ø
7. Johansson S
8. Myhre R
9. Sengpiel V
10. Njølstad PR
11. Jacobsson B
12. Muglia L
(2015) Assessing the causal relationship of maternal height on birth size and gestational age at birth: a mendelian randomization analysis
PLOS Medicine 12:e1001865.

https://doi.org/10.1371/journal.pmed.1001865
- PubMed
- Google Scholar
1. Zhou B
2. Bentham J
3. Di Cesare M
4. Bixby H
5. Danaei G
6. Cowan MJ
7. Paciorek CJ
8. Singh G
9. Hajifathalian K
10. Bennett JE
11. Taddei C
12. Bilano V
13. Carrillo-Larco RM
14. Djalalinia S
15. Khatibzadeh S
16. Lugero C
17. Peykari N
18. Zhang WZ
19. Lu Y
20. Stevens GA
21. Riley LM
22. Bovet P
23. Elliott P
24. Gu D
25. Ikeda N
26. Jackson RT
27. Joffres M
28. Kengne AP
29. Laatikainen T
30. Lam TH
31. Laxmaiah A
32. Liu J
33. Miranda JJ
34. Mondo CK
35. Neuhauser HK
36. Sundström J
37. Smeeth L
38. Soric M
39. Woodward M
40. Ezzati M
41. Abarca-Gómez L
42. Abdeen ZA
43. Rahim HA
44. Abu-Rmeileh NM
45. Acosta-Cazares B
46. Adams R
47. Aekplakorn W
48. Afsana K
49. Aguilar-Salinas CA
50. Agyemang C
51. Ahmadvand A
52. Ahrens W
53. Al Raddadi R
54. Al Woyatan R
55. Ali MM
56. Alkerwi Ala'a
57. Aly E
58. Amouyel P
59. Amuzu A
60. Andersen LB
61. Anderssen SA
62. Ängquist L
63. Anjana RM
64. Ansong D
65. Aounallah-Skhiri H
66. Araújo J
67. Ariansen I
68. Aris T
69. Arlappa N
70. Aryal K
71. Arveiler D
72. Assah FK
73. Assunção MCF
74. Avdicová M
75. Azevedo A
76. Azizi F
77. Babu BV
78. Bahijri S
79. Balakrishna N
80. Bandosz P
81. Banegas JR
82. Barbagallo CM
83. Barceló A
84. Barkat A
85. Barros AJD
86. Barros MV
87. Bata I
88. Batieha AM
89. Baur LA
90. Beaglehole R
91. Romdhane HB
92. Benet M
93. Benson LS
94. Bernabe-Ortiz A
95. Bernotiene G
96. Bettiol H
97. Bhagyalaxmi A
98. Bharadwaj S
99. Bhargava SK
100. Bi Y
101. Bikbov M
102. Bjerregaard P
103. Bjertness E
104. Björkelund C
105. Blokstra A
106. Bo S
107. Bobak M
108. Boeing H
109. Boggia JG
110. Boissonnet CP
111. Bongard V
112. Braeckman L
113. Brajkovich I
114. Branca F
115. Breckenkamp J
116. Brenner H
117. Brewster LM
118. Bruno G
119. Bueno-de-Mesquita H B(as)
120. Bugge A
121. Burns C
122. Bursztyn M
123. de León AC
124. Cacciottolo J
125. Cameron C
126. Can G
127. Cândido APC
128. Capuano V
129. Cardoso VC
130. Carlsson AC
131. Carvalho MJ
132. Casanueva FF
133. Casas J-P
134. Caserta CA
135. Chamukuttan S
136. Chan AW
137. Chan Q
138. Chaturvedi HK
139. Chaturvedi N
140. Chen C-J
141. Chen F
142. Chen H
143. Chen S
144. Chen Z
145. Cheng C-Y
146. Dekkaki IC
147. Chetrit A
148. Chiolero A
149. Chiou S-T
150. Chirita-Emandi A
151. Cho B
152. Cho Y
153. Chudek J
154. Cifkova R
155. Claessens F
156. Clays E
157. Concin H
158. Cooper C
159. Cooper R
160. Coppinger TC
161. Costanzo S
162. Cottel D
163. Cowell C
164. Craig CL
165. Crujeiras AB
166. Cruz JJ
167. D'Arrigo G
168. d'Orsi E
169. Dallongeville J
170. Damasceno A
171. Dankner R
172. Dantoft TM
173. Dauchet L
174. De Backer G
175. De Bacquer D
176. de Gaetano G
177. De Henauw S
178. De Smedt D
179. Deepa M
180. Dehghan A
181. Delisle H
182. Deschamps V
183. Dhana K
184. Di Castelnuovo AF
185. Dias-da-Costa JS
186. Diaz A
187. Dickerson TT
188. Do HTP
189. Dobson AJ
190. Donfrancesco C
191. Donoso SP
192. Döring A
193. Doua K
194. Drygas W
195. Dulskiene V
196. Džakula A
197. Dzerve V
198. Dziankowska-Zaborszczyk E
199. Eggertsen R
200. Ekelund U
201. El Ati J
202. Ellert U
203. Elliott P
204. Elosua R
205. Erasmus RT
206. Erem C
207. Eriksen L
208. de la Peña JE
209. Evans A
210. Faeh D
211. Fall CH
212. Farzadfar F
213. Felix-Redondo FJ
214. Ferguson TS
215. Fernández-Bergés D
216. Ferrante D
217. Ferrari M
218. Ferreccio C
219. Ferrieres J
220. Finn JD
221. Fischer K
222. Föger B
223. Foo LH
224. Forslund A-S
225. Forsner M
226. Fortmann SP
227. Fouad HM
228. Francis DK
229. Franco MdoC
230. Franco OH
231. Frontera G
232. Fuchs FD
233. Fuchs SC
234. Fujita Y
235. Furusawa T
236. Gaciong Z
237. Gareta D
238. Garnett SP
239. Gaspoz J-M
240. Gasull M
241. Gates L
242. Gavrila D
243. Geleijnse JM
244. Ghasemian A
245. Ghimire A
246. Giampaoli S
247. Gianfagna F
248. Giovannelli J
249. Goldsmith RA
250. Gonçalves H
251. Gross MG
252. Rivas JPG
253. Gottrand F
254. Graff-Iversen S
255. Grafnetter D
256. Grajda A
257. Gregor RD
258. Grodzicki T
259. Grøntved A
260. Gruden G
261. Grujic V
262. Gu D
263. Guan OP
264. Gudnason V
265. Guerrero R
266. Guessous I
267. Guimaraes AL
268. Gulliford MC
269. Gunnlaugsdottir J
270. Gunter M
271. Gupta PC
272. Gureje O
273. Gurzkowska B
274. Gutierrez L
275. Gutzwiller F
276. Hadaegh F
277. Halkjær J
278. Hambleton IR
279. Hardy R
280. Harikumar R
281. Hata J
282. Hayes AJ
283. He J
284. Hendriks ME
285. Henriques A
286. Cadena LH
287. Herrala S
288. Heshmat R
289. Hihtaniemi IT
290. Ho SY
291. Ho SC
292. Hobbs M
293. Hofman A
294. Dinc GH
295. Hormiga CM
296. Horta BL
297. Houti L
298. Howitt C
299. Htay TT
300. Htet AS
301. Hu Y
302. Huerta JM
303. Husseini AS
304. Huybrechts I
305. Hwalla N
306. Iacoviello L
307. Iannone AG
308. Ibrahim MM
309. Ikram MA
310. Irazola VE
311. Islam M
312. Ivkovic V
313. Iwasaki M
314. Jackson RT
315. Jacobs JM
316. Jafar T
317. Jamrozik K
318. Janszky I
319. Jasienska G
320. Jelakovic B
321. Jiang CQ
322. Joffres M
323. Johansson M
324. Jonas JB
325. Jørgensen T
326. Joshi P
327. Juolevi A
328. Jurak G
329. Jureša V
330. Kaaks R
331. Kafatos A
332. Kalter-Leibovici O
333. Kamaruddin NA
334. Kasaeian A
335. Katz J
336. Kauhanen J
337. Kaur P
338. Kavousi M
339. Kazakbaeva G
340. Keil U
341. Boker LK
342. Keinänen-Kiukaanniemi S
343. Kelishadi R
344. Kemper HCG
345. Kengne AP
346. Kersting M
347. Key T
348. Khader YS
349. Khalili D
350. Khang Y-H
351. Khaw K-T
352. Kiechl S
353. Killewo J
354. Kim J
355. Klumbiene J
356. Kolle E
357. Kolsteren P
358. Korrovits P
359. Koskinen S
360. Kouda K
361. Koziel S
362. Kristensen PL
363. Krokstad S
364. Kromhout D
365. Kruger HS
366. Kubinova R
367. Kuciene R
368. Kuh D
369. Kujala UM
370. Kula K
371. Kulaga Z
372. Kumar RK
373. Kurjata P
374. Kusuma YS
375. Kuulasmaa K
376. Kyobutungi C
377. Laatikainen T
378. Lachat C
379. Lam TH
380. Landrove O
381. Lanska V
382. Lappas G
383. Larijani B
384. Laugsand LE
385. Laxmaiah A
386. Bao KLN
387. Le TD
388. Leclercq C
389. Lee J
390. Lee J
391. Lehtimäki T
392. Lekhraj R
393. León-Muñoz LM
394. Levitt NS
395. Li Y
396. Lilly CL
397. Lim W-Y
398. Lima-Costa MF
399. Lin H-H
400. Lin X
401. Linneberg A
402. Lissner L
403. Litwin M
404. Lorbeer R
405. Lotufo PA
406. Lozano JE
407. Luksiene D
408. Lundqvist A
409. Lunet N
410. Lytsy P
411. Ma G
412. Ma J
413. Machado-Coelho GLL
414. Machi S
415. Maggi S
416. Magliano DJ
417. Majer M
418. Makdisse M
419. Malekzadeh R
420. Malhotra R
421. Rao KM
422. Malyutina S
423. Manios Y
424. Mann JI
425. Manzato E
426. Margozzini P
427. Marques-Vidal P
428. Marrugat J
429. Martorell R
430. Mathiesen EB
431. Matijasevich A
432. Matsha TE
433. Mbanya JCN
434. Posso AJMD
435. McFarlane SR
436. McGarvey ST
437. McLachlan S
438. McLean RM
439. McNulty BA
440. Khir ASM
441. Mediene-Benchekor S
442. Medzioniene J
443. Meirhaeghe A
444. Meisinger C
445. Menezes AMB
446. Menon GR
447. Meshram II
448. Metspalu A
449. Mi J
450. Mikkel K
451. Miller JC
452. Miquel JF
453. Mišigoj-Durakovic M
454. Mohamed MK
455. Mohammad K
456. Mohammadifard N
457. Mohan V
458. Yusoff MFM
459. Møller NC
460. Molnár D
461. Momenan A
462. Mondo CK
463. Monyeki KDK
464. Moreira LB
465. Morejon A
466. Moreno LA
467. Morgan K
468. Moschonis G
469. Mossakowska M
470. Mostafa A
471. Mota J
472. Motlagh ME
473. Motta J
474. Muiesan ML
475. Müller-Nurasyid M
476. Murphy N
477. Mursu J
478. Musil V
479. Nagel G
480. Naidu BM
481. Nakamura H
482. Námešná J
483. Nang EEK
484. Nangia VB
485. Narake S
486. Navarrete-Muñoz EM
487. Ndiaye NC
488. Neal WA
489. Nenko I
490. Nervi F
491. Nguyen ND
492. Nguyen QN
493. Nieto-Martínez RE
494. Niiranen TJ
495. Ning G
496. Ninomiya T
497. Nishtar S
498. Noale M
499. Noboa OA
500. Noorbala AA
501. Noorbala T
502. Noto D
503. Al Nsour M
504. O'Reilly D
505. Oh K
506. Olinto MTA
507. Oliveira IO
508. Omar MA
509. Onat A
510. Ordunez P
511. Osmond C
512. Ostojic SM
513. Otero JA
514. Overvad K
515. Owusu-Dabo E
516. Paccaud FM
517. Padez C
518. Pahomova E
519. Pajak A
520. Palli D
521. Palmieri L
522. Panda-Jonas S
523. Panza F
524. Papandreou D
525. Parnell WR
526. Parsaeian M
527. Pecin I
528. Pednekar MS
529. Peer N
530. Peeters PH
531. Peixoto SV
532. Pelletier C
533. Peltonen M
534. Pereira AC
535. Pérez RM
536. Peters A
537. Petkeviciene J
538. Pham ST
539. Pigeot I
540. Pikhart H
541. Pilav A
542. Pilotto L
543. Pitakaka F
544. Plans-Rubió P
545. Polakowska M
546. Polašek O
547. Porta M
548. Portegies MLP
549. Pourshams A
550. Pradeepa R
551. Prashant M
552. Price JF
553. Puiu M
554. Punab M
555. Qasrawi RF
556. Qorbani M
557. Radic I
558. Radisauskas R
559. Rahman M
560. Raitakari O
561. Raj M
562. Rao SR
563. Ramachandran A
564. Ramos E
565. Rampal S
566. Reina DAR
567. Rasmussen F
568. Redon J
569. Reganit PFM
570. Ribeiro R
571. Riboli E
572. Rigo F
573. de Wit TFR
574. Ritti-Dias RM
575. Robinson SM
576. Robitaille C
577. Rodríguez-Artalejo F
578. Rodriguez-Perez del Cristo M
579. Rodríguez-Villamizar LA
580. Rojas-Martinez R
581. Rosengren A
582. Rubinstein A
583. Rui O
584. Ruiz-Betancourt BS
585. Horimoto ARVR
586. Rutkowski M
587. Sabanayagam C
588. Sachdev HS
589. Saidi O
590. Sakarya S
591. Salanave B
592. Salazar Martinez E
593. Salmerón D
594. Salomaa V
595. Salonen JT
596. Salvetti M
597. Sánchez-Abanto J
598. Sans S
599. Santos D
600. Santos IS
601. dos Santos RN
602. Santos R
603. Saramies JL
604. Sardinha LB
605. Margolis GS
606. Sarrafzadegan N
607. Saum K-U
608. Savva SC
609. Scazufca M
610. Schargrodsky H
611. Schneider IJ
612. Schultsz C
613. Schutte AE
614. Sen A
615. Senbanjo IO
616. Sepanlou SG
617. Sharma SK
618. Shaw JE
619. Shibuya K
620. Shin DW
621. Shin Y
622. Siantar R
623. Sibai AM
624. Silva DAS
625. Simon M
626. Simons J
627. Simons LA
628. Sjöström M
629. Skovbjerg S
630. Slowikowska-Hilczer J
631. Slusarczyk P
632. Smeeth L
633. Smith MC
634. Snijder MB
635. So H-K
636. Sobngwi E
637. Söderberg S
638. Solfrizzi V
639. Sonestedt E
640. Song Y
641. Sørensen TIA
642. Jérome CS
643. Soumare A
644. Staessen JA
645. Starc G
646. Stathopoulou MG
647. Stavreski B
648. Steene-Johannessen J
649. Stehle P
650. Stein AD
651. Stergiou GS
652. Stessman J
653. Stieber J
654. Stöckl D
655. Stocks T
656. Stokwiszewski J
657. Stronks K
658. Strufaldi MW
659. Sun C-A
660. Sundström J
661. Sung Y-T
662. Suriyawongpaisal P
663. Sy RG
664. Tai ES
665. Tammesoo M-L
666. Tamosiunas A
667. Tang L
668. Tang X
669. Tanser F
670. Tao Y
671. Tarawneh MR
672. Tarqui-Mamani CB
673. Taylor A
674. Theobald H
675. Thijs L
676. Thuesen BH
677. Tjonneland A
678. Tolonen HK
679. Tolstrup JS
680. Topbas M
681. Topór-Madry R
682. Tormo MJ
683. Torrent M
684. Traissac P
685. Trichopoulos D
686. Trichopoulou A
687. Trinh OTH
688. Trivedi A
689. Tshepo L
690. Tulloch-Reid MK
691. Tuomainen T-P
692. Tuomilehto J
693. Turley ML
694. Tynelius P
695. Tzourio C
696. Ueda P
697. Ugel E
698. Ulmer H
699. Uusitalo HMT
700. Valdivia G
701. Valvi D
702. van der Schouw YT
703. Van Herck K
704. van Rossem L
705. van Valkengoed IGM
706. Vanderschueren D
707. Vanuzzo D
708. Vatten L
709. Vega T
710. Velasquez-Melendez G
711. Veronesi G
712. Verschuren WMM
713. Verstraeten R
714. Victora CG
715. Viet L
716. Viikari-Juntura E
717. Vineis P
718. Vioque J
719. Virtanen JK
720. Visvikis-Siest S
721. Viswanathan B
722. Vollenweider P
723. Voutilainen S
724. Vrdoljak A
725. Vrijheid M
726. Wade AN
727. Wagner A
728. Walton J
729. Mohamud WNW
730. Wang M-D
731. Wang Q
732. Wang YX
733. Wannamethee SG
734. Wareham N
735. Wederkopp N
736. Weerasekera D
737. Whincup PH
738. Widhalm K
739. Widyahening IS
740. Wiecek A
741. Wijga AH
742. Wilks RJ
743. Willeit J
744. Willeit P
745. Williams EA
746. Wilsgaard T
747. Wojtyniak B
748. Wong TY
749. Wong-McClure RA
750. Woo J
751. Woodward M
752. Wu AG
753. Wu FC
754. Wu SL
755. Xu H
756. Yan W
757. Yang X
758. Ye X
759. Yiallouros PK
760. Yoshihara A
761. Younger-Coleman NO
762. Yusoff AF
763. Yusoff MFM
764. Zambon S
765. Zdrojewski T
766. Zeng Y
767. Zhao D
768. Zhao W
769. Zheng Y
770. Zhu D
771. Zimmermann E
772. Zuñiga Cisneros J
(2017) Worldwide trends in blood pressure from 1975 to 2015: a pooled analysis of 1479 population-based measurement studies with 19·1 million participants
The Lancet 389:37–55.

https://doi.org/10.1016/S0140-6736(16)31919-5
- Google Scholar

Article and author information

Author details

Hakhamanesh Mostafavi

Department of Biological Sciences, Columbia University, New York, United States

Present address
Department of Genetics, Stanford University, Stanford, United States

Contribution
Conceptualization, Data curation, Software, Formal analysis, Investigation, Visualization, Methodology

Contributed equally with
Arbel Harpak

For correspondence
hsm2137@columbia.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-1060-2844
Arbel Harpak

Department of Biological Sciences, Columbia University, New York, United States

Contribution
Conceptualization, Data curation, Software, Formal analysis, Investigation, Visualization, Methodology

Contributed equally with
Hakhamanesh Mostafavi

For correspondence
ah3586@columbia.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-3655-748X
Ipsita Agarwal

Department of Biological Sciences, Columbia University, New York, United States

Contribution
Formal analysis, Investigation, Visualization

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-8537-0008
Dalton Conley
1. Department of Sociology, Princeton University, Princeton, United States
2. Office of Population Research, Princeton University, Princeton, United States
Contribution
Conceptualization

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-5174-7222
Jonathan K Pritchard
1. Department of Genetics, Stanford University, Stanford, United States
2. Department of Biology, Stanford University, Stanford, United States
3. Howard Hughes Medical Institute, Stanford University, Stanford, United States
Contribution
Conceptualization

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-8828-5236
Molly Przeworski
1. Department of Biological Sciences, Columbia University, New York, United States
2. Department of Systems Biology, Columbia University, New York, United States
Contribution
Conceptualization, Resources, Supervision, Methodology, Project administration

For correspondence
mp3284@columbia.edu

Competing interests
Reviewing editor, eLife

"This ORCID iD identifies the author of this article:" 0000-0002-5369-9009

Funding

National Institute of General Medical Sciences (GM121372)

Molly Przeworski

National Human Genome Research Institute (HG008140)

Jonathan K Pritchard

Robert Wood Johnson Foundation (84337817)

Dalton Conley

Simons Foundation (633313)

Arbel Harpak

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This study has been conducted using the UK Biobank resource under application Number 11138, as approved by Columbia University Institutional Review Board, protocol AAAS2914. We are grateful to Daniel Belsky, Jeremy Berg, Graham Coop, Peter Donnelly, Doc Edge, Iain Mathieson, Augustine Kong, Magnus Nordborg, Vincent Plagnol, Guy Sella, Alex Young and members of the Przeworski and Sella labs for valuable discussions and to Doc Edge, Guy Sella, Graham Coop and Magnus Nordborg for comments on a draft of the manuscript. This work was funded by NIH GM121372 to MP, NIH HG008140 to JKP, a Robert Wood Johnson Foundation Pioneer Award (grant number 84337817) to DC and a Junior Fellowship from the Simons Society of Fellows (number 633313) to AH.

Ethics

Human subjects: This study has been conducted using the UK Biobank resource under application Number 11138, as approved by Columbia University Institutional Review Board, protocol AAAS2914.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

11,748

views
1,648

downloads
295

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Article PDF

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Hakhamanesh Mostafavi
Arbel Harpak
Ipsita Agarwal
Dalton Conley
Jonathan K Pritchard
Molly Przeworski

(2020)

Variable prediction accuracy of polygenic scores within an ancestry group

eLife 9:e48376.

https://doi.org/10.7554/eLife.48376

Categories and tags

Research organism

Human

Share this article

Cite this article

Variable prediction accuracy of polygenic scores within an ancestry group.

Differences in environmental variance alone do not explain the variable prediction accuracy.

Comparison of prediction accuracy of standard and sib-based polygenic scores.

Variable prediction accuracy within an ancestry group.

Variable prediction accuracy (measured as R2) within an ancestry group.

Dependence on the polygenic score model.

Estimating mean effect size across strata.

Variable prediction accuracy within an ancestry also seen using a linear mixed model.

Comparison of siblings and unrelated individuals in the UK Biobank with respect to age, SES, and sex ratio.

Comparison of siblings and unrelated individuals in the UK Biobank with respect to population structure.

Comparison of prediction accuracies of polygenic scores based on standard and sib-GWAS for simulated traits.

Simulation results for polygenic scores based on standard GWAS and sib-GWAS in the presence of indirect effects.

Simulation results for polygenic scores based on standard GWAS and sib-GWAS in the presence of assortative mating.

Comparison of prediction accuracies of polygenic scores based on standard and sib-GWAS matched for sex ratio.

Prediction accuracy of polygenic scores based on sib-and standard GWAS, for a range of traits.

Prediction accuracy for years of schooling, for individuals with 0 or 1 full sibling.

Variable prediction accuracy for binary traits, when measured as incremental AUC.

Variable prediction accuracy for binary disease phenotypes, measured as incremental AUC, in men versus women.

Comparison of prediction accuracies of polygenic scores (measured as R2) based on standard and sib-GWAS.

UK Biobank phenotype data used in this study and their corresponding data fields.

Genetic correlations across samples that vary by a study characteristic.

Sample sizes used for siblings and unrelated sets.

Qualifications to years of schooling conversion table.

Author details

Hakhamanesh Mostafavi

Present address

Contribution

Contributed equally with

For correspondence

Competing interests

Arbel Harpak

Contribution

Contributed equally with

For correspondence

Competing interests

Ipsita Agarwal

Contribution

Competing interests

Dalton Conley

Contribution

Competing interests

Jonathan K Pritchard

Contribution

Competing interests

Molly Przeworski

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Research organism

Further reading

Variable prediction accuracy (measured as $R^{2}$ ) within an ancestry group.

Comparison of prediction accuracies of polygenic scores (measured as $R^{2}$ ) based on standard and sib-GWAS.