Ethnicity-Stratified Analysis of the Association between TNF-α Genetic Polymorphisms and Acute Kidney Injury: A Systematic Review and Meta-Analysis
Background. Several studies have reported conflicting findings regarding the association between tumor necrosis factor-alpha (TNF-α) genetic polymorphisms and acute kidney injury (AKI). Therefore, we performed this meta-analysis to further investigate whether TNF-α variants are related to AKI susceptibility. Methods. A comprehensive search of observational studies on the association of TNF-α polymorphism with AKI susceptibility was conducted in the PubMed, Cochrane, and Embase databases through February 10, 2020. Pooled odds ratios (ORs) and 95% corresponding confidence intervals (95% CIs) were analyzed to evaluate the strength of the relationship. Results. A total of 8 studies involving 6694 patients (2559 cases and 4135 controls) were included. Pooled analysis showed a trend of increased risk between the TNF-α rs1800629 variant and AKI (A vs. G: ) among the overall population. Ethnicity-stratified analysis indicated that the TNF-α rs1800629 variant was a risk factor for Asians () while it is not for Caucasians (). Additionally, we also found that TNF-α rs1799964 polymorphism was observed to have a significant relationship with AKI risk in Asian patients (C vs. T, ). Conclusions. The TNF rs1800629 polymorphism exhibited a trend toward AKI susceptibility with ethnic differences. The relationship was found to be significant among the Asian population, but not among those of Caucasian origin. Additionally, the TNF-α rs1799964 polymorphism was also related to a significantly increased risk of AKI in Asians.
Acute kidney injury (AKI), characterized by a rapid decline in kidney function within two days, has become one of the most common syndromes in hospital settings worldwide . AKI occurs in approximately 10–30% of patients admitted to the hospital. The rate is even as high as 50% when it comes to the high-risk patients [2, 3]. Furthermore, it is also associated with high mortality (10-30%) and a high occurrence of chronic kidney disease (24.6%) [2, 4]. The key step for preventing AKI is to identify high-risk patients . Besides clinical risk factors [6, 7], genetic factors have been reported relative to interindividual differences in susceptibility to AKI [8, 9].
Tumor necrosis factor-α (TNF-α) could initiate the inflammatory cascade and induces the production of inflammatory mediators [10, 11]. Several TNF-α single nucleotide polymorphisms (SNPs), including rs1800629 and rs1799964, had been reported to associate with altered TNF-α expression [12, 13]. Some studies had analyzed the association between two TNF-α promoter SNPs (rs1800629 and rs1799964) and AKI susceptibility [14–18]. Dalboni et al.  and McBride et al. [14, 15] conducted studies in Caucasians and found that TNF-α rs1800629 had no association with susceptibility to AKI, while studies conducted in Asians indicated that both TNF-α rs1800629 and rs1799964 had a positive association with AKI [16, 18].
Even though several studies made an attempt to explore the genetic association between TNF-α and AKI risk [14–21], these studies had inconsistent results or small sample size. Therefore, we aimed to conduct this study to assess the genetic association as well as the ethnic difference of two TNF-α SNPs (rs1800629 and rs1799964) with AKI risk.
This work adheres to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines, and the protocol of this meta-analysis has been published on PROSPERO (https://www.crd.york.ac.uk/prospero/), of which the registration number is CRD42020191747.
2.1. Literature Search
A comprehensive strategy was used to search for studies about the correlation between TNF-α polymorphisms and AKI from the PubMed, Embase, and Cochrane databases. Searching terms or keywords were “Polymorphism, Genetic,” “Genetic,” “TNF-α,” “Tumor necrosis factor,” “AKI,” and “acute kidney injury.” The searched articles and abstracts were limited to those written in English before February 10, 2020 (for details, see Additional file 1).
2.2. Inclusion and Exclusion Criteria
Studies were included according to the following criteria: (1) evaluation of the association between TNF polymorphisms and AKI, (2) studies focusing on humans, and (3) studies with detailed genotype data. Studies were excluded if they were classified as (1) duplication of publications; (2) commentaries, reviews, editorials, animal studies, or case reports; (3) studies not relevant to TNF-α polymorphism and AKI risk; (4) studies about kidney transplant; or (5) studies lacking control groups or without detailed genotype data.
2.3. Quality Assessment
We used the Newcastle-Ottawa Scale (NOS) to evaluate the quality of the included studies. According to the scoring results, the included articles were classified as low, medium, and high quality, with scores of “1-3,” “4-5,” and “6-9,” respectively.
2.4. Data Extraction from the Selected Article
According to the inclusion and exclusion criteria, two reviewers selected articles by screening the title, abstract, and full text independently. Then, these two reviewers analyzed each paper and extracted detailed information. The following information was extracted: name of first author, publication year, SNPs, country, ethnicity, total sample size, numbers of cases and controls, and numbers of genotype frequencies. Any divergence was settled by a consensus of all authors.
2.5. Statistical Analysis
The Hardy–Weinberg equilibrium (HWE) was calculated for control groups of each study by using the test. The study population was considered a significant departure from genetic balance if value < 0.05. Pooled ORs with 95% CIs were calculated for the allelic model (rs1800629: A vs. G, rs1799964: T vs. C), heterozygous model (rs1800629: GA vs. GG, rs1799964: TC vs. TT), homozygous model (rs1800629: AA vs. GG, rs1799964: CC vs. TT), dominant model (rs1800629: AA+GA vs. GG, rs1799964: CC+TC vs. TT), and recessive model (rs1800629: AA vs. GG+GA, rs1799964: CC vs. TT+TC). Heterogeneity was evaluated by the statistic and statistic . When was >50%, a random-effect model was adopted. Otherwise, the fixed-effect model was used . We also conducted a meta-regression analysis to determine the source of heterogeneity using a restricted maximum likelihood model with the Knapp-Hartung method [24, 25]. To assess the effect of ethnicity difference, we performed a subgroup analysis by ethnicity, categorized as Caucasian and Asian. Influence analysis was performed to recompute the pooled risk estimates for the remaining studies by removing each study successively.
Both Begg’s test and Egger’s test were used to assess publication bias . Publication bias was considered with a significant difference when . STATA 12.0 software (Stata Corporation, TX, USA) was performed to do all statistical analyses.
3.1. Characteristics of the Studies
The literature search strategy identified 198 related studies from the PubMed, Embase, and Cochrane databases. The literature selection process is shown in the schematic representation in Figure 1. A total of 181 articles were excluded by identifying the titles and abstracts, of which 23 were duplicates, 97 had no relation to this topic, 19 were related to kidney transplant, and 42 were commentaries, reviews, editorials, animal studies, or case reports. The remaining 17 studies were then reviewed in full text, and 9 studies were excluded, among which 2 lacked control groups and the others lacked detailed genotype data. Finally, 8 eligible studies [14, 15–21] including 6694 patients (2559 cases and 4135 controls) were chosen to be analyzed (Table 1). Within all of the included studies, there were 7 studies related to TNF-α rs1800629 polymorphism and 3 studies related to TNF-α rs1799964 polymorphism. Regarding the TNF-α rs1800629 studies, the subjects in 5 studies [14, 15, 19–21] were Caucasians and the subjects in the other 2 were Asians [16, 17]. Regarding the TNF-α rs1799964 studies, all the 3 included studies were conducted in the Asian population [16–18]. All the included studies were assessed as high quality (, see Additional file 2), and the genotyping distributions were genetically balanced ().
3.2. Meta-Analysis of TNF-α rs1800629
In the pooled analysis of 7 eligible studies, significant heterogeneity was found in the allelic model (A vs. G), heterozygous model (AG vs. GG), and dominant model (AG+AA vs. GG); therefore, the random-effect model was used for these genetic models. Although there was no statistical significance found in the allelic model (), an increasing trend of AKI risk could be seen, as shown in Figure 2, while a significant association can be seen in the heterozygous comparison () and dominant comparison (). Moreover, homozygous comparison () and recessive comparison () had no significance (see Table 2).
3.3. Meta-regression of TNF-α rs1800629
In order to explore the possible sources of heterogeneity, meta-regression analysis was performed. The result indicated that ethnicity had a significant association with outcomes at posttest (; ; , see Figure 3). To sum up, nearly most of heterogeneity could be explained by the ethnicity variance between included studies. However, sample size, case ratio, and sexual ratio were not significantly related to outcomes according to the results of the meta-regression (data not shown).
3.4. Ethnicity-Stratified Analysis of TNF-α rs1800629
We conducted a subgroup analysis by ethnicity. The ethnicity-stratified analysis had no significant heterogeneity for both Asian and Caucasian population analyses; thus, a fixed-effect model was employed. In Asians, all the five genetic models could be detected to significantly associate with AKI susceptibility (see Table 2). However, in Caucasians, there was no significant correlation in any genetic model (see Table 2).
3.5. Meta-Analysis of TNF-α rs1799964
A total of 3 studies were taken into the meta-analysis of the correlation between the TNF-α rs1799964 and AKI susceptibility. All the included studies were carried out in Asian populations; therefore, no significant heterogeneity was identified. The pooled analysis showed that the rs1799964 T allele was associated with the increased risk of AKI in all the genetic models (see Figure 4, Table 2, allelic comparison: ; heterozygous comparison: ; homozygous comparison: ; dominant comparison: ; and recessive comparison: ).
3.6. Influence Analysis
Influence analysis recomputing for the remaining studies by removing each study successively contributed to little change on the pooled results, from to (see Figure 5).
3.7. Publication Bias
Publication bias was evaluated by Begg’s funnel plot and Egger’s regression test. The results showed that no significant publication bias was found in the TNF-α rs1800629 variant (Begg’s test: ; Egger’s test: ). Symmetrical funnel plots are shown in Figure 6. Owing to the limited included studies, no publication bias test was conducted for the TNF-α rs1799964.
3.8. SNP Allele Frequency in Populations
As shown above, TNF-α rs1800629 was significantly associated with AKI susceptibility in Asians and but not in Caucasians. We next investigated the allele frequencies of TNF-α rs1800629 between ethnicity (see Table 3). In our meta-analysis, the rs1800629 A allele frequency in the control group (0.149 for Caucasians, 0.063 for Asians) was similar with that in the 1000 Genomes Project (0.134 for Europeans, 0.059 for Asians). However, TNF-α rs1799964 C allele frequency in the control group (0.133 for Asians) was different from that in the 1000 Genomes Project (0.196 for Asians).
To our knowledge, the present study is the first global systematic meta-analysis to evaluate the genetic relationship between TNF-α polymorphisms and AKI susceptibility. The results indicated that TNF-α rs1800629 had a trend of increased risk of AKI for all patients. What is more, this association was significant in Asian patients, but was not significant in Caucasian patients by subgroup analysis. Additionally, because of the limited number of included studies, TNF-α rs1799964 was shown to be related to an increased risk for AKI only in Asian patients.
In the previous systematic reviews, both Vilander et al.  and Larach et al.  indicated that the relationship between TNF-α polymorphisms and AKI susceptibility was controversial. In comparison, both of them did not assess the effect of TNF-α polymorphisms or conduct a detailed analysis of the conflicting results. Consequently, we performed this meta-analysis to quantificationally evaluate the relationship between TNF-α polymorphisms and AKI susceptibility. More importantly, we found that TNF-α polymorphisms had a trend of increasing the risk of AKI susceptibility with ethnic differences. In the present meta-analysis, the association between TNF-α rs1800629 and AKI susceptibility was not significant for the overall populations. Meta-regression and ethnicity-stratified analyses were performed to investigate the main source of heterogeneity. The results showed that the association of TNF-α rs1800629 with AKI susceptibility varied by ethnicity. What is more, the subgroup results indicated that the rs1800629 GA/AA genotype had a higher risk of AKI among Asians but not among Caucasians. Several possible explanations may account for these findings. One of the reasons is that the studies conducted in Caucasians mainly included critically ill patients from the intensive care unit with severe infection [14, 19–21] or underwent cardiac surgery . The patients included in Asian studies underwent coronary artery intervention  or were pediatric patients , with a lower incidence of comorbidity or better baseline renal function. Another possible explanation is that the distribution of TNF-α rs1800629 alleles is population dependent. In the present meta-analysis, the frequencies of the rs1800629 A allele in the control group were 6.26% for Asians [16, 17] and 14.90% for Caucasians [14, 15, 19–21]. This result was similar to that of the 1000 Genomes Project (as shown in Table 3) and that reported by Zhang et al. (8.20% for Asians and 15.70% for Caucasians) . Moreover, several meta-analysis studies also indicated that the rs1800629 A allele had a significant association with susceptibility to osteoarthritis , dilated cardiomyopathy , and chronic obstructive pulmonary disease  among the Asian population, but not among the Caucasian population. These findings indicated that the dependent distribution of alleles could partially explain the susceptibility differences in the ethnic background.
Regarding TNF-α rs1799964, a total of 3 eligible studies conducted in the Asian population including 1204 cases and 1881 controls were analyzed. Although there was a difference of the TNF-α rs1799964 allele frequency observed between our meta-analysis and the 1000 Genomes Project data, there was no deviation from HWE in the control group of all included studies. The difference may be due to the small sample size or the heterogeneity between the included studies. The relationship between TNF-α rs1799964 and AKI susceptibility was found in all the genetic comparisons, suggesting that the rs1799964 variant might lead to a higher risk for AKI in the Asian population, while the relationship needs to be verified in a larger sample size.
Nevertheless, some limitations of this study should not be ignored. First, it is better to further evaluate the effect of TNF-α polymorphisms on AKI susceptibility by screening in the online database, while it is difficult for us to get access to the raw data of AKI genomics to make the validation. Second, AKI is a multifactorial disease. Some confounding factors such as gender, history of nephrotoxic drugs, and heart function may be addressed across the included studies. However, we have no access to get these clinical data to estimate the effect of confounding factors. Third, significant heterogeneity was found in the rs1800629 analysis for the overall population. Thus, we performed meta-regression and ethnicity subgroup analyses to investigate the source of the heterogeneity. Last, the included studies for TNF-α rs1800629 subgroup analysis or rs1799964 analysis were relatively limited; the results of these analyses should be interpreted prudently.
In conclusion, TNF rs1800629 polymorphism had a trend of increasing the risk of AKI with ethnic differences. The trend was significant among Asian patients, but was not among Caucasian patients. Additionally, the TNF-α rs1799964 polymorphism was related to a significantly increased risk of AKI in Asians. However, as AKI is a multifactorial disease, our findings require verification in larger samples and functional studies in the future.
The data used and/or analyzed during the current study are available from the corresponding author on reasonable request.
The funding body played no direct role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Conflicts of Interest
The authors declare no conflict of interest.
Guanzhong Chen, Bowen Liu, Huanqiang Li, and Ziling Mai contributed equally to this work.
This work was supported by the National Science Foundation of China (Grant Numbers: 81670339 and 81970311).
Additional file 1: detailed strategy of the literature search. Additional file 2: quality assessment of the included studies using the Newcastle-Ottawa Scale. (Supplementary Materials)
C. Ronco, R. Bellomo, and J. A. Kellum, “Acute kidney injury,” Lancet, vol. 394, no. 10212, pp. 1949–1964, 2019.View at: Publisher Site | Google Scholar
E. A. J. Hoste, J. A. Kellum, N. M. Selby et al., “Global epidemiology and outcomes of acute kidney injury,” Nature Reviews Nephrology, vol. 14, no. 10, pp. 607–625, 2018.View at: Publisher Site | Google Scholar
E. A. Hoste, S. M. Bagshaw, R. Bellomo et al., “Epidemiology of acute kidney injury in critically ill patients: the multinational AKI-EPI study,” Intensive Care Medicine, vol. 41, no. 8, pp. 1411–1423, 2015.View at: Publisher Site | Google Scholar
K. L. Horne, R. Packington, J. Monaghan, T. Reilly, and N. M. Selby, “Three-year outcomes after acute kidney injury: results of a prospective parallel group cohort study,” BMJ Open, vol. 7, no. 3, article e015316, 2017.View at: Publisher Site | Google Scholar
K. Disease and I. G. Outcomes, “Acute Kidney Injury Work Group: KDIGO clinical practice guideline for acute kidney injury,” Kidney International Supplement, vol. 2, no. 1, pp. 1–138, 2012.View at: Google Scholar
M. E. Matheny, R. A. Miller, T. A. Ikizler et al., “Development of inpatient risk stratification models of acute kidney injury for use in electronic health records,” Medical Decision Making : an international journal of the Society for Medical Decision Making, vol. 30, no. 6, pp. 639–650, 2010.View at: Publisher Site | Google Scholar
S. Kheterpal, K. K. Tremper, M. Heung et al., “Development and validation of an acute kidney injury risk index for patients undergoing general surgery: results from a national data set,” Anesthesiology, vol. 110, no. 3, pp. 505–515, 2009.View at: Publisher Site | Google Scholar
L. M. Vilander, M. A. Kaunisto, and V. Pettila, “Genetic predisposition to acute kidney injury--a systematic review,” BMC Nephrology, vol. 16, no. 1, 2015.View at: Publisher Site | Google Scholar
D. B. Larach, M. C. Engoren, E. M. Schmidt, and M. Heung, “Genetic variants and acute kidney injury: a review of the literature,” Journal of Critical Care, vol. 44, pp. 203–211, 2018.View at: Publisher Site | Google Scholar
F. Balkwill, “Tumour necrosis factor and cancer,” Nature Reviews Cancer, vol. 9, no. 5, pp. 361–371, 2009.View at: Publisher Site | Google Scholar
V. Baud and M. Karin, “Signal transduction by tumor necrosis factor and its relatives,” Trends in Cell Biology, vol. 11, no. 9, pp. 372–377, 2001.View at: Publisher Site | Google Scholar
M. Szkup, E. Chelmecka, A. Lubkowska, A. J. Owczarek, and E. Grochans, “The influence of the TNFα rs1800629 polymorphism on some inflammatory biomarkers in 45-60-year-old women with metabolic syndrome,” Aging, vol. 10, no. 10, pp. 2935–2943, 2018.View at: Publisher Site | Google Scholar
H. Sato, L. Silveira, T. Fingerlin et al., “TNF polymorphism and bronchoalveolar lavage cell TNF-alpha levels in chronic beryllium disease and beryllium sensitization,” The Journal of Allergy and Clinical Immunology, vol. 119, no. 3, pp. 687–696, 2007.View at: Publisher Site | Google Scholar
M. A. Dalboni, B. M. R. Quinto, C. C. Grabulosa et al., “Tumour necrosis factor-α plus interleukin-10 low producer phenotype predicts acute kidney injury and death in intensive care unit patients,” Clinical & Experimental Immunology, vol. 173, no. 2, pp. 242–249, 2013.View at: Publisher Site | Google Scholar
W. T. McBride, P. S. Prasad, M. Armstrong et al., “Cytokine phenotype, genotype, and renal outcomes at cardiac surgery,” Cytokine, vol. 61, no. 1, pp. 275–284, 2013.View at: Publisher Site | Google Scholar
J. He, G. Xie, H. Wu et al., “Association between inflammatory-response gene polymorphisms and risk of acute kidney injury in children,” Bioscience Reports, vol. 38, no. 6, 2018.View at: Publisher Site | Google Scholar
C. F. Chang, T. M. Lu, W. C. Yang, S. J. Lin, C. C. Lin, and M. Y. Chung, “Gene polymorphisms of interleukin-10 and tumor necrosis factor-α are associated with contrast-induced nephropathy,” American Journal of Nephrology, vol. 37, no. 2, pp. 110–117, 2013.View at: Publisher Site | Google Scholar
H. Kamei, Y. Onishi, T. Nakamura, M. Ishigami, and N. Hamajima, “Role of cytokine gene polymorphisms in acute and chronic kidney disease following liver transplantation,” Hepatology International, vol. 10, no. 4, pp. 665–672, 2016.View at: Publisher Site | Google Scholar
P. Cardinal-Fernández, A. Ferruelo, M. el-Assar et al., “Genetic predisposition to acute kidney injury induced by severe sepsis,” Journal of Critical Care, vol. 28, no. 4, pp. 365–370, 2013.View at: Publisher Site | Google Scholar
A. Treszl, P. Tóth-Heyn, I. Kocsis et al., “Interleukin genetic variants and the risk of renal failure in infants with infection,” Pediatric Nephrology, vol. 17, no. 9, pp. 713–717, 2002.View at: Publisher Site | Google Scholar
L. M. Vilander, S. T. Vaara, M. A. Kaunisto, V. Pettilä, and The FINNAKI Study Group, “Common inflammation-related candidate gene variants and acute kidney injury in 2647 critically ill Finnish patients,” Journal of Clinical Medicine, vol. 8, no. 3, p. 342, 2019.View at: Publisher Site | Google Scholar
J. P. T. Higgins and S. G. Thompson, “Quantifying heterogeneity in a meta-analysis,” Statistics in Medicine, vol. 21, no. 11, pp. 1539–1558, 2002.View at: Publisher Site | Google Scholar
R. DerSimonian and N. Laird, “Meta-analysis in clinical trials,” Controlled Clinical Trials, vol. 7, no. 3, pp. 177–188, 1986.View at: Publisher Site | Google Scholar
T. D. Stanley and S. B. Jarrell, “Meta-regression analysis: a quantitative method of literature surveys,” Journal of Economic Surveys, vol. 19, no. 3, pp. 299–308, 2005.View at: Publisher Site | Google Scholar
D. Jackson, J. Bowden, and R. Baker, “Approximate confidence intervals for moment-based estimators of the between-study variance in random effects meta-analysis,” Research Synthesis Methods, vol. 6, no. 4, pp. 372–382, 2015.View at: Publisher Site | Google Scholar
M. Egger, G. D. Smith, M. Schneider, and C. Minder, “Bias in meta-analysis detected by a simple, graphical test,” BMJ, vol. 315, no. 7109, pp. 629–634, 1997.View at: Publisher Site | Google Scholar
H. F. Zhang, S. L. Xie, J. F. Wang, Y. X. Chen, Y. Wang, and T. C. Huang, “Tumor necrosis factor-alpha G-308A gene polymorphism and coronary heart disease susceptibility: an updated meta-analysis,” Thrombosis Research, vol. 127, no. 5, pp. 400–405, 2011.View at: Publisher Site | Google Scholar
J. Chen, Y. Wu, J. Yu, and J. Shen, “Association between tumor necrosis factor alpha rs1800629 polymorphism and risk of osteoarthritis in a Chinese population,” Brazilian Journal of Medical and Biological Research, vol. 51, no. 8, article e7311, 2018.View at: Publisher Site | Google Scholar
M. Salimi Asl, A. Ahmadi, J. Salimian, S. Shohani, S. Azimzadeh Jamalkandi, and M. Ghanei, “TNF-α -308 G/A variant and susceptibility to chronic obstructive pulmonary disease: a systematic review and meta-analysis,” Cytokine, vol. 123, p. 154763, 2019.View at: Publisher Site | Google Scholar