Characterization of the Newly Established Homoharringtonine- (HHT-) Resistant Cell Lines and Mechanisms of Resistance

Li, Fenglin; Ling, Qing; Hu, Chao; Wang, Huafeng; Ye, Wenle; Li, Xia; Zhang, Xiang; Lin, Xiangjie; Wei, Wenwen; Huang, Xin; Qian, Yu; Zhuang, Haihui; Jin, Jie; Lu, Ying

doi:https://doi.org/10.1155/2022/2813938

Journal of Oncology

On this page

Abstract Introduction Methods Results Discussion Conclusion Abbreviations Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 2813938 | https://doi.org/10.1155/2022/2813938

Characterization of the Newly Established Homoharringtonine- (HHT-) Resistant Cell Lines and Mechanisms of Resistance

Fenglin Li,^1,2Qing Ling,²Chao Hu,^2,3Huafeng Wang,^2,3Wenle Ye,²Xia Li,^2,3Xiang Zhang,^2,3Xiangjie Lin,^2,3Wenwen Wei,²Xin Huang,^2,3Yu Qian,²and Haihui Zhuang¹ et al.

Academic Editor: Subhadeep Roy

Received20 Feb 2022

Revised22 Jun 2022

Accepted07 Jul 2022

Published30 Aug 2022

Abstract

Homoharringtonine- (HHT-) based HHT, aclarubicin, and cytarabine (HAA) induction regimen is the first-line therapy for nonelder acute myeloid leukemia (AML) patients in China. However, drug resistance is a new challenge, and little attention has been devoted to excavating resistant mechanisms. This study used the classic method to construct six HHT-resistant cell lines with a gradually increasing resistance index (RI) to discover HHT drug resistance mechanisms dynamically. After HHT resistance, the cell growth rate decreased, cell cycle delayed, and P-glycoprotein (p-gp, CD243) expression levels increased. Furthermore, we explored the changes in transcriptomics between HHT-sensitive and HHT-resistant cells using RNA-sequence. Through Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and hub gene analyses, we found that immune activity, especially G-protein coupled receptor (GPR) and related molecules, may mediate HHT resistance. Moreover, Calcitonin Receptor-Like (CALCRL) and G Protein Subunit Alpha I1 (GNAI1), which belong to GPRs, were stimulated in HHT-resistant cell strains in vitro and vivo, indicating that they may play a critical role in HHT resistance. In addition, these two genes have prognostic significance for AML patients. Taken together, we successfully constructed HHT-resistant cell lines with dynamic RIs and explored the resistance mechanisms, which will help identify new drugs for HHT-resistant AML patients.

1. Introduction

Acute myeloid leukemia (AML) is the most common form of acute leukemia in adults and is characterized by stem cell transcriptional dysregulation and leads to a block in the differentiation [1]. Nowadays, advances in molecular biology and genetics have spurred advancements in diagnosis and prognosis in AML (2). However, the traditional treatment for AML, the DA or IA program, has not changed significantly over the past 40 years and is accompanied by a poor prognosis. Although many researchers have attempted to optimize traditional chemotherapy, the overall survival is still not ideal [2–4]. The main reasons for poor prognosis in AML patients are primary drug resistance and disease relapse. Moreover, the presence of drug-resistant leukemia cells is the key cause of relapse [5]. Therefore, characterizing drug-resistant leukemia cells and exploring resistance mechanisms are effective ways to improve AML prognosis.

Homoharringtonine (HHT) is a cytotoxic alkaloid that was initially extracted from Cephalotaxus hainanensis, and it has significant antileukemic effects in AML, chronic myeloid leukemia (CML), and myelodysplastic syndrome (MDS) [6]. One phase 3 clinical study showed that the complete remission (CR) (73% versus 61%, ) and the 3-year event-free survival (EFS) (35.4% versus 23.1%, ) in the HAA (HHT plus aclarubicin and cytarabine) treatment group were better than those in the DA (daunorubicin plus cytarabine) group [7]. The results suggest that the HAA regimen could be an alternative induction treatment for untreated AML patients, particularly those with favorable and intermediate cytogenetics. According to these results, HHT-based “HAA” program is the first-line therapy for nonelder AML patients in China [8]. The effectiveness of HAD (HHT + cytarabine + daunorubicin) program was also confirmed by another team in China [9]. However, the RFS of patients with HAA regimen was similar to that of DA regimen, and drug resistance has also become a challenge with the widespread use of HHT in AML. Exploring the resistance mechanism and related targets is crucial for refractory and relapsed AML patients who rely on HHT-based treatment.

Studies on the mechanisms of cytotoxicity and resistance to chemotherapy have been based on the development and analysis of resistant cancer cell lines [10]. For hematological malignancies, drug-resistant cell lines have been successfully constructed by increasing the concentration of a specific drug. The first successfully constructed drug-resistant cell line in hematology was K562 in 1983, which was exposed to increasing vincristine doses from 3 nM to 60 nM (12). Subsequently, drug-resistant cell lines for mantle cell lymphoma [11], multiple myeloma [12], and AML [13–15] have been successfully constructed.

