Influence of interleukin-18 polymorphisms on kidney transplantation outcomes: a meta-analysis

Objective: Kidney transplantation (KT) procedures are confronted with adverse outcomes that include allograft failure. Allograft survival are in large part attributed to genetics, which render the recipient susceptible or protected from allograft rejection. The genetics of KT outcomes point to single nucleotide polymorphisms (SNPs) where studies have reported the role of cytokines in allograft survival, one of which is interleukin-18 ( IL-18 ). Reported associations of IL-18 with KT outcomes have been inconsistent. This prompted a meta-analysis to obtain more precise estimates. Methods: From four included articles, we posed two hypotheses about IL-18 SNPs: (1) they are either high in patients (hp) /controls (hc) based on genotype distribution (GD) and (2) they either increase or decrease the risks of allograft rejection. To this end, we compared the IL-18 genotypes to estimate odds ratios [ORs] and 95% confidence intervals using standard genetic models (homozygous, recessive, dominant and codominant). Subgrouping was ethnicity-based. Heterogeneous (random-effects) associations were subjected to outlier treatment which split the outcomes as pre- (PRO) and post- (PSO) outlier. Stability and robustness of the outcomes were analyzed by Bonferroni-correction and sensitivity treatment, respectively. Results: Our results revealed two core outcomes based on significance (P a < 0.05): (1) genotype frequency was hp than hc (OR 1.34, P a = 0.0007) in the codominant model (PSO) based on stability and robustness and (2) protection from allograft rejection (OR 0.74, P a = 0.04) in the dominant model (PRO) based on homogeneity. Subgroup analysis showed that Caucasian and Asian outcomes validated the GD and


Conclusions:
The IL-18 SNPs showed associations (hp) with KT up to 1.3-fold and protected KT recipients from allograft rejection (26%). Subgroup outcomes delineated the Asian and Caucasian effects. Enabled by outlier treatment, these findings were supported by non-heterogeneity. More studies should confirm or counter our findings.

Background
The end-stage of renal failure resulting from kidney disease [1] points to kidney transplantation (KT) as the optimal therapeutic choice [1,2]. The transplanted material (allograft) in the recipient is successful only if it is not rejected [3].
Unrejected allografts are expected to perform the functions as normal kidneys.
Normal post-KT graft outcomes depend on immunology where variation in immune responses of the recipient is genetically influenced [4]. This variation may help individualize immunosuppressive regimens by identifying alleles that could increase risk or confer protection for immune-mediated complications [5]. Cytokine proteins modulate and mediate the immune response [6] and the gene that encode these proteins influence transcription, yielding differences in cytokine production [7]. Of

Selection of studies
Three databases (PubMed, Google Scholar and Science Direct) were searched for genetic association articles as of September 24, 2019. Search terms included: "interleukin", "IL-18 ", "cytokine ", "polymorphisms", "allograft" and "renal transplantation". Where duplicate articles were encountered, the later dated one was selected. Inclusion criteria were (1) studies that associated IL-18 SNPs with KT outcomes; (2) IL-18 genotype frequencies that compare KT patients and healthy controls, NR and NRJ. (3) genotype frequency data that allowed calculation of the odds ratios (ORs) and 95% confidence intervals (CIs). Exclusion criteria were studies that (1) did not examine renal allografts or KT outcomes; (2) were reviews; (3) were not about the IL-18 SNPs and (4) had unusable genotype or allele frequencies.

Linkage disequilibrium and data extraction
The included articles examined two IL-18 SNPs, -137G/C (rs187238) and -607A/C (rs1946518), each presented with genotype data (Table 1 and S1 Table). Observed phenotypic associations have been attributed to the proximity of two SNPs [11,12].
NCI LDLINK (https://ldlink.nci.nih.gov/) results shows that the two SNPS are in linkage equilibrium (LD) based on both European (CEU) and Han (CHB) genotypes [13]. LD is defined as the correlation between alleles located near each other [14] which is measured in terms of D′ with a value of 1 indicating complete LD [15]. Therefore, IL-18 SNPs with D′ values of 1.00 in this study were reported to be in LD (S1 Table) and combined in the analysis (S2/S3 Table). Given these conditions, rs187238 and rs1946518 SNPs in IL-18 were combined. This combination allowed analysis by GD and allograft as well as subgrouping by ethnicity (Tables 2 and 3).
Two investigators (TE and NP) independently extracted data and arrived at a consensus. The following information was obtained from each publication: first author's name, year of the study, country of origin, ethnicity, age of the subjects, IL-18 SNPs (rs number) and the Clark-Baudouin (CB) score (Table 1). Sample sizes as well as genotype data between the RJ and NRJ were also extracted along with calculated outcome of the minor allele frequency (maf) (S2 and S3 Tables).

Power calculations and HWE assessment
Using the G*Power program [16], we evaluated statistical power as its adequacy bolsters the level of associative evidence. Assuming an OR of 1.5 at a genotypic risk level of α = 0.05 (two-sided), power was considered adequate at ≥ 80%. The Hardy-Weinberg Equilibrium (HWE) was assessed using the application in https://ihg.gsf.de/cgi-bin/hw/hwa1.pl.

Methodological quality of the studies
We used the CB scale to evaluate methodological quality of the included studies [17]. The CB criteria include P-values, statistical power, correction for multiplicity, comparative sample sizes between cases and controls, genotyping methods and the HWE. In this scale, low, moderate and high have scores of < 5, 5-6 and ³ 7, respectively.

Meta-analysis
We estimated odds ratios (ORs) and 95 % confidence intervals (CIs) of association using two overall approaches: (i) genotype distribution (GD) between cases and healthy controls and (ii) allograft wherein RJ were compared with NRJ. Calculated pooled ORs for GD were either higher in patients (hp) or higher in controls (hc); in allograft, they were either increased (in) or decreased (de), indicating risk for rejection. Standard genetic modeling was used, wherein we compared the following, and codominant (Co: wt versus var) effects. Heterogeneity between studies was estimated with the c 2 -based Q test [18], with threshold of significance set at P b < 0.10. Heterogeneity was also quantified with the I 2 statistic which measures variability between studies [19]. I 2 values of > 50% indicate more variability than those £ 50% with 0% indicating zero heterogeneity. Evidence of functional similarities in population features of the studies warranted using the fixed-effects model [20], otherwise the random-effects model [21] was used. Sources of heterogeneity were detected with the Galbraith plot [22] followed by re-analysis  Figure 1 outlines the study selection process in a PRISMA-sanctioned flowchart (Preferred Reporting Items for Systematic Reviews and Meta-Analyses). Initial search resulted in 39 citations, followed by a series of omissions that eventually yielded four articles for inclusion [23][24][25][26].  (Tables   2 and 3). Table 1 shows that the methodological quality of the component studies was moderate based on mean and SD (6.37 ± 1.24) of the normally distributed CB scores (SW test: P = 0.33). This meta-analysis followed the PRISMA guidelines (S5 Table). Table 2 delineates the overall pooled ORs by direction of effect, where GDs were hp (OR > 1.00) and decreased risk in the allograft analysis (OR < 1.00). The results

Meta-analysis outcomes
show five statistically significant ORs (P a < 0.05), three (all PSO) and two outcomes in GD and allograft, respectively. Of the five significant overall pooled ORs, only one survived the Bonferroni correction. This Bonferroni-surviving pooled OR became our core finding in GD, the PSO-derived Co model outcome indicating hp (OR 1.34, 95% CI 1.13-1.58, P a = 0.0007). This finding was validated in Caucasian subgroup (OR 1.32, 95% CI 1.04 to 1.66, P a = 0.02). Based on homogeneity (I 2 = 0%) and initial fixed-effects features, the other core outcome was found in the Do model of the allograft analysis (OR 0.74, 95% CI 0.55-0.98, P a = 0.04). This finding was validated sources of heterogeneity (outliers), located above the +2 confidence limit (Fig 3). In  Table 2. Table 3 shows the outcomes from sensitivity analysis where two of the nine (22%) comparisons were robust, one of which was the core finding (GD overall PSO).

Summary of associations
The main findings of this study point to a dichotomized pattern of significant effects where all GD ORs favored patients (hp) over that of controls (up to 1.6-fold). In contrast, all outcomes in the allograft analysis indicated reduced risks of rejection (as much as 30%). Subgroup validation of these effects indicate consistency.
Bonferroni correction and sensitivity analysis indicate stability and robustness, respectively, of the core outcomes.
In this meta-analysis, subgroup and outlier treatments have unraveled significant and non-heterogeneous and homogeneous associations that were not present in the component single-study outcomes. Conflicting outcomes between primary studies may be attributed to their lack of power and small sample sizes. Underpowered outcomes appear to be common in candidate gene studies [27] and are prone to the risk of Type 1 error.

IL-18 SNPs in allografts
The crucial role of IL-18 in kidney physiology lies in its involvement in the filtration,

Strengths and limitations
Interpreting our findings should consider its limitations and strengths. Limitations include: (i) All the component studies were underpowered and (ii) most significant outcomes were non-robust. On the other hand, the strengths comprise of the following: (i) the combined sample sizes translated to adequate statistical power (80%); (ii) none of the articles had control frequencies that deviated from HWE.
Confining the analysis to studies in HWE did not materially alter the pooled ORs; in fact, HWE-analysis validated the overall pooled effects. Given this outcome, the risk of genotyping errors appears to be a minor issue which minimizes methodological weakness in our study. Outlier treatment was key to generating significance and reducing heterogeneity. This demonstrates the utility of this meta-analysis tool in elevating the level of evidence for associations.

Conclusions
We have shown that associations of the IL-18 SNPs with KT outcomes were genetic     Tables   Table 1 Characteristics of the included studies in IL-18      (Table 3); û significant outcomes that did not survive the Bonferroni correction; significant outcome that did P    (Table   3); û significant outcomes that did not survive the Bonferroni correction  Forest plot outcome of outlier treatment in the allograft analysis of the Co model. Diamond d