|
|
||||||||
1From the Glaucoma Service, The Rotterdam Eye Hospital, Rotterdam, The Netherlands; the 4Department of Ophthalmology and the 2Departments of Epidemiology and Biostatistics and 3Clinical Genetics, Genetic Epidemiology Unit, Erasmus Medical Center, Rotterdam, The Netherlands; and the 5Department of Molecular and Clinical Ophthalmogenetics, The Netherlands Institute for Neuroscience, Amsterdam, The Netherlands.
| Abstract |
|---|
|
|
|---|
METHODS. A total of 2620 subjects (mean age, 48 years; range 18–86) from extended pedigrees living in a small town in The Netherlands underwent an extensive ophthalmic examination. Their IOP was measured by Goldmann applanation tonometry, their RNFL thickness by scanning laser polarimetry (GDx VCC), and their optic disc parameters by confocal scanning laser ophthalmoscopy (HRT II). Risk associations were explored by linear regression analyses and heritability estimates by variance component methods.
RESULTS. Inbreeding was present in 2042 (81%) participants, and was significantly associated with a higher IOP (P < 0.001). The heritability estimate for IOP was 0.35 (95% confidence interval [CI], 0.27–0.43); for RNFL thickness, 0.48 (95% CI, 0.35–0.60); and for neuroretinal rim area, 0.39 (95% CI, 0.20–0.58). Nongenetic factors accounted for only a small proportion (
0.13) of the variance in all three traits.
CONCLUSIONS. Early, continuous markers of POAG are strongly determined by additive genetic effects. The results support a quantitative trait linkage strategy to discover new genes for POAG.
Gene-finding in POAG has been hampered by etiological and clinical heterogeneity, partly due to nonuniformity of diagnostic criteria. Its insidious onset and slow natural course impede a definite diagnosis at an early stage, whereas including only late-stage, outright POAG will greatly limit genetic studies because relatively few people would be available to participate. These difficulties fuel the question of whether studying heritable, continuous markers of POAG may improve the chance of success in gene finding.
Quantitative markers of POAG, apart from IOP, are retinal nerve fiber layer (RNFL) thickness and optic disc rim area, both indicators of the number of existing retinal ganglion cells, or indirectly, of any of their loss, typical of glaucoma.20 21 22 Previous studies, mostly performed in twins or nuclear families, reported heritability estimates for IOP ranging from 0.29 to 0.50 and for cup-to-disc ratio ranging from 0.48 to 0.80.23 24 25 26 The ranges of these heritability estimates are relatively large, and studies based on more extended pedigrees may provide more precise figures.
We performed a family-based cohort study in a genetically isolated population in The Netherlands and thereby had the opportunity to study large extended pedigrees. As a first investigation into the genetic etiology of early markers of POAG, we explored the heritability estimates for IOP, RNFL thickness, and optic disc rim area. We used the imaging techniques scanning laser polarimetry (SLP) and confocal scanning laser ophthalmoscopy (CSLO) to obtain objective measurements. This study also allowed us to demonstrate the effect of inbreeding and nongenetic factors on these POAG markers.
| Methods |
|---|
|
|
|---|
Eligibility for participation in the study was determined by genealogical background, not by any phenotypes of interest. Twenty-two families were selected who had at least six children baptized in the community church between 1880 and 1900. All living descendants of these families aged 18 years and older, as well as their spouses, were invited to attend a series of clinical examinations. A total of 2620 subjects underwent ophthalmic examination.
Based on genotyping half of the sample (n = 1430, 437 nuclear families with one or both tested parents and one or more tested offspring) with 5964 single nucleotide polymorphisms (SNPs), we identified only one nuclear family in which the father was not compatible with any of the three offspring (whereas the siblings were compatible with each other). Exclusion of this family did not alter the heritability estimates. As the percentage of nonpaternities is low and is known for only half the sample, the presence of nonpaternities was ignored in the data presented.
Clinical Examination
All data were collected between June 2002 and February 2005. Nonophthalmic examinations included anthropometric measurements, cardiovascular and endocrine assessments, neuropsychological tests, fasting blood samples, and interviews regarding medical history, medication, and putative risk factors.
The ophthalmic examination comprised the assessment of best-corrected visual acuity, refraction, and IOP. Keratometry was determined by an automatic refractometer, and the eyes axial lengths by an intraocular lens calculator (IOL Master; Carl Zeiss Meditec, Inc., Dublin, CA). Scanning laser polarimetry (SLP) was subsequently performed with the commercially available GDx VCC (Carl Zeiss Meditec, Inc.). In mydriasis, participants underwent fundus photography centered on the optic disc (20°, TRC-50XT retinal camera; Topcon Medical Systems, Inc., Paramus, NJ) and confocal scanning laser ophthalmoscopy (CSLO) measurements (Heidelberg Retina Tomograph II; [HRT II] Heidelberg Engineering GmbH, Dossenheim, Germany).
Intraocular Pressure
One of five trained examiners performed bilateral IOP measurements with Goldmann applanation tonometry. A drop of fluorescein sodium was instilled in each eye. The tonometer was set at 10, and the prism was carefully applied to the corneal surface of the right eye. Without looking at the scale, the examiner rotated the dial until the inner margins of the two semicircles touched each other. The examiner then moved the slit lamp away from the eye and read the IOP. The tonometer was set at 10, and the measurement was repeated. If the two measurements differed, a third measurement was performed, and the median value was recorded. The procedure was repeated for the left eye.
RNFL Thickness
After any refractive error was entered into the GDx VCC software, the birefringence of the anterior segment of each participant was assessed by means of the method described by Zhou and Weinreb.30 Subsequently, each eye was scanned with adjusted anterior segment compensation to estimate peripapillary RNFL thickness as described by Reus and Lemij.31 The cutoff for the quality of the image was a GDx VCC scan quality score of 8 or above. Images with lower scores were excluded.
After the boundaries of the optic disc were manually marked, the GDx VCC software positioned two circles centered on the disc: The first had a diameter of
2.5 mm (54 pixels), the second a diameter of
3.3 mm (70 pixels). The parameters of RNFL thickness that we used in this study were based on the RNFL thickness measurements between the circles: TSNIT average (temporal-superior-nasal-inferior-temporal), superior average (25°–144°), inferior average (215°–334°), and the worst hemifield average (the lowest average value of the RNFL thickness of the hemifields of both eyes).
Neuroretinal Rim
Details of the CSLO technique have been described elsewhere.32 Briefly, imaging was performed after entering the participants keratometry data into the software and after adjusting the settings in accordance with the refractive error. Only images with a standard deviation of height measurements below 50 µm were accepted. The optic disc margin was manually marked at the inner edge of Elschnigs ring by one observer (LMEvK). The HRT II software then calculated stereometric parameters of the optic disc and neuroretinal rim. The parameters that we used were disc area, rim area, rim area superotemporally (45°–90°), rim area inferotemporally (270°–315°), rim-to-disc area ratio, vertical cup-to-disc ratio, and cup shape measure, an index of depth variation and steepness of the cup walls.
Statistical Analysis
The inbreeding coefficient, which represents the probability that the two alleles at a given locus are identical by descent (i.e., derived from the same ancestral chromosome), was calculated, based on all available genealogical information, by means of PEDIG software.33 This coefficient was analyzed in quartiles, since its distribution was skewed toward zero.
Other putative covariates of glaucoma pathogenesis that were studied included age, sex, height, body mass index, systolic blood pressure, pulse rate, fasting blood glucose level, blood cholesterol level, IOP, time of IOP measurement, axial length of the eye, spherical equivalent of refractive error, and mean corneal curvature. The glaucoma markers that were studied were IOP, TSNIT average, superior average, inferior average, worst hemifield average, disc area, rim area, superotemporal rim area, inferotemporal rim area, rim-to-disc area ratio, vertical cup-to-disc ratio, and cup shape measure. These markers were based on the eye most representative of glaucoma: the eye with the lower RNFL thickness parameters, rim area parameters, and rim-to-disc area ratio and the eye with the higher IOP, vertical cup-to-disc ratio, and cup shape measurement. We calculated the mean of both eyes for the analyses of disc area. If a measurement could be obtained in only one eye, the parameters of this eye were included in the analyses.
Associations were explored by univariate and multivariate linear regression (SPSS version 11.0 for Windows; SPSS, Chicago, IL). All determinants below the 0.10 significance level in the multivariate analyses were retained in the final model for heritability estimation. The distribution of the multivariate regression residuals in the final model was tested for normality with the nonparametric, one-sample Kolmogorov-Smirnov test. To reduce kurtosis of the distribution in parameters describing RNFL thickness or rim area, we excluded the upper and lower 0.5 percentile values of these traits. We further transformed traits that were skewed using natural logarithm (disc area, rim area), or exponential function (vertical cup-to-disc ratio).
We estimated the heritability by means of a variance component maximum likelihood analysis, as implemented in Solar software (ver. 2.1.2).34 35 A variance component analysis separates the observed phenotypic variance into components that are attributable to different causes. Heritability describes the relative importance of the component that is attributable to heredity. This component is called the additive genetic variance and represents the cumulative effects of alleles. Heritability can be estimated from the resemblance between family members. We first examined the proportion of the phenotypic variance associated with the covariates. Subsequently, we estimated the proportion of the remaining phenotypic variance explained by additive genetic effects. Finally, the heritability of each glaucoma marker was calculated as the proportion of the total phenotypic variance explained by additive genetic effects. In addition, we investigated the genetic correlation between rim area and disc area with a linear bivariate analysis.
| Results |
|---|
|
|
|---|
|
Intraocular Pressure
Tonometry was successfully performed in 2457 (93.8%) subjects. Twenty-three subjects (0.9%) received IOP-lowering therapy or had a history of these medications and were excluded from the IOP analyses.
IOP values of the right and the left eye were significantly correlated (Pearson correlation coefficient 0.855; P < 0.001). IOP was significantly higher in the right eye than in the left eye (paired samples t-test; mean difference 0.10 mm Hg; 95% CI, 0.03–0.16).
Data have been presented for the eye with the higher IOP. The distribution of the IOP in the study population is given in Table 2 , and the results of the linear regression and variance component analyses are in Table 3 . Age, inbreeding, and fasting glucose were significantly associated with IOP, and corneal curvature was inversely related. These covariates accounted for 0.05 of the total phenotypic variance of IOP. Of the remaining variance, the proportion explained by additive genetic effects was 0.37 (95% CI, 0.29–0.45). The heritability estimate, calculated as the contribution of genetic factors to the total phenotypic variance, was 0.35 (95% CI, 0.27–0.43).
|
|
The distributions of the RNFL thickness parameters are presented in Table 2 . The nerve fiber indicator (NFI) was
40 (i.e., suggestive of glaucoma) in 20 (1.8%) of 1085 subjects younger than 55 years, and 41 (9.2%) of 446 subjects 55 years of age or older.
The results of the linear regression and variance component analyses are reported in Table 4 . Inbreeding and IOP were not significantly related to any RNFL thickness parameter. Age and systolic blood pressure were inversely associated with all parameters, although the relation between systolic blood pressure and inferior average did not reach statistical significance (P = 0.144). Axial length, spherical equivalent, and corneal curvature, all of which are covariates of refractive error, were significantly associated with most RNFL thickness parameters.
|
Neuroretinal Rim and Optic Disc Morphology
Data of the first 750 participants were included in this analysis. Of these, eight (1.1%) subjects were excluded because of poor quality of the topographic images.
The distributions of all optic disc parameters that were measured have been presented in Table 2 . The Moorfields regression classification was outside normal limits in 13 (2.8%) of 458 subjects less than 55 years of age, and 29 (10.2%) of 284 subjects more than 55 years of age.
Table 5 shows the results of the linear regression and variance component analyses. Inbreeding and IOP were not associated with any optic disc parameter. Age showed a statistically significant association with disc area, rim-to-disc area ratio, and vertical cup-to-disc ratio. Determinants of refractive error (axial length, spherical equivalent, and corneal curvature) were significantly related to most parameters.
|
We analyzed coaggregation of rim and disc area and found a low genetic correlation (rgenetic = 0.16; SE = 0.13, P = 0.26).
| Discussion |
|---|
|
|
|---|
The design of our study had several limitations. First, genetically isolated populations may exhibit genetic drift, and their genetic composition may therefore deviate from the general population. We performed simulation studies in the ERF population, which showed that the effects of genetic drift on the frequency of common alleles were negligible.28 Thus, we believe it is valid to generalize our results to an outbred population. A second limitation is that we studied a relatively young and healthy cohort, which may reduce the clinical relevance of our findings. This problem appeared to be small, since the proportion of subjects aged 60+ years was considerable (20%), and the range of outcome variables was representative of the clinical spectrum. Third, measurements of RNFL thickness and neuroretinal rim area were for logistic reasons not performed on the total study population. We do not think that this affected the outcome of our results, because the subsets were chosen randomly and had sufficient statistical power. Fourth, we did not have the opportunity to study central corneal thickness (CCT), a potential confounder of IOP measurements as well as an important risk factor for POAG.36 37 CCT has previously been reported to account for between 1% and 6% of the total variance in IOP measured with Goldmann applanation tonometry.38 39 40 41 The heritability of CCT has been estimated to be 0.95.42 Because we did not include CCT into the variance component analysis of IOP, CCT-determining genes may have been incorrectly considered to be IOP-determining genes, thus inflating our heritability estimate of IOP.
RNFL thickness and optic disc morphology were not associated with IOP in these relatively young individuals not suspected of having glaucoma. An elevated IOP has been recognized as an important risk factor for glaucoma.37 43 IOP levels have been shown to correlate with optic disc characteristics in eyes of subjects with ocular hypertension.44 This may suggest that IOP affects the RNFL and optic disc pathophysiologically rather than physiologically. The results of Chang et al.25 support this hypothesis. Other studies, however, reported a significant effect of IOP on optic disc morphology in normal subjects.45 46
How do our findings relate to other studies of POAG heritability? Our estimate of 0.35 for IOP is remarkably similar to that in the two most recent studies of the heritability of IOP. Klein et al.24 estimated a heritability of 0.35 from the parent–child correlation in the Beaver Dam Eye Study. Chang et al.25 reported a heritability estimate of 0.36 in a population-based cohort of Caucasian sibships 65 to 84 years of age. The heritability estimate of RNFL thickness in our study (0.48) was much lower than the estimate reported by Hougaard et al.47 (0.78–0.82), who studied monozygotic and dizygotic twins. Heritability may be population specific. Even populations with similar genetic backgrounds may show different heritability estimates due to different environmental variances and different study designs. Heritability studies based on twins assume that the environmental correlations among monozygotic and dizygotic twins are equal. However, if this "equal environment assumption" does not hold, the heritability estimated from twin data alone may be higher than that derived from extended families, as in our study.48 Different measurement methods (OCT in Hougaard et al.47 versus GDx VCC in our study) and different adjustments for covariates may also contribute to the heritability differences. There are no former studies that assessed the heritability of optic disc rim area, but three studies estimated the heritability of cup-to-disc ratio as a proxy. Our estimate for this parameter of 0.64 was in their reported range of 0.48–0.80.24 25 26
We were able to study optic disc morphology objectively with the HRT II, a technique that had not been used before in heritability studies. We found heritabilities of 0.52 for disc area, 0.39 for total rim area, and 0.47 for inferotemporal rim area. We found a markedly higher estimate (0.79) for the rim-to-disc ratio, indicating a larger genetic contribution for the combination than for any of the parameters separately. Speculating that different sets of genes may determine rim and disc area, we performed a bivariate analysis of these parameters and found no evidence of any genetic correlation.
As a general rule, inbreeding increases the probability that the gene profile comprises two identical alleles. In our study, inbreeding was significantly associated with a higher IOP, suggesting the presence of causative recessive alleles in the genetic background of this trait. Evidence from previous studies supports this notion. Recessive mutations in the CYP1B1 gene not only link to congenital glaucoma and anterior segment dysgenesis, but also play a role in high-pressure POAG with juvenile or adult onset.8 9 Inbreeding was not associated with RNFL thickness or neuroretinal rim area. Therefore, our study does not support a recessive inheritance of these POAG markers.
In conclusion, we have demonstrated that IOP, RNFL thickness, and neuroretinal rim area are continuous POAG markers that are strongly determined by genetic effects. Genome-wide association methods have been successfully applied to map genes for other complex disorders. A quantitative trait analysis greatly enhances the statistical power of this technology. The high heritabilities that we found in the current study encourage us to use this approach for identifying new POAG genes.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication December 21, 2006; revised March 18, 2007; accepted May 16, 2007.
Disclosure: L.M.E. van Koolwijk, None; D.D.G. Despriet, None; C.M. van Duijn, None; L.M. Pardo Cortes, None; J.R. Vingerling, None; Y.S. Aulchenko, None; B.A. Oostra, None; C.C.W. Klaver, None; H.G. Lemij, Carl Zeiss Meditec, Inc. (C, F)
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: Leonieke M. E. van Koolwijk, Glaucoma Service, The Rotterdam Eye Hospital, P.O. Box 70030, 3000 LM Rotterdam, The Netherlands; koolwijk{at}oogziekenhuis.nl.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. K. Lee, S. J. Woo, J.-i. Kim, S.-i. Cho, H. Kim, J. Sung, J.-S. Seo, and D. M. Kim Replication of a Glaucoma Candidate Gene on 5q22.1 for Intraocular Pressure in Mongolian Populations: The GENDISCAN Project Invest. Ophthalmol. Vis. Sci., March 1, 2010; 51(3): 1335 - 1340. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Carbonaro, T. Andrew, D. A. Mackey, T. L. Young, T. D. Spector, and C. J. Hammond Repeated Measures of Intraocular Pressure Result in Higher Heritability and Greater Power in Genetic Linkage Studies Invest. Ophthalmol. Vis. Sci., November 1, 2009; 50(11): 5115 - 5119. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. M. E. van Koolwijk, P. R. Healey, R. A. Hitchings, P. Mitchell, P. C. Sham, P. McGuffin, and A. C. Viswanathan Major Genetic Effects in Glaucoma: Commingling Analysis of Optic Disc Parameters in an Older Australian Population Invest. Ophthalmol. Vis. Sci., November 1, 2009; 50(11): 5275 - 5280. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Zheng, F. Xiang, W. Huang, G. Huang, Q. Yin, and M. He Distribution and Heritability of Intraocular Pressure in Chinese Children: The Guangzhou Twin Eye Study Invest. Ophthalmol. Vis. Sci., May 1, 2009; 50(5): 2040 - 2043. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Zheng, J. Ge, G. Huang, J. Zhang, B. Liu, Y.-M. Hur, and M. He Heritability of Central Corneal Thickness in Chinese: The Guangzhou Twin Eye Study Invest. Ophthalmol. Vis. Sci., October 1, 2008; 49(10): 4303 - 4307. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Amerasinghe, T. Y. Wong, W.-L. Wong, P. Mitchell, S. Y. Shen, S.-C. Loon, S.-M. Saw, P. J. Foster, T. Aung, and for the SiMES Study Group Determinants of the Optic Cup to Disc Ratio in an Asian Population: The Singapore Malay Eye Study (SiMES) Arch Ophthalmol, August 1, 2008; 126(8): 1101 - 1108. [Abstract] [Full Text] [PDF] |
||||
![]() |
F Carbonaro, T Andrew, D A Mackey, T D Spector, and C J Hammond Heritability of intraocular pressure: a classical twin study Br J Ophthalmol, August 1, 2008; 92(8): 1125 - 1128. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |