Research Article

The nature of nurture: Effects of parental genotypes

See allHide authors and affiliations

Science  26 Jan 2018:
Vol. 359, Issue 6374, pp. 424-428
DOI: 10.1126/science.aan6877

Genetic variants provide a nurturing environment

Genetic variants in parents may affect the fitness of their offspring, even if the child does not carry the allele. This indirect effect is referred to as “genetic nurture.” Kong et al. used data from genome-wide association studies of educational attainment to construct polygenic scores for parents that only considered the nontransmitted alleles (see the Perspective by Koellinger and Harden). The findings suggest that genetic nurture is ultimately due to genetic variation in the population and is mediated by the environment that parents create for their children.

Science, this issue p. 424; see also p. 386


Sequence variants in the parental genomes that are not transmitted to a child (the proband) are often ignored in genetic studies. Here we show that nontransmitted alleles can affect a child through their impacts on the parents and other relatives, a phenomenon we call “genetic nurture.” Using results from a meta-analysis of educational attainment, we find that the polygenic score computed for the nontransmitted alleles of 21,637 probands with at least one parent genotyped has an estimated effect on the educational attainment of the proband that is 29.9% (P = 1.6 × 10−14) of that of the transmitted polygenic score. Genetic nurturing effects of this polygenic score extend to other traits. Paternal and maternal polygenic scores have similar effects on educational attainment, but mothers contribute more than fathers to nutrition- and heath-related traits.

How the human genome (nature) and the environment (nurture) work together to shape members of our species is a fundamental question, and any insights into this topic would be an important milestone. One challenge encountered by those who aspire to shed light on this matter is the lack of independence between the genome and the environment; thus, models that fail to account for this limitation are incomplete. Here we demonstrate how the genomes of close relatives—parents and siblings—can affect the proband through their contributions to the environment.

In animal studies, it is well established that alleles in a parent that are not transmitted to the offspring can nonetheless influence the offspring’s phenotypes (1, 2). Most examples involve effects manifested at the fetal stage, at which only the nontransmitted maternal alleles are relevant. In humans, the nontransmitted maternal alleles have been used to examine the potential causal relationships between the state of the mother during pregnancy and the outcomes of the child (3, 4). Here, for humans, we consider an alternative causal path where both paternal and maternal nontransmitted alleles can have effects that are mostly manifested after birth. A sequence variant that affects the phenotype of an individual is also likely to affect the parent from whom it was inherited (Fig. 1A). For some phenotypes, the state of a parent can influence the state of its child. This gives rise to a situation in which a child’s phenotype is influenced not only by the transmitted paternal and maternal alleles (TP and TM) (Fig. 1A) but also by the alleles that were not transmitted (NTP and NTM). A good example is educational attainment (EA) (5, 6): The EA of parents provides an environmental effect for children, but one that has a genetic component (7, 8). We call this phenomenon “genetic nurture.” The transmitted and nontransmitted alleles (Fig. 1A) both exert effects on the parents, and thus both induce genetic nurturing effects. The effect of the transmitted allele includes both its direct effect on the proband and its effect manifested through nurturing from blood relatives. Because the amount of trait variance explained is proportional to the square of effect size, genetic nurture could have a larger impact on variance explained through the transmitted alleles (by magnifying the direct effect) than the nontransmitted alleles. However, data on the nontransmitted alleles are needed to separate the genetic nurturing effects from the direct effects of the transmitted alleles. Specifically, Embedded Image (transmitted) and Embedded Image (nontransmitted) denote the respective estimated effects of the alleles when the paternal and maternal alleles are grouped together. Denoting the direct effect as δ, we propose to estimate it by Embedded Image. By calculating the difference, genetic nurturing effects and other potential confounding effects induced by population structure and assortative mating (9, 10) (see below) are cancelled out. Even though the implementations are different, this approach is related to the transmission-disequilibrium test (TDT) (11, 12), as both use nontransmitted alleles as controls (13). However, the potential effects of the nontransmitted alleles are ignored in the TDT. Mathematically, genetic nurture is a form of associative (or indirect) genetic effect, as defined by the animal-breeding literature (2). Genetic nurture is not limited to effects manifested through the phenotypes of the parents, as additional contributions (albeit probably substantially smaller ones) may go through grandparents and great-grandparents, for example (Fig. 1B). This study takes advantage of our human data to empirically examine the magnitudes of such effects for traits such as EA.

