Molecular Vision 2013; 19:526-535 <http://www.molvis.org/molvis/v19/526>
Received 11 November 2012 | Accepted 25 February 2013 | Published 27 February 2013

Roles of tumor necrosis factor alpha gene polymorphisms, tumor necrosis factor alpha level in aqueous humor, and the risks of open angle glaucoma: A meta-analysis

Xiangyang Xin,1 Lili Gao,2 Tong Wu,3 Fengyuan Sun3

1First Center Clinical College, Tianjin Medical University, Tianjin, P. R. China; 2Department of Ophthalmology, Inner Mongolia Baogang Hospital, Inner Mongolia, P. R. China; 3Department of Ophthalmology, Tianjin First Center Hospital, Tianjin Medical University, Tianjin, P. R. China

Correspondence to: Fengyuan Sun, Department of Ophthalmology, Tianjin First Center Hospital, Tianjin Medical University, Tianjin, 300192, P. R. China; email: address: fysuntianjin@126.com

Abstract

Purpose: This meta-analysis was performed to clarify the association between tumor necrosis factor alpha (TNF-α) gene polymorphisms and open angle glaucoma (OAG) risks, and the association between the TNF-α level in aqueous humor (AH) and the risks of glaucoma.

Methods: A computerized literature search was performed for the relevant available studies from three databases including PubMed, ISI Web of Science, and Embase. The fixed or random effect model was selected based on the heterogeneity test using the Q test and the I2 statistic. The associations between TNF-α gene polymorphisms and OAG risks were estimated by calculating pooled odds ratios (ORs) and the 95% confidence interval (CI), while a pooled standardized mean difference with 95% CI was used for the comparison of TNF-α levels in AH between patients with OAG and controls. Publication bias was estimated using Begg’s funnel plots and Egger’s regression test.

Results: A total of 14 (1,182 cases and 3,003 controls), five (808 cases and 1,039 controls), three (645 cases and 666 controls), and three studies (404 cases and 625 controls) were finally included in the analyses for the associations between TNF-α −308G/A, −857C/T, −863C/A, and −238G/A polymorphisms and the risks of OAG, respectively. The combined results showed that the TNF-α −308G/A gene polymorphism was significantly associated with risks of high-tension glaucoma (A versus G: OR=1.660, 95% CI=1.033–2.667; AA/AG versus GG: OR=1.713, 95% CI=1.10–2.651), but not with normal tension glaucoma or exfoliation glaucoma. Ethnicity-stratified analysis revealed that a significant association also existed in Asians (A versus G: OR=1.947, 95% CI=1.097–3.456; AA/AG versus GG: OR=1.949, 95% CI=1.140–3.332). None of the other polymorphisms was significantly associated with OAG risks. Furthermore, the pooled results of six studies showed that the TNF-α levels in the AH of patients with OAG was higher than that of the control subjects (standardized mean difference=0.517, 95% CI=0.207–0.826, p=0.001). Probability of publication bias was low across all comparisons illustrated by the funnel plots and Egger’s test.

Conclusions: This meta-analysis suggests that patients with OAG may have higher TNF-α levels compared with the control subjects, and the TNF-α −308G/A polymorphism is significantly associated with the risks of high-tension glaucoma. Since potential confounders could not be ruled out completely, further studies are needed to confirm these results.

Introduction

Glaucoma, characterized by optic neuropathy and a progressive loss of the retinal ganglion cells (RGCs), is the second leading cause of blindness worldwide [1]. Glaucoma is classified into primary and secondary glaucoma according to the etiology and aqueous humor dynamics. Based on the anatomy of the anterior chamber, primary glaucoma is further classified as primary open angle glaucoma (POAG) and primary angle closure glaucoma (PACG). High-tension glaucoma (HTG) and normal tension glaucoma (NTG) are types of POAG, while exfoliation glaucoma (EXG) is a type of secondary glaucoma. Open angle glaucoma (OAG), including HTG, NTG, and EXG, is the most common form of glaucomas, accounting for about 90% of glaucoma cases in the United States [2]. The etiology of OAG is complex, and not fully understood at present. Although elevated intraocular pressure (IOP) has been recognized as the major risk factor for OAG, elevated IOP is not a diagnostic factor [2]. Furthermore, currently available drugs primarily aim to lower the IOP; however, the disease progression may continue despite significant IOP reduction [3]. Apparently, other factors other than IOP may also play important roles in the pathogenesis of glaucomatous optic neuropathy.

