|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1From the State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China; the 2UCL Institute of Ophthalmology, London, United Kingdom; and the 3Chonnam National University, Gwangju, South Korea.
| Abstract |
|---|
|
|
|---|
METHODS. The study participants were recruited from the Guangzhou Twin Registry. Anterior segment optical coherence tomography and custom software were used to quantify the angle opening distance (AOD) at the location 500 µm anterior to the scleral spur. Anterior chamber depth (ACD) and axial length (AL) were measured using laser interferometry. Cross-trait, cross-twin correlations for monozygotic (MZ) and dizygotic (DZ) twins and the Cholesky model were used to quantify shared genetic and environmental effects for AL, ACD, and AOD after adjusting for age and sex.
RESULTS. A group of 459 pairs of twins (304 MZ and 155 DZ) aged 8 to 16 years were available for analysis. The phenotypic correlations among AL, ACD, and AOD ranged from 0.39 to 0.64. Cross-twin, cross-trait correlations for these three phenotypes for MZ twins were consistently greater than the corresponding correlations for DZ twins. The results of the Cholesky model-fitting analyses can be summarized as follows: first, of 70% of additive genetic factors for AOD, 23% and 13% were those shared with ACD and AL, respectively, whereas the remaining 34% were those unique to AOD. Second, of 89% of additive genetic factors for ACD, 25% were those shared with AL, whereas 64% were those unique to ACD. Third, random environmental influences on covariances among AL, ACD, and AOD were very small.
CONCLUSIONS. Analyses of Chinese children twin data suggest that shared genes are responsible for the significant phenotypic correlations found for AL ACD, and AOD.
Twin studies are widely used to disentangle genetic and environmental contributions to variations in phenotypes. Monozygotic (MZ) twins are genetically identical, whereas dizygotic (DZ) twins have, on average, 50% of their genes in common. The heritability—the proportion of the total phenotypic variance attributable to genetic effects—can be estimated by the use of a univariate model. Furthermore, Cholesky models allow the estimation of genetic and environmental effects on the phenotypic associations among traits.11 The MZ and DZ twin correlations between trait A in one twin and trait B in his/her co-twin (termed cross-twin, cross-trait correlations) provide information on the extent to which genetic and environmental factors contribute to the phenotypic correlations between traits A and B. If genetic factors influence the association between the two traits, the MZ cross-twin, cross-trait correlation is expected to be larger than the DZ cross-twin, cross-trait correlation. On the other hand, if shared environmental factors are important for the relationship between the two traits, then the DZ cross-twin, cross-trait correlation is expected to be larger than
the MZ cross-twin, cross-trait correlation. If the association between the two traits is purely due to random environmental factors, then the MZ and DZ cross-twin, cross-trait correlations would not be significantly different from 0.
The use of intermediate phenotype (also called endophenotype) has been increasingly adopted for familial aggregation and gene mapping of complex diseases.12 Understanding the etiology of intermediate phenotypes may provide critical evidence for genetic and environmental pathways to the complex diseases of interest. As the anterior chamber depth (ACD) and drainage angle width (quantified as angle opening distance [AOD] in anterior segment optical coherence tomography [ASOCT]) are the anatomic basis of angle closure, they were chosen as intermediate phenotypes for angle closure in the present study. Axial length has been considered to be a phenotype closely correlated with myopia. Thus, AL, ACD, and AOD were selected as endophenotypes of myopia and angle closure in the present study. Our previous studies based on the Guangzhou Twin Registry confirmed that ACD, AL, and AOD are all highly heritable traits.13 14 We have shown that additive genetic effects explain 70% to 90% of the total variances of the three phenotypes. Given the substantial phenotypic associations among ACD, AOD, and AL and the high heritability observed for these traits, it is possible that shared genes may be responsible for the phenotypic covariation among ACD, AOD, and AL. The purposes of the present study were to determine whether common genetic and/or environmental effects influence the relationships among AL, ACD, and AOD and to estimate these common genetic and environmental effects.
| Materials and Methods |
|---|
|
|
|---|
In this twin cohort, the zygosity of all same-sex twin pairs was determined by 16 multiplex short-tandem repeats (STRs; PowerPlex 16 system, Promega, Madison, WI)16 at the Forensic Medicine Department of Sun Yat-Sen University in 2006. Opposite-sex twin pairs were assigned to be dizygotic (DZ) without genotyping.
Twin pairs were excluded from data analyses if one or both twins had pathologic changes (e.g., retinopathy of prematurity), recent orthokeratology contact lenses correction, or previous laser treatment for myopia.
Examination and Measurement
The right eye was arbitrarily selected to represent the phenotypic characteristics of the specific individual in the data analysis.
AL was measured by noncontact partial-coherence laser interferometry (IOLMaster; Carl Zeiss Meditec, Oberkochen, Germany) in a dark room (<5 lux illumination) before pharmacologic dilation of the pupils. The mean result of 10 continuous measurements was used. Poor measurements with signal-to-noise ratio (SNR) < 2.0 (displayed as "Borderline SNR," or "Error") or measurements with one result differing by more than 0.1 mm from the others (displayed as "Evaluation!") were deleted and remeasured. ACD was measured with the same interferometry device. The ACD measurement was made after automatic adjusting for the corneal radius. A mean of five measurements was taken. Measurements with displayed the "Error" message were deleted and remeasured. The ACD was measured as the distance between the anterior lens surface and corneal epithelium illumination and therefore included the thickness of the cornea.
ASOCT imaging (Visante; Carl Zeiss Meditec, Dublin, CA) was performed with the participants in a seated position without pupil dilation in a dark room with the same illumination as for the ACD and AL measurements. One scan, centered over the pupil, was taken on the horizontal meridian (between 0° and 180°). The scan direction was properly aligned until a full corneal reflex was achieved (indicated by an interference flare on the axis of anterior chamber). The fixation angle was adjusted to align and enable the image horizontal. One image with best quality was selected and saved for each eye of each person.
The ASOCT images were subsequently analyzed with custom software, the Zhongshan Angle Assessment Program (ZAAP, Guangzhou, China). This program performed noise and contrast reconditioning and allowed automatic measurement after the scleral spurs were identified. AOD, calculated at 500 µm from the scleral spur, was used in the present analysis. AOD in the present study was defined as the length of a line drawn from the anterior iris to the corneal endothelium perpendicular to a line along the trabecular meshwork 500 µm anterior to the scleral spur.17
Data Analysis and Genetic Modeling
To determine the causes of covariances among AL, ACD, and AOD, we computed twin correlations and cross-trait, cross-twin correlations for MZ and DZ twins and performed model-fitting analyses by using the Cholesky model.
As variances can be divided into genetic and environmental sources, covariances can be decomposed into additive genetic (A), dominance genetic (D), and random environmental (E) factors.11 The A factors represent the sum of the average effect of all alleles that influence a trait. The A factors correlate at 1.0 for MZ and at 0.5 for DZ twins. The D factors reflect the genetic factors that do not add up across alleles because of interaction. The D factors correlate at 1.0 for MZ and at 0.25 for DZ twins. Finally, the E factors, environmental factors that are unique to each member of a twin pair (e.g., accident, virus infection) and measurement error, do not contribute to the twin similarity. The shared family environmental factors were not considered in the present study, because the patterns of MZ and DZ twin correlations for AL, ACD, and AOD suggested that shared environmental factors are negligible (Table 1) , and because our previous studies have shown no significant shared environmental influences on AL, ACD, and AOD.13 14
|
|
2 – 2 df) and the likelihood-ratio
2 test (LRT) to select the best-fitting model. If competing models were not nested, the model that produced the lowest AIC was considered the best-fitting model. If alternative models were nested, then the LRT was used because the change in
2 is distributed as
2. A significant increase in
2 in the reduced model compared with the full model would suggest that the reduced model fit the data less well than the full model. A nonsignificant change in
2 would indicate that the reduction of the model parameter is acceptable. Our model-fitting analyses were based on observed variance and covariance matrices computed separately for MZ and DZ twins among AL, ACD, and AOD. We used Mx to fit the Cholesky decomposition model to the observed variance and covariance matrices.18
| Results |
|---|
|
|
|---|
Correlational Analyses
Table 1 presents phenotypic correlations, twin correlations (the elements in the rectangles on the diagonals), and cross-twin, cross-trait correlations (the off-diagonal elements) for AL, ACD, and AOD. All correlations were corrected for sex and age. The phenotypic correlations among AL, ACD, and AOD were high and significant, ranging from 0.39 to 0.64. MZ twin correlations were significantly greater than DZ twin correlations for all three phenotypes, confirming substantial genetic influences and no shared environmental influences on AL, ACD, and AOD. Cross-twin, cross-trait correlations in MZ were all greater than the corresponding DZ correlations, indicating that shared genetic factors played an important role in covariation among the three phenotypes. The impressions gained from these correlational results were tested using model-fitting analyses.
Model-Fitting Analyses
Cholesky Model-Fitting.
Table 2 provides the results of model-fitting analyses. The full Cholesky model fit the data well (
224 = 23.04, P = 0.52). When we removed nonadditive genetic variances and covariances among AL, ACD, and AOD from the full model, a nonsignificant change in
2 occurred (model 2; 
26 = 0.9, P = 0.99), suggesting that genetic effects on the three phenotypes are largely of the additive rather than the nonadditive kind. Eliminating additive genetic variances unique to AL, ACD, and AOD, respectively, significantly worsened the fit (models, 3, 4, and 5). Next, we removed additive genetic covariances between AL and ACD, between AL and AOD, and between ACD and AOD, respectively (models 6, 7, and 8). These procedures also produced significant changes in
2, suggesting the importance of common genetic effects on the covariances of AL, ACD, and AOD. Finally, random environmental covariances between AL and ACD, between AL and AOD, and between ACD and AOD were dropped from model 2, which yielded significant changes in
2 (models 9, 10, and 11). We did not attempt to remove random environmental variances unique to AL, ACD, and AOD, as these variances were confounded with measurement error. Taken together, these results suggested that the best-fitting model was model 2 where additive genetic and random environmental factors exert significant influences on variation and covariation among AL, ACD, and AOD.
|
|
Random environmental influences were modest: 12%, 11%, and 30% for AL, ACD, and AOD, respectively. When we decomposed these estimates, of 30% of random environmental influences on AOD, 26% were those unique to AOD (E3 in Fig. 2 ), 3% were those shared with ACD (E2 in Fig. 2 ), and 1% were those in common with AL (E1 in Fig. 2 ). Likewise, 11% of random environmental influences on ACD consisted of 2% of those shared with AL (E1 in Fig. 2 ) and 9% of those unique to ACD (E2 in Fig. 2 ).
| Discussion |
|---|
|
|
|---|
In our study, as we expected, we found high phenotypic associations among AL, ACD, and AOD. ACD and AL are both the distance measurement on the axial direction of the eye globe. AL is equivalent to a sum of ACD, lens thickness, and vitreous length and therefore, is closely associated with ACD. Using a case control design, Lowe5 showed that 35% of the variance of ACD was attributable to lens thickness and the remaining 65%, to the lens position. Population studies have consistently suggested that the anterior chamber tends to be deeper in myopic eyes.6 On the other hand, angle width, usually quantified using either gonioscopy or even better by ASOCT, is an estimation of the proximity between the peripheral iris and trabecular meshwork in the drainage angles of the eye. It is associated with ACD but also reflects a combination of other anatomic components, including iris thickness, iris curvature, and corneal curvature.
In most quantitative traits, genes are known to influence multiple phenotypes (pleiotropic effects of genes). The cross-twin, cross-trait correlations for AL, ACD, and AOD in the MZ pairs were greater than those in the DZ pairs, suggesting the importance of common genetic effects on the relationship among the three phenotypes. Cholesky model-fitting analyses confirmed the results of twin correlational analyses: Shared genetic variances for AL, ACD, and AOD were significant and high, whereas shared random environmental factors were very modest. Taken together, these results suggest that pleiotropic effects of genes may be operating in AL, ACD, and AOD and that the pleiotropic actions of genes probably contribute to the associations between angle closure and myopia-related traits.
Molecular genetic findings on angle closure and myopia remain limited and controversial. Pedigree-based studies have reported linkage loci on chromosome 11 for nanophthalmos but were not able to replicate linkage for angle-closure glaucoma or narrow angles.19 One recent study reported the MMP-9 gene to be associated with acute angle closure but yet to be confirmed.20 None of these angle closure–related chromosome regions have been reported to be overlapped with the loci related to myopia, although the linkage loci of myopia are mainly for high myopia and from pedigree-based studies.21 However, one has to be cautious that both angle closure and common myopia are likely polygenic disorders, in which the effect size at any single locus is likely to be small. The pedigree study or even population- and family-based association studies may lack the power to detect genes of such small effect size. Therefore, it would not be surprising, even if no chromosomal overlap has been detected in previous studies. However, the association between angle closure and myopia is biologically plausible. Myopia is characterized as excessive growth of the eye globe (longer eye), whereas angle closure occurs more often in shorter eyes, presumably representing no growth of the eye globe. Therefore, it is possible that there is a common set of genes underlying these two phenotypes that express in different pathways
In this study, the ACD and AOD are both quantitative traits and chosen as the intermediate phenotypes for the angle closure. However, these two traits are also closely associated with myopia—myopes tend to have longer eyes and deeper anterior chambers. Therefore, one may argue that the ACD and AOD may be phenotypes related to both angle closure and myopia. In our previous study,14 we divided the ACD by AL in an attempt to control the influence of myopia on the ACD trait. However, the heritability of ACD remained high even with adjustment for myopic effects. This finding underscores the need to consider the effect of myopia when we intend to use ACD or AOD as quantitative trait loci (QTL) traits of interest for the molecular genetic study of angle closure.
As our study participants were healthy young people, angle-closure glaucoma was extremely uncommon in the sample. Thus, the results of the present study are only applicable to normal children in China. Further inference on the etiology of angle closure or angle-closure glaucoma may depend on the validity of the assumption that shallow ACD in childhood indicates a propensity to a shallow ACD in later life. However, intrapair environmental factors are likely to be similar in young people, whereas adult or elderly twins may incur more diversified intrapair environmental influences. This may confer advantages for using relatively young twin for the heritability study.
In summary, we confirm strong, shared additive genetic effects on the traits related to angle closure and myopia. Only very small proportions of the covariances were explained by random environmental factors including measurement errors. This study underscores the need to consider the possible influence of pleiotropic effects of genes in angle closure and myopia in QTL-based association studies.
| Footnotes |
|---|
Submitted for publication April 4, 2008; revised June 6, 2008; accepted August 27, 2008.
Disclosure: M. He, None; Y.-M. Hur, None; J. Zhang, None; X. Ding, None; W. Huang, None; D. Wang, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Mingguang He, Department of Preventive Ophthalmology, Zhongshan Ophthalmic Center, Guangzhou 510060, Peoples Republic of China; mingguang_he{at}yahoo.com.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |