- Research article
- Open Access
- Open Peer Review
Association of Interleukin-1 gene clusters polymorphisms with primary open-angle glaucoma: a meta-analysis
BMC Ophthalmology volume 17, Article number: 218 (2017)
Previous studies have associated the Interleukin-1 (IL-1) gene clusters polymorphisms with the risk of primary open-angle glaucoma (POAG). However, the results were not consistent. Here, we performed a meta-analysis to evaluate the role of IL-1 gene clusters polymorphisms in POAG susceptibility.
PubMed, EMBASE and Cochrane Library (up to July 15, 2017) were searched by two independent investigators. All case-control studies investigating the association between single-nucleotide polymorphisms (SNPs) of IL-1 gene clusters and POAG risk were included. Odds ratios (ORs) with 95% confidence intervals (CIs) were calculated for quantifying the strength of association that has been involved in at least two studies.
Five studies on IL-1β rs16944 (c. -511C > T) (1053 cases and 986 controls), 4 studies on IL-1α rs1800587 (c. -889C > T) (822 cases and 714 controls), and 4 studies on IL-1β rs1143634 (c. +3953C > T) (798 cases and 730 controls) were included. The results suggest that all three SNPs were not associated with POAG risk. Stratification analyses indicated that the rs1143634 has a suggestive associated with high tension glaucoma (HTG) under dominant (P = 0.03), heterozygote (P = 0.04) and allelic models (P = 0.02), however, the weak association was nullified after Bonferroni adjustments for multiple tests.
Based on current meta-analysis, we indicated that there is lack of association between the three SNPs of IL-1 and POAG. However, this conclusion should be interpreted with caution and further well designed studies with large sample-size are required to validate the conclusion as low statistical powers.
Glaucoma is the second largest cause of bilateral blindness worldwide that affects approximately 60 million people . Primary open-angle glaucoma (POAG), the major type of glaucoma, is characterized by progressive damage of the optic nerve and retinal ganglion cells (RGC), and subsequent irrevocable vision loss . Although the etiology is not entirely clear, genetic factors are thought to be a potential risk to POAG patients, including both those with normal and elevated intraocular pressure (IOP). Recently, the relationship between genetic polymorphisms of the immune system and were not perceived as being related to glaucoma, however, their relationship has been highly suspected and investigated recently .
Interleukin-1 (IL-1), an important mediator of inflammation, has been suggested to play a crucial role in neuro-degeneration such as Alzheimer’s disease . In view of the potential similarities in cellular mechanism leading to neuro-degenerative disorder between glaucoma and Alzheimer’s disease (AD), some investigators speculated that the polymorphism in IL-1 gene clusters may be a genetic predisposing factor for glaucoma. It has been reported that IL-1 was expressed endogenously in trabecular meshwork (MT) cells that controls a stress response specific to the aqueous outflow pathway and confers a protective response against glaucoma . Moreover, experiments in vivo showed that IL-1 promotes glutamate uptake by Müller cells and increases the number of surviving of RGCs .
Recently, several single nucleotide polymorphisms (SNPs) in the IL-1 gene cluster [i.e. IL-1β rs16944 (c. -511C > T), IL-1α rs1800587 (c. -889C > T) and IL-1β rs1143634 (c. +3953C > T)] have been investigated for association with POAG [7,8,9,10,11,12,13,14]. However, the results remain controversial due to some studies identified there was an association between IL-1 gene SPN and POA while others failed. Furthermore, to our knowledge, no relevant genome-wide association study (GWAS) or meta-analysis has been published on this subject. In the current study, we conducted a systematic review and meta-analysis to summarize the contradictory results from these relevant studies and to clarify the relationship between IL-1 cluster polymorphisms and POAG.
This meta-analysis was followed the the Preferred Reporting Items for Systematic Review and Meta-Analyses (PRISMA) statement .
To obtain eligible literature, a comprehensive search without a language restriction was conducted in the PubMed, EMBASE and Cochrane Library covering publications up to July 15, 2017, by using of the following keywords: (“glaucoma” or “open angle” or “POAG” or “intraocular hypertension”) AND (“interleukin-1” or “IL-1” or “interleukin 1”) AND (“polymorphism” or “SNP” or “single nucleotide polymorphism” or “variation” or “mutation”) (Additional file 1). The titles, abstracts, and full texts were screened by two independent investigators (L.J.H., F.Y.F.) to determine inclusion and disagreements were resolved by discussion between all authors until a consensus was determined. Review articles and bibliographies of other relevant studies were searched manually to identify additional eligible studies.
Inclusion and exclusion criteria
Literature selection had to meet the the following criteria: (1) case-control or cohort studies on the relationship of the SNPs of IL-1 with POAG; (2) all patients in the eligible studies meeting the diagnostic criteria for POAG; (3) had sufficient information for estimating odds ratios (ORs) with corresponding 95% confidence interval (CI). The exclusion criteria were: (1) case only studies; (2) abstract, case report and review papers; (3) repeated or overlapped publication. Only the most recent or complete study was used in this meta-analysis in case of the same patient population was included in several publications.
The methodological quality of eligible studies was evaluated independently by two authors (J.H.L. and Y.F.F.), according to a modified version of the Newcastle–Ottawa scale (NOS) for genetic association studies  (Additional file 2). NOS quality scores ranged between 0 and 9 stars. Studies with a score of five stars or greater were considered high quality.
According to the PRISMA guidance, the necessary information from eligible studies was extracted independently by two reviewers (L.J.H. and S.M.S.) using a standardized form. The following data was retrieved: surname of first author, publication date, region, ethnicity, number of cases and controls, age, subtypes of POAG, genotyping method and genotypes frequency, evidence of Hardy-Weinberg equilibrium (HWE) in control group, etc. When the allele or genotype counts were not given specially in some articles, they were calculated from the frequencies and then rounded to the nearest integer. Disagreements were resolved by discussion between all authors and subsequent consensus.
Statistical analysis was undertaken using RevMan software (version 5.0; Cochrane Collaboration, Oxford, United Kingdom) and STATA 11.0 software (Stata Corporation, Texas, USA). The strength of association between IL-1 SNPs and POAG was estimated by ORs with 95% confidence intervals (CIs). The pooled ORs were performed for the allele model (T versus C), homozygote model (CC versus TT), heterozygote model (TC versus CC), dominant model (TT + CT versus CC) and recessive model (TT versus CT + CC), respectively.
I2value was applied to evaluate the between-study heterogeneity, with <25%, 25%–50%, and >50% to represent low, moderate, and high degree of heterogeneity, respectively . In addition, a Q-statistic test was conducted and the heterogeneity was considered significant at P < 0.10. The summary OR and 95% CI for each polymorphism were pooled by using the fixed-effect model (Mantel-Haenszel) when no significant heterogeneity was observed among studies. Otherwise, the random effect model (DerSimonian and Laird) was adopted.
In association analysis, a pooled P value of less than 0.05 was considered as suggestive evidence for a genetic association. Aim to control false positive error rate, the Bonferroni method was used to adjust for multiple comparisons . We performed 20 times comparisons for IL-1β rs16944, 15 times comparisons for IL-1α rs1800587 and IL-1β rs1143634 in POAG respectively, therefore, the P values that were less than 0.05/20 in IL-1β rs16944 and 0.05/15 in IL-1α rs1800587 and IL-1β rs1143634 showed significance.
Power analysis regarding the association of each IL-1 SNP with POAG was performed using the statistical software Power and Sample Size Calculation (PS) version 3.1.2 (http://biostat.mc.vanderbilt.edu/wiki/Main/PowerSampleSize) and a value greater than 0.8 meant high statistical power. The HWE was assessed by Fisher’s exact test. Subgroup analyses were performed based on ethnicity, as well as type of POAG including high tension glaucoma (HTG) (IOP > 21 mmHg) and non-HTG (IOP < 21 mmHg, with or without a strict diurnal testing). Sensitivity analysis was performed to assess the stability of individual studies, and the potential of publication bias was assessed by Begg’s funnel plot  and further evaluated by Egger’s linear regression test .
The flow chart for the study selection is summarized in Fig. 1. The initial search yielded 39 articles and 9 duplicate publications were excluded (Additional file 1 ). Based on the title and content of the abstract, 22 studies were further excluded. Nine articles with full-texts were further reviewed and 1 publication was excluded because it did not conform to our inclusion criteria. Thus, 8 case-control studies [7,8,9,10,11,12,13,14] were finally included in this meta-analysis and the main characteristics of each study are listed in Additional file 3: Table S1. The publication year of the included studies ranged from 2003 to 2016. All of the 8 studies had scores ≥5* and the mean score was 5.6*. The quality of these studies is summarized in Additional file 4: Table S2. In total, 64 SNPs have been investigated at least once in these 8 studies, among which 3 SNPs (IL-1β rs16944, IL-1β rs1143634, IL-1α rs1800587) were tested in at least 2 studies and included in the present meta-analysis. In all studies, the distributions of genotypes of all SNP in the control groups were consistent with HWE (Table 1).
Association between SNP rs16944 of IL-1β gene and POAG
SNP rs16944 was assessed in 5 case-control studies with a total of 1053 POAG cases and 986 controls (Table 1). No association between this SNP and POAG risk was found in the overall populations in dominant (OR = 0.97, 95% CI 0.71 to 1.31, P = 0.83), recessive (OR = 1.07, 95% CI 0.71 to 1.61, P = 0.74), homozygote (OR = 1.03, 95% CI 0.60 to 1.76, P = 0.91), heterozygote (OR = 0.97, 95% CI 0.79 to 1.20, P = 0.91) and allelic comparison (OR = 1.00, 95% CI 0.78 to 1.28, P = 1.00) models (Table 2; Fig. 2). In subgroup analysis stratified by type of POAG, the meta-analyses indicated no significant association of this SNP with POAG in any of the genetic models in patients with HTG and non-HTG, respectively. In subgroup analysis stratified by ethnicity, the meta-analyses indicated no significant association of this SNP with POAG in any of the genetic models in Asians.
Association between SNP rs1800587 of IL-1α gene and POAG
Four case-control studies had investigated the relationship between SNP rs1800587 and susceptibility to POAG with a total of 822 cases and 714 control subjects (Table 1). No association between this SNP and POAG risk was found in the overall populations in dominant (OR = 0.99, 95% CI 0.64 to 1.54, P = 0.98), recessive (OR = 1.01, 95% CI 0.64 to 1.61, P = 0.96), homozygote (OR = 1.19, 95% CI 0.44 to 3.17, P = 0.73), heterozygote (OR = 1.02, 95% CI 0.81 to 1.29, P = 0.84) and allelic comparison (OR = 1.08, 95% CI 0.75 to 1.55, P = 0.67) models (Table 2; Fig. 3). In subgroup analysis stratified by type of POAG, the meta-analyses indicated no significant association of this SNP with POAG in any of the genetic models in patients with HTG and non-NTG, respectively (Table 2 ).
Association between SNP rs1143634 of IL-1β gene and POAG
Three case-control studies had investigated the relationship between SNP rs1143634 and susceptibility to POAG with a total of 798 cases and 730 control subjects (Table 1). No association between this SNP and POAG risk was found in the overall populations in dominant (OR = 1.34, 95% CI 0.98 to 1.82, P = 0.06), recessive (OR = 1.56, 95% CI 0.50 to 4.84, P = 0.44), homozygote (OR = 1.64, 95% CI 0.53 to 5.11, P = 0.39), heterozygote (OR = 1.32, 95% CI 0.96 to 1.81, P = 0.09) and allelic comparison (OR = 1.32, 95% CI 0.99 to 1.75, P = 0.06) models (Table 2; Fig. 4). In subgroup analysis stratified by type of POAG, the meta-analyses indicated no significant association of this SNP with POAG in any of the genetic models in non-HTG. Although two study involving 216 HTG cases and 256 controls showed significant association with POGA risk under dominant (OR = 1.69, 95%CI: 1.07 to 2.12, P = 0.03), heterozygote (OR = 1.65, 95% CI: 1.03 to 2.65, P = 0.04) and allelic comparison (OR = 1.61, 95% CI: 1.06 to 2.42, P = 0.02), there was no statistical significance when Bonferroni correction was considered (Table 2 ).
Power calculation on the pooled frequencies showed that the statistical powers were all lower than 80% for allele models of the three SNPs (Table 2).
To examine the stability of pooled results, sensitivity analysis was conducted by sequentially excluding individual studies and calculating the pooled ORs for the remaining studies. As shown in Fig. 5, we found no individual study affected the pooled OR in the dominant models of all three SNPs.
Funnel plots (Fig. 6) and Egger’s linear regression test found no evidence of publication bias for the dominant models of SNP rs16944 (PEgger’s = 0.57), SNP rs1800587 (PEgger’s = 0.98) and SNP rs1143634 (PEgger’s = 0.56).
This is the first meta-analysis to summarize the evidence of associations between IL-1 gene cluster polymorphisms and susceptibility of POAG. Our data revealed that the SNPs of IL-1β rs16944, IL-1α rs1800587 and IL-1β rs1143634 were not associated with the POAG risk. Stratification analyses showed that IL-1β rs1143634 has a suggestive associated with the risk of HTG, however, this association was feeble after Bonferroni adjustments.
Recently, polymorphisms in the IL-1 gene clusters have been shown in chronic neurodegenerative disease, such as multiple sclerosis , Parkinson’s disease  and Alzheimer’s disease . The C allele of rs16944 and rs1800587, and the T allele of rs1143634 has been speculated to influence the risk of POAG by altering the expression level of the respective proteins, because the sites of rs1800587 and rs16944 are located within the transcriptional promoter regions of the IL-1α and IL-1β genes, respectively, and the site of rs1143634 is within the coding region of the IL-1β gene [24, 25]. IL-1 is a key mediator of immune and inflammatory responses, which act directly to mediate a number of cellular responses. Laboratory studies have showed that IL-1 could promote ganglion cell loss and optic nerve damage by increasing matrix mellanoproteinase-9 (MMP-9) synthesis in experimental models [26, 27]. IL-1 has also been reported to induce nitric oxide synthesis and reactive oxygen species (ROS), which involved in RGC damage leading to neurodegeneration . However, some studies have put forward argument and their results pointed to certain neuro-protective role of IL-1 [5, 6]. Therefore, we undertake the present meta-analysis on all the available data to establish a more robust estimate of the association between IL-1 gene clusters polymorphism and POAG.
Although the meta-analysis for overall studies showed that all the three IL-1 clusters polymorphisms were not associated with the POAG risk, the results of subgroup analyses should arouse our attention. On the one hand, genetic heterogeneity for POAG may exist in different populations. Stratified analyses in the current meta-analysis based on ethnicity showed that there is no statistical evidence of significant association between the rs16944 and POAG in Asian population. However, one study in Caucasians showed that the SNP rs16944 is important risk factors associated with POAG . Moreover, of the included studies for the Asian population, one studies  found a possible association of the SNP rs1143634 with POAG in Taiwan population but this association could not be further replicated in Singaporean  Chinese and Indian subjects . Thus, more studies with POAG cohort from different ethnic background are needed to further define this association.
On the other hand, genetic heterogeneity for POAG may exist in different subtypes of this disease. For example, Wang et al. found an increased risk for individuals carrying T allele of IL-1α rs1800587 in HTG patients with an IOP > 21 mmHg  but not for normal-tension glaucoma (NTG) cases with IOP < 21 mmHg . Similarly, the SNP rs1143634 seem to be associated with HTG rather than NTG. Glaucomatous damage to the retina and optic nerve is often accompanied by pathological IOP elevation. These findings of our meta-analysis suggested that the pathogenesis and effect of IL-1 may be different between HTG and NTG. Laboratory and clinical studies are needed to further clarify the molecular pathogenesis of the two subtypes of POAG as well as the role of IL-1 on IOP elevation.
Publication bias should be considered as it is an important factor that affects the reliability of the results of meta-analyses. In the current study, both Begg’s and Egger’s test were used to assess the potential publication bias and failed to detected a significant bias in all genetic models (Fig. 5), demonstrating the robustness and credibility of the present meta-analysis. However, some limitations of this meta-analysis should be taken into careful consideration. Firstly, some unescapable bias may exist in the results as only published studies with full-text were included in our analysis. Secondly, haplotype analysis for the possibility of linkage disequilibrium between SNPs was not performed in the present meta-analysis. How et al.  found that The TT haplotype of rs1143634 and rs16944 together and the TTT haplotype of rs16944, rs1800587 and rs1143634 was significantly more common in normal control subjects than in those with POAG. Further investigations of the haplotypic effect of a gene and the study of multiple polymorphisms in different genes are needed. Thirdly, the results of power calculations indicated that the combined sample sizes (overall and subgroups) in the current meta-analyses were still inadequate and underpowered to detect the association of IL-1 gene SNPs with POAG, due to limited availability of published data. Finally, the strength of association between IL-1 gene SNPs and POAG was carried out by unadjusted estimate, we didn’t adjust pool results by the factors like the age, gender, disease severity and genotyping procedure.
In conclusion, this meta-analysis did not demonstrate an association between the IL-1 SNPs (rs16944, rs1800587 and rs1143634) and the risk of POAG. Nevertheless, this conclusion should be interpreted with caution and well-designed studies with larger cohorts are urgently warranted to verify the present findings as low statistical powers.
High tension glaucoma
Normal tension glaucoma
Primary open-angle glaucoma
Retinal ganglion cells
Reactive oxygen species
Science Citation Index
Single nucleotide polymorphism
Quigley HA, Broman AT. The number of people with glaucoma worldwide in 2010 and 2020. Br J Ophthalmol. 2006;90:262–7.
Kwon YH, Fingert JH, Kuehn MH, Alward WL. Primary open-angle glaucoma. N Engl J Med. 2009;360:1113–24.
Wax MB. Is there a role for the immune system in glaucomatous optic neuropathy? Curr Opin Ophthalmol. 2000;11:145–50.
Mun MJ, Kim JH, Choi JY, Jang WC. Genetic polymorphisms of interleukin genes and the risk of Alzheimer's disease: an update meta-analysis. Meta Gene. 2016;8:1–10.
Wang N, Chintala SK, Fini ME, Schuman JS. Activation of a tissue-specific stress response in the aqueous outflow pathway of the eye defines the glaucoma disease phenotype. Nat Med. 2001;7:304–9.
Namekata K, Harada C, Guo X, Kikushima K, Kimura A, Fuse N, et al. Interleukin-1 attenuates normal tension glaucoma-like retinal degeneration in EAAC1-deficient mice. Neurosci Lett. 2009;465:160–4.
Lin HJ, Tsai SC, Tsai FJ, Chen WC, Tsai JJ, Hsu CD. Association of interleukin 1beta and receptor antagonist gene polymorphisms with primary open-angle glaucoma. Ophthalmologica. 2003;217:358–64.
Wang CY, Shen YC, Lo FY, Su CH, Lee SH, Lin KH, et al. Polymorphism in the IL-1alpha (−889) locus associated with elevated risk of primary open angle glaucoma. Mol Vis. 2006;12:1380–5.
How AC, Aung T, Chew X, Yong VH, Lim MC, Lee KY, et al. Lack of association between interleukin-1 gene cluster polymorphisms and glaucoma in Chinese subjects. Invest Ophthalmol Vis Sci. 2007;48:2123–6.
Wang CY, Shen YC, Lo FY, Su CH, Lee SH, Tsai HY, et al. Normal tension glaucoma is not associated with the interleukin -1alpha (−889) genetic polymorphism. J Glaucoma. 2007;16:230–3.
Wang CY, Shen YC, Su CH, Lo FY, Lee SH, Tsai HY, et al. Investigation of the association between interleukin-1beta polymorphism and normal tension glaucoma. Mol Vis. 2007;13:719–23.
Mookherjee S, Banerjee D, Chakraborty S, Banerjee A, Mukhopadhyay I, Sen A, et al. Association of IL1A and IL1B loci with primary open angle glaucoma. BMC Med Genet. 2010;11:99.
Markiewicz L, Majsterek I, Przybylowska K, Dziki L, Waszczyk M, Gacek M, et al. Gene polymorphisms of the MMP1, MMP9, MMP12, IL-1beta and TIMP1 and the risk of primary open-angle glaucoma. Acta Ophthalmol. 2013;91:e516–23.
Mookherjee S, Banerjee D, Chakraborty S, Mukhopadhyay I, Sen A, Ray K. Evaluation of the IL1 gene cluster single nucleotide polymorphisms in primary open-angle glaucoma pathogenesis. Genet Test Mol Biomarkers. 2016;20:633–6.
Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gotzsche PC, Ioannidis JP, et al. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate health care interventions: explanation and elaboration. J Clin Epidemiol. 2009;62:e1–e34.
Lv YB, Wang Y, Ma WG, Yan DY, Zheng WL, Chu C, et al. Association of Renalase SNPs rs2296545 and rs2576178 with the risk of hypertension: a meta-analysis. PLoS One. 2016;11:e158880.
Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003;327:557–60.
Boehringer S, Epplen JT, Krawczak M. Genetic association studies of bronchial asthma--a need for Bonferroni correction? Hum Genet. 2000;107:197.
Stuck AE, Rubenstein LZ, Wieland D. Bias in meta-analysis detected by a simple, graphical test. Asymmetry detected in funnel plot was probably due to true heterogeneity. BMJ. 1998;316(469):470–1.
Egger M, Davey SG, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315:629–34.
Huang J, Xie ZK, Lu RB, Xie ZF. Association of interleukin-1 gene polymorphisms with multiple sclerosis: a meta-analysis. Inflamm Res. 2013;62:97–106.
Arman A, Isik N, Coker A, Candan F, Becit KS, List EO. Association between sporadic Parkinson disease and interleukin-1 beta −511 gene polymorphisms in the Turkish population. Eur Cytokine Netw. 2010;21:116–21.
Nicoll JA, Mrak RE, Graham DI, Stewart J, Wilcock G, MacGowan S, et al. Association of interleukin-1 gene polymorphisms with Alzheimer's disease. Ann Neurol. 2000;47:365–8.
Hall SK, Perregaux DG, Gabel CA, Woodworth T, Durham LK, Huizinga TW, et al. Correlation of polymorphic variation in the promoter region of the interleukin-1 beta gene with secretion of interleukin-1 beta protein. Arthritis Rheum. 2004;50:1976–83.
Dominici R, Cattaneo M, Malferrari G, Archi D, Mariani C, Grimaldi LM, et al. Cloning and functional analysis of the allelic polymorphism in the transcription regulatory region of interleukin-1 alpha. Immunogenetics. 2002;54:82–6.
Zhang X, Chintala SK. Influence of interleukin-1 beta induction and mitogen-activated protein kinase phosphorylation on optic nerve ligation-induced matrix metalloproteinase-9 activation in the retina. Exp Eye Res. 2004;78:849–60.
Chintala SK, Zhang X, Austin JS, Fini ME. Deficiency in matrix metalloproteinase gelatinase B (MMP-9) protects against retinal ganglion cell death after optic nerve ligation. J Biol Chem. 2002;277:47461–8.
Neufeld AH, Sawada A, Becker B. Inhibition of nitric-oxide synthase 2 by aminoguanidine provides neuroprotection of retinal ganglion cells in a rat model of chronic glaucoma. Proc Natl Acad Sci U S A. 1999;96:9944–8.
Availability of data and materials
All data has been shared in the Figures and Tables.
Ethics approval and consent to participate
Consent for publication
All authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The full details of databases searching terms. (DOCX 36 kb)s (DOCX 36 kb)
Modified Newcastle-Ottawa Scale for studies of genetic association. (DOCX 16 kb) (DOCX 16 kb)
Characteristics of included studies. (DOC 67 kb) (DOC 67 kb)
Quality of included studies. (DOC 68 kb) (DOC 68 kb)