Tumor necrosis factor-alpha (TNF-α) is a cytokine that belongs to the TNF superfamily of 19 different protein ligands [4]. TNF-α is a proinflammatory cytokine with multiple functions in the immune response. Accumulating studies strongly support the involvement of TNF-α in the etiology of glaucoma [3,5]. Ischemic or pressure-loaded glial cells could produce TNF-α, which results in oligodendrocyte death and the subsequent apoptosis of RGCs [6]. The expression of TNF-α and TNF-α receptor-1 (TNF-R1) was upregulated in the retina and the optic nerve head, and the expression of TNF-α and TNF-R1 appeared to parallel the progression of optic nerve degeneration [6-8]. Furthermore, TNF-α was upregulated as a consequence of increasing IOP, and exogenous TNF-α could also lead to loss of oligodendrocytes and a delayed loss of RGCs [9]. The same study also demonstrated that functional blockade of TNF-α with an anti-TNF-α blocking antibody or deletion of the gene encoding TNF-α in a genetically altered mouse model prevented ocular hypertension (OH)-induced oligodendrocyte degeneration and the secondary loss of RGCs [9]. This research suggests that TNF-α may play a key role in glaucomatous neurodegeneration.

The human TNF-α gene is located within the highly polymorphic major histocompatibility complex class III region on chromosome 6p21.3 [10]. Several functional polymorphisms in the promoter region of the TNF-α gene have been identified, and have been related to the risks of glaucoma. The intensively studied polymorphism is characterized by a G to A substitution at position −308 (−308G/A, rs1800629), which could increase transcription six- to sevenfold [11]. Other polymorphisms including −238G/A (rs361525), −863C/A (rs645836), and −857C/T (rs1799724) have also been studied. However, the results have been highly inconsistent. Some studies suggested that the TNF-α −308 A allele was significantly associated with the increased risk of OAG [12-15], while other studies showed no significant association [16-18], and even reduced risks of OAG [19,20]. Furthermore, to investigate the roles of TNF-α and risks of glaucoma, several recent studies detected the TNF-α level in the aqueous humor (AH) in patients with glaucoma and the control subjects. Similarly, some studies showed that the TNF-α concentrations in AH and/or the ratio of TNF-α positive AH samples in patients with glaucoma was higher than those in controls [21,22], which was disclaimed by other studies [23].

In the current study, we performed a meta-analysis of all eligible studies, to provide a more accurate estimate of the association of the TNF-α gene (−308G/A, −238G/A, −863C/A, and −857C/T) polymorphisms and the risks of OAG. We also performed a meta-analysis to compare the TNF-α level in AH between patients with glaucoma and the control subjects.

Methods

Literature and search strategy

A computerized literature search was performed for the relevant available studies from three databases including PubMed, ISI Web of Science, and Embase. The search strategy for identifying the relevant studies involved combinations of the following keywords: “tumor necrosis factor-α” or “TNF-α” and “glaucoma.” The reference lists of review articles, clinical trials, and meta-analyses were also hand-searched to collect other relevant studies. If more than one article was published using the same case series, only the study with the largest sample size was selected. The literature search was updated on August 30, 2012.

The inclusion criteria were as follows: (1) evaluating the associations between TNF-α (−308G/A, −238G/A, −863C/A, or −857C/T) polymorphisms and the risk of OAG or comparing the TNF-α level in AH of patients with glaucoma and that of the control subjects, (2) case-control design, and (3) gene polymorphisms studies must provide sufficient data for calculating the odds ratio (OR) and the corresponding 95% confidence interval (95% CI). All identified studies were carefully reviewed independently by two investigators to determine whether an individual study was eligible for inclusion in this meta-analysis

Data extraction

Data were extracted independently by two investigators who reached a consensus on all of the items. The following information was extracted from each study: (1) name of the first author, (2) year of publication, (3) country of origin, (4) ethnicity of the study population, (5) number of cases and controls, (6) gender and age of enrolled subjects, and (7) number of genotypes in cases and controls.

Statistical analysis

