Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Abstract 


The m6A mRNA methylation involves in mRNA splicing, degradation and translation. Recent studies have revealed that reduced m6A mRNA methylation might promote cancer development. However, the role of m6A mRNA methylation in cervical cancer development remains unknown. Therefore, we investigated the role of m6A methylation in cervical cancer in the current study. We first evaluated the m6A mRNA methylation level in 286 pairs of cervical cancer samples and their adjacent normal tissues by dot blot assay. Then the role of m6A on patient survival rates and cervical cancer progression were assessed. The m6A level was significantly reduced in the cervical cancer when comparing with the adjacent normal tissue. The m6A level reduction was significantly correlated with the FIGO stage, tumor size, differentiation, lymph invasion and cancer recurrence. It was also shown to be an independent prognostic indicator of disease-free survival and overall survival for patients with cervical cancer. Reducing m6A level via manipulating the m6A regulators expression promoted cervical cancer cell proliferation. And increasing m6A level significantly suppressed tumor development both in vitro and in vivo. Our results showed that the reduced m6A level is tightly associated with cervical cancer development and m6A mRNA methylation might be a potential therapeutic target in cervical cancer.

Free full text 


Logo of oncotargetLink to Publisher's site
Oncotarget. 2017 Nov 17; 8(58): 98918–98930.
Published online 2017 Oct 24. https://doi.org/10.18632/oncotarget.22041
PMCID: PMC5716777
PMID: 29228737

Reduced m6A mRNA methylation is correlated with the progression of human cervical cancer

Abstract

The m6A mRNA methylation involves in mRNA splicing, degradation and translation. Recent studies have revealed that reduced m6A mRNA methylation might promote cancer development. However, the role of m6A mRNA methylation in cervical cancer development remains unknown. Therefore, we investigated the role of m6A methylation in cervical cancer in the current study. We first evaluated the m6A mRNA methylation level in 286 pairs of cervical cancer samples and their adjacent normal tissues by dot blot assay. Then the role of m6A on patient survival rates and cervical cancer progression were assessed. The m6A level was significantly reduced in the cervical cancer when comparing with the adjacent normal tissue. The m6A level reduction was significantly correlated with the FIGO stage, tumor size, differentiation, lymph invasion and cancer recurrence. It was also shown to be an independent prognostic indicator of disease-free survival and overall survival for patients with cervical cancer. Reducing m6A level via manipulating the m6A regulators expression promoted cervical cancer cell proliferation. And increasing m6A level significantly suppressed tumor development both in vitro and in vivo. Our results showed that the reduced m6A level is tightly associated with cervical cancer development and m6A mRNA methylation might be a potential therapeutic target in cervical cancer.

Keywords: m6A RNA methylation, cervical cancer, cancer

INTRODUCTION

Cervical cancer is one of the most prevalent gynecological malignancies in the world [1, 2]. Although major progresses have been made in cancer detection and treatment during the past decades, the 5-year survival rate remains low. Thus, it is important to determine the molecular mechanisms of cervical cancer development, to identify effective prognostic markers and develop novel therapeutic strategies.

Recent studies have uncovered the important role of RNA methylation in cell fate determination, self-renewal and cancer development, indicating a new and promising therapeutic target for investigation [3]. Chemical modified RNA has been demonstrated in recent decades [4]. Among them, the N6-methyladenosine (m6A) modification of mRNA/lncRNA is the most prevalent one and plays an important role in gene expression [5, 6].

The m6A methylation is reversible and regulated by adenosine methyltransferases [79] and demethylases [10, 11]. The m6A methyltransferases include METTL3 (Methyltransferase like 3), METTL14 (Methyltransferase like 14), WTAP (Wilms tumor 1-associated protein), RBM15 (RNA binding motif protein 15) and KIAA1429 [79, 12]. And the m6A demethylases include FTO (Fat mass and obesity-associated protein) and ALKBH5 (AlkB family member 5) [10, 11, 13]. Furthermore, the m6A could be selectively recognized by proteins such as HNRNPC (Heterogeneous nuclear ribonucleoprotein C), HNRNPA2B1 (Heterogeneous nuclear ribonucleoprotein A2/B1), YTHDF2 (YTH N6-methyladenosine RNA binding protein 2), YTHDF1 (YTH N6-methyladenosine RNA binding protein 1) and eIF3 (Eukaryotic initiation factor 3) [5, 1419].

The m6A mRNA methylation involves in mRNA splicing, degradation and translation [11, 14, 1619]. However, its roles in biological processes have just begun to be uncovered [20]. Previous reports have demonstrated that m6A mRNA methylation is implicated in stem cell maintenance and differentiation [2123], and dysregulation of this modification might contribute to development abnormalities, obesity, and other diseases [24, 25]. Recent studies revealed that the m6A mRNA methylation also plays an important role in cancer development, such as breast cancer [26, 27], hematologic malignancies [28, 29] and glioblastoma [30, 31]. However, whether m6A mRNA methylation also plays an important role in cervical cancer development remains unknown. Therefore, we investigated the role of m6A methylation in cervical cancer in the current study.

RESULTS

Decreased m6A level in cervical cancer biopsies

A total of 286 cases with cervical cancer were followed. All these patients had received no pre-operation chemotherapy. They were given the same radical operation and underwent the same adjuvant chemotherapy after the surgery. The m6A RNA methylation level was firstly analyzed and the data showed that the m6A level reduced significantly in cervical cancer tissues comparing with its paired adjacent non-cancerous tissues (Figure (Figure1A1A and and1B).1B). Furthermore, the mRNA level of the m6A methyltransferases (METTL3 and METTL14) were reduced while the demethylases (FTO and ALKBH5) were increased (Supplementary Figure 1). The reduced m6A level was associated with FIGO stage (P = 0.008), tumor size (P = 0.001), differentiation (P = 0.026), lymph invasion (P = 0.009) and recurrence (P = 0.005, Table Table11).

An external file that holds a picture, illustration, etc.
Object name is oncotarget-08-98918-g001.jpg
Reduced m6A is correlated with the poor prognosis of cervical cancer

(A) m6A levels in cancer tissues and adjacent normal tissues were determined by dot blot (n=286). The relative dot density was measured by Image J. * P<0.05. (B) Representative figure for m6A levels in cancer tissues and adjacent normal tissues. N: adjacent normal tissues; C: cancer tissues. (C) Kaplan-Meier survival curve of patients with high or low level of m6A. (D) Kaplan-Meier survival curve of patients with or without mutation of m6A regulators.

Table 1

Association between clinicopathological features and m6A mRNA methylation level
Nm6A mRNA methylation levelP value
Low (n=210)High (n=76)
Age, years 0.913
 <5021415955
 ≥50725121
FIGO stage0.008*
 < IIA1105852
 >IIB17615224
Tumor size, cm0.001*
 <41246460
 ≥416214616
Histological type0.058
 Squamous25819060
 others28208
Differentiation0.026*
 G1 (High)31247
 G2 (Medium)1227844
 G3 (Low)13310825
Lymph invasion0.009*
 Yes887216
 No19813860
Vaginal invasion
 Yes11887310.902
 No16812345
Recurrence0.005*
 Yes32284
 No25418272

*P < 0.05 indicates a significant association among the variables.

During the follow-up for all patients, 52 patients had died and 75 experienced recurrence. Disease-free survival (DFS) and overall survival (OS) was conducted to assess the predictive role of m6A level for distant metastasis. Both DFS and OS were significantly higher in the high m6A level group than the low m6A level group (Figure (Figure1C).1C). The low m6A level group subsequently developed more recurrence than the high m6A level group (P < 0.01, Table Table11).

