- Open Access
Connecting METTL3 and intratumoural CD33+ MDSCs in predicting clinical outcome in cervical cancer
Journal of Translational Medicine volume 18, Article number: 393 (2020)
Methyltransferase-like 3 (METTL3) is a member of the m6A methyltransferase family and acts as an oncogene in cancers. Recent studies suggest that host innate immunity is regulated by the enzymes controlling m6A epitranscriptomic changes. Here, we aim to explore the associations between the levels of METTL3 and CD33+ myeloid-derived suppressor cells (MDSCs) in tumour tissues and the survival of patients with cervical cancer (CC).
Specimens of paraffin embedded tumour from 197 CC patients were collected. The expression levels of METTL3 and CD33 were measured by immunohistochemical (IHC) staining. The clinical associations of the IHC variants were analysed by Pearson’s or Spearman’s chi-square tests. Overall survival (OS) and disease-free survival (DFS) were estimated by the Kaplan–Meier method and log-rank test. Hazard ratios (HRs) and independent significance were obtained via Cox proportional hazards models for multivariate analyses. METTL3 in CD33+ cells or CC-derived cells was knocked down by METTL3-specific siRNA, and MDSC induction in vitro was performed in a co-culture system in the presence of METTL3-siRNA and METTL3-knockdown-CC-derived cells compared with that of the corresponding controls.
We found that tumour tissues displayed increased levels of METTL3 and CD33+ MDSCs compared with tumour-adjacent tissues from the same CC patients. Importantly, METTL3 expression was positively related to the density of CD33+ cells in tumour tissues (P = 0.011). We further found that the direct CD33+CD11b+HLA-DR− MDSC induction and tumour-derived MDSC induction in vitro were decreased in the absence of METTL3. The level of METTL3 in tumour microenvironments was significantly related to advanced tumour stage. The levels of METTL3 and CD33+ MDSCs in tumour tissues were notably associated with reduced DFS or OS. Cox model analysis revealed that the level of METTL3 in tumour cells was an independent factor for patient survival, specifically for DFS (HR = 3.157, P = 0.022) and OS (HR = 3.271, P = 0.012), while the CD33+ MDSC number was an independent predictor for DFS (HR: 3.958, P = 0.031). Interestingly, in patients with advanced-disease stages (II–IV), METTL3 in tumour cells was an independent factor for DFS (HR = 6.725, P = 0.010) and OS (HR = 5.140, P = 0.021), while CD33+ MDSC density was an independent factor for OS (HR = 8.802, P = 0.037).
Our findings suggest that CD33+ MDSC expansion is linked to high levels of METTL3 and that METTL3 and CD33+ MDSCs are independent prognostic factors in CC.
Cervical cancer (CC) is one of the most common tumours, ranking fourth for both incidence and mortality in women worldwide [1,2,3]. CC is the result of continuous infection with some strains of human papillomavirus (HPV), such as HPV16 and HPV18 [4, 5]. Though there are abundant measures of prevention and cure, cervical cancer continues to exhibit high invasion and poor prognosis . In the past decade, researchers worldwide have found that the expression levels of molecular markers in the tumour microenvironment could be an essential factor for cervical cancer (CC) growth and metastasis [7, 8]. In addition to traditional prognostic factors, including age, WHO grade, TNM stage and clinical status, some of the molecular markers could be new predictors of CC prognosis [6, 9, 10]. However, there are no confirmed molecular markers for tumour progression or prognosis in CC patients. The related viral proteins E6 and E7 have been the focal points of research for the past several years [11, 12]. In other words, easily detectable and meaningful molecular markers need to be confirmed.
Methyltransferase-like 3 (METTL3) is associated with N6-methyladenosine (m6A) RNA methylation, which is the most abundant modification ubiquitously occurring in eukaryotic mRNAs [13, 14]. This modification regulates mRNA stability or translation and can affect many functions, such as immune cell differentiation, cell development, circadian periods and tumour growth [15, 16]. In previous studies, METTL3 was found to have an adverse influence on acute myeloblastic leukaemia (AML), breast cancer (BC), ovarian carcinoma, bladder cancer (BC) and gastric cancer (GC) [17,18,19,20,21,22]. Additionally, m6A modifications are carried out by a combination of m6A methyltransferases (also named writers: METTL3, METTL14 and WTAP), m6A demethylases (also named erasers: FTO and ALKBH5) and specific RNA-binding proteins (also named readers: YTHDF1/2/3, HNRNPA2B1, IGF2BP1/2/3, eIF3 [22,23,24].
CD33-positive cells are usually defined as myeloid-derived suppressor cells (MDSCs) with suppressive influence on human tumour tissues [25, 26]. MDSCs in the tumour environment were confirmed to be an independent indicator of poor prognosis in patients with many solid tumours [25, 27, 28]. In our previous studies, the MDSC proportion was expanded in the tumour microenvironment and showed extensive negative regulatory function for antitumour immunity in malignancies [29,30,31]. Recent studies have indicated that the differentiation of myeloid cells is regulated by m6A methyltransferases, including METTL3 [22, 32, 33]. We hypothesized that MDSC expansion may be linked to the level of METTL3 in the microenvironment of CC.
In the present study, we detected the levels of METTL3 and CD33+ MDSCs in tumour specimens from 197 CC patients by immunohistochemical (IHC) staining. We observed increased levels of METTL3 and CD33+ MDSCs in tumour tissues and positive associations between the levels of METTL3 and CD33+ MDSCs. The high levels of METTL3 and CD33+ MDSCs in the CC tumour microenvironment were significantly associated with poor disease-free survival (DFS) and overall survival (OS) in CC patients. Importantly, METTL3 and CD33+ MDSCs were independent prognostic predictors for CC patients. These findings suggest that METTL3 and MDSCs contribute to the development of disease and that METTL3 may respond to MDSC expansion in tumour microenvironments in CC.
Patients and tissue samples
A total of one hundred ninety-seven CC patients who received therapy at Sun Yat-Sen University Cancer Center in Guangzhou, China, and who accepted medical follow-up that continued until 2019 were included. Paraffin tumour specimens from 197 CC patients were collected at Sun Yat-Sen University Cancer Center between 2008 and 2010. In this retrospective study, none of the patients received antitumour treatment before tumour tissue was obtained, and all 197 patients were histologically confirmed as having primary CC. As shown in Table 1, according to the WHO and the International Federation of Gynecology and Obstetrics classification criteria, 127 (64.5%) patients had stage I disease, and 70 (35.5%) had stage II–IV disease; 26 (13.2%) died. Forty-five (22.8%) patients received surgery alone, 41 (20.8%) patients received radiation therapy alone, 27 (13.7%) patients received surgery + radiation therapy, and 84 (42.6%) patients received surgery + radiation therapy + chemotherapy. The tumour specimens and clinical information were provided by the Pathology Department of Sun Yat-Sen University Cancer Center. The study was approved by the Research Ethics Committee of the Sun Yat-sen University Cancer Center, and written informed consent was obtained from all 197 patients.
The HeLa cervical cancer cell line is maintained in our laboratory. HeLa cells were cultured RPMI 1640 (Gibco, Grand Island, NY, USA) containing 10% foetal bovine serum (FBS; EXCELL BIO, Florida, USA).
Antibodies and reagents
The antibodies and reagents used in this study were as follows. For immunoblotting, immunohistochemistry (IHC) and immunofluorescence staining, rabbit anti-METTL3 antibody (ab195352; Abcam, Cambridge, UK), rabbit anti-CD33 antibody (ab199432; Abcam, Cambridge, UK), rabbit mAb IgG control (ab172730; Abcam, Cambridge, UK), and rabbit anti-GAPDH antibody (10494; Proteintech, Wuhan, China)) were used. Rabbit anti-human IgG (H + L) secondary antibody, HRP (PV-6001-6.0, ZSGB Bio, Beijing, China), a DAB Horseradish Peroxidase Colour Development Kit (ZLI-9017, ZSGB Bio, Beijing, China), and RIPA cell lysis buffer (CW2333S CWBIO, Beijing, China) were also applied. Milk (1706404, Bio-Rad, Hercules, USA) and West dura extended duration substrate (34075, Thermo Scientific, Carlsbad, USA), haematoxylin (DH0001, Solarbio, Beijing, China) and DAPI (GB1012, Servicebio, Wuhan, China); Lipofectamine™ LTX Reagent with PLUS™ Reagent (15338100, Invitrogen, Carlsbad, USA) were used for the siRNA assay. For CD33+ cell isolation, human CD33-antibody-linked magnetic beads (23227, Invitrogen, Carlsbad, USA) were applied. Antibodies for flow cytometry were as follows: PerCP-Cyanine5.5 anti-human CD33 antibody (45-0338-42, eBioscience, San Diego, USA), APC anti-human HLA-DR antibody (17-9956-42, eBioscience, San Diego, USA), and PE anti-human CD11b antibody (12-0118-42, eBioscience, San Diego, USA). Fixable Viability Stain 700 (564997, BD, San Jose, USA) was used to distinguish the dead cells.
Immunohistochemistry (IHC) and immunofluorescence staining
Paraffin-embedded tissues were continuously sectioned at a thickness of 4 μm, and the immunohistochemistry kit was used according to the manufacturer’s instructions. In brief, tissue sections were deparaffinized by xylene, rehydrated in graded alcohols and immersed in EDTA (PH 8.0). Microwave (95 °C 12 min) was applied for antigen retrieval, and samples were cooled to room temperature. The endogenous enzyme block reagent was used to block the activities of endogenous peroxidase. The goat serum was applied to block nonspecific binding sites at room temperature for 30 min. Primary antibodies, including rabbit anti-METTL3 antibody (1:400), rabbit anti-CD33 antibody (1:200), and rabbit mAb IgG control (1:200) were incubated at 37 °C for 1 h and developed with peroxidase. After staining by haematoxylin, images were taken under a microscope (NIKON ECLIPSE 80i). The expression of METTL3 on CD33+ cells was measured by immunofluorescence staining; DAPI was used to stain the nuclei. The images were taken with a fluorescence microscope (NIKON ECLIPSE C1).
The METTL3 expression level was scored in tumour cells in five to ten separate × 400 high-power fields (HPFs). We scored METTL3 expression in the tumour cells of each specimen using a semiquantitative immunoreactivity scoring system, which ranged from 0 to 12 and was equal to multiplication of the intensity of immunohistochemical staining (zero: no staining; one: weak staining; two: moderate staining; and three: strong staining) and the percentage of positive tumour cells (one: less than 25%; two: 25–50%; three: 50–75%; and four: more than 75%). The expression of CD33 was determined by counting CD33-positive cells from five to ten separate × 400 HPFs from the same patient. METTL3 expression level in tumour-infiltrating cells (TILs) was evaluated based on the mean percentage from five to ten separate × 400 HPFs from the same patient. These METTL3- and CD33-positive scores were determined separately by two pathologists. An isotype control IgG antibody was included.
Knockdown of METTL3 by siRNA
To knock down METTL3 in HeLa cells or CD33+ cells, we generated METTL3-specific siRNA (siMETTL3) with the help of RiboBio; a control-siRNA vector was also generated. The siRNAs were transiently transfected into HeLa or CD33+ cells by Lipofectamine™ LTX Reagent with PLUS™ Reagent according to the manufacturer's instructions. After 48 h, the cells were harvested for immunoblotting and MDSC induction.
The sequences of siMETTL3 siRNAs were as follows:
siMETTL3_002 5′-GACTGCTCTTTCCTTAATA -3′; and.
MDSCs induction in vitro and fluorescence-activated cell sort (FACS) staining
Peripheral blood mononuclear cells (PBMCs) were derived from the peripheral blood of healthy donors by gradient centrifugation separation. The CD33+ cells were sorted by human CD33-antibody-linked magnetic beads. After isolation, the 1 × 10^6 CD33+ cells were seeded in a 24-well plate (outer well) and co-cultured with HeLa cells (inner well, at 1:2 ratio) with or without METTL3 knockdown in a Transwell system (3421, Corning, New York, NY, USA) for 48 h. The cells were harvested for FACS staining and detected by cytometry CytExpert software (Beckman Coulter, San Jose, USA) or immunoblotting. Cells were pipetted into single-cell suspensions and incubated with corresponding fluorescence-labelled antibodies according to the manufacturer’s instructions. The flow cytometer used was a cytoFLEX (Beckman), and all data were analysed by the original analysis software provided with the flow cytometer (CytExpert). The CD33+CD11b+HLA-DR− cells were defined as peripheral MDSCs in this study.
The harvested cells were lysed with pre-cooled RIPA buffer, and the proteins were quantified by a BCA protein assay kit (23227, Invitrogen, Carlsbad, USA) and then separated using a 10% SDS-PAGE. Proteins were transferred onto polyvinylidene difluoride membrane (IPVH00010; Millipore, Massachusetts, USA). The membrane was blocked with 5% milk and incubated with the corresponding primary antibodies at 4 °C overnight. Next, the membrane was incubated with HRP-coupled secondary antibodies at room temperature and detected using a West dura extended duration substrate.
SPSS 19.0 software (SPSS Inc., Chicago, USA) was used to analyse all the data, and GraphPad Prism 7 software (La Jolla, USA) was used to obtain the curves. The median values were used as cut-off values to divide the patients into two groups (high level and low level). We used Pearson’s chi-square test or Spearman's chi-square test to analyse the relationships between immunohistochemical variants in different cell populations and patients’ clinical parameters. The relationships among the expression of METTL3 in tumour cells, METTL3 in tumour-infiltrating immune cells and CD33 in tumour-infiltrating immune cells were determined using Pearson's or Spearman’s correlation coefficient and linear regression analyses. Cut-off selection was based on X-tile (Version 3.6.1, New Haven, USA). Then, we evaluated prognostic factors in univariate and multivariate analyses using the Cox proportional hazards model. In our research, *P < 0.05 was considered significant. Raw data of this article have been deposited in the Research Data Deposit (RDD) (www.researchdata.org.cn) with accession number RDDB2020000943.
The level of METTL3 is positively linked to the number of CD33+ MDSCs and contributes to tumour development
In the present study, the levels of METTL3 and CD33+ MDSCs were examined in tumour specimens from 197 patients with CC by IHC. METTL3 was located in the nuclei of tumour cells and tumour-infiltrating immune cells, while CD33+ cells were scattered mainly in the tumour stroma (Fig. 1a, b); isotype IgG was used as a control (Fig. 1c). Importantly, we found that CD33 and METTL3 co-localized in some tumour-infiltrated immune cells (Fig. 1d). We further demonstrated that the percentage of CD33+CD11b+HLA-DR− peripheral MDSCs was increased in CC patients compared with healthy donors, as was the percentage of tumour-derived CD33+CD11b+HLA-DR− MDSCs compared with that of tumour-adjacent tissues (Fig. 1e, f, n = 3). Consistent with the increase in the MDSC population in CC patients, the level of METTL3 was increased in the peripheral and tumour-infiltrating immune cells compared with the corresponding controls (Fig. 1g, h). Among the 197 patients with CC, the median survival time was 96 months (range: 0–120 months), and the 10-year DFS and 10-year OS rates were 88.83 and 86.80%, respectively (Additional file 1: Figure S1A and S1B). Table 2 shows the results of the relationships between clinicopathological features and immunohistochemical variants in different cell types in the tumour microenvironment. High METTL3 expression in the tumour and in tumour-infiltrating immune cells was linked to tumour stage (P = 0.040 and 0.020, respectively).
In addition, we analysed the relationship between METTL3 expression in tumour cells and in tumour-infiltrating immune cells and the number of CD33+ MDSCs via Spearman’s correlation coefficient and linear regression. The expression of METTL3 in tumour cells was positively correlated with that in tumour-infiltrating immune cells (R = 0.264, P < 0.001) (Fig. 2a). The number of CD33+ cells was positively correlated with the expression of METTL3 in tumour cells (R = 0.145, P = 0.041) and tumour-infiltrating immune cells (R = 0.182, P = 0.011) (Fig. 2b, c). Moreover, we found that the METTL3 level in tumour cells was positively correlated with TILs in the early (R = 0.049, P = 0.012) and advanced stage (R = 0.129, P = 0.002) (Additional file 1: Figure S1E and S1F), while we found that in the advanced stage, the number of CD33+ cells was positively correlated with the METTL3 level in TILs (R = 0.088, P = 0.013) (Additional file 1: Figure S1J).
To further investigate the role of METTL3 in the regulation of MDSC expansion, we knocked down METL3 expression in CD33+ cells or HeLa cells. We found that CD33+CD11b+HLA-DR− MDSCs and tumour-derived MDSCs were decreased when METTL3 was knocked down in CD33+ cells or HeLa cells (Fig. 2d–g).
High METTL3 levels and CD33+ MDSC density are associated with poor outcomes
To evaluate the expression of METTL3 and CD33 as predictors for the prognosis of the 197 patients, Kaplan–Meier survival curves were used for analysis. The high level of METTL3 in tumour cells was significantly correlated with decreased DFS (P < 0.001, Fig. 3a) and OS (P < 0.001, Fig. 3b) in CC patients. Accordingly, a high level of METTL3 in tumour-infiltrating immune cells was negatively correlated with DFS (P = 0.002, Fig. 3c) and OS (P < 0.001, Fig. 3d) in CC patients. The high density of CD33+ MDSCs was obviously correlated with decreased DFS (P < 0.001, Fig. 3e) and OS (P < 0.001, Fig. 3f) in CC patients.
METTL3 and CD33+ MDSCs are independent factors for patient prognosis
As shown in Table 3, univariate analysis showed that in addition to lymph node involvement and clinical stage, high levels of METTL3 in tumour cells (HR: 4.244, P = 0.002) and in tumour-infiltrating immune cells (HR: 4.857, P = 0.004) and a high density of CD33+ MDSCs (HR: 6.579, P = 0.002) were noticeably correlated with reduced DFS. In addition, we found that high levels of METTL3 in tumour cells (HR: 5.502, P = 0.001) and in tumour-infiltrating immune cells (HR: 6.021, P = 0.001) and a high density of CD33+ MDSCs (HR: 5.755, P = 0.001) were associated with decreased OS. Moreover, clinicopathological parameters such as clinical stage (HR: 3.511, P = 0.005) and nodal status (HR: 2.798, P = 0.032) also had prognostic value with decreased DFS, and clinical stage (HR: 3.820, P = 0.001) was related to decreased OS. In addition, other clinical characteristics, such as age and tumour status, were not clearly related to DFS and OS (Table 3). When we performed multivariate Cox proportional hazards regression analysis in Table 4, we included all the significant univariate variables. For all 197 patients, in addition to clinical stage (HR: 3.827, P = 0.003), N status (HR: 3.219, P = 0.021) was an independent factor for DFS, clinical stage (HR: 4.248, P < 0.001) was an independent factor for OS, and METTL3 levels in tumour cells (HR: 3.157, P = 0.022) and in tumour-infiltrating immune cells (HR: 3.368, P = 0.036) and CD33+ MDSCs (HR: 3.958, P = 0.031) were independent factors for both DFS and OS (Table 4).
METTL3 and CD33+ MDSCs have predictive value for patients with early and advanced disease stages
We further divided the 197 patients into two subgroups based on the clinicopathological stage: 127 of the total patients were in early disease stage (stage I), while 70 of the total patients were in advanced disease stage (stage II–IV). Through the Kaplan–Meier method, we found that the high expression of METTL3 in tumour-infiltrating immune cells was significantly correlated with poor DFS (P = 0.033) and OS (P = 0.019) (Additional file 2: Figure S2C and S2D) in patients with early disease stage, while there was no significant association between the high expression of METTL3 in tumour cells (P = 0.400 vs P = 0.183) and the number of CD33+ MDSCs (P = 0.393 vs P = 0.227) with the DFS and OS of patients with early-stage disease (Additional file 2: Figure S2A, S2B, S2E and S2F). For patients with advanced-stage disease (n = 70), a high level of METTL3 in tumour cells was dramatically correlated with decreased DFS (P < 0.001, Fig. 4a) and OS (P < 0.001, Fig. 4b), and a high level of METTL3 in tumour-infiltrating immune cells was negatively correlated with DFS (P = 0.004, Fig. 4c) and OS (P < 0.001, Fig. 4d); the increased number of CD33+ MDSCs was dramatically correlated with poor DFS (P < 0.001, Fig. 4e) and OS (P < 0.001, Fig. 4f). Using multivariate Cox regression analysis in the 70 patients with advanced-stage disease, METTL3 expression in tumour cells (HR: 6.725, P = 0.010) was an independent prognostic factor for DFS, while METTL3 expression in tumour cells (HR: 5.140, P = 0.021) and CD33+ MDSCs (HR: 8.802, P = 0.037) were independent prognostic factors for OS (Table 4).
The combination of METTL3 levels and CD33+ MDSCs was associated with the survival of patients with CC
Finally, considering that METTL3 levels were positively correlated with high CD33+ MDSC infiltration, we calculated the significance of the combination of these two biomarkers for the survival of CC patients. All 197 patients were divided into three groups. Patients with low levels of both METTL3 in tumour-infiltrating immune cells and CD33+ MDSCs were included in the combined low expression group, those with high levels of only one of the two biomarkers were included in the combined medium expression group, and those with high levels of both were included in the combined high expression group. The high combination of METTL3 and intratumoural CD33+ MDSCs was associated with reduced DFS (P < 0.001, Fig. 5a) and OS (P < 0.001, Fig. 5b). In the patients (127) with early-stage disease, the high combination of METTL3 and CD33+ MDSCs was not related to DFS (P = 0.063, Fig. 5c) but was clearly negatively related to OS (P = 0.037, Fig. 5d). In the patients (70) with advanced-stage disease, the combination of high METTL3 levels and CD33+ MDSCs was clearly related to unfavourable DFS (P < 0.001, Fig. 5e) and OS (P < 0.001, Fig. 5f).
Indeed, compared to METTL3 or CD33+ MDSCs, the combination of METTL3 and CD33+ MDSCs can improve the prognostic stratification of survival for CC patients, especially those in advanced disease stages. Based on a total of 197 patients, the combination of high METTL3 levels and CD33+ MDSCs was a predictor of worse patient prognosis, including DFS [HR (95% CI): 4.672 (2.149–10.156), P < 0.001] and OS [HR (95% CI): 4.890 (2.369–10.093), P < 0.001]. In the patients with early-stage disease, we found that the combination of high METTL3 levels and intratumoural CD33+ MDSCs was negatively correlated with OS [HR (95% CI): 3.071 (1.056–8.931), P = 0.039], but there was no significant association with DFS. Importantly, we found that the combination of high tumour METTL3 levels and intratumoural CD33+ MDSCs was significantly correlated with poor DFS [HR (95% CI): 7.673 (2.420–24.324), P = 0.001] and OS [HR (95% CI): 7.286 (2.667–19.902), P < 0.001] in patients with advanced disease stages (Table 3), suggesting that the combination of high METTL3 levels and CD33+ MDSCs improved patient prognostic stratification in those with advanced disease.
The development of tumour cells depends on the tumour microenvironment, which includes tumour cells, various other cells and extracellular components . The immunosuppressive cells in the tumour microenvironment, such as Tregs and MDSCs, not only affect each other, but their changes in number and types will affect tumour development [34, 35]. METTL3 is one of the ‘writers’, and its role is to catalyse the m6A methylation of mRNA (and other nuclear RNAs); after the methylation of m6A, RNAs will nucleate and transport to the cytoplasm faster and then produce more proteins for function and proliferation. Some studies have shown that METTL3 expression can promote tumour cell proliferation, leading to poor patient prognosis. The tumour-infiltrated MDSC population usually induces antitumour immunity tolerance by inhibiting the proliferation and function of T cells, such as hindering antigen presentation by antigen-presenting cells . Increased METTL3 levels and CD33+ MDSCs have been found in tumour microenvironments and lead to a poor prognosis [37,38,39,40]. In this study, we focused on the distribution of METTL3 and CD33+ MDSCs in the tumour microenvironment of 197 patients with CC. The positive association between METTL3 levels and CD33+ MDSCs and the prognostic value of these two variants in CC patients were demonstrated. Importantly, we demonstrated that knockdown of METTL3 in CD33+ cells or HeLa cells could attenuate MDSC or tumour-associated MDSC differentiation in vitro.
M6A methyltransferases, especially METTL3, can affect many physiological and pathological diseases through p53 and other genes . At the nucleic acid level, silencing m6A methyltransferase significantly affects gene expression and mRNA splicing patterns, leading to changes in normal cell signalling pathways and apoptosis . In bladder cancer cells, m6A-modified direct targets IKBKB and RELA (two key regulators of the NF-κB pathway) mediated by METTL3 become factors that promote tumour development ; in glioblastoma stem cells (GSCs), knocking down METTL3 can induce changes in m6A-enriched mRNA and alter the mRNA expression of genes with key biological functions in GSCs (such as ADAM1937) . In recent studies, high METTL3 levels were related to tumour invasion and poor outcomes in breast cancer and acute myeloid leukaemia (AML) [21, 42]. Our results are consistent with the results of these studies, showing that high METTL3 expression results in poor prognosis in CC patients. METTL3 regulates haematopoietic stem cell differentiation and induces the development of leukaemic cells by upregulating MYC expression [42, 43]. Therefore, we wondered whether METTL3 expression may be linked to the density of tumour-infiltrated MDSCs. Our data identified a positive association between METTL3 expression in tumour cells and in tumour-infiltrating immune cells and intratumoural CD33+ MDSC density. The results indicate that METTL3 could directly induce CD33+CD11b + HLA-DR− MDSC differentiation or tumour-associated MDSC differentiation in vitro. Moreover, both METTL3 and CD33+ MDSCs were independent factors for the prognosis of CC patients, and the combination of METTL3 levels and CD33+ MDSC density displayed prognostic value for CC patients, including patients at early or late disease stages. The function, distribution and clinical relevance of the proportion of tumour-derived CD33+ MDSCs have been explored in recent years. MDSCs are generally elevated in tumour tissues and in the peripheral blood of cancer patients and are linked to antitumour immunity suppression, resulting in tumour growth and metastasis [25, 34, 44]. In our study, CC patients with a high infiltration of MDSCs in the cervical cancer microenvironment showed a poor prognosis, which is consistent with observations in other solid cancers. The tumour microenvironment is a main battleground between tumour cells and the host immune system. Tumour cells usually ‘educate’ infiltrated immune cells through many factors, such as cytokines or tumour-derived exosomes, to affect the proliferation, differentiation and function of tumour-infiltrating immune cells, resulting in the expansion of suppressive immune cells, including M2 macrophages, MDSCs and Tregs, and limiting the antitumour effect of cytotoxic T cells. Epigenetic modifications such as RNA modification, DNA methylation and histone modifications can rapidly regulate infiltrated immune cell differentiation and activities in tumour microenvironments . Here, our data suggest that METTL3-mediated m6A RNA modification is positively associated with the increase in MDSC expansion and affects tumour development and prognosis in CC and induces CD33 + cells to differentiate into MDSCs in the tumour microenvironment. We further demonstrate the prognostic value of the combination of the METTL3 level in tumour-infiltrating immune cells and CD33+ MDSC density in CC patients, especially for those in advanced disease stages. A mechanistic study to support the role of METTL3 in the regulation of tumour-derived MDSC differentiation is currently underway, and the underlying mechanisms will be clarified in the near future.
Our study demonstrated a comprehensive result of the relationship between METTL3 and CD33 in CC and revealed that METTL3 could induce direct MDSC and tumour-associated MDSC differentiation in vitro. The results showed that both biomarkers were adverse indicators for prognosis and may have significant relationships in the microenvironment of CC. Our research may offer clues for further research into the mechanism behind METTL3 in the regulation of MDSC-mediated immune suppression in the CC microenvironment.
Availability of data and materials
The authenticity of this article has been validated by uploading the key raw data to the Research Data Deposit (RDD) public platform (www. researchdata.org.cn).
Myeloid-derived suppressor cells
Tumour node metastasis
World Health Organization
- T stage:
- N status:
Lymph node metastasis
Peripheral blood mononuclear cells
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Mayadev JS, Enserro D, Lin YG, Da Silva DM, Lankes HA, Aghajanian C, et al. Sequential ipilimumab after chemoradiotherapy in curative-intent treatment of patients with node-positive cervical cancer. JAMA Oncol. 2019;6:92.
Morel A, Neuzillet C, Wack M, Lameiras S, Vacher S, Deloger M, et al. Mechanistic signatures of human papillomavirus insertions in anal squamous cell carcinomas. Cancers. 2019;11(12):1846.
Cohen AC, Roane BM, Leath CA 3rd. Novel Therapeutics for Recurrent Cervical Cancer: Moving Towards Personalized Therapy. Drugs. 2020;80:217.
Luo X, Donnelly CR, Gong W, Heath BR, Hao Y, Donnelly LA, et al. HPV16 drives cancer immune escape via NLRX1-mediated degradation of STING. J Clin Invest. 2019;130:1635.
Burger EA, de Kok I, Groene E, Killen J, Canfell K, Kulasingam S, et al. Estimating the Natural History of cervical carcinogenesis using simulation models: a CISNET comparative analysis. J Natl Cancer Inst. 2019;112:955.
Wu T, Dai Y. Tumour microenvironment and therapeutic response. Cancer Lett. 2017;387:61–8.
Nicoli F, Mantelli B, Gallerani E, Telatin V, Bonazzi I, Marconi P, et al. HPV-specific systemic antibody responses and memory B cells are independently maintained up to 6 years and in a vaccine-specific manner following immunization with Cervarix and Gardasil in adolescent and young adult women in vaccination programs in Italy. Vaccines. 2020;8(1):26.
Babion I, Jaspers A, van Splunter AP, van der Hoorn IAE, Wilting SM, Steenbergen RDM. miR-9–5p exerts a dual role in cervical cancer and targets transcription factor TWIST1. Cells. 2019;9(1):65.
Small W Jr, Bacon MA, Bajaj A, Chuang LT, Fisher BJ, Harkenrider MM, et al. Cervical cancer: a global health crisis. Cancer. 2017;123(13):2404–12.
Choi CH, Choi HJ, Lee JW, Kang ES, Cho D, Park BK, et al. Phase I study of a B cell-based and monocyte-based immunotherapeutic vaccine, BVAC-C in human papillomavirus type 16- or 18-positive recurrent cervical cancer. J Clin Med. 2020;9(1):147.
Bergman H, Buckley BS, Villanueva G, Petkovic J, Garritty C, Lutje V, et al. Comparison of different human papillomavirus (HPV) vaccine types and dose schedules for prevention of HPV-related disease in females and males. Cochrane Database Syst Rev. 2019;2019(11):CD013479.
Cheng M, Sheng L, Gao Q, Xiong Q, Zhang H, Wu M, et al. The m(6)A methyltransferase METTL3 promotes bladder cancer progression via AFF4/NF-kappaB/MYC signaling network. Oncogene. 2019;38(19):3667–800.
Deng X, Su R, Weng H, Huang H, Li Z, Chen J. RNA N(6)-methyladenosine modification in cancers: current status and perspectives. Cell Res. 2018;28(5):507–17.
Lan Q, Liu PY, Haase J, Bell JL, Huttelmaier S, Liu T. The Critical role of RNA m(6)A methylation in cancer. Can Res. 2019;79(7):1285–92.
Li X, Tang J, Huang W, Wang F, Li P, Qin C, et al. The M6A methyltransferase METTL3: acting as a tumour suppressor in renal cell carcinoma. Oncotarget. 2017;8(56):96103–16.
Feng Z, Li Q, Meng R, Yi B, Xu Q. METTL3 regulates alternative splicing of MyD88 upon the lipopolysaccharide-induced inflammatory response in human dental pulp cells. J Cell Mol Med. 2018;22(5):2558–68.
Hua W, Zhao Y, Jin X, Yu D, He J, Xie D, et al. METTL3 promotes ovarian carcinoma growth and invasion through the regulation of AXL translation and epithelial to mesenchymal transition. Gynecol Oncol. 2018;151(2):356–65.
Lin S, Liu J, Jiang W, Wang P, Sun C, Wang X, et al. METTL3 promotes the proliferation and mobility of gastric cancer cells. Open Med. 2019;14:25–31.
Lin S, Choe J, Du P, Triboulet R, Gregory RI. The m(6)A methyltransferase METTL3 promotes translation in human cancer cells. Mol Cell. 2016;62(3):335–45.
Cai X, Wang X, Cao C, Gao Y, Zhang S, Yang Z, et al. HBXIP-elevated methyltransferase METTL3 promotes the progression of breast cancer via inhibiting tumour suppressor let-7g. Cancer Lett. 2018;415:11–9.
Chen XY, Zhang J, Zhu JS. The role of m(6)A RNA methylation in human cancer. Mol Cancer. 2019;18(1):103.
Du Y, Hou G, Zhang H, Dou J, He J, Guo Y, et al. SUMOylation of the m6A-RNA methyltransferase METTL3 modulates its function. Nucleic Acids Res. 2018;46(10):5195–208.
Geula S, Moshitch-Moshkovitz S, Dominissini D, Mansour AA, Kol N, Salmon-Divon M, et al. Stem cells. m6A mRNA methylation facilitates resolution of naive pluripotency toward differentiation. Science. 2015;347(6225):1002–6.
Veglia F, Perego M, Gabrilovich D. Myeloid-derived suppressor cells coming of age. Nat Immunol. 2018;19(2):108–19.
Vasquez-Dunddel D, Pan F, Zeng Q, Gorbounov M, Albesiano E, Fu J, et al. STAT3 regulates arginase-I in myeloid-derived suppressor cells from cancer patients. J Clin Invest. 2013;123(4):1580–9.
Yang Z, Guo J, Weng L, Tang W, Jin S, Ma W. Myeloid-derived suppressor cells-new and exciting players in lung cancer. J Hematol Oncol. 2020;13(1):10.
Haas L, Obenauf AC. Allies or enemies-the multifaceted role of myeloid cells in the tumour microenvironment. Front Immunol. 2019;10:2746.
Cai TT, Ye SB, Liu YN, He J, Chen QY, Mai HQ, et al. LMP1-mediated glycolysis induces myeloid-derived suppressor cell expansion in nasopharyngeal carcinoma. PLoS Pathog. 2017;13(7):e1006503.
Li ZL, Ye SB, OuYang LY, Zhang H, Chen YS, He J, et al. COX-2 promotes metastasis in nasopharyngeal carcinoma by mediating interactions between cancer cells and myeloid-derived suppressor cells. Oncoimmunology. 2015;4(11):e1044712.
Zhang H, Li ZL, Ye SB, Ouyang LY, Chen YS, He J, et al. Myeloid-derived suppressor cells inhibit T cell proliferation in human extranodal NK/T cell lymphoma: a novel prognostic indicator. Cancer Immunol Immunother. 2015;64(12):1587–99.
The RNA methyltransferase METTL3 promotes oncogene translation. Cancer Discov. 2016;6(6):572.
Tong J, Flavell RA, Li HB. RNA m(6)A modification and its function in diseases. Front Med. 2018;12(4):481–9.
Wu L, Liu H, Guo H, Wu Q, Yu S, Qin Y, et al. Circulating and tumour-infiltrating myeloid-derived suppressor cells in cervical carcinoma patients. Oncol Lett. 2018;15(6):9507–15.
Liu Y, Zheng P. Preserving the CTLA-4 checkpoint for safer and more effective cancer immunotherapy. Trends Pharmacol Sci. 2020;41(1):4–12.
Tan Z, Liu L, Chiu MS, Cheung KW, Yan CW, Yu Z, et al. Virotherapy-recruited PMN-MDSC infiltration of mesothelioma blocks antitumour CTL by IL-10-mediated dendritic cell suppression. Oncoimmunology. 2019;8(1):e1518672.
Li Z, Weng H, Su R, Weng X, Zuo Z, Li C, et al. FTO plays an oncogenic role in acute myeloid leukemia as a N(6)-methyladenosine RNA demethylase. Cancer Cell. 2017;31(1):127–41.
Shang W, Gao Y, Tang Z, Zhang Y, Yang R. The Pseudogene Olfr29-ps1 promotes the suppressive function and differentiation of monocytic MDSCs. Cancer Immunol Res. 2019;7(5):813–27.
Su R, Dong L, Li C, Nachtergaele S, Wunderlich M, Qing Y, et al. R-2HG exhibits anti-tumour activity by targeting FTO/m(6)A/MYC/CEBPA signaling. Cell. 2018;172(1–2):90–105.
Vu LP, Pickering BF, Cheng Y, Zaccara S, Nguyen D, Minuesa G, et al. The N(6)-methyladenosine (m(6)A)-forming enzyme METTL3 controls myeloid differentiation of normal hematopoietic and leukemia cells. Nat Med. 2017;23(11):1369–76.
Cui Q, Shi H, Ye P, Li L, Qu Q, Sun G, et al. m(6)A RNA methylation regulates the self-renewal and tumourigenesis of glioblastoma stem cells. Cell Rep. 2017;18(11):2622–34.
Barbieri I, Tzelepis K, Pandolfini L, Shi J, Millan-Zambrano G, Robson SC, et al. Promoter-bound METTL3 maintains myeloid leukaemia by m(6)A-dependent translation control. Nature. 2017;552(7683):126–31.
Lee H, Bao S, Qian Y, Geula S, Leslie J, Zhang C, et al. Stage-specific requirement for Mettl3-dependent m(6)A mRNA methylation during haematopoietic stem cell differentiation. Nat Cell Biol. 2019;21(6):700–9.
Zhang CX, Ye SB, Ni JJ, Cai TT, Liu YN, Huang DJ, et al. STING signaling remodels the tumour microenvironment by antagonizing myeloid-derived suppressor cell expansion. Cell Death Differ. 2019;26(11):2314–28.
Liu M, Zhou J, Chen Z, Cheng AS. Understanding the epigenetic regulation of tumours and their microenvironments. J Pathol. 2017;241:10.
Thanks for the staff in pathology department of Sun Yat-sen University Cancer Center and the donors who provide materials to our research.
This work was supported by the National Natural Science Foundation of China (Grant No. 8197110786, 81773256, 81572982, 81372442 and 81172164); the Sci-Tech Key Program of the Guangzhou City Science Foundation (Grant No. 201802020001); the Guangdong Province Natural Science Foundation [grant no. 2018A030310260].
Ethics approval and consent to participate
Total research procedures were performed with the approval of Research Ethics Committee of the Sun Yat-sen University Cancer Center (GZR2013-040).
Consent for publication
The authors declare no conflicts of interest. All authors read and approved the final manuscript.
The authors have declared that they have no competing interest. The sources that funded this study played no role in the study design, data collection, data analysis, decision to publish, or preparation of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The DFS and OS curves of 197 CC patients, and expression of METTL3 and CD33 in tumours or in tumour-adjacent tissues. (A and B) Disease-free survival (DFS, A) and overall survival (OS, B) curves of 197 CC patients in this study. (C and D) Statistical analysis showing the comparison between the expression of METTL3 (P = 0.999, n = 18) and the number of CD33+ MDSCs (P = 0.295, n = 18) in tumour-adjacent tissues (Non-tumour) and tumour tissues (Tumour). (E and F) The association between METTL3 expression in tumour cells and the expression of METTL3 in TILs in early stage (E, R = 0.049, P = 0.012) or in advanced stage (F, R = 0.129, P = 0.002). (G and H) The association between METTL3 expression in tumour cells and intratumoural CD33+ MDSC number in early stage (G, R = 0.011, P = 0.222) or in advanced stage (H, R = 0.048, P = 0.070). (I and J) The association between METTL3 expression in TILs and intratumoural CD33+ MDSC number in early stage (I, R = 0.013, P = 0.195) or in advanced stage (J, R = 0.088, P = 0.013). Statistics were conducted with a paired Student’s t test in C and D. The correlation statistical analysis was performed using Spearman’s correlation and linear regression. R, Spearman’s correlation, is the correlation coefficient.
Kaplan–Meier curves of DFS and OS according to METTL3 expression in different cell populations and intratumoural CD33+ cells of early-stage patients. (A and B) Kaplan–Meier curves showing the relationship of DFS (P = 0.400, n = 127) and OS (P = 0.183, n = 127) of CC patients and METTL3 expression in tumour cells. (C and D) Kaplan–Meier curves showing the relationship of DFS (P = 0.033, n = 127) and OS (P = 0.019, n = 127) of patients and METTL3 expression in TILs. (E and F) Kaplan–Meier curves showing the relationship of DFS (P = 0.393, n = 127) and OS (P = 0.227, n = 127) of CC patients and intratumoural CD33+ MDSC density. The percentages of DFS and OS were calculated by the Kaplan–Meier method, and P values were calculated by the log-rank statistic. Cut-off selection was based on X-tile.
About this article
Cite this article
Ni, Hh., Zhang, L., Huang, H. et al. Connecting METTL3 and intratumoural CD33+ MDSCs in predicting clinical outcome in cervical cancer. J Transl Med 18, 393 (2020). https://doi.org/10.1186/s12967-020-02553-z
- Tumour microenvironment
- Cervical cancer