Based on the above-mentioned background, we aimed to establish HHT-resistance AML cell lines with different RIs via exposure to gradually increasing HHT concentrations. After 15 months cultivated with HHT, 6 HHT-resistant cell lines were stabilization-established. Then cell characteristics including cell growth, apoptosis, cell cycle, and antigen p-gp were investigated to explore the biological behaviors difference between resistant cells and sensitive cells. From RNA-sequence data, we found that immune activity, especially G-protein coupled receptor and related molecules, such as CALCRL and GNAI1, may mediate HHT resistance. However, we have not further explored how CALCRL and GNAI1 mediate HHT resistance and the mechanism of resistance, and further exploration is needed in the future. The role of HHT in AML drug resistance is also explored from metabolomics or proteomics in the future.

2. Methods

2.1. Cell Culture

AML cell lines THP-1 (CVCL_0006), HL-60 (CVCL_A794), NB4 (CVCL-0005), Kasumi-1 (CVCL_0589), MOLM-13 (CVCL_2119), MV4-11 (CVCL_0064), KG-1 (CVCL_0374), and U937 (CVCL_0007) were purchased from American Type Culture Collection (ATCC), and OCI-AML3 (CVCL_1844) was purchased from Deutsche Sammlung von Mikroorganismen und Zellkulturen (DSMZ). THP-1, HL-60, OCI-AML3, NB4, Kasumi-1, KG-1α, and U937 were cultured in RPMI-1640 medium (HyClone, UT, USA); MOLM-13 and MV4-11 were cultured in IMDM medium (Corning Cellgro, USA), and supplemented with 10% fetal bovine (Thermo Fisher Scientific, Gibco, USA) at 37°C with 5% CO₂. All cells were passaged less than 20 times and had cell identification certificates.

2.2. Cell Proliferation Assay

Proliferation and viability of AML cells were measured by CellTiter 96® AQueous One Solution Cell Proliferation Assay (Promega, USA). Briefly, THP-1, HL-60, OCI-AML3, NB4, Kasumi-1, KG-1, U937, MOLM-13, MV4-11, and resistant strains were inoculated into 24 wells at a final concentration of 1 × 10⁵ ml and treated with the indicated concentrations of HHT for 48 hours. Then 100 μL cell supernatant was transferred into 96-well plates and 10 μL CellTiter 96® Aqueous One Solution Cell Proliferation Assay was added per well. Absorbance was measured by a microplate reader (Thermo Scientific Varioskan® Flash, USA) at 490 nm to determine the OD value after 4 hours. The GraphPad software (version 6) calculated the cell proliferation rate and half-maximal inhibitory concentration (IC₅₀) value according to the OD.

2.3. Construction of HHT-Resistant Cells with Dynamic Resistance Index (RI)

MV4-11 and MOLM-13 cells were treated with gradually increasing HHT concentrations (1 nmol/L to 50 nmol/L). The concentration of HHT was increased at approximately 2-week to 4-week intervals and maintained until cells grew normally in a complete growth medium in the presence of different concentrations of HHT for two weeks. Cells that grew normally in medium with 10% fetal bovine containing 10 nmol/L, 30 nmol/L, and 50 nmol/L HHT were preserved. These cells were designated as MV4-11 R10, MV4-11 R30, MV4-11 R50, MOLM-13 R10, MOLM-13 R30, and MOLM-13 R50. Then, the IC₅₀ value of those cells was calculated and the RI of those cells was calculated with the following formula: IC₅₀ (resistance)/IC₅₀ (sensitive).

2.4. Flow Cytometry Analysis

To examine the cell apoptosis, MV4-11 and MOLM-13 sensitive and resistant cells were inoculated into 24 wells at a final density of 1 × 10⁵/ml. Then cells were collected and washed twice with 1× phosphate-buffered saline (1× PBS) and inoculated into 300 μL 1× binding buffer with 10 μL Propidium (PI) and 5 μL Annexin-V for 30 minutes (Cat: AP101, MULTISCIENCES, Hangzhou, China). For cell cycle analysis, cells were seeded in 24 wells at a final concentration of 3 × 10⁵/ml cells per well and collected and fixed with precooled 75% ethanol at 4°C overnight. The next day, cells were washed twice with 1× PBS. Flow cytometry was utilized to detect the DNA content after PI (Cat: CCS012, MULTISCIENCES, Hangzhou, China) staining for about 30 minutes. For p-gp expression analysis, cells were inoculated into 300 μL of 1× PBS with 10 μL of p-gp antibody (RRID: AB_396548, BD, USA) labeled with APC for 30 minutes. Flow cytometry was performed and analyzed on ACEA NovoCyte (ACEA, USA).

2.5. Mouse Xenograft Model

Male BALB/c nude mice (5 weeks old) were obtained from Shanghai Experimental Animal Center of the Chinese Academy of Sciences (Shanghai, China). After one week of acclimatization, mice were subcutaneously injected with 1 × 10⁶ MV4-11 and MV4-11 R50 cells (suspended in 100 μl 1× PBS), respectively. Tumor volumes were evaluated using the following formula: 1/2 × length × width². After 7∼11 days, tumor volume reached 100–200 mm³, mice were randomly divided into MV4-11 group (n = 4) and MV4-11R group (n = 5), and HHT was administered by intraperitoneal injection every day at 0.5 mg/kg body weight for 7 days. Tumor volumes were measured twice daily using the caliper. All animal experiments were reviewed and approved by the Ethics Committee for Laboratory Animals of the First Affiliated Hospital, College of Medicine, Zhejiang University (Hangzhou, China) (reference number: 20201522).

2.6. Clinical Patients

