|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1From the State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China; the 2Department of Psychology, Chonnam National University, Gwangju, Korea; and the 3UCL Institute of Ophthalmology, London, United Kingdom.
| Abstract |
|---|
|
|
|---|
METHODS. Twins aged 8 to 16 years were recruited from the Guangzhou Twin Registry. Pachymetry data were obtained by one operator using the same imaging system. Zygosity was confirmed by genotyping with 16 polymorphic markers in all same-sex twin pairs. The CCT of the right eyes was chosen as the trait of interest in the analysis. Heritability was assessed by a general sex-limitation model, using Mx software (University of Richmond, Virginia).
RESULTS. Four hundred forty-nine twin pairs were available for data analyses, including 131 pairs of monozygotic boys (MZM), 44 pairs of dizygotic boys (DZM), 166 pairs of monozygotic girls (MZF), 31 pairs of dizygotic girls (DZF), and 77 pairs of opposite-sex dizygotic (OSDZ) twins. Twin correlations for CCT were 0.90 for MZM, 0.92 for MZF, 0.56 for DZM, 0.61 for DZF, and 0.44 for OSDZ twins. A sex-limitation model combining additive genetic and unique environmental factors produced the best fit for the data. Heritability estimates for CCT were 0.88 (95% confidence interval [CI]: 0.84–0.91) in the boys and 0.91 (95% CI: 0.89–0.93) in the girls. Unique environmental effects explained only 0.12 (95% CI: 0.09–0.16) and 0.09 (95% CI: 0.07–0.11) of the variance in CCT in the boys and the girls, respectively.
CONCLUSIONS. Additive genetic effects appear to be the major contributor to the variation of CCT in Chinese population. Heritability of CCT appears to be slightly greater in the girls than in the boys.
To be functional as a validated endophenotype, the trait should be heritable.8 In the context of POAG, heritability estimates of several quantitative traits have been reported, such as intraocular pressure (IOP),10 11 optic disc parameters,10 and retinal nerve fiber layer (RNFL) thickness.11 Central corneal thickness (CCT) has been considered to be a powerful risk factor and thus has been hypothesized to be another quantitative trait (endophenotype) for POAG.12 13 To our knowledge, the heritability of CCT in European people has been reported in only one study.12 Population studies have consistently suggested ethnic differences in CCT: The cornea tends to be thicker in Caucasians, thinner in Africans, and between the two in East Asians.13 14 15 Ethnic differences in prevalence and pattern of glaucoma between East Asian and Caucasian populations have been noted as well.16 17 Therefore, it would be interesting to confirm the heritability findings of CCT in East Asian populations. In the present investigation, we investigated the heritability estimate for CCT in a twin cohort aged from 8 to16 years of age recruited through a population-based twin registry in southern China. All twins were of homogenous Han ethnicity and were healthy without corneal disease and elevated IOP. We used a single imaging system (Pentacam; Oculus, Wetzlar, Germany), which is a noncontact imaging technique, to obtain objective measurements.
| Materials and Methods |
|---|
|
|
|---|
Zygosity of all same-sex twin pairs was determined by 16 multiplex STRs (PowerPlex 16 System; Promega, Madison, WI)19 at the Forensic Medicine Department of Sun Yat-Sen University in 2006. Opposite-sex twin pairs were deemed dizygotic and did not require genotyping.
Examination and Measurement
The imaging system (Pentacam; Oculus) captured the image of the anterior segment by a rotating Scheimpflug camera and a monochromatic slit light source (light-emitting diode at 475 nm). Using the displayed arrows, one experienced examiner (GH) focused and aligned the real-time image of the subjects eye. The 25 images per scan option was selected and automatic-release mode was used. The measurement automatically started only when correct alignment and focus of the eyes were achieved. All measurements were made from 9 to 12 AM and from 2:30 to 5:30 PM. The pachymetric results were automatically calculated by the device and output to a data sheet (Excel; Microsoft, Redmond, WA). The corneal thickness defined as a measurement at the apex point was taken to be the CCT in the present study. Only the measurements with a quality factor (QS) of >95% were considered valid and were included in the analysis. No eye drop was used before the imaging examination. The CCT results on the right and left eyes were similar (correlation coefficient = 0.93; P < 0.0001). Thus, we arbitrarily chose the right eyes of each twin in the data analysis.
Data Analysis and Genetic Modeling
Maximum likelihood correlations were calculated for the five groups of twins (monozygotic male [MZM], dizygotic male [DZM], monozygotic female [MZF], dizygotic female [DZF], and opposite sex dizygotic [OSDZ] twins) and model-fitting analyses were performed with a special software package (Mx; Statistical Modeling, Richmond, VA).20
In classic twin studies, variation of a trait can be decomposed into sources of additive genetic (A), common environmental (C), and unique environmental (E) effects.21 MZ twins are derived from one fertilized egg and share 100% of their genes. DZ twins are derived from two distinct fertilized eggs and share, on average, 50% of their genes. Because the twins in our sample were reared together, both MZ and DZ twins shared 100% of common family environmental effects. Heritability is defined as the proportion of the total variance attributable to genetic variance.
The use of opposite-sex DZ twin pairs in twin studies provides an opportunity to detect sex-specific effects as well as gender difference in genetic and environmental influences on the phenotype under study.21 In the present study, opposite-sex DZ correlation was lower than same-sex DZ correlation, suggesting that sex-specific effects may play a role in CCT variation. Thus, sex-limitation models were fit to the data using the model-fitting analysis software (Mx; Statistical Modeling). Figure 1 depicts the general sex-limitation model we used for the data. In the full general sex-limitation model (model 1), the A, C, and E parameters were assumed to differ between the boys and the girls (Am
Af, Cm
Cf, and Em
Ef). The full general sex-limitation model also assumed the existence of sex-specific genes by allowing the additive genetic correlation for opposite-sex twins (raO) to vary between 0 and 0.5. Variations of the full sex-limitation model were made to identify the best-fitting model. Four steps were taken specifically. First, to determine whether the full general sex-limitation model is acceptable, we compared the fit in that model with that in a saturated model where means and variances of CCT were allowed to differ across zygosity as well as between the first- and the second-born twins. Second, we fixed raO at 0.5, to detect effects of sex-specific genes. Third, we constrained the A, C, and E parameters to be equal across the two sexes, to determine the difference in the magnitude of genetic and environmental influences between the two sexes. Finally, either or both A and C parameters were eliminated from the full model and from the two sets of reduced sex-limitation models to test the significance of its effects.
|
2, with the degree of freedom equal to the difference in degree of freedom between the two models. The selection of the best-fitting model was made using the log-likelihood ratio test (LRT): A significant change in
2 between the full and reduced models suggested that the reduction was not acceptable, whereas a nonsignificant change in
2 indicated that the reduced model was better than the full model. | Results |
|---|
|
|
|---|
|
|
|
2 between the saturated model and the full general sex-limitation model (model 1) was not significant (
2= 21.45, for 16 df; P = 0.162), indicating that the general sex-limitation model was acceptable. Sex and age effects in twin analyses can bias estimates of genetic and environmental influences on the trait under study.21 Therefore, we treated sex and age as covariates in the models, to control their main effects. Fixing the additive genetic correlation for OSDZ twins (raO) at 0.5 did not worsen the model fit (model 2), indicating that the same genes may have been active in the boys and the girls. When we equated A, C, and E across sexes, the resultant change in
2 attained a borderline significance (model 3). Eliminating the C parameter from both sexes in model 2 did not significantly worsen the fit (model 4), whereas removing the A parameter did (model 5). Thus, model 4 was chosen as the best-fitting, most parsimonious model. In this model, additive genetic factors were 87.8% (95% CI: 83.9%–90.7%) for the boys and 91.4% (95% CI: 88.8%–93.3%) for the girls. The remaining unique environmental influences including measurement error were 12.2% (95% CI: 9.3%–16.1%) for the boys, and 8.6% (95% CI: 6.7%–11.2%) for the girls.
|
| Discussion |
|---|
|
|
|---|
Previous studies on the heritability estimates of CCT are scarce. Based on 86 families and 187 sibs, Alsbirk22 suggested familial aggregation of CCT with heritability estimates of 0.6 to 0.7. However, this estimate was based on simple correlations rather than on genetic modeling. Moreover, family studies cannot distinguish between genetic factors and shared environmental factors within families. The heritability was substantially lower in the Alsbirk study than in the present study. This finding is explicable on the ground that Alsbirks study included both adults and children, while our study included only youths. Heterogeneity in age in the Alsbirk sample may have resulted in a low heritability estimate for CCT. Another possible reason is that Alsbirk used a simple and nonautomated optical pachymetry, which may to some extent lower the within-pair correlations of CCT.
To our knowledge, only one study has been conducted to estimate the heritability of CCT in a twin design: In 256 twin pairs (131 MZ and 125 DZ) from Australia and the United Kingdom, Toh et al.12 reported a heritability estimate of 0.95 (95% CI: 0.93–0.96) and near 0 common environmental effects for CCT (intraclass correlation was 0.95 for the MZ and 0.52 for the DZ twins). Similarly, our study yielded heritability estimates of 0.88 (95% CI: 0.84–0.91) in the boys and 0.91 (95% CI: 0.89–0.93) in the girls with a substantially larger sample size (297 pairs of MZ twins and 152 pairs of DZ twins). Our study demonstrates that the heritability of CCT is similar in Chinese and Europeans, suggesting consistency, with the possibility of a central role of genetic factors. It would be interesting to explore further and compare the effects and number of CCT-determining genes among ethnic groups, given that CCT variations have been shown in different ethnic populations.13 14 15
In our analysis, the most parsimonious model was the one with higher CCT heritability in the girls than in the boys. This finding must be treated with caution, as the confidence intervals of heritability in boys and girls widely overlapped, reflecting an insufficient power in the present sample. There was no previous evidence supporting a difference between the sexes in CCT, and it is difficult to provide biological explanations of this finding. Previous studies of sex variations of CCT have been controversial. Patterns of male over female,23 24 female over male,13 and no sex difference12 14 have been reported. Although hormonal change in women is noted as having a possible impact,25 this effect may not be applicable to our young twins (aged, 8–16 years). Future research with a larger sample would help to clarify whether the current result can be replicated. For a complex trait with sex-specific genetic architecture, one broad implication of the present finding is that genetic models incorporating sex effects (for example, a sex-limited linkage model) can increase the ability to detect signals of susceptibility loci during a genome-wide screening.26
There is not a simple explanation of the genetic effect of CCT, given that its biological determinants remain unknown. No specific genes or chromosomal regions have been reported to be linked to CCT. An embryonic study showed that development of the cornea starts as early as 5 weeks in utero.27 Ehlers et al.28 suggested that CCT in children could reach adult levels before 2 years of age. Some rare disorders have been found to have a distinct CCT value in clinical studies. It is noteworthy that some of them are glaucoma-associated developmental syndromes, such as iris hypoplasia,29 aniridia,30 dysgenetic lens,29 pseudoexfoliation syndrome,31 Peters anomaly,32 and Axenfeld-Rieger Syndrome.33
The use of objective measurement for phenotyping is an important issue in conducting twin studies. Twins, especially MZ twins tend to participate together in the study examination, often dressing alike. Examiners who use a subjective instrument (such as optical pachymetry) can be unintentionally affected by the appearance of the twins and may give similar readings for those who look alike. In our study, we were able to collect CCT data objectively based on an imaging system (Pentacam; Oculus) that includes a rotating Scheimpflug camera that yields CCT readings with high repeatability and interoperator reproducibility.34 This instrument provides measurements comparable or even exchangeable with those obtained by ultrasound pachymetry.34 Moreover, the noncontact feature of the system makes it more feasible for making measurements in the children.
Several limitations of this study should be noted. First, previous studies have shown that results from twin studies can be generalized to the singleton population for many complex traits and mortality.35 However, whether findings on CCT from twin studies can be extrapolated to the singleton population is uncertain. To date, no study has been undertaken to compare CCT between twins and singletons in a uniform methodologic protocol. Nevertheless, in comparison findings in the present study (CCT = 551.8 ± 32.8 µm), investigators in a study identified a very similar CCT distribution (553 ± 32.7µm) in 1233 Chinese children in rural China (age: 14.7 ± 0.8 years),36 and a similar level of CCT (546.0 ± 31.8 µm) was also found in Singaporean Chinese children (age, 9–11 years) when slit lamp optical pachymetry was used.23 Second, a heritability estimate is population-specific. Our results obtained in Chinese children may not be applicable to other populations or adults. Third, classic twin studies are based on an assumption of equivalent environment,37 which suggests that both types of twins share broadly the same environment. Although it could be difficult to determine whether this assumption is true, it is unlikely that twin studies on such ocular biometry measurement as CCT may have violated this assumption.
In this study, we have demonstrated that genetic effects play an important role in determining CCT variation in healthy Chinese children, with heritability estimates of 0.88 and 0.91 in the boys and the girls, respectively. We believe that our results will help guide molecular genetic investigations to search genes for CCT and enhance understanding of the etiology of glaucoma.
| Footnotes |
|---|
Submitted for publication February 26, 2008; revised March 27 and May 15, 2008; accepted August 6, 2008.
Disclosure: Y. Zheng, None; J. Ge, None; G. Huang, None; J. Zhang, None; B. Liu, None; Y.-M. Hur, None; M. He, 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 |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. A. Landers, A. W. Hewitt, D. P. Dimasi, J. C. Charlesworth, T. Straga, R. A. D. Mills, R. Savarirayan, D. A. Mackey, K. P. Burdon, and J. E. Craig Heritability of Central Corneal Thickness in Nuclear Families Invest. Ophthalmol. Vis. Sci., September 1, 2009; 50(9): 4087 - 4090. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |