Investigative Ophthalmology & Visual Science Cover Image for Volume 46, Issue 10
October 2005
Volume 46, Issue 10
Free
Glaucoma  |   October 2005
Linkage to 10q22 for Maximum Intraocular Pressure and 1p32 for Maximum Cup-to-Disc Ratio in an Extended Primary Open-Angle Glaucoma Pedigree
Author Affiliations
  • Jac C. Charlesworth
    From the Menzies Research Institute, University of Tasmania, Hobart, Australia; the
  • Thomas D. Dyer
    Department of Genetics, Southwest Foundation for Biomedical Research, San Antonio, Texas; the
  • Jim M. Stankovich
    From the Menzies Research Institute, University of Tasmania, Hobart, Australia; the
    Genetics and Bioinformatics Division, Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia; the
  • John Blangero
    Department of Genetics, Southwest Foundation for Biomedical Research, San Antonio, Texas; the
  • David A. Mackey
    Center for Eye Research Australia, Royal Victorian Eye and Ear Hospital, University of Melbourne, Melbourne, Australia; the
  • Jamie E. Craig
    Department of Ophthalmology, Flinders University, Adelaide, Australia; and the
  • Catherine M. Green
    Center for Eye Research Australia, Royal Victorian Eye and Ear Hospital, University of Melbourne, Melbourne, Australia; the
  • Simon J. Foote
    Genetics and Bioinformatics Division, Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia; the
  • Paul N. Baird
    Center for Eye Research Australia, Royal Victorian Eye and Ear Hospital, University of Melbourne, Melbourne, Australia; the
  • Michèle M. Sale
    From the Menzies Research Institute, University of Tasmania, Hobart, Australia; the
    Center for Human Genomics and the
    Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, North Carolina.
Investigative Ophthalmology & Visual Science October 2005, Vol.46, 3723-3729. doi:https://doi.org/10.1167/iovs.05-0312
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Jac C. Charlesworth, Thomas D. Dyer, Jim M. Stankovich, John Blangero, David A. Mackey, Jamie E. Craig, Catherine M. Green, Simon J. Foote, Paul N. Baird, Michèle M. Sale; Linkage to 10q22 for Maximum Intraocular Pressure and 1p32 for Maximum Cup-to-Disc Ratio in an Extended Primary Open-Angle Glaucoma Pedigree. Invest. Ophthalmol. Vis. Sci. 2005;46(10):3723-3729. https://doi.org/10.1167/iovs.05-0312.

      Download citation file:


      © ARVO (1962-2015); The Authors (2016-present)

      ×
  • Supplements
Abstract

purpose. The purpose of this study was to identify genetic contributions to primary open-angle glaucoma (POAG) through investigations of two quantitative components of the POAG phenotype.

methods. Genome-wide multipoint variance–components 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.

