2.1. Protocol and Registration
This systematic review and meta-analysis study was prospectively registered with the National Institute for Health Research (PROSPERO: CRD42020216086). We followed the Reporting Checklist for Meta-analyses of Observational Studies (MOOSE) (23) and Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) (24) statements to report this study.
2.2. Eligibility Criteria
2.2.1. Type of Study
All observational studies (such as cross-sectional, case-control, and cohort), as well as interventional studies (randomized or non-randomized), were included in this study. In contrast, case-report, case-series, narrative reviews, animal studies, and letters were excluded.
2.2.2. Type of Participants
Kidney transplant recipients with and without acute rejection were enrolled in this study, which their pre- or post-transplant 25-hydroxyvitamin D levels were assessed. Recipients features are as follows: (i) biopsy-proven or clinically suspected acute rejection, (ii) received living or deceased donor kidney transplant, (iii) received their transplant for the first time or more, (iv) have experienced acute rejection up to 12 months after transplantation, and (v) with any age and gender. We excluded studies that did not separate the acute and chronic rejected patients and that measured 1,25-dihydroxyvitamin D instead of 25-hydroxyvitamin D.
2.2.3. Outcome
The outcome is the occurrence of acute rejection in kidney transplant recipients. Acute rejection is any type of cell-mediated or antibody-mediated rejection with any score or grade.
2.3. Search Strategy
The electronic databases searched by authors were PubMed, Scopus, Embase, ProQuest, Web of Science, and Google Scholar. Moreover, the gray literature (conferences/congress papers and theses) and reference lists of included studies were searched. The primary studies are from 1 January 1990 up to 31 October 2020. We used the keywords “kidney transplantation,” “acute rejection,” “25-hydroxyvitamin D,” and their synonym for the search syntax. The Medical Subject Headings (MeSH) and free-text method were used to create search syntax. The search syntax of the PubMed database is shown in Supplemental File 1.
2.4. Study Selection
The records searched in all databases were exported to the EndNote software, and duplicated records were removed. Authors performed the screening process using titles and abstracts. Then, two reviewers independently read the full text of screened papers and selected those that had the main data for conducting meta-analysis. Any disagreement was resolved by consensus.
2.5 Risk of Bias Assessment
The quality assessment of included studies was independently done by two reviewers using the Newcastle-Ottawa Quality Assessment Scale (NOS) for observational studies (25). Moreover, the quality of a trial study was assessed using the Cochrane checklist (26). The modified version of NOS checklists for observational studies are shown in Supplemental File 2.
2.6. Data Extraction
Two reviewers independently extracted data from included studies using a data extraction form as follows: first author name, publication year, maintenance therapy, induction therapy (IT), donor type (living vs. deceased), estimated glomerular filtration rate (eGFR), and age and gender of study subjects. Any disagreement was resolved by consensus, and, for incomplete data, we contacted the corresponding authors of the studies.
2.7. Data Synthesis and Analyses
The main data for conducting the meta-analysis are as follows (i.e., four numbers): the numbers of non-acute and acute rejection patients in the VitD-deficient (VitD < 20 ng/ml) group and the numbers of non-acute and acute rejection patients in the VitD-sufficient (VitD ≥ 20 ng/ml) group. We used the term “sufficient” for VitD ≥ 20 ng/ml patients in the whole manuscript and figures. According to the type of included studies and main data for performing the meta-analysis, we employed OR as an effect size to be reported in this study. According to the methodological heterogeneity between included studies, we selected a random-effects model as an appropriate combination model (27). The pooled effect size was plotted using a forest plot.
Statistical heterogeneity between included studies was evaluated using the I2 index (28). The severe statistical heterogeneity was considered I2 > 50%.
2.8. Publication Bias Assessment
We used three methods to assess publication bias as follows: the funnel plot (which is the visual assessment of publication bias), Begg’s (29) and Egger’s (30) tests (whish statistically assess publication bias), and trim and fill method (which assesses probable missing papers) (31).
2.9. Additional Analyses
2.9.1. Subgroup Analyses
We performed some subgroup analyses to find out potential sources of statistical heterogeneity observed between included studies. The other purpose of subgroup analyses was to know whether the pooled effect size was influenced by different variables in each subgroup analysis.
Subgroup analyses were performed as follows:
Quality of Studies
Studies with high quality (i.e., ≥ 6) were compared with those with low quality (i.e., < 6).
Induction Therapy
Studies with induction therapy (IT) were compared with those without IT.
Age
Studies with patients ≤ 47 years were compared with those with patients > 47 years in each group (VitD < 20 ng/ml and VitD ≥ 20 ng/ml) separately.
Gender
Studies with more male patients were compared with those with fewer male patients.
Donor Type
Studies with living donor more than 50% were compared with those with living donor less than 50% in each group (VitD < 20 ng/ml and VitD ≥ 20 ng/ml) separately.
2.9.2. Meta-Regression Analysis
Regarding the low number of studies that reported the eGFR of their patients, we decided to do meta-regression to evaluate whether eGFR affect the pooled effect size.
2.9.3. Sensitivity Analysis
We utilized the one-out-remove method to do sensitivity analysis. In this method, the pooled effect size was recalculated after removing the effect of each study. This method was used to know whether there was any difference between the results after removing a study. If yes, how much does the recalculated pooled effect size change?
This study was analyzed using STATA Version 13 (StataCorp, College Station, TX, USA).