- Open Access
DNA Methylation profiles as predictors of recurrence in non muscle invasive bladder cancer: an MS-MLPA approach
Journal of Experimental & Clinical Cancer Researchvolume 32, Article number: 94 (2013)
Although non muscle invasive bladder cancer (NMIBC) generally has a good long-term prognosis, up to 80% of patients will nevertheless experience local recurrence after the primary tumor resection. The search for markers capable of accurately identifying patients at high risk of recurrence is ongoing. We retrospectively evaluated the methylation status of a panel of 24 tumor suppressor genes (TIMP3, APC, CDKN2A, MLH1, ATM, RARB, CDKN2B, HIC1, CHFR, BRCA1, CASP8, CDKN1B, PTEN, BRCA2, CD44, RASSF1, DAPK1, FHIT, VHL, ESR1, TP73, IGSF4, GSTP1 and CDH13) in primary lesions to obtain information about their role in predicting local recurrence in NMIBC.
Formaldehyde-fixed paraffin-embedded (FFPE) samples from 74 patients operated on for bladder cancer were analyzed by methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA): 36 patients had relapsed and 38 were disease-free at the 5-year follow up. Methylation status was considered as a dichotomous variable and genes showing methylation ≥20% were defined as “positive”.
Methylation frequencies were higher in non recurring than recurring tumors. A statistically significant difference was observed for HIC1 (P = 0.03), GSTP1 (P = 0.02) and RASSF1 (P = 0.03). The combination of the three genes showed 78% sensitivity and 66% specificity in identifying recurrent patients, with an overall accuracy of 72%.
Our preliminary data suggest a potential role of HIC1, GSTP1 and RASSF1 in predicting local recurrence in NMIBC. Such information could help clinicians to identify patients at high risk of recurrence who require close monitoring during follow up.
Although superficial bladder cancer generally has a good long-term prognosis, up to 80% of patients will have local recurrence within 5 years of the primary tumor resection . After transurethral resection of bladder cancer (TURB), standard follow up involves numerous cystoscopies with consequently high healthcare costs and low patient compliance. Multiplicity, tumor size and prior relapse rate are the only recurrence-related parameters currently available for monitoring patients with bladder cancer , but such information would not seem to be accurate enough to ensure an adequate follow-up of individuals with stage Ta-T1 non muscle invasive bladder cancer (NMIBC). It would thus be extremely useful for clinicians to have new biological markers that can predict recurrence more accurately.
The role of epigenetic alterations in the carcinogenesis of solid tumors has been intensively investigated over the last ten years [2, 3]. DNA methylation at CpG rich regions often occurs at tumor suppressor gene promoters, frequently producing a reduction in the expression of target genes. An increasing number of papers are being published on the role of gene methylation and its potential clinical application in human tumors . Methylation seems to be an early event in the development of a number of solid tumors including bladder cancer [5, 6] and can thus be regarded as an early sign of cancer before the disease becomes muscle-invasive. Methylated tumor suppressor genes such as APC, RARB2, BRCA1 have recently been indicated as valid diagnostic markers for NMIBC [7–10]. A number of papers have also focused on the role of methylation as a prognostic marker, but it is not clear which methylated genes can accurately predict recurrence. Some studies have hypothesized hypermethylation of tumor suppressor genes, such as TIMP3, as a good prognostic marker [11, 12], while others have indicated hypermethylated E-cadherin, p16, p14, RASSF1, DAPK, APC, alone or in different combinations, as potential markers of early recurrence and poor survival [13–15].
In the present study we evaluated the methylation status of a panel of 24 genes (TIMP3, APC, CDKN2A, MLH1, ATM, RARB, CDKN2B, HIC1, CHFR, BRCA1, CASP8, CDKN1B, PTEN, BRCA2, CD44, RASSF1, DAPK1, FHIT, VHL, ESR1, TP73, IGSF4, GSTP1 and CDH13) in superficial bladder cancer to determine their ability to predict recurrence. Although methylation of some of these genes has already been investigated in bladder cancer [11–15], its relevance as an indicator of recurrence has yet to be confirmed. We used the relatively new methodology of methylation specific multiplex ligation dependent probe amplification (MS-MLPA) to evaluate epigenetic gene profiles. This approach permits methylation analysis of multiple targets in a single experiment [16, 17] and has been successfully used to evaluate the diagnostic or prognostic relevance of different markers in several tumor types such as lung , rectal , breast  and recently, bladder cancers [7, 8].
Case series (retrospective cohort study)
Tissue samples from 74 patients (65 males, 9 females) submitted to transurethral resection of primary bladder cancer at the Department of Urology of Morgagni-Pierantoni Hospital in Forlì between 1997 and 2006 were used for the study. All samples were retrieved from the archives of the Pathology Unit of the same hospital. Median age of patients was 73 years (range 39–92): 31 were <70 years and 43 ≥70 years. On the basis of 2004 World Health Organization criteria, final diagnosis was low grade non muscle invasive bladder cancer (NMIBC) in 55 patients and high grade NMIBC in 19 patients. At a median follow up of 5 years 38 patients were still disease-free and 36 had experienced one or more episodes of local recurrence. In this retrospective study, the two subgroups (disease free or relapsed) of patients were equally distributed for sex, age, grade and stage (Table 1).
All patients gave written informed consent for biological samples to be used for research purposes. The study protocol was reviewed and approved by the ‘Area Vasta’ Istituto Scientifico Romagnolo per lo Studio e la Cura dei Tumori (IRST) Ethics Committee.
Macrodissection and DNA isolation
Five 5-μm-thick sections were obtained from each paraffin-embedded block. Macrodissection was performed on hematoxylin-eosin stained sections and only cancer tissue was used for DNA isolation. Genomic DNA was purified using QIAmp DNA FFPE Tissue (Qiagen, Milan), according to the manufacturer’s instructions.
DNA was also isolated from a human bladder cancer cell line (HT1376) using Qiamp DNA minikit (Qiagen, Milan, Italy), according to the manufacturer’s instructions.
Methylation specific multiple ligation probe amplification (MS-MPLA)
MS-MLPA was performed using at least 50 ng of genomic DNA dissolved in 1XTE buffer (Promega, Madison, WI, USA). DNA isolated from HT 1376 cell line was used as internal control for MS MLPA analysis (Figure 1). The methylation status of 24 tumor suppressor gene promoters was analyzed using the ME001C1 kit (MRC-Holland, Amsterdam, The Netherlands) (Table 2). Two different probes that recognize two different sites of the promoter region were used for genes RASSF1 and MLH. We excluded CDKN2B gene from the analysis because its probe is sensitive to improper Hha1 digestion in FFPE samples. In brief, DNA was denatured (10 min at 98°C) and cooled at 25°C, after which the probe mix was added to the samples and hybridization was performed by incubation at 60°C for 16–18 h. The reaction was divided equally in two vials, one for ligation and the other for ligation-digestion reaction for each tumor. We added a mix composed of Ligase-65 buffer, Ligase-65 enzyme and water to the first vial and a mix of Ligase-65 Buffer, Ligase 65 enzyme, Hha1 enzyme (Promega, UK) and water to the second. The samples were then incubated at 49°C for 30 min. At the end of the ligation and ligation-digestion reactions, samples were amplified by adding a mix of PCR buffer, dNTPs and Taq polymerase. The PCR reaction was performed under the following conditions: 37 cycles at 95°C for 30 sec, 60°C for 30 sec and 72°C for 60 sec. The final incubation was performed at 73°C for 20 min.
Amplification products were analyzed by ABI-3130 genetic Analyzer (Applied Biosystem, UK). Universally methylated and unmethylated genomic DNA was used as positive or negative control, respectively.
Electropherograms obtained were analyzed using Gene Mapper software (Applied Biosystem, UK) and the peak areas of each probe were exported to a home-made excel spreadsheet. In accordance with the manufacturer’s instructions, we carried out “intrasample data normalization” by dividing the signal of each probe by the signal of every reference probe in the sample, thus creating as many ratios per probe as there were reference probes. We then calculated the median value of all probe ratios per probe, obtaining the normalization constant (NC). Finally, the methylation status of each probe was calculated by dividing the NC of a probe in the digested sample by the NC of the same probe in the undigested sample, and by multiplying this ratio by 100 to have a percentage value, as follows:
MS-MLPA technique reproducibility was assessed by performing three independent methylation profile analyses on a bladder cell line (HT1376). The methylation level for each gene was found to be the same in each experiment.
We considered the promoters showing a ratio ≥0.20 as methylated, while those with a ratio <0.20 were regarded as unmethylated. The cut-off was chosen on the basis of experiments performed on the bladder cancer cell line (HT1376) and on data from the literature [21, 22]. We have also performed the analysis on some samples from healthy tissues, to confirm that the background noise was inferior to 0.20 cut-off, such excluding false positive results due to experimental procedure.
Fisher’s exact test was used to compare the frequency of promoter methylation in the two subgroups: recurrent tumors versus non recurrent tumors. Methylation status was considered as a dichotomic variable and genes showing methylation ≥ 20% were classified as positive. A difference was considered significant if it showed a two-tailed P value ≤0.05. The genes showing a significant p value in Fisher’s exact test were used to analyze the methylator phenotype. Study endpoints were sensitivity (the proportion of recurrent cancer patients who were correctly identified by the test or procedures) and specificity (the proportion of non recurrent cancer patients who were correctly identified), with their 95% confidence intervals (CIs). We also evaluated overall accuracy, defined as the proportion of the total number of patients correctly identified by the test.
The student’s T test was used to assess the methylation index (MI), which was considered as a continuous variable. Logistic regression analysis was performed using the Epicalc of R to evaluate the performance of a panel of gene promoters (HIC1, RASSF1 and GSTP1) in discriminating between recurrent and non recurrent patients. We created logistic regression models with methylation levels of the three gene promoters (HIC1, RASSF1 and GSTP1). Probabilities were calculated as follows: P = exp ((Σ(bixi) + c)/(1 + Σ(bixi) + c), where p is the probability of each case, i = 1 to n; b is the regression coefficient of a given gene, x is the log2-transformed methylation level and c is a constant generated by the model. The ROCR package was used to obtain the ROC curves of the models and area under the curve (AUC) values. Recurrence-free survival was analyzed with the Log-rank test using SAS 9.3 software. All the molecular analyses were performed in a blind manner.
MS-MLPA analysis was feasible in all samples. The methylation frequency in the overall series varied widely (1% to 50%) for the different genes (Table 3). A separate analysis as a function of recurrence showed lower gene methylation in recurring than non recurring tumors, with the exception of CDKN1B, FHIT and IGSF4 genes. However, a significant difference between recurrent and non recurrent tumors was only observed for GSTP1, HIC1 and RASSF1 locus 2 (Table 3), with lower methylation in relapsed than non relapsed patients (Figure 2). The methylation index (MI), evaluated as the number of methylated genes relative to the total number of analyzed genes, showed values from 0 to 0.68 in the overall series of 23 genes and a significantly lower median value in non recurrent (0.08) than recurrent (0.12) (P = 0.011) patients (Table 4). To reduce the complexity of the methodological approach, further analysis was limited to a series of 10 genes (GSTP1, HIC1, RASSF1-locus2, CD44,DAPK, RASSF1-locus1, TP73, BRCA1, ESR1, TIMP3) that proved significant or showed a trend towards significance (P values varying from 0.02 to 0.31). Again, a higher median MI was seen in patients who relapsed compared to those who did not (0 versus 0.2; P = 0.0007) (Table 4).
We constructed a prognostic algorithm with the 3 significant genes (GSTP1, HIC1 and RASSF1) considering two phenotypes: the “methylated phenotype” (MP) (samples with at least one of the three genes methylated), and the “unmethylated phenotype” (samples with none of the three genes methylated). Of the 33 patients with methylated phenotype, 25 (76%) were still disease-free and 8 (24%) had had at least one intravescical recurrence at a median follow up of 5 years (Figure 3). Conversely, of the 41 patients with unmethylated phenotype, 28 (68%) had relapsed within 5 years of surgery and 13 (32%) had remained disease-free. The three-gene panel showed 78% sensitivity in identifying recurrent tumors and 66% specificity, with an overall accuracy of 72%.
We also performed ROC curve analysis for the three significant genes, singly or in combination, considered as continuous variables. Resultant AUCs were 0.5917 for HIC1, 0.6725 for RASSF1 and 0.5409 for GSTP1, the best AUC (0.6959) reached for the combination of the three genes (Figure 4).
Recurrence-free survival analysis of patients with methylated or unmethylated tumors highlighted a significantly higher recurrence-free survival (P = 0.0019) for those whose tumors showed the methylated phenotype (Figure 5).
The recurrence free survival analysis performed considering only the recurrent patients, showed that patients with unmethylated tumors had a lower median recurrent free survival time (14.5 months), with the respect to patients with methylated ones (18 months). However, the two subgroups are not equal distributed to give a statistical significant result (P = 0.9392, data not shown).
Multivariable analysis considering clinical and biological parameters (patient age and sex; tumor grade, stage and size; tumor multiplicity, methylated phenotype) showed that only age and methylated phenotype were independent predictors of recurrence. Specifically, patients under 70 years of age showed a higher probability of relapsing than older ones (P = 0.028) and their methylation phenotype was significantly predictive of recurrence (P < 0.0001).
The present study focused on evaluating the methylation status of tumor suppressor genes and on verifying its role in predicting recurrence of non muscle invasive bladder cancer (NMIBC). The MS-MLPA technique has the advantage of requiring only a small quantity of DNA, is capable of rapidly determining the methylation status of numerous genes in the same experiment, and has also been shown to work well in formalin-fixed paraffin-embedded samples. However, an important limitation of our study was the lack of a sufficient quantity of cancer tissue to confirm the methylation results using a second technique such as methylation specific PCR (MS PCR) or gene expression analyses.
In agreement with results from other studies , we found a positive correlation between gene methylation and lack of recurrence, highlighting that putative tumor suppressor genes do not always act as tumor suppressors but may actually have different biological functions.
Statistical analysis revealed 3 genes (HIC1, GSTP1, and RASSF1) capable of significantly predicting tumor recurrence. Their methylation was significantly indicative of a lack of recurrence at the 5-year follow up. The combined analysis of the three genes showed 72% accuracy in predicting recurrence or non recurrence.
HIC1 is a new candidate tumor suppressor gene , but the relevance of its methylation in bladder cancer prognosis is still unknown. Although GSTP1 methylation is a well known event in the carcinogenesis of prostate cancer, its role in bladder carcinoma has yet to be defined. A recent study by Pljesa-Ercegovac and coworkers  revealed that high GSTP1 expression is associated with an altered apoptotic pathway and bladder cancer progression. As methylation reduces gene expression, our data are in agreement with those of Pljesa-Ercegovac, the absence of GSTP1 methylation observed in our study supporting the hypothesis of more aggressive behavior of bladder tumors and consequently of a higher relapse rate.
Although the role of RASSF1 in bladder cancer development is still unclear, Ha and coworkers reported that its methylation would seem to play a part in predicting recurrence in low grade and stage bladder tumors . Surprisingly, we observed lower methylation levels of RASSF1 in recurrent tumors than in non recurrent ones, the discordance possibly due to different techniques used.
The MS-MLPA approach only permitted us to analyze one CpG site per probe, whereas several CpG sites may have been evaluated by Ha using the MS PCR technique . For these reasons, we believe that further evaluation is needed to clarify the role of RASSF1 in bladder cancer, especially with regard to the correlation between its methylation status and protein expression.We also observed fairly low methylation frequencies for all the loci analyzed compared to those reported in other papers . Such disagreement could, again, be due to the different analytical techniques adopted and/or to the different case series analyzed. Methylation cannot be the only mechanism of recurrence of NMIBC because the behavior of bladder tumors is fairly heterogeneous, as shown by Serizawa and coworkers  who observed an inverse correlation between FGFR mutations and hypermethylation events. In their study of the mechanisms of NMIBC recurrence, Bryan and coworkers , identified four reasons for relapse: incomplete resection, tumor cell re-implantation, growth of microscopic tumors and new tumor formation. These mechanisms differ greatly from each other and the identification of a single marker that is common to all four mechanisms appears improbable. It is more likely that a molecular marker characterizes tumor recurrence as a result of the third or fourth mechanisms, which may involve molecular alterations. This might explain why accuracy in our study only reached 72%.
Our preliminary findings pave the way for in depth evaluation of the methylation levels of HIC1, GSTP1, and RASSF1 genes in larger case series to improve the clinical surveillance of patients with superficial bladder cancer.
Written informed consent was obtained from the patient for the publication of this report and any accompanying images.
Van Rhijn BW, Burger M, Lotan Y, Solsona E, Stief CG, Sylvester RJ, et al: Recurrence and progression of disease in non-muscle-invasive bladder cancer: from epidemiology to treatment strategy. Eur Urol. 2009, 56 (3): 430-442. 10.1016/j.eururo.2009.06.028.
Sharma S, Kelly TK, Jones PA: Epigenetics in cancer. Carcinogenesis. 2010, 31 (1): 27-36. 10.1093/carcin/bgp220.
Tao W, Hongli L, Yeshan C, Wei L, Jing Y, Gang W: Methylation associated inactivation of RASSF1A and its synergistic effect with activated K-Ras in nasopharyngeal carcinoma. J Exp Clin Cancer Res. 2009, 28: 160-10.1186/1756-9966-28-160.
Jian Z, Yuyan W, Jianchun D, Hua B, Zhijie W, Lai W: DNA Methylation status of Wnt antagonist SFRP5 can predict the response to the EGFR-tyrosine kinase inhibitor therapy in non-small cell lung cancer. J Exp Clin Cancer Res. 2012, 31: 80-10.1186/1756-9966-31-80.
Sánchez-Carbayo M: Hypermethylation in bladder cancer: biological pathways and translational applications. Tumour Biol. 2012, 33 (22): 347-361.
Kim WJ, Kim YJ: Epigenetics of bladder cancer. Methods Mol Biol. 2012, 863: 111-118.
Cabello MJ, Grau L, Franco N, Orenes E, Alvarez M, Blanca A, et al: Multiplexed methylation profiles of tumor suppressor genes in bladder cancer. J Mol Diagn. 2011, 13 (1): 29-40. 10.1016/j.jmoldx.2010.11.008.
Zuiverloon TC, Beukers W, van der Keur KA, Munoz JR, Bangma CH, Lingsma HF, et al: A methylation assay for the detection of non-muscle-invasive bladder cancer (NMIBC) recurrences in voided urine. BJU Int. 2012, 109 (6): 941-948. 10.1111/j.1464-410X.2011.10428.x.
Eissa S, Swellam M, El-Khouly IM, Kassim SK, Shehata H, Mansour A, et al: Aberrant methylation of RARbeta2 and APC genes in voided urine as molecular markers for early detection of bilharzial and nonbilharzial bladder cancer. Cancer Epidemiol Biomarkers Prev. 2011, 20 (8): 1657-1664. 10.1158/1055-9965.EPI-11-0237.
Negraes PD, Favaro FP, Camargo JL, Oliveira ML, Goldberg J, Rainho CA, et al: DNA methylation patterns in bladder cancer and washing cell sediments: a perspective for tumor recurrence detection. BMC Cancer. 2008, 8: 238-10.1186/1471-2407-8-238.
Hoque MO, Begum S, Brait M, Jeronimo C, Zahurak M, Ostrow KL, Rosenbaum E: Tissue inhibitor of metalloproteinases 3 promoter methylation is an independent prognostic factor for bladder cancer. J Urol. 2008, 179 (2): 743-747. 10.1016/j.juro.2007.09.019.
Friedrich MG, Chandrasoma S, Siegmund KD, Weisenberger DJ, Cheng JC, Toma MI, et al: Prognostic relevance of methylation markers in patients with non-muscle invasive bladder carcinoma. Eur J Cancer. 2005, 41 (17): 2769-2778. 10.1016/j.ejca.2005.07.019.
Tada Y, Wada M, Taguchi K, Mochida Y, Kinugawa N, Tsuneyoshi M, et al: The association of death-associated protein kinase hypermethylation with early recurrence in superficial bladder cancers. Cancer Res. 2002, 62 (14): 4048-4053.
Lin HH, Ke HL, Wu WJ, Lee YH, Chang LL: Hypermethylation of E-cadherin, p16, p14, and RASSF1A genes in pathologically normal urothelium predict bladder recurrence of bladder cancer after transurethral resection. Urol Oncol. 2012, 30 (2): 177-181. 10.1016/j.urolonc.2010.01.002.
Yates DR, Rehman I, Abbod MF, Meuth M, Cross SS, Linkens DA, et al: Promoter hypermethylation identifies progression risk in bladder cancer. Clin Cancer Res. 2007, 13 (7): 2046-2053. 10.1158/1078-0432.CCR-06-2476.
Schouten JP, McElgunn CJ, Waaijer R, Zwijnenburg D, Diepvens F, Pals G: Relative quantification of 40 nucleic acid sequences by multiplex ligation-dependent probe amplification. Nucleic Acids Res. 2002, 30 (12): e57-10.1093/nar/gnf056.
Nygren AO, Ameziane N, Duarte HM, Vijzelaar RN, Waisfisz Q, Hess CJ, et al: Methylation-specific MLPA (MS-MLPA): simultaneous detection of CpG methylation and copy number changes of up to 40 sequences. Nucleic Acids Res. 2005, 33 (14): e128-10.1093/nar/gni127.
Castro M, Grau L, Puerta P, Gimenez L, Venditti J, Quadrelli S, et al: Multiplexed methylation profiles of tumor suppressor genes and clinical outcome in lung cancer. J Transl Med. 2010, 8: 86-10.1186/1479-5876-8-86.
Leong KJ, Wei W, Tannahill LA, Caldwell GM, Jones CE, Morton DG, et al: Methylation profiling of rectal cancer identifies novel markers of early-stage disease. Br J Surg. 2011, 98 (5): 724-734. 10.1002/bjs.7422.
Moelans CB, Verschuur-Maes AH, Van Diest PJ: Frequent promoter hypermethylation of BRCA2, CDH13, MSH6, PAX5, PAX6 and WT1 in ductal carcinoma in situ and invasive breast cancer. J Pathol. 2011, 225 (2): 222-231. 10.1002/path.2930.
Pavicic W, Perkiö E, Kaur S, Peltomäki P: Altered methylation at microRNA-associated CpG islands in hereditary and sporadic carcinomas: a methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA)-based approach. Mol Med. 2011, 17 (7–8): 726-735.
Joensuu EI, Abdel-Rahman WM, Ollikainen M, Ruosaari S, Knuutila S, Peltomäki P: Epigenetic signatures of familial cancer are characteristic of tumor type and family category. Cancer Res. 2008, 68 (12): 4597-4605. 10.1158/0008-5472.CAN-07-6645.
Fleuriel C, Touka M, Boulay G, Guérardel C, Rood BR, Leprince D: HIC1 (Hypermethylated in Cancer 1) epigenetic silencing in tumors. Int J Biochem Cell Biol. 2009, 41 (1): 26-33. 10.1016/j.biocel.2008.05.028.
Pljesa-Ercegovac M, Savic-Radojevic A, Dragicevic D, Mimic-Oka J, Matic M, Sasic T, Pekmezovic T: Enhanced GSTP1 expression in transitional cell carcinoma of urinary bladder is associated with altered apoptotic pathways. Urol Oncol. 2011, 29 (1): 70-77. 10.1016/j.urolonc.2008.10.019.
Ha YS, Jeong P, Kim JS, Kwon WA, Kim IY, Yun SJ: Tumorigenic and prognostic significance of RASSF1A expression in Low-grade (WHO grade 1 and grade 2) nonmuscle-invasive bladder cancer. Urology. 2012, 79 (6): e1-e6. 1411
Kim YK, Kim WJ: Epigenetic markers as promising prognosticators for bladder cancer. Int J Urol. 2009, 16 (1): 17-22. 10.1111/j.1442-2042.2008.02143.x.
Serizawa RR, Ralfkiaer U, Steven K, Lam GW, Schmiedel S, Schüz J, et al: Integrated genetic and epigenetic analysis of bladder cancer reveals an additive diagnostic value of FGFR3 mutations and hypermethylation events. Int J Cancer. 2011, 129 (1): 78-87. 10.1002/ijc.25651.
Bryan RT, Collins SI, Daykin MC, Zeegers MP, Cheng KK, Wallace DM, et al: Mechanisms of recurrence of Ta/T1 bladder cancer. Ann R Coll Surg Engl. 2010, 92 (6): 519-524. 10.1308/003588410X12664192076935.
The authors thank Grainne Thierney and Ursula Elbling for editing the manuscript.
The authors declare that they have no competing interests.
VC carried out the molecular genetic studies and drafted the manuscript; CM, DC, MT carried out the molecular genetic studies; RG, LS, FF participated in recruitment of patients and collection and assembly of data; CZ performed statistical analysis; RS helped to draft the manuscript and participated in the design of the study; DA and WZ participated in the design of the study and coordination. All authors read and approved the final manuscript.