|
|
||||||||
1From the Menzies Research Institute, University of Tasmania, Hobart, Australia; the 2Department of Genetics, Southwest Foundation for Biomedical Research, San Antonio, Texas; the 3Genetics and Bioinformatics Division, Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia; the 4Center for Eye Research Australia, Royal Victorian Eye and Ear Hospital, University of Melbourne, Melbourne, Australia; the 5Department of Ophthalmology, Flinders University, Adelaide, Australia; and the 6Center for Human Genomics and the 7Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, North Carolina.
| Abstract |
|---|
|
|
|---|
METHODS. Genome-wide multipoint variancecomponents linkage analyses of maximum recorded intraocular pressure (IOP) and maximum vertical cup-to-disc ratio were conducted on data from a single, large Australian POAG pedigree that has been found to segregate the myocilin Q368X mutation in some individuals.
RESULTS. Multipoint linkage analysis of maximum recorded IOP produced a peak LOD score of 3.3 (P = 0.00015) near marker D10S537 on 10q22, whereas the maximum cup-to-disc ratio produced a peak LOD score of 2.3 (P = 0.00056) near markers D1S197 to D1S220 on 1p32. Inclusion of the myocilin Q368X mutation as a covariate provided evidence of an interaction between this mutation and the IOP and cup-to-disc ratio loci.
CONCLUSIONS. Significant linkage has been identified for maximum IOP and suggestive linkage for vertical cup-to-disc ratio. Identification of genes contributing to the variance of these traits will enhance understanding of the pathophysiology of POAG as a whole.
The clinical diagnosis of glaucoma is based on a combination of several main features, including specific changes to the appearance of the optic nerve head constituting glaucomatous optic neuropathy, characteristic visual field loss with a slow and often asymptomatic progression, and, in most cases, increased intraocular pressure (IOP). The complexity of the phenotypic definition of POAG6 7 has contributed to the difficulties in identifying genes involved in this disease.
Two genes have been identified for adult-onset POAG8 9 ; however, mutations in these genes account for only a fraction of POAG cases. Myocilin (MYOC) on 1q24.3 has been shown to account for approximately 3% of adult-onset POAG.9 10 11 12 Optineurin (OPTN) on 10p15-p14 was shown by Rezaie et al.8 to be involved in up to 17% of low-tension glaucoma pedigrees, although a recent study indicated a low prevalence of OPTN mutations (<0.1%) in unselected cases of both POAG and NTG.13 The WD40-repeat 36 gene (WDR36) was recently identified at the GLC1G locus on 5q22.1,14 although the impact of this gene on the wider POAG population is yet to be determined. In addition to these three genes, four other loci for POAG have been mapped: GLC1B,15 GLC1C,16 GLC1D,17 and GLC1F.18 The results of three genome-wide scans suggest that several additional chromosomal regions may also be involved in susceptibility to POAG19 20 21 ; however, only GLC1B22 and GLC1C23 have been replicated in published studies.
Given the limited success so far in identifying genes conferring susceptibility to glaucoma, one approach that may have greater success is quantitative trait linkage analysis of precursors of glaucoma, such as raised IOP and increased cupping of the optic nerve. Such traits may have simpler genetic architecture than diagnosis of glaucoma, making it easier to map causative loci.24 25 26 Furthermore, quantitative trait linkage analysis is inherently more powerful than dichotomous trait linkage analysis27 28 and is particularly powerful in large families.29 30
The heritabilities of IOP and vertical cup-to-disc ratio were estimated to be 0.36 and 0.48, respectively, in the Beaver Dam Eye Study, providing evidence for genetic determinants for these components.31 Commingling analysis of IOP and glaucoma by Viswanathan et al.32 suggested the existence of a major gene accounting for 18% of the variance of IOP in the Blue Mountains Eye Study population.32 Duggal et al.33 recently conducted a complex segregation and linkage analysis of IOP, identifying two potential regions of linkage on chromosomes 6 and 13.
The purpose of the present investigation was to identify regions of linkage that contribute to maximum IOP and maximum cup-to-disc ratio, using measures from an extended Australian pedigree. This family already has been analyzed for linkage to glaucoma.34 The discovery of genes contributing to the variance of maximum IOP and maximum cup-to-disc ratio is expected to provide significant insights into the pathophysiology of glaucoma.
| Methods |
|---|
|
|
|---|
GTas02, the family of interest in this investigation, is one of the largest pedigrees identified as part of the GIST. It consists of more than 1350 members, and ancestry can be traced back to a founder couple six generations ago. The "core" pedigree containing the POAG cases consists of 246 individuals when deceased linking members are included. One hundred thirty-nine family members consented to clinical examination and blood collection for genotyping and mutation analysis. The structure, clinical diagnosis, and MYOC mutation status of this family have been published.36
Clinical Examination
Clinical examination and diagnosis of patients involved in the GIST are documented elsewhere.37 In brief, POAG was clinically defined as an optic neuropathy that had at least two of the following features: (1) optic nerve head excavation with thinning of the neuroretinal rim, often with Drance-type nerve fiber layer hemorrhages, notching, pitting, significant focal loss or general loss of the retinal nerve fiber layer (generally measured by an enlarged vertical cup-to-disc ratio
0.7); (2) elevated IOP above a population-based normal range or above the average of the unaffected individuals within a pedigree (generally IOP >21 mm Hg or two standard deviations from the population mean); and (3) visual field defects consistent with the disc changes and with common descriptions of glaucomatous field loss.37 Glaucoma cases secondary to trauma or anterior segment dysgenesis were excluded. Of the 139 individuals available for examination, 24 had a diagnosis of POAG. The details of the clinical diagnoses for this family have been extensively reported.36
IOP was measured with a calibrated Goldmann applanation tonometer. Multiple IOP measures were available for each individual; hence, maximum IOP was selected, to reduce any bias introduced by using postmedication pressure values. No corrections were made for corneal thickness, because this information was not routinely collected at the time of patient ascertainment. Optic disc appearance was classified by two clinicians at the time of examination, with a slit lamp biomicroscope after pupil dilation. In addition, optic disc stereo photographs (Nidek, Gamagori, Japan) were obtained for future reference in all cases. When there was a discrepancy between the two examiners, the stereo disc photographs were independently assessed by a glaucoma specialist. The highest vertical cup-to-disc ratio in either eye at any clinical examination was used as the trait measure.
Quantitative measures collected as part of the clinical examination of patients with POAG and their relatives included maximum recorded intraocular pressure (IOP) without medication and maximum vertical cup-to-disc ratio from the highest-scoring eye.36
Myocilin mutation detection in this family has been reported.36 38 Of the 139 individuals from GTas02 screened for mutations in the MYOC gene, 19 were known to carry the Q368X mutation.
Genotyping
A 10-cM genome-wide scan was conducted with 401 microsatellite markers from fluorescence marker sets (vers. 1 and 2; Applied Biosystems, Inc. [ABI], Foster City, CA), run on sequencers (model 377; ABI) and analyzed (GeneScan and GenoTyper software; ABI).34 Inconsistencies in Mendelian inheritance of genotypes were detected with Pedcheck.39 Marker allele frequencies were estimated from 72 elderly glaucoma-free control individuals drawn from the same population.40
Multipoint identity by descent (IBD) files were created with the Markov chain Monte Carlo (MCMC)based program Loki (ver. 2.4.7)41 42 from within SOLAR (Sequential Oligogenic Linkage Analysis Routines, ver. 2.1.1; http//:www.sfbr.org/ provided in the public domain by the Southwest Foundation for Biomedical Research, San Antonio, TX).25
Ascertainment Correction
The population ascertainment correction data set was taken from the Melbourne Visual Impairment Project (MVIP) and included 3905 individuals from the general population with data on age, sex, maximum recorded IOP, and maximum cup-to-disc ratio.43 The MVIP clinical data were collected by the same methods as were used in the GIST. A comparison between the MVIP and family GTas02 data is shown in Table 1 . Mean age and trait variables in family GTas02 and the ascertainment correction sample were compared by two-tailed t-tests. The ascertainment correction dataset was not intended to be age matched; however, comparison of age and trait between the two groups revealed the relative position of family GTas02 within the distribution of the general population.
|
Variance Components Analysis
Heritability for each trait was estimated with genetic variance component modeling as implemented in SOLAR (ver. 2.1.1).25 The covariates were selected among age, sex, agesex interaction; however, agesex interaction was not significant in any dataset and was thus removed. Variance component linkage analysis was performed to detect and localize quantitative trait loci (QTLs) influencing variation in maximum IOP and maximum cup-to-disc ratio, by using SOLAR.25
The variance-component linkage method is based on specifying the expected genetic covariances between arbitrary relatives as a function of IBD relationships at a given marker locus.25 The method involves partitioning the total trait phenotypic variance (
2P) into components attributable to covariate effects, effects of a specific QTL (QTL (
2q), and residual additive genetic effects (
2h2r). The test for linkage compares the likelihood of this model with the likelihood of a null model (no linkage), where the QTL effect size
q2 is fixed to be zero. The difference between the two log10 likelihoods produces a LOD score that can be interpreted in a fashion similar to that of the classic LOD scores of parametric linkage analysis.44 The variance component quantitative genetic approach enables penetrance model-free multipoint linkage analysis of complex quantitative traits in pedigrees of arbitrary size and complexity25 28 45 46 47 and has been used successfully to localize QTLs that influence many important disease-related traits, including risk of alcoholism,47 serum leptin levels,48 49 and resting heart rate.49
Expected LOD scores and empiric locus-specific probabilities were determined with the "lodadj" command of SOLAR.44 50 A total of 100,000 replicates were simulated to build up the distribution of LOD scores expected under the null hypothesis of no linkage, using a fully informative, unlinked marker. The observed LOD scores were then regressed on those expected for a multivariate normal trait, with the inverse slope of the regression line providing the LOD correction constant. Final LOD scores were multiplied by the correction constant only if the constant was <1.
To test the impact of the MYOC Q368X mutation on the linkage results, Q368X status was included as a covariate in a second round of analyses for each trait. The regression coefficient for Q368X mutation status (ßMYOC), used in the analysis of each trait, was estimated from the family GTas02 trait and mutation data and not from the population ascertainment correction dataset. As in previous analyses, age and sex were included in all models, and adjusted LOD scores and empiric probabilities were calculated as described earlier. The fitted coefficient for MYOC mutation status (ßMYOC) in the analysis of IOP was 5.065, whereas ßMYOC for the analysis of maximum cup-to-disc ratio was 0.192.
To control for the overall false-positive rate in our linkage screens, we converted the nominal probabilities associated with our peak LOD scores to genome-wide probabilities, using an approach developed by Feingold et al.51 This method takes into account the mean recombination rate in our study population and the marker density of the linkage map used in our genome scan.
| Results |
|---|
|
|
|---|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
The putative trait locus for maximum recorded IOP produced a significant peak LOD of 3.3 on 10q22. The genome-wide significance level of this result (P = 0.0165) strongly suggests that this region contains a gene that contributes to the variance of IOP. We did not see any overlap with the IOP linkage regions identified by Duggal et al.33 Although our region of interest for IOP on 10q22 has not been reported previously for IOP, linkage to this region has also been found for systemic hypertension in a Japanese population.52 The association between systemic blood-pressure (systolic or diastolic) and IOP has been well documented.53 54 55 It is possible that systemic hypertension and IOP share a common QTL in this region on the long arm of chromosome 10. The region contains BMPR1A (MIM: 601299), a bone morphogenic protein receptorinteresting because Bmp4 has been implicated in increased IOP in mice.56 The peak region also contains RGR (MIM: 600342), an opsin-related gene associated with retinitis pigmentosa.57 The OPTN gene is located outside our 10p13 linkage peak, approximately 60 cM upstream.
Applanation tonometry measurements of IOP are known to be influenced by central corneal thickness (CCT).58 59 CCT has a positive and apparently linear correlation with IOP60 and is strongly genetically determined (Toh TY, et al. IOVS 2005;46:ARVO E-Abstract 1093). There has been a report of thick corneas segregating in a family with apparent ocular hypertension,61 which leads to the suggestion that perhaps the variation in IOP in family GTas02 compared with the population is due to CCT. However, there are many individuals in family GTas02 with moderate to advanced visual field loss, which would support the elevated IOPs being genuine rather than artifactual in nature. It is possible that measurement of CCT and subsequent adjustment of Goldmann applanation IOP readings for GTas02 may have the effect of strengthening the linkage at the IOP locus on 10q22 by providing a more accurate approximation of true IOP (as opposed to measured IOP). This approach is worthy of future investigation.
This study identified a putative trait locus for maximum cup-to-disc ratio on 1p32the first reported locus for this traitwith a peak LOD score of 2.3. This corresponds to suggestive linkage, because there is a 20% chance of a LOD score of this size or greater occurring by chance in a genome-wide scan. The region covered by the linkage peak includes the gene POMGnT1 (MIM: 606822), mutant forms of which are responsible for muscle-eye-brain disease (MEB; MIM: 253280), a congenital muscular dystrophybased disorder with many additional features, including early-onset glaucoma, optic nerve atrophy, severe congenital myopia and retinal hypoplasia. MEB is inherited as a loss of function of POMGnT1.62 The peak region also contains FOXE3 (MIM: 601094), known to be associated with anterior segment ocular dysgenesis63 and is approximately 25 cM distal to the primary congenital glaucoma locus GLC3B (MIM: 600975), on 1p36.64
The MYOC Q368X mutation is the most common mutation identified in patients with POAG to date, found to account for approximately 1.6% of POAG.11 In comparison to other MYOC mutations associated with the juvenile-onset form of POAG (JOAG), this mutation generally gives rise to a milder phenotype with late age of onset and increased IOP.65 Myocilin mutation status and related phenotype modification effects within family GTas02 have been reported.36 38 The MYOC locus was not apparent in either of the linkage analyses for the two quantitative traits investigated in this study. Because only 9 of the 19 Q368X mutation carriers in this family were clinically diagnosed with POAG, this mutation is unable to account for the presence of all cases of POAG in the entire family.38 Despite the absence of a linkage signal at MYOC, the Q368X mutation appeared to account for nearly half the genetic variance for the quantitative trait based on maximum cup-to-disc ratio, and nearly 20% of the genetic variance for maximum IOP. The differences between carrier and noncarrier-trait values (shown in Table 2 ) were significant, which may indicate epistasis or some other form of interaction between MYOC and the putative trait loci for IOP and cup-to-disc ratio. The reduction in peak LOD scores for both traits after the inclusion of the Q368X mutation as a covariate provides further evidence for the presence of an interaction. However, the sample size of family GTas02 does not provide sufficient power to determine the exact mechanism behind this interaction. Subsequent identification of the genes at these loci is anticipated to aid investigation into the nature of this interaction.
Baird et al.34 recently identified linkage of clinical POAG diagnosis to 3p21-p22 in this family, using MCMC-based linkage analysis. We did not find any evidence of overlap with this region using the traits maximum IOP and cup-to-disc ratio. There are likely to be many genes contributing to the complex POAG phenotype, and different analytical approaches will have variable power to detect certain loci. It is also possible that the gene at the 3p locus may contribute to a different POAG trait, such as progression from elevated IOP to optic nerve damage and subsequent clinical signs of visual field loss. Although we are unaware of any systemic biases in pedigree or family member ascertainment, genotyping methodology, or trait measurements, if any such biases exist they would impact both linkage studies and any subsequent analyses of this dataset. Finally, as with any genome scan, one must also be mindful that some of these linkage results could represent chance events.
The loci identified in this study are believed to contribute to the variance of IOP and cup-to-disc ratio in the general population. However, since more extreme values are present in a POAG pedigree, when combined with population-based ascertainment correction, the approach used in this investigation is expected to provide greater power to map genes for these traits. Maximum IOP and maximum vertical cup-to-disc ratio may not be ideal measures, given the influence of pressure spikes and diurnal variation on IOP, for example, however other traits such as mean IOP are likely to be equally problematic considering the probable inclusion of postmedication levels. In addition, any measurement error in IOP and cup-to-disc ratio would also be present in the population data used for ascertainment correction.
In this investigation, we were able to identify regions of linkage contributing to maximum recorded IOP and maximum cup-to-disc ratio, quantitative clinical contributors to POAG diagnosis that are collected as a routine part of clinical practice. The discovery of the genes involved in these components of POAG at these loci is anticipated to provide significant insights into glaucoma pathophysiology as a whole.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication March 11, 2005; revised June 10, 2005; accepted August 24, 2005.
Disclosure: J.C. Charlesworth, None; T.D. Dyer, None; J.M. Stankovich, None; J. Blangero, None; D.A. Mackey, None; J.E. Craig, None; C.M. Green, None; S.J. Foote, None; P.N. Baird, None; M.M. Sale, 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: Michèle M. Sale, Center for Human Genomics, Wake Forest University School of Medicine, Medical Center Boulevard, Winston-Salem, NC 27157; msale{at}wfubmc.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
P Founti, F Topouzis, L van Koolwijk, C E Traverso, N Pfeiffer, and A C Viswanathan Biobanks and the importance of detailed phenotyping: a case study--the European Glaucoma Society GlaucoGENE project Br. J. Ophthalmol., May 1, 2009; 93(5): 577 - 581. [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] |
||||
![]() |
Y. S. Aulchenko, D.-J. de Koning, and C. Haley Genomewide Rapid Association Using Mixed Model and Regression: A Fast and Simple Method For Genomewide Pedigree-Based Quantitative Trait Loci Association Analysis Genetics, September 1, 2007; 177(1): 577 - 585. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. W. Hewitt, S. L. Bennett, J. H. Fingert, R. L. Cooper, E. M. Stone, J. E. Craig, and D. A. Mackey The Optic Nerve Head in Myocilin Glaucoma Invest. Ophthalmol. Vis. Sci., January 1, 2007; 48(1): 238 - 243. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. McLaren, D. M. Reed, D. C. Musch, C. A. Downs, M. E. Higashi, C. Santiago, P. A. Radenbaugh, R. R. Allingham, J. E. Richards, and S. E. Moroi Evaluation of the beta2-Adrenergic Receptor Gene as a Candidate Glaucoma Gene in 2 Ancestral Populations Arch Ophthalmol, January 1, 2007; 125(1): 105 - 111. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Y. Wang, B. J. Fan, J. K. H. Chua, P. O. S. Tam, C. K. S. Leung, D. S. C. Lam, and C. P. Pang A Genome-wide Scan Maps a Novel Juvenile-Onset Primary Open-Angle Glaucoma Locus to 15q Invest. Ophthalmol. Vis. Sci., December 1, 2006; 47(12): 5315 - 5321. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. A. Radenbaugh, A. Goyal, N. C. McLaren, D. M. Reed, D. C. Musch, J. E. Richards, and S. E. Moroi Concordance of Aqueous Humor Flow in the Morning and at Night in Normal Humans Invest. Ophthalmol. Vis. Sci., November 1, 2006; 47(11): 4860 - 4864. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |