Study characteristics
After exclusion of studies that not met the inclusion criteria, finally, 16 studies, published from 2001 to 2015, included 1,766 samples: 1,178 samples from NPC samples, and 588 samples from non-cancerous samples, were enrolled in the meta-analysis. The characteristics of included studies of RASSF1A methylation and risk of NPC were summarized in Table 1.
The number of NPC samples included studies ranged from 6 to 220 (mean: 42.07 ± 38.90), and number of control samples included studies ranged from 2 to 50 (mean: 21.78 ± 17.38). The patients’ ethnicity in 16 studies, comprised of 14 studies from Asian countries: China, Hong Kong, Singapore, Taiwan (counting for 87.50%), and 2 studies from African country: Tunisia (counting for 12.5%). Regarding to the test-method of RASSF1A methylation, twelve studies used MSP (counting for 75.00%), one study used COBRA (counting for 6.25%), one study used RT-qPCR (counting for 6.25%), one study used MMSP (counting for 6.25%), one study used MS-HRM counting for 6.25%).
Meta-analysis: The frequency of RASSF1A promoter methylation, and the association between RASSF1A gene methylation and NPC
The differences of RASSF1A gene methylation between the case-group and control-group in 16 studies, including 1,178 samples from NPC samples, and 588 samples from non-cancerous samples, were accessed. Concerning to the hetero heterogeneity between studies, as there was heterogeneity across studies in the case-group (Q = 443.5, p < 0.0001, I2 = 93.91%, 95%CI for I2 = 92.23-95.23), and no heterogeneity across studies in the control-group (Q = 9.67, p = 1.00, I2 = 0.00%, 95%CI for I2 = 0.00-0.00), thus, the random-effects model and fix-effects model were applied to calculate the frequency of RASSF1A gene in case-group and control-group, respectively (Fig. 1, Fig. 2). As the results, the frequency of RASSF1A gene methylation in case-group and control-group were 55.98% (ranged from 0.00% to 100.00%; 95% CI = 42.28-67.35) and 1.70% (ranged from 0.00 to 4.88%; 95% CI = 0.83-3.06), respectively. The meta-analysis result also indicated that the frequency of RASSF1A gene methylation in case-group was significant higher than control group.
The methylation of RASSF1A gene was significantly associated with an increased NPC risk with a pooled of Odds ratio (OR) of 33.40 (95%CI = 20.22-55.16) based on the fix-effect model (Q = 29.56, p = 0.29, I2 = 12.05%, 95%CI = 0.00-44.45) (Fig. 3). The funnel plot of pooled analysis, which was small asymmetric, indicated that there was still small bias among included studies, therefore publication bias cannot be completely excluded as a factor of influence on the present meta-analysis. In additionally, the sensitivity analysis found that study by Chang et al. (2003) (OR = 1.54, 95%CI = 0.09-24.10), Chang et al. (2003) (OR = 0.46, 95%CI = 0.02-11.79) and Tian et al. (2013) (OR = 4.14, 95%CI = 0.80-21.29) were the relative poor-quality studies. Therefore, those studies were omitted to evaluate again the association between the methylation of RASSF1A gene and NPC risk through OR. As the result, the between heterogeneity decreased to I2 = 0.00% (p = 0.99), indicated that no heterogeneity among enrolled studies was observed. Additionally, the association between RASSF1A methylation and risk of NPC increased, which was indicated by the increased OR of 51.43 (95%OR = 28.12-94.08) in fix-effects model (Q = 10.63, p = 0.99, I2 = 0.00, 95% CI = 0.00-0.00) (Fig. 5) (compared to previously calculated OR of 33.40). Moreover, the funnel plot of pooled analysis, which was quite symmetric, indicated that there was no significant bias among the included studies (Fig. 6).
The subgroup analysis by ethnicity, test-method, as well as source of sample were performed. The results of subgroup analysis showed that the heterogeneity totally disappeared in all subgroup (Asian vs non-Asian countries, MPS vs other test-methods, tissue vs other source of samples, I2 = 0.00, p > 0.05) (Table 2). With respect to the subgroups categorized by ethnicity, the significantly association between RASSF1A methylation and risk of NPC was observed among the Asian country and non-Asian countries in the fix-effect model (Asian countries: OR = 47.94, 95%CI = 25.61-89.75; non-Asian countries: OR = 135.96, 95%CI = 16.03-1,153.41). The subgroup analysis by source of NPC samples, there was a strong association between RASSF1A methylation and NPC among the NPC biopsy tissue group and non-biopsy group, including nasopharyngeal swab, mouth and rinsing fluid, plasma, blood in the fix-effect model (Biopsy tissue: OR = 55.29, 95%CI = 27.25-112.208; Non-biopsy: OR = 45.76, 95%CI = 15.56-134.60). Additionally, significant association between RASSF1A methylation and NPC risk among test-method subgroup was found in fix-effects model (MSP: OR = 40.86, 95%CI = 20.41-81.80; Other methods: OR = 76.24, 95%CI = 23.24-250.11).
Sensitivity analysis and publication bias
After removing the relative poor-quality study, the quite symmetric funnel plot, suggested there was no significant bias among included studies, was observed (Fig. 6). Additionally, the sensitivity analysis was performed to evaluate the stability and reliability of the conclusions according to the leave-one-out method by excluding one study. As the results, the pooled OR was ranged from 41.95 (95%CI = 22.13-79.49) to 55.06 (95%CI = 29.63-102.32) under the fix-effects model within the I2 = 0.00 (p > 0.05) (Table 3), which confirmed the stability and reliability of the results. Therefore, the results, and conclusion of present meta-analysis, which was to evaluate the association between methylation of RASSF1A and NPC risk, were stable and reliable.