The association of CMTM6 expression with prognosis and PD-L1 expression in triple-negative breast cancer
Introduction
Breast cancer is the most common malignancy in women worldwide and its clinical and molecular heterogeneity is well-documented. Endocrine receptors for estrogen (ER) and progesterone (PR), and the aberrant expression of human epidermal growth factor receptor-2 (HER2) are the specific biomarkers for breast cancer most commonly used in clinical practice. Triple-negative breast cancer (TNBC), which lacks expression of ER, PR, and HER2, is found in 10–20% of all breast cancers. TNBC is an aggressive type of breast cancer that usually displays a higher grade and poorer outcome than other breast cancer subtypes (1,2). Therefore, effective therapeutic strategies for TNBC are urgently needed.
Programmed death-ligand 1 (PD-L1) is a ligand of the programmed cell death-1 (PD-1) receptor and can be expressed by the tumor cell surface as well as by tumor infiltration lymphocytes (TILs). Binding of PD-L1 to PD-1 can downregulate anti-tumor T cell responses which lead to tumor immune escape (3). To date, immune checkpoint inhibitors that block the interaction between PD-1 and PD-L1 have been used to treat many types of metastatic cancers (3-5). In primary breast cancer (PBC), the expression of PD-L1 is heterogeneous and associated with higher histological grades and more aggressive molecular subtypes [triple-negative (TN), basal, and HER2-driven] (6,7). However, the regulation mechanism of PD-L1 remains elusive.
Recently, chemokine-like factor (CKLF)-like Marvel Transmembrane Domain-containing 6 (CMTM6) has been identified as a key regulator of the PD-L1 protein, which is also thought to be involved in modulating tumor immunity (8,9). Dysfunction of CMTM6 impairs the expression of PD-L1 in many human tumor cell types, such as lung cancer, thyroid cancer, and melanoma (8). The depletion of CMTM6, via downregulated PD-L1, can reduce the suppression of tumor-specific T cell activity in vitro and in vivo assays (10). Previous studies have suggested the potential value of a therapeutic target that elicits an immune response and avoids the escape of immune surveillance. However, the role of CMTM6 in breast cancer remains unclear.
In the present study, we used immunohistochemistry (IHC) and real-time polymerase chain reaction (RT-PCR) to analyze the correlation between CMTM6 and PD-L1 in patients with TNBC. Furthermore, we assessed the association of CMTM6 or PD-L1 expression with patients’ prognosis. This study gives a hint that the expression of CMTM6 may become a potential biomarker of immunotherapies through supplement to PD-L1 expression and could predict the patients’ outcome. We present the following article in accordance with the REMARK reporting checklist (available at http://dx.doi.org/10.21037/atm-20-7616).
Methods
Patients
A total of 143 patients with breast cancer who underwent surgical resection in the Cancer Hospital of the University of Chinese Academy of Sciences in China between January 2008 and December 2015 were enrolled in the study. Of these, 83 patients who had paired primary breast cancer tissues and matched metastatic cancer tissues were selected and comprised 45 HER2-driven breast cancer and 38 TNBC patients. All of these patients were without neoadjuvant chemotherapy at the first diagnosis and had sufficient archival tissue in formalin-fixed, paraffin-embedded (FFPE) blocks to perform IHC analyses. The remaining 60 individuals comprised 30 HER2-driven and 30 TNBC PBC patients with sufficient frozen tissue to perform mRNA expression detection (Figure S1).
The eighth edition of the tumor node metastasis (TNM) classification of the American Joint Commission on Cancer (AJCC) was used as a reference for pathological features and clinical stage. Progression-free survival (PFS) was defined as the time from the date of the initial treatment until the date of diagnosis of the initial recurrence or death from any cause. The PFS data were locked on August 30, 2020.
This study was approved by the Medical Ethics Committee of Zhejiang Cancer Hospital (IRB-2020-275). The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). Because of the retrospective nature of the research, the requirement for informed consent was waived.
IHC staining of CMTM6 and PD-L1
FFPE tissue specimens were collected from 83 patients with paired PBC and metastatic cancer tissues, and an IHC assay was performed to detect the expression of CMTM6 and PD-L1. FFPE tissue blocks were cut into 4-µm-thick sections and attached to a positively charged glass slide. Immunohistochemical staining of CMTM6 was carried out with a Leica BOND-III automatic IHC staining device. The samples were incubated with an antibody against CMTM6 for 15 min at room temperature (dilution 1:200, recombinant monoclonal antibody, Abcam, EPR23015-45, US). The signal was subsequently detected with a Leica Bond Polymer Refine Detection Kit (Leica Biosystems, Buffalo Grove, US). Post-primary antibody incubation and polymer incubation were set at 8 min for CMTM6 and followed by immersion in diaminobenzidine (DAB) for signal visualization for 10 min at room temperature. CMTM6 expression was defined as any intensity in the cytomembrane on the tumor cells (Figure S2A,B).
For PD-L1 staining, the VENTANA PD-L1 (SP142) rabbit monoclonal primary antibody (Ventana Medical Systems Inc., Tucson, AZ, USA) was optimized for use as a fully automated IHC assay on the BenchMark ULTRA (Ventana Medical Systems Inc., Tucson, AZ, USA) staining platform using the OptiView DAB IHC Detection Kit and OptiView Amplification Kit (Ventana Medical Systems Inc., Tucson, AZ, USA) as previously described (11). PD-L1 expression was identified as any intensity in the cytomembrane on the tumor-infiltrating immune cells (Figure S2C,D).
RNA extraction and quantitative RT-PCR
Total RNA was isolated from 60 frozen cancer tissues using TRIzol (Invitrogen, USA), following the manufacturer’s protocol. The concentration of RNAs was measured through a microvolume spectrophotometer. The 500 nanograms of RNA was reverse transcribed using a PrimeScript™ RT Reagent Kit (Takara, Dalian, China). Quantitative PCR (qPCR) was carried out as previously described (12). The primer sequences were: CMTM6-F: 5'-GCAACAATATCAGCAACTTCGT-3' and CMTM6-R: 5'-TTGGTCCTTAGGTGTGGTATCA-3'; PD-L1-F: 5'-CACCACCACCAATTCCAAGAG-3'; PD-L1-R: 5'- AGGATGTGCCAGAGGTAGTTC-3'; β-actin-F: 5'-TGGCACCCAGCACAATGAA-3', β-actin-R: 5'-CTAAGTCATAGTCCGCCTAGAAGCA-3'.
Western blotting
Cells were collected and prepared as described earlier (13). Proteins were separated by 10% SDS-PAGE and transferred to polyvinylidene fluoride membrane (Millipore, Bedford, MA, USA). Membranes were blocked in 5% non-fat milk in TBS-Tween 20 for 2 hours prior to overnight incubation with primary antibodies at 4 °C, and were then incubated for 2 hours at room temperature with the secondary antibodies (Servicebio, Wuhan, China). Finally, the protein level was detected using an enhanced chemiluminescence (ECL) reagent (CwBio, Beijing, China). Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used for normalization. Primary antibodies were directed against CMTM6 (Abcam, USA, 1:1,000), PD-L1 (Proteintech, USA, 1:1,000), and GAPDH (Santa Cruz, CA, USA, 1: 5,000).
Web-based mRNA profiling
The RNA-sequencing data and clinical records of 30 HER2-driven PBC and 123 triple-negative PBC patients were downloaded from The Cancer Genome Atlas (TCGA) dataset (see URLs http://cancergenome.nih.gov/). The present study conforms to the publication guidelines.
Statistical analysis
Statistical analysis was performed using SPSS statistics version 22.0 (IBM Corp., Armonk, NY, USA). The continuous variables were described as means and standard deviations (SD), and the categorical variables were described as number (percentage). Any skewed distribution data were expressed with median and interquartile ranges. The Mann-Whitney test was adopted to evaluate the non-normal distribution dataset. Spearman’s correlation analysis was used to describe the correlation between quantitative variables with skewed distribution data, and the Pearson Chi-square and McNemar’s tests were used to compare categorical variables. Kaplan-Meier (K-M) analysis was used, and between-group differences were evaluated by the log-rank test. The multivariable Cox proportional hazard model was applied to calculate the adjusted hazard ratio (HR). A two-sided P value <0.05 was considered statistically significant.
Results
Patient characteristics
A total of 143 patients with breast cancer, including 75 (52.4%) HER2-driven PBCs and 68 (47.6%) triple-negative PBCs were included in the study. The clinical features are shown in Table 1. All patients were female, with a mean age of 49.97±10.63 years. Of these, 16 (11.2%) had stage I, 60 (42.0%) had stage II, 47 (32.9%) had stage III, and 20 (14.0%) had stage IV diagnoses of breast cancer. A total of 97 (67.8%) patients showed lymph node metastases. There were 17 (11.9%) patients with a tumor size smaller than 2 cm, 96 (67.1%) patients with a tumor size of 2–5 cm, and 30 (21.0%) patients with a tumor size larger than 5 cm. Of the total number of 143 patients, 83 had matched PBC and metastatic breast cancer (MBC). There was a 97.6% (81/83) concordance between ER, PR, and HER2 expression between the matched PBC and MBC patients. The median PFS of these 83 patients was 48 months (range, 29–76 months).
Full table
CMTM6 and PD-L1 expression in HER2-driven breast cancer and TNBC
IHC was performed to test the expression of CMTM6 and PD-L1 in the 83 paired breast cancers. The typical CMTM6 and PD-L1 staining images from breast cancer samples are shown in Figure S2. Western blot was used to confirm the expression detected by IHC (Figure S3). In PBC, 56 (67.5%) cases were negative for CMTM6, while 27 (32.5%) showed CMTM6-positive staining. The positive expression rate of CMTM6 in patients with triple-negative PBC was significantly higher than in those with HER2-driven PBC (63.0% vs. 37.0%, P=0.036, Table 2). A similar result was found in the MBC patients (61.1% vs. 38.9%, P=0.048, Table 2). Moreover, the mRNA expression level of CMTM6 showed the same result (the median CMTM6 mRNA expression was 0.037 in triple-negative PBC and 0.026 in HER2-driven PBC, P=0.036, Table S1, Figure S4A).
Full table
For the PD-L1 expression, 50 (60.2%) showed negative staining, while 33 (39.8%) showed positive staining. However, no significant difference in PD-L1 expression was found between HER2-driven breast cancer and TNBC, either in PBC or in MBC samples (Table 2). Moreover, a significant difference was found in the PD-L1 mRNA expression level between triple-negative PBC and HER2-driven PBC (the median expression level was 0.00026 and 0.00012, respectively, P=0.040, Table S1, Figure S4B).
To further verify our findings, we performed CMTM6 and PD-L1 mRNA expression analyses in HER2-driven PBC and triple-negative PBC using the TCGA dataset. Although the expression level of PD-L1 between HER-2 driven PBC and triple-negative PBC did not differ (P=0.283, Table S2), a higher expression level of CMTM6 was observed in triple-negative PBC (P=0.039, Table S2).
CMTM6 and PD-L1 expression in the matched PBC and MBC samples
As shown in Table 3, we evaluated whether there was a difference in the expression levels of CMTM6 and PD-L1 in the matched PBC and MBC samples. Across the 83 matched samples, patients with a PD-L1 positive expression in the MBC group were more likely to have a PD-L1 positive expression in the matched PBC group (P=0.031). A similar result was found in the triple-negative samples (P=0.039), while no significant difference was found for the HER2-driven samples. There was no correlation between the expression of CMTM6 in the matched breast cancer samples in this study.
Full table
Correlation of CMTM6 expression and PD-L1 expression in breast cancers
The correlation between the expression of CMTM6 and PD-L1 was further analyzed. In MBC, we found that CMTM6 protein expression was positively correlated with PD-L1 protein expression. McNemar’s test showed that the expression of CMTM6 and PD-L1 were significantly correlated across the entire MBC sample (P=0.012, Table 4), but especially in the triple-negative MBC sample (P=0.004, Table 4). However, no significant correlation was observed in the individual PBC or HER2-driven MBC groups (Table 4).
Full table
Next, we performed correlation analyses between CMTM6 and PD-L1 expression in the mRNA data using Spearman’s correlation coefficient analysis. The results showed that the mRNA expression level of CMTM6 was positively correlated with the mRNA expression level of PD-L1 in total breast cancers (r Spearman=0.419 and P=0.001, Figure 1A), in HER2-driven PBC (r Spearman=0.373 and P=0.042, Figure 1B), and in TNBC (r Spearman=0.421 and P=0.021, Figure 1C).
Clinicopathological characteristics by CMTM6 and PD-L1 expression
CMTM6 expression was not significantly associated with clinicopathological characteristics (Table 5) but was significantly associated with PFS. In the total sample, those with CMTM6 expression were at higher risk for disease progression compared with those without CMTM6 expression, especially in TNBC (HR =1.83, 95% CI: 1.13–2.96, P=0.014 for the total sample; and HR =2.22; 95% CI: 1.08–4.57, P=0.030 for TNBC, Table 6). Median PFS times for those with CMTM6 expression in the total sample and triple-negative cases were 46.8 (IQR, 38.3–50.0) and 47.0 (IQR, 39.5–52.7) months, respectively, while the median for those without CMTM6 expression was over 51 months. The Kaplan-Meier survival analysis demonstrated similar results (Figure 2A,B,C). After adjusting for age and TNM stage, significant associations were found in CMTM6 (HR =2.21, 95% CI: 1.32–3.71, P=0.003), especially in the TNBC patients (HR =2.40, 95% CI: 1.11–5.20, P=0.027, Table 6).
Full table
Full table
In the case of PD-L1, the absence of PD-L1 expression was correlated with lymph node metastasis (P=0.026, Table 5). No other clinicopathologic characteristics were associated with PD-L1 expression. Moreover, our data showed that the expression of PD-L1 was not significantly associated with disease progression (Figure 2D,E,F, Table 6).
Discussion
Our data suggested that the expression of CMTM6 was higher in TNBC than in HER2-driven breast cancer and that CMTM6 was also an independent risk factor of PFS, especially in TNBC. The expression level of CMTM6 was positively correlated with the expression of PD-L1. Moreover, the expression of PD-L1 in triple-negative metastatic breast cancer was positively correlated with that in primary breast cancer. And the expression of PD-L1 was negative with metastases in lymph nodes in PBC.
TNBC has a more aggressive biological behavior than other types of breast cancer. Patients with TNBC do not benefit from hormonal or trastuzumab-based therapy because of the loss of the target receptors HER2, ER, and PGR (2,14); chemotherapy and surgery appear to be the only available treatment modalities (15). Recent studies have shown that TNBC has a higher mutational burden than other subtypes and presents with tumor-infiltrating lymphocytes (TILs). Moreover, patients with TNBCs combined low-TILs and high PD-L1 expression showed unfavorable outcome, which may be benefit from immune therapy (1,7,16-18). According to these results, many clinical trials are currently evaluating the role of checkpoint inhibitors, PD-1, and/or PD-L1 in TNBC with encouraging results (19,20).
Although anticancer therapies based on immune checkpoint inhibitors targeting the PD-1/PD-L1 pathway have emerged as a critical treatment option, it seems that only a small subset of patients can benefit from these (20,21). To solve this challenge, some studies have focused on investigating the regulation of PD-L1 expression identifying a type-3 transmembrane protein, CMTM6, as a regulator of PD-L1 expression (8,10). In our study, the expression of CMTM6 was higher in TNBC than in HER2-driven breast cancer, according to the TCGA dataset. The expression of CMTM6 was positively correlated with PD-L1 in PBC and MBC in mRNA expression and IHC data, respectively. Koh et al. and Gao et al. reported similar results in non-small cell lung cancer (9,22). Moreover, our data revealed no significantly different expression of CMTM6 between paired triple-negative PBC and MBC. Considering the transient PD-L1 expression in tumor cells which tends to disappear rapidly, the intratumoral heterogeneity of PD-L1 expression is frequently observed (23). Thus, the expression of CMTM6 may be a new potential biomarker of immunotherapies in TNBC.
In our study, TNBC patients with a high expression of CMTM6 had a poorer PFS. CMTM6 could therefore serve as an independent risk factor to predict patient outcome. At the mRNA expression level, Mamessier et al. found that a CMTM6-high group had a shorter overall survival than a CMTM6-low group in pancreatic adenocarcinomas. They also found that CMTM6-high and PD-L1-high groups were associated with better metastasis-free survival in triple-negative PBC and that the CMTM6 expression enhanced the prognostic value of PD-L1 expression (24). Furthermore, some studies reported that individuals diagnosed with non-small cell lung cancer and hepatocellular carcinoma who showed a high expression of CMTM6 had better overall survival (22,25). On the other hand, some studies have provided evidence that a high expression of CMTM6 is associated with poor prognosis in gliomas and gastric cancer (26,27). Tumor heterogeneity and ethnic differences may be responsible for these conflicting results, and a larger cohort of a specific cancer type based on multicenter studies is recommended.
This study has several limitations. Firstly, this is a retrospective study with a relatively small sample size. A larger cohort of TNBC patients and multicenter studies are recommended for future studies. Secondly, some cases were from biopsy samples that showed few tumor-infiltrating immune cells, which may explain the lack of significant correlation between the CMTM6 and PD-L1 protein expression in primary TNBC. Moreover, it was not possible to evaluate the PD-L1 expression in cancer cells. Thirdly, the exact mechanism of CMTM6 regulation in TNBC remains unclear and requires further investigation.
In conclusion, and to the best of our knowledge, this is the first study that has reported a higher expression of CMTM6 in triple-negative compared to HER2-driven breast cancer. The expression of CMTM6 was correlated with PD-L1 in TNBC and could be used as an independent risk factor to predict PFS. The results from this study suggest that CMTM6 may be able to guide prognosis and would become a new potential biomarker of immunotherapies as a supplement to PD-L1 expression in TNBC.
Acknowledgments
Funding: This study was supported by the Medical Health Science and Technology Project of Zhejiang Provincial Health Commission (2018KY030) and the Natural Science Foundation of Zhejiang Province (LQ21H260001).
Footnote
Reporting Checklist: The authors have completed the REMARK reporting checklist. Available at http://dx.doi.org/10.21037/atm-20-7616
Data Sharing Statement: Available at http://dx.doi.org/10.21037/atm-20-7616
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at http://dx.doi.org/10.21037/atm-20-7616). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. This study was approved by the ethics committee of Zhejiang Cancer Hospital (No. IRB-2020-275). The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). Because of the retrospective nature of the research, the requirement for informed consent was waived.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Mittendorf EA, Philips AV, Meric-Bernstam F, et al. PD-L1 expression in triple-negative breast cancer. Cancer Immunol Res 2014;2:361-70. [Crossref] [PubMed]
- Dent R, Trudeau M, Pritchard KI, et al. Triple-negative breast cancer: clinical features and patterns of recurrence. Clin Cancer Res 2007;13:4429-34. [Crossref] [PubMed]
- Schutz F, Stefanovic S, Mayer L, et al. PD-1/PD-L1 Pathway in Breast Cancer. Oncol Res Treat 2017;40:294-7. [Crossref] [PubMed]
- Topalian SL, Hodi FS, Brahmer JR, et al. Safety, activity, and immune correlates of anti-PD-1 antibody in cancer. N Engl J Med 2012;366:2443-54. [Crossref] [PubMed]
- Brahmer JR, Tykodi SS, Chow LQ, et al. Safety and activity of anti-PD-L1 antibody in patients with advanced cancer. N Engl J Med 2012;366:2455-65. [Crossref] [PubMed]
- Bertucci F, Goncalves A. Immunotherapy in Breast Cancer: the Emerging Role of PD-1 and PD-L1. Curr Oncol Rep 2017;19:64. [Crossref] [PubMed]
- Zhang M, Sun H, Zhao S, et al. Expression of PD-L1 and prognosis in breast cancer: a meta-analysis. Oncotarget 2017;8:31347-54. [Crossref] [PubMed]
- Mezzadra R, Sun C, Jae LT, et al. Identification of CMTM6 and CMTM4 as PD-L1 protein regulators. Nature 2017;549:106-10. [Crossref] [PubMed]
- Gao F, Chen J, Wang J, et al. CMTM6, the newly identified PD-L1 regulator, correlates with PD-L1 expression in lung cancers. Biochem Biophys Rep 2019;20:100690. [Crossref] [PubMed]
- Burr ML, Sparbier CE, Chan YC, et al. CMTM6 maintains the expression of PD-L1 and regulates anti-tumour immunity. Nature 2017;549:101-5. [Crossref] [PubMed]
- Vennapusa B, Baker B, Kowanetz M, et al. Development of a PD-L1 Complementary Diagnostic Immunohistochemistry Assay (SP142) for Atezolizumab. Appl Immunohistochem Mol Morphol 2019;27:92-100. [Crossref] [PubMed]
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001;25:402-8. [Crossref] [PubMed]
- Zhang F, Ying L, Jin J, et al. GAP43, a novel metastasis promoter in non-small cell lung cancer. J Transl Med 2018;16:310. [Crossref] [PubMed]
- Lyons TG. Targeted Therapies for Triple-Negative Breast Cancer. Curr Treat Options Oncol 2019;20:82. [Crossref] [PubMed]
- Bianchini G, Balko JM, Mayer IA, et al. Triple-negative breast cancer: challenges and opportunities of a heterogeneous disease. Nat Rev Clin Oncol 2016;13:674-90. [Crossref] [PubMed]
- Tomioka N, Azuma M, Ikarashi M, et al. The therapeutic candidate for immune checkpoint inhibitors elucidated by the status of tumor-infiltrating lymphocytes (TILs) and programmed death ligand 1 (PD-L1) expression in triple negative breast cancer (TNBC). Breast Cancer 2018;25:34-42. [Crossref] [PubMed]
- Cimino-Mathews A, Thompson E, Taube JM, et al. PD-L1 (B7-H1) expression and the immune tumor microenvironment in primary and metastatic breast carcinomas. Hum Pathol 2016;47:52-63. [Crossref] [PubMed]
- Schmid P, Adams S, Rugo HS, et al. Atezolizumab and Nab-Paclitaxel in Advanced Triple-Negative Breast Cancer. N Engl J Med 2018;379:2108-21. [Crossref] [PubMed]
- Adams S, Schmid P, Rugo HS, et al. Phase 2 study of pembrolizumab (pembro) monotherapy for previously treated metastatic triple-negative breast cancer (mTNBC): KEYNOTE-086 cohort A. J Clin Oncol 2017;35:abstr 1008.
- Sharma P, Allison JP. The future of immune checkpoint therapy. Science 2015;348:56-61. [Crossref] [PubMed]
- Ribas A, Wolchok JD. Cancer immunotherapy using checkpoint blockade. Science 2018;359:1350-5. [Crossref] [PubMed]
- Koh YW, Han JH, Haam S, et al. Increased CMTM6 can predict the clinical response to PD-1 inhibitors in non-small cell lung cancer patients. Oncoimmunology 2019;8:e1629261. [Crossref] [PubMed]
- Nakamura S, Hayashi K, Imaoka Y, et al. Intratumoral heterogeneity of programmed cell death ligand-1 expression is common in lung cancer. PLoS One 2017;12:e0186192. [Crossref] [PubMed]
- Mamessier E, Birnbaum DJ, Finetti P, et al. CMTM6 stabilizes PD-L1 expression and refines its prognostic value in tumors. Ann Transl Med 2018;6:54. [Crossref] [PubMed]
- Zhu X, Qi G, Li C, et al. Expression and Clinical Significance of CMTM6 in Hepatocellular Carcinoma. DNA Cell Biol 2019;38:193-7. [Crossref] [PubMed]
- Guan X, Zhang C, Zhao J, et al. CMTM6 overexpression is associated with molecular and clinical characteristics of malignancy and predicts poor prognosis in gliomas. EBioMedicine 2018;35:233-43. [Crossref] [PubMed]
- Li X, Chen L, Gu C, et al. CMTM6 significantly relates to PD-L1 and predicts the prognosis of gastric cancer patients. PeerJ 2020;8:e9536. [Crossref] [PubMed]
(English Language Editor: D. Fitzgerald)