Systems biology-based analysis implicates a novel role for vitamin D metabolism in the pathogenesis of age-related macular degeneration.

Vitamin D has been shown to have anti-angiogenic properties and to play a protective role in several types of cancer, including breast, prostate and cutaneous melanoma. Similarly, vitamin D levels have been shown to be protective for risk of a number of conditions, including cardiovascular disease and chronic kidney disease, as well as numerous autoimmune disorders such as multiple sclerosis, inflammatory bowel diseases and type 1 diabetes mellitus. A study performed by Parekh et al. was the first to suggest a role for vitamin D in age-related macular degeneration (AMD) and showed a correlation between reduced serum vitamin D levels and risk for early AMD. Based on this study and the protective role of vitamin D in diseases with similar pathophysiology to AMD, we examined the role of vitamin D in a family-based cohort of 481 sibling pairs. Using extremely phenotypically discordant sibling pairs, initially we evaluated the association of neovascular AMD and vitamin D/sunlight-related epidemiological factors. After controlling for established AMD risk factors, including polymorphisms of the genes encoding complement factor H (CFH) and age-related maculopathy susceptibility 2/HtrA serine peptidase (ARMS2/HTRA1), and smoking history, we found that ultraviolet irradiance was protective for the development of neovascular AMD (p = 0.001). Although evaluation of serum vitamin D levels (25-hydroxyvitamin D [25(OH)D]) was higher in unaffected individuals than in their affected siblings, this finding did not reach statistical significance. Based on the relationship between ultraviolet irradiance and vitamin D production, we employed a candidate gene approach for evaluating common variation in key vitamin D pathway genes (the genes encoding the vitamin D receptor [VDR]; cytochrome P450, family 27, subfamily B, polypeptide 1 [CYP27B1]; cytochrome P450, family 24, subfamily A, polypeptide 1 [CYP24A1]; and CYP27A1) in this same family-based cohort. Initial findings were then validated and replicated in the extended family cohort, an unrelated case-control cohort from central Greece and a prospective nested case-control population from the Nurse's Health Study and Health Professionals Follow-Up Studies, which included patients with all subtypes of AMD for a total of 2,528 individuals. Single point variants in CYP24A1 (the gene encoding the catabolising enzyme of the vitamin D pathway) were demonstrated to influence AMD risk after controlling for smoking history, sex and age in all populations, both separately and, more importantly, in a meta-analysis. This is the first report demonstrating a genetic association between vitamin D metabolism and AMD risk. These findings were also supplemented with expression data from human donor eyes and human retinal cell lines. These data not only extend previous biological studies in the AMD field, but further emphasise common antecedents between several disorders with an inflammatory/immunogenic component such as cardiovascular disease, cancer and AMD.

Introduction

Several biological processes/cellular pathways involved in conditions such as cancer, cardiovascular disease and type 1 diabetes have also been implicated in the pathophysiology of the neovascular form of age-related macular degeneration (AMD). AMD, a disease characterised by the loss of ability to drive, recognise faces and read, is the leading cause of blindness in the US elderly population. Most prominent among these shared disease-associated processes is angiogenesis, the defining hallmark of neovascular AMD [1-3]. Additionally, immune cell involvement, inflammation and oxidative stress have all been associated with these diseases [2,4-19].

Vitamin D is a circulating steroid hormone which is thought to be protective in a variety of cancers, including breast, prostate and cutaneous mela-noma [20-23]. Although classically associated with bone mineralisation and the regulation of Ca2+ and phosphorus (Pi),[24] vitamin D has been demonstrated to have anti-angiogenic and anti-inflammatory effects [25-28]. In addition, it has been shown to play a role in cellular proliferation, differentiation and apoptosis, as well as protecting cell proteins and membranes from oxidative damage [20,29,30].

Vitamin D synthesis and metabolism begins upon cutaneous ultraviolet (UV) B exposure, resulting in the conversion of 7-dehydrocholesterol to vitamin D. In addition to UVB, sources of vitamin D include diet and supplements [20]. In the liver, cytochrome P450, family 27, subfamily B, polypeptide 1 (CYP27B1) catalyses the conversion of vitamin D to its pro-hormonal major circulating form, 25-hydroxyvitamin D (25(OH)D). CYP27A1 then converts 25(OH)D to the biologically active form of vitamin D [1,25(OH)2D3] [31]. The biologically active form of vitamin D binds to the nuclear hormone receptor vitamin D receptor gene (VDR). VDR heterodimerises with retinoid X receptor alpha to regulate transcription by binding vitamin D response elements in gene promoters and recruiting transcriptional co-factors [32]. The active form of vitamin D is eventually catabolised after processing in the kidney or liver by cytochrome P450, family 24, subfamily A, polypeptide 1 (CYP24A1) [33]. Parekh et al. demonstrated an inverse relationship between serum vitamin D levels and early AMD but could not extrapolate the findings to advanced AMD, potentially owing to the small number of patients studied in this category (n = 10) [34]. Based on this association and the involvement of vitamin D in processes underlying several diseases with an inflammatory or immune component, we hypothesised that vitamin D might play a role in AMD pathophysiology and, in particular, neovascular AMD. A special focus was placed on neovascular AMD, as it is the advanced form of AMD that accounts for the majority of cases of severe vision loss and it is necessarily preceded by an early form of the disease. To evaluate comprehensively the contribution of vitamin D and vitamin D pathways to neovascular AMD and, moreover, all subtypes of AMD, we conducted epidemiological, biological and genetic studies initially on a family-based cohort. This was followed by replication and validation of significant genetic findings in an unrelated case-control cohort from central Greece, as well as two prospective cohorts comprising AMD cases and controls from the Nurse's Health Study (NHS) and Health Professionals Follow-Up Studies (HPFS).

An epidemiological study was conducted to determine the association of dietary supplements, UV light and sun-related variables, including iris colour, self-reported sun exposure and skin cancer history, with neovascular AMD. Serum levels of vitamin D were also examined to assess further the contribution of vitamin D levels to advanced AMD risk. To assess if genetic variation within vitamin D metabolism genes might be associated with neovascular AMD, tagging single nucleotide polymorphisms (SNPs) within the vitamin D pathway genes CYP27B1, CYP27A1, CYP24A1 and VDR were genotyped to find significantly associated AMD risk variants and haplotypes within a discovery cohort comprising extremely discordant sibling pairs --that is, where one sibling, the index patient, has the neovascular form of AMD and the other sibling has no signs of AMD and is older than 65 years of age. Moreover, CYP27B1, CYP27A1, CYP24A1 and VDR were chosen because they are vitamin D pathway genes that are located in genomic regions previously implicated in AMD disease susceptibility [35,36]. Based on findings of genetic association within this extremely discordant sibpair discovery cohort, we replicated and validated findings in our other cohorts as described below.

Materials and methods

Family-based patient population

The protocol was reviewed and approved by the Institutional Review Boards at the Massachusetts Eye and Ear Infirmary (MEEI), Boston, MA, USA, and conformed to the tenets of the Declaration of Helsinki. Eligible patients were enrolled in this study after giving informed consent and before completing a standardised questionnaire and donating 10-50 ml of venous blood.

Details of the recruitment of the sibling pairs, comprising mainly individuals of European ancestry, are described in detail elsewhere [37,38]. Briefly, the disease status of every participant was confirmed by at least two investigators by evaluation of fundus photographs or fluorescein angiograms, except when one of the investigators directly examined an unaffected sibling during a home visit (n = 4 cases).

All index patients were approximately 50 years of age or older and had the neovascular form of AMD in at least one eye, defined by subretinal haemorrhage, fibrosis or fluorescein angiographic presence of neovascularisation. Patients whose only exudative finding was a retinal pigment epithelium (RPE) detachment were excluded because this finding may not represent definite neovascular AMD. Patients with signs of pathological myopia, presumed ocular histoplasmosis syndrome, angioid streaks, choroidal rupture, any hereditary retinal diseases other than AMD, and previous laser treatment due to retinal conditions other than AMD were also excluded.

Patients with an early/intermediate form of dry AMD were classified according to the Age-Related Eye Disease Study (AREDS)[39]as either category 2 (small [< 63 μm] drusen with total area ≥ 125 μm, or at least one intermediate druse [≥63 and < 125 μm], or presence of pigment abnormalities); or AREDS category 3 (intermediate drusen comprising total area ≥360 μm in the presence of soft drusen, or ≥656 μm in the absence of soft drusen, or at least one large druse [≥125 μm], or non-central geographical atrophy), as previously described for this cohort [40].

The unaffected siblings had normal maculae and were aged at least 65 years or older at the time at which the index patient was first diagnosed with neovascular AMD. These criteria were based on published epidemiological studies that indicated that elderly individuals with such maculae rarely go on to develop neovascular AMD during a ten-year follow-up [41]. Unaffected maculae fulfilled the following criteria: 0-5 small drusen (all less than 63 μm in diameter), no pigment abnormalities, no geographical atrophy and no neovascularisation (as defined previously; AMD 'category 1 or less' on the AREDS scale[39]).

Replication cohorts

Replication of significant findings was performed on three unrelated case-control cohorts. The central Greece cohort (as previously described[40]) included patients without AMD (controls; n = 198), with early and intermediate dry AMD (AREDS category 2 and 3; n = 84) and with neo-vascular AMD (n = 139) (Table S1 (Table 5)). These patients were recruited from the medical retina outpatient clinic at the University Hospital of Larissa, Greece. The diagnosis of macular degeneration was confirmed by optical coherence tomography and fluorescein angiography. Colour fundus photography and indocyanine green angiography were performed in some cases.

Table S1

Subject characteristics