Genotypic frequency for the TNF-α gene polymorphisms was tested for deviation from Hardy–Weinberg equilibrium (HWE) in the control subjects using the chi-square goodness of fit. The associations between the TNF-α (−308G/A, −238G/A, −863C/A, and −857C/T) polymorphisms and OAG risks were estimated by calculating pooled ORs and 95% CI. The significance of the pooled effect size was determined with the Z test. Heterogeneity among studies was assessed using the Q test as well as the I2 statistic, which was documented for the percentage of the observed between-study variability due to heterogeneity rather than chance [24]. A significant Q-statistic (p<0.10) indicated heterogeneity across studies. The DerSimonian and Laird random effect model (REM) was used if heterogeneity existed; otherwise, the Mantel-Haenszel fixed effect model (FEM) was used [24]. Subgroup analyses were stratified by ethnicity and the type of OAG. Sensitivity analysis was conducted by removing an individual study each time to check whether a single study could bias the overall estimate [25]. Begg’s funnel plots and Egger’s regression test were performed to assess the potential publication bias [26].

A pooled standardized mean difference (SMD), together with 95% CI, was used to compare the TNF level in AH between the patients with OAG and the controls. The SMD was chosen because the TNF-α concentrations were measured using different methods in the included studies. The heterogeneity test was also performed by using the Q test and the I2 statistic as performed in other studies [27,28]. In cases of no statistical heterogeneity, the FEM was used for SMD meta-analysis; otherwise, the REM was chosen as the appropriate choice. Begg’s funnel plots and Egger’s linear regression were used to assess evidence for publication bias.

All statistical analysis was performed using STATA version 11 (StataCorp LP, College Station, TX). Two-sided p values less than 0.05 were considered statistically significant.

Results

Characteristics of studies included for the association between tumor necrosis factor alpha gene polymorphisms and open angle glaucoma risks

Eleven studies investigating the association between TNF-α gene polymorphisms and OAG risks were retrieved. Among these studies, Fan et al. recruited three forms of patients with OAG [19], while Razeghinejad et al. and Funayama et al. investigated two forms of OAG [13,18]. Three studies investigated the association of OAG risks with TNF-α (−308G/A) and TNF-α (−238G/A) [16,17,29], one study investigated the association of OAG with TNF-α (−308G/A) and TNF-α (−857C/T) [19], and one study investigated the association of OAG with three single nucleotide polymorphisms, TNF-α (−308G/A), TNF-α (−857C/T), and TNF-α (−863C/A) [18]. All other studies investigated a single nucleotide polymorphism in one form of OAG [12,14,15,20,30]. Finally, 14 (1,182 cases and 3,003 controls), five (808 cases and 1,039 controls), three (645 cases and 666 controls), and three studies (404 cases and 625 controls) were included in the analyses for the association between TNF-α −308G/A, −857C/T, −863C/A, and −238G/A polymorphisms and risks of OAG, respectively. The studies used blood samples for DNA extraction, while polymerase chain reaction-restriction fragments length polymorphism, TaqMan, or DNA sequencing methods were used for genotyping. Genotype distribution in control groups were in HWE in most of the studies except two studies for the TNF-α (−308G/A) polymorphism [14,19]. Details of the studies are shown in Table 1.

Quantitative data synthesis

Results of pooled analysis on the associations between TNF-α −308G/A polymorphisms and OAG risks are shown in Table 2. As shown in Table 2, when all the studies were included in the analysis, no significant association was observed between the TNF-α −308G/A polymorphism and the risk of OAG (A versus G: OR=1.379, 95% CI=0.877–2.170; AA/AG versus GG: OR=1.421, 95% CI=0.907–2.226). Subgroup analysis stratified by type of OAG and ethnicity further confirmed the null association except the homozygous genotype comparison for EXG. However, when the two studies in which the genotype distribution in the control groups deviated from HWE were excluded [14,19], a significant association was observed between TNF-α −308G/A polymorphisms and the risks of HTG (A versus G: OR=1.660, 95% CI=1.033–2.667; AA/AG versus GG: OR=1.713, 95% CI=1.10–2.651), but not with NTG (A versus G: OR=1.005, 95% CI=0.321–3.140; AA/AG versus GG: OR=1.005, 95% CI=0.319–3.165) or EXG (A versus G: OR=1.181, 95% CI=0.373–3.732; AA/AG versus GG: OR=1.130, 95% CI=0.346–3.689). Furthermore, ethnicity-stratified analysis revealed that a significant association also existed in Asians (A versus G: OR=1.947, 95% CI=1.097–3.456; AA/AG versus GG: OR=1.949, 95% CI=1.140–3.332) (Table 2 and Figure 1).

