Age-related macular degeneration (AMD) is a disease of the retina characterized by lipid-rich deposits called drusen which accumulate beneath the basal lamina of the retinal pigment epithelium (RPE)[
1], causing progressive loss of central vision[
2]. Clinically, it is subdivided into early stage AMD and late stage AMD[
2], and late stage can be further subdivided into nonneovascular and neovascular AMD[
2]. It is a leading cause of blindness affecting approximately 8.4 million individuals worldwide, especially in regions of the world with high proportions of elderly people where up to 14% of blindness cases are caused by the disease[
3]. However, other than age, studies have shown that genetic risk factors contribute significantly to risk of AMD[
4].
Genetic risk factors for AMD include genes involved in lipid metabolism including the cholesteryl ester transfer protein gene
(CETP)[
5], the hepatic lipase gene
(LIPC), the apolipoprotein E gene
(APOE), and the adenosine triphosphate-binding cassette transporter A1 (
ABCA1) gene [
4,
6]. These genes are all involved in the metabolism of high-density lipoprotein (HDL) cholesterol, however, the exact mechanism underlying the involvement of these genes in the pathogenesis of AMD remains unknown. ABCA1 is a cholesterol transporter that transports cholesterol onto lipid-poor apolipoproteins[
7,
8] and is expressed in various tissues including the eye, brain, and liver. Homozygotes for loss-of-function mutations in
ABCA1 have Tangier disease, a disease characterized by extremely low plasma concentrations of HDL cholesterol and accumulation of cholesterol in the tonsils, liver, spleen, lymph nodes, and cornea[
9]. Intronic variants in
ABCA1 are reported to be associated with risk of AMD in genome-wide association studies (GWAS)[
4,
6], but whether this also applies to amino acid-changing
ABCA1 variants in the general population is unknown.
Using two large population cohorts including 90,556 individuals, we tested the hypothesis that amino acid-changing genetic variants in ABCA1 associated with high HDL cholesterol concentrations are also associated with risk of AMD in the general population.
Methods
Setting, study population, and outcome
We included 90,556 individuals from two similar prospective cohort studies of the Danish general population, the Copenhagen City Heart Study (CCHS, n = 9,584) and the Copenhagen General Population Study (CGPS, n = 80,972). We excluded individuals with prevalent AMD (n = 50) at baseline. Institutional Review Board and Ethics Committee approval was obtained (KF-100-2039/91 and H-KF-01-144/01) and the study was conducted according to the Declaration of Helsinki. Written informed consent was obtained from all participants. All participants were White and of Danish descent. For an overview of the study design see Supplementary Fig. 1.
CGPS
The CGPS was initiated in 2003 with ongoing enrolment. Individuals were selected based on the national Danish Civil Registration System to reflect the adult Danish population aged 20–100 years or older. Data were obtained from a questionnaire, a physical examination, and a blood sample. We included 80,972 consecutive participants from this study in the present analysis. A total number of 1,370 AMD patients were diagnosed by the end of follow-up in the CGPS. Median follow-up was 10 years (< 1 to 15 years).
CCHS
The CCHS was initiated in 1976-78, with follow-up examinations in 1981-83, 1991-94, and 2001-03. Participants were recruited and examined as in the CGPS. We included all 9,584 consecutive participants in the 1991-94 and 2001-03 examinations in the present analysis. A total number of 142 AMD patients were diagnosed by the end of follow-up in the CCHS. Median follow-up was 18 years (< 1 to 27 years).
Endpoint definition
A diagnosis of AMD (World Health Organization, International Classification of Diseases, Tenth Revision [ICD-10]: DH353C, DH353E, DH353J, DH353K, DH353L, and DH353M, including senile macular degeneration (atrophic or exudative, with and without subretinal neovascularization and/or vitreomacular traction) was ascertained using the national Danish Patient Registry, which registers all contacts and diagnoses from public and private hospitals, from study entry until December 18th, 2018; not a single person was lost to follow-up due to the unique Danish registries. Danish national guidelines direct private practicing ophthalmologists to refer patients suspected of having exudative or neovascular AMD with active choroidal neovascularization to anti-VEGF treatment centralized at hospitals. This secures that most patients with AMD are registered with a diagnosis at a hospital, while all diagnoses of AMD in Denmark are given by medical doctors specialized in ophthalmology. The referral decision is (as a minimum) based on (1) visual acuity measurement, (2) slit lamp examination, ophthalmoscopy, or color fundus photography after pupil dilation, and (3) optical coherence tomography scan. Neovascular AMD was ICD10: DH353C, DH353J, DH353K, and DH353M. Nonneovascular AMD was ICD10: DH353E and DH353L. Individuals diagnosed with both nonneovascular and neovascular AMD were included in analyses of both endpoints. There were 391 individuals in CGPS and 35 individuals in CCHS who were diagnosed with both nonneovascular and neovascular AMD. The diagnoses of skin cancer were classified according to the World Health Organization international classification of diseases 10th edition (ICD-10): non-melanoma skin cancer as C44 and melanoma as C43. None of the non-melanoma skin cancer or melanoma diagnoses were based on self-report, and all were from the Danish Cancer Registry.
Genotyping
All 90,556 individuals were genotyped for fourteen amino acid-changing variants in
ABCA1: rs145183203, rs2230806, rs9282543, rs2066718, rs35819696, rs138880920, rs2066715, rs2066714, rs1422877738, rs33918808, rs76881554, rs2230808, rs146292819, and rs142688906. Briefly, these fourteen variants were selected because they were amino acid-changing and had a minor allele frequency above 0.001. From these fourteen variants nine were chosen for further analyses because of a significant association with HDL cholesterol concentrations in plasma (rs145183203, rs2066718, rs35819696, rs2066715, rs2066714, rs33918808, rs76881554, rs2230808, rs146292819). Genotyping was by TaqMan-based assays (Applied Biosystems, Foster City, CA, USA), or by an allele-specific PCR (Fluidigm; BioXpedia, Aarhus, Denmark). See Supplementary Fig. 2 for an overview of the genotyping and selection of variants for further analyses. Rs429358 and rs7412 in
APOE, rs708272 in
CETP, and rs1800588 in
LIPC were genotyped as previously described[
5,
10,
11].
Other covariates
Covariates for baseline characteristics were risk factors for atherosclerotic cardiovascular disease and AMD: age, sex, smoking status, alcohol consumption, use of lipid-lowering therapy, diabetes, hypertension, body mass index, physical inactivity, education, postmenopausal status, and hormonal replacement therapy. These risk factors were chosen since they have all been shown to be associated with risk of AMD[
2,
12‐
16]. Hypertension was use of anti-hypertensive medication and/or a systolic blood pressure of ≥ 140 mmHg, and/or a diastolic blood pressure of ≥ 90 mmHg. The reference was a systolic blood pressure < 140 mmHg or a diastolic blood pressure < 90 mmHg. Diabetes mellitus was self-reported disease, use of insulin or oral hypoglycaemic agents, and/or non-fasting plasma glucose levels of > 11 mmol/L (198 mg/dL). The reference was no self-reported disease or non-fasting plasma glucose levels ≤ 11 mmol/L (198 mg/dL). Smoking was current smoking. The reference was non-smokers. High alcohol consumption was > 14/21 units per week for women/men (1 unit = 12 g alcohol, equivalent to one glass of wine, one shot of spirit, or one beer (33 cL)). The reference was ≤ 14/21 units per week for women/men. Physical inactivity was ≤ four hours per week of light physical activity in leisure time. The reference was > four hours per week of light physical activity in leisure time. Women reported menopausal status (the reference was premenopausal women) and use of hormonal replacement therapy (in women who were postmenopausal). The reference was no use of hormonal replacement therapy (in women who were postmenopausal). Lipid-lowering therapy was primarily statins (yes/no). The reference was no use of lipid-lowering therapy.
Statistical analyses
R version 3.6.1 or Stata SE 17 was used. Chi-squared tests evaluated Hardy-Weinberg equilibrium. Kruskal-Wallis test or Cuzick’s test for trend was used to compare continuous covariates and Pearson’s chi-squared test to compare categorical covariates by genotype or weighted allele score. Missing data on covariates (0.4%) were imputed based on all other covariates in the study and their distributions using multivariable linear regression for continuous variables. Categorical variables where assigned a separate category. The mi impute command in Stata was used. If only individuals with complete data were included (no use of imputed data), results were like those reported.
In genetic analyses, we defined the exposure allele for each of the nine variants as the allele associated with the highest HDL cholesterol. For an in-depth description of the construction of the genetic score see Supplementary Table 1. For each individual, the nine variants were combined in an allele score weighted on higher HDL cholesterol (Supplementary Table 1) by multiplying the adjusted β-coefficient (see below) for each variant by the frequency of that variant and adding the adjusted β-coefficients across all variants. The adjusted β-coefficient for a specific variant takes the linkage structure from all nine variants into account by adjusting for the effects of all other variants in a linear regression analysis. The weighted allele scores across all nine variants for the different combined genotypes observed were then divided into groups of approximately equal size (called tertiles in the following): adjusted betascore from − 0.44 to < 0 (first tertile = reference), 0 (second tertile), and from > 0 to 0.13 (third tertile). External weights from
https://app.genebass.org/ were used to create a weighted allele score in a sensitivity analysis (Supplementary Table 1).
To examine the possible causal association between lifelong amino acid-changing genetic variants in
ABCA1 and risk of AMD we used Cox proportional hazard models with age as time scale and delayed entry at examination (left truncation). For Cox regression models, proportionality of hazards over time was assessed by plotting -ln(-ln[survival]) versus ln(analysis time). There was no suspicion of non-proportionality. All models were age- and sex adjusted and multivariable adjusted for risk factors for cardiovascular disease (other than lipids) and AMD: age (as time scale), sex, cohort, date of birth, body mass index, hypertension, diabetes, physical activity, smoking, alcohol consumption, postmenopausal status (in women), use of hormone replacement therapy (in women), education, and use of lipid-lowering therapy. In sensitivity analyses we further adjusted for
APOE, CETP, and
LIPC genotype, for LDL- and total cholesterol, as well as for skin cancer as a proxy for extensive sun exposure. The association between the
ABCA1 weighted allele score and risk of AMD was examined on a continuous scale using restricted cubic splines incorporated in the Cox regression model. The number of knots was chosen according to Akaike information criteria[
17]. The reference was set to a weighted allele score of 0 to have the largest possible reference group (Supplementary Table 1, panel B).
Interactions were tested using an interaction term between the ABCA1 weighted allele score and covariates, one at a time, in a multivariable adjusted Cox model. When investigating the interaction between age and the ABCA1 weighted allele score on AMD, time since study entry was used as the time scale.
In mediation analyses, the percent of the total risk of AMD conferred by genetic variants in
ABCA1 which was associated with genetically higher HDL cholesterol, was determined by the Karlson-Holm-Breen method[
18], using bootstrapping to calculate confidence intervals.
Associations of plasma HDL cholesterol (mmol/L) with LDL- and total cholesterol concentrations (mmol/L) were investigated using multiple linear regression and graphically displayed using kernel-weighted local polynomial smoothing and geometric means with 95% confidence intervals (CIs).
Discussion
The principal findings are that amino acid-changing genetic variants in ABCA1 which were associated with higher concentrations of HDL cholesterol, were associated with higher risk of AMD including both nonneovascular and neovascular AMD in the Danish general population. We observed these findings of genetically determined HDL cholesterol concentrations mediated by ABCA1 both on a continuous scale and in categories. These findings are novel, because we only use amino acid-changing variants in ABCA1, variation that only affects HDL cholesterol concentrations and no other lipoproteins, thereby suggesting a likely causal relationship from ABCA1 over HDL cholesterol to risk of AMD.
In support of our findings, the concentrations of various plasma lipids have been associated with risk of AMD, but the most robust association is observed between high plasma HDL cholesterol concentrations and higher risk of AMD in both observational and genetic analyses[
5,
19‐
24], as also confirmed in the present study. Genes involved in HDL cholesterol metabolism that have also been associated with risk of AMD include
CETP, LIPC, APOE, and
ABCA1[
5,
23]. The
ABCA1 gene is a significant signal in GWAS for AMD[
4,
6], and in smaller studies, single variants in
ABCA1 are associated with risk of AMD[
25] and progression from one stage of the disease to another[
26]. While Mendelian randomization studies investigating the association between genetically determined HDL cholesterol concentration and risk of AMD have included
ABCA1 variants, no study has, to our knowledge, investigated the association between a weighted allele score based solely on amino acid-changing variants in
ABCA1 and risk of AMD[
22‐
24]. This is important, as it may suggest direct causality from
ABCA1 over HDL to risk of AMD, and therefore highlights the possibility of using ABCA1 as a future drug target for preventing or maybe treating AMD.
The mechanism behind the association between genetically determined high concentrations of HDL cholesterol due to variation in
ABCA1 and risk of AMD is speculative, but below we describe a likely scenario. According to the “Oil spill hypothesis”, HDL particles are involved in the delivery of lipophilic nutrients to the retinal pigment epithelium (RPE) which is responsible for the subsequent delivery to the photoreceptors[
1,
27]. The lipophilic nutrients are taken up by scavenger receptor class B type I (SR-BI) situated at Bruch’s membrane. Proteins involved in HDL metabolism are localized in the subretinal space and are thought to supply unesterified cholesterol to the retina. However, excess unesterified cholesterol is cytotoxic and must be offloaded from the RPE[
1]. Mechanisms for offloading this excess unesterified cholesterol could include ABCA1-mediated transfer to circulating HDL particles[
1] and to HDL particles within the RPE. The cholesterol content in these particles will be esterified, and over time clearance of lipoprotein particles across Bruch’s membrane is hindered by the accumulation of esterified cholesterol causing lipoproteins to be captured in Bruch’s membrane and forming a lipid wall external to the RPE basal lamina[
1]. Subsequently, lipoproteins are degraded forming lipoprotein-derived debris. This eventually leads to the formation of the AMD-specific lesions called drusen[
1]. Peroxidised lipids in drusen are thought to cause choroidal neovascularization through inflammation[
1], which is characteristic of the type of AMD called “neovascular” AMD. Drusen will also lead to the blocking of transport of necessary nutrients from plasma to photoreceptors[
1], and with time this will lead to neurodegeneration of photoreceptors and ultimately blindness characteristic of the type of AMD called “nonneovascular” AMD (Supplementary Fig. 10). Although speculative, it is thus biologically plausible that
ABCA1 variants associated with high HDL concentration may also be involved directly in the pathogenesis of AMD.
Strengths of the current study are the inclusion of a large homogeneous sample from the general population, the thorough characterization of risk factors for AMD on an individual participant level, and the prospective study design. Another strength is the adjustment for important risk factors for AMD and the additional adjustment for other important risk genes. Further, unlike previous studies, we investigated several amino acid-changing variants in ABCA1 spanning the whole protein coding part of the gene.
Limitations include studying only White individuals of Danish descent; however, we are not aware of data suggesting that the present results should not apply to other ethnicities. Another limitation is the lack of functional studies of these structural genetic variants in
ABCA1; however, the association with HDL cholesterol and apolipoprotein A1 concentrations in individuals in the general population is a reasonable proxy for functionality. It should be noted that there are other indicators of blood lipids besides HDL cholesterol. These include LDL cholesterol and total cholesterol that are only weakly correlated with HDL cholesterol. It is, however, unlikely that parts of the effect of
ABCA1 variants on risk of AMD are mediated through other lipids since the direct associations between
ABCA1 variants and other lipids and lipoproteins are negliable. Another limitation is that family history of AMD and known AMD loci as complement factor H, ARMS2, and VEGF are not included in the multivariable adjustment used throughout this study. Complement factor H, ARMS2, and VEGF are, however, encoded by genes located on other chromosomes than
ABCA1, and are thus unlikely to bias the current results through linkage disequilibrium with
ABCA1 variants. Further, genetic results are very unlikely to be confounded[
28] which is also clearly demonstrated by the similar results for three different adjustments in Fig.
1 and in Supplementary Figs. 5 and 6.
In conclusion, amino acid-changing genetic variants in ABCA1 were associated with higher concentrations of plasma HDL cholesterol and higher risk of both nonneovascular and neovascular AMD in the Danish general population. These results may suggest causality from ABCA1 over HDL cholesterol to risk of AMD and highlights the possibility of using ABCA1 as a future drug target for preventing or maybe treating AMD.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.