DOI: 10.2337/db06-1299 © 2007 by the American Diabetes Association Genome-Wide Linkage Analyses to Identify Loci for Diabetic Retinopathy
1 Phoenix Epidemiology and Clinical Research Branch, National Institute of Diabetes and Digestive and Kidney Diseases, Phoenix, Arizona Address correspondence to Robert L. Hanson, PECRB, 1550 E. Indian School Rd., Phoenix, AZ 85014. E-mail: rhanson{at}mail.nih.gov
Abbreviations:
ACIBD, identity by descent; LOD, logarithm of odds; PADI, peptidyl arginine deiminase; WESDR, Wisconsin Epidemiology of Diabetic Retinopathy
Hyperglycemia and long duration of diabetes are widely recognized risk factors for diabetic retinopathy, but inherited susceptibility may also play a role because retinopathy aggregates in families. A genome-wide linkage analysis was conducted in 211 sibships in which 2 siblings had diabetes and retinal photographs were available from a longitudinal study. These sibships were a subset of 322 sibships who had participated in a previous linkage study of diabetes and related traits; they comprised 607 diabetic individuals in 725 sibpairs. Retinal photographs were graded for presence and severity of diabetic retinopathy according to a modification of the Airlie House classification system. The grade for the worse eye was adjusted for age, sex, and diabetes duration and analyzed as a quantitative trait. Heritability of diabetic retinopathy in this group was 18% (95% CI 2–36). A genome-wide linkage analysis using variance components modeling found evidence of linkage on chromosome 1p. Using single-point analysis, the peak logarithm of odds (LOD) was 3.1 for marker D1S3669 (34.2 cM), whereas with multipoint analysis the peak LOD was 2.58 at 35 cM. No other areas of suggestive linkage were found. We propose that an area on chromosome 1 may harbor a gene or genes conferring susceptibility to diabetic retinopathy. Diabetic retinopathy, a common and serious complication of diabetes (1), is one of the leading causes of blindness (2). Many risk factors for diabetic retinopathy have been identified, including poor glycemic control, hypertension (3), albuminuria, and diabetes duration (4–7). Retinopathy is also familial. A study of identical twins found a concordance for diabetic retinopathy in 21 of 31 pairs with type 1 diabetes (68%) and in 35 of 37 pairs with type 2 diabetes (95%) (8). In the Diabetes Control and Complications Trial cohort, the odds ratio (OR) for severe retinopathy when a relative had retinopathy was 3.1 and the pairwise correlation for retinopathy grade was 0.187 for family members (9). In a cohort with type 2 diabetes, the OR for any grade of retinopathy for siblings of subjects with retinopathy was 4.3 (10). More recently, in a study of retinopathy among Mexican-American families, no significant difference was found for the presence of any retinopathy according to the presence of retinopathy in a sibling, but there was a significant OR for the presence of severe retinopathy (11). To date, however, no specific genetic loci influencing susceptibility to diabetic retinopathy have been identified. Candidate gene studies have, for the most part, reported conflicting results (12), and there have been few genome-wide studies. Our group previously conducted genome-wide linkage analysis for retinopathy in diabetic Pima Indians who had participated in a linkage study designed to identify loci for diabetes and obesity (13,14). The results of this analysis showed only weak evidence for linkage on chromosomes 3 and 9 (13). However, the analysis was conducted by assessment of allele sharing among sibling pairs affected by retinopathy, and quantitative information on retinopathy severity was not used. Furthermore, at the time this preliminary analysis was conducted, retinal photographs were only available on a relatively small subset of participants. In the present study, genome-wide linkage analyses were conducted for participants with updated retinal photographs of quantitative retinopathy severity scores.
Data for the present analysis come from an ongoing longitudinal study among the Pima Indians of the Gila River Indian Community in central Arizona. Community residents are invited to undergo examinations at intervals of 2 years from the age of 5 years onward. These examinations include measures of BMI, blood pressure, plasma glucose, and A1C. Albumin-to-creatinine ratio was measured in a spot urine specimen. Diabetes was diagnosed on the basis of a 75-g oral glucose tolerance test according to World Health Organization criteria (15) or a previously documented clinical diagnosis. All participants 15 years old have retinal photographs taken through dilated pupils (two standard 45° fields for each eye) using a Canon CR4–45NM fundus camera, centered at the optic disk and at the macula. When subjects attended multiple examinations over the course of the study, the examination closest to a person having a 10-year duration of diabetes was selected. A genome-wide linkage scan was previously performed on 1,311 individuals selected from the longitudinal study in 332 nuclear families potentially informative for linkage with diabetes and related traits (14). For the current analyses, siblings from the genome-wide linkage group with a diagnosis of diabetes and graded retinal photographs (n = 607) were selected. The retinal photographs were graded in a standardized manner without knowledge of clinical details according to a modification of the Airlie House classification system and a modification of the Early Treatment of Diabetic Retinopathy Study severity system of diabetic retinopathy, which assigns a numeric score to each eye (Table 1) (16). A mean retinopathy score was calculated for each subject [(score in right eye + score in left eye)/2]; when only one eye was gradable, the score for that eye was taken as the mean score. In addition to the mean score, the score for the worse eye was also analyzed. Although these traits were highly correlated (r = 0.98), we analyzed both traits because the worse eye score may represent an individual's susceptibility to retinopathy, whereas the mean score may represent the precise estimate of the current level of retinopathy.
The traits (mean retinal score and worse eye score) for each individual were adjusted for age, sex, and diabetes duration using linear regression. To ensure that the assumption of multivariate normality was met, which is required for the genetic analyses, the residuals from the regression analyses were then normalized with an inverse Gaussian transformation before analysis. Variance components methods were used to calculate heritability (i.e., the proportion of variance potentially due to the additive effects of genetic factors) (17,18). CIs for the heritability were calculated by a systematic exploration of the likelihood surface (19). To further assess the extent of familial aggregation according to severity of retinopathy, the prevalence of retinopathy or varying severity was assessed in siblings of an "index" individual according to the presence or absence of retinopathy in the index sibling (the oldest [i.e., first born] diabetic sibling in each sibship was taken as the index). Thus, ORs of retinopathy in siblings were calculated using logistic regression models for varying levels of retinopathy using the retinopathy level (worse eye grade 21, "any retinopathy"; worse eye grade 31, "moderate nonproliferative diabetic retinopathy or more"; or worse eye grade 60, "proliferative diabetic retinopathy") for the older sibling as the risk factor for a corresponding degree of retinopathy in the younger siblings. These models were adjusted for age, sex, and diabetes duration of the younger sibling and were fit by generalized estimating equations to account for correlation among siblings from the same family.
Genome-wide linkage analysis.
From 211 families, 607 subjects were identified, giving a total of 725 sibpairs. The majority of families (n = 110) contributed one sibpair to the analysis, but up to nine siblings from a single family were included (Table 2). Excluding the larger families ( 7 siblings) from the analyses did not appreciably change the results, so the findings for the whole set are presented. Subject characteristics are shown in Table 3, and the distribution of retinopathy grades are shown in Table 4. Retinopathy grade was strongly correlated with diabetes duration (r = 0.49 for the worse eye score, r = 0.50 for the mean eye score) and with albumin-to-creatine ratio (0.40 for the worse eye score, r = 0.41 for the mean eye score), A1C (r = 0.23 for both worse eye score and mean eye score), diastolic blood pressure (r = 0.22 for the worse eye score, r = 0.23 for the mean eye score), systolic blood pressure (r = 0.17 for the worse eye score, r = 0.18 for the mean eye score), and age (r = 0.13 for the worse eye score, r = 0.14 for the mean eye score). BMI was negatively correlated with retinopathy grade (r = –0.30 for the worse eye score, r = –0.31 for the mean eye score). All correlation coefficients shown are Spearman correlations and P < 0.01 for each.
The heritability of retinopathy in this sample for the worse eye score was 18% (95% CI 2–36), and for the mean eye score, the heritability was 16% (1–33). For each family, the retinopathy grade in the first-born sibling (n = 211, 123 with worse eye grade of 10, 15 with worse eye grade of 21, 61 with worse eye grade of 31–41, and 12 with worse eye grade of 51) was used as a risk factor for younger siblings developing retinopathy. When the first-born sibling had any retinopathy (worse eye score 21), the risk for a younger sibling to have any retinopathy was modestly elevated but was not statistically significant (OR 1.10 [95% CI 0.61–1.97], P = 0.75 adjusted for age, sex, and diabetes duration of the younger sibling and for diabetes duration of the index sibling). The OR was similar for moderate nonproliferative retinopathy or more severe retinopathy (worse eye score 31) (1.11 [0.61–2], P = 0.73), but there was a greater risk for proliferative retinopathy (worse eye score 60) (2.46 [0.85–7.15], P = 0.10). All points where the LOD score exceeded 1 (locus-specific P < 0.03) are shown in Table 5 for the single-point linkage analyses. The strongest evidence for linkage was on chromosome 1 with marker D1S3699 for both the worse eye score and the mean eye score. Results of single-point and multipoint analyses for all chromosomes are shown in Fig. 1 for the worse eye phenotype. Findings were similar for the mean eye score phenotype. In the multipoint analysis, the peak LOD score was on chromosome 1 at 35 cM for both the worse eye and mean eye score phenotypes. The peak multipoint LOD was 2.58 for the worse eye and 1.94 for the mean eye score. No other chromosomal regions showed evidence of linkage for retinopathy with LOD >1.
Both twin studies (8) and familial aggregation studies (9,10) show evidence for genetic susceptibility to diabetic retinopathy. The heritability estimate in the present study shows a modest level of familial aggregation in this population that was significantly different from 0, suggesting a genetic role in the development and/or progression of retinopathy. The majority of studies on the genetics of diabetic retinopathy have focused on a candidate gene approach. The results of such studies have thus far been somewhat inconsistent, which may reflect, in part, the differing definitions for the phenotype used and the different populations included in the studies. The present genome-wide linkage analysis used variance components modeling to identify the areas of linkage for diabetic retinopathy. This method uses quantitative information about the degree of retinopathy from all siblings. The results showed suggestive evidence for linkage to chromosome 1p36. Genes that have shown evidence of association with any degree of retinopathy among candidate gene studies include inducible nitric oxide synthase gene (25), intracellular adhesion molecule-1 gene (26), and aldose reductase gene (27). The transforming growth factor-ß gene has shown evidence of association only with proliferative retinopathy (28). A previous candidate gene study in the Pimas reported a modest association between retinopathy (as assessed by fundoscopy) and a polymorphism in the promoter region of the plasminogen activator inhibitor 1 gene on chromosome 7 (29). None of these genes, however, map to the area of linkage identified in the current study. In linkage studies, it is customary to require particularly stringent levels of statistical significance. A LOD score >3 (locus-specific P < 0.0001) is conventionally considered statistically significant linkage, whereas a LOD score >2 (locus-specific P < 0.0012) is considered suggestive but not definitive evidence for linkage (30). A previous affected sibpair analysis in a subset of the individuals in the present study showed tentative linkage on chromosomes 3 and 9 (13); however, these regions did not show linkage in the current study. The previous study included a fairly small number of affected sibling pairs (n = 103, comprising 156 individuals, all of whom were included in the present study), and the LOD scores did not achieve levels that were suggestive of linkage by conventional criteria. The present analysis is presumably more powerful because it includes data on a larger number (n = 607) of individuals and because the LOD score on chromosome 1p is considered suggestive of linkage. The present analyses also incorporated information about severity of retinopathy based on fundus photographs, the more severe levels of which tend to be more familial than simply the presence of any retinopathy, which was the phenotype analyzed in the previous study. This is in keeping with results from Mexican-American families in Texas, where higher ORs were found for more severe retinopathy within siblings than for milder disease (11). It is possible that the differences between the two analyses are explained by interactions with other factors (such as duration of diabetes) or, given the relatively small sample size, by chance. Ultimately, replication studies in other populations will be required to more firmly establish susceptibility loci for diabetic retinopathy. For variance components modeling, a normal distribution of the phenotypic measure is assumed, although the method is robust to moderate violations of normality (17,18). Although retinopathy was treated as a continuous variable in this study, the distribution of retinopathy scores was highly skewed, because the majority of subjects had no retinopathy at their examination. Because of the severely skewed distribution of retinopathy scores, data were normalized (after adjustment for age, sex, and duration of diabetes) before linkage analysis. The normalization procedure ensures that the assumptions of normality are met, although it may compromise the power somewhat compared with analysis without normalization. The procedure allows for adjustment for important potential confounders such as duration of diabetes. The use of variance components analysis assumes that the same gene(s) are responsible for retinopathy in general. However, it is quite plausible that there are genes that cause a susceptibility to proliferative retinopathy, for instance, that are not involved in the susceptibility to nonproliferative retinopathy. The OR analysis certainly suggests a stronger familial risk for severe retinopathy compared with any grade of retinopathy. Because retinopathy is also highly dependent on diabetes duration, subjects studied at short durations without retinopathy may still have the genes for retinopathy susceptibility. An affected sibpair analysis might be a preferable method if this is the case rather than adjustment of the retinopathy grade for diabetes duration as is done here. However, we do not present data here for affected sibpair analysis because the number of affected (any retinopathy) sibpairs was low (n = 98), and thus the analysis lacked adequate power.
To assess how diabetic retinopathy in the Pima Indians is representative of diabetic retinopathy in the wider population, we can look at previous reports of prevalence, incidence, and progression. Estimates of the prevalence of diabetic retinopathy vary between studies. A recent analysis of pooled studies estimated a crude prevalence of retinopathy of 40.3% in adults aged
The incidence of diabetic retinopathy is strongly related to diabetes duration (4,5,37), and, therefore, the prevalence is higher when subjects are examined after longer durations of diabetes. The longitudinal study design allowed for subjects to be examined at varying durations of diabetes. Confounding by duration was minimized by selecting examinations nearest to 10-year duration, but a wide range of durations were included because many participants only had retinal photographs at other durations of diabetes. Differences in duration were accounted for by adjusting for duration by linear regression. It is possible that genetic factors have a stronger influence on early predisposition to retinopathy or the development of more severe retinopathy, but there were insufficient members with early-onset severe retinopathy to assess this in the present study. Because only The peak of linkage on the short arm of chromosome 1 has not been studied previously for genes associated with diabetic retinopathy. Although the 1p36 region is rich in genes, many have not yet been fully characterized. Figure 2 shows a diagram of identified genes that fall between the markers D1S228 and D1S522, which defined our linkage peak in the single-point analysis with the peak of linkage found at marker D1S3669 (diagram created using data from the National Center for Biotechnology Information Web site at http://www.ncbi.nlm.nih.gov as of August 2006).
Of the identified genes in the area, those coding for the peptidyl arginine deiminases ([PADI] 1, 2, 3, 4, and 6; Fig. 2, D–H), CASP-9 (Fig. 2, h), CLCN-Ka, and CLCN-Kb (Fig. 2, r and s) are close to our peak of linkage and offer possible mechanisms for diabetic retinopathy. PADI converts arginine residues to citrulline residues. PADI 2 is expressed in the brain, is found in retinal tissue (38), and may play a role in the pathogenesis of multiple sclerosis (39). PADI 4 is expressed primarily in neutrophils and eosinophils (40) and is associated with rheumatoid arthritis (41). Recently, PADI 4 was reported to regulate histone argine methylation and thus could potentially play a broader role in gene expression (42). CASP-9 codes for caspase-9, which is activated by the release of cytochrome c from the mitochondria and itself activates other caspases to induce cell apoptosis (43). In vitro, advanced glycosylation end products induce apoptosis in bovine retinal pericytes, and this effect was inhibited by the use of caspase inhibitors, including one specific to caspase-9 (44). The CLCN-Ka and CLCN-Kb genes code for a kidney-specific chloride channel. A functional mutation in the CLCN-Kb is associated with hypertension (45). However, in this study, the correlation between retinopathy score and blood pressure is modest. This study suggests that diabetic retinopathy is modestly heritable in Pima Indians with type 2 diabetes. It also demonstrates that an area of chromosome 1, which has not previously been studied in relation to diabetic retinopathy, has suggestive evidence of linkage. Further linkage studies of diabetic retinopathy may help to define the role of this locus in other populations. Association studies of candidate genes in the region, such as CASP-9, PADI 4, and CLCN-Kb, may also be informative.
This research has received support from the Intramural Research Program of the National Institute of Diabetes and Digestive and Kidney Diseases and by the National Eye Institute. We thank Susan Vitale at the National Eye Institute for her help during the planning of this project and Leslie Baier for her advice in the preparation of this manuscript. We are grateful to the participants of the Gila River Indian Community for their participation in this study, to all at the National Institutes of Health clinic (Sacaton, AZ), and to the staff at the Ocular Epidemiology Reading Center (Madison, WI).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Received for publication September 15, 2006 and accepted in revised form December 26, 2006
This article has been cited by other articles:
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||