178 AML patients between March 2010 and June 2014 with detailed diagnosis and treatment information were included in this study. Patients with acute promyelocytic leukemia (APL) and patients undergoing bone marrow transplantation were excluded. The clinical characteristics of the patients are summarized in Table 1. WHO classification, conventional cytogenetic banding assay, and molecular analyses were performed centrally as previously described in AML diagnosis [17]. Chromosomal abnormalities and gene mutations of NPM1, FLT3-ITD, CEBPA, DNMT3A, IDH1, and IDH2 were conducted using the previously described methods [18]. We performed the above cytogenetic analyses blinded to both the levels of CALCRL expression and clinical outcome; the clinical characteristics of the patients are summarized in Table 1. All patients provided written informed consent to participate in the study. The Research Ethics Committee of the First Affiliated Hospital, College of Medicine, Zhejiang University, approved this study (reference number: 20201280).

2.6.1. Quantitative Reverse Transcriptase-PCR

Total RNA was isolated from cells using the TRIzol reagent (Cat: 9109, Invitrogen, USA). cDNA was prepared using the reverse transcription kit following the instructions (Thermo, USA). Real-time quantitative PCR (RT-qPCR) was performed in duplicate with the fluorescent DNA binding dye SYBR green (Cat: RR420, Takara, JPN). PCR reactions were performed in a total volume of 25 μl containing 1 μl of 100 ng/μl sample cDNA, 12.5 μl 2× PCR mix, 1 μl 0.5 μM of each primer, and 10.5 μl ddH2O. mRNA levels were normalized to GAPDH. The primers used for quantitative PCR are presented in Table 2.

2.7. Western Blot

Total protein was extracted using RIPA buffer (REF: 89900, Sigma-Aldrich, USA) supplemented with protease inhibitor and phosphatase inhibitor cocktail (78420, Thermo Fisher, USA). The protein content was quantified using the BCA protein assay kit (23225, Thermo Fisher). An estimated 40 μg protein was loaded per well on 4–12% SDS-PAGE gel and transferred onto a PVDF membrane (Millipore, Billerica, MA, USA) preactivated with methanol. Membranes were blocked with 5% milk and incubated with antibodies against specific proteins. After overnight, membranes were washed with 1× TBST and incubated with horseradish peroxidase labeled secondary antibodies. The protein bands were visualized using enhanced chemiluminescence reagent (BIO-RAD Gel Dox TM XR+, USA). The antibody of GNAI1 12617-1-AP was purchased from Proteintech (Rosemont, USA), CALCRL (ab84467) was purchased from Abcam (Cambridge, MA, USA), and GAPDH ^#5174S and the Rabbit (DA1E) mAb IgG secondary antibody were purchased from CST (Danvers, MA) (^#3900).

2.7.1. RNA-Sequence Data Analysis

The mRNA expression profiles of MV4-11 sensitive cells and MV4-11 resistant strains (MV4-11 R10/R30/R50) were obtained using high-throughput sequencing (RNA-sequence). The sequenced data were obtained from Illumina HiSeq X Ten platforms, and the data were uploaded to NCBI’s SRA and access to cite for these SRA data: PRJNA664675. For RNA-seq data, mRNA expression levels were calculated as RPKM (Reads Per Kilo-base per Million reads). We used DE Seq (1.18.0) R package to analyze the differentially expressed mRNAs.

2.8. Statistical Analysis

2.8.1. Clinical Data Analysis

Overall survival (OS) was defined as the time from diagnosis to death due to any cause. Event-free survival (EFS) was defined as the time from diagnosis to removal from the study because of the absence of CR, relapse, or death. In the association between CALCRL expression and OS, EFS was estimated by the Kaplan-Meier method and log-rank test. In prognostic analysis of other genes like GNAI1, CALCRL was calculated via online database “GEPIA.”

Experimental analysis, unpaired two-tailed Student’s t-test, or ANOVA and Chi-square test were performed using IBM SPSS Statistics 20 software or GraphPad Prism 5. was considered statistically significant. , , and .

2.9. Hub Gene Analysis

Protein-protein interaction (PPI) network was created by online database “STRING” [19] and visualized by Cytoscape software [20]. Top 10 hub genes were analyzed using cytoHubba APP (22). Biological processes, molecule functions, and cell type analysis were performed using the “Enrichment” online database.

3. Results

3.1. MV4-11 and MOLM-13 Were the Most Sensitive AML Cell Lines to HHT

To explore the cytotoxicity of HHT on AML cells, we detected the antileukemic effect of HHT on nine AML cell lines with different genetic backgrounds via IC₅₀ value measurement. The IC₅₀ values of MOLM-13 and MV4-11 were 6.858 nM and 7.207 nM, respectively, while other cell lines exceeded 10 nM (Figure 1(a)). The inhibitory effect of HHT on FLT3-ITD mutant cell lines MOLM-13 and MV4-11 was significantly higher than that on the nonmutant cell lines (Figure 1(b)).

(a)

(b)

3.2. Establishment and Verification of the HHT-Resistant Cell Lines

