KIAA1429 expression in CML-BC patients is upregulated.
The expression of m6A related methyltransferase KIAA1429 is dysregulated in a variety of tumors (Fig. S1a) and is closely related to tumor progression and poor prognosis (Fig. S1b). Importantly, it is a key regulator of m6A modification (Fig. S1c). To explore the function of KIAA1429 in CML, the expression data of 16 cancers, documented in the Human Protein Atlas (HPA) database, were analyzed. The findings showed that in comparison with other cancers, the expression of KIAA1429 in CML cell line K562 was the highest (Fig. 1a). In addition, of all human tissues, bone marrow had the highest expression of KIAA1429 (Fig. S1d). Our team also analyzed the sequencing data in CML patients (GSE100026). Our findings showed that compared with CML-CP, the expression of m6A writer KIAA1429 in CML-BC was upregulated (Fig. 1b). The m6A modification level and KIAA1429 expression were measured in 40 patients with CML-CP and 18 patients with CML-BC included in this study. We found that the total m6A modification level (Fig. 1c) of KIAA1429 expression (Fig. 1d and e) (Fig. S1e) was significantly more upregulated in the CML-BC group than that in the CML-CP group, indicating that KIAA1429 could be associated with CML progression.
Imatinib is a first-line drug for the treatment of CML, and the progression of CML is generally associated with imatinib resistance. Therefore, the levels of m6A modification level and KIAA1429 expression were measured in the CML cell lines K562 and KCL22, as well as in the imatinib-resistant K562 cell line (K562/G01). The results obtained revealed that compared with the K562 and KCL22 cells, the total m6A modification level (Fig. 1f), as well as the mRNA and protein levels of KIAA1429 (Fig. 1g and h) were significantly upregulated in K562/G01 cells. In summary, KIAA1429 could upregulate the m6A modification to increase the imatinib resistance and promote tumor progression in CML.
KIAA1429 expression upregulation promoted the transformation of the malignant biological characteristics of CML cells.
To investigate the biological effects of KIAA1429 in CML, a KIAA1429 overexpression plasmid Lv-KIAA1429 was constructed and used to transfect CML cell lines (K562, K562/G01, and KCL22) (Fig. 2a and b). Our results showed that the total m6A modification level in the CML cells the of Lv-KIAA1429 group was significantly higher than that in the cells of the control group (Lv-NC) (Fig. 2c). KIAA1429 overexpression enhanced the proliferation and migration capabilities of CML cells (Fig. 2d, e, and g), while the cell apoptosis rate decreased (Fig. 2f). The Wright-Giemsa staining of cells before and after the overexpression revealed a lower number of primitive cells in the overexpressed cell lines than in the cells without KIAA1429 overexpression. In addition, the percentages of promyelocytes and myelocytes were increased, and the morphological characteristics were similar to the bone marrow characteristics of transformation from chronic phase to acute phase (Fig. 2h). The half-inhibitory concentration (IC50) assay showed that the resistance to Imatinib of K562/G01 cells was significantly upregulated (Fig. 2i), indicating that the increased expression of KIAA1429 could be associated with the imatinib resistance of CML cells. In summary, these findings showed that elevated expression of KIAA1429 reduced the sensitivity of CML cells to imatinib and inhibited the differentiation of CML cells, which could be a potential mechanism underlying the acute transformation of CML.
Knocking down KIAA1429 expression inhibited the transformation of malignant biological characteristics of CML cells.
The effect of KIAA1429 knockdown on CML progression was further investigated. First, a stable cell line of CML with KIAA1429 knockdown (Sh-KIAA1429) was constructed (Fig. 3a and b). In contrast to the overexpression, the KIAA1429 knockdown reduced the total level of RNA m6A modification in CML cells (Fig. 3c). CML cell proliferation was inhibited (Fig. 3d and e), cell apoptosis was increased (Fig. 3f), and cell migration was suppressed (Fig. 3g). Cell morphological observation showed that the KIAA1429 knockdown increased the percentage of lobulated or stab-nuclear cells, and the cells were stabilized in the chronic phase (Fig. 3h). Moreover, KIAA1429 knockdown significantly inhibited the Imatinib resistance of K562/G01 cells (Fig. 3i). These findings demonstrated that KIAA1429 knockdown could inhibit the malignant transformation of biological characteristics of CML.
RAB27B is the downstream target gene of KIAA1429
KIAA1429 is the critical enzyme of m6A methylation modification, the dysregulation of which induces the activation of relevant cancer pathways. In this study, the KIAA1429 expression was knocked down by shRNA in K562 cells, and K562-Sh1 and K562 cells were subjected to mRNA-seq analysis. Screening was then performed using |logFC| > 1 and P-value < 0.05 as the criteria, which retrieved 1587 differential expression genes (DEG) (Fig. 4a) (Fig. S2a and S2b). To clarify the association between DEGs and m6A modification, the KIAA1429 related MeRIP-seq data in the m6a2target database (http://m6a2target.canceromics.org/#/perturbation) was differentially analyzed (|logFC| > 1 and P-value < 0.05), and retrieved 1499 m6A related DEGs. OVERLAP analysis was performed for the two groups of DEGs, and 120 co-expressed DEGs (cDEGs) were identified (Table S1, Fig. 4b). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed for cDEGs, which showed the biological functions of the genes were mainly enriched in various aspects, including enzyme activity inhibitor, focal adhesion, and exocytosis (Fig. 4c). Network analysis of protein-protein interaction (PPI) (https://cn.string-db.org/) showed that RAB27B is the Hub gene of the network (Fig. S2c), indicating that RAB27B could be a potential key target gene regulated by KIAA1429. RT-qPCR analysis showed that the expression of RAB27B mRNA in K562/G01 cells was significantly higher than in K562 and KCL22 cells (Fig. 4d); while in K562 cells with KIAA1429 knockdown, RAB27B level also decreased with the reduction of KIAA1429 level (Fig. 4e and 4f). In clinical samples of CML, the expression of RAB27B mRNA was significantly higher in the CML-BC group than in the CML-CP group (Fig. 4g), which was significantly positively associated with KIAA1429 (r = 0.338, P = 0.010) (Fig. 4h). These findings evidenced that RAB27B could be an important downstream target gene for KIAA1429.
KIAA1429 regulates RAB27B mRNA stability in m6A-dependent manner and could be recognized by YTHDF1
To verify our speculation, RIP-qPCR was utilized to assess the relationship between KIAA1429 and RAB27B. After enrichment by KIAA1429-antibody, the expression of RAB27B mRNA increased significantly, indicating that RAB27B is a downstream mRNA binding to KIAA1429 (Fig. 4j). After the knockdown of KIAA1429, CML cells were treated by actinomycin D. We found that the RAB27B mRNA expression level decreased gradually over time (Fig. 4i), indicating that KIAA1429 regulated the stability of RAB27B mRNA and consequently influenced the expression and function of RAB27B.
The m6A methylation sites of RAB27B mRNA were predicted using the SRAMP website (http://www.cuilab.cn/sramp) (Fig. 4k), which showed that five m6A modification sites in the sequences of RAB27B mRNA, i.e. A116, A1466, A1508, A1518, A2300, A2981, A3478, A4641, had relatively high possibility of modification (Fig. 4k)(Fig. S2e). The SELECT assay was performed to measure the m6A modification levels at the sites before and after knockdown, which showed that after KIAA1429 knockdown, the m6A modifications at A1116, A2981, and A3478 sites decreased significantly (Fig. 4l and m) (Fig. S2f), while the modification at other sites did not change significantly, indicating that KIAA1429 knockdown downregulates the m6A modification level at these sites.
Previous studies have revealed that m6A reader proteins could recognize the information of m6A modification in mRNA of transcripts and consequently regulate the expression of target mRNA, influencing the biological processes in cancer. To further explore the molecular mechanisms of RAB27B m6A modification, the relevant reader proteins were searched by bioinformatics measurements. Of the various m6A reader proteins, YTHDF1 was positively correlated with the levels of RAB27B mRNA (GEO: GSE100026) (Fig. 4n), indicating that YTHDF1 could be a key reader protein of RAB27B.
This finding was further confirmed by in vitro experiments. The K562 sh-YTHDF1 cells were successfully constructed by using a lentivirus (Fig. 4p). Our results showed that after knockdown of YTHDF1 in cells, the mRNA and protein levels of RAB27B also decreased correspondingly (Fig. 4q and 4r), and the stability of RAB27B mRNA also decreased (Fig. 4o). RIP-qPCR findings showed that RAB27B mRNA could be abundantly enriched by YTHDF1 protein, and thus influence the stability of RAB27B mRNA (Fig. 4s). These findings indicated that YTHDF1 could bind to RAB27B mRNA in CML cells.
Inhibiting RAB27B expression rescued the promotion effects of KIAA1429 on CML cell proliferation and drug resistance
Previous studies have shown that RAB27B could regulate cell signaling transduction and exosome release, and thus influence tumor proliferation and drug resistance. However, the role of RAB27B in CML is still unclear. To clarify the biological functions of RAB27B, the lentivirus of RAB27B knockdown was used to transfect K562 and K562/G01 cells, as well as KIAA1429 overexpression cells (Fig 5a and 5b). The results we obtained revealed that after the knockdown of RAB27B, the proliferation capability of K562, K562/G01, and KIAA1429 overexpression cells decreased (Fig. 5c and 5d), while KIAA1429 overexpression partially rescued the inhibitory effects of RAB27B knockdown on cell proliferation capability. Cells were arrested in the S phase after RAB27B knockdown (Fig. 5e). Pathway array results showed that the expression levels of genes associated with S cell phase, such as ABL1 and MCM2, were lower, whereas those of negative regulation genes, such as CDKN2B and TP53, were higher (Fig. 5f). RAB27B knockdown induced the upregulation of cell apoptosis and suppressed the inhibitory effects of KIAA1429 overexpression on CML cell apoptosis (Fig. 5g). In addition, RAB27B knockdown increased the sensitivity of K562/G01 cells to imatinib (Fig. 5h), while KIAA1429 overexpression enhanced the resistance of the cells with RAB27B knockdown to the drug. To verify whether the drug resistance of CML cells was associated with RAB27B mediated drug transportation, K562/G01 cells in Sh-NC group and sh-RAB27B group were treated by 2 µmol/L imatinib for 48 h, and then culture medium with no imatinib was utilized to further treat cells for 24 h. Afterwards, LC-MS was used to measure the imatinib concentration in the cells and culture medium. As can be seen in Fig. 5i, the level of imatinib in the culture medium was lower, whereas its level in the cells in the Sh-RNA27B group was higher. The obtained findings also showed that the RAB27B knockdown reduced the number of exosomes in the culture medium (Fig. S3). These findings indicated that RAB27B knockdown inhibited the efflux of imatinib in cells, leading to drug enrichment in the cells, which could be associated with the inhibition of the exosome secretion.
Inhibiting KIAA1429 axis inhibited the tumorigenesis capability of CML cell lines in vivo
In vivo experiments were further performed. The K562 cells in the negative-control group, as well as K562 cells with KIAA1429, RAB27B, or YTHDF1 knockdown, were inoculated into BALB/C-nu nude mice. The results showed that the sizes of the tumors in the nude mice inoculated with KIAA1429, RAB27B, or YTHDF1 knockdown cells were smaller than those in the nude mice in the control group (Fig. 5j). Additionally, the tumor growth rate in the former group was lower (Fig. 5k), indicating that the KIAA1429, RAB27B, or YTHDF1 knockdown significantly reduced the tumorigenesis capability of K562 cells. The body weight curve of the nude mice was also used to assess their overall nutrition status. We found that the body weight of the nude mice in the control group was reduced, whereas that of the nude mice in the KIAA1429, RAB27B, and YTHDF1 groups remained stable (Fig. 5l), which was significantly improved compared with that in the control group.
Rucaparib inhibited the KIAA1429 expression and the proliferation and drug resistance of CML cells
Our previous studies have revealed the role of KIAA1429 in promoting the progression of CML. We also aimed to treat CML by targeting this gene. To screen the treatment drugs, the Genomics of Drug Sensitivity in Cancer (GDSC) database and mRNA-seq data were analyzed comprehensively, and IC50 was used to assess the sensitivity of the cells in each group to chemotherapy drugs. The results obtained (Fig. 6a) showed that the KIAA1429 knockdown group had higher IC50 values than those in the ponatinib (P = 0.002), rucaparib (P = 0.026), axitinib group (P = 0.026), and ATRA groups (P = 0.026). These findings indicated that cells with KIAA1429 overexpression are more sensitive to these anti-tumor drugs. Ponatinib, axitinib, and ATRA have already been applied for CML treatment in clinical practice. Thus, in this study, we focused on the treatment with rucaparib.
Interestingly, the 48-h treatment with concentration gradients of rucaparib decreased the cell expression of KIAA1429 mRNA and protein (Fig. 6b and c). The degrees of this reduction were positively associated with drug concentrations, indicating that rucaparib could inhibit the expression of KIAA1429. In addition, rucaparib also suppressed the proliferation of CML cells (Fig. 6d), while the IC50 values of K562 and K562/G01 were 101.3 μM and 87.81 μM (Fig. 6e), respectively. The 48-h rucaparib treatment dose-dependently promoted the apoptosis of CML cells (Fig. 6f). These findings demonstrated that rucaparib suppressed the proliferation and promoted the apoptosis of CML cells. These effects were more pronounced in imatinib-resistant CML cells than in sensitive cells. We further investigated whether rucaparib could enhance the drug resistant CML cells to imatinib. Low concentrations of rucaparib and imatinib of approximately 50% of IC50 were used to treat cells independently or in combination, and then the proliferation and apoptosis of cells were observed. We established that the proliferation of cells in the rucaparib + imatinib group was significantly lower than that in the control, rucaparib, and imatinib groups (Fig. 6g). The apoptosis rate in the rucaparib + imatinib group was significantly higher than those in the other three groups (Fig. 6h). Rucaparib treatment decreased significantly the IC50 of imatinib (Fig. 6i). These findings indicated that the treatment with rucaparib significantly increased the sensitivity of CML cell lines to imatinib.
We further investigated whether rucaparib could be used to treat CML in vivo. Firstly, toxicological experiments in BALB/C-nu mice were performed for verification. The body weight (Fig. 6j), blood routine parameters (Fig. 6k), and biochemical indicators such as ALT and UREA (Fig. 6l) in the treatment group were not significantly different from those in the control group, indicating that this dose of rucaparib exerts no liver or renal toxicity in vivo. No tissue injury or inflammatory lesions were found in the treatment group (Fig. 6m), indicating the high biocompatibility and safety of rucaparib at this dose.
The nude mice models of subcutaneous tumorigenesis models of CML were treated by rucaparib to assess its therapeutic effects. The tumor size (Fig. 6n), weight (Fig. 6o), and growth weight (Fig. 6p) in the treatment group were significantly lower than those in the control group, indicating that rucaparib inhibited the tumorigenesis in the CML cells in vivo, suggesting its potential therapeutic value.