Cholesky model for axial length (AL) and height. Measured phenotypes for twin 1 and twin 2 are shown as rectangles, and latent additive genetic and random environmental factors are shown as circles. A indicates additive genetic factor; a, additive genetic determinants of AL and height; E, random environmental factor; e, additive environmental determinants of AL and height; and 1.0, additive genetic correlation in monozygotic twins; and 0.5, additive genetic correlation in dizygotic twins.
Standardized variable estimates (95% confidence intervals) in the best-fitting Cholesky model for axial length (AL) and height. Measured phenotypes are shown as rectangles, and latent additive genetic and random environmental factors are shown as circles. A indicates additive genetic factor; E, random environmental factor.
Zhang J, Hur Y, Huang W, Ding X, Feng K, He M. Shared Genetic Determinants of Axial Length and Height in ChildrenThe Guangzhou Twin Eye Study. Arch Ophthalmol. 2011;129(1):63-68. doi:10.1001/archophthalmol.2010.323
To describe the association between axial length (AL) and height and to estimate the extent to which shared genetic or environmental factors influence this covariance.
Study participants were recruited from the Guangzhou Twin Registry. Axial length was measured using partial coherence laser interferometry. Height was measured with the participants standing without shoes. We computed twin pairwise correlations and cross-twin cross-trait correlations between AL and height for monozygotic and dizygotic twins and performed model-fitting analyses using a multivariate Cholesky model. The right eye was arbitrarily selected to represent AL of participants.
Five hundred sixty-five twin pairs (359 monozygotic and 206 dizygotic) aged 7 to 15 years were available for analysis. Phenotypic correlation between AL and height was 0.46 but decreased to 0.19 after adjusting for age, sex, and age × sex interaction. Bivariate Cholesky model–fitting analyses revealed that 89% of phenotypic correlation was due to shared genetic factors and 11% was due to shared random environmental factors, which includes measurement error.
Covariance of AL and height is largely attributable to shared genes. Given that AL is a key determinant of myopia, further work is needed to confirm gene sharing between myopia and stature.
Results of population-based studies1- 9 have suggested that taller adults and children have greater axial length (AL), although the association with refraction was inconsistent. The eye globe is increased 0.25 mm per 0.1-mm increase in height among adults.5 Significant positive correlations between eye size growth and overall physical development have also been reported.10- 12 The fact that cessation in growth of AL and height may occur at a similar age6 suggests that these 2 traits may share similar developmental pathways. Although no convincing explanation for these associations has been proposed to date, 2 previous studies13,14 found that correlations between height and refraction remained significant after data were corrected for socioeconomic status variables, such as education and income, suggesting that stature may be a risk factor for myopic development independent of socioeconomic status. Given that stature15 and myopia16,17 are significant genetically influenced traits, we speculate that shared genes or shared environmental factors may explain the association between height and refractive error.
Refraction is a complex phenotype that involves various biometric variables of the eye, including AL, corneal curvature, anterior chamber depth, lens thickness, vitreous chamber depth, and refractive power of the structure. Among these, AL has consistently been shown to be a genetically influenced trait that is closely related to myopia. Therefore, we chose AL in the present study to explore underlying genetic relationships between AL and stature. To the best of our knowledge, this issue has not been investigated before.
Monozygotic (MZ) twins are genetically identical, whereas dizygotic (DZ) twins have, on average, 50% of their genes in common. Quantitative genetic analysis of twins can be used to estimate genetic and environmental contributions to trait variations. In addition, biometrical models based on twin data have been used to estimate shared genetic determinants of multiple-trait covariance. Cross-twin cross-trait correlations, defined as correlations between trait A in one twin and trait B in the other twin, can be compared in MZ and DZ twins to investigate whether shared genetic or environmental factors influence the association between traits A and B.
In the present study, our objective was to assess the association between AL and height. We further aimed to estimate the extent to which shared genetic or environmental factors influence this covariance.
Study participants were recruited from the Guangzhou Twin Registry, which has been described in depth elsewhere.18 In brief, this registry was established in Guangzhou City, China, in 2005-2006. All twins born between January 1, 1987, and December 31, 2000, were identified using the official Household Registry of Guangzhou, followed by door-to-door verification. In July and August 2006, we invited twins aged 7 to 15 years (as of July 1, 2006) living in 2 districts for baseline data collection. Among 705 pairs invited, 565 pairs (359 MZ and 206 DZ twins) participated in the study, yielding a participation rate of approximately 80%. For all invited participants, written informed consent was obtained from parents or legal guardians of the twins following an in-depth explanation of the study. Ethical approval was obtained from the Sun Yat-sen University ethical review board and the Zhongshan Ophthalmic Center ethics committee, and the study was conducted in accord with the tenets of the World Medical Association's Declaration of Helsinki.
In this twin cohort, zygosity of all same-sex twin pairs was determined by 16 multiplex short tandem repeats (PowerPlex 16; Promega, Madison, Wisconsin)19 at the Forensic Medicine Department of Sun Yat-sen University in 2006. Opposite-sex twin pairs were considered DZ without genotyping.
Axial length was measured using a noncontact partial coherence laser interferometry technique (IOLMaster; Carl Zeiss Meditec, Oberkochen, Germany) in a dark room (illumination, <5 lux) before pharmacologic dilation of pupils. The mean of 5 continuous measurements was used. Measurements with a ratio of signal to noise less than 2.0 (displayed as “Borderline SNR” or “Error”) or with more than a 0.1-mm difference from the others (displayed as “Evaluation”) were deleted and remeasured. The right eye was arbitrarily selected to represent phenotypic characteristics of participants in the data analysis. Height was measured with the participants standing without shoes following a standard protocol. The median values of 3 consecutive readings were then recorded to the nearest 0.1 cm.
To assess gene sharing between AL and height, we first examined twin pairwise correlations and cross-twin cross-trait correlations for MZ and DZ twins. We then performed model-fitting analyses using a multivariate Cholesky model.
In the classic twin design, the following 4 components comprise the total phenotype variance: additive genetic (A), dominance genetic (D), common environmental (C), and random environmental (E), which includes measurement error. Given the confounding effects of C and D, it is inappropriate to put these 2 factors in the same model. Therefore, we fit ADE and ACE models separately.
The bivariate Cholesky model is shown in Figure 1. Measured phenotypes for twin 1 and twin 2 are shown as rectangles, and latent additive genetic and random environmental factors are shown as circles. The first additive genetic (A1) and random environmental (E1) factors exert influences on AL and height. The second Cholesky factors (ie, A2 and E2) are those unique to AL. Although ADE and ACE models were tested, C and D are not shown in the figure for ease of interpretation. The a21 path reflects shared additive genetic determinants of AL and height. The A, C/D, and E variance and covariance matrices were computed by the product of their respective Cholesky factor-loading matrix and its transposition. Genetic and environmental correlations, derived from variances and covariances, allowed us to determine the percentages of phenotypic correlation attributable to shared genetic factors (ie, bivariate heritability) and to shared environmental factors (ie, bivariate environmentality). In addition to bivariate heritability and environmentality, we computed genetic and environmental variances unique to AL, as well as those shared with height.
We used the Mx program (Virginia Commonwealth University, Richmond) to fit the Cholesky model to variance and covariance matrices computed from MZ and DZ twins.20 Akaike information criterion (χ2 − 2 df) and likelihood ratio χ2 test were used to select the best-fitting model. Lower Akaike information criterion was considered a better fit. For likelihood ratio χ2 test, a significant χ2 increase in the reduced model compared with the full model suggests that the reduced model fit the data worse than the full model. An insignificant χ2 change indicates that the reduced model is acceptable. Before model-fitting analyses, AL and height were adjusted for age, sex, and age × sex interaction. Predicted values were calculated for each observation based on a linear regression model that included age, sex, and their interaction as independent variables. Residual values were then calculated as differences between observed and predicted values, and these standardized residuals were then used in variance and covariance matrices for heritability estimation.
The study included 565 twin pairs (359 MZ and 206 DZ), with AL and height data available for all. The mean (SD) age was 11.7 (2.5) years, and 50.6% were female. The mean (SD) AL and height in this cohort were 23.50 (1.10) mm (range, 20.27-27.11 mm) and 144.3 (14.6) cm (range, 106.5-183.4 cm), respectively. Estimated regression coefficients for AL and height with age, sex, birth order, and zygosity are given in Table 1. Multivariate linear regression models revealed that age and sex were statistically significant for AL (R2 = 0.24) and height (R2 = 0.76). However, age × sex interaction was significant for height but not for AL; therefore, correlation coefficients were adjusted for this interaction for AL and height. The regression coefficient for birth order was not significant for AL (−0.03, P = .57) or for height (−0.05, P = .90). For AL and height, no significant differences in regression coefficients were found between MZ and DZ twins, fulfilling the basic assumptions for twin analyses.
Table 2 gives phenotypic correlations, twin correlations, and cross-twin cross-trait correlations between AL and height. Unadjusted phenotypic correlation was 0.46, which was highly statistically significant. This correlation remained significant when it decreased to 0.19 after adjusting for age, sex, and their interaction for AL and height. Monozygotic twin correlations exceeded DZ twin correlations for both phenotypes. Cross-twin cross-trait correlations for AL and height were similar in MZ and DZ twins.
Table 3 gives results of model-fitting analyses. Although ADE and ACE models fit the data well, we selected ADE as the Cholesky full model because the Akaike information criterion was lower in ADE than in ACE. Results suggested that the best-fitting model was model 3, in which additive genetic and random environmental factors exert significant influences on AL and height variations and covariance.
Figure 2 shows standardized Cholesky factor loadings and their 95% confidence intervals in the best-fitting model. Estimates of shared genetic factors computed from the best-fitting model were 89% for AL and 92% for height, confirming substantial genetic influences on these traits.15,21 Estimates of shared random environmental factors were 11% for AL and 8% for height. Correlations between AL and height were 0.19 for genetic factors and 0.22 for random environmental factors. Bivariate heritability (ie, proportion of phenotypic correlation between AL and height [r = 0.19] attributable to shared genetic factors) was 89% ([0.19 × √0.92 × √0.89] / 0.19). Bivariate environmentality (ie, proportion of phenotypic correlation between AL and height attributable to shared random environmental factors) was 11% ([0.22 × √0.08 × √0.11] / 0.19). These results suggested that phenotypic correlation between AL and height, although modest (r = 0.19), was largely determined by genetic rather than environmental commonality.
Recently, associations between ocular variables and anthropometric measurements of height, weight, and body mass index have received considerable attention.1- 5 These investigations were unable to quantify shared genetic or environmental contributions to associations among phenotypes (ie, AL, height, and refraction) because results were primarily based on unrelated singletons. To our knowledge, the present study is the first to explore shared genetic and environmental determinants of associations between AL and height in a population-based sample of young Chinese twins.
The association between stature and ocular biometry observed herein is consistent with some findings and inconsistent with others.1- 5 In a Finnish study8 and a Danish study,22 men with myopia were taller than those without myopia, a difference that was not observed among the women. Significant associations between height and refraction were also found in children.4,9 Eye globe and body stature are in developmental phases in children, whereas they are stable in adults. Lens opacity may also alter refraction in adults. Although the age range of our twins is large, we conducted model-fitting analyses based on the total sample because such analyses require a large sample. However, when we divided our total sample into 2 age groups (<10 years and ≥10 years) and computed correlations between AL and height in each, no significant difference in correlation was found between the 2 groups, suggesting that the pattern of association between AL and height is similar in younger and older samples in our study. Nevertheless, it will be important to replicate our findings in adults and in large population samples of children with a smaller age range.
In our study, we found a significant association between AL and height (r = 0.45). When we adjusted data for age and sex, the correlation decreased (r = 0.19) but remained significant. Further modeling results suggest that 89% of adjusted correlation between AL and height (r = 0.19) was due to shared genetic factors, and the remaining 11% was due to shared environmental factors. In contrast, analyses of genetic variance of AL suggest that only about 3% was in common with height, and the remaining genetic variance was unique to AL. Axial length and height are likely determined by many genes with small effect sizes. Our findings suggest that pleiotropic effects of these genes may influence the relationship between eye size and stature.
The present study is limited because our conclusions were derived from model-fitting analyses rather than from gene identification studies. Model-fitting analyses test hypotheses indirectly; therefore, only inferences of a genetic association can be made by comparing alternative models. However, the Online Mendelian Inheritance in Man database23 shows overlap of chromosome regions between loci identified for myopia and stature, supporting our results of shared genetic associations between AL and height. Among loci identified for myopia (MYP1 through MYP14), MYP6 and MYP14 were identified by genome-wide linkage in common myopia among Jewish families.24,25MYP7 through MYP10 were identified by genome-wide linkage in common myopia among DZ twins.26 The remaining loci were mainly identified by linkage analysis in familial high -myopia cases. Based on a twin cohort, a recent genome-wide linkage study27 identified the quantitative trait locus for AL on chromosome 5. There have been many reports of quantitative trait loci for stature.28 Notably, STQTL6 (Xq24-q25) and MYP13 (Xq23-q25) overlapped, whereas STQTL10 (3q23) and MYP8 (3q26), STQTL (7q31.3) and MYP4 (7q36), and STQTL12 (4q28-q32) and MYP11 (4q22-q27) were close to each other. Given that GPR126,29HMGA2,30ZBTB38,31CDK6,31LCORL,31UQCC,32 and GDF532 genes have been recently identified within quantitative trait loci for stature, the finding of significant shared genes for AL and height in the present study suggests that identified genes for stature may be considered candidate genes for myopia as well.
Axial length may be considered a possible endophenotype for myopia. However, myopia is a complex phenotype involving various biometric variables, and AL is only one of these components. It is possible that genetic and environmental factors related to other myopia-related variables such as corneal refraction and lens refraction differ for AL. Therefore, the genetic covariance identified herein may be applicable only to the association between AL and height, not between myopia and height. On the other hand, it is recognized that myopia may develop with reduction of corneal and lens refractive power in tandem with axial elongation during early development of emmetropia.33 Results of a population-based study34 among 12-year-old Australian children suggested that AL accounted for about half of the refraction variation. A 2008 twin study35 estimated that 50% of genetic covariance in refraction was due to genetic factors influencing AL. Therefore, AL represents one of the most important ocular dimensions related to development of refractive error.
In summary, we confirm that there is a significant association between AL and height in children and that most of this correlation is attributable to genetic covariance between the 2 traits. Given the existence of significant shared genetic factors, we may need to consider pleiotrophic effects of myopia and height in future quantitative trait locus–based association studies and gene-searching efforts.
Correspondence: Mingguang He, MD, PhD, State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou 510060, China (email@example.com).
Submitted for Publication: December 10, 2008; final revision received December 14, 2009; accepted January 24, 2010.
Financial Disclosure: None reported.
Funding/Support: This study was supported by grant 30772393 from the National Natural Science Foundation of China (Dr He); by grant NCET-06-0720 from the Program for New Century Excellent Talents in University, National Ministry of Education (Dr He); by special fund 30371513 from the Sino-Australia Chinese State National Science Fund Committee–Australian Department of Education, Science, and Training (Dr He); and by grant 2006Z3-E0061 from the Guangzhou Science and Technology Development Fund (Dr He).