Molecular Vision 2013; 19:822-828
Received 13 March 2012 | Accepted 03 April 2013 | Published 05 April 2013
Inderjeet Kaur,1 Stuart Cantsilieris,2 Saritha Katta,1 Andrea J. Richardson,2 Maria Schache,2 Rajeev R. Pappuru,3 Raja Narayanan,3 Annie Mathai,3 Ajit B. Majji,3 Nicole Tindill,2 Robyn H. Guymer,2 Subhabrata Chakrabarti,1 Paul N. Baird2
First three and last two authors contributed equally to this work.
1Hyderabad Eye Research Foundation, L.V. Prasad Eye Institute, Hyderabad, India; 2Centre for Eye Research Australia, University of Melbourne, Royal Victorian Eye and Ear Hospital, Melbourne, Victoria, Australia; 3Hyderabad Eye Institute, L.V. Prasad Eye Institute, Hyderabad, India
Correspondence to: Subhabrata Chakrabarti, Kallam Anji Reddy Molecular Genetics Laboratory, Champalimaud Translational Centre, Brien Holden Eye Research Centre, L.V. Prasad Eye Institute, Hyderabad 500034, India; Phone: +91-40-30612524; FAX: +91-40-23548271; email: email@example.com
Purpose: The ARMS2/HTRA1 genes at the 10q26 locus have been associated with risk of age-related macular degeneration (AMD), with the most significantly associated variants being A69S (rs10490924), del443ins54 (EU427539) and rs11200638. We wished to explore the association of the del443ins54 in two ethnically different populations from India and Australia.
Methods: The del443ins54 was screened in a large cohort of ~1500 subjects from these two populations by a combination of PCR-based agarose gel electrophoresis and validated by resequencing. Statistical analysis comprised the calculations of allele, genotype and haplotype frequencies along with their p values and corresponding odds ratios (OR), and 95% confidence intervals (95% CI) and measures of linkage disequilibrium (LD).
Results: The del443ins54 was significantly associated with AMD in both the Indian (p=1.74×10−13; OR=2.80, 95%CI, 2.12–3.70) and Australian cohorts (p=2.78×10−30; OR=3.15, 95%CI, 2.58–3.86). These associations were similar to those previously identified for the A69S and the rs11200638 variant in these populations that also exhibited high degrees of LD (D’ of 0.87-0.99). A major risk haplotype of “T-indel-A” (p=5.7×10−16; OR=3.16, 95%CI, 2.34–4.19 and p=6.33×10−30; OR=3.15, 95%CI, 2.57–3.85) and a protective haplotype of “G-wild type-G” (p=2.35×10−11; OR=0.39, 95%CI, 0.29–0.52 and p=1.02×10−30; OR=0.31, 95%CI, 0.25–0.38) were identified in the Indian and Australian cohorts, respectively.
Conclusions: These data provide an independent replication of the association of del443ins54 variant in two different ethnicities, despite differences in allele and haplotype frequencies between them. High levels of LD in both populations limit further genetic dissection of this region in AMD.
Age-related macular degeneration (AMD) is a complex multifactorial disease and a leading cause of irreversible blindness in the world . A major AMD susceptibility locus on 10q26 has been found to harbor risk associated variants in ARMS2 (rs10490924) and HTRA1 (rs11200638) in multiple populations worldwide [1-4]. An insertion-deletion (indel) polymorphism (EU427539) that affects the stability of ARMS2 mRNA by the removal of a polyadenylation signal (443 bases) and insertion of a 54bp AU rich element in the 3′-UTR (del443ins54), has also been identified in the 3′ end of the ARMS2 gene as increasing risk of AMD by several fold in both Caucasian  as well as Asian populations. Previous associations for this indel vary (from p=3.5×10−5 to p=8.4×10−34) in individuals of European origin [6-9], whereas, this has only been reported in two non-European cohorts consisting of Han Chinese  and Japanese populations .
We have previously reported a significant association of the A69S (rs10490924) and rs11200638 variants with AMD in both South Indian and Australian cohorts [11,12]. We now wished to determine the risk conferred by del443ins54 and its combined effect with both the A69S and rs11200638 variants with AMD susceptibility in two ethnically different cohorts from Southern India and Australia.
The study protocols adhered to the tenets of the Declaration of Helsinki and were approved by the Institutional Review Boards of the Royal Victorian Eye and Ear Hospital, Melbourne, Australia and L.V. Prasad Eye Institute, Hyderabad, India. The del443ins54 was screened in end stage AMD cases (mainly choroidal neovascular) and normal controls from cohorts in India (n=433) and Australia (n=1054). The detailed methods of clinical diagnosis along with the inclusion and exclusion criteria have been previously reported [11,12]. Amplification was performed using forward (5′-TCT GTG CAG CTG GTG AAA TC-3′) and reverse (5′-TCC AGG GTG GTG TAA TCC AT-3′) primers at an annealing temperature of 61 °C. Amplicons were visualized on a 2% agarose gel and genotypes directly scored from the gels. Subsets of samples were further validated by bi-directional sequencing on an automated DNA sequencer (ABI 3100), using the BigDye chemistry as per manufacturer’s guidelines (both from Applied Biosystems, Foster City, CA). Genotyping results were independently validated by a second investigator who was masked to the phenotype data.
Allele and genotype frequencies were determined by the gene counting method and estimates of Hardy–Weinberg equilibrium were assessed. Odds ratios (OR) and 95% confidence intervals (95%CI) were calculated to assess the risk conferred for each variant using the PLINK software . Haplotypes were generated using various combinations of the A69S, the indel and rs11200638 variants and the estimated haplotype frequencies and linkage disequilibrium (LD) were assessed with the Haploview software (version 4.2) that uses the EM algorithm .
All statistical analyses were based on samples where genotyping was successful across all three A69S, del443ins54 and rs11200638 genetic variants. Each variant was in Hardy–Weinberg equilibrium in both the Indian and Australian cohorts (p>0.05). Allele frequencies for each of the risk variants (T allele of A69S, presence of the indel in del443ins54 and the A allele of rs11200638) exhibited a relatively higher frequency (0.60–0.63) in the Indian cohort compared to the European cohorts (0.36–0.53) but lesser than that reported in Han Chinese (0.73–0.77) or the in the Japanese (0.86–0.88) populations (Table 1). The allele frequencies in the Australian cohort were similar (0.44) to those previously reported for other European cohorts (Table 1).
The frequency of the del443ins54 variant was significantly higher among AMD cases than controls in both the Indian (p=1.74×10−13; OR=2.80, 95%CI, 2.12–3.70) and Australian (p=2.78×10−30; OR=3.15, 95%CI, 2.58–3.86), cohorts (Table 1). This increased risk was also observed with respect to the del443ins54 genotypes. These findings were similar for the A69S (ARMS2) and the rs11200638 (HTRA1) variants in both the Indian and Australian cohorts, respectively (Table 2).
Homozygosity of the indel and the other variants were strongly associated with an increased risk of AMD in both the Indian and Australian cohorts. Combined homozygosities at the A69S and the rs11200638 along with the indel variant did not alter the risk of AMD significantly either in the Indian (OR=7.69, 95%CI, 4.07–14.51) or Australian cohorts (OR=10.61, 95%CI, 7.05–15.96).
The measure of linkage disequilibrium (LD) between the A69S, del443ins54 and rs11200638 variants were remarkably high across this 10q26 region with relatively higher values in the Australian (D’=0.99; r2=0.98) compared to the Indian (D’=0.87; r2=0.71) cohorts (Figure 1). Two major haplotypes (frequency >5%) were identified across the three variants with “T-Indel-A” being the risk haplotype and “G-WT (wild-type)-G” being protective in the Indian and Australian cohorts, respectively (Table 3). Different pairwise haplotype combinations with either the ‘T’ or ‘A’ allele at A69S in presence of the del443ins54 or its wild-type form along with the ‘G’ or ‘A’ allele of rs11200638, did not substantially alter the p values or ORs observed for either the risk or protective haplotypes as opposed to when all three variants were assessed together, reinforcing the observation of the high LD between variants in these two genes (Table 3).
These data provide an independent replication of the association of the ARMS2 del443ins54 variant in two cohorts, and to the best of our knowledge for the first time among South Indians with AMD. The strong association of the del443ins54 along with the A69S and rs11200638 variants in the Indian and Australian cohorts were consistent with that observed in other populations [6-10]. Haplotype analysis with these three variants indicated that inclusion of del443ins54 in the haplotype neither increased nor decreased the risk of AMD in either cohort (Table 3).
The ARMS2 and HTRA1 genes are in high LD in European populations and thus dissecting out the role of one gene over the other has proved difficult. The advantage of undertaking a comparative analysis of genetic variants in populations of differing ethnicities expands the genetic diversity available and may provide the opportunity of identifying a more defined but associated region for further study. The current study highlighted similar degrees of associations across these three variants despite a relatively lower LD between the A69S, del443ins54 and rs11200638 variants in the Indian compared to the Australian cohort. However, it did confirm the presence of stratification differences between ethnicities with the allele frequency of the Indel of del443ins54 in South Indians being higher at 0.63 in cases compared to that in European populations (0.36–0.53) but lower than other Asian populations (0.73–0.88). The allele frequency in cases is similar to that previously shown by us in the assessment of the A69S (0.63) and the rs11200638 (0.60) variants of the ARMS2 and HTRA1 genes respectively, in the South Indian cohort. Evidence of population stratification has also been observed in AMD studies of the protective CFHR3–1 deletion, with the highest frequencies of the deleted allele being present in African populations (16%–20%) compared to Asians (<2%) [15,16].
The potential role of the ARMS2 and HTRA1 gene in AMD is still unclear but functional dissection of the effect of the rs11200638 promoter variant in the HTRA1 gene has revealed that this variant resides within a putative transcription binding site for the factors AP2α and SRF (serum response factor) [3,4]. Initial investigation of the influence of the homozygous risk genotype on HTRA1 expression levels revealed consistently higher levels of expression with the AA genotype compared to the GG genotype [3,4]. In contrast, other studies of the rs11200638 variant have revealed no functional effect on HTRA1 expression [2,6,17].
Analysis of the chromosome 10q26 risk haplotype inclusive of the ARMS2 del443ins54 indel found decreased ARMS2 expression and almost 3.0 fold increase in HTRA1 expression . Interestingly, a subsequent study has shown that while the ARMS2 risk del443ins54 results in decrease in mRNA transcription levels of the ARMS2 gene, a non risk associated variant (rs2736911) also leads to significantly reduced ARMS2 transcript levels suggesting that ARMS2 protein deficiency alone is unlikely to be pathogenic in AMD . A functional role for ARMS2 in mitochondrial homeostasis has also been suggested and the biology concerning mitochondrial dysfunction and the effects on age supports this notion [2,5]. However, subsequent immunofluorescence and immunoblot experiments localized ARMS2 in retinal epithelial ARPE-19 cells and COS7 transfected cells to the cytosol rather than the mitochondria suggesting that ARMS2 may not confer risk to AMD through the mitochondrial pathway . Studies concerning the effects of AMD risk variants on HTRA1 expression are equivocal and further investigations on the functional role of these variants are required.
In conclusion, we provide convincing evidence for the association of the del443ins54 variant with AMD, despite differences in allele, genotype and haplotype frequencies and LD across the 10q26 region reflecting population stratification differences in two different ethnicities. AMD, a complex multi-factorial disease is associated with multiple genomic regions with varying magnitudes of effect and the relevance of genetic associations differ between populations. Further, elucidation of the genetic basis of this disease through the analysis of individuals from different ethnic groups has the potential to provide useful insights into the genetic diversity of risk and protective variants within a gene as well as their contributions to disease. Also, meaningful genetic dissection of the ARMS2 and HTRA1 gene in this region will require much larger patient cohorts than have currently been assessed, or through the identification of other ethnic populations which show relatively lower levels of LD over this 10q26 region.
The authors thank all the Australian and Indian patients and the normal volunteers for their participation and Drs. Nazimul Hussain, Taraprasad Das, Anjli Hussain and Avinash Pathangay for their help in collecting some of the earlier patients’ data. RHG and PNB acknowledge an NHMRC Practitioner Award and NHMRC Fellowship, respectively. SK and Stuart Cantsilieris acknowledge a Senior Research Fellowship of the Council of Scientific and Industrial Research (CSIR), Government of India, and an Australian Postgraduate Award, respectively. CERA receives Operational Infrastructure Support from the Victorian Government, Australia. Grant support: This work was supported by the Australia-India Strategic Research Fund (AISRF) jointly funded through the Department of Biotechnology (DBT), Government of India (SC and IK) and the Department of Innovation, Industry, Science and Research (DIISR), Government of Australia (PNB and RHG), the National Health and Medical Research Council (NHMRC) Centre for Clinical Research Excellence #529923 - Translational Clinical Research in Major Eye Diseases.