The pooled results on the associations between the TNF-α (−238G/A, −863C/A, and −857C/T) polymorphisms and the OAG risks are shown in Table 3. The combined results showed that none of these gene polymorphisms was significantly associated with the risks of OAG (for −857C/T, T versus C, OR=1.055, 95% CI=0.931 −1.216; TT/CT versus CC, OR=0.984, 95% CI=0.794–1.219; for −863C/A, A versus C, OR=0.928, 95% CI=0.636–1.353; AA/AC versus CC, OR=0.871, 95% CI=0.618–1.227; for −238G/A, OR=1.088, 95% CI=0.691–1.941; AA/AG versus GG, OR=1.076, 95% CI=0.681–1.701). Due to the limited studies available, we did not conduct a sub-group analysis for the TNF-α (−238G/A, −863C/A, −857C/T) polymorphisms.

Sensitivity analysis and publication bias

Sensitivity analysis was performed by sequential omission of individual studies in every comparison, and the data showed that no study significantly influenced the pooled effects by omitting any study. After studies that deviated from HWE in controls were excluded, no other studies were found to significantly influence the pooled effects in each genetic model. Begg’s funnel plots were generated to assess publication bias. Egger’s test was performed to statistically evaluate funnel plot symmetry. The results suggested no publication bias for the association of the TNF-α (−308G/A, −238G/A, −863C/A, and −857C/T) polymorphisms and the OAG risks (PEgger test>0.05; Figure 2).

Pooled results of tumor necrosis factor alpha level in aqueous humor of patients with glaucoma and the control subjects

Five studies in which the TNF-α levels in AH between patients with glaucoma and control subjects were compared were retrieved. Enzyme-linked immunosorbent assay, singleplex bead immunoassay analysis, multiplex bead immunoassay, and microparticle-based immunoassays were used to measure the TNF-α concentration in the AH in different studies. The TNF-α level was below the detection limit in one study [31], while Takai et al. recruited two forms of patients with OAG [23], and Sawada et al. enrolled three forms of patients with OAG [32]. Some studies provided the mean and standard deviation of TNF-α levels in patients with OAG and controls, while other studies reported the TNF-α positive and negative AH samples in patients with OAG and controls. Details of the retrieved studies are shown in Table 4.

The results of the meta-analysis for the comparisons of TNF-α levels in AH between the patients with glaucoma and the controls are shown in Figure 3. The pooled results of six studies demonstrated that patients with OAG had higher TNF-α levels in AH than the control subjects (SMD=0.517, 95% CI=0.207–0.826, p=0.001). No significant heterogeneity was observed (I2=43.3%, p=0.116). Egger’s test and Begg’s funnel plot were applied for comparison to assess the publication bias of the literature, and no possibility of publication bias for this test was observed (p=0.191).

Discussion

Glaucoma is the second leading cause of blindness worldwide. Growing evidence supports the important role of TNF-α as a mediator of neurodegeneration in OAG [3,5]. Some of these studies focused on the association between the TNF-α gene polymorphisms and OAG risks, in which several gene polymorphisms in the TNF-α gene promoter region including TNF-α (−308G/A, −238G/A, −863C/A, and −857C/T) have been investigated. However, these epidemiological studies obtained contradictory results, which might be related to factors such as the sample size of the study, the ethnicity of the enrolled subjects, and the type of OAG. Considering the potential important roles of TNF-α in the etiology of OAG, a timely meta-analysis was conducted to provide an appropriate approach to obtain a more definitive conclusion.

In the current study, 14, five, three, and three studies were included in the analyses for the associations between TNF-α −308G/A, −857C/T, −863C/A, and −238G/A polymorphisms and the risks of OAG, respectively. When all the studies were included, the combined results showed no significant association between the TNF-α −308G/A polymorphism and the risks of OAG. However, when the studies in which the genotype distribution in the control groups deviated from HWE were excluded [14,19], the pooled results clearly showed that the TNF-α −308G/A polymorphisms were significantly associated with the risks of HTG, but not with the risks of NTG or EXG. Furthermore, ethnicity-stratified subgroup analysis revealed that significant association also existed in Asians. Regarding the other genetic polymorphisms including −238G/A, −863C/A, and −857C/T, the results of the current meta-analysis demonstrated no significant association.