Fig. 1 Direct genetic effect and genetic nurturing effect.

(A) Alleles at an autosomal site carried by a parents-offspring trio are labeled with respect to the offspring (proband). TP and TM denote, respectively, the alleles transmitted from the father and the mother to the proband, and NTP and NTM denote the paternal and maternal alleles that are not transmitted. The transmitted alleles can influence the phenotype of the offspring, XO, through a direct path. The alleles of the parents, both transmitted and nontransmitted, can influence the parents’ phenotypes, YP and YM, and through them may have a nurturing effect on XO. This pathway combines a genetic effect (TP, NTP, TM, and NTM) on YP and YM with a nurturing effect (YP and YM) on XO. Note that although XO is often an individual trait of interest, Y would include a much broader set of phenotypes and is not completely known. (B) Red diamonds denote phenotypes of relatives; the blue diamond denotes the phenotype of the proband. Using the maternally transmitted allele as an example (denoted by T), we highlight that, in addition to the parents, the genetic nurturing effect can be manifested through the phenotypes of older ancestors and nonancestors such as siblings.

Estimating direct effects

To maximize the power to detect the effects of the nontransmitted alleles, we used 618,762 single-nucleotide polymorphisms (SNPs) spanning the genome to construct polygenic scores (14). The per-locus allele-specific weightings for the polygenic scores were derived from applying LDpred (15) to the results of a large genome-wide association study (8) (GWAS) of EA measured in years of education, with Icelandic data removed (13). The first analysis focused on 21,637 Icelandic probands, born between 1940 and 1983 (9139 males, 12,498 females), with EA data and at least one parent genotyped (Table 1). Because we could establish the parent of origin of the transmitted alleles (16), the nontransmitted allele from a genotyped parent was easily determined. polyTP and polyTM represent the polygenic scores computed from the transmitted paternal and maternal alleles, respectively, and polyNTP and polyNTM denote the corresponding polygenic scores for the nontransmitted alleles. To maximize power, we start by providing the results for polyT = polyTP + polyTM and polyNT = polyNTP + polyNTM. Here, polyTP and polyTM are scaled so that polyT has a mean of 0 and a variance of 1, and the trait EA is standardized to have a variance of 1. polyNTP and polyNTM were similarly calculated, and a 0 was imputed when the parent was not genotyped (13). Associations between EA and the polygenic scores computed from a joint analysis of polyT and polyNT that adjusts for sex, year of birth (yob) up to the cubic term, interactions between sex and yob, and 100 principal components (PCs) (13) are presented in Table 1. The estimated effect of polyT, Embedded Image, is 0.223 and significant [P = 1.6 × 10−174, calculated with genomic control adjustment (13, 17)]. Because both polyT and EA are standardized, the estimated fraction of the trait variance explained by polyT is Embedded Image (R2 in Table 1). However, the estimated effect of polyNT, Embedded Image, is also significant (P = 1.6 × 10−14). Thus, the estimated direct effect of polyT, Embedded Image, explains only Embedded Image of the trait variance, approximately one-half of R2. Noting that Embedded Image, the value of Embedded Image is determined (Table 1). In addition to the polygenic scores, individual results for 120 genome-wide significant SNPs (P < 5 × 10−8) in the Iceland-excluded meta-analysis are provided (table S1). Fifteen of the 120 SNPs (12.5%) have a one-tailed P value that is <0.05 for the nontransmitted alleles, which is more than that expected from noise [P = 1.5 × 10−3 (13)]. The average estimated effect of the nontransmitted alleles is 34.2% of that of the transmitted alleles. These results are consistent with previous observations that within-family EA effects calculated for dizygotic twins tended to be smaller than the standard GWAS effect estimates (7, 8).

Table 1 Decomposition of the observed effect of the polygenic score into direct, genetic nurturing, and confounding effects.

