Aurora kinase A (AURKA) promotes the progression and imatinib resistance of advanced gastrointestinal stromal tumors
Cancer Cell International volume 21, Article number: 407 (2021)
Gastrointestinal stromal tumor (GIST) is a common tumor that originates from the alimentary system mesenchyme. Compared to typical gastrointestinal carcinomas, GISTs exhibit unique malignant behaviors. Bioinformatic tools and subsequent experiments were applied to investigate novel targets involved in GIST progression and imatinib resistance.
Differences in gene expression profiles between advanced and nonadvanced GISTs were comprehensively analyzed based on the Gene Expression Omnibus (GEO) dataset GSE136755. A protein–protein interaction (PPI) network was constructed to identify the potential target gene. Gene set enrichment analysis (GSEA) was used to elucidate relevant biological events related to the target gene based on the GSE47911 dataset. Subsequently, immunohistochemistry and Kaplan–Meier analysis were performed to validate the prognostic value of the target gene in GISTs. Overexpression of the target gene was conducted to analyze its function in the proliferation, apoptosis, and imatinib resistance of GIST/T1 cells.
In the current study, a total of 606 differentially expressed genes (DEGs) were screened based on the GSE136755 dataset, and the upregulated DEGs in advanced GISTs were mainly involved in cell division through functional annotations. The intersecting hub gene, Aurora kinase A (AURKA), was identified by degree and bottleneck algorithms. GSEA revealed that AURKA was involved in cell cycle-related biological processes. Analysis of the Oncomine and GEPIA databases revealed a pattern of elevated AURKA expression in most human malignances. Clinical assays demonstrated that AURKA could be an independent prognostic factor for GISTs. Additionally, overexpression of AURKA was experimentally demonstrated to promote cell proliferation, inhibit cell apoptosis, and enhance imatinib resistance in GIST/T1 cells.
These findings indicated that overexpression of AURKA promoted GIST progression and enhanced imatinib resistance, implying that AURKA is a potential therapeutic target for GISTs.
Gastrointestinal stromal tumor (GIST) is the most common mesenchymal tumor of the alimentary system and originates from the interstitial cells of Cajal (ICC) . Approximately 71% of GISTs present with KIT (71%) or PDGFRα (14%) mutations, and 10–15% of GISTs do not have KIT or PDGFRα mutations, referred to as KIT/PDGFRα wild type GISTs [2, 3]. The malignant potential of GISTs is stratified based on tumor size, mitotic index, and location according to the modified NIH criteria . Based on these criteria, GISTs are classified as high-risk, intermediate-risk, low-risk, and very low-risk. The 5-year survival rate of patients with advanced GIST is between 35 and 65% . The main therapeutic option for primary localized GISTs is surgical resection. However, the recurrence rate for GISTs, even after complete surgical resection, is 40–80% . The median time to recurrence for most patients is approximately 12–16 months . Treatment with imatinib, a tyrosine kinase inhibitor (TKI) that targets KIT and PDGFRα, has improved the prognosis of GIST patients. However, when used to eliminate mature GIST cells, imatinib has limited efficacy, and studies have revealed that GIST persists with prolonged TKI therapy . Due to acquired resistance to imatinib, approximately 85–90% of patients with GIST experience disease progression within 20–24 months [2, 3, 7]. Further research is urgently needed to reveal the mechanism of GIST progression and explore novel therapeutic targets to combat imatinib resistance.
Currently, bioinformatic tools are being used to evaluate the molecular signatures associated with progression and clinical outcomes in several types of malignancies [8,9,10]. In this study, the GSE136755 and GSE47911 datasets were downloaded from the Gene Expression Omnibus (GEO) and used to evaluate the potential target genes involved in GIST progression. Among several hub genes in the GSE136755 dataset, AURKA was considered a key hub gene in GIST progression. Gene set enrichment analysis (GSEA) based on the GSE47911 dataset indicated that AURKA promotes GIST progression by regulating cell cycle processes. Subsequent clinical data analyses demonstrated the value of AURKA as a prognostic factor for GISTs. Furthermore, overexpression of AURKA was experimentally demonstrated to significantly promote GIST/T1 cell proliferation, inhibit apoptosis, and enhance their resistance to imatinib.
Materials and methods
Data acquisition and differentially expressed gene (DEG) identification
The GSE136755 dataset was downloaded from the Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) . GSE136755 is based on the GPL17077 platform (Agilent-039494 SurePrint G3 Human GE v2 8 × 60 K Microarray 039381) and includes clinicopathological information for 65 human GIST tumor samples without preoperative imatinib treatment. Primary GISTs with KIT mutations were selected to screen the differentially expressed genes (DEGs). Advanced GIST samples were defined as samples from patients with high-risk GISTs (16 samples), while nonadvanced GIST samples were defined as samples from patients with low-risk and very low-risk GISTs (31 samples). Intermediate-risk GISTs were not included because they do not have clear distinct biological behaviors in comparison with high-risk and low-risk GISTs. GEO2R, an R-associated web tool from the National Center for Biotechnology Information, was used to screen the DEGs between advanced and nonadvanced GISTs . The DEGs were identified using the cutoff values of |log2FoldChange| > 1 and adju. p < 0.05. For hierarchical clustering analysis in Morpheus (https://software.broadinstitute.org/morpheus), the DEGs were downloaded in text format.
Functional enrichment analysis of DEGs
The Database for Annotation, Visualization, and Integrated Discovery (DAVID) (http://david.ncifcrf.gov, version 6.8) is a web-based bioinformatics resource that is used to extract genes functional annotation information . Gene ontology (GO) is a major bioinformatics tool for gene annotation and analysis [14, 15]. The Kyoto Encyclopedia of Genes and Genomes (KEGG) is a popular database for the analysis of advanced gene functions and potential signaling pathways in large-scale molecular data [16, 17]. DAVID was used to perform GO and KEGG enrichment analyses of DEGs. The cutoff criterion was a false discovery rate (FDR) < 0.05.
Protein–protein interaction (PPI) network and module analysis
The Search Tool for the Retrieval Interacting Genes (STRING) (http://string-db.org) is an online database that is used to identify interactions among DEGs [18, 19]. A confidence score ≥ 0.7 was set for conducting the PPI network. Cytoscape (version 3.7.1) is an open-source bioinformatics software platform that is used for visualizing the PPI network and for further analyses [20, 21]. The Molecular Complex Detection (MCODE) plugin in Cytoscape was used to identify significant modules based on the PPI network topology. The criteria were degree cutoff = 2, node score cutoff = 0.2, K-core = 2, and max. depth = 100. The plug-in app ClueGO was used to analyze and visualize the biological processes and pathways in significant modules.
Hub gene identification and analysis
The cytoHubba plugin in Cytoscape was used to identify hub genes based on the degree and bottleneck algorithms. The Oncomine database is an online platform that computes gene expression signatures, clusters, and gene-set modules . The Gene Expression Profiling Interactive Analysis (GEPIA) database is a newly developed web server for cancer, normal gene expression profiling and interactive analysis . CytoHubba was used to identify the hub genes whose expression patterns have been evaluated in common human malignancies using the Oncomine and GEPIA databases.
Clinicopathological features and KIT/PDGFRα mutation types were extracted from the GSE136755 dataset and from the raw data provided by Lagarde et al. . Correlations between key gene expression patterns and clinicopathological features as well as KIT/PDGFRα mutation types were statistically determined, with P < 0.05 set as the threshold for statistical significance.
Gene set enrichment analysis (GSEA)
GSEA is a computational method that determines whether an a priori defined gene set shows statistically significant and concordant differences between two biological states [25, 26]. GSEA computes biological information from different perspectives and further elucidates relevant biological events. The GSE47911 dataset, based on the GPL6480 platform (Agilent-014850 Whole Human Genome Microarray 4 × 44 K G4112F), was also downloaded from the GEO database. This dataset has 15 gastric GIST samples . Based on the key gene expression level, GIST samples were divided into two groups, after which GSEA was subsequently performed. Annotated gene sets (c2. cp.kegg.v7.2.symbols.gmt [Curated], c2.cp.ractome.v7.2.symbols.gmt [Curated], and c5.bp.v7.2.symbols.gmt. [Gene Oncology]) were chosen as the reference gene sets. Gene size > 20, FDR < 0.05 and normalized enrichment score (NES) > 2.00 were set as the cutoff criteria.
Immunohistochemistry and survival analysis
Between 2001 and 2015, a total of 49 patients admitted to the First Affiliated Hospital of Zhejiang University (Zhejiang Province, China) who were diagnosed with GISTs were enrolled in this study. The GIST patients who had incomplete resection, neoadjuvant or adjuvant imatinib treatment, or a family history of GIST, were excluded in this study. Clinical stratification of GISTs was based on the modified NIH criteria . Paraffin-embedded GIST samples were obtained from the study participants and analyzed by immunohistochemical (IHC) staining and were used for the survival analysis. All GIST tissue samples were provided by the Department of Pathology within the First Affiliated Hospital of Zhejiang University. IHC staining was performed as previously described . Briefly, tissue sections were incubated at 4 °C overnight with anti-human AURKA rabbit polyclonal antibody diluted 1:500 (NOVUS Biologicals, USA). A total of five adjacent fields using 400× magnification in areas with the highest density of positive staining were scored according to the summation of the percentage of staining intensity. The immunostaining percentage was defined as 0 (< 5%), 1 (< 20%), 2 (20–50%), and 3 (> 50%). Staining intensity was defined from 0 (no staining) to 3 (strongest staining). The maximum score of IHC staining was 6, in which > 50% of the cells had the strongest staining intensity. Staining scores lower than the mean value were considered low expression, while scores higher than the mean value were considered high expression. The follow-up time for all the patients was calculated from the date of surgery to the date of disease recurrence or last visit. The use of human tumor samples and clinical data in this study was approved by the Ethical Committee of the First Affiliated Hospital of Zhejiang University. Study participants were required to sign written informed consent before enrollment.
The GIST cell line (GIST/T1) donated by Prof. Wenbin Chen (Zhejiang University, Hangzhou, China) was obtained from the Cell Bank of the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China). Cells were maintained in RPMI-1640 medium (GIBCO) supplemented with 10% heat-inactivated fetal bovine serum (GIBCO) in a humidified atmosphere with 5% CO2 at 37 °C.
Construction and transfection of lentiviral vectors for AURKA overexpression
The DNA fragment encoding the AURKA sequence was synthesized and inserted into the lentivirus expression vector pLVX-IRES-tdtomato (TaKaRa, China). The resulting vector was identified as pLVX-AURKA-IRES-tdtomato. Lentiviral plasmids were transfected into HEK 293T cells with psPAX2 and pMD2. G plasmids at a ratio of 4:3:1 using Lipofectamine 2000 (Invitrogen, USA). Then, 48 h after transfection, the virus was isolated. GIST/T1 cells were infected with lentivirus for 48 h, and the transfection efficiency was measured by PCR and western blotting.
RNA extraction and real-time quantitative PCR (RT-qPCR)
Total RNA was extracted using TRIzol reagent (Generay Biotech, China) according to the manufacturer’s instructions. The extracted RNA was then treated with RQ1 RNase-Free DNase (Promega, USA). Then, reverse transcription was performed using PrimeScriptTM RT Master Mix (Takara, China) according to the manufacturer’s instructions. RT-qPCR analysis was performed to measure the expression levels of AURKA using the CFX Connect Real-Time System (BIO-RAD, USA) with the SuperReal PreMix Color SYBR Green kit (Tiangen, China). The primer sets used for RT-qPCR are shown in Table 1. Gene expression levels were normalized against the internal control using the 2−ΔΔCT method.
Protein lysates from each sample were separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to a polyvinylidene difluoride (PVDF) membrane (Merck Millipore, USA). The membrane was blocked with 5% nonfat milk or 5% BSA in TBST (tris-buffered saline with 0.1% Tween 20) for 1 h at room temperature. Then, the membrane was incubated with a primary antibody against AURKA (1:1000 dilution, NOVUS Biologicals, USA) overnight at 4 °C, washed three times using 0.1% TBST buffer for 30 min, probed with goat anti-rabbit IgG-HRP secondary antibody (1:2000 dilution; Santa Cruz, USA), and washed three times using 0.1% TBST buffer for 30 min. The HRP-conjugated secondary antibody was detected and visualized using an enhanced chemiluminescence detection system (GE Healthcare, USA). Band intensity was quantified by densitometry using ImageJ software (version 1.49; National Institutes of Health, USA) .
Cell proliferation was assessed by the CCK-8 assay (Tongren Chemical Society, Japan). Briefly, 1 × 105 cells per well were seeded in 96-well plates and incubated in a 5% CO2 atmosphere at 37 °C for 24 h. Imatinib (3 µM) was then added to the culture medium for 48 h to evaluate the effect of AURKA on imatinib resistance. The medium was discarded, replaced with serum-free medium and CCK8 (10 µl), and incubated for 2 h. The Biokinetics Reader (MD corporate, USA) was used to detect absorbance at 450 nm.
Cell apoptosis was assessed using the Annexin V-APC/7-AAD apoptosis kit (MultiSciences, China) according to the manufacturer’s instructions. Briefly, 3 × 105 cells per well were seeded in 6-well plates and incubated in a 5% CO2 atmosphere at 37 °C for 24 h. Imatinib (3 µM) or solvent was added to the culture medium for 48 h, after which the cells were collected. After being washed twice in PBS at 4 °C, cells were resuspended in binding buffer (500 µl). Annexin V-APC (5 µl) and 7-AAD (10 µl) were added to the suspension and the mixture was incubated for 5 min at 4 °C. The apoptosis index was examined by flow cytometry (ACCURI C6; BD, USA).
All statistical analyses were performed using SPSS 25.0 software (SPSS Inc., USA). Descriptive data were expressed as the mean ± SD. Comparison of more than two mean values was performed by one-way analysis of variance (ANOVA), while Student’s t-test was used to compare two mean values. Kaplan–Meier analysis was performed to establish disease-free survival (DFS) curves, while the log-rank test was used for survival curve comparison. The Cox proportional hazards model was then used to perform multivariable analysis, while the forward likelihood ratio method was used to identify independent variables. P < 0.05 was set as the threshold for statistical significance.
Identification of differentially expressed genes (DEGs)
A total of 65 human GIST samples without preoperative imatinib treatment were included in the GSE13675 dataset. A cohort of 47 samples with KIT mutations, consisting of 16 patients with advanced (high-risk) GIST and 31 with nonadvanced (low-risk and very low-risk) GIST, was used to screen the DEGs using GEO2R. A total of 606 genes (244 with upregulated expression and 362 with downregulated expression) were identified using the cutoff criteria of adj. p < 0.05 and log2FoldChange > 1. The top 50 genes with up-and downregulated expression are presented as a heatmap (Fig. 1A, Additional file 1).
Functional annotation of DEGs
To evaluate the biological clustering of DEGs, GO and KEGG analyses for the up- and downregulated DEGs were performed using DAVID. Based on GO analysis, upregulated DEGs were found to be significantly enriched in cell division, sister chromatid cohesion, mitotic nuclear division, microtubule-based movement, and mitotic metaphase plate congression. Cellular components (CCs) of the upregulated DEGs were significantly enriched in the midbody, kinesin complex, spindle, spindle microtubule, and condensed chromosome kinetochore. The molecular function (MF) of the upregulated DEGs was significantly enriched in microtubule motor activity and microtubule binding. KEGG analysis showed that the upregulated DEGs were mainly involved in the cell cycle and oocyte meiosis. GO analysis revealed that the downregulated DEGs were significantly enriched in the interferon-gamma-mediated signaling pathway, type I interferon signaling pathway, and antigen processing and presentation of antigens. The CC of the downregulated DEGs was significantly enriched in the integral component of the luminal side of the endoplasmic reticulum membrane and ER to Golgi transport vesicle membrane, while the MF of the downregulated DEGs was significantly enriched in peptide antigen binding and MHC class II receptor activity. KEGG analysis further revealed that the downregulated DEGs were mainly involved in graft-versus-host disease, type I diabetes mellitus, allograft rejection, antigen processing and presentation, and viral myocarditis. These results are presented in Fig. 1B (Additional file 2).
Module analysis through the PPI network of DEGs
To clarify the DEGs functionally, STRING was used to construct a PPI network, which was composed of 603 nodes and 1768 edges. The PPI enrichment p-value was < 1.0E−16. The PPI network was visualized by Cytoscape and further analyzed by the MCODE plugin. The two most significant modules (module 1 and module 2) were identified and analyzed using GO, KEGG and REACTOME annotations to infer their biological functions. Module 1 (MCODE score = 36.667) was mainly involved in cell cycle-related biological processes and signaling pathways, while module 2 (MCODE score = 18.759) was mainly involved in immunological processes and signaling pathways. These results are presented in Fig. 2A (Additional file 2).
Identification of hub genes
To reveal the crucial genes underlying the regulation of GIST progression, we filtered hub genes among DEGs using the cytoHubba plugin of Cytoscape. Two algorithms, degree and bottleneck, were applied to weight the DEGs. The degree algorithm calculates the relevance and abundance of genes, while the bottleneck algorithm evaluates key gene positions in an entire regulatory network. According to the degree algorithm, the top 15 hub genes were CDK1, KIF11, KIF2C, CENPE, KIF20A, BUB1, CCNA2, CCNB1, AURKA, MAD2L1, CDCA8, KIF4A, CENPF, NDC80, and KIF23, with scores ranging from 58 to 48. According to the bottleneck algorithm, the top 15 hub genes were AURKA, FN1, CD44, VEGFA, IL6, HLA-DQB1, HLA-DPA1, CXCL8, NT5E, ANK2, FOXM1, CHEK1, STAT1, CDC25A, and IFIH1, with scores ranging from 64 to 8. A Venn diagram was used to identify the intersection of the key hub genes between the two hub gene cohorts. The results showed that AURKA was the only overlapping hub gene (Fig. 2B).
Gene set enrichment analysis (GSEA)
There was a total of 15 human gastric GIST samples in the GSE47911 dataset, comprised of 6 high-risk patients, 1 intermediate-risk patient, 3 low-risk patients, and 5 very low-risk patients. To further verify significant biological processes associated with AURKA expression, GSE47911 gene profiles were divided into two groups, after which GSEA was performed based on the AURKA expression level. Samples with the highest (25%, 4 samples) and lowest (25%, 4 samples) expression levels were selected for further analysis using GSEA. Cell cycle-related gene sets were associated with elevated AURKA expression (Fig. 2C, Additional file 3).
Correlation between AURKA expression and the clinicopathological features of GISTs
To evaluate the clinical significance of AURKA expression in GISTs, AURKA expression levels in 49 GIST tissues were assessed by IHC staining (Fig. 3A). The correlations between AURKA expression and clinicopathological features (age, sex, location and risk stratification) were determined (Table 2, Additional file 4). AURKA expression was closely associated with tumor risk stratification (Fig. 3B; P < 0.001). The clinical significance of AURKA expression in GISTs was also evaluated using the data from GSE136755 and the raw data provided by Lagarde et al.  (Tables 3 and 4, and Fig. 4, Additional file 5). Findings from the GSE136755 dataset analysis revealed significant associations between AURKA expression and tumor risk stratification (P < 0.001) as well as tumor stage (P < 0.001) and the analyses of the raw data from Lagarde et al.  also showed a significant association between AURKA expression and tumor risk stratification (P < 0.001) as well as tumor recurrence (P < 0.001) and metastasis (P < 0.001). However, apart from GSE136755, which revealed a significant association between AURKA expression and tumor location (P = 0.018), the data provided by Lagarde et al. did not establish a significant association between AURKA expression and tumor location (P = 0.156).
To determine the prognostic value of AURKA expression in GISTs, Kaplan-Meier survival analysis was performed. The range of observation time was 9–79 months. As shown in Fig. 3C, patients with GIST with elevated AURKA expression exhibited poorer DFS than those with low AURKA expression levels (43.25 ± 6.94 months vs. 98.48 ± 3.44 months, P < 0.001). We further included the statistically significant variables (tumor size, mitotic index, risk stratification, and AURKA score) into Cox regression model. The Cox proportional hazards model showed that AURKA could be used as an independent prognostic marker for GISTs (P = 0.002) and the hazard ratio (HR) was 0.087.
Gene mutation types can predict the responsiveness of GISTs to imatinib. GISTs with KIT exon 11, PDGFRα exon 12 and PDGFRα exon 14 mutations were considered sensitive to imatinib. GISTs with other mutations, such as KIT exon 9, KIT exon 13, KIT exon 14, KIT exon 17, KIT exon 18, PDGFRα exon 18 and D842V mutations and KIT/PDGFRα wild-type GISTs, were insensitive/resistant to imatinib [1, 3, 30]. In the GSE136755 dataset, 56 samples were from imatinib-sensitive GISTs and 7 were from imatinib-resistant GISTs; there was a weak association between AURKA expression and imatinib-resistant gene mutations (Fig. 4A, P = 0.074). Through analysis of the raw data provided by Lagarde et al. , comprising 45 samples from imatinib-sensitive GISTs and 15 from imatinib-resistant GISTs, it was shown that there was a strong association between AURKA expression and imatinib-resistant gene mutations (Fig. 4B, P = 0.018).
AURKA expression patterns in common human malignancies
To determine whether elevated AURKA expression is common in human digestive malignancies, mRNA expression levels of AURKA in stomach carcinoma, liver hepatocellular carcinoma, and colorectal carcinoma were evaluated using data from the GEPIA database. AURKA expression was found to be significantly upregulated in all the above malignancies compared to normal tissues. Findings from the Oncomine database also indicated that AURKA expression is upregulated in most human malignancies (Fig. 3D).
AURKA overexpression promotes GIST/T1 cell proliferation and anti-apoptosis
To assess the biological effects of AURKA expression in GISTs, AURKA was overexpressed in GIST/T1 cells by transfection with an AURKA-expressing virus; these cells were defined as the AURKA overexpression group (AURKA group). Normal GIST/T1 cells (blank group) and GIST/T1 cells transfected with vacant plasmids (vector group) were considered the control groups. The transfection efficiency was determined by observing the red fluorescence from the tdtomato reporter and quantified by RT-qPCR and western blotting. Figure 5 shows that compared to the blank and vector groups, AURKA was overexpressed in the AURKA group (Additional file 6).
The CCK-8 assay was performed to assess the effect of AURKA overexpression on cell proliferation. Compared to the blank and vector groups, the overexpression of AURKA in the AUKRA group significantly enhanced GIST/T1 cell proliferation (P = 0.018) (Fig. 6A). Imatinib treatment significantly inhibited cell proliferation in all three groups. However, compared to cells in the control groups, AUKRA group cells still showed a relatively higher proliferation rate in the presence of imatinib (P < 0.001, Fig. 6A; Additional file 6).
We also established that AURKA overexpression markedly suppressed the apoptotic process in GIST/T1 cells (P < 0.001, Fig. 6B). A similar result was observed with imatinib administration. Compared to the blank and vector groups, AURKA overexpression obviously inhibited cell apoptosis after imatinib administration (P < 0.001, Fig. 6B). The results suggested that AURKA overexpression enhanced the resistance of GIST cells to imatinib (Additional file 6).
GIST is a common mesenchymal malignancy of the human alimentary system. Compared to gastrointestinal carcinomas, GISTs are known to possess unique biological features. For example, lymph node metastasis is not common in GIST and preferentially occurs in patients at a young age [31, 32]. A 1- or 2-cm macroscopic margin may be sufficient to achieve microscopically negative margins . GISTs respond poorly to conventional chemotherapy and radiotherapy . As such, to better understand GIST biological behavior and inform the development of therapeutic strategies, it is important to establish the crucial genes that regulate the malignant behavior of GISTs. Bioinformatics advances have been useful for exploring molecular targets indicating the progression and prognosis of GIST [11, 35, 36].
In this study, gene profiles of 47 GIST samples from the GSE136755 dataset were selected for further analyses. Comparison of the gene profiles between advanced and nonadvanced GISTs generated 244 upregulated DEGs and 362 downregulated DEGs. Functional annotation based on GO and KEGG analyses showed that upregulated DEGs were mainly enriched in cell cycle-related biological processes and signaling pathways, while the downregulated DEGs were mainly enriched in immune-related biological processes and signaling pathways. The STRING database and Cytoscape software were used for further exploration of the DEGs. Two important modules were extracted and visualized. Module 1 consisted of upregulated DEGs and was mainly involved in cell cycle-related biological processes and signaling pathways, while module 2 consisted of downregulated DEGs and was mainly involved in immunological processes and signaling pathways. This indicates that the difference in gene profiles between advanced and nonadvanced GISTs is mainly reflected in the cell cycle and tumor immunity.
Based on the degree and bottleneck algorithms, the cytoHubba plugin in Cytoscape software was used to screen for novel key genes associated with GIST progression. The degree algorithm calculates the relevance and abundance of genes, while the bottleneck algorithm evaluates key gene positions in an entire regulatory network. In this study, a significant key gene, AURKA, was identified using a Venn diagram. AURKA is a protein-coding gene that encodes a cell cycle-regulated kinase involved in microtubule formation and/or stabilization at the spindle pole during chromosomal segregation. It has been documented that AURKA promotes tumor progression by enhancing cell cycle progression, cell survival, genomic instability, epithelial-mesenchymal transition (EMT) and stem-like properties of cancer cells . In most solid tumors, AURKA regulates cell cycle checkpoints and promotes the cell cycle process . GSEA based on GSE47911 gene profiles further validated the association between AURKA overexpression and cell cycle progression in GISTs.
To confirm the importance of AURKA expression in GISTs, we performed IHC staining to establish the associations between AURKA expression and the clinicopathological characteristics of the 49 enrolled patients with GIST. In advanced GISTs, the expression level of AURKA was found to be elevated. This result is consistent with the analyses of the data provided by GSE136755 and Lagarde et al. . Survival analysis further showed that AURKA overexpression was a potential independent prognostic factor for patients with GIST. Furthermore, a series of in vitro experiments demonstrated that overexpressing AUKRA in GIST cells promoted cell proliferation and was antiapoptotic, indicating enhanced malignancy. These findings validated the results of the bioinformatics analyses.
Drug resistance is a major obstacle in cancer chemotherapy and greatly affects a patient’s prognosis. Adjuvant imatinib has been widely used as a first-line therapeutic option for patients with advanced GIST [38, 39]. However, the prevalence of imatinib resistance has increased in recent years. By analyzing the raw data provided by Lagarde et al. , a significant relationship was found between AURKA overexpression and gene mutations causing imatinib resistance. The GSE136755 dataset showed a similar result, but the result was not significant, which could be attributable to the small sample size. In addition, in vitro experiments showed that AURKA overexpression enhanced the resistance of GIST cells to imatinib by promoting cell proliferation and inhibiting cell apoptosis.
The role of AURKA overexpression in tumor progression has been reported in a variety of human malignancies. AURKA phosphorylates RPS6KB1 and promotes cell proliferation and anti-apoptosis . AURKA also stabilizes the transcription factor N-MYC, thereby promoting G1/S cell cycle transition and tumor cell proliferation . Pharmacological inhibition of AURKA promotes the chemosensitivity of cervical cancer cells . Compounds targeting AURKA, particularly alisertib, have been extensively studied in preclinical models, where they have shown synergistic effects with other targeted therapies, leading to tumor regression in a variety of cancer models . Yeh et al. confirmed the contribution of the AURKA inhibitor MLN8237 to the suppression of metastatic GISTs . Findings from the GEPIA and Oncomine databases also supported the contribution of AURKA overexpression to tumorigenesis.
This study had the following limitations. First, the case number in each GIST cohort was not large. Therefore, we made up for this deficiency to some extent by incorporating different GIST cohorts for comprehensive analysis. Second, this study is a preliminary exploration and certification of AURKA as a therapeutic target. We did not carry out experiments to investigate the potential of AURKA-targeted therapy. Further validation based on various in vitro and in vivo experiments is required.
In conclusion, our findings demonstrated a significant overexpression of AURKA in advanced GISTs by bioinformatics analyses, which predicts poor patient prognosis. Overexpression of AURKA was experimentally demonstrated to promote the proliferation of GIST cells and inhibit GIST cell apoptosis, which contributes to imatinib resistance, implying that AURKA is a potential therapeutic target for GISTs.
Availability of data and materials
All data and materials generated and analyzed during this study are included in the main paper.
Corless CL, Barnett CM, Heinrich MC. Gastrointestinal stromal tumours: origin and molecular oncology. Nat Rev Cancer. 2011;11:865–78.
Wozniak A, Gebreyohannes YK, Debiec-Rychter M, Schoffski P. New targets and therapies for gastrointestinal stromal tumors. Expert Rev Anticancer Ther. 2017;17:1117–29.
Liu P, Tan F, Liu H, Li B, Lei T, Zhao X. The use of molecular subtypes for precision therapy of recurrent and metastatic gastrointestinal stromal tumor. Onco Targets Ther. 2020;13:2433–47.
Joensuu H. Risk stratification of patients diagnosed with gastrointestinal stromal tumor. Hum Pathol. 2008;39:1411–9.
Iorio N, Sawaya RA, Friedenberg FK. Review article: the biology, diagnosis and management of gastrointestinal stromal tumours. Aliment Pharmacol Ther. 2014;39:1376–86.
Valsangkar N, Sehdev A, Misra S, Zimmers TA, O’Neil BH, Koniaris LG. Current management of gastrointestinal stromal tumors: surgery, current biomarkers, mutations, and therapy. Surgery. 2015;158:1149–64.
Serrano C, George S, Valverde C, Olivares D, Garcia-Valverde A, Suarez C, et al. Novel insights into the treatment of imatinib-resistant gastrointestinal stromal tumors. Target Oncol. 2017;12:277–88.
Zhao X, Xu M, Cai Z, Yuan W, Cui W, Li MD. Identification of LIFR, PIK3R1, and MMP12 as novel prognostic signatures in gallbladder cancer using network-based module analysis. Front Oncol. 2019;9:325.
Yang JF, Shi SN, Xu WH, Qiu YH, Zheng JZ, Yu K, et al. Screening, identification and validation of CCND1 and PECAM1/CD31 for predicting prognosis in renal cell carcinoma patients. Aging. 2019;11:12057–79.
Zheng H, Zhang G, Zhang L, Wang Q, Li H, Han Y, et al. Comprehensive review of web servers and bioinformatics tools for cancer prognosis analysis. Front Oncol. 2020;10:68.
Ohshima K, Fujiya K, Nagashima T, Ohnami S, Hatakeyama K, Urakami K, et al. Driver gene alterations and activated signaling pathways toward malignant progression of gastrointestinal stromal tumors. Cancer Sci. 2019;110:3821–33.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41:D991–5.
Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.
Gene Ontology C. The Gene Ontology resource: enriching a GOld mine. Nucleic Acids Res. 2021;49:D325–34.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49:D545–51.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45:D362–8.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607–13.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Ramos PIP, Arge LWP, Lima NCB, Fukutani KF, de Queiroz ATL. Leveraging user-friendly network approaches to extract knowledge from high-throughput omics datasets. Front Genet. 2019;10:1120.
Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007;9:166–80.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45:W98–102.
Lagarde P, Perot G, Kauffmann A, Brulard C, Dapremont V, Hostein I, et al. Mitotic checkpoints and chromosome instability are strong predictors of clinical outcome in gastrointestinal stromal tumors. Clin Cancer Res. 2012;18:826–38.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102:15545–50.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34:267–73.
Lee EJ, Kang G, Kang SW, Jang KT, Lee J, Park JO, et al. GSTT1 copy number gain and ZNF overexpression are predictors of poor response to imatinib in gastrointestinal stromal tumors. PLoS ONE. 2013;8:e77219.
Chen WB, Cheng XB, Ding W, Wang YJ, Chen D, Wang JH, et al. Centromere protein F and survivin are associated with high risk and a poor prognosis in colorectal gastrointestinal stromal tumours. J Clin Pathol. 2011;64:751–5.
Schneider CA, Rasband WS, Eliceiri KW. NIH image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.
Corless CL, Schroeder A, Griffith D, Town A, McGreevey L, Harrell P, et al. PDGFRA mutations in gastrointestinal stromal tumors: frequency, spectrum and in vitro sensitivity to imatinib. J Clin Oncol. 2005;23:5357–54364.
von Mehren M, Joensuu H. Gastrointestinal stromal tumors. J Clin Oncol. 2018;36:136–43.
Agaimy A, Wunsch PH. Lymph node metastasis in gastrointestinal stromal tumours (GIST) occurs preferentially in young patients < or = 40 years: an overview based on our case material and the literature. Langenbecks Arch Surg. 2009;394:375–81.
Joensuu H, Hohenberger P, Corless CL. Gastrointestinal stromal tumour. Lancet. 2013;382:973–83.
Keung EZ, Fairweather M, Raut CP. The role of surgery in metastatic gastrointestinal stromal tumors. Curr Treat Options Oncol. 2016;17:8.
Chen JS, Yeh CN, Cheng CT, Yen CC, Chen YY, Huang SC, et al. Role of PLK1 signaling pathway genes in gastrointestinal stromal tumors. Oncol Lett. 2018;16:3070–82.
Amirnasr A, Gits CMM, van Kuijk PF, Smid M, Vriends ALM, Rutkowski P, et al. Molecular comparison of imatinib-Naive and resistant gastrointestinal stromal tumors: differentially expressed microRNAs and mRNAs. Cancers. 2019;11:882.
Yan M, Wang C, He B, Yang M, Tong M, Long Z, et al. Aurora-A kinase: a potent oncogene and target for cancer therapy. Med Res Rev. 2016;36:1036–79.
Dematteo RP, Ballman KV, Antonescu CR, Maki RG, Pisters PW, Demetri GD, et al. Adjuvant imatinib mesylate after resection of localised, primary gastrointestinal stromal tumour: a randomised, double-blind, placebo-controlled trial. Lancet. 2009;373:1097–104.
Joensuu H, Eriksson M, Sundby Hall K, Reichardt A, Hartmann JT, Pink D, et al. Adjuvant imatinib for high-risk GI stromal tumor: analysis of a randomized trial. J Clin Oncol. 2016;34:244–50.
Wang-Bishop L, Chen Z, Gomaa A, Lockhart AC, Salaria S, Wang J, et al. Inhibition of AURKA reduces proliferation and survival of gastrointestinal cancer cells with activated KRAS by preventing activation of RPS6KB1. Gastroenterology. 2019;156:662-675.e7.
Otto T, Horn S, Brockmann M, Eilers U, Schuttrumpf L, Popov N, et al. Stabilization of N-Myc is a critical function of Aurora A in human neuroblastoma. Cancer Cell. 2009;15:67–78.
Sun JM, Yang LN, Xu H, Chang B, Wang HY, Yang G. Inhibition of Aurora A promotes chemosensitivity via inducing cell cycle arrest and apoptosis in cervical cancer cells. Am J Cancer Res. 2015;5:1133–45.
Otto T, Sicinski P. Cell cycle proteins as promising targets in cancer therapy. Nat Rev Cancer. 2017;17:93–115.
Yeh CN, Yen CC, Chen YY, Cheng CT, Huang SC, Chang TW, et al. Identification of aurora kinase A as an unfavorable prognostic factor and potential treatment target for metastatic gastrointestinal stromal tumors. Oncotarget. 2014;5:4071–86.
We acknowledge and appreciate Prof. Wenbin Chen for donating the GIST/T1 cell line.
This research was supported by the Zhejiang Provincial Natural Science Foundation of China under Grant No. LBY21H160005.
Ethics approval and consent to participate
This study was approved by the Ethical Committee of the First Affiliated Hospital, Zhejiang University School of Medicine.
Consent for publication
The authors declare no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The expression matrix of the top 50 up-regulated and down-regulated DEGs.
Functional annotation of DEGs based on GSE136755 dataset.
GSEA based on GSE47911 dataset.
Clinicopathological features, AURKA expression and follow-up of GIST patients.
Characteristics of GIST samples in GSE136755 and raw data provided by Lagarde et al. (PMID:22167411).
Raw data from in vitro experiments.
About this article
Cite this article
Cheng, X., Wang, J., Lu, S. et al. Aurora kinase A (AURKA) promotes the progression and imatinib resistance of advanced gastrointestinal stromal tumors. Cancer Cell Int 21, 407 (2021). https://doi.org/10.1186/s12935-021-02111-7