As MV4-11 and MOLM-13 were the most sensitive AML cell lines to HHT, these two cell lines were selected to construct HHT-resistant cell lines by exposing them to gradually increasing HHT concentrations for 15 months. Six resistant strains were successfully built, and HHT concentrations were induced at 10 nM, 30 nM, and 50 nM, respectively. The RIs of the three MOLM-13 resistant cell lines (R10, R30, and R50) were 7.13 nM, 23.75 nM, and 109.9 nM, respectively, which were 1.82-fold, 6.07-fold, and 28.11-fold higher than those of the sensitive cells. Contemporarily, the RIs of MV4-11 resistance cell lines (R10, R30, and R50) were 18.31 nM, 40.67 nM, and 58.82 nM, respectively, which were 4.92-fold, 13.08-fold, and 15.81-fold higher than those of sensitive cells (Figure 2(a)). Furthermore, there was no obvious apoptosis of HHT-sensitive cells and HHT-resistant cells under natural culture conditions (Figure S1A). To verify that HHT had no cytotoxic effects on these six cell lines, we evaluated cell growth in vitro and in vivo. The growth rate of MV4-11 R50 and MOLM-13 R50 was normally or partly inhibited with 10 nM HHT treatment, whereas sensitive cells (MV4-11S and MOLM-13S) were utterly suppressed (Figure 2(b)). In addition, the growth rate of HHT-resistant cells was lower than that of the HHT-sensitive cells (Figure 2(b)). Then we found that G0/G1 phase was delayed, and S phase replication was gradually declined in both MOLM-13 and MV4-11 cells with a progressively increasing HHT RI (Figure 2(c)). Additionally, results from flow cytometry showed that original HHT-resistant strains and constructed resistant cell lines had a higher expression of multidrug resistance protein p-gp (Figure 2(d)).

(a)

(b)

(c)

(d)

Next, we used the nude mouse subcutaneous tumor-bearing model to observe the drug resistance of HHT-resistant strains in vivo. Subcutaneous tumor mass in mice is mainly infiltrated by leukemia cells (Figure S1B). The animal experimental results also showed that tumor volume, growth rate, and weight were not significantly decreased after HHT treatment in MV4-11 R50 mice, while they were significantly decreased in MV4-11 sensitive mice (Figures 3(a)-3(b)). After HHT treatment, tumor volume was reduced more significantly in the MV4-11 sensitive group than in the resistant group (Figure S1C). These results revealed six HHT-resistant cell lines with different RIs were successfully established.

(a)

(b)

3.3. Differentially Expressed Genes (DEGs) and Pathways between Sensitive and Resistant Strains

As the HHT-resistant cell lines were successfully established, we aimed to explore the dynamic changes in gene expression with gradually increasing HHT RIs. Gene expression profiles between different RIs of HHT-sensitive and HHT-resistant cells were analyzed by RNA-sequence. There were 342, 485, and 489 upregulated genes and 338, 450, and 668 downregulated genes in the three sets MV4-11 R10 versus MV4-11 S, MV4-11 R30 versus MV4-11 S, and MV4-11 R50 versus MV4-11 S, respectively (Figure 4(a)). The RT-qPCR results from three upregulated genes and three downregulated genes verified the reliability of RNA-sequence (Figure S2). Further intersection analyses showed that 227 genes were all upregulated, and 243 were downregulated in three MV4-11 resistant strains by Venn methods (Figure S3).

(a)

(b)

(c)

(d)

To better understand the biological functions of these DEGs, we performed enrichment analysis. Upregulated genes in resistant groups took part in several stable pathways related to the immune activity, such as neutrophil activation involved in immune response and neutrophil-mediated immunity (Figures 4(b)–4(d)). The downregulated genes participated in several metabolic and biosynthetic processes, and the differences increased as the resistance index increased (Figures S4A–S4C).

3.4. Key Regulated DEGs in HHT-Resistant Cells

To explore the critical proteins in HHT resistance, PPI nodes from the online database “STRING” were obtained to discover the interactions between proteins. In addition, the top 10 hub genes were analyzed using cytoHubba software from both upregulated and downregulated DEGs among three sets; the top 10 hub genes among upregulated DEGs are shown in Figures 5(a)–5(c) with increasing RI, and the top 10 hub genes among downregulated DEGs are shown in Figures 5(d)–5(f). From biological function (BP) analysis, upregulated hub genes in MV4-11 resistant strains took part in dendritic cell chemotaxis; G-protein coupled receptor signaling pathway, coupled to cyclic nucleotide second messenger; positive regulation of T-helper 1 type immune response; and so forth (Figure S5). In those processes, the biological processes “G-protein coupled receptor (GPR) signaling pathway, coupled to cyclic nucleotide second messenger” were significantly upregulated as the RI increased (Figure S5). These results indicated that immune biology process, especially GPR signaling, may mediate HHT resistance. The biological processes of downregulated genes included a variety of biological processes, while the metabolic and biosynthetic processes were not stable as the resistance index decreased (Figure S6).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5

The hub genes, respectively, in both upregulated and downregulated differential genes among three sets. (a) The top 10 hub genes among upregulated genes in MV4-11 R10 versus MV4-11 sets. (b) The top 10 hub genes among upregulated genes in MV4-11 R30 versus MV4-11 sets. (c) The top 10 hub genes among upregulated genes in MV4-11 R50 versus MV4-11 sets. (d) The top 10 hub genes among downregulated genes in MV4-11 R10 versus MV4-11 sets. (e) The top 10 hub genes among downregulated genes in MV4-11 R30 versus MV4-11 sets. (f) The top 10 hub genes among downregulated genes in MV4-11 R50 versus MV4-11 sets.

3.5. CALCRL and GNAI1 May Mediate HHT Resistance and Have Clinical Significance in AML

Given that immune biological processes, especially GPR signaling, may mediate HHT resistance, we further focused on four upregulated DEGs (GPR183, CNR2, CALCR, and GNAI1), which belong to the G-protein family. First, the mRNA expression levels of these genes were detected in the HHT innate resistant cell lines THP-1 and KG-1 and sensitive cell lines MV4-11 and MOLM-13. Only CALCRL and GNAI1 had higher expression levels in innate resistant cells THP-1and KG-1 (Figure 6(a)). Contemporarily, the protein levels of CALCRL and GNAI1 were also higher in innate resistant cells and mice tumor (Figure 6(b)). In addition, compared to the MV4-11 sensitive mice group, the MV4-11 R50 mice had higher protein levels of CALCRL and GNAI1 (Figure 6(b)). Then, we explored the clinical significance of CARCRL and GNAI1 in AML. First, higher expression levels of CALCRL and GNAI1 were found in AML cells than in the normal cells (Figure 6(c)). Moreover, patients with high GNAI1 expression tended to have a poor prognosis according to the public database GEPIA (Figure 6(c)). Patients with high CALCRL expression had also worse OS from GEPIA (Figure 6(c)). We also verified the clinical significance of CALCRL by our database, and the results showed that AML patients with high CALCRL expression had poor OS and EFS (Figure 6(d)). Interestingly, CALCRL and GNAI1 had significant positive correlation by analyzing the TCGA database (Figure 6(e)). These results indicated that CALCRL and GNAI1 were highly expressed in HHT-resistant cells and may mediate HHT resistance and might have interaction.

(a)

(b)

(c)

(d)

(e)

Figure 6

Verification of the hub genes (CALCRL and GNAI1) and the prognostic significance of CALCRL and GNAI1 for AML. (a) The mRNA level of GPR183, CNR2, CALCRL, and GNAI1 in four AML cell lines by RT-qPCR. (b) The protein level of CALCRL and GNAI1 in four AML cell lines and the mice tumor mass from HHT-sensitive and HHT-resistant mice groups. (c) The expression level and OS of GNAI1 and CALCRL in online database GEPIA in AML. (d) The OS and EFS of CALCRL in our cohorts. (e) The correlation of GNAI1 and CALCRL through correlation analysis. and .

4. Discussion

With the widespread use of HHT in AML treatment, drug resistance has become a great challenge, and few attempts have been made to explore the resistance mechanisms. In this study, we successfully constructed six HHT-resistant cell lines in two HHT-sensitive AML cell lines with different RIs. The growth rate of HHT-resistant cell lines was decreased and accompanied by a delay in the cell cycle and an increased p-gp (CD243) expression. GPRs mediate immune activity, for molecular mechanisms may play an essential role in HHT resistance. The GPR signaling, especially two GPRs, CALCRL and GNAI1, may mediate HHT resistance and have prognostic significance in AML.

Drug-resistant cell lines can be used to study biological resistance changes. In this study, we constructed different RIs of cell lines to observe the resistance characteristics and mechanism dynamically. Recent studies have shown that the cell cycle plays a critical role in chemotherapeutic resistance in cancer. When cancer cells exhibit Palbociclib resistance, cells are delayed in the G1 phase [21]. Further studies revealed that cancer cells with delayed G1 phase could affect intracellular platinum concentrations and protect them from chemotherapy-induced damage [22]. Another study found a lower percentage of S phase in paclitaxel-resistant cell lines of ovarian carcinoma [23]. These studies were consistent with our results in HHT-resistant cell lines, in which G0/G1 phase was delayed and S phase was decreased in HHT-resistant cells.

The increase of p-gp expression in resistant cells was another crucial difference between resistant and sensitive cells. Overexpression of p-gp is considered the primary cause of multidrug resistance (MDR) in AML (26, 27). p-gp functions as an ATP-dependent membrane protein and is involved in the absorption, distribution, and excretion of drugs, xenobiotics, and endogenous compounds [24]. These results indicate that cell cycle delay reduces cell growth rate, and increasing p-gp expression may induce HHT resistance in AML cells.

After KEGG and hub gene analyses, we found that immune activity may play an important role in HHT resistance, and top 10 upregulated hub genes were focused on GPRs. GPRs can regulate many fundamental biological processes, such as cardiovascular functions, immune responses, and hormones [25]. Dysregulated expression and aberrant function of GPRs have been linked to many human diseases including cancer [26]. Cancer immunotherapy targeting GPRs is a promising treatment strategy [27]. In this study, the critical DEGs after HHT resistance belonged to GPRs, indicating that immune regulation by GRPs may play an essential role in HHT resistance. We focused on upregulated DEGs which belonged to GPRs and found that two upregulated DEGs, CALCRL (31, 32) and GNAI1 (33, 34), were highly expressed in primary HHT-resistant cell lines and newly constructed HHT-resistant cell lines. These indicated that CALCRL and GNAI1 may mediate HHT resistance in AML cells. GNAI1 has been confirmed as an oncogene and therapeutic target in serous ovarian cancer [28] and hepatocellular carcinoma [29]. However, its role in AML remains poorly understood. CALCRL is a potential therapeutic target for AML [30]. Another study found that CALCRL increased chemotherapy resistance (Ara-c and daunorubicin) and stem cell properties in AML (32). These confirmed the important role of the two genes in cancer therapy including AML. From UniProt database, we knew that CALCRL was localized in the cell membrane, but GNAI1 was widely localized in various parts of the cell. As we found that they had significant positive correlation, whether they interact to mediate HHT resistance in AML is unclear. So, in this study, we did not further explore the exact role and mechanism of CALCRL and GNAI1 in regulating HHT resistance in AML; we also did not reveal the interaction between CALCRL and GNAI1; this requires future detailed studies. Further research using targeting knocked down or overexpression methods is needed.