Population	Mean age	Range	Standard deviation	% Male
Massachusetts Eye and Ear Infirmary Extremely Discordant Sibling Pair Discovery Cohort (n = 270 individuals from 135 families)
Neovascular subjects (n = 135)	71.7	49.0-85.5	8.0	47%
Control subjects (n = 135)	76.0	50.3-93.9	7.6	39%
Massachusetts Eye and Ear Infirmary Cohort Extended Sibling Pair Cohort (n = 657 individuals from 284 families)
Population	Mean age	Range	Standard deviation	% Male
Neovascular subjects (n = 353)	73.9	49.0-94.0	7.8	41%
Dry subjects (n = 106)	77.6	58.2-100.6	7.2	35%
Control subjects (n = 198)	75.4	50.3-94.3	8.3	44%
Central Greece cohort (n = 344 individuals)
Population	Mean age	Range	Standard deviation	% Male
Neovascular subjects (n = 139)	76.23	49.0-94.0	7.5	45%
Dry subjects (n = 84)	74.52	52.0-91.0	7.8	44%
Control subjects (n = 121)	73.53	48.0-88.0	7.3	51%
Prospective-based cohort derived from the Nurses' Health Study/Health Professionals Study (n = 1,527 individuals)
Population	Mean age	Range	Standard deviation	% Male
Neovascular subjects (n = 164)	61.1	42.7-74.0	6.0	46%
Dry subjects (n = 293)	59.5	41.3-74.0	5.7	18%
Control subjects (n = 1070)	60.2	41.3-74.0	5.9	36%

A second replication cohort was a prospective, nested case-control sample of participants from two cohorts: the NHS and HPFS. Details of these studies and ascertainment of AMD cases have been previously described [40,42,43]. This combined NHS/HPFS AMD cohort contained 1,527 individuals, including 293 with the early or dry form of AMD (AREDS category 2 and 3) and associated vision loss to 20/30 or worse, 164 with the neovascular form of the disease (AREDS category 4) and 1,070 age-matched controls who remained free of AMD during the follow-up period (Table S1 (Table 5)).

Epidemiological data

Data were obtained by conducting telephone or in-person interviews to complete a standardised questionnaire. Data included iris colour as blue (light), green (medium) or brown (dark). The regular use of multivitamins or vitamin D supplements was ascertained if supplements were taken at least twice a week for at least six months prior to the index patient's diagnosis of AMD. Smoking data, as measured in pack years (as described previously[37]), was available for every participant. Briefly, a smoker was defined as smoking ≥10 pack-years and a non-smoker as smoking less than 10 pack-years in their lifetime.

To quantify sunlight exposure, average UV radiation (UVR) exposure was calculated for the location at which each subject lived for the majority of his or her lifetime on 133 sibpairs, for a total of 276 individuals. Latitude, longitude and elevation for each city and state were determined using the website http://www.topozone.com. Using the calculated longitude (range, -123.03-63.19), latitude (range, 25.99-46.49) and elevation (range, 0.001-1.598 km) for each participant's city and state, UVR exposure was obtained using the tropospheric UV-visible (TUV) calculator found at the website http://cprm.acd.ucar.edu/Models/TUV/Interactive_TUV/, using input option 1: latitude and longitude, and output 2: irradiance, weighted. The TUV calculator uses a model that corrects for the effects of cloud cover, based on data collected from November 1978 to June 2000, as described on the website. UVR index and UVA and UVB irradiances were reported in watts/m2 for noon on the 15th of each month for the year 2006 in order to assess UV exposure throughout the year. Additionally, UV measures were taken on two specific dates: summer solstice and winter solstice, as these two days may have different UV exposures, being the longest and shortest days of the year, respectively. Yearly UV exposure was calculated as the average of all the recorded exposures.

To validate our findings of association between UVexposure and AMD, we obtained additional UV data on our family-based cohort from a separate source. The monthly average UV index was determined using the website for the US Environmental Protection Agency's SunWise Program http://www.epa.gov/sunwise1/doc/uvimonth.html. The website reports average monthly exposure in terms of UV index, based on data from the years 1994-2008. Yearly exposure was calculated by averaging monthly exposures obtained for the location at which each subject lived for the majority of their lifetime.

Self-reported sun exposure was also obtained, to account for cumulative sun exposure. Self-reported sun exposure was recorded for both current and --separately -- 15 years ago, as one of four categories: (i) spends most of the day outside; (ii) spends most free time outside; (iii) occasionally spends a day in the sun; or (iv) almost never outside.

Measurement of circulating vitamin D

For the collection of sera, non-fasting blood was drawn from 50 sibling pairs (100 subjects) in the discovery cohort of extremely discordant sibling pairs. Blood was collected into BD Vacutainer SST tubes (BD, Franklin Lakes, NJ, USA) and samples were centrifuged at room temperature (approximately 20°C) at 1500 g for 10 minutes. Serum was immediately drawn off and stored at -80°C in 2.0 ml cryogenic tubes (Corning Life Sciences, Lowell, MA, USA).

Levels of circulating 25(OH)D were measured in ng/ml, according to previously described methods [44]. Patients were classified as deficient if serum 25(OH)D was less than 20 ng/ml; as insufficient if levels were between 21 and 29 ng/ml; and as sufficient if levels were between 30 and 150 ng/ml [24]. Serum 25(OH)D levels greater than 150 ng/ml were classified as toxic.

RPE-choroid expression profile assessment

Whole-transcriptome expression profiles were available from 126 RPE-choroid and 118 retina punches (each 6 mm in diameter) collected from the macular and extramacular regions of eyes derived from 66 human donors. These eyes were selected from a well-characterised repository derived from 3,903 donors and collected by G.S.H. over a 20-year period at the University of Iowa and St Louis University. Medical and ophthalmic histories, a family questionnaire and blood and sera were obtained from the majority of donors. Gross pathological features, as well as the corresponding fundus photographs and angiograms (when available) of all eyes in this repository were read and classified by retinal specialists. Fundi and/or posterior poles were graded using a slightly modified version of two standardised classification systems, as published previously [4,5,45-47]. The ages of the donors ranged from nine to 101 years; approximately 50 per cent had documented clinical histories of AMD. RNA expression profiles were assessed using two-colour, 44 K Agilent whole-genome in situ oligonucleotide microarray analysis and a universal reference RNA experimental design. The universal reference RNA consisted of a 1:1 pool of RPE-choroid and retinal RNA generated from donors with and without AMD. After correcting for dye effects using Lowess normalisation, the net intensity values were determined and expressed as a percentage of the total array intensity. The ratios of the experimental and reference RNA signals were then calculated and the normalised percentage total of each experimental value was calculated by multiplication, using the geometric mean of all determinations of each probe's reference RNA value. For those probes with replicates in the array, the average values were determined. Inter-array differences were further corrected by quantile normalisation and probes that did not have net intensity values greater than six times the standard deviation of the background in at least 5 per cent of the samples were omitted. This resulted in a final dataset comprising 28,127 unique probes. Expression of the vitamin D pathway-associated genes VDR, CYP27A1, CYP27B1 and CYP24A1 was examined.

Investigation of vitamin D in human retinal cells

To investigate vitamin D signalling in retinal cells, vitamin D was applied to retinoblastoma cells as a first iteration. Expression levels of mRNA were compared between untreated and treated cells by quantitative real-time polymerase chain reaction (qRT-PCR) analysis. Weri-RB1 (American Type Culture Collection [ATCC] catalogue number HTB 169) and Y79 (ATCC catalogue number HTB 18) were utilised for cell culture and subsequent treatment with vitamin D3. Cells were cultured in accordance with the manufacturer's specification/requirements. Cells were treated with 10 μM Vitamin D3 (final concentration) for 1 hour, 6 hours, 12 hours, 24 hours or 48 hours. Following the allotted treatment time, cells were then harvested. RNA isolated from these cells was used for qRT-PCR assays. For each gene tested, three technical replicates were used, with 1 μl of cDNA being used in each well, with 5 μl of SYBR® Green master mix (Invitrogen Corporation Carlsbad, CA, USA), 3 μl dH20 and 1 μl of a forward/reverse mix of primers for a final primer concentration of 1 μM. The final reaction volume for each sample was 10 μl. Samples were run using an ABI 9600 platform (Life Technologies Corp, Carlsbad, CA, USA), for a total of 40 cycles, with readings being taken every cycle. Triplicates were averaged and the delta-Ct, compared with β-actin, was determined for each gene. Each triplicate reaction was performed for three separate RNA samples. Standard deviation was calculated for each dataset. Comparisons between time points and controls were made for fold-change estimation and relative expression, using the following formula: 1000/2Delta-CtBactin-Delta-Cttestgene.

Genotyping analysis

Genotyping was performed on leukocyte DNA that was purified by using standard phenol-chloroform or DNAzol (Invitrogen) extraction protocols. Initially, genotyping was performed with tagging SNPs (tagSNPs) from the vitamin D pathway genes CYP27A1, CYP27B1, VDR and CYP24A1. These particular genes were chosen because they are in regions of the genome previously reported by linkage studies to harbour AMD susceptibility loci [36,48]. CYP27A1 is located on chromosome 2q33-qter and spans approximately 34 kilobases (kb). CYP27B1 is located on chromosome 12q13.1-q13.3 and spans approximately 5 kb. VDR is located on chromosome 12q13.11 and spans 63 kb. CYP24A1 is located on chromosome 20q13 and spans approximately 21 kb. To assess variation within these genes, tagSNPs were chosen to span each gene using data from HapMap http://www.hapmap.org/ and the following criteria: (i) the minor allele frequency was greater than 10 per cent; (ii) the r2 value was at least 0.8; and (iii) SNPs tagged for at least six other SNPs. The location of SNPs with respect to introns/exons was analysed in terms of the largest coding transcript for each gene, as listed in Ensembl http://www.ensembl.org/index.html (Table S2 (Table 6)). These SNPs were supplemented with additional SNPs, chosen to be evenly spaced along the length of the gene in an effort fully to capture variation within a given gene. In total, 59 SNPs were initially screened by Sequenom (San Diego, CA, USA) in the extremely discordant sibling pair discovery cohort (135 sibling pairs, 270 individuals): CYP27A1, n = 4; VDR, n = 33; CYP27B1, n = 2; CYP24A1, n = 20 (Table S2 (Table 6)).

