Ordered subset linkage analysis supports a susceptibility locus for age-related macular degeneration on chromosome 16p12
BMC Genetics volume 5, Article number: 18 (2004)
Age-related macular degeneration (AMD) is a complex disorder that is responsible for the majority of central vision loss in older adults living in developed countries. Phenotypic and genetic heterogeneity complicate the analysis of genome-wide scans for AMD susceptibility loci. The ordered subset analysis (OSA) method is an approach for reducing heterogeneity, increasing statistical power for detecting linkage, and helping to define the most informative data set for follow-up analysis. OSA assesses the linkage evidence in subsets of potentially more homogeneous families by rank-ordering family-specific lod scores with respect to trait-associated covariates or phenotypic features. Here, we present results of incorporating five continuous covariates into our genome-wide linkage analysis of 389 microsatellite markers in 62 multiplex families: Body mass index (BMI), systolic (SBP) and diastolic (DBP) blood pressure, intraocular pressure (IOP), and pack-years of cigarette smoking. Chromosome-wide significance of increases in nonparametric multipoint lod scores in covariate-defined subsets relative to the overall sample was assessed by permutation.
Using a correction for testing multiple covariates, statistically significant lod score increases were observed for two chromosomal regions: 14q13 with a lod score of 3.2 in 28 families with average IOP ≤ 15.5 (p = 0.002), and 6q14 with a lod score of 1.6 in eight families with average BMI ≥ 30.1 (p = 0.0004). On chromosome 16p12, nominally significant lod score increases (p ≤ 0.05), up to a lod score of 2.9 in 32 families, were observed with several covariate orderings. While less significant, this was the only region where linkage evidence was associated with multiple clinically meaningful covariates and the only nominally significant finding when analysis was restricted to advanced forms of AMD. Families with linkage to 16p12 had higher averages of SBP, IOP and BMI and were primarily affected with neovascular AMD. For all three regions, linkage signals at or very near the peak marker have previously been reported.
Our results suggest that a susceptibility gene on chromosome 16p12 may predispose to AMD, particularly to the neovascular form, and that further research into the previously suggested association of neovascular AMD and systemic hypertension is warranted.
Age-related macular degeneration (AMD) affects the central region of the retina (macula), which has the highest concentration of cone photoreceptors and is responsible for central visual acuity. In developed nations, it is the most common cause of irreversible blindness in older adults. Approximately 4% of individuals over 60 years of age and 10% of those over 75 years of age have advanced stages of the disorder, which include geographic atrophy (dry AMD) and neovascular (wet) AMD . While both forms lead to loss of central vision, it is more common and occurs more rapidly with neovascular AMD. It is unknown whether the two clinical subtypes have a distinct etiology, but longitudinal studies have shown that the presence of large soft drusen, which are extensive extracellular protein/lipid deposits in the macula, increases the risk of progressing to either form of advanced AMD [2–4].
AMD has a complex etiology likely to result from the interplay of several risk factors, both genetic and environmental. A contribution of genetic susceptibility is supported by epidemiologic [5–7] and twin studies [8, 9], as well as segregation analyses . Candidate gene association studies have examined many genes responsible for retinal disorders with Mendelian inheritance, with generally negative results. One of the most intensely studied genes has been the ABCA4 gene on chromosome 1p, which causes juvenile-onset autosomal-recessive Stargardt disease. An initial study reporting an increased risk of AMD for carriers of two particular ABCA4 sequence variants  was followed by only one positive replication study  and multiple studies reporting an absence of this association [13–19]. At this point, the ABCA4 gene is not believed to be a major susceptibility gene for AMD, although it may account for a small proportion of the disease, possibly only the dry form. A candidate gene that has consistently been reported to be associated with AMD in multiple independent studies is the apolipoprotein E (APOE) gene on chromosome 19q. The APOE-4 allele, or a nearby allele in linkage disequilibrium with APOE-4, appears to be protective for AMD [20–28]. To date, only two studies failed to confirm this finding in a Caucasian  and a Chinese sample . The APOE-2 allele may increase AMD risk in smokers .
Several research groups have collected multiplex families (2+ sampled family members with AMD) to perform a genome-wide screen for AMD susceptibility loci [32–37]. Despite variable phenotype definitions and different analysis approaches, these genome scans identified remarkably consistent regions of linkage on several chromosomes, including chromosome 1q25-31, 10q25-26, 12q21-23 and 16p11-12. None of the genes responsible for these linkage signals have thus far been identified. The hemicentin-1 (FIBL6) gene, located on 1q31, has been proposed as a rare cause of AMD , but awaits confirmation by other research groups.
We recently genotyped 62 multiplex AMD families ascertained through Duke University Medical Center (DUMC) and Vanderbilt University Medical Center (VUMC) for 389 microsatellite markers distributed at 10 cM density across the human genome. These families were screened for the first time and were not included in previously published genome screens for AMD [32–37]. For all individuals enrolled in our study, an extensive array of clinical, anthropometric, demographic and environmental covariates were collected. The ordered subset analysis method  is one approach for incorporating such covariate information into nonparametric linkage analysis. The goal of the method is to test whether the evidence for linkage is significantly influenced by a trait-related covariate, which may define a genetically more homogeneous subset of families. In addition to family history and increasing age, there are several well-established risk factors for AMD. Smoking is considered a major modifiable risk factor and appears to increase the risk of both the atrophic and neovascular disease type . Systemic hypertension is another risk factor, particularly for the neovascular form of AMD , and is potentially associated with increased intraocular pressure in some racial groups [42, 43]. Increased systolic blood pressure (SBP) was shown to be a significant predictor of AMD incidence in two large prospective studies [44, 45]. Other factors associated with an increased risk of hypertension can be considered indirect risk factors for AMD, but some of them have also been implicated as independent predictors of risk, such as obesity . A recent study reported that overall obesity (measured as body mass index, BMI) and abdominal obesity (measured by waist-to-hip ratio and waist circumference) were the most significant variables associated with an increased risk of progressing from early to advanced stages of AMD . On the basis of these reported clinical associations, this study focused on incorporating the above covariates into a nonparametric linkage analysis of our multiplex families to reduce the phenotypic and genetic heterogeneity of AMD and potentially improve our ability to detect linkage.
OSA analysis of multiplex families with early or advanced AMD
Descriptive characteristics of the study population are shown in Table 1. OSA-defined family subsets with lod score increases meeting uncorrected statistical significance (p ≤ 0.05) are shown in Table 2. When a corrected significance level of 0.006 was used, only regions on chromosome 14q13 and 6q14 produced significant increases in lod scores. When ordering families by ascending average IOP, a maximum lod score of 3.2 was obtained between markers D14S608 and D14S599 (33 cM) in 28 families with average IOP ≤ 15.5. This was a significant increase from the baseline lod score for all families (0.31; p = 0.002). When ordering families by descending average BMI, a maximum lod score of 1.6 in eight families with average BMI ≥ 30.1 was obtained at marker D6S1031 (90 cM; p = 0.0004). While the third region of interest, chromosome 16p12, did not reach the corrected statistical significance level of 0.006, it met all of the other criteria for judging the clinical relevance of OSA results mentioned in the Methods section: Consistency across several clinically meaningful covariates, agreement with previously reported linkage signals, and continued significance when using a more stringent phenotype definition. A maximum lod score of 2.2 at marker D16S403 (44 cM) was obtained for a subset of six families with average SBP ≥ 153.5 (p = 0.04). Ordering families by descending average IOP produced a nominally significant lod score increase in a subset of 35 families with average IOP ≥ 15.4 very close to this map position (39 cM, p = 0.04). The peak marker (D16S403) is located only 7 cM distal to D16S769, which was reported as a significant marker (p = 0.009) in a genome-wide analysis of AMD as a quantitative trait in 102 nuclear families from Beaver Dam, Wisconsin , and was also significant (p = 0.005) in a second sample of 34 extended pedigrees ascertained from the same geographic region . In both studies, the Wisconsin Age-Related Maculopathy Grading System, which is a 16-level severity scale for the extent of macular damage, was used as the trait in a Haseman-Elston regression applied to all genotyped sibling pairs .
OSA analysis of multiplex families with advanced AMD
To implement a more stringent binary phenotype definition for further evaluation of our OSA results, we excluded individuals with early AMD (grade 3) from the set of affected individuals, which reduced the number of families with at least one sampled affected relative pair to 45. When the OSA procedure was repeated on the chromosomes where lod score increases met nominal statistical significance in the initial analysis (Table 2), only the result for chromosome 16p12 maintained nominal significance. The maximum lod score for the IOP-defined subset of 32 families increased to 2.9 at 39 cM (p = 0.008). The maximum lod score for the SBP-defined subset was 2.3 near the same map position (44 cM, p = 0.05), and the number of families in this subset increased from 6 in the previous analysis to 25 (with average SBP ≥ 133). In addition, a nominally significant lod score increase at 44 cM was observed when families were rank-ordered by average BMI (9 families with average BMI ≥ 28.6 and lod score 2.0 at D16S403, p = 0.04). Eight of the nine families in the BMI subset and 20 of the 25 families in the SBP subset were part of the IOP subset. Multipoint lod score curves for the BMI, SBP and IOP subsets with significantly increased linkage evidence to chromosome 16p12, along with the baseline lod score for the data set, are displayed in Figure 1. Similar plots for chromosomes 14 and 6 are shown in Figures 2 and 3.
Clinical features of families identified by OSA
The clinical features of family subsets identified by OSA are shown in Table 3. The eight families in the BMI-defined subset linked to chromosome 6q14 were characterized by higher mean BMI (31.4) and pack-years (PKYRS) of cigarette smoking (24.5), compared to the entire data set of 62 multiplex families (averages of 26.2 and 21.1, respectively). They also included a higher proportion of early AMD (33.3%, versus 19.7% in the entire data set). The nine families in the BMI-defined subset linked to 16p12 were also characterized by higher mean BMI (29.3) and PKYRS (25.5), compared with the entire data set of 45 multiplex families (averages of 25.7 and 22.5, respectively). Only three families were in the BMI-defined subset contributing to both linkage signals on 6q14 and 16p12. Due to the opposite rank ordering of the IOP covariate in OSA (lowest to highest for 14q13 and highest to lowest for 16p12), the 28 families with linkage to 14q13 and the 32 families with linkage to 16p12 had IOP averages below (13.9) and above (17.3) the population average of about 15.0, as expected. Seven families with IOP averages between 15 and 15.5 contributed to both linkage signals. Finally, the 25 families in the SBP-defined subset with linkage to 16p12 had higher mean SBP (145.0) than the entire data set of 45 multiplex families (137.8). As mentioned above, most of these families (20/25) were also part of the IOP-defined subset. In all three covariate-based subsets with linkage to 16p12, as well as the baseline data set of 45 multiplex families with advanced AMD, the proportion of neovascular AMD was much higher than the proportion of geographic atrophy (71–83%). While there is some overlap in families included in the various covariate-based subsets, the five variables used for OSA showed limited correlation in the AMD patient data set included in this analysis. Only the Pearson correlations of IOP and BMI (r = 0.20, p = 0.03), and, as expected, SBP and DBP (r = 0.51, p < 0.0001) were statistically significant. BMI was inversely correlated with age at exam (r = -0.26, p = 0.004), while IOP (r = -0.12, p = 0.18), SBP (r = 0.19, p = 0.06), DBP (r = -0.03, p = 0.77) and PKYRS (r = -0.11, p = 0.25) were not significantly correlated with age.
Our analysis supports the possibility of distinct AMD susceptibility loci on three chromosomal regions: 6q14, with a lod score of 1.6 in a subset of eight overweight families; 14q13, with a lod score of 3.2 in 28 families with lower-than-average IOP values; and 16p12, with a lod score of 2.9 in 32 families with higher-than-average IOP values, most of which also had above-average values of SBP and BMI. Of these three regions, we believe that the 16p12 linkage is the most interesting finding since there is consistency across several clinically meaningful covariates, agreement with prior studies, and increased statistical significance with a more stringent phenotype definition. Based on our results, a gene on 16p12 may be associated with an increased risk of primarily neovascular AMD. The peak marker, D16S403, is located very close to marker D16S769, which was implicated in two prior genome screens of the Beaver Dam Study based on different phenotype definitions and statistical analysis methods (p = 0.009 and p = 0.005, respectively). However, the families included in these prior screens were not described with respect to the vascular risk factors we have evaluated here. It is possible that the significant result for the 16p12 marker in the Beaver Dam Study is primarily due to a contrast in identity-by-descent sharing of individuals with neovascular versus other forms of AMD, although a more detailed analysis would be necessary to confirm this hypothesis.
The BMI-defined subgroup with linkage to 16p12 fits the definition of being overweight. Increased BMI has previously been associated with a greater risk of AMD [48, 49]. It was proposed that overweight and obese individuals may have lower macular pigment optical density (MPOD), which is a measure of retinal levels of the carotenoids lutein and zeaxanthin . Carotenoids may be protective for AMD, implying that lower MPOD may confer a greater risk of AMD . Our results suggest that individuals in the BMI-defined subset with linkage to 16p12 were also heavier smokers, as indicated by a greater average of pack-years of smoking. Epidemiologic studies have consistently reported an increased risk of both dry and wet AMD due to smoking , but little is known about the underlying mechanism. A recent study reported that nicotine increased the size and severity of experimental choroidal neovascularization in a mouse model of wet AMD and suggested that nicotinic receptor activation may mediate this harmful effect .
The SBP-defined subset of families meets the definition of systemic hypertension (SBP ≥ 140 mmHG or DBP ≥ 90 mmHG according to the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure, ). Elevated SBP is a more important predictor of cardiovascular risk than elevated DBP in persons older than 50 years. Consistent with our results, the Beaver Dam Study reported that higher SBP at baseline and longitudinal increase of SBP were significantly associated with the 10-year incidence of neovascular, but not atrophic AMD . The Rotterdam Study reported an increased incidence of early AMD with higher baseline SBP . The shorter follow-up time precluded the assessment of whether SBP specifically influenced the progression of early to neovascular, rather than atrophic, AMD. A large case-control study also reported a positive association of hypertension and neovascular, but not atrophic AMD . The 16p12 region harbors genes for monogenic forms of hypertension (SCNN1B [MIM 600760], SCNN1G [MIM 600761]). However, until the region responsible for the linkage signal in our data set can be substantially narrowed down by association mapping, we believe it would be premature to consider these particular genes as promising locational candidate genes for neovascular AMD.
The IOP-defined subgroup with linkage to 16p12 is more difficult to interpret, since there is little evidence to suggest that elevated IOP is a risk factor for AMD. In African-American samples, cross-sectional associations between SBP and IOP as well as positive associations of elevated SBP and DBP at baseline with longitudinal IOP increases have been reported [42, 43]. However, the general relationship of systemic hypertension and IOP is not well documented, and in the data set used here, IOP was only significantly correlated with BMI. It is unknown whether the mechanisms that contribute to hypertension or increased IOP have anything in common with the angiogenic processes observed in neovascular AMD. It is conceivable that both IOP and SBP values in the subgroup of families linked to 16p12 may be correlated with unmeasured markers of retinal vascular or choroidal blood flow characteristics. Correlations between elevated systemic blood pressure and retinal microvascular changes have been reported in multiple large-scale studies [54–56]. Retinal arteriolar narrowing in particular has been associated with increased risk of coronary heart disease  and with AMD progression in the Blue Mountains Eye Study . Due to the limited size of our study, further interpretation of the SBP- and IOP-related results from this analysis awaits replication in a larger study population.
A linkage signal on 14q13 near the peak markers in our data set (D14S608, D14S599) was previously reported by one other study with a two-point parametric lod score of 1.5 . While we observed the highest subset-based lod score in our data set (3.2) in this region, it is not clear whether the particular subset of 28 families with lower-than-average IOP values is defined by a clinically meaningful covariate. As mentioned above, little is known about the relationship of IOP and AMD risk. A statistically significant heritability of IOP was recently reported by the Beaver Dam Study , which gives some support to the use of IOP as a family covariate. In general, IOP is known to be influenced by many factors, including age, sex, refractive error, serum cholesterol, SBP, DBP, and BMI , although in our data set, IOP was only significantly correlated with BMI. The IOP-defined subset with linkage to 14q13 may represent individuals with a distinct profile of AMD risk factors that were not captured by the covariates considered here, but could be correlated with lower IOP. Alternatively, this result may be due to type I error, despite the correction for multiple testing we applied.
A linkage signal on 6q14 for the same peak marker as in our data set (D6S1031) was reported by two prior genome screen studies, with a p-value of 0.04 in Haseman-Elston regression  and an NPLpairs p-value of 0.00001 . However, since the number of families contributing to the linkage evidence in our data set is very small, it is difficult to speculate about the plausibility of our finding. The observation of a higher proportion of early AMD in the BMI-defined subset linked to 6q14, relative to the overall data set, is consistent with previous reports of an increased risk of early, but not late, AMD for both underweight and overweight individuals  and of a specific association between RPE abnormalities and BMI . Alternatively, the lod score in this region may reflect identity-by-descent sharing of a gene that predisposes to obesity, rather than to AMD itself. One of the most consistently reported regions linked to BMI as a quantitative trait is chromosome 6q23-25 , however, this region is located at least 50 cM away from our region. It has been suggested that the incorporation of continuous covariates into binary trait linkage analysis may identify linkage signals that are distinct from those detected by analyzing such covariates as quantitative traits [64, 65], but more methodological research is needed to further explore this question.
Our results, particularly those for chromosome 16p12, illustrate the utility of OSA for incorporating continuous covariates into linkage analysis of complex traits. The method is a conceptually and computationally simple approach to evaluating linkage evidence in subsets of families that are more homogeneous with respect to clinical features and/or non-genetic risk factors for the disease under study [64, 66].
Phenotypically more similar families may be genetically more homogeneous as well, in which case OSA can greatly improve the power of linkage analysis. In contrast to the admixture test for parametric lod score analysis implemented in the program HOMOG , which simply allows for a proportion of families to be unlinked to a particular region under study, OSA may provide insight into the reasons for the underlying heterogeneity. It can have better power than the admixture test to detect linkage in subsets of families when the overall genetic effect is low or when the families are small, while the admixture test tends to be more powerful when the families are larger and provide more variability in family-specific lod scores . It may also provide better localization of the putative susceptibility gene. The reduction of genetic heterogeneity is especially important in the analysis of complex disorders, for which AMD is a prime example. Limitations of OSA include the inability to incorporate more than one covariate at a time. The power of the method largely depends on the degree of correlation between the evidence for linkage and the levels of the OSA covariate, i.e., on the extent to which phenotypic heterogeneity between families, as measured by averaged covariate values of affected family members, reflects underlying genetic heterogeneity .
In summary, our data support the presence of an AMD susceptibility locus on chromosome 16p12 that may predispose primarily to neovascular AMD. While it would be premature to speculate about biological relationships between SBP, IOP, and genetic predisposition in AMD etiology, our findings suggest that further research into the previously suggested association of neovascular AMD and systemic hypertension is warranted. The question of whether wet and dry AMD may have a different pathogenesis has long been debated. Distinct risk factor associations for the two disease forms were recently reported by the Beaver Dam Study . At the genetic level, our data provide some support of a distinct etiology. Larger family data sets with sufficiently large proportions of both dry and wet AMD, particularly those for which baseline linkage results already exist, should be characterized with respect to the clinical variables considered here to further investigate the possibly distinct genetic basis of the two advanced stages of AMD.
Multiplex AMD families were ascertained at DUMC and VUMC via a proband with early or advanced AMD, as previously described . All individuals included in this analysis were white, and their demographic and clinical characteristics are shown in Table 1. The assignment of AMD affection status was based on the clinical evaluation of stereoscopic color fundus photographs of the macula (EAP, AA, MADLP), according to a system described previously [5, 68]. This system is a slight modification of the Age-Related Eye Disease Study (AREDS) grading system , using example slides from the Wisconsin Grading System  and the International Classification System  as guides. Briefly, individuals were assigned an AMD grade ranging from 1 through 5 based on the macular characteristics found within a 3000 μm-radius centered on the fovea. Eyes with extensive (≥ 15) small drusen (< 63 μm), non-extensive intermediate (≥ 63 μm) drusen or pigment abnormalities were assigned grade 2; eyes with extensive intermediate or any large (≥ 125 μm) drusen, with or without drusenoid (non-fluid) RPE detachments, were assigned grade 3; eyes with geographic atrophy were assigned grade 4, and eyes with serous or hemorrhagic RPE detachments, or choroidal neovascular membrane, were assigned grade 5. Eyes without any drusen and pigment abnormalities, or only small non-extensive drusen, were assigned grade 1.
For the purposes of performing a genome-wide screen for AMD susceptibility loci, multiplex families were defined as those with at least two sampled first- or second-degree relatives with grade 3 (early AMD/ARM), grade 4 (atrophic AMD), or grade 5 (neovascular AMD) in at least one eye. The resulting data set included 62 families, 147 affected individuals, 38 additional family members, and a total of 119 affected sibling pairs. Sixteen families had three or more affected siblings. Only three families had affected relative (avuncular) pairs other than sibling pairs. The number of affected siblings ranged from 2 to 5, with an average of 2.4.
Laboratory and statistical analysis
A total of 389 microsatellite markers spaced at an average 10 cM density across the human genome were genotyped on 185 individuals by the Center for Inherited Disease Research (CIDR). Prior to removal of genotypes that were inconsistent with Mendelian inheritance, pedigree relationships were verified with the programs RELPAIR  and PREST , both of which use multipoint identity-by-descent sharing estimates to infer the most likely relationship between pairs of individuals in the data set. No misspecified relationships were detected. For analysis, inconsistent genotypes detected by the program PEDCHECK  were removed. Inter-marker distances and marker order were obtained from the genetic linkage maps developed by the Marshfield Medical Research Foundation .
Motivated by successful applications of the ordered subset analysis (OSA) method in studies of other complex disorders [64, 66], the primary goal of the analyses presented here was to incorporate AMD-associated clinical covariates into the nonparametric linkage analysis of the genome screen data. Given the extensive phenotypic heterogeneity of AMD, this may help identify more homogeneous subsets of families for follow-up analysis and has the potential to replicate previously published linkage signals that could be obscured in an analysis of the entire data set. To this end, the OSA method  was applied as follows: First, families were rank-ordered by the average covariate value of affected family members. Family-specific multipoint lod scores, which can in principle be parametric or nonparametric, were added one at a time in the covariate-based rank order, at each position on the chromosome map. Since the vast majority (59 of 62, 95%) of our multiplex families were nuclear families with two or more affected siblings, but no other affected relative pairs, we used the nonparametric MLS method for affected sibling pair data  to compute family-specific lod scores. This method has been implemented in the program SIBLINK . Families with n affected siblings were weighted by a factor of n-1, and an additive model was assumed. For each ordered subset of families, the maximum lod score anywhere on the chromosome was determined, the next family was added, and the procedure was repeated until all families had been analyzed in this way. The maximum subset-based lod score for each covariate ordering, along with the map position at which it occurred, was obtained. To evaluate whether the covariate-based subset of families provided significantly increased evidence of linkage, the observed maximum OSA lod score was compared to an empirical distribution of lod scores. This distribution was generated by randomly permuting the order in which families with available covariate information were added and computing the maximum lod score for each permutation as described above. The empirical p-value thus computed indicates how likely it is to obtain a subset-based lod score of the same or greater size than the observed OSA maximum lod score. It corresponds to a test of the null hypothesis of no increase in linkage evidence by rank-ordering families with respect to their covariate values. For the results presented here, we used a minimum of 10,000 permutations to compute empirical p-values.
We applied the OSA procedure with the following covariates: Body mass index (BMI), defined as self-reported weight (in kilograms) divided by squared height (in meters); systolic (SBP) and diastolic (DBP) blood pressure, defined as the average of two sequential measurements taken with the Hawksley random zero sphygmomanometer at the time of the clinical exam; intraocular pressure (IOP), measured by Goldmann applanation tonometry; and pack-years of self-reported cigarette smoking (PKYRS). To compute this combined measure of duration and dosage of cigarette smoking prior to study enrollment, we asked study participants (i) whether they had ever smoked cigarettes at least once per week, (ii) at which age they started and, if applicable, stopped smoking, and (iii) how many cigarettes, on average, they smoked per day. From this information, pack-years of cigarettes were computed as the product of smoking duration (in years), relative to a reference age 10 years prior to study enrollment, and dosage (number of cigarettes per day divided by 20). For never-smokers, zero pack-years were used. For BMI, IOP and PKYRS, we applied OSA with two independent ranking orders, lowest to highest and highest to lowest. For the two blood pressure variables, it was not clear how to best correct for the potential use of anti-hypertensive medication. We felt that only the presence of higher blood pressure at the time of the clinical exam indicated possible hypertension, while lower pressures may be observed for true normotensives as well as hypertensives on blood pressure-lowering medication. Therefore, we only applied the highest to lowest covariate ordering for SBP and DBP in OSA. Thus, eight maximum nonparametric lod scores, corresponding to two covariates with one ranking order and three covariates with two ranking orders, were obtained for each of the 22 autosomes, and a corrected chromosome-wide significance level of 0.006 (= 0.05/8) was used. Results with p-values ≤ 0.0003 (= 0.006/22) may be considered as having genome-wide significance. In addition to correcting for multiple testing, we applied the following considerations to help identify the most promising results of our OSA analysis: (i) Consistency of results across multiple clinically plausible covariates; (ii) agreement with prior published genome screens of AMD; (iii) increase or persistence in statistical significance when using a more stringent phenotype definition. Clinical features for family subsets with a statistically significant increase in nonparametric lod score were analyzed with the Statistical Analysis System (SAS Institute, Cary, NC, version 8).
Klein R, Klein BE, Linton KL: Prevalence of age-related maculopathy. The Beaver Dam Eye Study. Ophthalmology. 1992, 99: 933-943.
Seddon JM, Cote J, Davis N, Rosner B: Progression of age-related macular degeneration: association with body mass index, waist circumference, and waist-hip ratio. Arch Ophthalmol. 2003, 121: 785-792. 10.1001/archopht.121.6.785.
Klaver CCW, Assink JJM, van Leeuwen R, Wolfs RCW, Vingerling JR, Stijnen T, Hofman A, De Jong PTVM: Incidence and progression rates of age-related maculopathy: The Rotterdam study. Invest Ophthalmol Vis Sci. 2001, 42: 2237-2241.
van Leeuwen R, Klaver CC, Vingerling JR, Hofman A, de Jong PT: The risk and natural course of age-related maculopathy: follow-up at 6 1/2 years in the Rotterdam study. Arch Ophthalmol. 2003, 121: 519-526. 10.1001/archopht.121.4.519.
Seddon JM, Ajani UA, Mitchell BD: Familial aggregation of age-related maculopathy. Am J Ophthalmol. 1997, 123: 199-206.
Klaver CC, Wolfs RC, Assink JJ, Van Duijn CM, Hofman A, de Jong PT: Genetic risk of age-related maculopathy. Population-based familial aggregation study. Arch Ophthalmol. 1998, 116: 1646-1651.
Klein BE, Klein R, Lee KE, Moore EL, Danforth L: Risk of incident age-related eye diseases in people with an affected sibling : The Beaver Dam Eye Study. Am J Epidemiol. 2001, 154: 207-211. 10.1093/aje/154.3.207.
Meyers SM, Greene T, Gutman FA: A twin study of age-related macular degeneration. Am J Ophthalmol. 1995, 120: 757-766.
Hammond CJ, Webster AR, Snieder H, Bird AC, Gilbert CE, Spector TD: Genetic influence on early age-related maculopathy: a twin study. Ophthalmology. 2002, 109: 730-736. 10.1016/S0161-6420(01)01049-1.
Heiba IM, Elston RC, Klein BE, Klein R: Sibling correlations and segregation analysis of age-related maculopathy: the Beaver Dam Eye Study [published erratum appears in Genet Epidemiol 1994;11(6):571]. Genet Epidemiol. 1994, 11: 51-67.
Allikmets R, Shroyer NF, Singh N, Seddon JM, Lewis RA, Bernstein PS, Peiffer A, Zabriskie NA, Li Y, Hutchinson A, Dean M, Lupski JR, Leppert M: Mutation of the Stargardt disease gene (ABCR) in age-related macular degeneration [see comments]. Science. 1997, 277: 1805-1807. 10.1126/science.277.5333.1805.
Allikmets R: Further evidence for an association of ABCR alleles with age-related macular degeneration. The International ABCR Screening Consortium. Am J Hum Genet. 2000, 67: 487-491. 10.1086/303018.
Webster AR, Heon E, Lotery AJ, Vandenburgh K, Casavant TL, Oh KT, Beck G, Fishman GA, Lam BL, Levin A, Heckenlively JR, Jacobson SG, Weleber RG, Sheffield VC, Stone EM: An analysis of allelic variation in the ABCA4 gene. Invest Ophthalmol Vis Sci. 2001, 42: 1179-1189.
Rivera A, White K, Stohr H, Steiner K, Hemmrich N, Grimm T, Jurklies B, Lorenz B, Scholl HP, Apfelstedt-Sylla E, Weber BH: A comprehensive survey of sequence variation in the ABCA4 (ABCR) gene in Stargardt disease and age-related macular degeneration. Am J Hum Genet. 2000, 67: 800-813. 10.1086/303090.
Stone EM, Webster AR, Vandenburgh K, Streb LM, Hockey RR, Lotery AJ, Sheffield VC: Allelic variation in ABCR associated with Stargardt disease but not age-related macular degeneration [letter]. Nat Genet. 1998, 20: 328-329. 10.1038/3798.
De La Paz MA, Guy VK, Abou-donia SM, Heinis R, Bracken B, Vance JM, Gilbert JR, Gass JDM, Haines JL, Pericak-Vance MA: Stargardt disease gene (ABCR) mutations in age-related macular degeneration. Ophthalmology. 1999, 106: 1531-1536. 10.1016/S0161-6420(99)90449-9.
Guymer RH, Heon E, Lotery AJ, Munier FL, Schorderet DF, Baird PN, McNeil RJ, Haines H, Sheffield VC, Stone EM: Variation of codons 1961 and 2177 of the Stargardt disease gene is not associated with age-related macular degeneration. Arch Ophthalmol. 2001, 119: 745-751.
Kuroiwa S, Kojima H, Kikuchi T, Yoshimura N: ATP binding cassette transporter retina genotypes and age related macular degeneration: an analysis on exudative non-familial Japanese patients. Br J Ophthalmol. 1999, 83: 613-615.
Schmidt S, Postel EA, Agarwal A, Allen I.C.,Jr., Walters SN, De La Paz MA, Scott WK, Haines JL, Pericak-Vance MA, Gilbert JR: Detailed analysis of allelic variation in the ABCA4 gene in age-related maculopathy. Invest Ophthalmol Vis Sci. 2003, 44: 2868-2875. 10.1167/iovs.02-0957.
Klaver CC, Kliffen M, Van Duijn CM, Hofman A, Cruts M, Grobbee DE, Van, Broeckhoven C, de Jong PT: Genetic association of apolipoprotein E with age-related macular degeneration [published erratum appears in Am J Hum Genet 1998 Oct;63(4):1252]. Am J Hum Genet. 1998, 63: 200-206. 10.1086/301901.
Souied EH, Benlian P, Amouyel P, Feingold J, Lagarde JP, Munnich A, Kaplan J, Coscas G, Soubrane G: The epsilon4 allele of the apolipoprotein E gene as a potential protective factor for exudative age-related macular degeneration. Am J Ophthalmol. 1998, 125: 353-359. 10.1016/S0002-9394(99)80146-9.
Schmidt S, Klaver C, Saunders A, Postel E, De La Paz M., Agarwal A, Small K, Udar N, Ong J, Chalukya M, Nesburn A, Kenney C, Domurath R, Hogan M, Mah T, Conley Y, Ferrell R, Weeks D, De Jong P, van Duijn C, Haines J, Pericak-Vance M, Gorin M: A pooled case-control study of the apolipoprotein E (APOE) gene in age- related maculopathy. Ophthalmic Genet. 2002, 23: 209-223. 10.1076/opge.18.104.22.16883.
Simonelli F, Margaglione M, Testa F, Cappucci G, Manitto MP, Brancato R, Rinaldi E: Apolipoprotein E polymorphisms in age-related macular degeneration in an Italian population. Ophthalmic Res. 2001, 33: 325-328. 10.1159/000055688.
Baird PN, Guida E, Cain M, Chu DT, Mukesh BN, Guymer RH: Association studies of the apolipoprotein E (ApoE) gene and Age related macular degeneration (AMD). Invest Ophthalmol Vis Sci. 2002, 43: 2833-
Magnusson KP, Sigurdsson H, Smarason S, Stefansson K, Group AMD Genetics Research: Genetic association of apolipoprotein E with severe exudative Age-related macular degeneration (AMD). Invest Ophthalmol Vis Sci. 2002, 43: 1843-
Schmidt S, Saunders AM, De La Paz MA, Postel EA, Heinis RM, Agarwal A, Scott WK, Gilbert JR, McDowell JG, Bazyk A, Gass JD, Haines JL, Pericak-Vance MA: Association of the apolipoprotein E gene with age-related macular degeneration: possible effect modification by family history, age, and gender. Mol Vis. 2000, 6: 287-293.
Zareparsi S, Reddick AC, Branham KE, Moore KB, Jessup L, Thoms S, Smith-Wheelock M, Yashar BM, Swaroop A: Association of apolipoprotein e alleles with susceptibility to age-related macular degeneration in a large cohort from a single center. Invest Ophthalmol Vis Sci. 2004, 45: 1306-1310. 10.1167/iovs.03-1253.
Baird PN, Guida E, Chu DT, Vu HT, Guymer RH: The epsilon2 and epsilon4 Alleles of the Apolipoprotein Gene Are Associated with Age-Related Macular Degeneration. Invest Ophthalmol Vis Sci. 2004, 45: 1311-1315. 10.1167/iovs.03-1121.
Schultz DW, Klein ML, Humpert A, Majewski J, Schain M, Weleber RG, Ott J, Acott TS: Lack of an association of apolipoprotein E gene polymorphisms with familial age-related macular degeneration. Arch Ophthalmol. 2003, 121: 679-683. 10.1001/archopht.121.5.679.
Pang CP, Baum L, Chan WM, Lau TC, Poon PM, Lam DS: The apolipoprotein E epsilon4 allele is unlikely to be a major risk factor of age-related macular degeneration in Chinese. Ophthalmologica. 2000, 214: 289-291. 10.1159/000027506.
Scott WK, Schmidt S, Fan Y-T, Postel EA, Agarwal A, Gass JDM, Gilbert JR, Haines JL, Pericak-Vance MA: Cigarette smoking and APOE genotype interaction in age-related macular degeneration. Invest Ophthalmol Vis Sci 2004, 45: 2302. 2004
Weeks DE, Conley YP, Tsai HJ, Mah TS, Rosenfeld PJ, Paul TO, Eller AW, Morse LS, Dailey JP, Ferrell RE, Gorin MB: Age-related maculopathy: an expanded genome-wide scan with evidence of susceptibility loci within the 1q31 and 17q25 regions. Am J Ophthalmol. 2001, 132: 682-692. 10.1016/S0002-9394(01)01214-4.
Majewski J, Schultz DW, Weleber RG, Schain MB, Edwards AO, Matise TC, Acott TS, Ott J, Klein ML: Age-related macular degeneration--a genome scan in extended families. Am J Hum Genet. 2003, 73: 540-550. 10.1086/377701.
Schick JH, Iyengar SK, Klein BE, Klein R, Reading K, Liptak R, Millard C, Lee KE, Tomany SC, Moore EL, Fijal BA, Elston RC: A whole-genome screen of a quantitative trait of age-related maculopathy in sibships from the Beaver Dam Eye Study. Am J Hum Genet. 2003, 72: 1412-1424. 10.1086/375500.
Seddon JM, Santangelo SL, Book K, Chong S, Cote J: A genomewide scan for age-related macular degeneration provides evidence for linkage to several chromosomal regions. Am J Hum Genet. 2003, 73: 780-790. 10.1086/378505.
Weeks DE, Conley YP, Mah TS, Paul TO, Morse L, Ngo-Chang J, Dailey JP, Ferrell RE, Gorin MB: A full genome scan for age-related maculopathy. Hum Mol Genet. 2000, 9: 1329-1349. 10.1093/hmg/9.9.1329.
Iyengar SK, Song D, Klein BE, Klein R, Schick JH, Humphrey J, Millard C, Liptak R, Russo K, Jun G, Lee KE, Fijal B, Elston RC: Dissection of genomewide-scan data in extended families reveals a major locus and oligogenic susceptibility for age-related macular degeneration. Am J Hum Genet. 2004, 74: 20-39. 10.1086/380912.
Schultz DW, Klein ML, Humpert AJ, Luzier CW, Persun V, Schain M, Mahan A, Runckel C, Cassera M, Vittal V, Doyle TM, Martin TM, Weleber RG, Francis PJ, Acott TS: Analysis of the ARMD1 locus: evidence that a mutation in HEMICENTIN-1 is associated with age-related macular degeneration in a large family. Hum Mol Genet. 2003, 12: 3315-3323. 10.1093/hmg/ddg348.
Hauser ER, Watanabe RM, Duren WL, Bass MP, C.D. Langefeld, Boehnke M: Ordered subset analysis in genetic linkage mapping of complex traits. Genet Epidemiol. 2004, 27: 53-63-10.1002/gepi.20000.
Smith W, Assink J, Klein R, Mitchell P, Klaver CC, Klein BE, Hofman A, Jensen S, Wang JJ, de Jong PT: Risk factors for age-related macular degeneration: Pooled findings from three continents. Ophthalmology. 2001, 108: 697-704. 10.1016/S0161-6420(00)00580-7.
Hyman L, Schachat AP, He Q, Leske MC: Hypertension, cardiovascular disease, and age-related macular degeneration. Age-Related Macular Degeneration Risk Factors Study Group. Arch Ophthalmol. 2000, 118: 351-358.
Hennis A, Wu SY, Nemesure B, Leske MC: Hypertension, diabetes, and longitudinal changes in intraocular pressure. Ophthalmology. 2003, 110: 908-914. 10.1016/S0161-6420(03)00075-7.
Nemesure B, Wu SY, Hennis A, Leske MC: Factors related to the 4-year risk of high intraocular pressure: the Barbados Eye Studies. Arch Ophthalmol. 2003, 121: 856-862. 10.1001/archopht.121.6.856.
van Leeuwen R, Ikram MK, Vingerling JR, Witteman JC, Hofman A, de Jong PT: Blood pressure, atherosclerosis, and the incidence of age-related maculopathy: the Rotterdam Study. Invest Ophthalmol Vis Sci. 2003, 44: 3771-3777. 10.1167/iovs.03-0121.
Klein R, Klein BE, Tomany SC, Cruickshanks KJ: The association of cardiovascular disease with the long-term incidence of age-related maculopathy: the Beaver Dam Eye Study. Ophthalmology. 2003, 110: 1273-1280. 10.1016/S0161-6420(03)00599-2.
Group Age-Related Eye Disease Study Research: Risk factors associated with age-related macular degeneration; a case-control study in the age-related eye disease study: Age Related Eye Disease Study Report number 3. Ophthalmology. 2000, 107: 2224-2232. 10.1016/S0161-6420(00)00409-7.
Elston RC, Buxbaum S, Jacobs KB, Olson JM: Haseman and Elston revisited. Genet Epidemiol. 2000, 19: 1-17. 10.1002/1098-2272(200007)19:1<1::AID-GEPI1>3.0.CO;2-E.
Schaumberg DA, Christen WG, Hankinson SE, Glynn RJ: Body mass index and the incidence of visually significant age-related maculopathy in men. Arch Ophthalmol. 2001, 119: 1259-1265.
Delcourt C, Michel F, Colvez A, Lacroux A, Delage M, Vernet MH: Associations of cardiovascular disease and its risk factors with age-related macular degeneration: the POLA study. Ophthalmic Epidemiol. 2001, 8: 237-249. 10.1076/opep.22.214.171.1243.
Hammond B.R.,Jr., Ciulla TA, Snodderly DM: Macular pigment density is reduced in obese subjects. Invest Ophthalmol Vis Sci. 2002, 43: 47-50.
Beatty S, Murray IJ, Henson DB, Carden D, Koh H-H, Boulton ME: Macular Pigment and Risk for Age-Related Macular Degeneration in Subjects from a Northern European Population. Invest Ophthalmol Vis Sci. 2001, 42: 439-446.
I.J. Suner, D.G. Epsinosa-Heidmann, M.E. Marin-Castano, E.P. Hernandez, Pereira-Simon S, Cousins SW: Nicotine Increases Size and Severity of Experimental Choroidal Neovascularization. Invest Ophthalmol Vis Sci. 2004, 45: 311-317. 10.1167/iovs.03-0733.
Chobanian AV, Bakris GL, Black HR, Cushman WC, Green LA, Izzo J.L.,Jr., Jones DW, Materson BJ, Oparil S, Wright J.T.,Jr., Roccella EJ: The Seventh Report of the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure: the JNC 7 report. JAMA. 2003, 289: 2560-2572. 10.1001/jama.289.19.2560.
Klein R, Klein BE, Moss SE: The relation of systemic hypertension to changes in the retinal vasculature: the Beaver Dam Eye Study. Trans Am Ophthalmol Soc. 1997, 95: 329-348.
Leung H, Wang JJ, Rochtchina E, Tan AG, Wong TY, Klein R, Hubbard LD, Mitchell P: Relationships between age, blood pressure, and retinal vessel diameters in an older population. Invest Ophthalmol Vis Sci. 2003, 44: 2900-2904. 10.1167/iovs.02-1114.
Hubbard LD, Brothers RJ, King WN, Clegg LX, Klein R, Cooper LS, Sharrett AR, Davis MD, Cai J: Methods for evaluation of retinal microvascular abnormalities associated with hypertension/sclerosis in the Atherosclerosis Risk in Communities Study. Ophthalmology. 1999, 106: 2269-2280. 10.1016/S0161-6420(99)90525-0.
Wong TY, Klein R, Sharrett AR, Duncan BB, Couper DJ, Tielsch JM, Klein BE, Hubbard LD: Retinal arteriolar narrowing and risk of coronary heart disease in men and women. The Atherosclerosis Risk in Communities Study. JAMA. 2002, 287: 1153-1159. 10.1001/jama.287.9.1153.
Wang JJ, Mitchell P, Rochtchina E, Tan AG, Wong TY, Klein R: Retinal vessel wall signs and the 5 year incidence of age related maculopathy: the Blue Mountains Eye Study. Br J Ophthalmol. 2004, 88: 104-109. 10.1136/bjo.88.1.104.
Klein BE, Klein R, Lee KE: Heritability of Risk Factors for Primary Open-Angle Glaucoma: The Beaver Dam Eye Study. Invest Ophthalmol Vis Sci. 2004, 45: 59-62. 10.1167/iovs.03-0516.
Klein BE, Klein R, Linton KL: Intraocular pressure in an American community. The Beaver Dam Eye Study. Invest Ophthalmol Vis Sci. 1992, 33: 2224-2228.
Smith W, Mitchell P, Leeder SR, Wang JJ: Plasma fibrinogen levels, other cardiovascular risk factors, and age-related maculopathy: The Blue Mountains Eye Study. Arch Ophthalmol. 1998, 116: 583-587.
Klein R, Klein BE, Jensen SC: The relation of cardiovascular disease and its risk factors to the 5-year incidence of age-related maculopathy: the Beaver Dam Eye Study. Ophthalmology. 1997, 104: 1804-1812.
Atwood LD, Heard-Costa NL, Cupples LA, Jaquish CE, Wilson PW, D'Agostino RB: Genomewide linkage analysis of body mass index across 28 years of the Framingham Heart Study. Am J Hum Genet. 2002, 71: 1044-1050. 10.1086/343822.
Scott WK, Hauser ER, Schmechel DE, Welsh-Bohmer KA, Small GW, Roses AD, Saunders AM, Gilbert JR, Vance JM, Haines JL, Pericak-Vance MA: Ordered-subsets linkage analysis detects novel Alzheimer disease Loci on chromosomes 2q34 and 15q22. Am J Hum Genet. 2003, 73: 1041-1051. 10.1086/379083.
Witte JS, Goddard KA, Conti DV, Elston RC, Lin J, Suarez BK, Broman KW, Burmester JK, Weber JL, Catalona WJ: Genomewide scan for prostate cancer-aggressiveness loci. Am J Hum Genet. 2000, 67: 92-99. 10.1086/302960.
Shao Y, Cuccaro ML, Hauser ER, Raiford KL, Menold MM, Wolpert CM, Ravan SA, Elston L, Decena K, Donnelly SL, Abramson RK, Wright HH, DeLong GR, Gilbert JR, Pericak-Vance MA: Fine mapping of autistic disorder to chromosome 15q11-q13 by use of phenotypic subtypes. Am J Hum Genet. 2003, 72: 539-548. 10.1086/367846.
Ott J: Linkage analysis and family classification under heterogeneity. Ann Hum Genet. 1983, 47: 311-320.
De La Paz MA, Pericak-Vance MA, Lennon F, Haines JL, Seddon JM: Exclusion of TIMP3 as a candidate locus in age-related macular degeneration. Invest Ophthalmol Vis Sci. 1997, 38: 1060-1065.
The Age-Related Eye Disease Study (AREDS): design implications AREDS report no. 1. The Age-Related Eye Disease Study Research Group. Control Clin Trials. 1999, 20: 573-600. 10.1016/S0197-2456(99)00031-8.
Klein R, Davis MD, Magli YL, Segal P, Klein BE, Hubbard L: The Wisconsin age-related maculopathy grading system. Ophthalmology. 1991, 98: 1128-1134.
Bird AC, Bressler NM, Bressler SB, Chisholm IH, Coscas G, Davis MD, de Jong PT, Klaver CC, Klein BE, Klein R: An international classification and grading system for age-related maculopathy and age-related macular degeneration. The International ARM Epidemiological Study Group. Survey of Ophthalmology. 1995, 39: 367-374.
Epstein MP, Duren WL, Boehnke M: Improved inference of relationship for pairs of individuals. Am J Hum Genet. 2000, 67: 1219-1231.
McPeek MS, Sun L: Statistical tests for detection of misspecified relationships by use of genome-screen data. Am J Hum Genet. 2000, 66: 1076-1094. 10.1086/302800.
O'Connell JR, Weeks DE: PedCheck: A program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet. 1998, 63: 259-266. 10.1086/301904.
Broman KW, Murray JC, Sheffield VC, White RL, Weber JL: Comprehensive human genetic maps: Individuals and sex-specific variation in recombination. Am J Hum Genet. 1998, 63: 861-869. 10.1086/302011.
Risch N: Linkage strategies for genetically complex traits. II. The power of affected relative pairs. Am J Hum Genet. 1990, 46: 229-241.
Hauser ER, Boehnke M: Genetic linkage analysis of complex genetic traits by using affected sibling pairs. Biometrics. 1998, 54: 1238-1246.
We thank the patients with AMD and their family members for participating in this research, and the personnel at the Duke Center for Human Genetics for assistance in patient and family ascertainment, data management, and preparation of this manuscript. Grant support from the National Eye Institute (EY12118, to MAP-V; EY09859, to MBG), the National Institute on Aging (AG11268, to Harvey Cohen, MD), The Eye and Ear Foundation, Pittsburgh, PA (to MBG), and Research to Prevent Blindness (to MBG) is gratefully acknowledged. The OSA software is available at http://wwwchg.duhs.duke.edu/software/osa.html.
SS participated in the design of the study, performed all statistical analyses and drafted the manuscript. WKS designed the study questionnaire used to collect environmental risk factor information (such as smoking history) and assisted in the interpretation of results. EAP, AA and MADLP were responsible for the clinical evaluation of all study participants. ERH developed the OSA method and assisted in the interpretation of results. JRG participated in the design of the study. DEW and MBG coordinated the genomic screening of study samples performed at CIDR. JLH and MAP-V conceived of the study and oversaw its overall design and coordination. All authors read and approved the final manuscript.
About this article
Cite this article
Schmidt, S., Scott, W.K., Postel, E.A. et al. Ordered subset linkage analysis supports a susceptibility locus for age-related macular degeneration on chromosome 16p12. BMC Genet 5, 18 (2004). https://doi.org/10.1186/1471-2156-5-18