Traits: educational attainment (EA), age at first child (AGFC), high-density lipoprotein level (HDL), body mass index (BMI), fasting glucose level (FG), height (HT), cigarettes per day for smokers (CPD), and composite health trait (HLTH). Traits are standardized to have a variance of 1. N: number of probands with at least one parent genotyped; NNTP: number with father genotyped; NNTM: number with mother genotyped. Embedded Image and Embedded Image: estimated effects of the polygenic scores computed for the transmitted and nontransmitted alleles, respectively, when they are analyzed jointly. Embedded Image: estimated direct effect of the polygenic score. R2: estimated variance accounted for by the transmitted polygenic score, which captures both the direct effect and the genetic nurturing effect. Embedded Image: estimated variance accounted for by the direct effect alone. These fractions of variance explained are for trait values adjusted for sex, yob (year of birth), and PCs. Corresponding values for unadjusted trait values would be somewhat smaller (13). Embedded Image, Embedded Image, and Embedded Image: estimates, respectively, of the assortative mating–induced confounding effect for the direct effect component, the genetic nurturing effect, and the confounding effect of the genetic nurturing component.

View this table:

Assortative mating and estimating the genetic nurturing effect

We designate η to denote the magnitude of the genetic nurturing effect. Even though our analyses have adjustment for 100 PCs, which should have eliminated much of the population stratification–induced confounding, Embedded Image can still capture effects other than η. When there is assortative mating with respect to the genetic component underlying EA (10), a subtle confounding effect may result. Figure 2 illustrates a simple scenario in which the phenotype is assumed to be influenced by two loci: A and B. If there is assortative mating in the parents’ generation, it would lead to correlation of alleles between partners; for instance, the A alleles of the father (A1 and A2 in Fig. 2) will be correlated with the B alleles of the mother (B3 and B4) and vice versa. Consequently, the paternally transmitted A allele AP will be positively correlated with the maternally transmitted B allele BM, and AM will be correlated with BP. This correlation between alleles inherited from different parents is referred to as trans correlation, whereas the correlation between alleles inherited from the same parent (e.g., AP and BP) is referred to as cis correlation. This assortative mating–induced correlation differs from correlation between markers that are close physically, that is, within the same linkage-disequilibrium block. The latter correlation is mainly driven by the cis component, whereas the assortative mating–induced correlation could be dominated by the trans component. If trait association is calculated for locus A individually, the observed effect will capture both the effect of locus A and part of the effect of locus B. We let Embedded Image denote this added confounding effect. Similarly, assortative mating would also lead the A alleles to capture some of the nurturing effect of locus B, an effect denoted by Embedded Image. Under our model assumptions (13) Embedded ImageThe factor of 2 arises because the nontransmitted alleles have the same nurturing effects as the transmitted ones, and thus the transmitted and nontransmitted A alleles are capturing, through correlation, the nurturing effects of both the transmitted and nontransmitted B alleles. Additionally, we have the decompositionsEmbedded ImageandEmbedded Imagewhere E[ ] denotes expectation. Because both the transmitted and nontransmitted A alleles capture the confounding effects, Embedded Image remains an appropriate estimate of the direct effect δ. Locus A and locus B in Fig. 2 can be generalized to represent two nonoverlapping sets of loci. For our study, we think of locus A as the EA polygenic score, whereas locus B represents the genetic component of EA that is statistically orthogonal to locus A (under a scenario of no assortative mating). The mathematical relationships highlighted above continue to hold for the polygenic scores, either exactly or approximately. By using a method for estimating heritability that also incorporates data on the nontransmitted alleles (18), we estimate the full genetic component of EA to have a direct effect that explains 17.0% of the variance of EA. In other words, polyT is estimated to be 2.45/17.0 = 14.4% of the full genetic component, whereas the remaining 85.6% corresponds to the B components. From this estimate, we extrapolate the correlations observed between the paternal polygenic scores (polyTP and polyNTP) and the maternal polygenic scores (polyTM and polyNTM) to estimate the correlations between them and the unobserved B components (13). From the latter, Embedded Image and Embedded Image are estimated as 0.065 and 0.130, respectively. For this calculation, we avoided making the assumption that assortative mating between parents was manifested only through correlation of their EAs, which would have led to lower estimates for the Embedded Image values (13). From these estimates and the above equations, Embedded Image, Embedded Image, and Embedded Image were computed and presented in Table 1 as fractions of Embedded Image. For EA, Embedded Image accounts for ~75% of the value of Embedded Image and Embedded Image is 31.9% of Embedded Image. Finally, we note that assortative mating occurring before the parents’ generation could lead to additional confounding. However, this effect appears to be negligible in our study, as after adjustment for 100 PCs, the within-parent correlation of the transmitted and nontransmitted polygenic scores is actually negative (but P > 0.05) (13).