It has been demonstrated that the genetic alterations of m6A regulators could predict poorer survival in acute myeloid leukemia [29]. Therefore, we also analyzed the gene mutations of m6A regulators, including METTL3, METTL14, WTAP, RBM15, KIAA1429, FTO, ALKBH5, HNRNPC, HNRNPA2B1, YTHDF2, YTHDF1 and eIF3 (Supplementary Figure 2). Mutations in individual genes did not show any significant association with clinicopathological features (data no shown). But when all the mutations were pooled together, the m6A regulators mutations showed significantly association with FIGO stage (P = 0.005), tumor size (P = 0.001), differentiation (P = 0.016), lymph invasion (P = 0.002) and recurrence (P = 0.000, Table Table2).2). Both DFS and OS were significantly higher in the patients without m6A regulators mutations (Figure (Figure1D).1D). Univariate and multivariate analysis showed that patients with reduced m6A level or containing gene mutations of m6A regulators had a significantly reduced OS and DFS (Table (Table3,3, Table Table44).

Table 2

Association between clinicopathological features and m6A regulators mutation
Nm6A regulators mutationP value
No (n=240, %)Yes (n=46, %)
Age, years0.913
 <5021416440
 ≥5072666
FIGO stage0.005*
 ≤IIA1107337
 ≥IIB1761679
Tumor size, cm0.001*
 <41247945
 ≥41621611
Histological type0.168
 Squamous25821345
 others28271
Differentiation0.016*
 G1 (High)31301
 G2 (Medium)1228735
 G3 (Low)13312310
Lymph invasion0.002*
 Yes88871
 No19815345
Vaginal invasion
 Yes118102160.912
 No16813830
Recurrence0.000*
 Yes32320
 No25420846

*P < 0.05 indicates a significant association among the variables.

Table 3

Univariate Cox proportional hazards model for disease-free survival (DFS) and overall survival (OS)
DFSOS
HR95% CIP valueHR95% CIP value
Age, years
 <50
 ≥501.0090.659–1.8410.7140.9340.530–1.6440.812
FIGO stage
 ≤IIA
 ≥IIB3.0130.617–7.6620.031*3.5460.860–8.7760.015*
Differentiation
 G1
 G20.8160.309–2.1550.6810.7890.267–2.3340.669
 G31.2510.718–2.1760.4311.1450.614–2.1350.671
Tumor size, cm
 <4
 ≥45.8873.025–11.456<0.001*4.1572.009–8.602<0.001*
Histological type
 Squamous
 others1.1080.305–4.0270.8770.6670.201–2.2150.508
Lymph invasion
 No
 Yes4.3150.750–12.3060.041*4.4580.764–12.7800.023*
m6A level
 Low4.9012.469–9.721<0.001*4.6382.152–9.997<0.001*
 High
m6A regulator mutation
 No
 Yes6.1183.004–12.462<0.001*6.3482.875–14.014<0.001*

*P < 0.05 indicates a significant association among the variables.

Table 4

Multivariate Cox proportional hazards model for DFS and OS
DFSOS
HR95 % CIP valueHR95 % CIP value
Lymph invasion2.7961.919–4.161<0.001*2.6591.711–4.223<0.001*
m6A level1.7011.129–2.5410.008*3.9811.854–9.173<0.001*
m6A regulator mutation4.4021.299–14.5510.011*8.0012.403–26.815<0.001*

*P < 0.05 indicates a significant association among the variables.

Taken together, the data showed that the reduced m6A level might contribute to cervical cancer development and the poor outcomes.

Reducing m6A level promotes cervical cancer cell proliferation

To uncover the important role of m6A in cervical cancer development, cell proliferation was analyzed in human cervical cancer cell line SiHa. The m6A methylation is regulated by adenosine methyltransferases and demethylases [711]. To reduce the m6A level, the adenosine methyltransferases (METTL3 and METTL14) were knocked down through two different shRNA. The knock-down efficiency was evaluated with qPCR (Figure 2A-2C, Supplementary Figure 3A). Knocking down both METTL3 and METTL14 promoted SiHa cell proliferation (Figure 2B, 2D). On the other hand, overexpressing adenosine demethylases (FTO and ALKBH5) also decreases m6A level. The FTO and ALKBH5 overexpression were validated by western blot (Figure (Figure2E,2E, Supplementary Figure 3B). And their overexpression significantly promoted cell proliferation (Figure (Figure2F).2F). Furthermore, the cell motility was also increased by knocking down METTL3 and METTL14 or overexpressing FTO and ALKBH5 (Supplementary Figure 4).

An external file that holds a picture, illustration, etc.
Object name is oncotarget-08-98918-g002.jpg
Reducing m6A level promoted cervical cancer cell proliferation

(A) METTL3 shRNA knocking-down efficiency was determined by qPCR in cervical cancer cell line SiHa. n = 3. * P<0.05. (B) Knocking down METTL3 promoted cell proliferation. n = 3. * P<0.05. (C) METTL14 shRNA knocking-down efficiency was determined by qPCR. n = 3. * P<0.05. (D) Knocking down METTL14 promoted cell proliferation. n = 3. * P<0.05. (E) Overexpressing FTO and ALKBH5 were validated via western blot. (F) Overexpressing FTO and ALKBH5 promoted cell proliferation. n = 3. * P<0.05.

Thus, reducing m6A level could promote cervical cancer cell proliferation, indicating that increasing m6A level might have anti-cancer effects in cervical cancer.

Cervical cancer cell proliferation could be suppressed by increasing m6A level

As the reduced m6A level was tightly associated with poor outcomes and reducing m6A in cervical cancer cell lines could promote cell proliferation, we wondering whether increasing m6A would have anti anti-cancer effects or not in cervical cancer. Indeed, increasing m6A through knocking down adenosine demethylases (FTO and ALKBH5) or overexpressing adenosine methyltransferases (METTL3 and METTL14) suppressed the cancer cell proliferation (Figure 3A-3F, Supplementary Figure 3C-3D). Furthermore, the cell motility was also decreased by knocking down FTO and ALKBH5 or overexpressing METTL3 and METTL14 (Supplementary Figure 5). These data demonstrated that increasing m6A might inhibit cervical cancer development.

An external file that holds a picture, illustration, etc.
Object name is oncotarget-08-98918-g003.jpg
Increasing m6A level suppressed cervical cancer cell proliferation

(A) FTO shRNA knocking-down efficiency was determined by qPCR in cervical cancer cell line SiHa. n = 3. * P<0.05. (B) Knocking down FTO suppressed cell proliferation. n = 3. * P<0.05. (C) ALKBH5 shRNA knocking-down efficiency was determined by qPCR. n = 3. * P<0.05. (D) Knocking down ALKBH5 suppressed cell proliferation. n = 3. * P<0.05. (E) Overexpressing METTL3 and METTL14 was validated via western blot. (F) Overexpressing METTL3 and METTL14 suppressed cell proliferation. n = 3. * P<0.05.

Increasing m6A level inhibits cervical cancer development in vivo

To further confirm the important role of m6A level on cervical cancer development, human cervical cancer cell line SiHa was first infected with lentivirus expressing shRNA to knock down adenosine methyltransferases (METTL3 and METTL14) or overexpressing adenosine demethylases (FTO and ALKBH5). Then the cells were transplanted subcutaneously into the dorsal scapula region of the NOD/SCID mice. Data showed that decreasing m6A level promoted cervical cancer development in vivo (Figure (Figure4A).4A). On the other hand, knocking-down adenosine demethylases (FTO and ALKBH5) or overexpressing adenosine methyltransferases (METTL3 and METTL14) suppressed cervical cancer development in vivo (Figure (Figure4B).4B). Furthermore, when treating the mice with FTO Inhibitor MA2 [32], the tumor size was significantly reduced (Figure 4C-4E). These data showed that increasing m6A level inhibits cervical cancer development in vivo.

An external file that holds a picture, illustration, etc.
Object name is oncotarget-08-98918-g004.jpg
Increasing m6A level suppressed cervical cancer development in vivo

(A) Knocking down METTL3 and METTL14 or overexpressing FTO and ALKBH5 promoted cervical cancer development in vivo. n = 12. * P<0.05. (B) Knocking down FTO and ALKBH5 or overexpressing METTL3 and METTL14 suppressed cervical cancer development in vivo. n = 12. * P<0.05. (C) FTO inhibitor MA2 treatment suppressed cervical cancer development in vivo. n = 12. * P<0.05. (D) Representative luciferase signal density measured via IVIS 50 imaging system. (E) Representative tumors isolated from the mice treated with MA2 or control.

DISCUSSION

Despite the current treatments for cervical cancer having improved the survival significantly, the development of drug resistance still occurs in a great number of patients determining recurrence [1, 2]. Therefore, more efforts should be made to uncover the underlying mechanisms and develop novel therapeutic targets.

In the present study, data showed that the m6A level was significantly decreased in cervical cancer tissues comparing with tumor paired normal tissues. Further analysis showed that the reduced m6A level tightly associated with cancer progression and poor survival, and could be developed as novel prognostic marker to predict tumor recurrence.

Then the m6A level was reduced in human cervical cancer cell lines through knocking down the adenosine methyltransferases or overexpressing adenosine demethylases, both of which could promote cell proliferation. On the other hand, increasing m6A level suppressed the cancer cell proliferation. The potential therapeutic application of increasing m6A was further validated in vivo.

In conclusion, the data showed that cervical cancer cells had very low level of m6A. Increasing m6A level impeded the cervical cancer development. Therefore, m6A might be a potential target for cervical cancer treatment.

MATERIALS AND METHODS

Patients

A total of 286 pairs of samples were obtained from patients with primary cervical cancer who had undergone surgery without any preoperative therapy at Qilu Hospital of Shandong University between January 2005 and December 2009. Patients staged according to the International Federation of Gynecology and Obstetrics (FIGO) criteria. The histological subtype was assigned according to the criteria of the World Health Organization classification. The clinical and pathologic parameters were reviewed from impatient medical records and presented in Table Table1.1. Samples were collected from the patients consecutively with the standardized protocol. Every patient specimen included two matched pairs, namely, cervical cancer tissues and adjacent normal tissues (≥5 cm away from the tumor). Surgically resected specimens were collected immediately after tumor removal and divided into two aliquots: half were immediately flash-frozen in liquid nitrogen and then frozen at −80°C until RNA and DNA extraction was performed; the remainder was fixed with formalin for histopathological analysis. The diagnosis was confirmed by at least two pathologists. The study was conducted according to the Declaration of Helsinki and approved by the Ethics Committee of Qilu Hospital of Shandong University. Written informed consent was obtained from all patients.

RNA extraction and m6A dot blot assay

Total RNA was extracted from tissues with Trizol (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. mRNA was prepared from total RNA using the Dynabeads mRNA purification kit (Ambion, catalog no. 61006). RNA samples were quantified using UV spectrophotometry, and equal amounts were mixed 1:1 with glyoxal loading dye (Ambion) and denatured for 20 min at 50°C. 500ng mRNAs were spotted onto a nylon membrane (GE Healthcare). RNA was UV crosslinked to the membrane, and the membranes were blocked for 1 hr in 5% nonfat dry milk in 0.1% PBST (0.1% Tween-20 in 1x PBS, pH 7.4) (Blocking Buffer). Anti-m6A antibody (Synaptic Systems; catalog no. 202003) was diluted 1:1000 in 0.1% PBST and incubated on the membranes for 1 hr (25°C) to overnight (4°C). Following extensive washing with 0.1% PBST, HRP-conjugated donkey anti-rabbit IgG (GE Healthcare) was diluted 1:2500 in Blocking Buffer and added to the membranes for 1 hr at 25°C. Membranes were washed again in 0.1% PBST and developed with enhanced chemiluminescence (ECL; GE Healthcare). The intensity of dot blot signal was quantified by ImageJ.

Gene mutation analysis

Genomic DNA was extracted from tissues with QIAamp DNA Blood Mini Kit (Qiagen) according to the manufacturer's instructions. Exons were amplified with high fidelity DNA polymerase (NEB) and sequenced with ABI 3730. Sequence deletion, truncation and missense mutations were characterized as gene mutation in the present study.

Real-time PCR

cDNA was prepared by using the iScript™ cDNA Synthesis kit (Bio-Rad, USA). PCR primers (Generay, Shanghai, China) used for RT-PCR were as follows: for METLL3, sense: 5′-TCCTGACTGACCTTCTTGCTC-3′ and anti-sense: 5′- TCAGCATCGGAA CCAGCAAAG-3′; METLL14, sense: 5′- GTTGGAACA TGGATAGCCGC -3′ and anti-sense: 5′-CAATGCTGTCGGCACTTTCA -3′; FTO, sense: 5′- AACACCAGGCTCTTTACGGTC -3′ and anti-sense: 5′- TGTCCGTTGTAGGATGAACCC -3′; ALKBH5, sense: 5′- ATGCACCCCGGTTGGAAAC -3′ and anti-sense: 5′- GACTTGCGCCAGTAGTTCTCA-3′;β-actin, sense: 5′-CCTGACTGACTACCTCATGAAG-3′ and anti-sense: 5′-GACGTAGCACAGCTTCTCCTTA-3′. RT-PCR amplification reaction was prepared with the SYBR Green PCR kit (Bio-rad, USA) and performed using the 7500 fast Real-Time PCR system (Applied Biosystems, USA). PCR products were verified by melting curve analysis. Relative mRNA levels of target genes were calculated by the 2-ΔΔct method.

Western blotting

Total protein from cultured cells were lysed in RIPA buffer with protease inhibitor (Beyotime, Shanghai, China). The protein was quantified using a BCA assay kit (Beyotime, Shanghai, China). A total of 20 μg of total protein were separated by 10% SDS-PAGE, transferred onto polyvinylidene fluoride membranes, and then reacted with primary antibodies against METLL3, METLL14, FTO, ALKBH5 and β-actin (all from Abcam, Cambridge, UK). After being extensively washed with PBS containing 0.1 % Triton X-100, the membranes were incubated with alkaline phosphatase-conjugated goat anti-rabbit antibody for 30 min at room temperature. The bands were visualized using 1-step TM NBT/BCIP reagents (Thermo Fisher Scientific, Rockford, IL, USA) and detected by an Alpha Imager (Alpha Innotech, San Leandro, CA, USA).

Cell culture

The cervical cancer cell lines SiHa was obtained from the American Type Culture Collection (ATCC; Rockville, MD, USA) and cultured in DMEM (GIBCO, Shanghai, China) supplemented with 10 % FBS.

Plasmid DNA and virus preparation

shRNAs were cloned into lentiviral pLKO.1-puro vector. The following shRNA sequences were used: control shRNA, 5`-ACTCAAAAGGAAGTG ACAAGA-30; METTL3 shRNA-1, 50-GCTGCACTTCAGACG AATT-30; METTL3 shRNA-2, 5`-CCACCTCAGTGGATC TGTT-3`; METTL14 shRNA-1, 5`-GCTAAAGGATGAG TTAAT-3`; METTL14 shRNA-2, 5`-GGACTTGGGATG ATATTAT- 3`; FTO shRNA-1, 5`-GGAAGATTTAAACT CCATGAAG-3`; FTO shRNA-2, 5`-CAAAGTGTTCAA TGGATGCAAC-3`; ALKBH5 shRNA-1, 5`-AGGTTCTC ATATTCTTGGTATC-3`; ALKBH5 shRNA-2, 5`-GATGAAATCACTCACTGCATAC-3`. The METTL3, METTL14, FTO and ALKBH5 expressing lentiviral vector was prepared by cloning the human full coding sequences into the pLVX-puro lentiviral vector.

Lentiviruses were prepared using 293T cells according to the manufacturer's instructions. Cells were incubated with lentivirus and 4 mg/mL polybrene (AmericanBio) for 24 hr.

Cell growth assay

Cells were seeded at 5×104 cells per well in 24-well plates and cultured for 3 days. Cell number was counted using a hemocytometer.

Animal study

6–8 weeks old NOD/SCID mice (Charles River Laboratories, Beijing, China) were housed in specific pathogen-free conditions. The study was approved by the Research Ethics Committee of Shandong University. Mice were housed in the pathogen free region and monitored daily during the experiments and the mice would be sacrificed when the weight loss is more than 20%. For evaluation of the tumor growth in vivo, 5 × 106 cells were suspended in 200 μl PBS and injected subcutaneously into the dorsal scapula region of the mice. For FTO Inhibitor MA2 treatment, the mice were treated with MA2 (5 nmol in 100μL PBS) or vehicle control by intratumoral injection once a week for 4 weeks. Tumor size was measured with fine digital calipers and calculated by the following formula: tumor volume =0.5 × width2 × length. For imaging the tumor clearance in vivo, cells were transduced with a lentiviral construct that drives the expression of the click beetle red luciferase (CBR). Bioluminescence imaging was performed 4 weeks after injection. For imaging, mice were injected i.p. with D-luciferin (150 μg/g body weight), anesthetized (2% isoflurane), and placed in an IVIS 50 imaging system (Xenogen). Regions of interest (ROI) were defined manually over the whole body with Living-Image software (Igor Wavemetrics) for determining tumor burden signal intensities.

Statistical analysis

Data were expressed as mean (±SE) and analyzed by a SPSS software package (SPSS Standard version 13.0, SPSS Inc, USA). Differences between variables were assessed by the Chi-square test. Survival analysis of patients with cervical cancer was calculated by Kaplan-Meier analysis. A log rank test was used to compare different survival curves. A Cox proportional hazards model was used to calculate univariate and multivariate hazard ratios for the variables. Unpaired Student's t test and one way ANOVA were used as appropriate to assess the statistical significant of difference. P values under 0.05 were considered statistically significant.

SUPPLEMENTARY MATERIALS FIGURES

Footnotes

CONFLICTS OF INTEREST

The authors declare that they have no competing interests.

FUNDING

This work was supported by the Seeding Fund to Beihua Kong (81672427).

REFERENCES

1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016;66:7–30. https://doi.org/10.3322/caac.21332 [Abstract] [Google Scholar]
2. Siegel R, Naishadham D, Jemal A. Cancer statistics, 2012. CA Cancer J Clin. 2012;62:10–29. https://doi.org/10.3322/caac.20138 [Abstract] [Google Scholar]
3. Jaffrey SR, Kharas MG. Emerging links between m6A and misregulated mRNA methylation in cancer. Genome Med. 2017;9:2. https://doi.org/10.1186/s13073-016-0395-8 [Europe PMC free article] [Abstract] [Google Scholar]
4. Machnicka MA, Milanowska K, Osman Oglou O, Purta E, Kurkowska M, Olchowik A, Januszewski W, Kalinowski S, Dunin-Horkawicz S, Rother KM, Helm M, Bujnicki JM, Grosjean H. MODOMICS: a database of RNA modification pathways--2013 update. Nucleic Acids Res. 2013;41:D262–7. https://doi.org/10.1093/nar/gks1007 [Europe PMC free article] [Abstract] [Google Scholar]
5. Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, Cesarkas K, Jacob-Hirsch J, Amariglio N, Kupiec M, Sorek R, Rechavi G. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 2012;485:201–6. https://doi.org/10.1038/nature11112 [Abstract] [Google Scholar]
6. Meyer KD, Saletore Y, Zumbo P, Elemento O, Mason CE, Jaffrey SR. Comprehensive analysis of mRNA methylation reveals enrichment in 3’ UTRs and near stop codons. Cell. 2012;149:1635–46. https://doi.org/10.1016/j.cell.2012.05.003 [Europe PMC free article] [Abstract] [Google Scholar]
7. Liu J, Yue Y, Han D, Wang X, Fu Y, Zhang L, Jia G, Yu M, Lu Z, Deng X, Dai Q, Chen W, He C. A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat Chem Biol. 2014;10:93–5. https://doi.org/10.1038/nchembio.1432 [Europe PMC free article] [Abstract] [Google Scholar]
8. Ping XL, Sun BF, Wang L, Xiao W, Yang X, Wang WJ, Adhikari S, Shi Y, Lv Y, Chen YS, Zhao X, Li A, Yang Y, et al. Mammalian WTAP is a regulatory subunit of the RNA N6-methyladenosine methyltransferase. Cell Res. 2014;24:177–89. https://doi.org/10.1038/cr.2014.3 [Europe PMC free article] [Abstract] [Google Scholar]
9. Schwartz S, Mumbach MR, Jovanovic M, Wang T, Maciag K, Bushkin GG, Mertins P, Ter-Ovanesyan D, Habib N, Cacchiarelli D, Sanjana NE, Freinkman E, Pacold ME, et al. Perturbation of m6A writers reveals two distinct classes of mRNA methylation at internal and 5’ sites. Cell Rep. 2014;8:284–96. https://doi.org/10.1016/j.celrep.2014.05.048 [Europe PMC free article] [Abstract] [Google Scholar]
10. Jia G, Fu Y, Zhao X, Dai Q, Zheng G, Yang Y, Yi C, Lindahl T, Pan T, Yang YG, He C. N6-methyladenosine in nuclear RNA is a major substrate of the obesity-associated FTO. Nat Chem Biol. 2011;7:885–7. https://doi.org/10.1038/nchembio.687 [Europe PMC free article] [Abstract] [Google Scholar]
11. Zheng G, Dahl JA, Niu Y, Fedorcsak P, Huang CM, Li CJ, Vagbo CB, Shi Y, Wang WL, Song SH, Lu Z, Bosmans RP, Dai Q, et al. ALKBH5 is a mammalian RNA demethylase that impacts RNA metabolism and mouse fertility. Mol Cell. 2013;49:18–29. https://doi.org/10.1016/j.molcel.2012.10.015 [Europe PMC free article] [Abstract] [Google Scholar]
12. Patil DP, Chen CK, Pickering BF, Chow A, Jackson C, Guttman M, Jaffrey SR. m(6)A RNA methylation promotes XIST-mediated transcriptional repression. Nature. 2016;537:369–73. https://doi.org/10.1038/nature19342 [Europe PMC free article] [Abstract] [Google Scholar]
13. Mauer J, Luo X, Blanjoie A, Jiao X, Grozhik AV, Patil DP, Linder B, Pickering BF, Vasseur JJ, Chen Q, Gross SS, Elemento O, Debart F, et al. Reversible methylation of m6Am in the 5′ cap controls mRNA stability. Nature. 2017;541:371–5. https://doi.org/10.1038/nature21022 [Europe PMC free article] [Abstract] [Google Scholar]
14. Liu N, Dai Q, Zheng G, He C, Parisien M, Pan T. N(6)-methyladenosine-dependent RNA structural switches regulate RNA-protein interactions. Nature. 2015;518:560–4. https://doi.org/10.1038/nature14234 [Europe PMC free article] [Abstract] [Google Scholar]
15. Alarcon CR, Goodarzi H, Lee H, Liu X, Tavazoie S, Tavazoie SF. HNRNPA2B1 Is a Mediator of m(6)A-Dependent Nuclear RNA Processing Events. Cell. 2015;162:1299–308. https://doi.org/10.1016/j.cell.2015.08.011 [Europe PMC free article] [Abstract] [Google Scholar]
16. Wang X, Lu Z, Gomez A, Hon GC, Yue Y, Han D, Fu Y, Parisien M, Dai Q, Jia G, Ren B, Pan T, He C. N6-methyladenosine-dependent regulation of messenger RNA stability. Nature. 2014;505:117–20. https://doi.org/10.1038/nature12730 [Europe PMC free article] [Abstract] [Google Scholar]
17. Zhou J, Wan J, Gao X, Zhang X, Jaffrey SR, Qian SB. Dynamic m(6)A mRNA methylation directs translational control of heat shock response. Nature. 2015;526:591–4. https://doi.org/10.1038/nature15377 [Europe PMC free article] [Abstract] [Google Scholar]
18. Wang X, Zhao BS, Roundtree IA, Lu Z, Han D, Ma H, Weng X, Chen K, Shi H, He C. N(6)-methyladenosine Modulates Messenger RNA Translation Efficiency. Cell. 2015;161:1388–99. https://doi.org/10.1016/j.cell.2015.05.014 [Europe PMC free article] [Abstract] [Google Scholar]
19. Meyer KD, Patil DP, Zhou J, Zinoviev A, Skabkin MA, Elemento O, Pestova TV, Qian SB, Jaffrey SR. 5′ UTR m(6)A Promotes Cap-Independent Translation. Cell. 2015;163:999–1010. https://doi.org/10.1016/j.cell.2015.10.012 [Europe PMC free article] [Abstract] [Google Scholar]
20. Jia G, Fu Y, He C. Reversible RNA adenosine methylation in biological regulation. Trends Genet. 2013;29:108–15. https://doi.org/10.1016/j.tig.2012.11.003 [Europe PMC free article] [Abstract] [Google Scholar]
21. Batista PJ, Molinie B, Wang J, Qu K, Zhang J, Li L, Bouley DM, Lujan E, Haddad B, Daneshvar K, Carter AC, Flynn RA, Zhou C, et al. m(6)A RNA modification controls cell fate transition in mammalian embryonic stem cells. Cell Stem Cell. 2014;15:707–19. https://doi.org/10.1016/j.stem.2014.09.019 [Europe PMC free article] [Abstract] [Google Scholar]
22. Wang Y, Li Y, Toth JI, Petroski MD, Zhang Z, Zhao JC. N6-methyladenosine modification destabilizes developmental regulators in embryonic stem cells. Nat Cell Biol. 2014;16:191–8. https://doi.org/10.1038/ncb2902 [Europe PMC free article] [Abstract] [Google Scholar]
23. Geula S, Moshitch-Moshkovitz S, Dominissini D, Mansour AA, Kol N, Salmon-Divon M, Hershkovitz V, Peer E, Mor N, Manor YS, Ben-Haim MS, Eyal E, Yunger S, et al. Stem cells. m6A mRNA methylation facilitates resolution of naive pluripotency toward differentiation. Science. 2015;347:1002–6. https://doi.org/10.1126/science.1261417 [Abstract] [Google Scholar]
24. Klungland A, Dahl JA. Dynamic RNA modifications in disease. Curr Opin Genet Dev. 2014;26:47–52. https://doi.org/10.1016/j.gde.2014.05.006 [Abstract] [Google Scholar]
25. Blanco S, Frye M. Role of RNA methyltransferases in tissue renewal and pathology. Curr Opin Cell Biol. 2014;31:1–7. https://doi.org/10.1016/j.ceb.2014.06.006 [Europe PMC free article] [Abstract] [Google Scholar]
26. Zhang C, Zhi WI, Lu H, Samanta D, Chen I, Gabrielson E, Semenza GL. Hypoxia-inducible factors regulate pluripotency factor expression by ZNF217- and ALKBH5-mediated modulation of RNA methylation in breast cancer cells. Oncotarget. 2016;7:64527–42. https://doi.org/10.18632/oncotarget.11743 [Europe PMC free article] [Abstract] [Google Scholar]
27. Zhang C, Samanta D, Lu H, Bullen JW, Zhang H, Chen I, He X, Semenza GL. Hypoxia induces the breast cancer stem cell phenotype by HIF-dependent and ALKBH5-mediated m(6)A-demethylation of NANOG mRNA. Proc Natl Acad Sci U S A. 2016;113:E2047–56. https://doi.org/10.1073/pnas.1602883113 [Europe PMC free article] [Abstract] [Google Scholar]
28. Bansal H, Yihua Q, Iyer SP, Ganapathy S, Proia DA, Penalva LO, Uren PJ, Suresh U, Carew JS, Karnad AB, Weitman S, Tomlinson GE, Rao MK, et al. WTAP is a novel oncogenic protein in acute myeloid leukemia. Leukemia. 2014;28:1171–4. https://doi.org/10.1038/leu.2014.16 [Europe PMC free article] [Abstract] [Google Scholar]
29. Kwok CT, Marshall AD, Rasko JE, Wong JJ. Genetic alterations of m6A regulators predict poorer survival in acute myeloid leukemia. J Hematol Oncol. 2017;10:39. https://doi.org/10.1186/s13045-017-0410-6 [Europe PMC free article] [Abstract] [Google Scholar]
30. Zhang S, Zhao BS, Zhou A, Lin K, Zheng S, Lu Z, Chen Y, Sulman EP, Xie K, Bogler O, Majumder S, He C, Huang S. m6A Demethylase ALKBH5 Maintains Tumorigenicity of Glioblastoma Stem-like Cells by Sustaining FOXM1 Expression and Cell Proliferation Program. Cancer Cell. 2017;31:591–606. https://doi.org/10.1016/j.ccell.2017.02.013 e6. [Europe PMC free article] [Abstract] [Google Scholar]
31. Cui Q, Shi H, Ye P, Li L, Qu Q, Sun G, Sun G, Lu Z, Huang Y, Yang CG, Riggs AD, He C, Shi Y. m6A RNA Methylation Regulates the Self-Renewal and Tumorigenesis of Glioblastoma Stem Cells. Cell Rep. 2017;18:2622–34. https://doi.org/10.1016/j.celrep.2017.02.059 [Europe PMC free article] [Abstract] [Google Scholar]
32. Huang Y, Yan J, Li Q, Li J, Gong S, Zhou H, Gan J, Jiang H, Jia GF, Luo C, Yang CG. Meclofenamic acid selectively inhibits FTO demethylation of m6A over ALKBH5. Nucleic Acids Res. 2015;43:373–84. https://doi.org/10.1093/nar/gku1276 [Europe PMC free article] [Abstract] [Google Scholar]

Articles from Oncotarget are provided here courtesy of Impact Journals, LLC

Citations & impact 


Impact metrics

Jump to Citations

Citations of article over time

Alternative metrics

Altmetric item for https://www.altmetric.com/details/99564937
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/99564937

Article citations


Go to all (64) article citations

Data 


Data behind the article

This data has been text mined from the article, or deposited into data resources.

Similar Articles 


To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.