5. Conclusion

In summary, we successfully constructed six HHT-resistant AML cell lines with different RIs; these will be specific cell tools to study HHT resistance. The cell cycle delay and increase in p-gp expression may be the biological reasons for HHT resistance. At the transcriptional level, we found that GPRs, especially CALCRL and GNAI1, may mediate HHT resistance. There are some limitations in this study; although we found that GPRs may mediate HHT, we did not study how to mediate HHT resistance, and we need to further explore how CALCRL and GNAI1 regulate HHT resistance in AML and even chemotherapy resistance in AML. The interrelationship and role of CALCRL and GNAI1 also require further research in the future. We only explored the possible mechanism of HHT resistance in AML at the transcriptional level. Further, we need to use the constructed HHT-resistant cell lines to comprehensively study the HHT resistance mechanism at the metabolic and proteomic levels and to explore whether each omics level there is a cross.

Abbreviations

HHT:	Homoharringtonine
HAA:	HHT + aclarubicin + cytarabine
AML:	Acute myeloid leukemia
RI:	Resistance index
p-gp:	P-glycoprotein (CD243)
KEGG:	Kyoto Encyclopedia of Genes and Genomes
GO:	Gene Ontology
GPR:	G-protein coupled receptor
CALCRL:	Calcitonin Receptor-Like
GNAI1:	G Protein Subunit Alpha I1
CML:	Chronic myeloid leukemia
MDS:	Myelodysplastic syndrome
CR:	Complete remission
EFS:	Event-free survival
DA:	Daunorubicin plus cytarabine
IA:	Idarubicin plus cytarabine
IC₅₀:	Half-maximal inhibitory concentration
APL:	Acute promyelocytic leukemia
PI:	Propidium
RT-qPCR:	Real-time quantitative PCR
OS:	Overall survival
DEGs:	Differentially expressed genes
MDR:	Multidrug resistance.

Data Availability

The datasets used and analyzed during this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Fenglin Li, Ying Lu, and Jie Jin contributed to research design, data methodology, data interpretation, and final manuscript writing and approval. Qing Ling, Chao Hu, and Huafeng Wang contributed to literature search, data integration, molecular genetic studies, and approval of the final manuscript. Wenle Ye, Xia Li, Xiangjie Lin, Wenwen Wei, Xin Huang, Yu Qian, and Haihui Zhuang contributed to data verification, as well as final manuscript discussion and approval. Fenglin Li and Qing Ling contributed equally to this work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant 81820108004), Basic Public Welfare Research Project of Zhejiang Province (LGF22H080007), Medical Health Science and Technology Plan Project of Zhejiang (2020KY912), and Natural Science Foundation of Ningbo (202003N4292).

Supplementary Materials

Figure S1: the phenotype of HHT-resistant and HHT-sensitive strains. A: the apoptosis of HHT-resistant cell lines and HHT-sensitive cell lines without HHT cultivate. B: the HE stains of mice tumor mass between HHT-sensitive and HHT-resistant groups. C: the fold change of tumor volume after HHT treatment. NS: not significant; . Figure S2: the verification of DEGs in three MV4-11 resistant cell lines compared with MV4-11 sensitive cell lines. A: fold change of upregulated gene expression in three MV4-11 resistant cell lines via RT-qPCR. B: fold change of downregulated gene expression in three MV4-11 resistant cell lines via RT-qPCR. Figure S3: intersection between MV4-11 R10 versus MV4-11, MV4-11 R30 versus MV4-11, and MV4-11 R50 versus MV4-11. A: in these three sets, there were 227 upregulated genes. B: in these three sets, there were 243 downregulated genes. Figure S4: BP about downregulated genes among MV4-11 resistant strains and sensitive ones. A: BP analysis among MV4-11 R10 and MV4-11. B: BP analysis among MV4-11 R30 and MV4-11. C: BP analysis among MV4-11 R50 and MV4-11. Figure S5: BP analysis of upregulated hub genes in MV4-11 resistant strains. A: BP analysis of upregulated hub genes between MV4-11 R10 resistant strains and MV4-11 sensitive ones. B: BP analysis of upregulated hub genes between MV4-11 R30 resistant strains and MV4-11 sensitive ones. C: BP analysis of upregulated hub genes between MV4-11 R50 resistant strains and MV4-11 sensitive ones. Figure S6: BP analysis of downregulated hub genes in MV4-11 resistant strains. A: BP analysis of downregulated hub genes between MV4-11 R10 resistant strains and MV4-11 sensitive ones. B: BP analysis of downregulated hub genes between MV4-11 R30 resistant strains and MV4-11 sensitive ones. C: BP analysis of downregulated hub genes between MV4-11 R50 resistant strains and MV4-11 sensitive ones. (Supplementary Materials)

References