Table S2

Genomic location of SNPs†

Gene	SNP	Exon/intron‡	Position (bp)
CYP27A1	rs11677711*	Intron 1	chr2: 219,669,375
	rs7566656*	Intron 1	chr2: 219,664,519
	rs7594289*	Intron 1	chr2: 219,652,286
	rs933994*	Intron 1	chr2: 219,650,616
CYP27B1	rs8176345*	Exon 5	chr12: 58,158,558
	rs4646536*	Intron 6	chr12: 58,157,988
VDR	rs4516035*	Promoter	chr12: 48,299,826
	rs7299460*	Intron 1	chr12: 48,296,268
	rs10783219*	Intron 1	chr12: 48,295,488
	rs7136534*	Intron 1	chr12: 48,294,626
	rs3890733*	Intron 2	chr12: 48,289,373
	rs11574026*	Intron 2	chr12: 48,288,246
	rs11574027*	Intron 2	chr12: 48,287,373
	rs4237855*	Intron 2	chr12: 48,287,203
	rs4328262*	Intron 2	chr12: 48,285,648
	rs11168287*	Intron 2	chr12: 48,285,414
	rs4760648*	Intron 2	chr12: 48,280,665
	rs2853564*	Intron 2	chr12: 48,278,487
	rs2238136*	Intron 2	chr12: 48,277,713
	rs2254210*	Intron 3	chr12: 48,273,714
	rs2228570§,*	Exon 4	chr12: 48,272,895
	rs10875693*	Intron 4	chr12: 48,269,650
	rs2239186*	Intron 4	chr12: 48,269,410
	rs3782905*	Intron 4	chr12: 48,266,167
	rs3819545*	Intron 4	chr12: 48,265,006
	rs2189480*	Intron 4	chr12: 48,263,828
	rs12721370*	Intron 4	chr12: 48,262,073
	rs12717991*	Intron 4	chr12: 48,259,126
	rs2239179*	Intron 5	chr12: 48,257,766
	rs2283342*	Intron 5	chr12: 48,255,859
	rs2107301*	Intron 5	chr12: 48,255,570
	rs2239182*	Intron 5	chr12: 48,255,411
	rs11574077*	Intron 5	chr12: 48,252,927
	rs11168267*	Intron 5	chr12: 48,251,542
	rs7975128*	Intron 8	chr12: 48,245,828
	rs1544410*	Intron 10	chr12: 48,239,835
	rs7975232*	Intron 10	chr12: 48,238,837
	rs731236*	Exon 11	chr12: 48,238,757
	rs3847987*	Exon 11	chr12: 48,238,068
CYP24A1	rs2248359*	Promoter	chr20: 52791518
	rs6022999*	Intron 3	chr20: 52788013
	rs2585428*	Intron 3	chr20: 52786897
	rs2245153*	Intron 3	chr20: 52786406
	rs2296241	Exon 4	chr20: 52786219
	rs4809960*	Intron 4	chr20: 52786073
	rs4809959*	Intron 4	chr20: 52785859
	rs2181874*	Intron 4	chr20: 52784478
	rs3787557*	Intron 4	chr20: 52783135
	rs3787555*	Intron 4	chr20: 52782691
	rs3787554*	Intron 4	chr20: 52782680
	rs4809958*	Intron 4	chr20: 52782438
	rs2762939*	Intron 5	chr20: 52781251
	rs6068816*	Exon 6	chr20: 52781091
	rs6022993	Intron 7	chr20: 52779247
	rs6091828	Intron 7	chr20: 52779210
	rs6127119	Intron 7	chr20: 52779082
	c.990 + 189 G > A	Intron 7	chr20: 52779066
	rs6127118*	Intron 7	chr20: 52779052
	c.990 + 222_223ins99	Intron 7	chr20: 52779033
	rs2209314*	Intron 7	chr20: 52778961
	rs2296239	Exon 8	chr20: 52775528
	rs2274130	Intron 9	chr20: 52774601
	rs1570670	Intron 9	chr20: 52774579
	rs1977297	Intron 9	chr20: 52774479
	rs1570669*	Intron 9	chr20: 52774427
	rs927650*	Intron 11	chr20: 52772741
	rs2762934*	3'UTR	chr20: 52771261
	rs4809957	3'UTR	chr20: 52771171
	rs6022987	3'UTR	chr20: 52770596
	rs11907350	3'UTR	chr 20: 52770439
	rs10623012	3'UTR	chr20: 52770055-52770054
	rs4811494	3'UTR	chr20: 52769984
	rs2762933	3'UTR	chr20: 52769903
	rs6097807*	3'UTR	chr20: 52769455
	rs6068810*	3'UTR	chr20: 52769351

Abbreviations: chr, chromosomes; UTR, untranslated region; SNP, single nucleotide polymorphism.

†SNP chromosomal location according to Ensembl.

‡Exon/intron description based on largest transcript listed in Ensembl.

§rs2228570 was formerly rs10735810.

*SNPs investigated in the initial discovery cohort.

Multiplex PCR assays were designed using Sequenom SpectroDESIGNER software (version 3.0.0.3) by inputting sequence containing the SNP site and 100 base pairs (bp) of flanking sequence on either side of the SNP. Briefly, 10 ng of genomic DNA was amplified in a 5 μl reaction containing 1 × HotStar Taq PCR buffer (Qiagen, Valencia, CA, USA), 1.625 mM MgCl2, 500 μl each dNTP, 100 nM each PCR primer and 0.5 U HotStar Taq (Qiagen). The reaction was incubated at 94°C for 15 minutes, followed by 45 cycles at 94°C for 20 seconds, 56°C for 30 seconds and 72°C for 1 minute, followed by another 3 minutes at 72°C. Excess dNTPs were then removed from the reaction mixture by incubation with 0.3 U shrimp alkaline phosphatase (USB Corporation, Cleveland, OH, USA) at 37°C for 40 minutes, followed by 5 minutes at 85°C to deactivate the enzyme. Single primer extension over the SNP was carried out in a final concentration of between 0.625 μM and 1.5 μM for each extension primer (depending on the mass of the probe), iPLEX termination mix (Sequenom) and 1.35 U iPLEX enzyme (Sequenom) and cycled using a two-step 200 short cycles programme; 94°C for 30 seconds, followed by 40 cycles at 94°C for 5 seconds, five cycles at 52°C for 5 seconds, 80°C for 5 seconds and then 72°C for 3 minutes. The reaction was then desalted by the addition of 6 mg cation exchange resin, followed by mixing and centrifugation to settle the contents of the tube. The extension product was then spotted onto a 384-well SpectroCHIP before being flown in the matrix-assisted laser desorption/ionisation time-of-flight (MALDI-TOF) mass spectrometer. Data were collected, in real time, using SpectroTYPER Analyzer 3.3.0.15, SpectraAQUIRE 3.3.1.1 and SpectroCALLER 3.3.0.14 (Sequenom). Additionally, to ensure data quality, genotypes for each subject were checked manually.

In order to discover novel variation, a combination of direct sequencing and NextGen exon capture were performed. Oligonucleotide primers were selected using the Primer3 program http://primer3.sourceforge.net/ to encompass the SNP and flanking sequences (Table S3). All PCR assays were performed using genomic DNA fragments from 20 ng of leukocyte DNA in a solution of 10× PCR buffer containing 25 mM MgCl2, 0.2 mM each of dATP, dTTP, dGTP and dCTP, and 0.5 units of Taq DNA polymerase (USB). 5 M betaine was added to the reaction mix for rs2414687 (Sigma-Aldrich, St Louis, MO, USA). The temperatures used during the PCR were as follows: 95°C for 5 minutes followed by 35 cycles of 58°C for 30 seconds, 72°C for 30 seconds and 95°C for 30 seconds, with a final annealing at 58°C for 1.5 minutes and extension of 72°C for 5 minutes. For sequencing reactions, PCR products were digested according to the manufacturer's protocol using ExoSAP-IT (USB) and were then subjected to a cycle sequencing reaction using the Big Dye Terminator v3.1 Cycle Sequencing kit (Applied Biosystems, Foster City, CA, USA), according to the manufacturer's protocol. Products were purified using Performa DTR Ultra 96-well plates (Edge Biosystems, Gaithersburg, MD, USA) in order to remove excess dye terminators. NextGen sequencing was performed on exonic regions using methods that have been previously described [49]. For direct sequencing, samples were sequenced on an ABI Prism 3100 DNA sequencer (Applied Biosystems). Electropherograms generated from the ABI Prism 3100 were analysed using the Lasergene DNA and protein analysis software (DNASTAR, Inc., Madison, WI, USA). Electropherograms were read by two independent evaluators without knowledge of the subject's disease status. All patients were sequenced in the forward direction (5'to 3'), unless variants or polymorphisms were identified, in which case confirmation was obtained in some cases by sequencing in the reverse direction.

For the genotyping of the variation discovered by direct and NextGen sequencing, a combination of direct sequencing and Taqman genotyping was employed to supplement data previously acquired by Sequenom or to investigate further variation. Direct sequencing was used to genotype the rs6022993, rs6091828, rs6127119, c.990 +189 G > A, rs6127118 and c.990 +222_223ins99 SNPs using a single set of primers: forward primer TGGACTGAGGTTTTCTTCTGTG and reverse primer GGCACATGAAGTCTCCTCCT. For all other SNPs, Taqman genotyping was employed. Briefly, either TaqMan Pre-Designed SNP Genotyping Assays or Custom TaqMan SNP Genotyping Assays kits (Applied Biosystems) were designed (for listing of SNPs and probes, see Table S4). The 40× stock of the probes were diluted to 16× with 0.5× trisethylenediamine tetra-acetic acid (EDTA) and stored at -20°C. The amplification reaction was carried out in a total reaction volume of 16.25 μl containing 2.5 μl DNA (10 ng), 1.25 μl of 16× probe and 12.5 μl of TaqMan Genotyping Master Mix. Sample DNA was amplified using the following reaction: 1 minute at 60°C, 10 minutes at 95°C and 40 cycles of 15 seconds at 92°C and 1 minute at 60°C. The amplification reaction was carried out on thermo-cyclers and then fluorescence was measured on the ABI 7500 Real-Time PCR System by which the genotypes are analysed with the accompanying software or, in some cases, manually.

