Abstract
Aim: The aim of the study was to identify novel biomarkers that are vital for improving management of patients with gastric cancer (GC). Materials and Methods: An RNA-sequencing analysis was conducted using gastric tissue from patients with metastatic GC. In vitro cell functions were evaluated by siRNA-mediated knockdown assays. A total of 230 pairs of gastric tissue were subjected to expression analysis of mRNA and protein in situ. The serum levels of the candidate biomarker were determined by ELISA. Results: MELTF was identified as a candidate biomarker. Inhibition of MELTF expression suppressed the invasion ability of GC cells. Increased tissue MELTF mRNA expression was associated with shorter survival. Furthermore, staining intensity of tissue MELTF protein was linked to recurrence rates. Serum MELTF levels gradually were increased from healthy controls to advanced GC. Patients with high serum MELTF levels had poor prognosis. Conclusion: Both tissue and serum MELTF levels may serve as biomarkers of GC progression.
Despite important advances in multidisciplinary treatment, prognosis of patients with gastric cancer (GC) remains poor (1, 2). Discovery of biomarkers for accurate staging, monitoring of recurrence and prognostication contribute to improvement of patient management.
Since GC is recognized as a heterogeneous disease both phenotypically and genetically, currently available biomarkers are inadequate (3). Particularly, biomarkers detectable in blood samples are of practical value, because they can be measured repeatedly and non-invasively (4-7). However, conventional serum tumor markers, such as carbohydrate antigen (CA) 19-9, CA125 and carcinoembryonic antigen, may be widely used in clinical practice but are not ideal markers due to insufficient sensitivity and specificity in the risk stratification of advanced GC (8-10).
In this context, we sought to identify novel candidate genes involved in GC progression. Using an RNA-sequencing analysis, we identified melanotransferrin (MELTF), which is a cell-surface glycoprotein (11). The aim of this study was to determine the oncological role of MELTF and significance of its expression in gastric tissue and sera as a biomarker of GC progression.
Materials and Methods
RNA-sequencing. To identify candidate genes associated with the metastatic potential of GC cells, surgically resected gastric tissue from four patients with distant metastasis of GC was subjected to RNA-sequencing analysis (HiSeq platform; Illumina, San Diego, CA, USA). We compared the expression levels of 57,749 genes in primary GC tissues with those of the corresponding normal gastric mucosa (12).
Cell lines and clinical samples. Fourteen GC cell lines (AGS, GCIY, IM95, KATOIII, MKN1, MKN7, MKN45, MKN74, N87, NUGC2, NUGC3, NUGC4, OCUM1 and SC-6-JCK) and a control, nonmalignant small intestinal epithelial cell line (FHs74) were obtained from the American Type Culture Collection (ATCC; Manassas, VA, USA) or the Japanese Collection of Research Bioresources Cell Bank (JCRB; Osaka, Japan). A total of 230 pairs of gastric tissue (tumor and the corresponding normal mucosa) were obtained from patients who underwent complete gastrectomy without neoadjuvant therapy between 2001 and 2014 as previously described (7). Since 2006, postoperative adjuvant therapy using S-1 has been recommended for stage II or III GC patients unless contraindicated (13). As an external validation cohort, a public dataset comprising 1065 patients with GC was used (14, 15).
List of genes overexpressed in primary cancerous tissues from patients with metastatic gastric cancer.
Transcription of MELTF and PCR array analysis. The levels of MELTF mRNA were analyzed using a real-time quantitative reverse-transcription PCR (RT-PCR) analysis. As a reference gene, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used (16). To investigate the involvement of MELTF by identifying cancer-related genes expressed simultaneously with MELTF, PCR array analysis was performed using the human epithelial to mesenchymal transition (EMT) RT2 Profiler PCR Array (Qiagen, Hilden, Germany) (17).
Knockdown of MELTF expression. MELTF was silenced by transfection with small interfering RNA (siRNA) as described previously using LipoTrust EX Oligo (Hokkaido System Science, Sapporo, Japan) (17). To confirm MELTF knockdown, both mRNA and protein expression of MELTF were analyzed. A rabbit polyclonal antibody against MELTF (NBP1-85777PEP; Novus Biologicals, Centennial, CO, USA; dilution 1:1000) was used for western blot analysis (18).
Cell functional analyses. Cell proliferation was evaluated with the Cell Counting Kit-8 (Dojindo Molecular Technologies, Inc., Kumamoto, Japan) (19, 20). We determined the invasive ability of GC cells using BioCoat Matrigel invasion chambers (BD Biosciences, Bedford, MA, USA). Wound-healing assay was performed to evaluate cell migration (17). The CytoSelect 48-Well Cell Adhesion Assay (Cell Biolabs, Inc., San Diego, CA) was employed to analyze adherence of GC cells to the extracellular matrix components (18).
Immunohistochemistry (IHC). In situ MELTF protein expression was determined with IHC in 84 representative patients who had well-preserved sections of GC tissues. Formalin-fixed paraffin-embedded sections were incubated for 1 h at room temperature with a rabbit polyclonal antibody against MELTF (NBP1-85777PEP, Novus Biologicals) diluted 1:100. We graded the patients according to the staining intensity of MELTF in GC component: none, weak, moderate, and strong staining as previously described (21).
Enzyme-linked immunosorbent assay (ELISA). The serum levels of MELTF were determined in 20 healthy controls and 68 patients with GC. Adults (n=20, median age 42 years) in the healthy controls had annual medical checkups (chest X-rays, blood tests and physical examination) and were not under disease treatment. All patients were diagnosed as adenocarcinoma of the stomach by upper gastrointestinal endoscopy before blood collection for the ELISA assay. We collected blood samples within 7 days before gastrectomy. Serum was immediately separated and stored at −80°C. MELTF ELISA Kit (CSB-EL013754HU, CUSABIO, Wuhan, PR China) was used and all samples were assayed in duplicate (22).
Statistical analyses. Differences in the data between the two groups were analyzed by the Mann–Whitney test or χ2 test. The receiver-operating characteristic (ROC) curve was used to assess goodness-of-fit by calculating the area under the curve (AUC). Survival data was analyzed with Kaplan–Meier curves. To identify prognostic factors, multivariable regression analysis was performed using the Cox proportional hazards model, and differences between variables with p<0.05 were entered in the final model. JMP 13 software (SAS Institute Inc., Cary, NC, USA) was used for all statistical analyses. A p-value of less than 0.05 was considered to indicate statistically significant differences.
Results
Detection of MELTF. As a result of an RNA-sequencing analysis, 25 candidates identified that were highly expressed in GC tissues compared to the noncancerous tissue and comparably expressed in primary and metastatic tissue (Table I). Of them, we selected MELTF for further investigation, because it encodes for a cell-surface glycoprotein that could be released into the circulatory system, and no previous reports were found regarding GC.
Expression in GC cell lines. The MELTF mRNA levels were fivefold higher in MKN7, NUGC4, OCUM1, AGS, MKN1, and N87 cells than in FHs74 cells. No clear trends in MELTF mRNA levels were indicated based on the differentiation of the GC cell lines. Three epithelial mesenchymal transition-related genes encoding vacuolar protein sorting 13 homolog A (VPS13A), snail family transcriptional repressor 3 (SNAI3), and Erb-B2 receptor tyrosine kinase 3 (ERBB3) were expressed at levels with correlation coefficient >0.6 compared with MELTF (Figure 1a).
Effect of MELTF knockdown on biological activities of gastric cancer cells. N87, AGS, and MKN1 were differentiated-type cells displaying high levels of MELTF mRNA expression. NUGC4 and OCUM1 were undifferentiated-type cells also displaying high MELTF mRNA expression. Based on the characteristics of the cells (e.g. invasion and migration abilities), we selected MKN1 and AGS from the differentiated-type cells, and NUGC4 from the undifferentiated-type cells for the knockdown experiments. To evaluate the oncological functions of MELTF in GC cells, MELTF was silenced in MKN1, AGS, and NUGC4 cells by specific siRNAs. After confirming knockdown efficiency of mRNA and protein expression (Figure 1b), we compared cell proliferation, invasiveness, migration, and adhesion of MKN1 and AGS cells with and without MELTF knockdown. Inhibition of MELTF expression significantly decreased the proliferation of MKN1 cells (Figure 1c) and invasion of MKN1 and AGS cells (Figure 2a). Furthermore, the migration of AGS cells was reduced by siMELTF (Figure 2b). However, the silencing of MELTF had little effect on the adhesive ability (Figure 2c).
Prognostic value of tissue MELTF mRNA expression. There were 165 males and 65 females aged 26 to 96 years [66.3±10.8 years, mean±standard deviation (SD)]. With respect to pathological stage, 50, 71, and 109 patients were in stages I, II, and III, respectively, indicating that all 230 patients underwent curative gastrectomy. GC tissue samples had significantly higher mean levels of MELTF mRNA than the adjacent normal tissues. The AUC value of MELTF levels was 0.749 for detection of cancer-related death within 2 years after gastrectomy, and the optimal cutoff value was 0.00248 (sensitivity 91%, specificity 55%). Patients were classified into the low MELTF (below the cut-off value) or high MELTF (equal to or above the cut-off value) groups. The high MELTF group had a significantly lower overall survival rate than the low MELTF group (5-year survival rates, 78% and 89%, respectively; Figure 3a). Furthermore, the disease-free survival rate was significantly lower in the high MELTF group than in the low MELTF group (3-year disease-free survival rates, 64% and 89%, respectively, p<0.001). The frequencies of peritoneal and lymph node recurrences were significantly higher in the high MELTF group than in the low MELTF group, whereas there was no significant difference in that of hematogenous recurrences. Consistent results were obtained in the analysis of the extra-validation cohort (Figure 3b).
Next, we investigated the correlation between MELTF expression and clinicopathological characteristics in patients. High MELTF expression in GC tissue was significantly correlated with differentiated tumor and vascular invasion. High MELTF mRNA expression in GC tissue was found to be a significant prognostic factor together with CA19-9 >37 IU/ml, larger tumor size, pT4, invasive growth, lymphatic involvement, vascular invasion and nodal metastasis in the univariate analysis of overall survival. Multivariate analysis identified lymph nodal metastasis and high MELTF mRNA levels as independent prognostic factors for overall survival (Table II).
Expression of MELTF protein in situ using immunohistochemistry. Examples of each staining intensity are shown in Figure 3a. Forty-three, 20, 18, and 3 patients were categorized as none, weak, moderate, and strong staining, respectively. Prevalence of recurrence after surgery gradually increased as the intensity of staining of MELTF increased (Figure 3c).
Prognostic value of serum MELTF levels. Fifty-one males and 17 females were included for the serum analysis and they were aged 24 to 88 years (65.1±12.6 years, mean±SD). Based on the UICC TNM Classification of Malignant Tumors, 8th edition, 17, 15, 16, and 20 patients were in stage I, II, III, and IV, respectively. Twenty-seven and 41 patients were diagnosed with differentiated and undifferentiated GC. Serum MELTF levels of the 20 healthy controls ranged from 0 to 105.2 pg/ml (15.9±28.0 pg/ml, mean±SD). Patients with GC (n=68) had significantly higher preoperative serum MELTF levels (388±1834 pg/ml, mean±SD) than the healthy controls. Serum MELTF levels gradually increased along with disease stages (Figure 3d). We set the cut-off of serum MELTF levels at 105.2 pg/ml, which were the highest levels in the healthy controls. The sensitivity, specificity, false positive, and false negative at the cutoff value of MELTF 105.2 pg/ml were 50%, 100%, 0% and 50%, respectively. Using this cut-off, GC patients were classified into the MELTF-positive (equal to or above the cut-off) and MELTF-negative (less than the cut-off) groups. Patients in the MELTF-positive group had significantly shorter survival than those in the MELTF-negative group (Figure 3e).
Expression analysis of cell lines and effects of siRNA-mediated knockdown of MELTF expression. (a) Expression of MELTF and genes at similar differential levels were identified using polymerase chain reaction (PCR) array analysis. MELTF mRNA levels positively correlated with those of SNAI3, VPS13A, and ERBB3. (b) Silencing of MELTF in MKN1 and AGS cells was examined using PCR and western blot analysis. (c) Cell proliferation of MKN1 and AGS cells transfected with siMELTF, control siRNA, or untreated. *p<0.05. Error bars indicate standard deviation. siRNA, Small interfering RNA.
Prognostic factors for overall survival in 230 patients with gastric cancer.
Discussion
RNA-sequencing analysis using tissue from patients with metastatic GC, identified MELTF as a candidate biomarker. Expression levels of MELTF were correlated with those of three EMT-related genes, SNAI3, VPS13A, and ERBB3, and knockdown of MELTF expression attenuated proliferation, invasion, and migration abilities of GC cells. Tissue and serum MELTF levels were associated with disease recurrence and prognosis.
MELTF is a cell-surface glycoprotein initially found on melanoma cells, and it belongs to the transferrin superfamily (11, 23). MELTF is related to transferrin in tertiary structure and has the physiological role of iron cellular uptake (11, 24, 25). In malignancy, it has been reported that MELTF is significantly overexpressed in melanoma and associated with increased proliferation of melanoma cells (11, 25, 26). Additionally, recent reports have confirmed MELTF as a potential serological marker for colorectal cancer (27, 28). However, the role of MELTF in GC is poorly understood.
PCR array analysis showed that SNAI3, VPS13A, and ERBB3 were expressed together with MELTF with high correlation coefficient. SNAI3 plays a role in mesodermal formation during embryogenesis and acts as a transcriptional repressor implicated in the EMT process (29, 30). VPS13A is reported to control the process of cycling through the trans-Golgi network and the plasma membrane, lysosomes and endosomes (16, 31). It has been reported that VPS13A is expressed in various cancer cells (32). ERBB3 is a member of the epidermal growth factor receptor (EGFR) family and has been reported that this gene is overexpressed in GC (33-35). We showed here that MELTF contributes to the proliferation, invasion, and migration ability of GC cells. Our present findings support the conclusion that MELTF acts as a putative oncogene in GC.
In our cohort, high expression of tissue MELTF mRNA was significantly associated with shorter disease-free and overall survival after curative resection. These findings are consistent with the findings from our analysis of a validation cohort (14). Multivariate analysis showed that high mRNA levels of MELTF are an independent prognostic factor for survival. Moreover, expression status of tissue MELTF can be evaluated by immunohistochemical staining. These results highlighted the value of MELTF expression as a biomarker for GC progression.
Next, we detected MELTF in serum samples, because blood sampling is less invasive and readily accessible (6, 36). MELTF levels in serum samples of GC patients were significantly higher than that in healthy controls. Because preoperative serum MELTF levels gradually increased along with higher UICC stage and were associated with poor overall survival, examining preoperative MELTF serum levels may help predict patient prognosis. Shin et al. have conducted quantitative proteomic analysis to identify serum biomarkers for colorectal cancer and confirmed MELTF as a potential serological marker (27). They measured levels of MELTF by ELISA in 325 samples from 77 healthy controls and 228 patients with colorectal cancer and found that serum MELTF levels were elevated in patients with colorectal cancer (27). Detection of MELTF in serum samples is an ideal method for screening and emphasizes the clinical practicality of evaluating MELTF expression.
Effect of inhibition of MELTF on MKN1 and AGS cells. (a) Cell invasion assay of MKN1 and AGS cells transfected with siMELTF, control siRNA, or untreated (200× magnification). (b) Migration assays in MKN1 and AGS cells transfected as indicated. Graph shows mean length of migration at the indicated times. (c) Cell adhesion assay in MKN1 and AGS cells transfected as indicated. *p<0.05. Error bars indicate standard deviation. siRNA, Small interfering RNA.
Overall survival of patients with stage I-III gastric cancer according to tissue mRNA levels of MELTF in our institutional cohort (a) and the extra-validation cohort (b). (c) Immunohistochemical staining. Representative data showing none, weak, moderate, and strong staining for MELTF expression (400× magnification). The results show a positive correlation between the staining intensity of MELTF in gastric cancer component and frequency of recurrence. (d) Measurement of serum MELTF levels. Stepwise elevation of serum MELTF levels in healthy controls and patients with gastric cancer. *p<0.05. Error bars indicate standard deviation. (e) Patients in the MELTF-positive group had poorer prognosis compared with those in the MELTF-negative group.
Several limitations of this study should be acknowledged. Further pathway analysis is desirable to understand the oncological roles of MELTF in GC. The relatively small sample size limited the statistical power. Unfortunately, very limited cases experienced disease recurrences after measurement of serum MELTF levels. Therefore, it is difficult to analyze correlations between serum MELTF levels and recurrent patterns this time.
To further evaluate the diagnostic performance of serum MELTF levels, a prospective study to determine if serum MELTF levels normalize after resection of primary gastric cancer and rise prior to radiographic evidence of disease recurrences is needed.
In conclusion, MELTF acts as a putative oncogene in GC and its levels in primary GC tissues and serum samples may be a potential biomarker.
Acknowledgements
The Authors thank Edanz Group (www.edanzediting.com/ac) for editing a draft of this article.
Footnotes
Authors' Contributions
KS and MK conceived the study concept and design, analyzed data and wrote the manuscript. MK, SU, TM, CT, DK, MH, SY, GN, KO, MK and YK contributed to data acquisition and interpretation. KO contributed to statistical analysis. CT, DK, MH, SY, GN, MK, YK and KO revised the draft. All authors have read and approved the final version of the manuscript.
Conflicts of Interest
The Authors have no conflicts of interest directly relevant to the content of this article.
- Received October 9, 2019.
- Revision received October 21, 2019.
- Accepted October 22, 2019.
- Copyright© 2019, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved