Abstract
Parents pass on both their genes and environment to offspring, prompting debate about the relative importance of nature versus nurture in the inheritance of complex traits. Advances in molecular genetics now make it possible to quantify an individual’s genetic predisposition to a trait via his or her ‘polygenic score’. However, part of the risk captured by an individual’s polygenic score may actually be attributed to the genotype of their parents. In the most well-studied example of this indirect ‘genetic nurture’ effect, about half the genetic contribution to educational attainment was found to be attributed to parental alleles, even if those alleles were not inherited by the child. Refractive errors, such as myopia, are a common cause of visual impairment and pose high economic and quality-of-life costs. Despite strong evidence that refractive errors are highly heritable, the extent to which genetic risk is conferred directly via transmitted risk alleles or indirectly via the environment that parents create for their children is entirely unknown. Here, an instrumental variable analysis in 1944 pairs of adult siblings from the United Kingdom was used to quantify the proportion of the genetic risk (‘single nucleotide polymorphism (SNP) heritability’) of refractive error contributed by genetic nurture. We found no evidence of a contribution from genetic nurture: non-within-family SNP-heritability estimate = 0.213 (95% confidence interval 0.134–0.310) and within-family SNP-heritability estimate = 0.250 (0.152–0.372). Our findings imply the genetic contribution to refractive error is principally an intrinsic effect from alleles transmitted from parents to offspring.
Similar content being viewed by others
Introduction
A mismatch between the focal and axial lengths of the eye causes blurry vision, known as myopia or hyperopia (Fig. 1). Approximately 30–50% of children in the USA and Europe and 50–90% of children in East and South-East Asia develop myopia by early adulthood [1, 2]. Refractive errors impose a heavy economic burden on society in terms of the need for sight tests and visual correction (glasses, contact lenses or refractive surgery). Myopia also imposes a healthcare and societal burden via an increased risk of visual impairment and blindness, through myopic macular degeneration, retinal detachment and a heightened risk of glaucoma and cataract [3].
Refractive error is a highly heritable trait [4,5,6,7], with most cases of myopia developing when genetically susceptible individuals are exposed to environmental risk factors [8]. To date, insufficient time outdoors and intensive education are regarded as major risk factors for myopia [1, 9,10,11,12]. The interplay between genetic and lifestyle risk factors is a major focus of current research efforts in refractive error [13,14,15].
Genetic nurture
The genetic risk transmitted from parents to offspring generally occurs via the transmission of parental alleles to children. However, alleles present in the parents can potentially influence the phenotype of a child through the child’s environment, independently of their transmission to the child (‘genetic nurture’) [16]. Accordingly, children may ‘inherit both phenotype-associated SNPs and phenotype-associated environments from parents’ [17]. Figure 2 illustrates how the family environment of a child (the proband) can be influenced by the proband’s own genotype and the genotypes of the two parents and sibling(s).
A genetic variant displaying a significant genotype-phenotype association provides evidence of a causal path affecting the phenotype [18]. However, in the context of genetic nurture, it is apparent that this causal path could act partly or fully via the environment of a proband (which is affected by the genetics of his/her relatives). The presence of genetic nurture may have profound implications. For instance, when aiming to translate discoveries from genome-wide association studies (GWASs) into treatments, it is critical to determine whether disease-associated SNPs act within a proband or through his/her environment. If the former is true, then an intervention (such as a drug) would need to be targeted at the proband, while if the latter is true, an intervention may need to be targeted at the proband’s environment.
To date, the role of genetic nurture has only been examined for a few traits. Most studies have focused on either educational attainment or birth weight [16, 19,20,21,22,23,24,25,26,27,28]. These studies have demonstrated that as much as half of the genetic contribution to educational attainment may arise via genetic nurture rather than via transmitted alleles [16, 20, 22,23,24,25,26]. This is suggestive, for example, that alleles carried by the parents and influencing their own educational attainment [29] may lead them to influence the development of educational attainment in children by placing more or less than average emphasis on their child’s academic environment, e.g. via encouragement to study or provision of resources, such as books and private tuition. In the case of birth weight, several genome-wide significantly associated variants have been shown to act via the foetal environment [30].
Polygenic scores and SNP-heritability
A polygenic score (PGS) quantifies a person’s genetic predisposition for a specific trait [31] and can be harnessed to provide an unbiased estimate of the SNP-heritability by implementing a ‘split-sample’ GWAS approach [32, 33]. This involves dividing the available GWAS sample at random into two equally sized groups of participants and running a GWAS in each subsample, thereby yielding two independent estimates of the SNP effect sizes (beta coefficients). An unbiased estimate of the SNP-heritability can then be obtained by using one of the resultant polygenic scores as an instrumental variable for the other [32, 33].
When a GWAS is carried out in a sample of unrelated individuals, the SNP beta coefficients will not only capture effects of transmitted alleles, but also any effects due to genetic nurture, population stratification and assortative mating. By contrast, an appropriately designed family-based GWAS will only capture the genetic effects of transmitted alleles. Here, we carried out analyses to assess the evidence for genetic nurture by comparing SNP-heritability estimates obtained in within-family vs. non-within-family analyses [21]. Given the evidence for a causal role of education in predisposing children to myopia [11, 12, 34], we hypothesized that genetic nurture would make a major contribution to the SNP-heritability of refractive error. We also performed the same analysis for a trait known to show genetic nurture, educational attainment.
Methods
Participants and phenotypes
UK Biobank is a prospective study of health and well-being of adults living in the United Kingdom. Ethical approval for the study was obtained from the NHS Research Ethics Committee (Reference: 11/NW/0382) and all participants provided signed, informed consent. Recruitment occurred during 2006–2010, when approximately 500,000 individuals aged 37–73 years attended the baseline assessment visit [35]. At this visit, participants completed a structured interview that included the question, “Which of the following qualifications do you have? (You can select more than one)”. A follow-up question was asked of those who did not report having a university or college degree, “At what age did you complete your continuous full-time education?”. Educational attainment (EduYears) was defined as age participants reported completing their full-time education (minus 5 years, which is the age the participants started formal schooling), except that those who reported leaving school before age 15 or after age 21 years-old were assigned a school leaving age of 15 or 21 (thus, an EduYears of 10 or 16 years), respectively [12]. Participants with a university or college degree were assigned a school leaving age of 21 (an EduYears of 16 years). An ophthalmic assessment was introduced late in the UK Biobank recruitment period, with only about a quarter of participants undergoing this assessment. The mean spherical equivalent (sphere plus 0.5 × cylinder) for each eye was calculated from the autorefraction readings. The average mean spherical equivalent (avMSE) of the two eyes was taken as the participant’s refractive error phenotype value [12].
Selection of siblings for genetic nurture analysis and participants for split-sample genome-wide association studies
An overview of the selection scheme is shown in Fig. 3. All analyses were restricted to participants of European ancestry. Siblings were identified as described by Bycroft et al. [36]. Two outlier pairs whose kinship vs. identity-by-state allele sharing pattern did not cluster with the other sibling pairs were excluded (Supplementary Fig. S1). For sibships with more than 2 individuals with information on the traits of interest, we randomly selected one pair. We refer to this sample of siblings as the ‘genetic nurture analysis sample’ (Fig. 3). For the split-sample GWASs for each trait-of-interest, we selected participants who were unrelated to any individual in the genetic nurture analysis sample and who were also unrelated to one another (Fig. 3), using the R package igraph (Bycroft et al. [36]). The pool of available participants for the split-sample GWAS analyses was divided at random into two equally sized groups. The sample size for the split-sample GWAS for avMSE was 2 × n = 43,168 and for EduYears was 2 × n = 147,813.
Genome-wide association (GWA) analyses and polygenic scores
A list of autosomal genetic variants in the HapMap3 reference panel was downloaded from the LDSC website [37] (https://alkesgroup.broadinstitute.org/LDSCORE/w_hm3.snplist.bz2). Using PLINK v2.00a2LM [38], imputed genotype data from UK Biobank release 3 were filtered to retain only variants present on the HapMap3 list that had a missing genotype rate <0.02, restricted to participants of European ancestry with a missing genotyping rate <0.05. The variants were sorted in descending order of minor allele frequency (MAF) and the top 999,998 variants were selected for inclusion in the genome-wide association studies (these retained variants all had MAF > 0.02). GWA analyses were performed separately for split sample 1 and split sample 2 for each trait of interest (n = 43,168 for avMSE and n = 147,813 for EduYears) using the --predBetasFile option in BOLT v2.3.5 [39]. Age, age-squared, sex, genotyping array and the first 10 genetic principal components (PCs) were included as covariates. The above GWA analyses were performed for an infinitesimal model. Polygenic scores for the genetic nurture samples were calculated using the PLINK v1.9 --score function, with the weights obtained from BOLT. In total, four polygenic scores were calculated: two polygenic scores for avMSE and two polygenic scores for EduYears. The incremental R2 was defined as the variance in the trait-of-interest explained by the polygenic score in a linear regression model, over and above the variance explained by a baseline linear regression model that included the predictors age, age-squared, sex, genotyping array and the first 10 genetic PCs.
Assessment of genetic nurture
Each participant’s genetic propensity for the trait of interest was quantified using a polygenic score (PGS), using Eq. (1).
where, βj is a vector of weights for a set of J genetic variants indexed by j = 1, 2, … J and k = 0, 1 or 2 is the number of effect alleles for variant j carried by the participant. The weights βj are obtained in a GWAS for the trait of interest, after accounting for linkage disequilibrium (LD). The aim of the current study was to examine the extent to which the genetic risk captured by a polygenic score for refractive error is mediated by genetic nurture effects. In this context, use of an imprecise polygenic score will tend to underestimate the contribution of the mediating pathway [31]. One method to avoid this source of bias is to estimate \(\widehat {\beta _j}\) using a split-sample GWAS [32, 33]. This yields two independent estimates of the genetic variant weights \(\widehat {\beta _j}\) and an unbiased estimate of PGS can be obtained by using one of the \(\widehat {\beta _j}\) estimates as an instrumental variable for the other [32, 33].
We fitted Eq. (2), following Brumpton et al. [40].
where, α0 is the mean phenotype level in the sample, yfi is the phenotype of sibling i = A or B from family f who has an instrumented polygenic score of \(\widehat {{{{{\rm{PGS}}}}}}_{(m = 2)fi}\), and Wf is the mean level of the phenotype in family f. Terms \({{{{{\rm{PGS}}}}}}_{\left( {m = 1} \right)fi}\)and PGS(m=2)fi are the polygenic scores derived from the first and second tranches of a split-sample GWA, in which the first polygenic score is used as an instrumental variable for the second. Age, age-squared, sex, genotyping array and the first 10 principal components (PCs) derived by Bycroft et al. [36] were included in Eq. (2). The parameter α1 in Eq. (2) is a within-family measure of the genetic contribution to the trait of interest, which we refer to as: \(\alpha _1^{wf}\). An equivalent model can also be fitted for a sample of unrelated individuals comprising of one sibling selected at random from each family (without including terms for family effects). This estimates the total contribution of the polygenic score to the phenotype (\(\alpha _1^{nwf}\)), by capturing both within-family and between-family effects. A difference in the magnitude of \(\alpha _1^{wf}\) vs. \(\alpha _1^{nwf}\) provides evidence for genetic nurture (and/or assortative mating and population stratification). Models were fitted in R v3.6.3, using the packages plm and ivmodel (code in Supplementary Methods).
Results
Table 1 presents the demographic characteristics of the two samples of siblings used to assess the role of genetic nurture in educational attainment and refractive error. Table 1 also presents the demographic characteristics of the split-sample GWAS samples used to derive polygenic scores for each trait. A much larger sample size was available for the educational attainment analyses compared to the refractive error analyses, as autorefraction was performed only at the latter stages of UK Biobank recruitment.
Contribution of genetic nurture to educational attainment
Split-sample GWAS analyses for EduYears were carried out in two non-overlapping samples of UK Biobank participants (Fig. 3). From the two sets of GWAS summary statistics, two independent polygenic scores for EduYears were derived. The variance in EduYears explained by each of these polygenic scores was assessed in a sample comprising of one randomly chosen sibling (sibling ‘A’) from the genetic nurture analysis sample for EduYears, which comprised of 18,146 sibling pairs who were unrelated to the participants in the GWAS samples. Each polygenic score explained approximately 8% of the variance in EduYears. Specifically, the incremental R2 for the first polygenic score PGS(m=1) was R2 = 0.077 and for the second polygenic score PGS(m=2) it was R2 = 0.080. The correlation of the two polygenic scores (PGS(m=1) vs. PGS(m=2) in sibling A) in this sample was 0.509 (95% CI 0.498–0.519).
For the first polygenic score PGS(m=1) the correlation of the polygenic score for sibling A vs. sibling B in each sibship, i.e. the within-sibling-pair correlation, was 0.556 (95% CI 0.546–0.566). For PGS(m=2) the corresponding within-sibling-pair correlation was 0.553 (95% CI 0.542–0.563). These within-sibling-pair correlations were above the level of 0.5 expected under random mating, suggesting the presence of assortative mating, as reported previously for the trait EduYears [41].
The results of the genetic nurture analysis for EduYears are shown in Table 2, Fig. 4 and Supplementary Table S1. For probands in a sample composed of one randomly selected sibling from each pair, the instrumented polygenic score yielded an estimate of the SNP-heritability of \(\alpha _1^{nwf}\)= 0.174 (95% CI 0.158–0.191). By contrast, a within-family analysis provided a lower estimate, \(\alpha _1^{wf}\) = 0.040 (95% CI 0.030–0.051). The relatively lower magnitude of \(\alpha _1^{wf}\) compared to \(\alpha _1^{nwf}\) (0.040 vs. 0.174) indicated that the contribution to the SNP-heritability of EduYears from genetic nurture and assortative mating combined was larger than the contribution from transmitted alleles. This is consistent with previous work, which has suggested that as much as half of the variance in educational attainment captured by a polygenic score acts via alleles present in the parents irrespective of transmission to the proband or reflects assortative mating.
Contribution of genetic nurture to refractive error
Following the same methodology described above for studying the contribution of genetic nurture to EduYears, split-sample GWAS analyses for avMSE were carried out in two further non-overlapping samples of UK Biobank participants (Fig. 3). From the two sets of GWAS summary statistics, two independent polygenic scores for avMSE were derived. The variance in avMSE explained by each of the polygenic scores was assessed in a sample comprising of one randomly chosen sibling from the genetic nurture analysis sample for avMSE, which comprised of 1944 sibling pairs. Each polygenic score for avMSE explained ~6% of the variance in avMSE in this independent sample. Specifically, the incremental R2 = 0.064 for PGS(m=1) and R2 = 0.059 for PGS(m=2) (n = 1944). The correlation of the two polygenic scores for avMSE (PGS(m=1) vs. PGS(m=2) in sibling A) in this sample of unrelated individuals was 0.316 (95% CI 0.275–0.355).
In the sample of n = 1944 sibling pairs with information on avMSE, the within-sibling-pair correlation of PGS(m=1) was 0.477 (95% CI 0.442–0.511) and the within-sibling-pair correlation for PGS(m=2) was 0.499 (95% CI 0.465–0.532). These correlations were close to the value of 0.5 expected under random mating, although the confidence intervals were wide. To provide more precise estimates, the within-sibling-pair analyses were repeated in the larger sample of siblings (n = 18,146 pairs) used for the genetic nurture analysis of EduYears. (Only genetic data were required to calculate the within-sibling correlation in polygenic scores for avMSE. Therefore, this calculation was possible even though information for the avMSE phenotype was mostly not available in this sample of 18,146 sibling pairs). Calculated for this sample of 18,146 sibling pairs, the within-sibling-pair correlation in polygenic scores for avMSE was 0.526 (95% CI 0.516–0.537) for PGS(m=1) and 0.525 (95% CI 0.514 –0.535) for PGS(m=2). These more precise estimates were indicative of assortative mating, although to a lesser degree than for EduYears.
An instrumented polygenic score for avMSE yielded a non-within-family SNP-heritability estimate of \(\alpha _1^{nwf}\)= 0.213 (95% CI 0.134–0.310) and a within-family estimate of \(\alpha _1^{wf}\)= 0.250 (95% CI 0.152–0.372). Notably, the within-family SNP-heritability estimate for avMSE was not lower than the non-within-family estimate (0.213 vs. 2.50; Table 2, Fig. 4 and Supplementary Table S2). This lack of a reduction in the within-family estimate suggested that genetic nurture (and assortative mating) made little contribution to the SNP-heritability of avMSE.
Discussion
Time spent outdoors, reading or playing video games is associated with the risk of myopia and is also correlated within siblings, suggesting a familial element to these behaviours [42]. Furthermore, a child’s refractive error PGS and the number of parents with myopia were reported to be independent predictors of myopia in children, indicating that parents may pass on myopia-predisposing environments as well as myopia-predisposing alleles to their children [43, 44]. Gene-environment correlation (expected in genetic nurture) has also been shown in the development of myopia [44]. A history of smoking by paternal grandmothers has been reported to be associated with early-onset myopia [45]. Epidemiological studies have consistently demonstrated a strong association between educational attainment and refractive error [46]. Indeed, Mendelian randomization and regression discontinuity analyses have suggested that this relationship may be causal [11, 12, 34]. These prior links between family environment, education and myopia led us to hypothesize that genetic nurture would play an important role in refractive error development, and that, if so, genetic nurture may upwardly bias estimates of the SNP-heritability for refractive error [47]. The key finding in this study was a lack of support for a contribution from genetic nurture to the genetics of refractive error. The lower bound of the 95% confidence interval for the within-family genetic contribution \(\alpha _1^{wf}\) was 0.152, which, in comparison to the point estimate of 0.250, suggests genetic nurture contributed at most 40% (0.250–0.152/0.250) of the SNP-heritability of refractive error in the UK Biobank sample. Analyses in larger samples will be needed to gain more precise estimates.
The above-chance correlation of polygenic scores within pairs of siblings suggested that assortative mating may be present for both educational attainment and refractive error, which supports past findings [41, 48, 49]. Our estimate of the SNP-heritability for educational attainment (\(h_{{{{{\rm{SNP}}}}}}^2\) ≈ 0.17) is in line with previous estimates [16, 50]. Our SNP-heritability estimate for refractive error (\(h_{{{{{\rm{SNP}}}}}}^2\) ≈ 0.22) is lower than previous estimates [7, 47], but the precision of the current assessment was limited by the small sample size. It was notable, in this study, that only about 25% of the SNP-heritability of educational attainment was attributed to alleles transmitted from parents to children (\(\alpha _1^{wf} = 0.04\) vs. \(\alpha _1^{nwf} = 0.17\)), with the rest being attributed to genetic nurture or assortative mating. Previous studies have suggested the contribution from transmitted alleles to educational attainment to be approximately 50% rather than 25% [16, 50]. One potential reason for this is that previous studies may not have fully accounted for attenuation bias resulting when performing a PGS-based mediation analysis [31], whereas we used an instrumented PGS to minimize this bias. Additional factors, such as differences between study populations and phenotype definitions, may also be relevant.
Studies of siblings offer distinct benefits compared to studying genetic risk in samples of unrelated individuals. Case-control association studies in unrelated samples rely on the principle that disease risk-increasing alleles will be over-represented in cases compared to controls. However, three other phenomena can also cause a significant over-representation of alleles in cases: population stratification, assortative mating and genetic nurture [40, 48]. Analyses of siblings, on the other hand, can be designed to be free from these biases [41]. In the current study, our finding that \(\alpha _1^{nwf}\) ≈ \(\alpha _1^{wf}\) for refractive error suggested that population stratification, assortative mating and genetic nurture made little contribution to the SNP-heritability of refractive error.
Strengths of the current work were the highly standardized method for measuring refractive error, use of powerful PGS constructs implemented in an instrumental variable framework [31, 33]. The main weakness of the current work was the limited sample available for the analysis of refractive error (n = 1944 pairs of siblings), leading to wide confidence intervals for the key parameters of \(\alpha _1^{wf}\)and \(\alpha _1^{nwf}\) and to insufficient power to examine genetic nurture effects for specific SNPs strongly associated with refractive error–as previously applied to investigate specific variants associated with birth weight [30]. Although autorefraction is considered the gold standard method for measuring refractive error in research studies [51], differences between autorefraction readings and manifest refraction would have contributed to the imprecision of the polygenic scores for refractive error.
When considering these new results, it is important to bear in mind that our analyses were conducted in a group of participants who grew up in the United Kingdom during the 1940s to 1980s. As such, our findings do not rule out a greater role for genetic nurture in more recent birth cohorts, or in communities with differing lifestyle-related risk factors, such as urban regions in East and South-East Asia where the prevalence of myopia is much higher than the United Kingdom [1]. Given these caveats, the current observation that genetic nurture makes a limited contribution to refractive error has several implications. First, it suggests that genetic variants identified in GWASs for refractive error are likely to exert their effects in probands rather in their parent, highlighting the genes that are potential therapeutic targets in children; this point has previously been assumed, but without any supporting evidence [5, 7]. Second, the evidence that genetic variants act directly in probands suggests that polygenic scores for predicting children at high risk of myopia may have greater clinical utility than they would otherwise. For instance, selecting a child’s treatment based on his/her PGS may be both more effective and more ethically justified, if genetic risk for myopia is intrinsic rather than a function of the environment. More widely, it will be of interest to investigate if genetic nurture contributes to other heritable ocular traits such as astigmatism and glaucoma [52, 53].
In summary, in this study, we obtained minimally-biased estimates of the genetic contribution to refractive error and educational attainment; specifically, free from bias due to genetic nurture and assortative mating. Building on prior work, we estimated that ~25% of the SNP-heritability of educational attainment was attributed to alleles transmitted from parents to children and that the remaining ~75% occurred as a result of genetic nurture and assortative mating. In stark contrast, our results suggested little contribution from genetic nurture to the SNP-heritability of refractive error. While we did find evidence of assortative mating for refractive error, there was no evidence that this appreciably inflated its SNP-heritability estimate. We conclude that, whereas the genetic contribution to educational attainment is in large part inherited through the environment, the genetic contribution to refractive error occurs mainly through direct parent-to-child transmission of alleles. Our findings validate the assumption that genetic variants associated with refractive error highlight potential therapeutic targets.
Data availability
Individual-level data from UK Biobank can be accessed by applying to the UK Biobank (https://www.ukbiobank.ac.uk/enable-your-research/apply-for-access).
Code availability
R code for implementing the statistical analysis is provided in Supplementary Methods.
References
Morgan IG, Wu P-C, Ostrin LA, Tideman JWL, Yam JC, Lan W, et al. IMI risk factors for myopia. Investig Ophthalmol Vis Sci. 2021;62:3.
Baird PN, Saw S-M, Lanca C, Guggenheim JA, Smith EL, Zhou X, et al. Myopia. Nat Rev Dis Prim. 2020;6:99.
Ohno-Matsui K, Wu P-C, Yamashiro K, Vutipongsatorn K, Fang Y, Cheung CMG, et al. IMI Pathologic Myopia. Investig Ophthalmol Vis Sci. 2021;62:5–5.
Sanfilippo PG, Hewitt AW, Hammond CJ, Mackey DA. The heritability of ocular traits. Surv Ophthalmol. 2010;55:561–83.
Kiefer AK, Tung JY, Do CB, Hinds DA, Mountain JL, Francke U, et al. Genome-wide analysis points to roles for extracellular matrix remodeling, the visual cycle, and neuronal development in myopia. PLoS Genet. 2013;9:e1003299.
Pickrell JK, Berisa T, Liu JZ, Segurel L, Tung JY, Hinds DA. Detection and interpretation of shared genetic influences on 42 human traits. Nat Genet. 2016;48:709–17.
Hysi PG, Choquet H, Khawaja AP, Wojciechowski R, Tedja MS, Yin J, et al. Meta-analysis of 542,934 subjects of European ancestry identifies new genes and mechanisms predisposing to refractive error and myopia. Nat Genet. 2020;52:401–7.
Wojciechowski R. Nature and nurture: the complex genetics of myopia and refractive error. Clin Genet. 2011;79:301–20.
Wu PC, Tsai CL, Wu HL, Yang YH, Kuo HK. Outdoor activity during class recess reduces myopia onset and progression in school children. Ophthalmology. 2013;120:1080–5.
He M, Xiang F, Zeng Y, Mai J, Chen Q, Zhang J, et al. Effect of time spent outdoors at school on the development of myopia among children in China: a randomized clinical trial. JAMA. 2015;314:1142–8.
Cuellar-Partida G, Lu Y, Kho PF, Hewitt AW, Wichmann HE, Yazar S, et al. Assessing the Genetic Predisposition of Education on Myopia: A Mendelian Randomization Study. Genet Epidemiol. 2016;40:66–72.
Mountjoy E, Davies NM, Plotnikov D, Davey Smith G, Rodriguez S, Williams CE, et al. Education and myopia: assessing the direction of causality by mendelian randomisation. BMJ. 2018;361:k2022.
Summers JA, Schaeffel F, Marcos S, Wu H, Tkatchenko AV. Functional integration of eye tissues and refractive eye development: Mechanisms and pathways. Exp Eye Res. 2021;209:108693.
Tkatchenko TV, Tkatchenko AV. Genome-wide analysis of retinal transcriptome reveals common genetic network underlying perception of contrast and optical defocus detection. BMC Med Genomics. 2021;14:153.
Jong M, Jonas JB, Wolffsohn JS, Berntsen DA, Cho P, Clarkson-Townsend D, et al. IMI 2021 Yearly Digest. Invest Ophthalmol Vis Sci. 2021;62:7.
Kong A, Thorleifsson G, Frigge ML, Vilhjalmsson BJ, Young AI, Thorgeirsson TE, et al. The nature of nurture: Effects of parental genotypes. Science. 2018;359:424–8.
Morris TT, Davies NM, Hemani G, Davey Smith G. Population phenomena inflate genetic associations of complex social traits. Sci Adv. 2020;6:eaay0328.
McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JP, et al. Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008;9:356–69.
Tubbs JD, Hwang LD, Luong J, Evans DM, Sham PC. Modeling Parent-Specific Genetic Nurture in Families with Missing Parental Genotypes: Application to Birthweight and BMI. Behav Genet. 2021;51:289–300.
Wu Y, Zhong X, Lin Y, Zhao Z, Chen J, Zheng B, et al. Estimating genetic nurture with summary statistics of multigenerational genome-wide association studies. Proc Natl Acad Sci USA. 2021;118:e2023184118.
Cawley J, Han E, Kim J, Norton EC. Testing for family influences on obesity: the role of genetic nurture. Health Econ. 2019;28:937–52.
Bates TC, Maher BS, Medland SE, McAloney K, Wright MJ, Hansell NK, et al. The nature of nurture: using a virtual-parent design to test parenting effects on children’s educational attainment in genotyped families. Twin Res Hum Genet. 2018;21:73–83.
Willoughby EA, McGue M, Iacono WG, Rustichini A, Lee JJ. The role of parental genotype in predicting offspring years of education: evidence for genetic nurture. Mol Psychiatry. 2019;26. https://doi.org/10.1038/s41380-41019-40494-41381.
de Zeeuw EL, Hottenga J-J, Ouwens KG, Dolan CV, Ehli EA, Davies GE, et al. Intergenerational transmission of education and ADHD: effects of parental genotypes. Behav Genet. 2020;50:221–32.
Young AI, Nehzati SM, Lee C, Benonisdottir S, Cesarini D, Benjamin DJ, et al. Mendelian imputation of parental genotypes for genome-wide estimation of direct and indirect genetic effects. bioRxiv. 2020. https://doi.org/10.1101/2020.07.02.185199.
Hwang LD, Tubbs JD, Luong J, Lundberg M, Moen GH, Wang G, et al. Estimating indirect parental genetic effects on offspring phenotypes using virtual parental genotypes derived from sibling and half sibling pairs. PLoS Genet. 2020;16:e1009154.
Conley D, Johnson R, Domingue B, Dawes C, Boardman J, Siegal M. A sibling method for identifying vQTLs. PLoS ONE. 2018;13:e0194541.
Juliusdottir T, Steinthorsdottir V, Stefansdottir L, Sveinbjornsson G, Ivarsdottir EV, Thorolfsdottir RB, et al. Distinction between the effects of parental and fetal genomes on fetal growth. Nat Genet. 2021;53:1135–42.
Okbay A, Beauchamp JP, Fontana MA, Lee JJ, Pers TH, Rietveld CA, et al. Genome-wide association study identifies 74 loci associated with educational attainment. Nature. 2016;533:539–42.
Warrington NM, Beaumont RN, Horikoshi M, Day FR, Helgeland Ø, Laurin C, et al. Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors. Nat Genet. 2019;51:804–14.
Becker J, Burik CAP, Goldman G, Wang N, Jayashankar H, Bennett M, et al. Resource profile and user guide of the Polygenic Index Repository. Nat Hum Behav. 2021;5. https://doi.org/10.1038/s41562-41021-01119-41563.
Tucker-Drob EM. Measurement error correction of genome-wide polygenic scores in prediction samples. bioRxiv. 2017. https://doi.org/10.1101/165472.
DiPrete TA, Burik CAP, Koellinger PD. Genetic instrumental variable regression: Explaining socioeconomic and health outcomes in nonexperimental data. Proc Natl Acad Sci USA. 2018;115:E4970–9.
Plotnikov D, Williams C, Atan D, Davies NM, Ghorbani Mojarrad N, Guggenheim JA, et al. Effect of education on Myopia: evidence from the United Kingdom ROSLA 1972 reform. Investig Ophthalmol Vis Sci. 2020;61:7.
Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK Biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12:e1001779.
Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562:203–9.
Bulik-Sullivan BK, Loh P-R, Finucane HK, Ripke S, Yang J. Schizophrenia Working Group of the Psychiatric Genomics C, et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet. 2015;47:291–5.
Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015;4:7.
Loh P-R, Tucker G, Bulik-Sullivan BK, Vilhjalmsson BJ, Finucane HK, Salem RM, et al. Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nat Genet. 2015;47:284–90.
Brumpton B, Sanderson E, Heilbron K, Hartwig FP, Harrison S, Vie G, et al. Avoiding dynastic, assortative mating, and population stratification biases in Mendelian randomization through within-family analyses. Nat. Commun. 2020;11:3519.
Selzam S, Ritchie SJ, Pingault J-B, Reynolds CA, O’Reilly PF, Plomin R. Comparing within- and between-family polygenic score prediction. Am J Hum Genet. 2019;105:351–63.
Guggenheim JA, Pong-Wong R, Haley CS, Gazzard G, Saw SM. Correlations in refractive errors between siblings in the Singapore Cohort Study of Risk-factors for Myopia. Br J Ophthalmol. 2007;91:781–4.
Ghorbani Mojarrad N, Williams C, Guggenheim JA. A genetic risk score and number of myopic parents independently predict myopia. Ophthalmic Physiol Opt. 2018;38:492–502.
Enthoven CA, Tideman JWL, Polling JR, Tedja MS, Raat H, Iglesias AI, et al. Interaction between lifestyle and genetic susceptibility in myopia: the Generation R study. Eur J Epidemiol. 2019;34:777–84.
Williams C, Suderman M, Guggenheim JA, Ellis G, Gregory S, Iles-Caven Y, et al. Grandmothers’ smoking in pregnancy is associated with a reduced prevalence of early-onset myopia. Sci Rep. 2019;9:15413.
Morgan I, Rose K. How genetic is school myopia? Prog Retin Eye Res. 2005;24:1–38.
Guggenheim JA, St Pourcain B, McMahon G, Timpson NJ, Evans DM, Williams C. Assumption-free estimation of the genetic contribution to refractive error across childhood. Mol Vis. 2015;21:621–32.
Robinson MR, Kleinman A, Graff M, Vinkhuyzen AAE, Couper D, Miller MB, et al. Genetic evidence of assortative mating in humans. Nat Hum Behav. 2017;1:0016.
Xiang F, He M, Morgan IG. The impact of parental myopia on myopia in chinese children: population-based evidence. Optom Vis Sci. 2012;89:1487–96.
Wang B, Baldwin JR, Schoeler T, Cheesman R, Barkhuizen W, Dudbridge F, et al. Robust genetic nurture effects on education: A systematic review and meta-analysis based on 38,654 families across 8 cohorts. Am J Hum Genet. 2021;108:1780–91.
Wolffsohn JS, Kollbaum PS, Berntsen DA, Atchison DA, Benavente A, Bradley A, et al. IMI - clinical myopia control trials and instrumentation report. Investig Ophthalmol Vis Sci. 2019;60:M132–60.
Parssinen O, Kauppinen M, Kaprio J, Koskenvuo M, Rantanen T. Heritability of refractive astigmatism: a population-based twin study among 63- to 75-year-old female twins. Investig Ophthalmol Vis Sci. 2013;54:6063–7.
Springelkamp H, Iglesias AI, Mishra A, Höhn R, Wojciechowski R, Khawaja AP, et al. New insights into the genetics of primary open-angle glaucoma based on meta-analyses of intraocular pressure and optic disc characteristics. Hum Mol Genet. 2017;26:438–53.
Acknowledgements
This research has been conducted using the UK Biobank Resource (application #17351). UK Biobank was established by the Wellcome Trust; the UK Medical Research Council; the Department for Health (London, UK); Scottish Government (Edinburgh, UK); and the Northwest Regional Development Agency (Warrington, UK). It also received funding from the Welsh Assembly Government (Cardiff, UK); the British Heart Foundation; and Diabetes UK. Collection of eye and vision data was supported by The Department for Health through an award made by the NIHR to the Biomedical Research Centre at Moorfields Eye Hospital NHS Foundation Trust, and UCL Institute of Ophthalmology, London, United Kingdom (grant no. BRC2_009). Additional support was provided by The Special Trustees of Moorfields Eye Hospital, London, United Kingdom (grant no. ST 12 09). Data analysis was carried out using the HAWK computing cluster, maintained by Supercomputing Wales and Cardiff University ARCCA. This work was supported by the Welsh Government and Fight for Sight (Ref: 24WG201) to JAG and CW; NIH Grant R37MH107649–07S1 and by Research Council of Norway (NRC) Grant No. 288083 to TZ. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Author information
Authors and Affiliations
Consortia
Contributions
JAG and CW designed the research and obtained funding. JAG and RC performed the research. All authors contributed to interpretation of the results and drafting of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Ethical approval
Ethical approval for the study was obtained from the NHS Research Ethics Committee (Reference: 11/NW/0382). All participants provided signed, informed consent.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Guggenheim, J.A., Clark, R., Zayats, T. et al. Assessing the contribution of genetic nurture to refractive error. Eur J Hum Genet 30, 1226–1232 (2022). https://doi.org/10.1038/s41431-022-01126-6
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41431-022-01126-6