Fig. 2 Correlation and confounding induced by assortative mating.

An example of two loci, A and B, contributing to the phenotype. Through assortative mating, alleles in the father become correlated with alleles in the mother. Consequently, the transmitted paternal alleles (AP and BP) become correlated with the maternally transmitted alleles (AM and BM). This correlation between alleles with different parental origins is referred to as trans correlation, whereas correlation between alleles with the same parental origins (e.g., AP and BP) is referred to as cis correlation. When AP/AM and BM/BP are correlated, association analysis between the phenotype and A alone will also capture part of the effect of B.

Direct and nurturing effects on other traits

The EA polygenic score is associated with other quantitative traits in our database. Among them, those with the strongest statistical significance (Table 1) are age at first child (AGFC) (19), high-density lipoprotein level (HDL) (20), body mass index (BMI) (21), fasting glucose level (FG) (22), height (HT) (23), and cigarettes smoked per day by smokers (CPD) (24). The effects of the transmitted and nontransmitted EA polygenic scores on these phenotypes were estimated as before for the EA phenotype (Table 1). Although the fraction of variance explained by polyT (R2) is smaller than that for EA, the effect of polyNT is statistically significant. Moreover, except for BMI, the ratio Embedded Image is higher for these traits than for EA and exceeds 1 for HT.

Parent of origin

Table 2 provides the estimated effects of polyTP, polyTM, polyNTP, and polyNTM separately (13). For EA, Embedded Image, the estimated effect of polyNTP, is significant (P = 5.2 × 10−7), and its value is nearly identical to that of Embedded Image (the higher P value for polyNTP is due to fewer fathers genotyped than mothers). This indicates that the effect observed for polyNT is not driven by epigenetic effects such as imprinting or genetic interactions between fetus and mother in the womb and does capture a genetic nurturing effect [also see tables S2 and S3, which have results for polygenic scores calculated without SNPs in imprinted regions (25)]. However, even with both parents contributing to genetic nurture, the magnitude of the effect can differ between fathers and mothers. We designate Embedded Image and Embedded Image to denote the paternal and maternal genetic nurturing effects, respectively. Because the transmitted alleles also contribute to the nurturing effect, we use a weighted average of Embedded Image and Embedded Image, with weights proportional to the inverse square of the standard error (13), to estimate Embedded Image (Table 2). Combining this estimate with Embedded Image from Table 1, considered as an estimate of a weighted average of Embedded Image and Embedded Image with weights proportional to the numbers of fathers and mothers genotyped, we calculated individual estimates of Embedded Image and Embedded Image (13), denoted by Embedded Image and Embedded Image, and the ratio Embedded Image (Table 2). For EA, Embedded Imageis estimated to be 0.011, but it is not significantly different from zero (P = 0.31)—that is, the ratio Embedded Image is not significantly different from 1. For all of the other six traits, Embedded Image but was significant only for HT (Embedded Image, P = 1.1 × 10−2). HDL and FG have P values between 0.05 and 0.10. To increase power, for individuals for whom we had data for one or more of the five health- and nutrition-related traits (HDL, BMI, FG, HT, and CPD), a composite health trait (HLTH) was constructed by taking the sum of the standardized values of the available traits (positive signs for HDL and HT; negative signs for BMI, FG, and CPD) and dividing it by the square root of the number of trait values summed. It was then standardized to have a variance of 1. For HLTH, Embedded Image has a larger value than that for the individual health- and nutrition-related traits and is highly significant (P = 8.9 × 10−11) (Table 1). Both Embedded Image and Embedded Image are significant, but Embedded Image with a P value of 4.8 × 10−3 (Table 2). This supports the notion that mothers have a stronger nurturing effect than fathers on the health of the child.

Table 2 Parent-of-origin–specific effects of the polygenic scores.

Traits are as defined in Table 1. Embedded Image, Embedded Image, Embedded Image, and Embedded Image: estimates of the effect of the polygenic scores polyTP, polyTM, polyNTP, and polyNTM, respectively. Embedded Image and Embedded Image: estimates of the paternal and maternal genetic nurturing effects, respectively.

View this table:

Variance explained and effects of siblings

The existence of genetic nurture complicates the estimation and interpretation of heritability (18). For example, maternal effects have been shown to affect heritability estimates from animal-breeding data (26). Though distinct from the direct effect of inherited genetic variants, we demonstrate here how genetic nurture can be measured and taken into consideration. If there are two uncorrelated variants of the same frequency, one having a direct effect δ only and the other having a nurturing effect η only, then the variance explained is proportional to δ2 + η2. By comparison, if one variant has both effects, then the variance explained is proportional to (δ + η)2 = δ2 + 2δη + η2 (Fig. 3), with the extra 2δη term. Moreover, (δ + η)2 captures the effect of the transmitted allele(s) only; the phenotypic variance accounted for by the transmitted and nontransmitted alleles together is proportional to (δ + η)2 + η2 (Fig. 3). With EA, Embedded Image, Embedded Image, and Embedded Image. Assuming that the direct effect alone explains 17.0% of the variance, the variance explained by the transmitted alleles with the nurturing effects included increases to 17.0% × 1.74 = 29.6%. Additionally including the nontransmitted alleles would increase the variance explained to 17.0% × 1.84 = 31.3%. The genetic nurturing effect not only magnifies the variance explained, it also induces an even larger amplification of the phenotypic correlations of parents and offspring and of siblings (Fig. 3) (13). Also worth noting is that the 2δη term highlighted above does not exist for adopted children, as then both alleles of a parent would be nontransmitted.

Fig. 3 Variance explained and induced correlations between phenotypes of parent-offspring and siblings.

Results are displayed as a function of the ratio η/δ. The y axis is the relative amplification; that is, various measures relative to what can be accounted for by the direct effect alone, the latter proportional to δ2. The total variance explained by the transmitted alleles is proportional to (δ + η)2 [the plotted curve is hence (δ + η)2/ δ2], whereas the total variance explained by the transmitted alleles plus the nontransmitted alleles is proportional to (δ + η)2 + η2. Formulas for the induced parent-offspring and sibling correlations are derived (13). η, magnitude of genetic nurturing effect; δ, direct effect; T, transmitted; NT, nontransmitted.

Genetic nurture could go through a sibling (Fig. 1B) if, as proposed (27), the phenotypes of the proband are influenced by the phenotypes or behavior of a sibling. On the basis of the genealogy, for each EA proband who has at least one sibling, the sibling most likely to have the biggest effect on the proband was identified as follows. If the proband has older siblings, the older sibling with a yob closest to the proband was selected (monozygotic twins were excluded, but we count a dizygotic twin of the proband as an older sibling). If the proband is the eldest child, a younger sibling with the closest yob was chosen. There are 7798 probands whose chosen sibling is genotyped and whose parents are both genotyped. A polygenic score, denoted by polyTS, was computed using the alleles transmitted from the parents to the sibling. The EA of the proband was then regressed on polyT, polyNT, and polyTS jointly. The effect of polyTS is significant (P = 0.015) and is estimated to be 24.1% (95% confidence interval: 4.7 to 43.6%) of the direct effect. The uncertainty is large because polyTS is strongly correlated with polyT and polyNT. One compensation is that, having adjusted for both polyT and polyNT, the estimated effect of polyTS is free of confounding from assortative mating.

Heritability is defined as the fraction of phenotypic variance explained by direct effects alone. The presence of parental genetic nurture introduces bias to estimates of heritability from GREML (genomic relatedness–based restricted maximum likelihood)–type methods (28), such as those embodied in the software package GCTA (29), that use correlations due to transmitted alleles without distinction between direct genetic effects and genetic nurturing effects (18). By contrast, heritability estimates from comparing correlations between monozygotic versus dizygotic twins (30) are unaffected as the effects of parental genetic nurture are cancelled out. However, when genetic nurturing effects that go through the phenotypes of a sibling or twin are present, both twin-based heritability estimates (31) and estimates from GREML-type methods are affected.