Statistical analysis

Epidemiological data

Conditional logistic regression (CLR) (SAS 9.0, SAS Institute Inc., Cary, NC, USA) was performed to identify factors associated with neovascular AMD. Potential risk factors of interest, as defined above, were evaluated initially in single factor analysis. Numerical variables such as serum 25(OH)D levels and UV irradiance were tested as continuous variables, while other variables--including iris colour, skin cancer history, self-reported sun exposure and vitamin supplementation--were treated as categorical. In addition to smoking, known genetic risk variants in both CFH and ARMS2/HTRA1 were controlled for [38,50-52]. For these genetic variables, each SNP (n = 18) was tested under an additive genetic model, along with UV irradiance. The backward selection procedure in conditional logistic regression was used to create the most parsimonious statistical model predictive of neovascular AMD. Specifically, each variable was tested for its association with AMD and entered into the model if its p-value was less than 0.15. Those variables that met these criteria were then kept in the model if the covariate-adjusted p-value was less than 0.1. Using the variables selected in backward selection, a multiple variable model was constructed and tested for its association with neovascular AMD.

Genotyping data

SNPs were tested for association using the minor allele, as defined by the allele occurring less frequently in the normal subjects for all cohorts examined. Initial testing of association between the variation and AMD in the extremely discordant sibling pair was done using the family-based association test (FBAT) http://biosun1.harvard.edu/~fbat/fbat.htm and generalised estimating equations (GEE), assuming an additive genetic model in SAS. For validation in the extended sibling cohort, allelic associations were performed in SAS (v9.1) using conditional logistic regression (CLR). For replication in the unrelated case-control cohorts, allelic associations were performed using unconditional logistic regression (UCLR) in SAS. Both analyses were performed assuming an additive genetic model and controlling for covariates (age, sex and smoking). All haplotypic associations were performed initially using FBAT in the discovery cohort and then for validation and replication using the program UNPHASED http://www.mrc-bsu.cam.ac.uk/personal/frank/software/unphased/. For the family data, the model for sibships was used. UNPHASED uses likelihood ratio association analysis, and for sibships uses nuclear families with missing parental genotype data [53]. Covariates--specifically age, sex and smoking history--were controlled for in analysis using residuals, as calculated in R Statistical Software http://www.r-project.org/.