Glaucoma is a major cause of visual impairment and the second leading cause of blindness worldwide. 1 The most common form is adult-onset primary open-angle glaucoma (POAG), which has a strong genetic component, with family history of the disease an acknowledged risk factor. 2 3 4 A 10-fold increase in risk of POAG has been documented in first-degree relatives of affected individuals, 3 whereas underreporting of family history suggests that the genetic component of POAG may be even greater than is generally acknowledged. 5  
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 POAG 6 7 has contributed to the difficulties in identifying genes involved in this disease. 
Two genes have been identified for adult-onset POAG 8 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 POAG 19 20 21 ; however, only GLC1B 22 and GLC1C 23 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 analysis 27 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
The Glaucoma Inheritance Study in Tasmania
This investigation was conducted as part of the Glaucoma Inheritance Study in Tasmania 35 (GIST), a large population study of glaucoma-affected families in Tasmania, Australia. Ethics approval was obtained from the Human Research Ethics Committees of the Royal Children’s Hospital, the Royal Victorian Eye and Ear Hospital, the Royal Hobart Hospital, and the University of Tasmania, and the study was conducted in accordance with the tenets of the Declaration of Helsinki. Written informed consent was obtained from all participants. 
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. 
The ascertainment correction was conducted with SOLAR (ver. 2.1.1) 25 and the population-based ascertainment correction dataset. The mean effects for each trait were estimated from this dataset, including the regression coefficients βAge and βSex, which were used to adjust the trait mean for the effect of the covariates age and sex. The resultant mean, SD, βAge, and βSex obtained from the population-based dataset were then used in the model to analyze the pedigree data. By constraining these mean effects and the total variance to that of the general population, trait values for individuals in family GTas02 were given the appropriate weight. 
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, age–sex interaction; however, age–sex 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 (σ2 P) into components attributable to covariate effects, effects of a specific QTL (QTL (σ2 q), and residual additive genetic effects (σ2 h2r). The test for linkage compares the likelihood of this model with the likelihood of a null model (no linkage), where the QTL effect size σq 2 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 complexity 25 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
Maximum recorded IOP had a heritability of 0.55 for family GTas02. The covariates sex and age were statistically controlled for and included in the analysis, although in the ascertainment correction population, neither showed effects that were substantial (P = 0.073 for sex, P = 0.53 for age), and accounted for only 0.1% of the variance. The mean value for maximum recorded IOP in family GTas02 (17.9 ± 4.2 mm Hg [SD]) was significantly higher (P < 0.0001) than the mean IOP in the general population (15.2 ± 3.3 mm Hg). The empiric LOD correction constant was greater than 1; hence, no adjustment was applied to the linkage results. The highest LOD score for IOP in family GTas02 was 3.3 (locus-specific P = 0.00015) near marker D10S537 (Fig. 1) , with the LOD-1 interval spanning approximately 20 cM on 10q22 (Fig. 2) . The genome-wide probability for this result was 0.0165. 
Maximum cup-to-disc ratio had a heritability of 0.39 for family GTas02. The effect of the covariates sex (P = 3.1 × 10−5) and age (P = 1.6 × 10−6) were highly significant in the ascertainment correction population, accounting for 1.4% of the variance; hence, both were statistically controlled for and included in the analysis of family GTas02. The mean ± SD cup-to-disc ratio in family GTas02 (0.47 ± 0.20) was not significantly different (P > 0.05) from the mean cup-to-disc ratio in the general population (0.43 ± 0.21). The empiric LOD correction constant was greater than 1; hence, no adjustment was applied. The highest LOD score for maximum cup-to-disc ratio in family GTas02 was 2.3 (locus-specific P = 0.00,056) near markers D1S197 and D1S220 (Fig. 3) , with the LOD-1 interval spanning approximately 20 cM on 1p32 (Fig. 4) . The genome-wide probability for this LOD score was determined to be 0.208. 
The MYOC Q368X mutation was present in 9 (37.5%) of 24 individuals in family GTas02 with diagnosed POAG and in 10 individuals without a diagnosis of POAG at examination. Whereas the mean age of individuals with the Q368X mutation in family GTas02 was not significantly different from those without the mutation, mean maximum IOP and mean maximum cup-to-disc ratio were significantly higher in Q368X mutation carriers, regardless of POAG affection status (Table 2) . MYOC Q368X mutation status (βMYOC) accounted for 16% of the genetic variance of the maximum IOP (heritability reduced from 0.55 to 0.46) and 49% of the genetic variance of the maximum cup-to-disc ratio (heritability reduced from 0.39 to 0.20). Addition of MYOC Q368X mutation status as a covariate in the analysis of the quantitative traits resulted in the peak LOD score for maximum recorded IOP (near marker D10S537 on 10q22) decreasing from 3.3 to 1.9 (Fig. 2) , and the peak LOD score for maximum cup-to-disc ratio (near markers D1S197 and D1S220 on 1p32) was reduced from 2.3 to 0.9 (Fig. 4)
Discussion
In this study, we investigated two individual disease components of POAG—maximum recorded IOP and maximum vertical cup-to-disc ratio—as quantitative traits using variance components linkage analysis. Genome scan analyses of these traits in an extended pedigree revealed one region of significant or suggestive linkage for each trait. Multipoint linkage analysis of maximum recorded IOP identified a peak LOD of 3.3 (locus-specific P = 0.00015) near marker D10S537 on 10q22, whereas analysis of the maximum cup-to-disc ratio produced a peak LOD score of 2.3 (locus-specific P = 0.00056) near markers D1S197 to D1S220 on 1p32. 
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 receptor—interesting 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 IOP 60 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 IOP’s 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 1p32—the first reported locus for this trait—with 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 dystrophy–based 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 dysgenesis 63 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 non–carrier-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. 
 
Table 1.
 
Age and Trait Data Distribution Statistics from the MVIP Dataset and Family GTas02
Table 1.
 
Age and Trait Data Distribution Statistics from the MVIP Dataset and Family GTas02
Trait MVIP Family GTas02 Fitted Coefficients
n Mean SD n Mean SD βAge βSex
Age (y) 3892 59.4 11.6 139 53.1 14.4
Maximum cup-to-disc ratio 2986 0.43 0.21 135 0.47 0.20 0.0016 −0.032
Maximum IOP (mm Hg) 3221 15.2 3.32 137 17.9 4.15 0.0033 −0.21
Figure 1.
 
Genome-wide multipoint variance components linkage results for maximum recorded IOP.
Figure 1.
 
Genome-wide multipoint variance components linkage results for maximum recorded IOP.
Figure 2.
 
Multipoint variance–components linkage results for maximum recorded IOP for chromosome 10 (solid line) and linkage signal after the inclusion of Q368X status as a covariate (dashed line). Covariates are indicated on the plots. The location of the optineurin (OPTN) gene is indicated.
Figure 2.
 
Multipoint variance–components linkage results for maximum recorded IOP for chromosome 10 (solid line) and linkage signal after the inclusion of Q368X status as a covariate (dashed line). Covariates are indicated on the plots. The location of the optineurin (OPTN) gene is indicated.
Figure 3.
 
Genome-wide multipoint variance components linkage results for maximum vertical cup-to-disc ratio.
Figure 3.
 
Genome-wide multipoint variance components linkage results for maximum vertical cup-to-disc ratio.
Figure 4.
 
Multipoint variance components linkage results for maximum vertical cup-to-disc ratio for chromosome 1 (solid line), and linkage signal after the inclusion of Q368X status as a covariate (dashed line). Covariates are indicated on the plots. The location of the myocilin (MYOC) gene is indicated.
Figure 4.
 
Multipoint variance components linkage results for maximum vertical cup-to-disc ratio for chromosome 1 (solid line), and linkage signal after the inclusion of Q368X status as a covariate (dashed line). Covariates are indicated on the plots. The location of the myocilin (MYOC) gene is indicated.
Table 2.
 
Comparison between Myocilin Q368X Mutation Carriers and Mutation-Free Individuals from Family GTas02
Table 2.
 
Comparison between Myocilin Q368X Mutation Carriers and Mutation-Free Individuals from Family GTas02
Trait Q368X Mutation Carriers Q368X Mutation Free P *
n Mean SD n Mean SD
Maximum IOP 19 22.21 5.32 118 17.17 3.48 7.70 × 10−8
Maximum cup-to-disc ratio 19 0.64 0.23 116 0.44 0.18 2.26 × 10−4
The authors thank the GTas02 family members, the population control subjects, and the MVIP participants for their participation; Hien Vu (Population Health Division of the University of Melbourne) for providing the MVIP population data; and Danielle Healy, Susan Stanwix, Tiffany Wong, Maree Ring, Julie Barbour, Robin Wilkinson, Colleen Wilkinson, Robert Buttery, Andrew McNaught, Michael Coote, and Julian Rait for assistance with examination and data collection. 
QuigleyHA. Number of people with glaucoma worldwide. Br J Ophthalmol. 1996;80:389–393. [CrossRef] [PubMed]
WeihLM, NanjanM, McCartyCA, TaylorHR. Prevalence and predictors of open-angle glaucoma: results from the visual impairment project. Ophthalmology. 2001;108:1966–1972. [CrossRef] [PubMed]
WolfsRC, KlaverCC, RamrattanRS, van DuijnCM, HofmanA, de JongPT. Genetic risk of primary open-angle glaucoma. Population-based familial aggregation study. Arch Ophthalmol.. 1998;116:1640–1645. [CrossRef] [PubMed]
RosenthalAR, PerkinsES. Family studies in glaucoma. Br J Ophthalmol. 1985;69:664–667. [CrossRef] [PubMed]
McNaughtAI, AllenJG, HealeyDL, et al. Accuracy and implications of a reported family history of glaucoma: experience from the Glaucoma Inheritance Study in Tasmania. Arch Ophthalmol. 2000;118:900–904. [PubMed]
WolfsRC, BorgerPH, RamrattanRS, et al. Changing views on open-angle glaucoma: definitions and prevalences–The Rotterdam Study. Invest Ophthalmol Vis Sci. 2000;41:3309–3321. [PubMed]
KahnHA, MiltonRC. Alternative definitions of open-angle glaucoma. Effect on prevalence and associations in the Framingham eye study. Arch Ophthalmol.. 1980;98:2172–21777. [CrossRef] [PubMed]
RezaieT, ChildA, HitchingsR, et al. Adult-onset primary open-angle glaucoma caused by mutations in optineurin. Science. 2002;295:1077–1079. [CrossRef] [PubMed]
StoneEM, FingertJH, AlwardWL, et al. Identification of a gene that causes primary open angle glaucoma. Science. 1997;275:668–670. [CrossRef] [PubMed]
AlwardWL, KwonY, KhannaC, et al. Variations in the myocilin gene in patients with open-angle glaucoma. Arch Ophthalmol. 2002;120:1189–1197. [CrossRef] [PubMed]
FingertJH, HeonE, LiebmannJM, et al. Analysis of myocilin mutations in 1703 glaucoma patients from five different populations. Hum Mol Genet. 1999;8:899–905. [CrossRef] [PubMed]
AlwardWL, FingertJH, CooteMA, et al. Clinical features associated with mutations in the chromosome 1 open- angle glaucoma gene (GLC1A). N Engl J Med. 1998;338:1022–1027. [CrossRef] [PubMed]
AlwardWL, KwonYH, KawaseK, et al. Evaluation of optineurin sequence variations in 1,048 patients with open-angle glaucoma. Am J Ophthalmol. 2003;136:904–910. [CrossRef] [PubMed]
MonemiS, SpaethG, DasilvaA, et al. Identification of a novel adult-onset primary open angle glaucoma (POAG) gene on 5q22.1. Hum Mol Genet. 2005;14:725–733. [CrossRef] [PubMed]
StoilovaD, ChildA, TrifanOC, CrickRP, CoakesRL, SarfaraziM. Localization of a locus (GLC1B) for adult-onset primary open angle glaucoma to the 2cen-q13 region. Genomics. 1996;36:142–150. [CrossRef] [PubMed]
WirtzMK, SamplesJR, KramerPL, et al. Mapping a gene for adult-onset primary open-angle glaucoma to chromosome 3q. Am J Hum Genet. 1997;60:296–304. [PubMed]
TrifanOC, TraboulsiEI, StoilovaD, et al. A third locus (GLC1D) for adult-onset primary open-angle glaucoma maps to the 8q23 region. Am J Ophthalmol. 1998;126:17–28. [CrossRef] [PubMed]
WirtzMK, SamplesJR, RustK, et al. GLC1F, a new primary open-angle glaucoma locus, maps to 7q35–q36. Arch Ophthalmol. 1999;117:237–241. [CrossRef] [PubMed]
WiggsJL, AllinghamRR, HossainA, et al. Genome-wide scan for adult onset primary open angle glaucoma. Hum Mol Genet. 2000;9:1109–1117. [CrossRef] [PubMed]
NemesureB, JiaoX, HeQ, et al. A genome-wide scan for primary open-angle glaucoma (POAG): the Barbados Family Study of Open-Angle Glaucoma. Hum Genet. 2003;112:600–609. [PubMed]
WiggsJL, LynchS, YnagiG, et al. A genomewide scan identifies novel early-onset primary open-angle glaucoma loci on 9q22 and 20p12. Am J Hum Genet. 2004;74:1314–1320. [CrossRef] [PubMed]
CharlesworthJC, StankovichJM, MackeyDA, et al. Confirmation of the adult-onset primary open angle glaucoma locus GLC1B at 2cen-q13 in an Australian family. Ophthalmologica.. .In press
KitsosG, EibergH, Economou-PetersenE, et al. Genetic linkage of autosomal dominant primary open angle glaucoma to chromosome 3q in a Greek pedigree. Eur J Hum Genet. 2001;9:452–457. [CrossRef] [PubMed]
FreimerN, SabattiC. The use of pedigree, sib-pair and association studies of common diseases for genetic mapping and epidemiology. Nat Genet. 2004;36:1045–1051. [CrossRef] [PubMed]
AlmasyL, BlangeroJ. Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998;62:1198–1211. [CrossRef] [PubMed]
AmosCI, de AndradeM. Genetic linkage methods for quantitative traits. Stat Methods Med Res. 2001;10:3–25. [CrossRef] [PubMed]
WilliamsJT, BlangeroJ. Comparison of variance components and sibpair-based approaches to quantitative trait linkage analysis in unselected samples. Genet Epidemiol. 1999;16:113–134. [CrossRef] [PubMed]
DuggiralaR, WilliamsJT, Williams-BlangeroS, BlangeroJ. A variance component approach to dichotomous trait linkage analysis using a threshold model. Genet Epidemiol. 1997;14:987–992. [CrossRef] [PubMed]
WilliamsJT, BlangeroJ. Power of variance component linkage analysis to detect quantitative trait loci. Ann Hum Genet. 1999;63:545–563. [CrossRef] [PubMed]
WijsmanEM, AmosCI. Genetic analysis of simulated oligogenic traits in nuclear and extended pedigrees: summary of GAW10 contributions. Genet Epidemiol. 1997;14:719–735. [CrossRef] [PubMed]
KleinBE, KleinR, LeeKE. Heritability of risk factors for primary open-angle glaucoma: the Beaver Dam Eye Study. Invest Ophthalmol Vis Sci. 2004;45:59–62. [CrossRef] [PubMed]
ViswanathanAC, HitchingsRA, IndarA, et al. Commingling analysis of intraocular pressure and glaucoma in an older Australian population. Ann Hum Genet. 2004;68:489–497. [CrossRef] [PubMed]
DuggalP, KleinAP, LeeKE, et al. A genetic contribution to intraocular pressure: the Beaver Dam Eye Study. Invest Ophthalmol Vis Sci. 2005;46:555–560. [CrossRef] [PubMed]
BairdPN, FooteSJ, MackeyDA, CraigJE, SpeedTP, BureauA. Evidence for a novel glaucoma locus at chromosome 3p21-22. Hum Genet. 2005;117:249–257. [CrossRef] [PubMed]
MackeyDA, CraigJE. Glaucoma Inheritance Study in Tasmania: An International Collaboration. Basic and Clinical Science Course Section 13. 2002;265–269.American Academy of Ophthalmology San Francisco.
CraigJE, BairdPN, HealeyDL, et al. Evidence for genetic heterogeneity within eight glaucoma families, with the GLC1A Gln368STOP mutation being an important phenotypic modifier. Ophthalmology. 2001;108:1607–1620. [CrossRef] [PubMed]
CooteMA, McCartneyPJ, WilkinsonRM, MackeyDA. The ‘GIST’ score: ranking glaucoma for genetic studies. Glaucoma Inheritance Study of Tasmania. Ophthalmic Genet.. 1996;17:199–208. [CrossRef] [PubMed]
BairdPN, CraigJE, RichardsonAJ, et al. Analysis of 15 primary open-angle glaucoma families from Australia identifies a founder effect for the Q368STOP mutation of myocilin. Hum Genet. 2003;112:110–116. [PubMed]
O’ConnellJR, WeeksDE. PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet. 1998;63:259–266. [CrossRef] [PubMed]
VickersJC, CraigJE, StankovichJ, et al. The apolipoprotein epsilon4 gene is associated with elevated risk of normal tension glaucoma. Mol Vis. 2002;8:389–393. [PubMed]
HeathSC. Markov chain Monte Carlo segregation and linkage analysis for oligogenic models. Am J Hum Genet. 1997;61:748–760. [CrossRef] [PubMed]
HeathSC, SnowGL, ThompsonEA, TsengC, WijsmanEM. MCMC segregation and linkage analysis. Genet Epidemiol. 1997;14:1011–1016. [CrossRef] [PubMed]
WensorMD, McCartyCA, StanislavskyYL, LivingstonPM, TaylorHR. The prevalence of glaucoma in the Melbourne Visual Impairment Project. Ophthalmology. 1998;105:733–739. [CrossRef] [PubMed]
BlangeroJ, WilliamsJT, AlmasyL. Variance component methods for detecting complex trait loci. Adv Genet. 2001;42:151–181. [PubMed]
BlangeroJ, AlmasyL. Multipoint oligogenic linkage analysis of quantitative traits. Genet Epidemiol. 1997;14:959–964. [CrossRef] [PubMed]
WilliamsJT, Van EerdeweghP, AlmasyL, BlangeroJ. Joint multipoint linkage analysis of multivariate qualitative and quantitative traits. I. Likelihood formulation and simulation results. Am J Hum Genet.. 1999;65:1134–1147. [CrossRef] [PubMed]
WilliamsJT, BegleiterH, PorjeszB, et al. Joint multipoint linkage analysis of multivariate qualitative and quantitative traits. II. Alcoholism and event-related potentials. Am J Hum Genet.. 1999;65:1148–1160. [CrossRef] [PubMed]
MartinLJ, MahaneyMC, AlmasyL, et al. A quantitative trait locus on chromosome 22 for serum leptin levels adjusted for serum testosterone. Obes Res. 2002;10:602–607. [CrossRef] [PubMed]
ComuzzieAG, HixsonJE, AlmasyL, et al. A major quantitative trait locus determining serum leptin levels and fat mass is located on human chromosome 2. Nat Genet. 1997;15:273–276. [CrossRef] [PubMed]
BlangeroJ, WilliamsJT, AlmasyL. Robust LOD scores for variance component-based linkage analysis. Genet Epidemiol. 2000;19(suppl 1)S8–S14. [CrossRef] [PubMed]
FeingoldE, BrownPO, SiegmundD. Gaussian models for genetic linkage analysis using complete high-resolution maps of identity by descent. Am J Hum Genet. 1993;53:234–251. [PubMed]
WuZ, NakuraJ, AbeM, et al. Genome-wide linkage disequilibrium mapping of hypertension in Japan. Hypertens Res. 2003;26:533–540. [CrossRef] [PubMed]
DielemansI, VingerlingJR, AlgraD, HofmanA, GrobbeeDE, de JongPT. Primary open-angle glaucoma, intraocular pressure, and systemic blood pressure in the general elderly population. The Rotterdam Study. Ophthalmology.. 1995;102:54–60. [CrossRef] [PubMed]
TielschJM, KatzJ, SommerA, QuigleyHA, JavittJC. Hypertension, perfusion pressure, and primary open-angle glaucoma: a population-based assessment. Arch Ophthalmol. 1995;113:216–221. [CrossRef] [PubMed]
McLeodSD, WestSK, QuigleyHA, FozardJL. A longitudinal study of the relationship between intraocular and blood pressures. Invest Ophthalmol Vis Sci. 1990;31:2361–2366. [PubMed]
ChangB, SmithRS, PetersM, et al. Haploinsufficient Bmp4 ocular phenotypes include anterior segment dysgenesis with elevated intraocular pressure. BMC Genet. 2001;2:18. [CrossRef] [PubMed]
BernalS, CalafM, Garcia-HoyosM, et al. Study of the involvement of the RGR, CRPB1, and CRB1 genes in the pathogenesis of autosomal recessive retinitis pigmentosa. J Med Genet. 2003;40:E89. [CrossRef] [PubMed]
HansenFK, EhlersN. Elevated tonometer readings caused by a thick cornea. Acta Ophthalmol (Copenh). 1971;49:775–778. [PubMed]
StodtmeisterR. Applanation tonometry and correction according to corneal thickness. Acta Ophthalmol Scand. 1998;76:319–324. [CrossRef] [PubMed]
WolfsRC, KlaverCC, VingerlingJR, GrobbeeDE, HofmanA, de JongPT. Distribution of central corneal thickness and its association with intraocular pressure: The Rotterdam Study. Am J Ophthalmol. 1997;123:767–772. [CrossRef] [PubMed]
DohadwalaAA, DamjiKF. Familial occurrence of artefactual ocular hypertension from thick corneas and of primary open angle glaucoma in a French Canadian kindred. Ophthalmic Genet. 2000;21:1–7. [CrossRef] [PubMed]
ManyaH, SakaiK, KobayashiK, et al. Loss-of-function of an N-acetylglucosaminyltransferase, POMGnT1, in muscle-eye-brain disease. Biochem Biophys Res Commun. 2003;306:93–97. [CrossRef] [PubMed]
SeminaEV, BrownellI, Mintz-HittnerHA, MurrayJC, JamrichM. Mutations in the human forkhead transcription factor FOXE3 associated with anterior segment ocular dysgenesis and cataracts. Hum Mol Genet. 2001;10:231–236. [CrossRef] [PubMed]
AkarsuAN, TuracliME, AktanSG, et al. A second locus (GLC3B) for primary congenital glaucoma (Buphthalmos) maps to the 1p36 region. Hum Mol Genet. 1996;5:1199–1203. [CrossRef] [PubMed]
AllinghamRR, WiggsJL, De La PazMA, et al. Gln368STOP myocilin mutation in families with late-onset primary open- angle glaucoma. Invest Ophthalmol Vis Sci. 1998;39:2288–2295. [PubMed]
Figure 1.
 
Genome-wide multipoint variance components linkage results for maximum recorded IOP.
Figure 1.
 
Genome-wide multipoint variance components linkage results for maximum recorded IOP.
Figure 2.
 
Multipoint variance–components linkage results for maximum recorded IOP for chromosome 10 (solid line) and linkage signal after the inclusion of Q368X status as a covariate (dashed line). Covariates are indicated on the plots. The location of the optineurin (OPTN) gene is indicated.
Figure 2.
 
Multipoint variance–components linkage results for maximum recorded IOP for chromosome 10 (solid line) and linkage signal after the inclusion of Q368X status as a covariate (dashed line). Covariates are indicated on the plots. The location of the optineurin (OPTN) gene is indicated.
Figure 3.
 
Genome-wide multipoint variance components linkage results for maximum vertical cup-to-disc ratio.
Figure 3.
 
Genome-wide multipoint variance components linkage results for maximum vertical cup-to-disc ratio.
Figure 4.
 
Multipoint variance components linkage results for maximum vertical cup-to-disc ratio for chromosome 1 (solid line), and linkage signal after the inclusion of Q368X status as a covariate (dashed line). Covariates are indicated on the plots. The location of the myocilin (MYOC) gene is indicated.
Figure 4.
 
Multipoint variance components linkage results for maximum vertical cup-to-disc ratio for chromosome 1 (solid line), and linkage signal after the inclusion of Q368X status as a covariate (dashed line). Covariates are indicated on the plots. The location of the myocilin (MYOC) gene is indicated.
Table 1.
 
Age and Trait Data Distribution Statistics from the MVIP Dataset and Family GTas02
Table 1.
 
Age and Trait Data Distribution Statistics from the MVIP Dataset and Family GTas02
Trait MVIP Family GTas02 Fitted Coefficients
n Mean SD n Mean SD βAge βSex
Age (y) 3892 59.4 11.6 139 53.1 14.4
Maximum cup-to-disc ratio 2986 0.43 0.21 135 0.47 0.20 0.0016 −0.032
Maximum IOP (mm Hg) 3221 15.2 3.32 137 17.9 4.15 0.0033 −0.21
Table 2.
 
Comparison between Myocilin Q368X Mutation Carriers and Mutation-Free Individuals from Family GTas02
Table 2.
 
Comparison between Myocilin Q368X Mutation Carriers and Mutation-Free Individuals from Family GTas02
Trait Q368X Mutation Carriers Q368X Mutation Free P *
n Mean SD n Mean SD
Maximum IOP 19 22.21 5.32 118 17.17 3.48 7.70 × 10−8
Maximum cup-to-disc ratio 19 0.64 0.23 116 0.44 0.18 2.26 × 10−4
×
×

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.

×