|
|
||||||||
1From the National Human Genome Center and the 2Department of Ophthalmology, Howard University, College of Medicine, Washington, DC; the 3Department of Medicine, University of Science and Technology, Kumasi, Ghana; the 4College of Medicine, Endocrine and Metabolic Unit, University of Lagos, Lagos, Nigeria; the 5Department of Medicine and Ophthalmology, University of Nigeria Teaching Hospital, Enugu, Nigeria; and the 6Department of Medicine and Surgery, University of Ghana Medical School, Accra, Ghana.
| Abstract |
|---|
|
|
|---|
METHODS. Genomewide linkage analysis was conducted for loci linked to IOP (measured by applanation tonometry) in 244 affected sibling pairs with T2D using 372 autosomal short-tandem repeat markers at an average spacing of 9 cM.
RESULTS. Multipoint variance components linkage analyses revealed suggestive linkage on chromosome 5 (5q22) with a logarithm of odds (LOD) score of 2.50 (nominal P = 0.0003; empiric P = 0.0004) and on chromosome 14 (14q22) with an LOD score of 2.95 (nominal P = 0.0001; empiric P = 0.0003). Fine mapping at a marker density of 2 cM in the 5q region confirmed the linkage signal, with an increase in peak LOD score to 4.91.
CONCLUSIONS. The strong signal on chromosome 5 lies in the region in which a novel gene, WDR36, in the GLC1G locus was recently identified as causative for adult-onset primary open-angle glaucoma and provides additional evidence that chromosome 5 contains susceptibility loci for glaucoma in multiple human populations. The evidence provided in this study is particularly important given the evolutionary history of these West African populations and the recent ancestral relationship to African Americansa population with one of the highest rates of diabetes and associated complications (including glaucoma) in the world.
Elevated intraocular pressure (IOP) is the strongest known risk factor for glaucoma. Experimental elevation of IOP can induce glaucomatous neuropathy.5 A recent meta-analysis of randomized, controlled trials showed that lowering IOP in patients with ocular hypertension or glaucoma helps to reduce the long-term risk of visual field loss.6 Given that glaucoma can and does occur without high IOP (called normal tension glaucoma)7 and that most persons with high IOP do not develop glaucoma,8 it may be that there is considerable individual variation in disease susceptibility. Several studies also show variation in disease progression and response to treatment in patients with open-angle glaucoma.9 10
The genetic factors determining IOP level and susceptibility to pressure-induced damage are largely unknown, as are the specific processes that lead to retinal ganglion cell death. Several studies provide evidence for a significant genetic contribution.5 In fact, positive family history has long been recognized as a major risk factor for glaucoma, with the risk for glaucoma in first-degree relatives estimated to be as high as 7 to 10 times that in the general population.11 12 13 14 15 Several studies are directed at identifying the chromosomal regions and genes and the associated risk factors, including IOP, that contribute to glaucoma.4 15 16 17 18 19 At least seven genetic lociGLC1A,20 GLC1B,16 GLC1C,18 GLC1D,21 GLC1E,22 GLC1F,23 and another chromosome 2 locus15 have been reported for POAG. However, only three genesmyocilin (MYOC),17 optineurin (OPTN),24 and WD40-repeat 36 (WDR36)4 have been reported to be causally linked to glaucoma in these loci.
To our knowledge, no genomewide linkage study for POAG has been conducted in sub-Saharan African populations, and only one25 has been conducted in any population of African descent living on other continents.19 In this report, we provide evidence of quantitative trait loci (QTL) for IOP in West Africa (a major source population for present-day African Americans) and show strong evidence of replication of a recently reported POAG locus.4
| Materials and Methods |
|---|
|
|
|---|
All subjects were affected sibling pairs with T2D and were enrolled from five centers in two countries in West Africa: Enugu, Ibadan, and Lagos in Nigeria and Accra and Kumasi in Ghana. Diabetes diagnosis was based on the criteria established by the America Diabetes Association Expert Committee, as follows: fasting plasma glucose (FPG) concentration greater than 126 mg/dL (7.0 mM) or a 2-hour postload value in oral glucose tolerance test (OGTT) greater than 200 mg/dL (11.1 mM) on more than one occasion. The detection of autoantibodies to glutamic acid decarboxylase (GAD) antibody or a fasting C-peptide level less than 0.03 nM was used to exclude subjects with type 1 diabetes. Eye examination was part of a comprehensive physical examination of each participant in the study, as previously described.26 Measurement of IOP was performed with a Goldmann applanation tonometer after instillation of a drop of fluorescein in each eye. Before measuring IOP in the right eye, the tonometer was set to 10 mm Hg, and 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. Measurements were taken in duplicate on each eye, and the average was recorded. None of the subjects was taking glaucoma medication at the time of the measurement. Those found to have elevated IOP were referred for management.
Genotyping was performed at the Center for Inherited Disease Research (CIDR). The CIDR microsatellite marker set is composed primarily of trinucleotide and tetranucleotide repeats, with an average spacing of 8.9 cM and no gaps in the map larger than 18 cM. Average marker heterozygosity was 0.76. Approximately 10% of the marker loci are different between the current CIDR marker set and the Marshfield Genetics screening set version 8. For this study, 372 autosomal short tandem repeat markers were genotyped for an average sex-equal distance of 9 cM, with no gaps greater than 18 cM. The error rate was 0.1% per genotype, and the inconsistency rate was 0.11%. Extensive quality checks were carried out to verify consistency of marker genotyping and stated pedigree relationships, as previously described.27
Descriptive statistics and regressive statistical values were calculated using the Statistical Analysis System (SAS Institute, Cary, NC). The potential confounding influences of age and sex on the distribution of IOP were removed using a multiple regression model. Because IOP displayed a non-Gaussian distribution with a right skew, we performed a logarithmic transformation to obtain a normal (Gaussian) distribution. Transformed data were used in all subsequent analyses. QTL linkage analysis was performed with software (SOLAR, Southwest Foundation for Biomedical Research, San Antonio, TX) through the multipoint variance components approach.28 In variance components linkage analysis, the variance of a trait is decomposed into locus-specific effects determined by the identity-by-descent (IBD) relationships (additive QTL variance), the residual additive genetic effects (additive polygenic variance), and individual specific random environmental effects (random environmental variance).28 29 The null hypothesis is that the additive QTL variance equals zero (no linkage), and this was tested by comparing the likelihood of the restricted model with that of a model in which the additive QTL variance is estimated. The difference between the two log likelihoods produces a logarithm of odds (LOD) score. Twice the difference between the two log likelihoods of these models yields a test statistic that is asymptotically distributed as a 50:50 mixture of a
2 variable and a point mass at zero. One-LOD unit support intervals were obtained by identifying the peaks of the maximum LOD score on the plot of the multipoint results, dropping down 1 LOD unit and finding the chromosomal region defined by the shoulders of the curve.30 An LOD score
3.3 was taken as evidence of significant linkage, and an LOD score
1.9 but <3.3 was taken as evidence of suggestive linkage.31 Marshfield age- and sex-averaged maps were used in the linkage analyses.
To estimate the probability of obtaining false-positive evidence of linkage, we conducted gene-dropping simulations using MERLIN.32 Marker data were simulated under the null hypothesis of no linkage or of association to observed phenotypes while retaining the same pedigree structures, maps, marker allele frequencies, and missing data patterns. We simulated 10,000 replicates and conducted the same linkage analyses described earlier. The probability of obtaining a false-positive result was defined as the proportion of replicates for which we obtained a specified LOD or higher score.
| Results |
|---|
|
|
|---|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
The 5q QTL in this study, the first genome scan for IOP in West Africans, overlaps with a recently described adult-onset POAG locus on 5q22.1,4 designated GLC1G by the HUGO Gene Nomenclature Committee (HGNC). For several reasons, we believe the 5q QTL identified in the present study is a replication of the GLC1G locus. First, the locus identified in the first-pass genome scan in the present study overlaps with the GLC1G locus.4 Second, fine mapping of the locus confirmed and increased the LOD score to greater than 4. Third, the 2-Mb critical interval for GLC1G corresponded to one of our fine-mapping linkage peaks at 116 cM. The present study, in replicating linkage evidence for IOP at a 5q locus for POAG, suggests that this QTL may be important in several geographic populations. On the other hand, a recent study33 mapped a juvenile-onset POAG (JPAG) locus to 5q22.1-q32; this locus did not overlap with the GLC1G minimal interval, and the authors ruled out WDR36 coding sequence mutations in the pedigree studied.
The story of how WDR36 was identified from within the GLC1G locus is instructive. A mutation screen of seven candidate genes from the GLC1G critical region, which spans approximately 2 Mb (D5S1466-D5S2051), identified a significant alteration (D658G) in the WDR36 (WD40-repeat 36) gene.4 Additional screening of WDR36 in 130 POAG families revealed 24 DNA variations with four mutations (N355S, A449T, R529Q, D658G) identified in 17 unrelated POAG subjects but absent in more than 200 healthy control chromosomes.4 These mutations were conserved between WDR36 orthologs in mouse, rat, dog, chimp, and human. WDR36, a novel gene with 23 exons encoding 951 amino acids, is expressed in lens, iris, sclera, ciliary muscles, ciliary body, trabecular meshwork, retina, and optic nerve, as established by RT-PCR.4 However, several lines of evidence suggest that other POAG genes may exist in this region. For example, the region is large and often has several linkage peaks (as in the present study); moreover, previous studies have identified POAG families that map to this region but show no WDR36 mutations. An examination of the annotated human genome sequence (build 35.1) of this region shows other potential candidate genes for POAG apart from the WDR36 gene. These genes include SEMA6A and TGFBI (transforming growth factor ß-induced), which share the "visual perception" gene ontology (GO) annotation with MYOC; DMXL1 and WDR55 (formerly NP_060176), which share the G-protein ß WD-40 repeat Interpro domain with WDR36; and ACSL6, which shares the adenosine monophosphate (AMP)dependent synthetase and ligase Interpro domain with WDR36. However, it should be noted that these are just a few of the potential candidate genes in the region. We are in the process of screening some of these genes in our West African patients for the potential identification of novel variants for subsequent functional analysis.
The 14q22 linkage signal (LOD = 2.5; P = 0.000038) is also noteworthy in that other groups have identified evidence for POAG susceptibility loci in this region.15 25 For example, Wiggs et al.15 reported a maximum LOD score greater than 2.0 for five regions, including chromosome 14. This finding was subsequently supported, though weakly, by the genome scan for POAG in the Barbados Family Study of Open-Angle Glaucoma.25 Interestingly, the 14q22 region is the location for the methylenetetrahydrofolate dehydrogenase gene, MTHFD1, a trifunctional enzyme involved in de novo purine and pyrimidine biosynthesis and in homocysteine metabolism.35 It has been reported that mild hyperhomocysteinemia was associated with POAG and secondary open-angle glaucoma.36 37 38 Furthermore, a polymorphism of another enzyme directly involved in homocysteine metabolism, methylenetetrahydrofolate reductase (MTHFR), showed significant association with hypercysteinemia and POAG.39
Other important glaucoma loci, not observed in the present study, have been reported.5 These loci include 1q23 to 25 (MYOC); 2cen-q13 (GLC1B); 3q21 to 3q24 (GLC1C); 8q23 (GLC1D); 10p15 to 10p14 (OPTN); 7q35 (GLC1F); 9q22 (GLC1); and 20q12 (GLC1).4 16 17 18 19 21 22 23 In the only other study with participants from a population exclusively of recent African descent, Nemesure et al.25 found no support for the myocilin (MYOC) or optineurin (OPTN) genes as causative genes for glaucoma in an Afro-Caribbean population in Barbados.
The present study has some limitations. First, study participants were affected sibling pairs with T2D, but it is uncertain whether the linkage regions tracked IOP, the underlying T2D, or both. However, our genome scan for T2D in this sample27 found no chromosome 5 linkage regions but obtained suggestive linkage to T2D on chromosomes 12 and 20. Therefore, we think it unlikely that the chromosome 5 QTL tracked T2D. Second, we focused on IOP as the phenotype rather than POAG or glaucoma, reasoning that the continuous variable provides more power to detect linkage than a binary variable such as the presence or absence of glaucoma. However, because IOP is an intermediate physiologic phenotype whose elevation does not always correlate with disease (American Academy of Ophthalmology, 2001), it is uncertain whether a study of IOP can be directly compared with one of POAG. Third, we did not measure central corneal thickness, which is a major determinant of IOP. Central corneal thickness itself is highly heritable, and our linkage regions may track (at least in part) this phenotype rather than IOP. Despite these limitations, however, we believe that the data presented here provide an additional important step in the understanding of the genetic susceptibility to glaucoma in various world populations.
Based on our findings and the support provided by other studies, we plan to conduct fine mapping of our linkage signals on chromosomes 5 and 14 with densely placed single nucleotide polymorphisms (SNPs). We are of the opinion that this data set of West Africans with T2DM and associated complications, including IOP, is likely to be highly useful for fine mapping and gene localization given that the region of linkage disequilibrium around a susceptibility allele is likely to be narrower in Africans than in Europeans or Asians.40 Because glaucoma is the leading cause of blindness in African Americans who share an ancestral relationship with West African populations, we anticipate that findings from this project will be directly relevant to the ongoing effort to help reduce the substantial ethnic disparity in glaucoma rates in the United States.
| Footnotes |
|---|
Submitted for publication December 2, 2005; revised March 2 and 29, 2006; accepted June 21, 2006.
Disclosure: C.N. Rotimi, None; G. Chen, None; A.A. Adeyemo, None; L.S. Jones, None; K. Agyenim-Boateng, None; B.A. Eghan, Jr., None; J. Zhou, None; A. Doumatey, None; K. Lashley, None; H. Huang, None; O. Fasanmade, None; F.B. Akinsola, None; F. Ezepue, None; A. Amoah, None; S. Akafo, None; Y. Chen, None; J. Oli, None; T. Johnson, 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: Charles N. Rotimi, National Human Genome Center, Howard University, College of Medicine, 2216 6th Street, NW, Washington, DC 20059; crotimi{at}howard.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] |
||||
![]() |
F. Pasutto, C. Y. Mardin, K. Michels-Rautenstrauss, B. H. F. Weber, H. Sticht, G. Chavarria-Soley, B. Rautenstrauss, F. Kruse, and A. Reis Profiling of WDR36 Missense Variants in German Patients with Glaucoma Invest. Ophthalmol. Vis. Sci., January 1, 2008; 49(1): 270 - 274. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. K. Iyengar The Quest for Genes Causing Complex Traits in Ocular Medicine: Successes, Interpretations, and Challenges Arch Ophthalmol, January 1, 2007; 125(1): 11 - 18. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L. Wiggs Genetic Etiologies of Glaucoma Arch Ophthalmol, January 1, 2007; 125(1): 30 - 37. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |