Knockdown of RBM15 inhibits tumor progression and the JAK-STAT signaling pathway in cervical cancer
BMC Cancer volume 23, Article number: 684 (2023)
RNA binding motif protein 15 (RBM15), a writer of N6-methyladenosine (m6A) methylation, contributes significantly to the development of various tumors. However, the function of RBM15 in cervical cancer (CC) has not been determined.
Based on the GSE9750, GSE63514, and m6A datasets, m6A-related differentially expressed genes (DEGs) were screened out. The hub genes were identified by generating a Protein-Protein Interaction (PPI) network. RT-qPCR was conducted to assess the mRNA expression of hub genes. CCK8, scratch wound healing, and transwell assays were utilized to examine the influence of RBM15 on HeLa and SiHa cells. Tumor xenograft models were used to assess the effects of RBM15 on tumorigenesis. A mechanistic analysis of RBM15 in CC tumors was conducted using the GeneCards and Coxpresdb databases, followed by a Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, and the pathway-related genes were subsequently validated using Western blotting.
Five DEGs were screened, including WTAP, RBM15, CBLL1, and YTHDC2. Among them, WTAP, RBM15, CBLL1, and YTHDC2 were hub genes and can be used as biomarkers for CC. RBM15 expression was considerably increased, while WTAP, CBLL1, and YTHDC2 were significantly downregulated. Knockdown of RBM15 significantly suppressed the proliferation, invasion, and migration of CC cells and tumorigenesis. Moreover, knockdown of RBM15 significantly reduced the expression levels of proteins related to the JAK-STAT pathway.
Knockdown of RBM15 inhibited the progression of CC cells, which probably by inhibiting the JAK-STAT pathway pathway.
Cervical cancer (CC) ranks second in the prevalence of female malignancy . The infection of high-risk human papillomavirus , early age at first intercourse , and smoking  are all risk factors for CC. CC tumorigenesis is a complex process involving transcriptomic, genetic, and epigenetic alterations . The stage at diagnosis is a crucial factor in determining the prognosis of CC, but most patients are at an intermediate to advanced stage  accompanied by high recurrence and metastasis , making treatment more difficult. Therefore, much more research is required to clarify the molecular mechanisms of CC.
N6-methyladenosine (m6A) is an epigenetic alteration, which is methylated at the N6 position of adenosine . Being the most prevalent RNA modification in mammals, m6A regulates multiple processes in RNA metabolism, including splicing, export, and translation . Several studies have demonstrated that m6A affects the progression of various tumors, including CC. However, research on the function and mechanism of m6A in CC is not been fully expounded. RBM15, the m6A-modified catalytic protein, combines with METTL3, WTAP, and METTL14 to form the m6A methyltransferase complex, which is accountable for the installation of m6A . Previous reports suggest that RBM15 plays a significant role in tumor development. The m6A site modified by RBM15 is recognized by IGF2BP3, enhancing laryngeal squamous cell carcinoma progression . RBM15 accelerates the development of hepatocellular carcinoma in an IGF2BP1-dependent way . Knockdown of RBM15 suppresses the growth and metastasis of colorectal cancer . However, the function of RBM15 in CC has not been elucidated.
The Janus kinase/signal transducer and activator of transcription (JAK/STAT) signaling pathway is widely recognized as a pivotal hub in cellular function, serving as a central node for intercellular communication . The JAK/STAT signaling pathway plays an important role in tumor development. The activation of JAK/STAT signaling leads to the upregulation of various proteins involved in cell proliferation, cell survival, stemness, self-renewal, evasion of immune surveillance mechanisms, and overall tumor progression . Enhancement of the epithelial-mesenchymal transition through the activation of the JAK/STAT3 signaling pathway results in increased tumorigenic and metastatic potential, facilitates the transition of cancer stem cells, and contributes to chemoresistance in cancer . In CC, the JAK/STAT pathway has important correlation with HPV E6/E7 oncoprotein . And the activation of the JAK-STAT pathway by CXCL10 and CXCR3 leads to heightened expression of exosomal PD-L1, which may contribute to HPV evasion of the immune response and subsequent development of cancer . In addition, studies have shown that the JAK/STAT signaling pathway is involved in the progression of multiple tumors by regulating m6A-related genes. The cell cycle regulation in human glioma cells is influenced by YTHDC1-mediated VPS25, which specifically targets the JAK-STAT signaling pathway . ZC3H13 has been identified as a factor influencing the progression of hepatocellular carcinoma, which potentially impacts the JAK-STAT signaling pathway . However, the correlation between the JAK-STAT signaling pathway and m6A has rarely been reported in CC.
Bioinformatics allows for faster and comprehensive insights into disease mechanisms and treatments. For example, over the course of just a year, from 2021 to 2022, the Ren et al. team make a series of significant discoveries in Kidney disease. They identifies a biological signature and predictive model for biopsy-proven acute rejection after kidney transplantation, one o six stemness-associated genes predicting overall survival in renal clear cell carcinoma: AC010973.2, and RNA-binding protein genes and drug candidates associated with prognosis in renal papillary cell carcinoma, as well as also examined the DNA damage and prostate cancer correlation, which are investigated by bioinformatic methods [21,22,23,24]. And research on the molecular mechanisms of m6A modification in tumors has been greatly facilitated by the rapid development of bioinformatics. In this study, we selected GSE9750 and GSE63514 datasets from the GEO database to screen DEGs in CC. DEGs from each dataset were then intersected with the m6A data to screen for key genes related to m6A. RBM15 was finally identified to be further studied. The function of RBM15 in CC progression and its specific mechanism were studied.
Materials and methods
The GSE9750 and GSE63514 datasets were acquired from the Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/). The GSE9750 dataset (GPL96, [HG-U133A] Affymetrix Human Genome U133A Array) contains 33 primary tumors and 24 normal cervical epithelia, and the GSE63514 dataset (GPL570, [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array) contains 28 tumor tissue samples and 24 normal tissue samples.
Screening for DEGs
The GSE9750 and GSE63514 datasets were subjected to analysis by GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r) to obtain gene expression data. DEGs were screened according to the p < 0.05 and |logFC| > 0.5. By reviewing the m6A2Target database (http://rm2target.canceromics.org/), we collected m6A-related genes. The Venn diagram tool (http://bioinformatics.psb.ugent.be/webtools/Venn/) was employed to identify the m6A-related DEGs across three datasets: GSE9750, GSE63514, and the m6A2Target datasets.GO and KEGG enrichment analysis of DEGs.
Functional enrichment analysis of the DEGs was analyzed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) in the DAVID (Database for Annotation, Visualization and Integrated Discovery) database (https://david.ncifcrf.gov/), which is a comprehensive database encompassing genomic, chemical, and systemic functional information .
PPI network construction and m6A-related hub DEGs determination
Using the STRING database (https://string-db.org/) with a confidence score of 0.15, a PPI (protein-protein interaction interaction) network of the DEGs was established. We then screened out the linkage data for proteins associated with m6A and loaded them into the Cytoscape software (version 3. 9. 1) for visualization. Hub modules were identified using the Molecular Complex Detection (MCODE) plug-in. The plug-in of CytoHubba was used to calculate each protein degree of nodes and filtered out hub genes based on connectivity.
Analysis of m6A-related hub DEGs
Receiver operating characteristic (ROC) curves, Principal component analysis (PCA), and a ridge plot of the hub genes expression abundance were performed by the pROC package, factoextra package, and ggridges package, respectively in R software (version 3.6.0). In the GEPIA2 database (http://gepia2.cancer-pku.cn/), the mRNA expression levels of hub genes were explored in CC tissues. In the Kaplan-Meier Plotter database (https://kmplot.com/analysis/), the effect of RBM15 expression on the survival of CC patients was evaluated. The correlation between RBM15 and CD4 + T cells was analyzed by the TIMER2.0 database (http://timer.cistrome.org/).
Normal cervical epithelial cell line: HcerEpic, and CC cell line: SiHa and Hela (Icellbioscience, Shanghai, China) were grown in DMEM medium (Gibco, NY, USA) containing 1% penicillin and streptomycin, and 10% fetal bovine serum (Gibco). Cells were placed in an incubator (Thermo Fisher Scientific, Massachusetts, USA) in a humid environment, and the temperature of the incubator was set to 37 °C with a 5% CO2 condition.
Reverse transcription quantitative PCR (RT-qPCR)
Cells were subjected to RNA extraction using TRIzol Reagent (Thermo Fisher Scientific) to obtain the complete RNA. Following the evaluation of RNA samples for quantification and quality using a NanoDrop spectrophotometer (Thermo Fisher Scientific) and standard agarose gel electrophoresis, respectively, the RNA was subjected to reverse transcription into cDNA using a reverse transcription kit (Thermo Fisher Scientific). The PCR quantification was performed using the 7500 real-time PCR System (Thermo Fisher Scientific) with the SYBR Premix Ex Taq reagent (Takara, Dalian, China). The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used to normalize the results, and the relative expression levels of the mRNAs were determined using the 2−ΔΔCT method. Supplementary Table 1 contains the primer sequences used in the study.
The siRNA transfection was utilized to transiently knock down RBM15 expression. Three types of the small interfering RNA (siRNA) aimed at human RBM15 and the non-targeting control siRNA, were procured from Shanghai Integrated Biotech Solutions Co., Ltd. (Shanghai, China). A total of 1 × 105 cells/well SiHa and Hela cells were incubated in a 24-well plate until they attained a confluence of 70–90%. Subsequently, they were subjected to transfection with 50 nmol siRBM15 and negative control (NC), with the Lipofectamine 3000 reagent (Thermo Fisher Scientific), according to the manufacturer’s instructions. The effectiveness of silencing RBM15 was validated by RT-qPCR analysis. Supplementary Table 1 displays the sequences of the siRNAs used.
Cell counting kit-8 (CCK8) assays
The transfected 100 µL cell suspension (The density was 1 × 105 cells/well) was seeded into 96 well plates. Subsequently, 10 µL of CCK8 solution (Solarbio, Beijing, China) was added to each well at 0, 24, 48, 72, and 96 h, respectively. Following a 2-hour incubation, the absorbance at 450 nm was measured with a microplate reader (DALB, Shanghai, China).
Scratch wound healing
The density of cells in each group was adjusted to 5 × 105/well and then inoculated into six-well plates. When the cells had reached 80% confluence, a 10 µL pipette was used to draw a horizontal line on the cell surface. Images were captured immediately after scratching 24 h later. The migration distances and migration rates were calculated using the Image J program.
Following transfection for 48 h, we used serum-free medium to resuspend the cells. After coating the transwell chambers with matrigel gel (Becton, Dickinson and Company, New York, USA), in the upper chamber, cells were added, and in the lower chamber, 600 µL of medium supplemented with 15% fetal bovine serum was added. Following a 24-hour incubation, cells were fixed in 4% paraformaldehyde at room temperature for 15 min before being stained with 0.1% crystal violet at 37 °C for 30 min. Finally, cells were observed and photographed using a light microscopy (OLYMPUS, Tokyo, Japan).
Mice tumor xenograft model experiments
Six-week-old, female BALB/c nude mice (SPF (Beijing) Biotechnology Co, Ltd, Beijing, China) were randomly allocated into two groups (n = 6 in each group). To establish a tumor xenograft model, 0.1 ml Hela cell (transfected with RBM15 shRNA or control shRNA) suspension at a concentration of 1 × 107/ml was injected subcutaneously into the right axillary region of mice. The growth of the xenograft tumor was checked weekly after the first administration, and the tumour diameter was determined by a vernier caliper. The calculation of the tumor volume was performed by the formula V = 1/2ab2 (a is the maximum tumour diameter, and b is the minimum tumour diameter) and continuously measured for 4 weeks. Mice were executed after the last measurement, and the tumors were weighed. Our Animal Ethics Committee approved all the above experiments.
Western blot analysis
Proteins extracted from tumor tissues and cells were lysed using RIPA lysis buffer containing a protein inhibitor (Solarbio). The protein concentrations of the samples were determined by the BCA assay (Thermo Fisher Scientific). The protein was loaded onto a sodium dodecyl sulfate–polyacrylamide gel to complete separation. Then they were transferred to the polyvinylidene difluoride (PVDF) membranes (Roche, Basel, Switzerland), which was conducted at a current of 200 mA for approximately 50–60 min, depending on the size of the proteins. The PVDF membrane was subjected to blocking in 5% bovine serum albumin (BSA; Thermo Fisher Scientific) at room temperature for 1 h and then was incubated overnight at 4 °C with primary antibody. Following three washes in TBST (10 min each), the membrane was incubated with the secondary antibody (1:5000, Abcam) at room temperature for 2 h. Subsequently, the membrane was washed three times with TBST. The immunoblots were visualized using an ECL luminescent solution (Amersham, Little Chalfont, UK) as per the instructions provided by the manufacturer. The resulting immunoblots were then analyzed using Image J software. And the antibodies utilized were as listed below: anti-RBM15 (1:2000, Abcam),, P-JAK2 (1:2000, Abcam), JAK2 (1:2000, Abcam), P-STAT3 (1:2000, Abcam), STAT3 (1:2000, Abcam), P-STAT5 (1:2000, Abcam), STAT5 (1:2000, Abcam), and anti-β-actin (1:2000, Abcam).
Mechanisms of RBM15 affecting CC progression
We performed a screening of the top 1000 genes related to CC using the GeneCards database (https://www.genecards.org), taking into consideration their relevance scores. Furthermore, the Coxpresdb database (https://coxpresdb.jp/) was utilized to select the top 1000 genes that were co-expressed with the RBM15 gene. By creating a Venn diagram, we identified overlapping genes between the two databases. Then KEGG enrichment analysis of overlapping genes was performed by DAVID database.
Statistical data analysis
The data obtained from the experiment were processed using Graphpad Prism 8.0 statistical software, and the plots were generated using R software (version 3.6.1). The measurement data were expressed as mean ± standard deviation (SD). For the comparison of two groups, the T-test was used, while analysis of variance (ANOVA) was employed to compare multiple groups. Following ANOVA analysis, Tukey’s multiple comparisons test was used to compare two groups. A p-value less than 0.05 was considered statistically significant.
Identification of DEGs in CC
In the GEO database, the GSE9750 and GSE63514 datasets were selected to obtain DEGs by the GEO2R tool, and mapped the volcano (Fig. 1A and B). The read counts for each sample were normalized. The results revealed that the median values for each sample were near to identical, indicating that the GSE7670 and GSE63514 datasets fitted the criteria for further study (Supplementary Fig. 1). A total of 7644 DEGs were identified in the GSE63514 dataset and 5929 DEGs in the GSE9750 dataset (Supplementary Tables 2 and 3). The 25 overlapping most highly up- and down-regulated DEGs from the GSE9750 and GSE63514 datasets (Supplementary Table 4) were filtered out and plotted as a Heatmap using R software (Fig. 1C). According to the Venn diagram, there were 2682 DEGs that were commonly identified in both the GSE9750 and GSE63514 datasets. By reviewing the m6A2Target database, we collected 24 m6A-related genes (Supplementary Table 5) to the intersection with the overlapping DEGs for screening m6A-related DEGs, and five m6A-related DEGs were selected (Fig. 1D; Supplementary Table 6).
GO and KEGG analysis of the overlapping DEGs
We uploaded the overlapping DEGs to the DAVID database for GO and KEGG pathway analyses. GO functional enrichment analysis is divided into 3 categories, including MF (molecular function), BP (biological process), and CC (cellular component). The bar chart demonstrated the enrichment analysis results that ranked in the top 10 for every individual GO term, according to the p-value. In the MF group, the DEGs were predominantly involved in pathways related to cancer, human papillomavirus infection, and proteoglycans in cancer. In the BP, the DEGs were mainly involved in the membrane, nucleus, and cytoplasm. Regarding the CC category, the DEGs were found to be predominantly enriched in processes such as positive regulation of transcription from RNA polymerase II promoter, regulation of cell cycle, and DNA repair (Fig. 2A). Using the p-values as a basis, we filtered for the top 10 KEGG pathways and visualized them in a bubble plot (Fig. 2B). The findings demonstrated that the DEGs were primarily enriched in functions related to protein binding, identical protein binding, and ATP binding.
Construction of a protein-protein interaction (PPI) network and analysis of modules comprising m6A-related genes
A PPI network was created using the m6A-related DEGs (Fig. 3A). Then a highly interconnected cluster from the PPI network was determined using the MCODE in the Cytoscape software and four hub genes were obtained: WTAP, RBM15, CBLL1, and YTHDC2 (Fig. 3B).
Bioinformatics analysis was performed to identify hub genes
ROC curves were utilized to assess the precision of the four hub genes as CC biomarkers (Fig. 4A). The results showed that RBM15 had the highest AUC value (0.824), meanwhile, the AUC values for WTAP, CBLL1, and YTHDC2 were 0.653, 0.749, and 0.682, respectively. The ridge plot of four hub genes expression showed that RBM15 and WTAP had high expression levels (Supplementary Fig. 2). In addition, hub genes expressions in CC were explored in the GEPIA2 database. The findings indicated a significant increase in RBM15 expression in tumor tissues compared to non-tumor tissues. Conversely, the results demonstrated a significant decrease in YTHDC2 expression in tumor tissues compared to normal tissues. (Fig. 4B). Based on the above analysis, we finally determined RBM15 as the object of follow-up research.
Survival and immune infiltration analysis of RBM15
The results indicated a negative correlation between RBM15 expression and the survival of CC patients. The higher the expression of RBM15, the shorter the survival period of patients (Supplement Fig. 3A). In the TIMER2.0 database, there was a positive association between RBM15 and CD4 + T-cell infiltration in uterine corpus endometrial carcinoma (CESC), but the correlation coefficient was only 0.156 (Supplement Fig. 3B). Besides, there was a positive correlation between RBM15 expression and CD4 + T cell infiltration in a variety of tumors by different algorithms such as TIMER, EPIC, MCPCOUNTER, CIBERSORTT and XCELL (Supplement Fig. 3C).
Hub gene expression validation
We then examined the WTAP, RBM15, CBLL1, and YTHDC2 expressions in CC cells by RT-qPCR. According to the results, RBM15 expression was substantially increased in SiHa and Hela cells compared with HCerEpic cells. In contrast, WTAP, CBLL1, and YTHDC2 expressions were significantly downregulated (Fig. 5).
Knockdown of RBM15 inhibits proliferation, invasion and migration of CC cells
To further verify the role of RBM15 on CC cells, we knocked down RBM15 using siRNA. Figure 6 A showed the RBM15 expression in the si-RBM15-1, si-RBM15-2, and si-RBM15-3 groups. The si-RBM15-3 groups had the highest knockdown efficiency, so we applied si-RBM15-3 for our subsequent experiments. The results of the CCK-8 assay showed that the proliferation ability of SiHa and Hela cells was significantly reduced in the si-RBM15 group (Fig. 6B). The results from the Transwell and wound healing assays indicated that SiHa and Hela cells in the si-RBM15 group showed significantly reduced invasion and migration abilities compared to those in the si-NC group (Fig. 6C-D).
Knockdown of RBM15 inhibits the growth of CC in vivo
Then we verified the of RBM15 in CC growth by establishing mouse xenograft tumor models. The results indicated that the sh-RBM15 group exhibited significantly lower volume and weight of transplanted tumors compared to the sh-NC group (Fig. 7A-C). Figure 7D showed a significant reduction in RBM15 protein expression.
Silencing RBM15 inhibits the expression levels of proteins related to the JAK-STAT signaling pathway
Then, the molecular mechanisms of RMB15 affecting CC progression was explored. We screened CC-related top 1000 genes using the GeneCards database based on their relevance score. Additionally, the Coxpresdb database was employed to select the top 1000 genes co-expressed with the RBM15 gene. By drawing a Venn diagram, we identified 37 overlapping genes between the two databases (Fig. 8A). Then a KEGG enrichment analysis performed on the 37 overlapping genes, and the results revealed that they were significantly enriched in the JAK-STAT signaling pathway (Fig. 8B). Next, Western blot was used to detect the expression levels of JAK-STAT-signaling pathway-related proteins in SiHa and Hela cells transfected with siRBM15. The results showed that knockdown of RBM15 significantly decreased the proteins expression of P-JAK2/JAK2, P-STAT3/STAT3, and P-STAT5/STAT5 in SiHa and Hela cells (Fig. 8C and D).
Approximately 84% of all CC cases occur in countries with poor human development indexes. The majority of patients are diagnosed with CC at an advanced stage, and the prognosis is extremely poor, even in developed countries . Therefore, it is essential to find biomarkers and therapeutic targets for CC. In the present study, RBM15 was identified as an independent prognostic factor for CC, and the knockdown of RBM15 significantly inhibited the proliferation, invasion, and migration of CC cells and the JAK-STAT signaling pathway.M6A modification appears to be dynamic and reversible as it is regulated by methyltransferases, and demethylases. M6A recognition proteins play a crucial role in the function of m6A modifications . The m6A-associated proteins regulate multiple tumor progression. Upregulation of KIAA1429 expression in hepatocellular carcinoma induces m6A modification on pre-mRNA of GATA3, leading to increased tumor cell proliferation and metastasis . YTHDF1 is overexpressed in ovarian cancer and knockdown of YTHDF1 inhibits tumor progression . Several studies have elucidated the important role of m6A-associated proteins in CC development. METTL3 is overexpressed in CC and silencing of METTL3 can reduce cancer cell viability, promote apoptosis and increase cisplatin sensitivity . METTL14 expression is upregulated in CC, and knocking down METTL14 significantly inhibited tumor growth and invasion . YTHDF1 regulates RANBP2 translation efficiency in an m6A-dependent manner, promoting CC cell growth, migration and invasion . In this study, our results showed that RBM15 was noticeably upregulated in SiHa and Hela cells, while WTAP, CBLL1, and YTHDC2 expressions were significantly downregulated. Research has shown that WTAP is located in the nuclei and has a significant correlation with overall survival in CC patients . The prognosis of UCEC is influenced by missense mutations in YTHDC2 . The high expression of RBM15 in CC suggested a potential involvement as an oncogene, which had sparked significant interest among us. So we subsequently focused on the role of RBM15 in CC.
RBM15 belongs to the SPEN (split-end) protein family, which can be found on chromosome 1p13.3 . It has been reported that RBM15 directly binds to U-rich regions on mRNA to attract the methyltransferase complex to its target transcripts . Based on the current study, RBM15 was found to be an oncogenic factor in most tumors. RBM15-MKL1 fusion resulting from chromosomal translocation t(1;22)(p13;q13) can lead to acute megakaryocytic leukemia . RBM15 enhances the growth, metastasis, and infiltration of clear cell renal cell carcinoma . In addition, RBM15 interacts with circ-CTNNB1 to enhance the expression of HK2, GPI and PGK1 in an m6A modification-dependent manner, leading to the promotion of the glycolytic process and the development of osteosarcoma . In our study, we discovered that RBM15 functions as a tumor oncogene in CC, and knockdown of RBM15 inhibits cell proliferation, migration, and invasion, as well as tumor growth in vivo.
Subsequently, we conducted further investigations to understand the molecular mechanisms underlying the impact of RBM15 on the progression of CC. Our findings revealed a significant correlation between RBM15 and the JAK/STAT signaling pathway in CC. A previous bioinformatics study has shown that immune-related genes of CC are actively involved in the JAK-STAT pathway, which is consistent with our study , suggesting that the JAK-STAT signaling pathway plays an important role in CC progression. Moreover, the study has shown that HPV exploits the JAK/STAT signaling pathway to circumvent the immune system and stimulate cellular proliferation, thereby establishing viral persistence and propelling the progression of cancer . The progression of CC is impeded through the regulation of the JAK/STAT3 signaling pathway by silencing ITGB6 . And the involvement of STAT3 and its associated key genes, highlighting a significant pathogenic role of the JAK/STAT3 pathway in HPV-mediated CC . In our study, we found that knockdown of RBM15 significantly inhibited the protein expression of P-JAK2, JAK2, P-STAT3, STAT3, P-STAT5, and STAT5 in CC through further functional analysis. The evolutionary conservation of the JAK/STAT signaling pathway comprises essential constituents, namely ligand-receptor complexes, JAKs, and STATs . Upon receptor-ligand coupling, which leads to receptor dimerization/oligomerization, the associated JAK proteins phosphorylate themselves and trans-phosphorylate the tyrosine residues on the receptors, providing docking sites for the SH2 domains of the latent cytoplasmic STAT molecules . JAK2 is a non-receptor tyrosine kinase, which can be phosphorylated by members of both the gp130 receptor family and the class II cytokine-receptor family . Synergistic suppression of the survival and proliferation of CC cells is regulated through the combined inhibition of JAK1/2 and DNMT1 using newly identified small-molecule compounds . STAT3 and STAT5 are both members of the STAT family . STAT3 and STAT5 activity is frequently elevated aggressive subtypes of cancer, making them valuable prognostic indicators [49,50,51]. The pivotal role of STAT3 in tumor development stems from its capacity to regulate the transcription of genes involved in various key processes, such as the cell cycle (CCND1, CCNE1), metabolism (OCT1, HIF1A), and cell survival (BCL2, BCLXL, and HSP70) . Interference with STAT3 using siRNA leads to an enhancement in the sensitivity of CC cells to cisplatin . And in 90% of the CC patients, HPV infection is detected and shows a significant association with overexpression of STAT5 . Therefore, combined with the above findings, we hypothesized that RMB15 was likely to influence the development of CC by regulating the JAK-STAT signaling pathway, and that the RBM15/JAK2/STAT3/STAT5 axis hold promise as a potential therapeutic target for CC treatment.
Finally, the present study aimed to identify m6A-related DEGs that may contribute to the development of CC. And RBM15 has the potential to be developed as a diagnostic biomarker for CC. Knockdown of RBM15 in CC has been shown to suppress tumor cell proliferation, invasion, and migration as well as the JAK-STAT signaling pathway. Our study revealed an important correlation between the the JAK-STAT signaling pathway and m6A in CC.
All data in the manuscript is available through the responsible corresponding author.
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Ca-Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.
Gupta SM, Mania-Pramanik J. Molecular mechanisms in progression of HPV-associated cervical carcinogenesis. J Biomed Sci. 2019;26(1):28. https://doi.org/10.1186/s12929-019-0520-2.
Plummer M, Peto J, Franceschi S, International Collaboration of Epidemiological Studies of Cervical C. Time since first sexual intercourse and the risk of cervical cancer. Int J Cancer. 2012;130(11):2638–44. https://doi.org/10.1002/ijc.26250.
Hoover DS, Spears CA, Vidrine DJ, Walker JL, Shih Y-CT, Wetter DW. Smoking Cessation Treatment needs of low SES Cervical Cancer Survivors. Am J Health Behav. 2019;43(3):606–20. https://doi.org/10.5993/ajhb.43.3.14.
Balasubramaniam SD, Balakrishnan V, Oon CE, Kaur G. Key molecular events in Cervical Cancer Development. Med (Kaunas). 2019;55(7). https://doi.org/10.3390/medicina55070384.
Jain A, Ganesh B, Bobdey SC, Sathwara JA, Saoba S. Sociodemographic and Clinical Profile of Cervical Cancer Patients visiting in a Tertiary Care Hospital in India. Indian J Med Paediatr Oncol. 2017;38(3):291–5. https://doi.org/10.4103/ijmpo.ijmpo_20_16.
Rosen VM, Guerra I, McCormack M, Nogueira-Rodrigues A, Sasse A, Munk VC, Shang A. Systematic review and network Meta-analysis of Bevacizumab Plus First-Line Topotecan-Paclitaxel or Cisplatin-Paclitaxel Versus Non-Bevacizumab-Containing Therapies in Persistent, recurrent, or metastatic cervical Cancer. Int J Gynecol Cancer. 2017;27(6):1237–46. https://doi.org/10.1097/IGC.0000000000001000.
Desrosiers R, Friderici K, Rottman F. Identification of methylated nucleosides in messenger RNA from Novikoff hepatoma cells. Proc Natl Acad Sci U S A. 1974;71(10):3971–5. https://doi.org/10.1073/pnas.71.10.3971.
Zhang H, Shi X, Huang T, Zhao X, Chen W, Gu N, Zhang R. Dynamic landscape and evolution of m6A methylation in human. Nucleic Acids Res. 2020;48(11):6251–64. https://doi.org/10.1093/nar/gkaa347.
Garcias Morales D, Reyes JL. A birds’-eye view of the activity and specificity of themRNA m6Amethyltransferase complex. WIREs RNA. 2020;12(1). https://doi.org/10.1002/wrna.1618.
Wang X, Tian L, Li Y, Wang J, Yan B, Yang L, Li Q, Zhao R, Liu M, Wang P, Sun Y. RBM15 facilitates laryngeal squamous cell carcinoma progression by regulating TMBIM6 stability through IGF2BP3 dependent. J Exp Clin Cancer Res. 2021;40(1):80. https://doi.org/10.1186/s13046-021-01871-4.
Cai X, Chen Y, Man D, Yang B, Feng X, Zhang D, Chen J, Wu J. RBM15 promotes hepatocellular carcinoma progression by regulating N6-methyladenosine modification of YES1 mRNA in an IGF2BP1-dependent manner. Cell Death Discov. 2021;7(1):315. https://doi.org/10.1038/s41420-021-00703-w.
Zhang Z, Mei Y, Hou M. Knockdown RBM15 inhibits Colorectal Cancer Cell Proliferation and Metastasis Via N6-Methyladenosine (m6A) modification of MyD88 mRNA. Cancer Biother Radiopharm. 2022;37(10):976–86. https://doi.org/10.1089/cbr.2021.0226.
Darnell JE Jr. STATs and gene regulation. Science. 1997;277(5332):1630–5. https://doi.org/10.1126/science.277.5332.1630.
Sabaawy HE, Ryan BM, Khiabanian H, Pine SR. JAK/STAT of all trades: linking inflammation with cancer development, tumor progression and therapy resistance. Carcinogenesis. 2021;42(12):1411–9. https://doi.org/10.1093/carcin/bgab075.
Jin W. Role of JAK/STAT3 signaling in the regulation of Metastasis, the transition of Cancer Stem cells, and Chemoresistance of Cancer by epithelial-mesenchymal transition. Cells. 2020;9(1). https://doi.org/10.3390/cells9010217.
Gutierrez-Hoya A, Soto-Cruz I. Role of the JAK/STAT pathway in Cervical Cancer: its relationship with HPV E6/E7 oncoproteins. Cells. 2020;9(10). https://doi.org/10.3390/cells9102297.
Chen X, He H, Xiao Y, Hasim A, Yuan J, Ye M, Li X, Hao Y, Guo X. CXCL10 produced by HPV-Positive cervical Cancer cells stimulates exosomal PDL1 expression by fibroblasts via CXCR3 and JAK-STAT pathways. Front Oncol. 2021;11:629350. https://doi.org/10.3389/fonc.2021.629350.
Zhu X, Yang H, Zhang M, Wu X, Jiang L, Liu X, Lv K. YTHDC1-mediated VPS25 regulates cell cycle by targeting JAK-STAT signaling in human glioma cells. Cancer Cell Int. 2021;21(1):645. https://doi.org/10.1186/s12935-021-02304-0.
Wu S, He G, Liu S, Cao Y, Geng C, Pan H. Identification and validation of the N6-methyladenosine RNA methylation regulator ZC3H13 as a novel prognostic marker and potential target for hepatocellular carcinoma. Int J Med Sci. 2022;19(4):618–30. https://doi.org/10.7150/ijms.69645.
Han Q, Zhang X, Ren X, Hang Z, Yin Y, Wang Z, Chen H, Sun L, Tao J, Han Z, Tan R, Gu M, Ju X. Biological characteristics and predictive model of Biopsy-Proven Acute rejection (BPAR) after kidney transplantation: evidences of Multi-Omics Analysis. Front Genet. 2022;13:844709. https://doi.org/10.3389/fgene.2022.844709.
Liu Y, Wang J, Li L, Qin H, Wei Y, Zhang X, Ren X, Ding W, Shen X, Li G, Lu Z, Zhang D, Qin C, Tao L, Chen X. AC010973.2 promotes cell proliferation and is one of six stemness-related genes that predict overall survival of renal clear cell carcinoma. Sci Rep. 2022;12(1):4272. https://doi.org/10.1038/s41598-022-07070-1.
Jiang S, Ren X, Liu S, Lu Z, Xu A, Qin C, Wang Z. Integrated Analysis of the Prognosis-Associated RNA-Binding protein genes and candidate drugs in renal papillary cell carcinoma. Front Genet. 2021;12:627508. https://doi.org/10.3389/fgene.2021.627508.
Zhang D, Xu X, Wei Y, Chen X, Li G, Lu Z, Zhang X, Ren X, Wang S, Qin CJCcjotMCC. (2022) Prognostic role of DNA damage response genes mutations and their Association with the sensitivity of olaparib in prostate Cancer patients. 29:10732748221129451. https://doi.org/10.1177/10732748221129451.
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and Integrated Discovery. Genome Biol. 2003;4(5):P3.
Pei J, Li M, Wu C, Xu M, Shu T, Zhu C. Worldwide trends in cervical cancer incidence and mortality. Cancer. 2022;128(5):1141. https://doi.org/10.1002/cncr.34029.
Ran Y, Yan Z, Jiang B, Liang P. N6-methyladenosine functions and its role in skin cancer. Exp Dermatol. 2023;32(1):4–12. https://doi.org/10.1111/exd.14696.
Lan T, Li H, Zhang D, Xu L, Liu H, Hao X, Yan X, Liao H, Chen X, Xie K, Li J, Liao M, Huang J, Yuan K, Zeng Y, Wu H. KIAA1429 contributes to liver cancer progression through N6-methyladenosine-dependent post-transcriptional modification of GATA3. Mol Cancer. 2019;18(1):186. https://doi.org/10.1186/s12943-019-1106-z.
Liu T, Wei Q, Jin J, Luo Q, Liu Y, Yang Y, Cheng C, Li L, Pi J, Si Y, Xiao H, Li L, Rao S, Wang F, Yu J, Yu J, Zou D, Yi P. The m6A reader YTHDF1 promotes ovarian cancer progression via augmenting EIF3C translation. Nucleic Acids Res. 2020;48(7):3816–31. https://doi.org/10.1093/nar/gkaa048.
Li R, Song Y, Chen X, Chu M, Wang ZW, Zhu X. METTL3 increases cisplatin chemosensitivity of cervical cancer cells via downregulation of the activity of RAGE. Mol Ther Oncolytics. 2021;22:245–55. https://doi.org/10.1016/j.omto.2021.05.013.
Geng F, Fan MJ, Li J, Liang SM, Li CY, Li N. Knockdown of METTL14 inhibits the growth and invasion of cervical cancer. Transl Cancer Res. 2019;8(6):2307–15. https://doi.org/10.21037/tcr.2019.09.48.
Wang H, Luo Q, Kang J, Wei Q, Yang Y, Yang D, Liu X, Liu T, Yi P. YTHDF1 aggravates the progression of Cervical Cancer through m(6)A-Mediated Up-Regulation of RANBP2. Front Oncol. 2021;11:650383. https://doi.org/10.3389/fonc.2021.650383.
Condic M, Ralser DJ, Klumper N, Ellinger J, Qureischi M, Egger EK, Kristiansen G, Mustea A, Thiesler T. Comprehensive analysis of N6-Methyladenosine (m6A) Writers, Erasers, and readers in Cervical Cancer. Int J Mol Sci. 2022;23(13). https://doi.org/10.3390/ijms23137165.
Zhang C, Guo C, Li Y, Ouyang L, Zhao Q, Liu K. The role of YTH domain containing 2 in epigenetic modification and immune infiltration of pan-cancer. J Cell Mol Med. 2021;25(18):8615–27. https://doi.org/10.1111/jcmm.16818.
Hu M, Yang Y, Ji Z, Luo J. RBM15 Functions in Blood Diseases. Curr Cancer Drug Targets. 2016;16(7):579–85. https://doi.org/10.2174/1568009616666160112105706.
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(7620):369–73. https://doi.org/10.1038/nature19342.
Torres L, Lisboa S, Vieira J, Cerveira N, Santos J, Pinheiro M, Correia C, Bizarro S, Almeida M, Teixeira MR. Acute megakaryoblastic leukemia with a four-way variant translocation originating the RBM15-MKL1 fusion gene. Pediatr Blood Cancer. 2011;56(5):846–9. https://doi.org/10.1002/pbc.22765.
Zeng X, Chen K, Li L, Tian J, Ruan W, Hu Z, Peng D, Chen Z. Epigenetic activation of RBM15 promotes clear cell renal cell carcinoma growth, metastasis and macrophage infiltration by regulating the m6A modification of CXCL11. Free Radic Biol Med. 2022;184:135–47. https://doi.org/10.1016/j.freeradbiomed.2022.03.031.
Yang F, Liu Y, Xiao J, Li B, Chen Y, Hu A, Zeng J, Liu Z, Liu H. Circ-CTNNB1 drives aerobic glycolysis and osteosarcoma progression via m6A modification through interacting with RBM15. Cell Prolif. 2023;56(1):e13344. https://doi.org/10.1111/cpr.13344.
Wang Q, Vattai A, Vilsmaier T, Kaltofen T, Steger A, Mayr D, Mahner S, Jeschke U, Hildegard Heidegger H. Immunogenomic Identification for Predicting the prognosis of Cervical Cancer Patients. Int J Mol Sci. 2021;22(5). https://doi.org/10.3390/ijms22052442.
Morgan EL, Macdonald A. Manipulation of JAK/STAT signalling by high-risk HPVs: potential therapeutic targets for HPV-Associated Malignancies. Viruses. 2020;12(9). https://doi.org/10.3390/v12090977.
Zheng X, Zhu Y, Wang X, Hou Y, Fang Y. Silencing of ITGB6 inhibits the progression of cervical carcinoma via regulating JAK/STAT3 signaling pathway. Ann Transl Med. 2021;9(9):803. https://doi.org/10.21037/atm-21-1669.
Thakur K, Janjua D, Shishodia G, Chhokar A, Aggarwal N, Yadav J, Tripathi T, Chaudhary A, Senrung A, Bharti AC. Investigation of molecular mechanisms underlying JAK/STAT signaling pathway in HPV-induced cervical carcinogenesis using ‘omics’ approach. Med Oncol. 2022;39(12):255. https://doi.org/10.1007/s12032-022-01854-1.
Valentino L, Pierre J. JAK/STAT signal transduction: regulators and implication in hematological malignancies. Biochem Pharmacol. 2006;71(6):713–21. https://doi.org/10.1016/j.bcp.2005.12.017.
Verhoeven Y, Tilborghs S, Jacobs J, De Waele J, Quatannens D, Deben C, Prenen H, Pauwels P, Trinh XB, Wouters A, Smits ELJ, Lardon F, van Dam PA. The potential and controversy of targeting STAT family members in cancer. Semin Cancer Biol. 2020;60:41–56. https://doi.org/10.1016/j.semcancer.2019.10.002.
Schindler C, Strehlow I. Cytokines and STAT signaling. Adv Pharmacol. 2000;47:113–74. https://doi.org/10.1016/s1054-3589(08)60111-8.
She S, Zhao Y, Kang B, Chen C, Chen X, Zhang X, Chen W, Dan S, Wang H, Wang YJ, Zhao J. Combined inhibition of JAK1/2 and DNMT1 by newly identified small-molecule compounds synergistically suppresses the survival and proliferation of cervical cancer cells. Cell Death Dis. 2020;11(9):724. https://doi.org/10.1038/s41419-020-02934-8.
Polak KL, Chernosky NM, Smigiel JM, Tamagno I, Jackson MW. Balancing STAT activity as a therapeutic strategy. Cancers (Basel). 2019;11(11). https://doi.org/10.3390/cancers11111716.
Wagner KU, Rui H. Jak2/Stat5 signaling in mammogenesis, breast cancer initiation and progression. J Mammary Gland Biol Neoplasia. 2008;13(1):93–103. https://doi.org/10.1007/s10911-008-9062-z.
Ahonen TJ, Xie J, LeBaron MJ, Zhu J, Nurmi M, Alanen K, Rui H, Nevalainen MT. Inhibition of transcription factor Stat5 induces cell death of human prostate cancer cells. J Biol Chem. 2003;278(29):27287–92. https://doi.org/10.1074/jbc.M304307200.
Wu P, Wu D, Zhao L, Huang L, Shen G, Huang J, Chai Y. Prognostic role of STAT3 in solid tumors: a systematic review and meta-analysis. Oncotarget. 2016;7(15):19863–83. https://doi.org/10.18632/oncotarget.7887.
Carpenter RL, Lo HW. STAT3 target genes relevant to human cancers. Cancers (Basel). 2014;6(2):897–925. https://doi.org/10.3390/cancers6020897.
Huang LL, Rao W. SiRNA interfering STAT3 enhances DDP sensitivity in cervical cancer cells. Eur Rev Med Pharmacol Sci. 2018;22(13):4098–106. https://doi.org/10.26355/eurrev_201807_15401.
Sobti RC, Singh N, Hussain S, Suri V, Bharadwaj M, Das BC. Deregulation of STAT-5 isoforms in the development of HPV-mediated cervical carcinogenesis. J Recept Signal Transduct Res. 2010;30(3):178–88. https://doi.org/10.3109/10799891003786218.
Program Category: Frontier and Critical Technology Innovation
Project Category: Ganzhou City 2022 “Science and Technology + National Regional Medical Center” joint plan project – General project of key research and development plan
Department in charge: Agricultural and Social Development Department
Ethical approval and consent to participate
The animal experiments conformed to the Guide for the Care and Use of Laboratory Animals. Animal study has been approved by the Animal Ethics Committee of Ganzhou People’s Hospital.
All experiments were performed in accordance with relevant guidelines and regulations.
The manuscript reporting adheres to the ARRIVE guidelines for the reporting of animal experiments.
Consent to participate
Consent for publication
The authors declare no competing interests.
Conflict of interest
The authors declare that they have no conflicts of interest to disclose.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhang, C., Gu, L., Xiao, J. et al. Knockdown of RBM15 inhibits tumor progression and the JAK-STAT signaling pathway in cervical cancer. BMC Cancer 23, 684 (2023). https://doi.org/10.1186/s12885-023-11163-z