The nature of genetic nurture and other polygenic scores

To further use the EA trait data, we performed analyses that treated the nontransmitted polygenic score of a genotyped parent as missing if the EA of that parent was unknown. For these data, (unadjusted) estimates of Embedded Image were calculated as before (table S4). Also given are estimates of Embedded Image adjusted for the EAs of the parents, obtained by adding the latter to the explanatory variables in the regressions. For EA, AGFC, HT, and HLTH, the adjusted estimate remains significant (P < 0.005), and the ratio of the adjusted versus unadjusted estimate is, respectively, 47.6, 63.0, 80.3, and 68.6%. This indicates that the EA of the parent is an important part of the parental phenotypes (Y in Fig. 1A) through which genetic nurture operates, but it is far from all of it. The EA polygenic score is likely associated with intelligence, conscientiousness, and future planning. Parents with a high score enhance the nurturing of their offspring through many behaviors, not exclusively through their own EA.

To contrast the results presented for the polygenic score constructed from a GWAS of EA (EA polygenic score), we examined polygenic scores constructed from GWASs of HT (32) (HT polygenic score) and BMI (33) (BMI polygenic score). (Results corresponding to Table 1 are in tables S5 and S6.) Noting that the HT and BMI polygenic scores are, respectively, positively (r = 0.087) and negatively correlated (r = −0.146) with the EA polygenic score, we computed HT and BMI polygenic scores adjusted for the EA polygenic score by regressing the former on the latter and calculating the residuals (tables S7 and S8). Whereas the unadjusted nontransmitted polygenic score has a few significant associations (tables S5 and S6), with adjustment (tables S7 and S8) the only significant effect of the nontransmitted polygenic score is between the HT trait and the nontransmitted HT polygenic score. Furthermore, most of this observed effect is estimated to be due to confounding from assortative mating.


Through the study of the nontransmitted alleles, we demonstrated that genetic nurturing effects exist and can have an impact on variance explained. These results also reveal that the observed effects from GWAS do not necessarily reflect direct effects alone. They can be amplified by genetic nurturing effects and, to a lesser extent, assortative mating–induced confounding. Owing to power considerations, we mostly studied variants as an aggregate. However, given the complexity of the EA trait (6) and our observed effects of the EA polygenic score on other traits, for individual variants, the ratio of the genetic nurturing effect versus the direct effect must have variations both between and within traits. Thus, we should aim to gather enough data to perform GWAS with the nontransmitted alleles. This would add insight into the pathway(s) through which the effect of an individual variant is manifested, as well as enable a better understanding of some pleiotropic effects (34).

Although genes have been shown to affect the environment (24, 35, 36), the contribution of a genetic effect manifested through nurturing has mostly been ignored in GWAS. Results here highlight the importance of family data.

Our focus has been on genetic nurture in one direction, but the effects are likely to be bidirectional. For a parent-offspring pair, the magnitude of the effect in the direction of parent to offspring is likely to dominate the effect in the opposite direction. However, with siblings and twins, the effects would be reciprocal.

Our analyses implicitly assume that direct genetic effects and genetic nurturing effects are additive, but interactive effects could certainly exist and further complicate the interpretation of observed effects. Moreover, alleles other than those in the parents can also have an effect; for example, the genetic makeup of the population of the probands could also be an important environmental contributor to their phenotypes.

Supplementary Materials

Materials and Methods

Tables S1 to S9

References (3743)

References and Notes

  1. Materials and methods are available as supplementary materials.
Acknowledgments: We thank A. Okbay for providing the EA meta-analysis results with Icelandic data removed; J. Hirshhorn for pointing out that the nontransmitted alleles could be capturing some confounding effects due to assortative mating; and the GIANT consortium, A. Wood, and A. Locke for assisting us to obtain meta-analysis results for HT and BMI with Icelandic data removed. B.V.H. is an associate professor in the School of Science and Engineering at Reykjavík University. A summary of the data used in this manuscript is in table S9. Icelandic law allows for unimpeded sharing of summary-level data. However, the law does not allow the sharing of individual-level data on genotypes and phenotypes outside of Iceland.
View Abstract

Navigate This Article