In a recently published meta-analysis, Yu et al. investigated the association between the TNF-α −308G/A polymorphism and the risks of glaucoma [33]. In that study, eight studies investigating the association between the TNF-α −308G/A polymorphism and the risks of glaucoma were included; those were included in the current study, in which five additional newly published studies were included [12,19,29]. Importantly, the meta-analysis by Yu et al. did not include subgroup analysis for HTG and NTG, which are apparently different, although both are types of POAG. HTG is characterized by high IOP, while the IOP of NTG is usually within normal levels. Findings of ongoing in vivo studies support that TNF-α and TNF-R1 are upregulated following experimental elevation of IOP, which indicates that a potential link may exist between TNF-α and IOP [5]. Furthermore, the only two studies in which the association between TNF-α −308G/A polymorphisms and the NTG risks was investigated reported null association [18,19]. Therefore, it is not surprising that the TNF-α −308G/A gene polymorphism was significantly associated with HTG risk, but not with NTG risks. EXG is another type of OAG with quite different characteristics from POAG. EXG is the most frequently reported type of secondary glaucoma, which is caused by deposition of exfoliation material and liberated iris pigment in the trabecular meshwork leading to elevation of IOP [34]. Thus, EXG may be significantly associated with other genetic polymorphisms such as the lysyl oxidase-like 1 (LOXL1) gene [35].

To investigate whether TNF-α might play important roles in the etiology of OAG, comparison of the TNF-α level between patients with OAG and controls would provide direct evidence. Therefore, in the present study, in addition to exploring the association between TNF-α gene polymorphisms and the risks of OAG, we summarized the studies in which the TNF-α levels in the AH were detected in patients with OAG as well as in control subjects. Some studies provided the mean and standard deviation of TNF-α levels in patients with OAG and controls, while other studies reported the positive and negative AH samples of patients with OAG and controls (Table 4). The pooled results of six individual studies showed that patients with OAG had higher TNF-α levels in AH than the control subjects (SMD=0.517, 95% CI=0.207–0.826, p=0.001; Figure 3). Other studies showed that patients with glaucoma were more likely to have detectable levels of TNF-α in their AH [21]. These data support the theory that TNF-α might play important roles in the etiology of OAG. However, most of these studies enrolled Japanese patients. Thus, whether similar results could be obtained in people in other countries remains unclear.

Despite the clear strengths of our study such as the larger sample size compared with the previous individual ones, it does have some limitations. First, the present meta-analysis was based primarily on unadjusted effect estimates and CIs for the association study between TNF-α gene polymorphisms and the OAG risks; thus, the effect estimates were relatively imprecise. Second, OAG is a multifactor disease, which may be related to genetic and environmental factors. However, the effects of gene–gene and gene–environment interactions were not addressed in this meta-analysis, and thus the potential roles of these gene polymorphisms may be masked or magnified by other gene–gene/gene–environment interactions. Third, although the funnel plot and Egger’s test showed no publication bias, selection bias may also exist because only published studies in the three selected databases were retrieved.

In summary, the current meta-analysis suggests that the TNF-α −308G/A gene polymorphism is significantly associated with the risks of HTG, but not with NTG or EXG, while the TNF-α (−238G/A, −863C/A, −857C/T) polymorphisms are not significantly associated with OAG risks. Furthermore, the pooled results of six individual studies indicate that patients with OAG may have higher TNF-α levels in AH than the control subjects, which supports the important roles of TNF-α in the development of OAG. Since potential confounders could not be ruled out completely, further studies are needed to confirm these results.