E. H. Estey, “Acute myeloid leukemia: 2019 update on risk-stratification and management,” American Journal of Hematology, vol. 93, no. 10, pp. 1267–1291, 2018.
View at: Publisher Site | Google Scholar
X. Li, Y. Dong, Y. Li et al., “Low-dose decitabine priming with intermediate-dose cytarabine followed by umbilical cord blood infusion as consolidation therapy for elderly patients with acute myeloid leukemia: a phase II single-arm study,” BMC Cancer, vol. 19, no. 1, p. 819, 2019.
View at: Publisher Site | Google Scholar
U. Krug, C. Röllig, A. Koschmieder et al., “Complete remission and early death after intensive chemotherapy in patients aged 60 years or older with acute myeloid leukaemia: a web-based application for prediction of outcomes,” The Lancet, vol. 376, no. 9757, pp. 2000–2008, 2010.
View at: Publisher Site | Google Scholar
B. Löwenberg, G. J. Ossenkoppele, W. van Putten et al., “High-dose daunorubicin in older patients with acute myeloid leukemia,” New England Journal of Medicine, vol. 361, no. 13, pp. 1235–1248, 2009.
View at: Publisher Site | Google Scholar
H. Hackl, K. Astanina, and R. Wieser, “Molecular and genetic alterations associated with therapy resistance and relapse of acute myeloid leukemia,” Journal of Hematology & Oncology, vol. 10, no. 1, p. 51, 2017.
View at: Publisher Site | Google Scholar
H. M. Kantarjian, M. Talpaz, V. Santini, A. Murgo, B. Cheson, and S. M. O’Brien, “Homoharringtonine: history, current research, and future direction,” Cancer, vol. 92, no. 6, pp. 1591–1605, 2001.
View at: Google Scholar
J. Jin, J. X. Wang, F. F. Chen et al., “Homoharringtonine-based induction regimens for patients with de-novo acute myeloid leukaemia: a multicentre, open-label, randomised, controlled phase 3 trial,” The Lancet Oncology, vol. 14, no. 7, pp. 599–608, 2013.
View at: Publisher Site | Google Scholar
L. Leukemia and Group CSoHCMA, “Chinese guidelines for diagnosis and treatment of adult acute myeloid leukemia (not APL) (2017),” Zhonghua Xue Ye Xue Za Zhi, vol. 38, no. 3, pp. 177–182, 2017.
View at: Google Scholar
Y. Zhang, X. Li, X. Weng et al., “Optimization of idarubicin and cytarabine induction regimen with homoharringtonine for newly diagnosed acute myeloid leukemia patients based on the peripheral blast clearance rate: a single-arm, phase 2 trial (RJ-AML 2014),” American Journal of Hematology, vol. 97, no. 1, pp. 43–51, 2022.
View at: Publisher Site | Google Scholar
M. McDermott, A. J. Eustace, S. Busschots et al., “In vitro development of chemotherapy and targeted therapy drug-resistant cancer cell lines: a practical guide with case studies,” Frontiers Oncology, vol. 4, no. 40, 2014.
View at: Publisher Site | Google Scholar
T. Takimoto-Shimomura, H. Nagoshi, S. Maegawa et al., “Establishment and characteristics of a novel mantle cell lymphoma-derived cell line and a bendamustine-resistant subline,” Cancer Genomics & Proteomics, vol. 15, no. 3, pp. 213–223, 2018.
View at: Publisher Site | Google Scholar
J. Park, E. K. Bae, C. Lee et al., “Establishment and characterization of bortezomib-resistant U266 cell line: constitutive activation of NF-κB-mediated cell signals and/or alterations of ubiquitylation-related genes reduce bortezomib-induced apoptosis,” BMB Reports, vol. 47, no. 5, pp. 274–279, 2014.
View at: Publisher Site | Google Scholar
O. Lindblad, E. Cordero, A. Puissant et al., “Aberrant activation of the PI3K/mTOR pathway promotes resistance to sorafenib in AML,” Oncogene, vol. 35, no. 39, pp. 5119–5131, 2016.
View at: Publisher Site | Google Scholar
Y. Saijo, N. Kumano, S. Suzuki et al., “Establishment of an etoposide (VP-16)-resistant subline of THP-1 human monocytic leukemia cell line,” Tohoku Journal of Experimental Medicine, vol. 157, no. 3, pp. 215–219, 1989.
View at: Publisher Site | Google Scholar
T. Ohnoshi, T. Ohnuma, I. Takahashi, K. Scanlon, B. A. Kamen, and J. F. Holland, “Establishment of methotrexate-resistant human acute lymphoblastic leukemia cells in culture and effects of folate antagonists,” Cancer Research, vol. 42, no. 5, pp. 1655–1660, 1982.
View at: Google Scholar
C. H. Chin, S. H. Chen, H. H. Wu, C. W. Ho, M. T. Ko, and C. Y. Lin, “cytoHubba: identifying hub objects and sub-networks from complex interactome,” BMC Systems Biology, vol. 8, no. S4, p. S11, 2014.
View at: Publisher Site | Google Scholar
M. Yu, J. Wang, Z. Zhu et al., “Prognostic impact of MYH9 expression on patients with acute myeloid leukemia,” Oncotarget, vol. 8, no. 1, pp. 156–163, 2017.
View at: Publisher Site | Google Scholar
W. L. Chen, J. H. Wang, A. H. Zhao et al., “A distinct glucose metabolism signature of acute myeloid leukemia with prognostic value,” Blood, vol. 124, no. 10, pp. 1645–1654, 2014.
View at: Publisher Site | Google Scholar
D. Szklarczyk, J. H. Morris, H. Cook et al., “The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible,” Nucleic Acids Research, vol. 45, no. D1, pp. D362–D368, 2017.
View at: Publisher Site | Google Scholar
N. T. Doncheva, J. H. Morris, J. Gorodkin, and L. J. Jensen, “Cytoscape StringApp: network analysis and visualization of proteomics data,” Journal of Proteome Research, vol. 18, no. 2, pp. 623–632, 2019.
View at: Publisher Site | Google Scholar
D. W. Fry, P. J. Harvey, P. R. Keller et al., “Specific inhibition of cyclin-dependent kinase 4/6 by PD 0332991 and associated antitumor activity in human tumor xenografts,” Molecular Cancer Therapeutics, vol. 3, no. 11, pp. 1427–1438, 2004.
View at: Publisher Site | Google Scholar
J. M. Replogle, W. Zhou, A. E. Amaro et al., “Aneuploidy increases resistance to chemotherapeutics by antagonizing cell division,” Proceedings of the National Academy of Sciences of the USA, vol. 117, no. 48, pp. 30566–30576, 2020.
View at: Publisher Site | Google Scholar
J. Zhang, J. Zhao, W. Zhang et al., “Establishment of paclitaxel-resistant cell line and the underlying mechanism on drug resistance,” International Journal of Gynecological Cancer, vol. 22, no. 9, pp. 1450–1456, 2012.
View at: Publisher Site | Google Scholar
M. Huls, F. G. M. Russel, and R. Masereeuw, “The role of ATP binding cassette transporters in tissue defense and organ regeneration,” Journal of Pharmacology and Experimental Therapeutics, vol. 328, no. 1, pp. 3–9, 2009.
View at: Publisher Site | Google Scholar
R. Fredriksson, M. C. Lagerstrom, L. G. Lundin, and H. B. Schioth, “The G-protein-coupled receptors in the human genome form five main families. Phylogenetic analysis, paralogon groups, and fingerprints,” Molecular Pharmacology, vol. 63, no. 6, pp. 1256–1272, 2003.
View at: Publisher Site | Google Scholar
R. Bar-Shavit, M. Maoz, A. Kancharla et al., “G protein-coupled receptors in cancer,” International Journal of Molecular Sciences, vol. 17, no. 8, Article ID E1320, 2016.
View at: Publisher Site | Google Scholar
V. Wu, H. Yeerna, N. Nohata et al., “Illuminating the Onco-GPCRome: novel G protein-coupled receptor-driven oncocrine networks and targets for cancer immunotherapy,” Journal of Biological Chemistry, vol. 294, no. 29, pp. 11062–11086, 2019.
View at: Publisher Site | Google Scholar
S. J. Zhan, B. Liu, and H. Linghu, “Identifying genes as potential prognostic indicators in patients with serous ovarian cancer resistant to carboplatin using integrated bioinformatics analysis,” Oncology Reports, vol. 39, no. 6, pp. 2653–2663, 2018.
View at: Publisher Site | Google Scholar
J. Yao, L. H. Liang, Y. Zhang et al., “GNAI1 suppresses tumor cell migration and invasion and is post-transcriptionally regulated by mir-320a/c/d in hepatocellular carcinoma,” Cancer Biol Med, vol. 9, no. 4, pp. 234–241, 2012.
View at: Publisher Site | Google Scholar
L. Angenendt, E. Bormann, C. Pabst et al., “The neuropeptide receptor calcitonin receptor-like (CALCRL) is a potential therapeutic target in acute myeloid leukemia,” Leukemia, vol. 33, no. 12, pp. 2830–2841, 2019.
View at: Publisher Site | Google Scholar
J. Prada-Arismendy, J. C. Arroyave, and S. Rothlisberger, “Molecular biomarkers in acute myeloid leukemia,” Blood Reviews, vol. 31, no. 1, pp. 63–76, 2017.
View at: Publisher Site | Google Scholar
T. Tsuruo, H. Iida, E. Ohkochi, S. Tsukagoshi, and Y. Sakurai, “Establishment and properties of vincristine-resistant human myelogenous leukemia K562,” Gann, vol. 74, no. 5, pp. 751–758, 1983.
View at: Google Scholar
F. Gao, W. Dong, W. Yang, J. Liu, Z. Zheng, and K. Sun, “Expression of P-gp in acute myeloid leukemia and the reversal function of As2O3 on drug resistance,” Oncology Letters, vol. 9, no. 1, pp. 177–182, 2015.
View at: Publisher Site | Google Scholar
T. Gluexam, A. M. Grandits, A. Schlerka et al., “CGRP signaling via CALCRL increases chemotherapy resistance and stem cell properties in acute myeloid leukemia,” International Journal of Molecular Sciences, vol. 20, no. 23, p. 5826, 2019.
View at: Publisher Site | Google Scholar
Z. W. Li, B. Sun, T. Gong et al., “GNAI1 and GNAI3 reduce colitis-associated tumorigenesis in mice by blocking IL6 signaling and down-regulating expression of GNAI2,” Gastroenterology, vol. 156, no. 8, pp. 2297–2312, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Fenglin Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

252

Downloads

536

Citations