|
|
||||||||
1From the Statistical Genetics Section, Inherited Disease Research Branch, National Human Genome Research Institute, National Institutes of Health, Baltimore, Maryland; the 2Department of Ophthalmology and Visual Sciences, University of Wisconsin Medical School, Madison, Wisconsin; and the 3Department of Epidemiology and Biostatistics, Case Western Reserve University, Cleveland, Ohio.
| Abstract |
|---|
|
|
|---|
METHODS. Segregation analysis is an analytical method that provides statistical evidence supporting the involvement of a major gene or polygenes in a particular phenotype. Detailed medical histories and eye examinations were performed on all participants. From the two eyes, the higher IOP measurement was used as a continuous trait after adjustment for covariates. A genome-wide scan (GWS) using affected sib pair linkage analysis was performed on 218 sibling pairs.
RESULTS. In this segregation analysis the model that allowed for an unmeasured major environmental effect plus a polygenic/multifactorial effect provided the best fit and was the most parsimonious model. The lack of an adequate fit for the Mendelian single-gene models is consistent with a multifactorial model of inheritance that may include multiple genes and environmental factors that contribute to IOP. The results of the GWS yielded two novel loci as potential linkage regions for IOP on chromosomes 6 (P = 0.008) and 13 (P = 0.0007). Neither of these regions has previously been identified in GWS of POAG.
CONCLUSIONS. The segregation and familial correlation analyses of IOP suggest a polygenetic component with environmental influences. The pilot linkage study further confirms the heterogeneity of IOP with the identification of two novel genetic loci.
2.0) to loci on chromosomes 2, 14, 17 (short arm), and 19, using 113 sibling pairs.7 Recently, a genome-wide scan of an Afro-Caribbean population provided evidence for linkage to POAG on 2q and 10p.8 Two loci, TIGR/myocilin (1q23) and OPTN/optineurin (10p14) have been identified9 10 and >43 mutations have been identified in myocilin in patients with POAG.11 These studies support POAG as a genetically heterogeneous disorder. Elevated intraocular pressure (IOP) is a principal risk factor for POAG and results in optic nerve damage. In the United States, it is estimated that 3 to 6 million people have elevated IOP and are at an increased risk for development of glaucoma.12 The goal of this study was to investigate a potential genetic contribution to IOP, a predictor of POAG. We performed a complex segregation analysis on 2336 individuals in 620 extended pedigrees ascertained through a population-based cohort, the Beaver Dam Eye Study (BDES) in Beaver Dam, Wisconsin. In addition, we performed a nonparametric genome-wide linkage analysis of 218 sibling pairs from the BDES.
| Methods |
|---|
|
|
|---|
Data Collection
IOP was measured with a Goldmann applanation tonometer after instillation of a drop of fluorescein (Fluoress; Barnes-Hind Armour Pharmaceutical Co., Kankakee, IL) in each eye. Before measuring the IOP in the right eye, the tonometer was set to 10 mm Hg and the pressure was recorded only after the tonometer was moved back from the cornea. The tonometer was then reset to 10 mm Hg for measurement of the left eye. Any measurement the examiner thought was unreliable was excluded from the analysis. In addition to a full eye examination each participant also provided a detailed medical history and examination, which included information about hypertension, diabetes, and other medical conditions as well as a history of medication use. Although longitudinal eye measurements are collected, only baseline examination IOP measurements were used in this analysis.
Statistical Methods
Preliminary Analysis.
Linear regression was used to assess potential confounding variables (SAS, ver. 8.1; SAS, Cary, NC). IOP was defined as the higher IOP measurement between each eye and was the dependent variable. Age, sex, systolic blood pressure (SBP), and treatment of IOP were considered potential confounders and were either included as covariates or adjusted for before further analysis. Familial correlations including parent-offspring, sibling, cousin, avuncular, and spouse pairs were estimated with the FCOR program of S.A.G.E. (1987; Statistical Analysis for Genetic Epidemiology). We estimated the heritability (h2) from the adjusted sibling correlations (r), by the equation h2 = 2r.
Segregation Analysis.
The REGC program of S.A.G.E. (1987) version 3.0 was used to perform complex segregation analysis. In addition, the REGCHUNT program15 which allows REGC to run multiple times with random starting values was used. Using the regressive models proposed by Bonney16 17 (S.A.G.E., 1994), REGC performs segregation analysis of continuous traits. Because the BDES is a population-based survey, no ascertainment correction was necessary.
Mendelian (autosomal dominant, autosomal codominant, and autosomal recessive), sporadic, and environmental models were used to test for a locus that influences IOP. Maximum-likelihood methods were used to estimate the parameters for each model. Each model was compared to a general, unrestricted model by using the likelihood ratio test. Models can include different discrete factors that may affect a persons phenotype or IOP. These are called "type"18 or "ousiotype."19 In the Mendelian models, type is represented by genotype, and three types are possible: AA (homozygote), AB (heterozygote), and BB (homozygote). However, in the non-Mendelian models, type may represent any transmission, including shared exposures to environmental risk factors within families. If the type frequencies are in Hardy-Weinberg equilibrium, then they are defined as the frequency of allele A, qA. The probability of a parents transmitting these discrete factors to an offspring is represented by a transmission parameter (
). Under Mendelian transmission the
s are interpreted as the probability that a parent transmits the A allele to an offspring. Thus,
AA = 1,
AB = 0.5, and
BB = 0. However, when these transmission factors are fixed equal to each other, they represent an environmental model in which parental type (environmental exposure) has no effect on offspring type.
Using the S.A.G.E. statistical program, we performed all analyses with class D models that assume, after conditioning on the parents, that the correlations among siblings are equal but that they are not completely due to common parentage. Additional familial correlations can therefore be estimated that may explain the effects of other genes (polygenic inheritance) or common environmental factors. These correlations may include spousal (
fm), parent-offspring (
po), and sibling (
ss). In addition, age, sex, treatment of IOP, and SBP measurements at baseline were included as covariates in the analysis because they were considered potential confounders. The IOP measurements were initially analyzed both with and without transformation to normality by using the Box-Cox transformation.20 However, the IOP measurements followed a normal distribution, and the data were analyzed without the need for transformation.
Six hypotheses were tested against the likelihood of a general or unrestricted model and likelihood ratio tests, and Akaikes Information Criterion (AIC) was used to select the most parsimonious model. The six hypotheses of transmission included (1) no major gene, (2) Mendelian dominant, (3) Mendelian codominant, (4) Mendelian recessive (arbitrary major gene), (5) environmental mixed model with a major environmental effect (no parent-offspring transmission) which estimates the transmission probabilities to be equal to the allele frequency, and (6) a semigeneral or
AB estimated model that estimates the parameter
AB when
AA and
BB are fixed to the Mendelian values of 1.0 and 0.0, respectively. Likelihood ratio tests were computed by 2 times the difference in lnLikelihood of the general model and a nested model. AIC is defined as 2lnL + 2(number of parameters estimated).21 The model with the minimum AIC is the most parsimonious model.
Pilot Linkage Analysis.
Previously, 325 individuals (263 sibling pairs) from 102 pedigrees in the BDES were selected for genotyping and linkage analysis of 345 genome-wide markers (Weber panel 8 marker set) because some members of these families were affected with age-related maculopathy (ARM).22 Although ARM is unrelated to IOP, many individuals genotyped (218 sibling pairs within 87 pedigrees) had complete IOP and covariate measurements and could be used in a pilot genome-wide linkage analysis.
We determined the allele frequencies of all the markers by using maximum-likelihood methods in the S.A.G.E. program FREQ (S.A.G.E., ver. 4.5). MARKERINFO (S.A.G.E., ver. 4.5) identified two sibships with Mendelian inheritance inconsistencies. These sibships were then excluded from the remaining analyses. RELTEST (S.A.G.E., ver. 4.5), RELCHECK,23 24 and PREST25 were used to identify potential misspecified family relationships. Based on these results, data from an additional 13 sibships were excluded. IOP was treated as a quantitative trait, and we used the higher untransformed IOP measurement in either eye, adjusted for age, sex, IOP treatment, and SBP as the dependent variable in a Haseman-Elston regression. Using GENIBD, we estimated the single point and multipoint identity by descent-sharing for the study subjects. Linkage analysis was performed with the modified Haseman-Elston regression models and the option that allows for the nonindependence of sibling pairs and uses a weighted combination of the squared trait difference and squared mean corrected trait sum (option W4)26 26 in the model-free linkage program SIBPAL (S.A.G.E., ver. 4.5). For each marker, P-values were obtained using the asymptotic distribution of the likelihood-ratio test statistics. Each chromosome was analyzed separately.
| Results |
|---|
|
|
|---|
|
AB free model 6; P < 0.0001). The environmental mixed model with
s equal to the allele frequency and multiple distributions (model 5) with polygenic effects provided a better fit to the data and could not be rejected (P = 0.12). However, both the environmental mixed model
s equal to the allele frequency and the
AB free model were highly parsimonious according to the AIC (11974.0469 and 11974.2469, respectively).
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
In the segregation analysis, the model that allowed for an unmeasured major environmental effect plus a polygenic/environmental effect provided the best fit to the data and was the most parsimonious model. Another version of the general model (semigeneral) commonly used to test for the presence of a major locus effect in a mixed model is one in which the
of the A allele for AA and BB genotypes are fixed to 1.0 and 0, and the remaining parameters are unrestricted. This model was marginally rejected (P = 0.02) and highly parsimonious and yielded an estimate of the
AB transmission of 0.70, close to its Mendelian value of 0.5. All the Mendelian models tested for only a single major gene to control IOP across the spectrum of IOP values. The lack of an adequate fit for these models is consistent with a multifactorial or complex model of inheritance that may include multiple genes (at least more than one gene) and environmental factors that contribute to IOP. The conclusions of the segregation analysis are not surprising, because IOP has been associated with numerous physiological and environmental factors including SBP and diastolic (DBP) blood pressure, body mass index, hematocrit, serum glucose, glycohemoglobin, cholesterol levels, pulse, nuclear sclerosis, and season and time of day of measurement.31 As with most genetic analyses it is difficult to account for all potential confounders because of the complexity of the genetic model and current limitations of analytical programs. However, we tried to account for those covariates that may have the greatest effect on IOP, but there may be some confounding by unaccounted for covariates.
The inclusion of SBP as a covariate is important, because a recent study of the BDES by Klein et al.32 determined that IOPs correlated with SBP and DBP at both baseline and 5-year follow-up. Specifically, they found that reduced systemic blood pressure was associated with reduced IOP.
In addition, whole-genome scans for hypertension including both SBP and DBP have also indicated a polygenic inheritance with heterogeneity.22 33 34 35 36 37 38 39 40 41 42 Our analysis of IOP appears to be consistent with blood pressure and is suggestive that we are looking at a similar pattern of inheritance. It is plausible that IOP, like blood pressure, is heterogeneous, and there may be loci that control both the pressure in the eye and blood.
The two regions we identified by our pilot linkage analysis, chromosomes 6 and 13, warrant further attention, especially since neither of these loci has been identified in previous genome-wide scans of POAG. The identification of two new regions for POAG may be due to the definition of IOP and the ascertainment of cases. This study analyzed the quantitative trait, IOP, compared with the dichotomous POAG previously studied in linkage analyses of glaucoma. Although elevated IOP is a key characteristic in the diagnosis of glaucoma, IOP may represent only one aspect of development of glaucoma. Further localization of these novel regions will provide greater insight into the mechanism of IOP and the contribution to POAG. In addition, the ascertainment of elevated IOP individuals in the Beaver Dam population is in contrast to some of the previous linkage studies involving individuals with POAG from highly aggregated families. The previous linkage studies selected families with well-defined cases of POAG and may have identified rare highly penetrant susceptibility alleles. The BDES involved a predominantly European-American community-based sample which ascertained all age-eligible individuals in Beaver Dam, Wisconsin. Thus, it is possible this sample lacked the power to replicate the initial linkage findings for rare and highly penetrant alleles, and instead we identified more common susceptibility alleles at two potential loci that may control the variation in IOP. Finally, it is important to note that many of the previous linkage findings have not yet been replicated, nor has a gene been identified for each identified locus, and so it is likely some of these regions are not true disease loci. However, the replication of a true linkage finding can be difficult if there are multiple genes that contribute to a disease43 as with POAG, which may help to explain the lack of replication in our study and others.
Although we identified a strong two-point linkage signal, the strength of the signal decreased in the multipoint analysis. The weakened signal for linkage in the multipoint analysis is not surprising, since, for both chromosome 6 and 13, there is a wide distance surrounding the marker that provided the greatest linkage signal (25 cM and 28 cM, respectively). Further fine mapping in this region may identify a narrower linkage peak. A whole-genome linkage analysis is currently under way and involves the complete family structure from the BDES (2231 individuals). This should give us additional power to identify new genes or further localize novel genes that contribute to the complex traits IOP and glaucoma.
| Footnotes |
|---|
Submitted for publication June 21, 2004; revised October 13, 2004; accepted October 19, 2004.
Disclosure: P. Duggal, None; A.P. Klein, None; K.E. Lee, None; S.K. Iyengar, None; R. Klein, None; J.E. Bailey-Wilson, None; B.E.K. Klein, 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: Barbara E. K. Klein, Department of Ophthalmology and Visual Sciences, University of Wisconsin-Madison, 610 North Walnut St., 4th Floor WARF, Madison, WI 53726-2336; kleinb{at}epi.ophth.wisc.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
P. Duggal, A. P. Klein, K. E. Lee, R. Klein, B. E. K. Klein, and J. E. Bailey-Wilson Identification of Novel Genetic Loci for Intraocular Pressure: A Genomewide Scan of the Beaver Dam Eye Study Arch Ophthalmol, January 1, 2007; 125(1): 74 - 79. [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] |
||||
![]() |
C. N. Rotimi, G. Chen, A. A. Adeyemo, L. S. Jones, K. Agyenim-Boateng, B. A. Eghan Jr, J. Zhou, A. Doumatey, K. Lashley, H. Huang, et al. Genomewide scan and fine mapping of quantitative trait Loci for intraocular pressure on 5q and 14q in west africans. Invest. Ophthalmol. Vis. Sci., August 1, 2006; 47(8): 3262 - 3267. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Klein and B. E. K. Klein The Epidemiology of Eye Disease: From Glycemia to Genetics The Friedenwald Lecture Invest. Ophthalmol. Vis. Sci., May 1, 2006; 47(5): 1747 - 1753. [Full Text] [PDF] |
||||
![]() |
J. C. Charlesworth, T. D. Dyer, J. M. Stankovich, J. Blangero, D. A. Mackey, J. E. Craig, C. M. Green, S. J. Foote, P. N. Baird, and M. 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., October 1, 2005; 46(10): 3723 - 3729. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |