Patients with low nicotinamide N-methyltransferase expression benefit significantly from bevacizumab treatment in ovarian cancer
BMC Cancer volume 21, Article number: 67 (2021)
The role of nicotinamide N-methyltransferase (NNMT) in ovarian cancer is still elusive. Our aim is to explore the expression of NNMT in ovarian cancer and to assess its association with patient prognosis and treatment response.
We first analyzed the differential expression of NNMT among fallopian tube epithelium, primary ovarian cancers, metastatic ovarian cancers, and recurrent ovarian cancers using Gene Expression Ominus (GEO) database (GSE10971, GSE30587, GSE44104 and TCGA datasets). Then, we assessed the association of NNMT expression with clinical and molecular parameters using CSIOVDB database and GSE28739 dataset. Next, we evaluate the association of NNMT expression with the prognosis of ovarian cancer patients in both GSE9891 dataset and TCGA dataset. Finally, GSE140082 dataset was used to explore the association of NNMT expression with bevacizumab response.
NNMT expression was significantly elevated in lymphovascular space invasion (LVSI)-positive ovarian cancers compared with that in LVSI-negative ovarian cancers (TCGA dataset, P < 0.05), Moreover, increased expression of NNMT was associated with increased tumor stage, grade, and mesenchymal molecular subtype (CSIOVDB database). Survival analysis indicated that increased expression of NNMT was associated with a reduced OS in both GSE9891 dataset (HR: 2.28, 95%CI: 1.51–3.43, Log-rank P < 0.001) and TCGA dataset (HR: 1.55, 95%CI: 1.02–2.36, Log-rank P = 0.039). Multivariate analysis further confirmed the negative impact of NNMT expression on OS in ovarian cancer patients in those two datasets. Furthermore, the NNMT-related nomogram showed that NNMT shared a larger contribution to OS, compared with debulking status. More interestingly, bevacizumab conferred significant improvements in OS for patients with low NNMT expression (HR: 0.56, 95%CI: 0.31–0.99, Log-rank P = 0.049). In contrast, patients with high NNMT expression didn’t benefit from bevacizumab treatment significantly (HR: 0.85, 95%CI: 0.48–1.49, Log-rank P = 0.561). NNMT expression was positively correlated with the expression of genes, LDHA and PGAM1, involved in Warburg effect.
In conclusion, NNMT expression is associated with the aggressive behavior of ovarian cancer, correlates with a poor prognosis, and is predictive of sensitivity to bevacizumab treatment.
Ovarian cancer, usually diagnosed at an advanced stage and characterized by metastatic bulky disease burden, has the highest mortality rate of all gynecologic cancer . After debulking surgery and combined chemotherapy, most of the patients relapse and die from drug resistance. This underscores the significant clinical need to decipher the molecular biology of ovarian cancer and to identify new biomarkers to predict prognosis, which may be translated into personalized treatment strategies and survival improvements.
Nicotinamide N-methyltransferase (NNMT) is an S-adenosylmethionine-dependent enzyme, which plays a critical role in the biotransformation and detoxification of many drugs and xenobiotic compounds . Overexpression of NNMT has been implicated in various cancers, including but not limited to colon, lung, hepatocellular, and bladder cancer [3,4,5,6]. In addition, increased expression of NNMT has been associated with tumor aggressiveness and demonstrated to facilitate the migration, invasion, viability, and proliferation of various cancer cells [3, 7,8,9,10,11,12]. Recently, NNMT has also been demonstrated to a master metabolic regulator of in several cancers and may be a therapeutically targeted [13, 14]. However, its role in ovarian cancer is still largely elusive [12, 13].
In present study, we analyzed the expression pattern of NNMT in ovarian cancer and evaluate its association with patient prognosis by taking advantage of the Gene Expression Omnibus (GEO) dataset and The Cancer Genome Altas (TCGA) dataset.
Database used in present study
The publicly available gene expression data used in our study is described in GEO database (GSE10971, GSE30587, GSE44104, GSE28739, GSE9891, and GSE140082). GSE10971 dataset was used to determine the differential expression of NNMT (probe ID:202237_at and 202238_s_at) between fallopian tube epithelium (FTE) (n = 24) and primary high grade serous tubal cancer/ovarian cancers (n = 13). GSE30587 dataset was used to determine the differential expression of NNMT (probe ID: 7943998) between paired primary ovarian cancer tissues(n = 9) and their omental metastasis counterparts(n = 9). GSE44104 dataset was used to determine the differential expression of NNMT (probe ID:202237_at and 202238_s_at) between ovarian cancer with recurrence (n = 20) and ovarian cancer with non-recurrence(n = 40). GSE28739 dataset was used to determine the differential expression of NNMT (probe ID: A_23_P127584) between drug sensitive ovarian cancer (n = 20) and drug resistant ovarian cancer (n = 30). If there are multiple probes all targeting NNMT, the average value of all probes is taken as an expression value of NNMT. Level 3 gene expression data with lymphovascular space invasion (LVSI) from TCGA dataset was used to explore the differential expression of NNMT between LVSI-positive ovarian cancers and LVSI-negative ovarian cancers.
GSE9891 dataset and TCGA dataset were used to develop and validated a NNMT-related nomograms to predict overall survival (OS) respectively. We only included patients with enough information (including age at initial diagnosis, histological grade, FIGO stage, debulking status, days to death, and survival status) in our final analysis. Totally, 242 patients in GSE9891 dataset and 462 patients in TCGA dataset were included. The expression value of NNMT expression together with the clinical information in GSE9891 dataset and TCGA dataset were extracted from the “curatedOvarianData” Bioconductor package (version 2.12 for R 3.0.3).
Kaplan-Meier analysis of OS for combination therapy of standard treatment and bevacizumab vs. standard treatment alone in ovarian cancer patients stratified by NNMT expression was performed in GSE140082 dataset.
Detailed methods of statistical analysis were described in our previous published paper . The expression of genes in GSE9891 and TCGA dataset were extracted from the curatedOvarianData R package. Since the present study was a single gene-centric study, we retained both the P value generated from student-t test and the adjusted P value generated from Limma R package for comparing gene expression between groups . When using limma R package, the Benjamini-Hochberg multiple comparison adjustment was applied and the adjusted P value < 0.05 was considered statistically significant. “Low” and “High” were classified according to the gene expression level. The optimal cutoff values for NNMT, lactate dehydrogenase A (LDHA), and phosphoglycerate mutase 1 (PGAM1) expression were determined by the maximally selected rank statistics using R package survminer (http://www.sthda.com/english/rpkgs/survminer/). Correlation analysis was performed by Pearson method. All the statistical analyses were performed using IBM SPSS Statistics and R. P values < 0.05 are considered statistically significant.
Differential expression of NNMT in fallopian tube epithelium (FTE), primary ovarian cancers (POCs)/primary tubal cancers (PTCs), metastatic ovarian cancers (MOCs), and recurrent ovarian cancers (ROCs)
The mean expression value of NNMT was elevated in POCs/PTCs compared with that in FTs (Fig. 1a, GSE10971), and was increased in MOCs and ROCs compared with POCs (Fig. 1b, GSE30587; Fig. 1c, GSE44104). But the differences were not statistically significant suggested by adjusted P value generated from Limma R package. However, NNMT expression was significantly elevated in LVSI-positive ovarian cancers compared with LVSI-negative ovarian cancers (Fig. 1d, TCGA dataset, P < 0.05).
The association of NNMT expression with clinical and molecular parameters
Next, we try to explore whether NNMT expression was correlated with clinical and molecular parameters using CSIOVDB database. We showed that increased expression of NNMT was associated with increased tumor stage (Fig. 2a), grade (Fig. 2b), and mesenchymal molecular subtype (Fig. 2c). The mean expression value of NNMT was elevated in drug-resistant ovarian cancers compared with drug-sensitive ovarian cancers (Fig. 2d), but the difference was not statistically significant suggested by adjusted P value generated from Limma R package.
The association of NNMT expression with OS
Next, we evaluate the association of NNMT expression with the prognosis of ovarian cancer patients. Interestingly, survival analysis indicated that increased expression of NNMT was associated with a reduced OS in both GSE9891 dataset (Fig. 3a, HR: 2.28, 95%CI: 1.51–3.43, Log-rank P < 0.001) and TCGA dataset (Fig. 3b, HR: 1.55, 95%CI: 1.02–2.36, Log-rank P = 0.039). Multivariate analysis further confirmed the negative impact of NNMT expression on OS in ovarian cancer patients in both GSE9891 dataset (Table 1) and TCGA dataset (Table 2).
To quantitatively predict the survival of ovarian cancer patients, we developed a nomogram using GSE9891 dataset (training dataset) based on age at initial diagnosis, debulking status, FIGO stage, and NNMT expression, and validated it using TCGA dataset (validation dataset). Among the included factors, age at initial diagnosis had the largest contribution to OS, followed by FIGO stage and NNMT expression (Fig. 4). The C-index for predicting OS was 0.70 (95%CI, 0.64–0.75) in GSE9891 dataset (training dataset) and 0.64 (95%CI, 0.60–0.68) in TCGA dataset (validation dataset), respectively. The calibration plots demonstrated optimal agreement between the predicted and observed probability (Fig. 5a-f).
The influence of NNMT expression on bevacizumab response
Then, we explore whether NNMT expression was predictive of response to bevacizumab usding GSE140082 dataset. Interestingly, it was found that bevacizumab conferred significant improvements in OS for patients with low NNMT expression (Fig. 6a, HR: 0.56, 95%CI: 0.31–0.99, Log-rank P = 0.049). In contrast, patients with high NNMT expression didn’t benefit from bevacizumab treatment significantly (Fig. 6b, HR: 0.85, 95%CI: 0.48–1.49, Log-rank P = 0.561).
The correlations of NNMT expression with the expression of genes involved in Warburg effect
Finally, we explore whether NNMT expression was correlated with the expression of genes (hexokinase 2 (HK2), LDHA, and PGAM1) involved in Warburg effect. Interestingly, our results suggested that NNMT expression was positively correlated with LDHA (Fig. 7a, r = 0.347, P < 0.001) and PGAM1 (Fig. 7b, r = 0.38, P < 0.001) expression, but not with HK2 expression (Fig. 7c, r = − 0.01, P = 0.905). The positive correlations of NNMT expression with the expression of LDHA (Fig. 7d, r = 0.15, P = 0.003) and PGAM1 (Fig. 7e, r = 0.30, P < 0.001) were confirmed in GSE140082 dataset. LDHA expression was not associated with bevacizumab response in ovarian cancer patients (Fig. 8a, Fig. 8b). However, bevacizumab conferred significant improvements in OS for patients with low PGAM1 expression (Fig. 8c, HR: 0.34, 95%CI: 0.12–0.94, Log-rank P = 0.037). In contrast, patients with high PGAM1 expression didn’t benefit from bevacizumab treatment significantly (Fig. 8d, HR: 0.80, 95%CI: 0.52–1.24, Log-rank P = 0.325).
The major findings of the present study can be summarized as follows. First, NNMT expression was significantly increased in LVSI-positive ovarian cancer tissues compared with LVSI-negative ovarian cancer tissues. Second, increased expression of NNMT was associated with increased tumor stage, grade, and mesenchymal molecular subtype. Third, higher NNMT expression correlated with unfavorable outcomes in ovarian cancer patients. Finally, NNMT expression was predictive of response to bevacizumab in ovarian cancer patients.
Recently, NNMT expression has been associated with a mesenchymal signature in bladder cancer and renal cancer [7, 8]. Consistently, our results also indicated that NNMT expression was increased in mesenchymal molecular subtype of ovarian cancer compared with other molecular subtypes of ovarian cancer. Functional studies revealed that knockdown of NNMT expression resulted in reduced migration and invasion of various cancer cells, suggesting an pro-metastatic role of NNMT [7, 8, 18, 19]. Consistent with these findings, we also found that the mean expression value of NNMT was increased in metastatic tumor tissues compared with primary counterparts, suggesting a potential pro-metastatic role of NNMT in ovarian cancer. But the differences were not statistically significant.
Previous studies have connected the increased expression of NNMT in cancer patients with chemoresistance . The inhibition of NNMT expression in colorectal cancer HT-29 cells diminishes 5-FU resistance, while upregulation of NNMT in SW480 cells enhances it . Our results also indicated that the mean expression value of NNMT was increased in drug resistant patients compared with drug sensitive patients, suggesting a potential link between NNMT elevation and drug resistance. But the differences were not statistically significant.
The prognostic value of NNMT has also been observed in various cancers. For example, increased expression of NNMT was associated with a worse prognosis in gastric cancer  and pancreatic cancer . Our data also showed that NNMT abundance was associated with aggressive disease characteristics such as advanced stage and poor differentiation. More importantly, increased NNMT expression was linked to a reduced OS. Moreover, our NNMT-related nomogram showed that NNMT shared a larger contribution to OS, compared with debulking status. Previously, NNMT was identified as a novel serum marker colorectal cancer. Evaluation of the possibility of NNMT as serum marker for ovarian cancer is also needed.
Our previous work indicated that the presence of LVSI was associated with a worse clinical outcome in ovarian cancer . LVSI is the first step of hematogenous metastasis of cancer cells , and we speculated that the expression of LVSI associated genes might be predictive of response to anti-angiogenesis therapy. Our work revealed that NNMT expression was higher in ovarian cancer tissues with LVSI than those without LVSI . Thus, we explored whether the expression of LVSI associated gene, NNMT, was predictive of response to bevacizumab, an anti-angiogenesis agent which has been added to standard first-line chemotherapy in ovarian cancer [25, 26]. Interestingly, our results indicated that bevacizumab conferred significant improvements in OS for patients with low NNMT expression. In contrast, patients with high NNMT expression didn’t benefit from bevacizumab treatment significantly. Therefore, more survival benefit from bevacizumab treatment in ovarian cancer may be expected if it is used in patients with lower NNMT expression. And in patients with higher NNMT expression, combined therapy with NNMT inhibitor and bevacizumab may be translated into survival improvements.
Previous studies have indicated that NNMT is a master regulator in cancer cell metabolism in general and Warburg effect in particular [13, 14]. NNMT was an up-regulator of PGAM1, which was involved in Warburg effect . Our data also demonstrated that NNMT expression was positively correlated with PGAM1 expression in ovarian cancer. Consistently, PGAM1 was also associated with bevacizumab response in ovarian cancer. This finding link NNMT expression to Warburg effect in particular. It may be helpful to clarify the mechanism of bevacizumab response in ovarian cancer.
In conclusion, NNMT expression is associated with the aggressive behavior of ovarian cancer, correlates with a poor prognosis, and is predictive of sensitivity to bevacizumab treatment. Future researches are needed to establish the role of NNMT in ovarian cancer development, and to explore whether combined therapy with NNMT inhibitor and bevacizumab could be translated into survival improvements.
Availability of data and materials
The datasets analysed during the current study are available in GEO database and TCGA database.
Lactate dehydrogenase A
Phosphoglycerate mutase 1
Fallopian tube epithelium
Primary ovarian cancers
Metastatic ovarian cancers
Recurrent ovarian cancers
Lymphovascular space invasion
Coleman RL, Monk BJ, Sood AK, Herzog TJ. Latest research and treatment of advanced-stage epithelial ovarian cancer. Nat Rev Clin Oncol. 2013;10(4):211–24.
Aksoy S, Szumlanski CL, Weinshilboum RM. Human liver nicotinamide N-methyltransferase. cDNA cloning, expression, and biochemical characterization. J Biol Chem. 1994;269(20):14835–40.
Roessler M, Rollinger W, Palme S, Hagmann ML, Berndt P, Engel AM, et al. Identification of nicotinamide N-methyltransferase as a novel serum tumor marker for colorectal cancer. Clin Cancer Res. 2005;11(18):6550–7.
Tomida M, Mikami I, Takeuchi S, Nishimura H, Akiyama H. Serum levels of nicotinamide N-methyltransferase in patients with lung cancer. J Cancer Res Clin Oncol. 2009;135(9):1223–9.
Kim J, Hong SJ, Lim EK, Yu YS, Kim SW, Roh JH, et al. Expression of nicotinamide N-methyltransferase in hepatocellular carcinoma is associated with poor prognosis. J Exp Clin Cancer Res. 2009;28:20.
Kassem H, Sangar V, Cowan R, Clarke N, Margison GP. A potential role of heat shock proteins and nicotinamide N-methyl transferase in predicting response to radiation in bladder cancer. Int J Cancer. 2002;101(5):454–60.
Wu Y, Siadaty MS, Berens ME, Hampton GM, Theodorescu D. Overlapping gene expression profiles of cell migration and tumor invasion in human bladder cancer identify metallothionein 1E and nicotinamide N-methyltransferase as novel regulators of cell migration. Oncogene. 2008;27(52):6679–89.
Tang SW, Yang TC, Lin WC, Chang WH, Wang CC, Lai MK, et al. Nicotinamide N-methyltransferase induces cellular invasion through activating matrix metalloproteinase-2 expression in clear cell renal cell carcinoma cells. Carcinogenesis. 2011;32(2):138–45.
Ulanovskaya OA, Zuhl AM, Cravatt BF. NNMT promotes epigenetic remodeling in cancer by creating a metabolic methylation sink. Nat Chem Biol. 2013;9(5):300–6.
Palanichamy K, Kanji S, Gordon N, Thirumoorthy K, Jacob JR, Litzenberg KT, et al. NNMT silencing activates tumor suppressor PP2A, inactivates oncogenic STKs, and inhibits tumor forming ability. Clin Cancer Res. 2017;23(9):2325–34.
Xie X, Yu H, Wang Y, Zhou Y, Li G, Ruan Z, et al. Nicotinamide N-methyltransferase enhances the capacity of tumorigenesis associated with the promotion of cell cycle progression in human colorectal cancer cells. Arch Biochem Biophys. 2014;564:52–66.
Kanska J, Aspuria PP, Taylor-Harding B, Spurka L, Funari V, Orsulic S, et al. Glucose deprivation elicits phenotypic plasticity via ZEB1-mediated expression of NNMT. Oncotarget. 2017;8(16):26200–20.
Eckert MA, Coscia F, Chryplewicz A, Chang JW, Hernandez KM, Pan S, et al. Proteomics reveals NNMT as a master metabolic regulator of cancer-associated fibroblasts. Nature. 2019;569(7758):723–8.
Cui Y, Yang D, Wang W, Zhang L, Liu H, Ma S, et al. Nicotinamide N-methyltransferase decreases 5-fluorouracil sensitivity in human esophageal squamous cell carcinoma through metabolic reprogramming and promoting the Warburg effect. Mol Carcinog. 2020;59(8):940–54.
Tan TZ, Yang H, Ye J, Low J, Choolani M, Tan DS, et al. CSIOVDB: a microarray gene expression database of epithelial ovarian cancer subtype. Oncotarget. 2015;6(41):43843–52.
Li J, Yue H, Yu H, Lu X, Xue X. Development and validation of SIRT3-related nomogram predictive of overall survival in patients with serous ovarian cancer. J Ovarian Res. 2019;12(1):47.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Mobley A, Zhang S, Bondaruk J, Wang Y, Majewski T, Caraway NP, et al. Aurora kinase a is a biomarker for bladder cancer detection and contributes to its aggressive behavior. Sci Rep. 2017;7:40714.
Chen C, Wang X, Huang X, Yong H, Shen J, Tang Q, et al. Nicotinamide N-methyltransferase: a potential biomarker for worse prognosis in gastric carcinoma. Am J Cancer Res. 2016;6(3):649–63.
Xie X, Liu H, Wang Y, Zhou Y, Yu H, Li G, et al. Nicotinamide N-methyltransferase enhances resistance to 5-fluorouracil in colorectal cancer cells through inhibition of the ASK1-p38 MAPK pathway. Oncotarget. 2016;7(29):45837–48.
Xu Y, Liu P, Zheng DH, Wu N, Zhu L, Xing C, et al. Expression profile and prognostic value of NNMT in patients with pancreatic cancer. Oncotarget. 2016;7(15):19975–81.
Li J, Li S, Chen R, Lu X. Increased risk of poor survival in ovarian cancer patients with high expression of SNAI2 and lymphovascular space invasion. Oncotarget. 2017;8(6):9672–85.
Yeung TL, Leung CS, Yip KP, Au YC, Wong ST, Mok SC. Cellular and molecular processes in ovarian cancer metastasis. A review in the theme: cell and molecular processes in Cancer metastasis. Am J Physiol Cell Physiol. 2015;309(7):C444–56.
Yue H, Wang J, Chen R, Hou X, Li J, Lu X. Gene signature characteristic of elevated stromal infiltration and activation is associated with increased risk of hematogenous and lymphatic metastasis in serous ovarian cancer. BMC Cancer. 2019;19(1):1266.
Perren TJ, Swart AM, Pfisterer J, Ledermann JA, Pujade-Lauraine E, Kristensen G, et al. A phase 3 trial of bevacizumab in ovarian cancer. N Engl J Med. 2011;365(26):2484–96.
Burger RA, Brady MF, Bookman MA, Fleming GF, Monk BJ, Huang H, et al. Incorporation of bevacizumab in the primary treatment of ovarian cancer. N Engl J Med. 2011;365(26):2473–83.
We thank all the investigators who shared their high throughput expression data.
This study was supported by funding from the National Natural Science Foundation of China (No.81802596), and Natural Science Foundation of Shanghai (20ZR1408900). The funders have no roles in writing and publication of the manuscript.
Ethics approval and consent to participate
All data in this article are obtained from published studies, which have received ethics approval and consent to participate. The study was performed in accordance with the Declaration of Helsinki.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, J., Yue, H., Yu, H. et al. Patients with low nicotinamide N-methyltransferase expression benefit significantly from bevacizumab treatment in ovarian cancer. BMC Cancer 21, 67 (2021). https://doi.org/10.1186/s12885-021-07785-w