Linkage disequilibrium (LD) (D' and r2) between each of the SNPs was determined using Haploview http://www.broad.mit.edu/mpg/hap-loview/. Deviation from Hardy-Weinberg equilibrium was tested on each SNP using the chi-square test in the unaffected population. Haplotype blocks were constructed in Haploview using the method proposed by Gabriel et al [54]. Haplotype blocks defined by the Gabriel rule were tested for association with all subtypes of AMD using UNPHASED, and additionally FBAT in the families. In addition to the use of residuals, multiple testing in FBAT was corrected for using the permutation test. Haplotype testing was also performed by two additional methods, backward selection and sliding window. Backward selection was performed in UNPHASED, beginning with all SNPs with association p value ≤0.12 in the model and taken out one at a time in order to determine the most significantly parsimonious haplotype [55]. Sliding window analysis was also performed in UNPHASED, using a window size of 2 and using all SNPs with association p ≤ 0.12. Genotype and allele frequencies for all SNPs identified as either individually significant or significant as part of a haplotype were calculated in the affected patients and separately in unaffected controls (Table S4).

Based on the variation found within the extended sibling cohort, replication was performed on a cohort of unrelated cases and controls from central Greece and the combined NHS and HPFS cohort. Meta-analysis of all four cohorts was performed for SNPs genotyped in all four cohorts and also significant haplotypes using the Comprehensive Meta-Analysis v2 (Biostat, Englewood NJ, USA) program. Z-scores and overall p-values for each SNP were calculated under a fixed model based on odds ratios obtained from CLR for the family-based cohorts and UCLR for the case-control cohorts calculated in SAS while controlling for the covariates age, sex and smoking history. To test for the effects of population heterogeneity, Cochran's Q was also calculated, which is the weighted sum of squared differences between individual study effects and the pooled effect across studies.

In order to test for interaction or epistasis between the genes investigated here and previously reported AMD-associated genes, an interaction term was added into the regression model, including both variables. SNPs included in statistical analyses were SNPs that had been genotyped in more than one cohort: VDR (rs2189480); CFH (rs572515 and rs1061170); ARMS2/HTRA1 (rs10490924, rs11200638 and rs1049331); and RAR-related Orphan receptor A gene (RORA; rs12900948 and rs4335725) [38,40,42,51,52]. An additional pooled case-only analysis was performed in logistic regression to test for interaction with greater power. Statistical interaction was also tested for between UV measures and the variants examined using an interaction term in the regression models.

Bioinformatic analysis of functional implications of genetic variants

Variants were analysed by computational methods to infer potential functional relevance. Conservation analysis was performed using the comparative genomics feature in Ensembl http://useast.ensembl.org/index.html to compare the human variant with nine eutherian mammals, including chimp, mouse and pig. To analyse the putative biological effects of each exonic or intronic variation, either the major or minor allele was input into predictive programs flanked on either side by 100 bp of sequence. Potential human transcription factor binding sites were analysed using default parameters of the transcription element search system from the University of Pennsylvania http://www.cbil.upenn.edu/cgi-bin/tess/tess. Additional analysis of putative splice sites was performed using the default parameters of the Alternative Splice Site Predictor http://www.es.embnet.org/~mwang/assp.html. Also, the publicly available program mRNA by SNP Browser v 1.0.1[56,57] was used to investigate the association of significant CYP24A1 SNPs and expression quantitative trait loci (eQTLs).

Pathway analyses

We used Ingenuity Pathway Analysis (IPA) software http://www.ingenuity.com/ to create a visual representation of the complex shared pathways involved in the aetiology of AMD, the metabolism/synthesis of vitamin D and the vascularisation processes. The likelihood that each gene/protein is randomly associated with a given function, pathway or other protein/gene is automatically calculated using the right-tailed Fisher's exact test. We accepted only those relationships with p > 0.0001, suggesting that their association with our gene/protein functions or diseases was non-random.

Results

Patient population

For our initial discovery cohort, we recruited a total of 135 extremely discordant sibling pairs (one sibling with neovascular AMD and one unaffected sibling -- AREDS category 1 or less). The mean ± SD age of the unaffected siblings was 76.0 ± 7.6 years (age range, 50.3-93.9) and the mean ± SD of the affected (neovascular) siblings was 71.7 ± 8.0 years (age range, 49.0-86.5) (Table S1 (Table 5)). Of the affected patients, 47 per cent were male and 39 per cent of the matching unaffected siblings were male. Most participants were of European descent.

The expanded sibling cohort comprised 481 sibling pairs and 657 individuals: 198 unaffected controls, 106 early and intermediate dry AMD patients (AREDS category 2 and 3), and 353 neo-vascular AMD patients (Table S1 (Table 5)). The mean ± SD age of the unaffected subjects was 75.4 ± 8.3 years (age range, 50.3-94.3 years), the mean age of the patients with early or intermediate AMD was 77.6 ± 7.2 years (age range, 58.2-100.6 years), and the age of those subjects with neovascular AMD was 73.9 ± 7.8 years (age range, 49.0-94.0 years). Forty-four percent of the unaffected subjects, 35 per cent of the AREDS category 2 and 3 AMD subjects and 41 per cent of the neovascular subjects were male. The cohort from central Greece consisted of a total of 344 patients: 121 unaffected controls, 84 early and intermediate dry AMD patients (AREDS category 2 and 3) and 139 neovascular AMD patients (Table S1 (Table 5)). The mean ± SD age of the unaffected subjects was 73.5 ± 7.3 years (age range, 48.0-88.0 years), the age of the patients with early or intermediate AMD was 74.5 ± 7.8 years (age range, 52.0-91.0 years) and the mean age of those subjects with neovascular AMD was 76.2 ± 7.5 years (age range, 49.0-94.0 years). Fifty-one per cent of the unaffected subjects, 44 per cent of the AREDS category 2 and 3 AMD subjects and 45 per cent of the neovascular subjects were male. All participants were white and from central Greece.

The prospective cohorts from the NHS and HPFS consisted of a total of 1,527 patients: 1,070 unaffected controls, 293 early and intermediate dry AMD patients (AREDS category 2 and 3) and 164 AMD neovascular patients. The mean ± SD age of the unaffected subjects was 60.2 ± 5.9 years (age range, 41.3-74.0 years), the age of the patients with early or intermediate AMD was 59.5 ± 5.7 years (age range, 41.3-74.0 years) and the mean age of those subjects with neovascular AMD was 61.1 ± 6.0 years (age range, 42.7-74.0 years). Thirty-six per cent of the unaffected subjects, 18 per cent of the AREDS category 2 and 3 AMD subjects and 46 per cent of the neovascular subjects were male.

Analysis of sunlight-associated epidemiological factors

The analysis of epidemiological data was performed on 133 extremely discordant sibling pairs, equalling 266 subjects. Twenty-six affected patients (19.5 per cent) had skin cancer and 25 (18.8 per cent) of their matched unaffected siblings also had skin cancer. Iris colour in the affected patients was 42.9 per cent blue (n = 151), 24.8 per cent green (n = 33) and 32.3 per cent brown (n = 34). In the matched unaffected patients iris colour was 41.6 per cent blue (n = 132), 25.0 per cent green (n = 33) and 33.4 per cent brown (n = 45). Of the affected siblings, 94 (70.6 per cent) reported taking multivitamin supplements at least twice a week for at least six months and 76 (57.1 per cent) of the unaffected patients reported the same. Specific vitamin D supplementation was reported in 11 affected patients (8.3 per cent) and 15 unaffected patients (11.3 per cent).

The average yearly UV index in the affected patients was 5.24 W/m2 (range, 4.42-7.49) and 5.43 W/m2 (range, 4.37-7.49) in the unaffected patients. Average yearly UVA in the affected patients was 40.22 W/m2 (range, 36.16-50.58) and 41.17 W/m2 (range, 36.17-50.57). Average yearly UVB in the affected patients was 0.99 W/m2 (range, 0.84-1.41) and 1.03 W/m2 (range, 0.84-1.41) in the unaffected patients.

Single CLR regression showed that, of these factors, UV, UVA and UVB irradiance were the risk factors most significantly associated with neo-vascular AMD (Table 1) (p = 0.0013 for all). Analysis of self-reported sun exposure suggested a protective effect for those who reported minimal sun exposure 15 years in the past (p = 0.04); however, this did not remain significant after correction for multiple testing. Skin cancer status, iris colour, vitamin D supplementation and various exposures of sunlight were not associated with neo-vascular AMD risk in CLR analysis.

Table 1

Single factor analysis of epidemiological factors using conditional logistic regression

			Distribution of continuous variables (if applicable)
Risk factor	Odds ratio (95% confidence interval)	p-Value	Average	Min	Max
Vitamin D supplements for ≥6 months	0.862 (0.505-1.472)	0.5866
Multivitamin supplements for ≥6 months	0.667 (0.237-1.873)	0.4417
Serum 25(OH)D level (ng/ml)	0.973 (0.934-1.013)	0.1859	26.10	5.10	67.30
Blue eyes	1.064 (0.654-1.731)	0.8039
Green eyes	1.000 (0.573-1.744)	1.0000
Brown eyes	0.966 (0.579-1.613)	0.8960
Skin cancer	1.050 (0.570-1.933)	0.8762
Sun exposure now 1†	1.071 (0.517-2.220)	0.8528
Sun exposure now 2†	1.350 (0.757-2.407)	0.3090
Sun exposure now 3†	0.929 (0.545-1.584)	0.7856
Sun exposure now 4†	0.750 (0.407-1.382)	0.3562
Sun exposure now 1 or 2†	1.381 (0.788-2.421)	0.2600
Sun exposure 15 years ago 1†	0.955 (0.525-1.736)	0.8788
Sun exposure 15 years ago 2†	1.333 (0.831-2.141)	0.2336
Sun exposure 15 years ago 3†	1.045 (0.583-1.876)	0.8815
Sun exposure 15 years ago 4†	0.375 (0.147-0.958)	0.0405
Sun exposure 15 years ago 1 or 2†	1.409 (0.816-2.433)	0.2186
Year UV index (TUV)‡	0.115 (0.031-0.432)	0.0013	5.33	4.37	7.49
Year UV-A (TUV)‡	0.619 (0.463-0.829)	0.0013	40.69	36.16	50.58
Year UV-B (TUV)‡	0.009 (0.008-0.010)	0.0013	1.01	0.84	1.41
Year UV index (EPA)§	0.104 (0.020-0.541)	0.0071	4.13	1.00	10.00

Abbreviations: EPA, Environmental Protection Agency; Min, minimum; Max, maximum; 25(OH)D, 25-hydroxyvitamin D.

†Self-reported sun exposure 1: spends most of the day outside; self-reported sun exposure 2: spends most of free time outside; self-reported sun exposure 3: occasionally spends a day in the sun; self-reported sun exposure 4: almost never outside.

‡As calculated using data from the TUV calculator http://cprm.acd.ucar.edu/Models/TUV/Interactive_TUV/.

§As calculated using data from the US EPA SunWise Program http://www.epa.gov/sunwise1/doc/uvimonth.html.

The backward selection model using previously reported SNPs (12 in CFH and six in ARMS2/HTRA1), smoking history and yearly UV irradiance showed that the model most predictive of neovascular AMD included smoking greater than ten pack-years, CFH rs572515, HTRA1 rs1049331, and total UV index, UVA irradiance or UVB irradiance (data not shown). Each UV irradiance measure was shown to be significant in the multiple model (p = 0.01 for all). Most significantly, for every one-unit increase in UV index between 4.73 and 7.49, AMD risk was reduced 66.6-fold after controlling for smoking, CFH rs572515 and HTRA1 rs1049331 (Table 2).

Table 2

Multiple conditional logistic regression for UV irradiance

Risk factor	Odds ratio (95% CI)	p-Value
Year UV index	0.015 (0.001-0.433)	0.0143
Smoking ≥10 pack-years	3.581 (1.590-8.063)	0.0021
CFH rs572515 (additive model)	7.971 (2.960-21.459)	3.99E-05
HTRA1 rs1049331 (additive model)	9.829 (3.589-26.920)	8.77E-06
Risk factor	Odds ratio (95% CI)	p-Value
Year UV-A	0.368 (0.163-0.828)	0.0157
Smoking ≥10 pack-years	3.705 (1.631-8.418)	0.0018
CFH rs572515 (additive model)	8.403 (3.057-23.098)	3.69E-05
HTRA1 rs1049331 (additive model)	9.969 (3.618-27.467)	8.72E-06
Risk factor	Odds ratio (95% CI)	p-Value
Year UV-B	0.009 (0.000-0.011)	0.0153
Smoking ≥ 10 pack-years	3.606 (1.597-8.141)	0.0020
CFH rs572515 (additive model)	8.084 (2.985-21.895)	3.93E-05
HTRA1 rs1049331 (additive model)	9.869 (3.598-27.073)	8.73E-06

Abbreviation: CI, confidence interval.

Analysis of serum levels of vitamin D

Serum 25(OH)D levels were measured for 100 patients comprising 50 extremely discordant sibling pairs. Of the 50 unaffected siblings, 30 took multi-vitamins, five took both multivitamins and vitamin D supplements and one took vitamin D supplements alone. The mean ± SD serum 25(OH)D levels in affected patients was 24.7 ± 10.5 ng/ml (concentration range, 5.1-49.6 ng/ml) and in unaffected patients was 27.5 ± 10.7 ng/ml (concentration range, 8.9-67.3 ng/ml). Serum 25(OH)D levels were not statistically different when comparing affected and unaffected siblings (p = 0.22) (Figure 1), or when controlling for the month in which serum was collected.

Figure 1

Analysis of serum levels of vitamin D in extremely discordant sibpairs.

RPE-choroid expression analysis

Of the vitamin D pathway-associated genes examined in a whole transcriptome analysis of ocular tissues derived from 66 human donors, no significant association as a function of age or disease status was observed; however, differences in expression levels between the macular and extramacular regions were observed. Expression of both VDR and CYP27A1 was detected in the RPE-choroid; CYP27A1 and CYP27B1 were detected in the retina. CYP27A1 was expressed at low to moderate levels, with significant differences detected between the extramacula and macula regions (Figure 2A). Similar differences in the extramacula:macula expression ratio of CYP27A1 were observed in both non-AMD (1.5-fold, p = 0.00012; one-tailed Mann-Whitney U test) and AMD (1.7-fold, p = 5 × 10-8) RPE-choroids. Similar macula:extramacula expression differences were observed for CYP27A1 (Figure 2C) and CYP27B1 (Figure 2D) in the retina. VDR expression in many RPE-choroid tissues was at or below the level of detection, but was observed at significantly higher levels in several donors, concomitant with an increased expression of genes involved in local inflammatory responses (M.J. Radeke, D.H. Anderson and L.V. Johnson, personal communications). There was a trend of slightly higher VDR expression in the macular region of both normal and diseased RPE-choroids (1.3- and 1.5-fold, respectively) which rose to a level of significance (p = 0.02) in the diseased RPE-choroid (Figure 2B), but not the retina (data not shown).

Figure 2

RPE-choroid expression analysis.

Investigation of vitamin D in retinal cells

Both cell lines were derived from a human retino-blastoma; however, they were isolated and deemed continuous at two separate institutions. As both cell lines were derived from a retinoblastoma, and with previous studies implicating vitamin D3 as a potential therapeutic agent for this cancer,[58] a subset of genes involved in the vitamin D3 pathway were selected for testing via qRT-PCR. Vitamin D treatment of retinoblastoma cells did not result in significant changes in mRNA expression of any of the genes investigated. The genes tested by qRT-PCR were CYP24A1 and the calreticulin gene (CALR). For these two genes, vitamin D3 treatment resulted in gene expression changes of varying degrees in both cell lines. A general increasing trend was noted for CYP24A1, with peak expression occurring at the 48-hour time point in Weri-RB1-RB1-RB1 cells. By contrast, the Y79 cells exhibited a similar increasing trend, with peak expression at six hours, but then followed by a drastic decrease which remained until the final time point at 48 hours. CALR exhibited a pattern of an initial increase in expression, followed by a subsequent decrease and then a gradual increase over the subsequent time points. In Weri-RB1 cells, this was accompanied by peak expression at the final time point, while in Y79 cells the subsequent decrease occurred later, with no subsequent increase being noted owing to the decrease in expression occurring at the final time point at 48 hours. It would have been interesting to see if expression began to increase at later time points, as it did in Weri cells (see Figure 3).

Figure 3

Vitamin D treatment.

Analysis of CYP27A1, CYP27B1, VDR and CYP24A1 SNPs in the extremely discordant sibling pair cohort: Discovery cohort

No significant deviations from Hardy-Weinberg equilibrium were observed for any of the variants studied in VDR, CYP27A1, CYP27B1 and CYP24A1 in unaffected siblings. Of the four genes tested, variation within both VDR and CYP24A1 were found to be significantly associated with neo-vascular AMD. In VDR, rs2189480 (intron 4) was found to be significantly associated with risk of neovascular AMD by FBAT (p = 0.02) (Table 3). In FBAT, the minor alleles of CYP24A1 rs6127118 (intron 7) and rs2762934 (exon 12) were found to be associated with risk of neovascular AMD (p = 0.03 and p = 0.01). Two haplotype blocks were shown to be significantly associated with neovascular AMD risk in FBAT (Figure 4). In block 4, comprising rs6068816 (exon 6) and rs6127118 (intron 7), haplotype h1 (GC) was modestly significant (p < 0.03) and haplotype h2 (GT) was significant under an additive model (p = 0.002). In block 5, comprising rs1570669 (intron 9), rs927650 (intron 11) and rs2762934 (exon 12), haplotype h3 (TGT) was significantly associated (p = 0.01). Due to the fact that these haplotypes were more significantly associated with AMD than any single SNP alone, and that there is more power to detect an association through haplotypes,[59] we used these regions to drill down in the replication and validation cohorts.

Table 3

Family-based association test analysis of individual CYP27A1, CYP27B1, VDR and CYP24A1 SNPs in discovery cohort

Marker	Allele	Frequency	Informative families	Z score	p-Value
CYP27A1
rs7594289	G	0.428	51	-0.808	0.4189
rs7566656	A	0.485	60	1.179	0.2383
CYP27B1
rs8176345	A	0.030	< 2	0.000	1.0000
rs4646536	C	0.324	< 2	0.000	1.0000
VDR
rs4516035	G	0.443	47	1.116	0.2643
rs7299460	A	0.262	36	1.333	0.1824
rs10783219	A	0.361	42	-0.630	0.5287
rs7136534	A	0.238	35	1.622	0.1048
rs11574026	T	0.120	23	0.626	0.5316
rs11574027	T	0.017	2	1.414	0.1573
rs4237855	T	0.490	47	-0.329	0.7419
rs4328262	A	0.568	49	0.956	0.3390
rs11168287	C	0.485	56	-0.570	0.5688
rs4760648	A	0.447	52	-0.640	0.5221
rs2853564	C	0.395	46	0.000	1.0000
rs2238136	A	0.283	43	-1.000	0.3173
rs2254210	T	0.369	54	0.843	0.3994
rs10735810	T	0.374	48	0.739	0.4602
rs10875693	A	0.679	46	0.000	1.0000
rs2239186	C	0.248	44	-0.137	0.8907
rs3782905	G	0.282	34	-0.152	0.8788
rs3819545	C	0.427	54	0.226	0.8208
rs2189480	C	0.593	45	-1.807	0.0707
rs12721370	T	0.112	15	1.414	0.1573
rs12717991	G	0.564	51	0.000	1.0000
rs2239179	G	0.349	47	-0.911	0.3621
rs2283342	C	0.205	38	-0.302	0.7630
rs2107301	T	0.340	46	0.405	0.6858
rs2239182	G	0.420	46	0.788	0.4308
rs11574077	G	0.028	10	-0.632	0.5271
rs11168267	T	0.116	17	0.728	0.4669
rs7975128	T	0.333	45	0.132	0.8946
rs1544410	A	0.337	46	-0.640	0.5221
rs7975232	T	0.494	52	0.239	0.8111
rs731236	C	0.320	45	-0.775	0.4386
rs3847987	T	0.150	28	1.029	0.3035
CYP24A1
rs2248359	A	0.425	49	-0.611	0.5413
rs6022999	C	0.250	42	-1.089	0.2763
rs2585428	A	0.441	56	-0.342	0.7324
rs2245153	G	0.201	48	1.089	0.2763
rs4809960	G	0.229	51	1.192	0.2332
rs4809959	C	0.500	46	0.122	0.9028
rs2181874	T	0.247	45	-1.633	0.1025
rs3787557	G	0.144	40	1.265	0.2059
rs3787555	T	0.263	42	1.260	0.2076
rs3787554	T	0.065	20	0.209	0.8348
rs4809958	C	0.128	25	1.134	0.2568
rs2762939	G	0.275	36	-0.816	0.4142
rs6068816	A	0.079	22	0.000	1.0000
rs6127118	T	0.157	22	2.132	0.0330
rs2209314	G	0.256	43	-1.387	0.1655
rs1570669	C	0.335	45	0.630	0.5287
rs927650	G	0.511	44	-0.391	0.6961
rs2762934	T	0.156	29	-2.535	0.0112
rs6097807	C	0.241	50	-1.131	0.2579
rs6068810	A	0.059	17	0.728	0.4669

Figure 4

Significant .

Sequencing of area surrounding significant variation

Based on the initial findings of significant variation within VDR and CYP24A1, a combination of direct sequencing and NextGen sequencing was performed on 95 neovascular patient samples in order to capture variation beyond the significantly associated tagging SNPs to include exons and surrounding exon/intron borders. Specifically, exons 4 and 5 of VDR and exons 4-12 of CYP24A1 were assayed for variation. No further variation was found within VDR, but the following 22 SNPs were found in the regions investigated in CYP24A1: rs2296241, rs4809960, rs4809959, rs2181874, rs6068816, rs6022993, rs6091828, rs6127119, rs6127118, rs2296239, rs2274130, rs1570670, rs1977297, rs1570669, rs927650, rs2762934, rs4809957, rs6022987, rs11907350, rs10623012, rs4811494 and rs2762933 (see Table S2 (Table 6)). Additionally, novel variations were found within CYP24A1 intron 7 at c.990 +189 (chromosome 20: 52,779,066; G > A) and c.990 + 222-223 (chromosome 20: 52,779,033; an insertion of 99 bp of the sequence GCGCACAGCATCTATACCACGGT GCCTCAGCGCACGGCCTCATACTAGACACCA CGGTGCCTCCGCGCACGGCCTCATACCAG ACGCCACGGTGCCTCC). Both variants were confirmed by sequencing in the reverse direction.

Further direct sequencing showed the 99 bp insert to be a repeated sequence. Patients were found to be both homozygous and heterozygous for this repeat. Interestingly, an alternative heterozygote was found that consisted of one 99 bp repeat and a second, smaller insert. This second insert had the same sequence as that described for rs3831062 but at a different location than that initially described (chromosome 20: 52,779,033; an insertion of 35 bp of the sequence GGAGGCACCGTGGCGTCTG GTATGAGGCCGTGCGC). Direct sequencing was performed for all subsequent variation of this sequence. This 99 bp insertion was not associated with AMD in either the sibling cohort or the Greek cohort.

Analysis of VDR and CYP24A1 SNPs in the replication cohorts

Based on the results of direct sequencing, 24 SNPs/structural variants were genotyped in the extended sibling cohort consisting of 657 individuals/481 sibling pairs and the unrelated case-control cohort from central Greece by a combination of Sequenom, Taqman and direct sequencing. Based on the preliminary findings of association in the discovery and expanded sibling cohorts, the following 12 SNPs were genotyped in the NHS/HPFS cohort: VDR rs2189480, CYP24A1 rs2181874, rs6127118, c.990 +222-223ins99, rs2296239, rs2274130, rs1570670, rs1570669, rs927650, rs2762934, rs4809957 and rs11907350.

Haplotype analysis of CYP24A1

Haploview, using haplotype blocks defined by the Gabriel rule, [54] demonstrated that in the extended sibling cohort there were four haplotype blocks spanning the more than 16,000 base pairs covered by genotyping (Figure 5). Of these haplotype blocks, none was shown to be significantly associated with any subtype of AMD risk by either FBAT analysis or UNPHASED analysis (data not shown). Similar results were found for the Greek cohort and the NHS/HPFS in UNPHASED (data not shown).

Figure 5

Linkage disequilibrium in the extended sibpair cohort. (A) shows linkage disequilibrium measures in D'and (B) shows linkage disequilibrium measures in r2

Those CYP24A1 SNPs that were shown to have a p-value ≤ 0.12 in the analysis of individual SNPs using UNPHASED were analysed in a backward selection fashion and, additionally, by using a sliding window approach for haplotype association using the same program. In each of the four cohorts, no method identified significant haplo-types (data not shown).

Meta-analysis of VDR and CYP24A1

Comprehensive meta-analysis v2 was used to analyse the variation which had been genotyped in all populations. VDR rs2189480 was not significantly associated with AMD. Several SNPs within CYP24A1 were shown to be significant by meta-analysis. When comparing all AMD with unaffected patients, rs1570669, rs1570670, rs2274130, rs2296239 and rs4809957 were all significantly associated with AMD (p = 0.0005, p = 0.02, p = 0.008, p = 0.02 and p = 0.03, respectively). For each of these SNPs, the overall p-value obtained by meta-analysis was more significant than that for each of the cohorts analysed where, in some cases, these SNPs were not individually significant (Table 4). Cochran's Q, the measure of heterogeneity, was significant for SNP rs927650 (p = 0.02). This was the SNP that the minor allele was defined as (A) and the Sibling cohort and the NHS/HPFUS, but was defined as (G) in the Greeks. Overall, the test for heterogeneity was not statistically significant (p > 0.1).

Table 4

Single SNP association results comparing all AMD subtypes with normal subjects

			Expanded sibpair cohort		Greeks		NHS/HPFS		Meta-analysis
SNP	Gene	Allele	Odds ratio(95% CI)	p-Value	Odds ratio(95% CI)	p-Value	Odds ratio(95% CI)	p-Value	Odds ratio(95% CI)	p-Value
rs2189480	VDR	A	1.801 (1.069-3.035)	0.0270	1.199 (0.846-1.700)	0.3080	0.903 (0.766-1.064)	0.2213	0.997 (0.865-1.151)	0.9710
rs2296241	CYP24A1	C	1.345 (0.829-2.182)	0.2294	1.185 (0.842-1.668)	0.3306			1.036 (0.814-1.319)	0.7711
rs4809960	CYP24A1	G	0.983 (0.567-1.703)	0.9515	1.027 (0.690-1.527)	0.8967			1.116 (0.843-1.477)	0.4425
rs4809959	CYP24A1	T	1.238 (0.770-1.989)	0.3781	1.058 (0.761-1.470)	0.7369			1.124 (0.885-1.427)	0.3391
rs2181874	CYP24A1	T	0.871 (0.501-1.515)	0.6242	0.687 (0.481-0.983)	0.0400	0.989 (0.812-1.204)	0.9116	0.894 (0.763-1.046)	0.1629
rs6068816	CYP24A1	A	2.081 (0.945-4.586)	0.0690	1.141 (0.695-1.875)	0.6012			1.243 (0.851-1.817)	0.2600
rs6091828	CYP24A1	A	1.338 (0.281-6.382)	0.7148	0.536 (0.206-1.393)	0.2006			0.745 (0.356-1.560)	0.4352
rs6127119	CYP24A1	A	1.284 (0.724-2.274)	0.3924	1.245 (0.831-1.865)	0.2881			1.257 (0.937-1.688)	0.1275
c.990 + 189G > A	CYP24A1	A	0.887 (0.138-5.700)	0.8993	1.699 (0.323-8.925)	0.5312			1.123 (0.347-3.630)	0.8464
rs6127118	CYP24A1	T	1.297 (0.725-2.319)	0.3816	1.274 (0.845-1.921)	0.2471	1.050 (0.864-1.276)	0.6230	1.111 (0.946-1.305)	0.1995
c.990 + 222-223 ins99	CYP24A1	ins	1.341 (0.744-2.417)	0.3288	1.297 (0.867-1.939)	0.2055	1.086 (0.874-1.350)	0.4564	1.145 (0.963-1.363)	0.1258
rs2296239	CYP24A1	A	1.457 (0.819-2.590)	0.2000	1.529 (0.976-2.396)	0.0638	1.205 (0.993-1.463)	0.0584	1.217 (1.035-1.430)	0.0174
rs2274130	CYP24A1	C	1.751 (0.975-3.147)	0.0610	1.512 (0.970-2.358)	0.0679	1.255 (1.029-1.529)	0.0247	1.248 (1.059-1.47)	0.0082
rs1570670	CYP24A1	C	1.790 (0.981-3.267)	0.0579	1.512 (0.970-2.358)	0.0679	1.195 (0.984-1.452)	0.0725	1.220 (1.037-1.436)	0.0165
rs1977297	CYP24A1	A	0.725 (0.293-1.791)	0.4856	1.142 (0.621-2.101)	0.6693			1.009 (0.643-1.582)	0.9694
rs1570669	CYP24A1	C	1.153 (0.690-1.927)	0.5875	1.287 (0.894-1.852)	0.1753	1.309 (1.104-1.551)	0.0019	1.285 (1.116-1.478)	0.0005
rs927650	CYP24A1	A	0.914 (0.596-1.402)	0.6804	1.468 (1.048-2.056)	0.0257	0.853 (0.728-0.999)	0.0487	0.969 (0.852-1.102)	0.6316
rs2762934	CYP24A1	T	0.743 (0.367-1.502)	0.4082	1.067 (0.704-1.619)	0.7595	0.930 (0.754-1.148)	0.5016	0.949 (0.798-1.130)	0.5574
rs4809957	CYP24A1	C	1.475 (0.822-2.644)	0.1925	1.462 (0.943-2.268)	0.0896	1.191 (0.981-1.445)	0.0772	1.196 (1.017-1.405)	0.0300
rs6022987	CYP24A1	G	0.755 (0.427-1.336)	0.3347	1.169 (0.804-1.699)	0.4136			1.109 (0.845-1.456)	0.4567
rs11907350	CYP24A1	T	0.363 (0.066-2.000)	0.2445	1.296 (0.438-3.836)	0.6397	1.246 (0.829-1.872)	0.2896	1.111 (0.772-1.598)	0.5711

Abbreviation: CI, confidence interval; SNP, single nucleotide polymorphism.

Bioinformatic analysis of potential function of variants

Of the SNPs found to be significantly associated with AMD, three SNPs were well conserved across species: rs2296239, rs2274130 and rs1570670, spanning from exon 8 to intron 9 of CYP24A1. Transcription Element Search System (TESS) analysis showed that, of our five significant SNPs, three SNPs changed or created transcription factor binding sites with the presence of the minor allele. Specifically, the minor allele of rs2274130 eliminated the binding site of B factor T00286 factor i T00557 NF-E2; rs1570670 changed the binding site from Gata-1 to GR; and the minor allele of rs1570669 changed the binding site from NF-BA1 to STE12. Using the publicly available eQTL database mRNA by SNP Browser v 1.0.1,[56,57] we investigated the association of the CYP24A1 SNPs that were significant alone or as part of a haplotype with eQTLs. Of our six SNPs, rs2181874 and rs1570669 were associated with the genes TSC22 domain family, member 1 (TSC22D1), unc-51-like kinase 2 (ULK2) and ring finger protein 213 (RNF213), respectively. These results suggest that further experimental analysis of these data would be worthwhile.

In order to identify causal variation, the creation of putative splice sites by genetic variants was analysed by inputting sequence into the alternative splice site predictor http://www.es.embnet.org/~mwang/assp.html (using the default parameters). Splice site analysis of the minor allele (G) of rs927650--the SNP most significantly associated with an increased risk of AMD in the Greek cohort--predicted the creation of an alternative isoform/cryptic donor at the bp immediately preceding this SNP (chromosome 20: 52772742) within intron 11 of the longest CYP24A1 transcript (CYP24A1-003) with a score of 10.967. Further studies should be conducted to test this prediction in a biological setting.

IPA analyses indicated that five genes (RXRA, VDR, PTH, CYP24 and CYP24A1) involved in the regulation of vitamin D metabolism, are connected by direct and indirect relationships to one AMD-associated gene (HTRA1) for which the pathway is unknown [10,38,52,60-62]. Additionally, some of these pathway genes directly or indirectly interact with genes involved in angiogenesis (the CALR, catenin [cadherin-associated protein] beta 1 (CTNNB1), bone morphogenetic protein 4 (BMP4) and interleukin 6 [IL6] genes), one of the hallmarks of neovascular AMD.

As the pathway map shows (Figure 6), some of the genes are implicated in both AMD and vitamin D pathways (VDR and CYP24A1), or both angiogenesis and AMD (BMP4). When the pathways linking the genes/proteins described above are grown to incorporate other related molecules, the retinoid X receptor-α gene (RXRA), encoding a ligand-dependent nuclear receptor which binds directly to VDR, is identified as a potentially relevant gene in the pathogenesis of AMD, given its role in the inflammatory response. Clearly, the genes illustrated in our pathway map are involved in multiple processes. Although it is a hypothetical pathway, the map may help to explain the protective anti-inflammatory and anti-angiogenic role of vitamin D in the likelihood of developing AMD.

Figure 6

IPA pathway analysis of vitamin D metabolism. Graphical representation of the relationships between genes/proteins (represented as nodes). The biological relationship between two nodes is represented as a line (segmented to illustrate indirect relationships and continuous to illustrate direct relationships). All lines are supported by at least one reference from the literature or from canonical information stored in the Ingenuity Pathways Knowledge Base. Nodes are displayed using various shapes to illustrate the functional class of the gene product.

Analysis of interactions between gene variants and environmental factors

For gene-gene interactions, all CYP24A1 SNPs were tested for interaction with well-established risk genotypes in ARMS2/HTRA1, CFH and RORA for which we had previous genotyping data and, furthermore, for interaction with VDR rs2189480. Analysis of all AMD subtypes compared with unaffected individuals showed several CYP24A1 SNPs and VDR rs2189480 to interact modestly, but these associations did not remain significant after adjusting for multiple comparisons. The pooled case-only analysis showed no significant interactions between CYP24A1 SNPs and VDR rs2189480.

Interaction between both CYP24A1 and VDR variants and sunlight exposure was also tested using CLR in SAS. No significant interaction was found between UV irradiation and any of the variants tested (data not shown).

Discussion

Many pathways and processes involved in the progression of systemic disorders such as cancer, cardiovascular disease and autoimmune diseases are also thought to be involved in AMD. It is therefore likely that valuable insights into the pathophysiology of AMD may be gained by studying pathways common to AMD and other well-characterised diseases that may overlap in pathogenesis. In fact, cross-talk between disease fields has already proven beneficial for the treatment of AMD: anti-vascular endothelial growth factor (anti-VEGF) antibodies initially developed for the treatment of cancer are currently the most effective treatment for neovascular AMD [3].

The important role of vitamin D in diseases of complex aetiology, its anti-angiogenic properties and an association between vitamin D serum levels and early AMD led us initially to explore the association between vitamin D and neovascular AMD by examining vitamin D-associated epidemiological factors and genetic variants of vitamin D metabolism genes in an extremely discordant sibling pair cohort. Here, we showed a protective effect of UV exposure on the development of neo-vascular AMD. Further, we identified several CYP24A1 SNPs significantly associated with all subtypes of AMD in four diverse cohorts, both retrospectively and prospectively. Moreover, these significant associations were demonstrated in a meta-analysis of the four cohorts. To the best of our knowledge, this was the first study to demonstrate a protective effective of UV exposure on the development of AMD and a genetic link between vitamin D metabolism and AMD.

The role of sunlight in AMD pathophysiology remains controversial. Human exposure to natural sunlight may be both beneficial and detrimental. One benefit of cutaneous sunlight exposure is the production of vitamin D. Studies have also shown that vitamin D produced as a result of sunlight exposure reduces the risk of autoimmune diseases such as multiple sclerosis (MS) and inflammatory bowel disease (IBD) [63]. By contrast, natural sunlight (ie UV exposure) has long been known to have harmful effects on the human body, including permanent DNA damage. For example, studies have linked sunlight exposure to skin cancer, for both melanomas [64,65] and carcinomas [66]. In dermatological malignancies, the most detrimental factor is UVB, but UVA may also play a role.

To distinguish between the effects of general sunlight exposure and specifically UV exposure, sunlight exposure was self-assessed by patients and UV exposure was quantified separately by examining UV radiation, as measured in irradiance of UVB (290 nm to 320 nm), UVA (320 nm to 400 nm) and UV index. In addition, we examined vitamin supplementation, as well as sun-related variables including iris colour, self-reported sun exposure and skin cancer history. We also measured 25(OH)D in sera as a biomarker for vitamin D levels. To the best of the authors' knowledge, this was the first examination of sunlight exposure in the context of established AMD risk genotypes and smoking.

Quantification of sunlight exposure is difficult, as it is an ongoing life-long process [67-79]. One of the first studies of ocular exposure to UV based on interview data and field measurements in Chesapeake Bay watermen demonstrated no association between UV exposure and AMD [79]. Further analysis of this same dataset showed that there was in fact an association between AMD and ocular exposure, but to blue light specifically [76,78]. Alternatively, various studies have tested the association of sunlight exposure and AMD by measuring various sun-related variables. For example, Khan et al. investigated the role of variables such as iris colour, place of residence and hair colour but found no association between AMD and these sun related factors [80]. Cruickshanks et al. measured lifetime number of sunburns, use of hats and sunglasses and number of hours spent outdoors, and showed an association between leisure time spent outdoors and risk of early AMD [69]. Our study demonstrated a protective effect of moderate levels of UV radiation exposure and AMD, independent of smoking history and CFH and HTRA1 genotype. While experimental studies have shown that UV light has damaging effects on the retina--primarily via oxidative damage and formation of reactive oxygen species[81]--previous epidemiological studies have shown inconsistent results regarding an association between UV light and risk of AMD [67-79]. UV light is part of the sunlight spectrum, but only UVA (320-400 nm) and UVB (280-320 nm) reach the Earth's surface. UV light has traditionally been considered harmful because it induces DNA damage and causes oxidative damage to RPE cells [82-84]. Furthermore, UV light has been shown to have pro-angiogenic properties. UVA-irradiated RPE cells have been shown to have increased levels of prostaglandin-endoperoxide synthase 2 (also known as cyclooxygenase-2 [COX2]),[85] which is implicated in choroidal neovascularisation [86]. UVB-induced tumours express elevated levels of factors associated with angiogenesis, including VEGF, matrix metallopeptidase 2 (MMP2) and matrix metallopeptidase 3 (MMP3) [87].

Emerging evidence shows that solar radiation may have beneficial effects, however. Epidemiological studies have demonstrated that sunlight reduces the incidence and mortality of various cancers, including non-Hodgkin's lymphoma and melanoma [88,89]. UV light is known to suppress the immune response[90] and growing evidence points to the role of inflammation in AMD pathogenesis. Since aberrant inflammation is believed to be a contributing factor to AMD, the immunosuppressive effects of UV radiation may provide some protective effect.

We specifically examined the role of UV light and AMD in this study, although previous studies have mainly looked at ambient sunlight exposure. Previous studies have been inconsistent about the role of sunlight exposure and AMD [67-79,91,92]. Some studies suggested a positive association,[69,76,91,92] while others found no correlation [71,78,93]. To our knowledge, our study was the first to show a protective influence of UV light on neovascular AMD. Our finding may contrast with previous studies as a result of using a larger number of neovascular cases (n = 133) than that used in most previous studies. Furthermore, previous studies measured ambient solar radiation, which includes light of various wavelengths (UV, infrared, visible), while our study distinguished the contribution of UV irradiance.

Although our gene expression studies on vitamin D metabolism genes in human donor eyes found no significant association as a function of age or disease status, differences in expression levels between the macular and extramacular regions were observed. Specifically, the expression of both VDR and CYP27A1 was detected in the RPE-choroid, and CYP27A1 and CYP27B1 were detected in the retina, thus suggesting that the vitamin D metabolic genes VDR and CYP27B1 may function differently between the maculae and extramaculae, whether diseased or non-diseased. Also, AMD may not be a localised disease but rather may manifest systemically, and this may explain why CYP24A1 was not detected in human retinal tissue regardless of disease status. The gene expression changes that we observed in retinal cell lines suggested that CYP24A1 activity may alter over time and this would support its role in a late-onset disease such as AMD.

There were several limitations to this study. First, there was the possibility of recall bias. There was no definitive method to confirm an individual's recall of personal sun exposure and thus data could have been prone to misclassification. We attempted to minimise this bias by using a quantitative measurement of UV exposure based on the geographical location where the participants had resided for the majority of their lifetime; however, assessment of UV radiation based on residential history represents potential rather than actual exposure and could be subject to error. Secondly, although the subjects were questioned about the amount of time they spent outdoors, no information was gathered about their use of protection against sunlight, such as hats and sunglasses, which shield against UV light. Further, the subjects' sera were collected at various times throughout the year, and thus serum levels of vitamin D measured in this study may not accurately have represented the true range among the patient population, since vitamin D levels are influenced by sunlight exposure, diet and age.

Although a previous study by Parekh et al. showed an inverse relationship between vitamin D levels and early AMD but not late AMD,[34] this study was limited by a small number of patients with an advanced form of AMD (n = 10). A more recent study by Millen et al. was able to show a statistically significant association between serum 25(OH)D concentrations and early AMD in women younger than 75 years; specifically, that high serum 25(OH)D concentrations may be protective [94]. In our study, we further investigated the association of serum levels of 25(OH)D, the circulating form of vitamin D, with neovascular AMD by analysing the serum levels of 50 extremely discordant sibling pairs (n = 100). Our results were similar to those found by Millen and colleagues: a trend towards higher levels of vitamin D in the serum of unaffected patients compared with neo-vascular patients, although the difference was not statistically significant. Since the ability to produce vitamin D is diminished with advanced age (ageing results in decreased amounts of cutaneous 7-dehydrocholesterol, and by 70 years of age vitamin D3 synthesis is reduced by approximately 75 per cent[95]), it is not unexpected that serum vitamin D levels were low, given that the average age of our unaffected patients was 75.4 years. This could explain why the association was not significant. Additionally, since this was measured only once, the serum values reflected the vitamin D production/intake over a limited amount of time, thus suggesting enhanced random measurement error.

In addition to epidemiological data, the association of AMD with genetic variation in vitamin D pathway genes was further suggested by existing linkage data for genomic regions associated with AMD. The vitamin D pathway genes assayed in this study all lie within regions that have been associated with AMD -- specifically, 2q33[48] (CYP27A1), 12q13[36] (CYP27B1 and VDR), and 20q13[36] (CYP24A1).

The CYP24A1 SNPs significantly associated with neovascular AMD in the discovery cohort (either individually [rs6127118 and rs2769234] or as part of a haplotype [rs6068816, rs6127118, rs1570699, rs92760 and rs2762934]) were concentrated in the area between exon 6 and exon 12. This area was further examined and refined by direct sequencing and analysed in an extended sibling cohort. After controlling for smoking history, sex and age, we were able to show significant variation in CYP24A1 in all populations, both separately and, more importantly, in a meta-analysis. Variation within CYP24A1 was available in the dbGAP dataset http://www.ncbi.nlm.nih.gov/projects/gap/ but, when tested for association with AMD, no significance was observed. This lack of association could be explained by the small number of subjects in the dbGAP dataset, as compared with the current study (n = 500 and 2,528, respectively). Additionally, these data were based on a 100,000 SNP chip, and the CYP24A1 variation was imputed with low imputation quality, which is another limitation of the dbGAP dataset.

Moreover, pathway analysis showed that, in a hypothetical pathway, CYP24A1 was directly linked to another AMD-associated gene, HTRA1, for which the exact pathway is unknown. In a complex disease such as AMD, one would expect the combination of multiple effects, including the modest effect of CYP24A1, to contribute to disease causality. Although there was no statistical interaction seen between CYP24A1 and HTRA1, there may have been binding between the two that differed between those with and those without AMD. Further molecular studies should be performed to investigate this association.

An anti-angiogenic role for vitamin D has been well documented in the cancer literature [96-98]. Therefore, a role for vitamin D may be protective of AMD by its anti-angiogenic properties [27]. For example, VEGF expression was downregulated after tumour cells were treated with vitamin D [25,26]. Vitamin D may exert anti-inflammatory properties by enhancing T suppressor cell activity and down-regulating T helper cells, T cytotoxic cells and natural killer cells [28]. At physiological concentrations, vitamin D has also been shown to protect cell proteins and membranes from oxidative damage [30]. In addition to playing a role in cancer, the vitamin D pathway has been implicated in several autoimmune diseases. Several studies have found that low levels of vitamin D are associated with increased autoimmunity.

A role for vitamin D in AMD is not only plausible in terms of the known biological roles of vitamin D, but also because of the diminished ability to produce vitamin D with advanced age [99]. The problem of vitamin D deficiency is not only prevalent among the elderly: as of 2005, approximately 40 per cent of men and 50 per cent of women aged over 18 from the USA were estimated to have inadequate levels of 25(OH)D [100].

Conclusion

In summary, using multiple lines of evidence, this study demonstrated an association between the vitamin D pathway and AMD. Specifically, there is a protective association between UV radiation and neovascular AMD which is independent of smoking and CFH and HTRA1 polymorphisms. Additionally, we identified an association between genetic variation in the major vitamin D catabolic gene, CYP24A1, and risk for all forms of AMD and, specifically, the neovascular form of AMD. Further experiments are needed to elucidate the role of vitamin D and UV light in the pathogenesis of AMD. Ongoing large-scale trials of supplemental vitamin D may provide a unique opportunity to examine the effects of vitamin D on AMD in a randomised trial setting. Such studies may have important implications for the prevention and/or treatment of neovascular AMD by the regulation of modifiable lifestyle factors that influence levels of vitamin D.