High expression of RBP7 correlates with better prognosis in breast cancer
Previous studies had indicated that high expression of RBP7 was an independent biomarker of poor cancer specific survival in breast cancer, colon cancer and bladder cancer. Through analysis the mRNA level of RBP7 using TCGA database, as shown in Fig. 1A, the expression of RBP7 was downregulated in breast cancer, colon cancer and bladder cancer tissues compared with their normal samples in TCGA. What’s more, the expression level of RBP7 was not related with tumor stage in breast cancer (Fig. 1B). However, the Kaplan-Meier curve indicated that the breast cancer patients with high level of RBP7 mRNA trends to have a longer survival time especially in the Lumina A subtype (Fig. 2).
Methylation status of RBP7 gene in breast cancer tissues and its association with clinicopathological features
We then examined promoter methylation of RBP7 using MSP approach in a large cohort of breast cancer samples. As shown in Table 1, the methylation rate of RBP7 was 59.1% (110/186) in the cohort. However, no significant difference of RBP7 methylation rate was seen between WHO grade, LNM, radiotherapy, chemotherapy and either hormone receptors.
In order to further examine the relationship of RBP7 methylation status with clinicopathological characteristics in this cohort. Logistic regression analysis was then explored. The univariate analyses showed that lower RBP7 methylation rate was found to be significantly associated with patients with chemotherapy (p = 0.03) and patients with lymph node metastasis (p = 0.02). However, higher RBP7 methylation rate was associated with Her2 positive group (p = 0.004), cancer relapsed (p = 0.04), and cancer related death (p = 0.01) (Table 2). In order to assess the independent association of gene methylation with molecular status (including ER, PR and Her2), lymph node metastasis, radiotherapy, chemotherapy and cancer-related death, we conducted multiple multivariable logistic regressions. As also shown in Table 2, lower RBP7 methylation rate remained significantly associated with lymph node metastasis. However, higher RBP7 methylation rate remained associated with Her2 status and cancer related death.
Table 2
Methylation of RBP7 in breast cancer: univariate and multivariate models with clinicopathological characteristics
Characteristics | Univariate | | Multivariate | |
OR* (95% CI) | P | | OR* (95% CI) | P | |
ER | 1.58 (0.87–2.86) | 0.13 | | 6.16 (1.97–19.27) | 0.02 | |
PR | 0.84 (0.47–1.51) | 0.56 | | 0.35 (0.12–1.03) | 0.06 | |
Her2 | 3.49 (1.06–11.44) | 0.004* | | 5.42 (1.91–15.43) | 0.002 | |
LNM | 0.50 (0.28–0.91) | 0.02 | | 0.38 (0.18–0.81) | 0.01 | |
Radiotherapy | 1.11 (0.53–2.33) | 0.79 | | 2.76 (1.08–7.09) | 0.04 | |
Chemotherapy | 0.38 (0.16–0.89) | 0.03 | | 0.64 (0.24–1.72) | 0.37 | |
Relapse | 2.25 (1.02–4.97) | 0.04 | | n/a | n/a | |
Survival status | 2.90 (1.24–6.78) | 0.01 | | 4.71 (1.66–13.35) | 0.04* | |
*OR: odds ratio with 95% confidence interval; Survival status (alive vs. dead). |
Given that RBP7 methylation was found to be associated with the key molecular status including ER, PR and Her2 in our cohort. We then divided the cohort based on this molecular status and examine the relationship of methylated RBP7 with clinicopathological characteristics again using logistic regression analysis. In the ER genitive group, as shown in Table S1, the univariate and multivariable logistic regressions showed that higher RBP7 methylation rate was significantly associated with Her2 positive group and cancer related death. Besides, in the PR genitive group, as shown in Table S2, we also found that higher RBP7 methylation rate was significantly associated with Her2 positive group and cancer related death. What’s more, in the Her2 genitive group, as shown in Table S3, we found that higher RBP7 methylation rate was significantly associated with ER positive group and cancer related death both using univariate and multivariable logistic regressions. The above results indicated that patients with high level of promoter methylation of RBP7 might have a worse survival time especially in those patients with some negative hormone receptor.
RBP7 gene promoter methylation predicted worse prognosis in breast cancer
Whether methylation of RBP7 is indeed associated with a worse survival in breast cancers, as suggested by its association with clinicopathological characteristics of breast cancer patients, was subsequently investigated in unvariate and multivariable survival analysis. As shown in Table 3, univariate Cox regression showed RBP7 methylation (p = 0.001), ER negative status (p = 0.02), Her2 positive status (p = 0.001), patients without radiotherapy (p = 0.04), patients without chemotherapy (p = 0.003), and cancer relapse (p < 0.001) acted as a predictor of poor survival for breast cancer patients. Cox multivariate regression showed that RBP7 methylation (p < 0.001) and lymph node metastasis (p < 0.001) is also a predictor of poor survival as an independently variable with respect to the ER, PR, Her2. And radiotherapy act as a predictor of better survival factor (Table 4).
Table 3
Prognostic value of clinicopathological factors and methylation of RBP7 using univariate Cox regression analysis
Characteristics | HR* (95% CI) | P |
RBP7 methylation | 5.67 (1.98–16.20) | 0.001 |
ER | 0.44 (0.22–0.86) | 0.02 |
PR | 0.55 (0.28–1.08) | 0.08 |
Her2 | 3.04 (1.55–5.98) | 0.001 |
LNM | 2.62 (1.31–5.24) | 0.07 |
Radiotherapy | 0.23 (0.05–0.94) | 0.04 |
Chemotherapy | 0.35 (0.17–0.70) | 0.003 |
Relapse | 127.132 (30.31-533.25) | < 0.001 |
*HR: hazard ratio with 95% confidence interval (CI) |
Table 4
Prognostic value of clinicopathological factors and methylation of RBP7 using multivariate Cox regression analysis
Characteristics | HR* (95% CI) | P |
RBP7 methylation | 9.52 (3.05–29.72) | < 0.001 |
ER | 0.39 (0.14–1.12) | 0.08 |
PR | 1.22 (0.46–3.26) | 0.69 |
Her2 | 1.06 (0.42–2.72) | 0.91 |
LNM | 4.27 (2.05–8.91) | < 0.001 |
Radiotherapy | 0.12 (0.03–0.54) | 0.006 |
*HR: hazard ratio with 95% confidence interval (CI) |
In order to understand the relationship of RBP7 methylation with the key molecular status including ER, PR and Her2 in the survival status of these breast cancer patients. We then again divided the cohort based on these molecular status and explored unvariate and multivariable survival analysis using Cox regression. As show in Table S4, in the ER genitive group, the univariate and multivariable Cox regression analysis showed that RBP7 methylation and lymph node metastasis acted as a predictor of poor survival. Besides, as show in Tables 5 and 6, in the PR and Her2 genitive group, we also found that RBP7 methylation acted as a predictor of poor survival.
The Kaplan-Meier estimator of the survivorship function was used to evaluate the impact of aberrant gene methylation on the survival of breast cancer patients. The survival of cancer patients with and without gene methylation was compared using the log-rank test. As shown in Fig. 3, methylation of RBP7 was significantly associated with worse survival of breast cancer patients (median survival time: 71.4 vs 90 months; p = 0.019). What’s more, we also found that methylation of RBP7 was significantly associated with worse survival of breast cancer patients both in the Her2 negative and ER negative subtypes (Fig. 3). Since lymph node metastasis, cancer relapse, radiotherapy and chemotherapy play an important role of survival time. We then analyze the survival times after surgery under those conditions. As also shown in Fig. 3, methylation of RBP7 was significantly associated with worse survival both in the patients with lymph node metastasis and patients without radiotherapy subtype.