References

  1. Quigley HA, Broman AT. The number of people with glaucoma worldwide in 2010 and 2020. Br J Ophthalmol. 2006; 90:262-7. [PMID: 16488940]
  2. Distelhorst JS, Hughes GM. Open-angle glaucoma. Am Fam Physician. 2003; 67:1937-44. [PMID: 12751655]
  3. Agarwal R, Agarwal P. Glaucomatous neurodegeneration: An eye on tumor necrosis factor-alpha. Indian J Ophthalmol. 2012; 60:255 [PMID: 22824592]
  4. Dempsey PW, Doyle SE, He JQ, Cheng G. The signaling adaptors and pathways activated by TNF superfamily. Cytokine Growth Factor Rev. 2003; 14:193-209. [PMID: 12787559]
  5. Tezel G. TNF-alpha signaling in glaucomatous neurodegeneration. Prog Brain Res. 2008; 173:409-21. [PMID: 18929124]
  6. Tezel G, Wax MB. Increased production of tumor necrosis factor-alpha by glial cells exposed to simulated ischemia or elevated hydrostatic pressure induces apoptosis in cocultured retinal ganglion cells. J Neurosci. 2000; 20:8693-700. [PMID: 11102475]
  7. Yuan L, Neufeld AH. Tumor necrosis factor-alpha: a potentially neurodestructive cytokine produced by glia in the human glaucomatous optic nerve head. Glia. 2000; 32:42-50. [PMID: 10975909]
  8. Yan X, Tezel G, Wax MB, Edward DP. Matrix metalloproteinases and tumor necrosis factor alpha in glaucomatous optic nerve head. Arch Ophthalmol. 2000; 118:666-73. [PMID: 10815159]
  9. Nakazawa T, Nakazawa C, Matsubara A, Noda K, Hisatomi T, She H, Michaud N, Hafezi-Moghadam A, Miller JW, Benowitz LI. Tumor necrosis factor-alpha mediates oligodendrocyte death and delayed retinal ganglion cell loss in a mouse model of glaucoma. J Neurosci. 2006; 26:12633-41. [PMID: 17151265]
  10. Hajeer AH, Hutchinson IV. Influence of TNFalpha gene polymorphisms on TNFalpha production and disease. Hum Immunol. 2001; 62:1191-9. [PMID: 11704281]
  11. Agarwal P, Oldenburg MC, Czarneski JE, Morse RM, Hameed MR, Cohen S, Fernandes H. Comparison study for identifying promoter allelic polymorphism in interleukin 10 and tumor necrosis factor alpha genes. Diagn Mol Pathol. 2000; 9:158-64. [PMID: 10976723]
  12. Al-Dabbagh NM, Al-Dohayan N, Al-Asmari A, Arfin M, Tariq M. Association of TNF-α and TNF-β Gene Polymorphisms with Primary Open Angle and Primary Angle Closure Glaucoma. In: Tomas Kubena, editor. The Mystery of Glaucoma. InTech; 2011. p.229-256.
  13. Razeghinejad MR, Rahat F, Kamali-Sarvestani E. Association of TNFA −308 G/A and TNFRI +36 A/G gene polymorphisms with glaucoma. Ophthalmic Res. 2009; 42:118-24. [PMID: 19556827]
  14. Khan MI, Micheal S, Rana N, Akhtar F, den Hollander AI, Ahmed A, Qamar R. Association of tumor necrosis factor alpha gene polymorphism G-308A with pseudoexfoliative glaucoma in the Pakistani population. Mol Vis. 2009; 15:2861-7. [PMID: 20029655]
  15. Lin HJ, Tsai FJ, Chen WC, Shi YR, Hsu Y, Tsai SW. Association of tumour necrosis factor alpha −308 gene polymorphism with primary open-angle glaucoma in Chinese. Eye (Lond). 2003; 17:31-4. [PMID: 12579167]
  16. Mossböck G, Renner W, El-Shabrawi Y, Faschinger C, Schmut O, Wedrich A, Zimmermann C, Weger M. TNF-alpha −308 G>A and −238 G>A polymorphisms are not major risk factors in Caucasian patients with exfoliation glaucoma. Mol Vis. 2009; 15:518-22. [PMID: 19279689]
  17. Mossböck G, Weger M, Moray M, Renner W, Haller-Schober EM, Mattes D, Schmut O, Wegscheider B, El-Shabrawi Y. TNF-alpha promoter polymorphisms and primary open-angle glaucoma. Eye (Lond). 2006; 20:1040-3. [PMID: 16138112]
  18. Funayama T, Ishikawa K, Ohtake Y, Tanino T, Kurosaka D, Kimura I, Suzuki K, Ideta H, Nakamoto K, Yasuda N, Fujimaki T, Murakami A, Asaoka R, Hotta Y, Tanihara H, Kanamoto T, Mishima H, Fukuchi T, Abe H, Iwata T, Shimada N, Kudoh J, Shimizu N, Mashima Y. Variants in optineurin gene and their association with tumor necrosis factor-alpha polymorphisms in Japanese patients with glaucoma. Invest Ophthalmol Vis Sci. 2004; 45:4359-67. [PMID: 15557444]
  19. Fan BJ, Liu K, Wang DY, Tham CC, Tam PO, Lam DS, Pang CP. Association of polymorphisms of tumor necrosis factor and tumor protein p53 with primary open-angle glaucoma. Invest Ophthalmol Vis Sci. 2010; 51:4110-6. [PMID: 20357201]
  20. Tekeli O, Turacli ME, Egin Y, Akar N, Elhan AH. Tumor necrosis factor alpha-308 gene polymorphism and pseudoexfoliation glaucoma. Mol Vis. 2008; 14:1815-8. [PMID: 18852869]
  21. Chua J, Vania M, Cheung CM, Ang M, Chee SP, Yang H, Li J, Wong TT. Expression profile of inflammatory cytokines in aqueous from glaucomatous eyes. Mol Vis. 2012; 18:431-8. [PMID: 22355254]
  22. Balaiya S, Edwards J, Tillis T, Khetpal V, Chalam KV. Tumor necrosis factor-alpha (TNF-alpha) levels in aqueous humor of primary open angle glaucoma. Clin Ophthalmol. 2011; 5:553-6. [PMID: 21607023]
  23. Takai Y, Tanito M, Ohira A. Multiplex cytokine analysis of aqueous humor in eyes with primary open-angle glaucoma, exfoliation glaucoma, and cataract. Invest Ophthalmol Vis Sci. 2012; 53:241-7. [PMID: 22159018]
  24. Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med. 2002; 21:1539-58. [PMID: 12111919]
  25. Tobias A. Assessing the influence of a single study in the meta-analysis estimate. Stata Tech Bull. 1999; 47:15-7.
  26. Harbord RM, Egger M, Sterne JA. A modified test for small-study effects in meta-analyses of controlled trials with binary endpoints. Stat Med. 2006; 25:3443-57. [PMID: 16345038]
  27. Zeng T, Guo FF, Zhang CL, Song FY, Zhao XL, Xie KQ. A meta-analysis of randomized, double-blind, placebo-controlled trials for the effects of garlic on serum lipid profiles. J Sci Food Agric. 2012; 92:1892-902. [PMID: 22234974]
  28. Cui G, Wang H, Li R, Zhang L, Li Z, Wang Y, Hui R, Ding H, Wang D. Polymorphism of tumor necrosis factor alpha (TNF-alpha) gene promoter, circulating TNF-alpha level, and cardiovascular risk factor for ischemic stroke. J Neuroinflammation. 2012; 9:235 [PMID: 23050663]
  29. Bozkurt B, Mesci L, Irkec M, Ozdag BB, Sanal O, Arslan U, Ersoy F, Tezcan I. Association of tumour necrosis factor-alpha −308 G/A polymorphism with primary open-angle glaucoma. Clin Experiment Ophthalmol. 2012; 40:e156-62. [PMID: 21575121]
  30. Wang CY, Shen YC, Wei LC, Lin KH, Feng SC, Yang YY, Chiu CH, Tsai HY. Polymorphism in the TNF-alpha(−863) locus associated with reduced risk of primary open angle glaucoma. Mol Vis. 2012; 18:779-85. [PMID: 22509108]
  31. Kuchtey J, Rezaei KA, Jaru-Ampornpan P, Sternberg P, , Jr Kuchtey RW. Multiplex cytokine analysis reveals elevated concentration of interleukin-8 in glaucomatous aqueous humor. Invest Ophthalmol Vis Sci. 2010; 51:6441-7. [PMID: 20592224]
  32. Sawada H, Fukuchi T, Tanaka T, Abe H. Tumor necrosis factor-alpha concentrations in the aqueous humor of patients with glaucoma. Invest Ophthalmol Vis Sci. 2010; 51:903-6. [PMID: 19737888]
  33. Yu QQ, Yao Y. A detailed meta-analysis shows no association between TNF-alpha −308G/A polymorphism and different forms of glaucoma. Ophthalmic Res. 2012; 47:47-51. [PMID: 21701244]
  34. Ritch R, Schlötzer-Schrehardt U, Konstas AGP. Why is glaucoma associated with exfoliation syndrome? Prog Retin Eye Res. 2003; 22:253-75. [PMID: 12852486]
  35. Schlötzer-Schrehardt U. Genetics and genomics of pseudoexfoliation syndrome/glaucoma. Middle East Afr J Ophthalmol.. 2011; 18:30-6. [PMID: 21572731]