Deciphering key genes involved in cisplatin resistance in kidney renal clear cell carcinoma through a combined in silico and in vitro approach
Services Hospital, Lahore, Pakistan
Sabzazar Medical Center, Lahore, Pakistan
BHU Pourmiana, Attock, Pakistan
Public Health Department, University of Health Sciences, Lahore, Pakistan
Department of Zoology, Shah Abdul Latif University, Khairpur, Pakistan
Department of Animal Sciences, University of Florida, Gainesville, USA
Pharmacology and Toxicology Department College of Pharmacy, Umm Al-Qura University, Makkah, Saudi Arabia
Department of Pharmacology, College of Pharmacy, Jouf University, Aljouf, Saudi Arabia
Regional Laboratory and Blood Bank (Total Laboratory Automation), Riyadh, Saudi Arabia
Health Affairs in Najran, New Najran General Hospital, Najran, Saudi Arabia
Department of Pharmaceutics, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
PARC Arid Zone Research Center, Dera Ismail Khan, Pakistan
Abstract
The low survival rate of Kidney renal clear cell carcinoma (KIRC) patients is largely attributed to cisplatin resistance. Rather than focusing solely on individual proteins, exploring protein-protein interactions could offer greater insight into drug resistance. To this end, a series of in silico and in vitro experiments were conducted to identify hub genes in the intricate network of cisplatin resistance-related genes in KIRC chemotherapy. The genes involved in cisplatin resistance across KIRC were retrieved from the National Center for Biotechnology Information (NCBI) database using search terms as “Kidney renal clear cell carcinoma” and “Cisplatin resistance”. The genes retrieved were analyzed for hub gene identification using the STRING database and Cytoscape tool. Expression and promoter methylation profiling of the hub genes was done using UALCAN, GEPIA, OncoDB, and HPA databases. Mutational, survival, functional enrichment, immune cell infiltration, and drug prediction analyses of the hub genes were performed using the cBioPortal, GEPIA, GSEA, TIMER, and DrugBank databases. Lastly, expression and methylation levels of the hub genes were validated on two cisplatin-resistant RCC cell lines (786-O and A-498) and a normal renal tubular epithelial cell line (HK-2) using two high throughput techniques, including targeted bisulfite sequencing (bisulfite-seq) and RT-qPCR. A total of 124 genes were identified as being associated with cisplatin resistance in KIRC. Out of these genes, MCL1, IGF1R, CCND1, and PTEN were identified as hub genes and were found to have significant (p < 0.05) variations in their mRNA and protein expressions and effects on the overall survival (OS) of the KIRC patients. Moreover, an aberrant promoter methylation pattern was found to be associated with the dysregulation of the hub genes. In addition to this, hub genes were also linked with different cisplatin resistance-causing pathways. Thus, hub genes can be targeted with Alvocidib, Estradiol, Tretinoin, Capsaicin, Dronabinol, Metribolone, Calcitriol, Acetaminophen, Acitretin, Cyclosporine, Azacitidine, Genistein, and Resveratrol drugs. As the pathogenesis of KIRC is complex, targeting hub genes and associated pathways involved in cisplatin resistance could bring a milestone change in the drug discovery and management of drug resistance, which might uplift overall survival among KIRC patients.
Untitled section
Keywords: KIRC, Cisplatin resistance, Chemotherapy, Overall survival
Article notes
Untitled section
Received 2023 Apr 21; Accepted 2023 May 29; Collection date 2023.
Introduction
Kidney renal clear cell carcinoma (KIRC) is a type of kidney cancer that originates from the cells lining the small tubules in the kidney that filter waste from the blood [1]. KIRC is the most common subtype of renal cell carcinoma, accounting for approximately 75% of all cases [2]. The incidence of KIRC has been increasing over the past few decades [3], and it is estimated that there will be approximately 76,080 new cases of kidney cancer in the United States in 2022, with KIRC accounting for the majority of cases [4].
The prognosis for KIRC varies based on factors such as Tumour, Node, Metastasis (TNM) staging, age, and gender, with the patient’s five-year survival rate heavily influenced by the tumor’s pathological grade [5,6]. KIRC patients with an early-stage tumor (grade I or II) tend to have better treatment outcomes and a longer survival period post-surgery, with a survival rate of 80%–90% [7]. Conversely, those diagnosed in the advanced stages of the disease (grade III or IV) have a poorer prognosis, with a survival rate of around 20% due to the high malignancy and probability of metastasis or recurrence [8]. Risk factors for KIRC include smoking, obesity, high blood pressure, and a family history of kidney cancer. KIRC is often asymptomatic in its early stages, and symptoms may not appear until the disease has advanced. Common symptoms of KIRC include blood in the urine, back pain, weight loss, and fatigue [9,10].
The management of KIRC involves a variety of approaches, including surgery, radiotherapy, chemotherapy, and immunotherapy [11–13]. Chemotherapy is considered a foundational aspect of cancer treatment and is frequently used alone or in conjunction with radiation therapy for locally advanced tumors [14]. Cisplatin is a commonly used chemotherapeutic drug for managing KIRC that operates by creating DNA adducts and disrupting cell cycles, leading to cell death [15,16]. For locally advanced tumors (stages III and IV), cisplatin is usually given during concurrent chemoradiation therapy (CRT), either alone or following surgery, or as induction therapy followed by CRT [17]. In cisplatin-based CRT, the drug is administered at a dosage of 100 mg/m2 IV on days 1, 22, and 43, or at a weekly dose of 30-40 mg/m2 IV for 6 to 7 weeks [18]. Despite all this treatment modality, only around 50% of locally advanced KIRC patients respond to therapy [19].
Earlier, several studies reported the development of cisplatin resistance in KIRC patients. In a study by Taniguchi et al. they found that cisplatin resistance was associated with the overexpression of the MDR1 gene in KIRC cells [20]. Similarly, Liu et al. reported that high level of miR-193a-3p expression was associated with cisplatin resistance in KIRC patients [21]. Other studies, such as the one by Sui et al. identified the role of autophagy in cisplatin resistance development in KIRC cells [22]. These findings suggest that cisplatin resistance in KIRC could be due to multiple factors, and the underlying mechanisms need to be further investigated to overcome this challenge.
Moreover, cancer patients who are resistant to cisplatin treatment typically have a bleak prognosis [23]. Proteins are vital macromolecules involved in various cellular operations, and understanding their expression and interactions can provide insight into the complicated molecular pathways involved in drug resistance [24]. Protein-protein interaction (PPI) is an indirect phenomenon that contributes to various cellular functions [25]. While transient PPI controls signaling pathways, permanent PPI forms a protein complex. Approximately 80% of protein functions occur through PPI [26], making the assessment of protein complexes or networks of proteins a more effective approach than targeting a single protein in defeating drug resistance. Furthermore, PPI interactions can help predict the role of previously unexplained proteins.
Hence, this study was initiated for identifying biomarker hub genes through in silico and in vitro methodologies, that are associated with cisplatin resistance in KIRC. In silico and in vitro combined approach for identifying cisplatin resistance-associated hub genes is the integration of computational and experimental methods. The in silico approach involves the utilization of bioinformatics tools and algorithms to analyze large-scale genomic, transcriptomic, and proteomic data sets. On the other hand, the in vitro approach involves the application of experimental methods to validate the identified hub genes, including their expression and methylation levels. The integration approach of in silico and in vitro methodologies can help improve our understanding of cisplatin resistance development in KIRC, leading to the development of more effective therapies.
Materials and Methods
Extraction of cisplatin resistance-associated genes in KIRC
The National Center for Biotechnology Information (NCBI, https://www.ncbi.nlm.nih.gov/) database [27] was searched to identify cisplatin resistance-associated genes in KIRC patients. The search term included “Kidney renal clear cell carcinoma” AND “Cisplatin resistance”. At the end of the search process, a total of 191 items related to Homo sapiens have appeared. Out of the appeared items, a total of 124 cisplatin resistance genes were used in the present work (Fig. 1A). Detail of the search further includes: (Kidney renal clear cell carcinoma [All Fields], (Cisplatin [All Fields], Resistance [All Fields]), “Homo sapiens” [porgn].
mRNA and protein expression profiling of hub genes
UALCAN (https://ualcan.path.uab.edu/analysis-prot.html) [31] is an online database that provides gene expression analysis of cancer data from The Cancer Genome Atlas (TCGA) project. The database is designed to be user-friendly and accessible for researchers and clinicians who want to study the expression of genes and proteins in different cancer types [31]. UALCAN allows users to search for a specific gene or protein of interest and visualize its expression levels across different cancer types, subtypes, and clinical characteristics [31]. UALCAN is a powerful tool for cancer research and personalized medicine, and its open-access nature makes it a valuable resource for the scientific community. UALCAN was used in this work for mRNA and protein expression profiling of the cisplatin resistance-associated hub genes across KIRC samples relative to controls.
mRNA expression validation and survival analysis of hub genes using additional TCGA detests
GEPIA (http://gepia.cancer-pku.cn/) [32] and OncoDB (https://oncodb.org/index.html) [33] are widely used web-based tools that allow researchers to perform gene expression analysis on different types of cancer as well as normal tissues. These tools provide an interactive interface for users to explore and visualize gene expression patterns, differential expression analysis, and survival analysis. These databases include RNA sequencing data from The Cancer Genome Atlas, providing comprehensive resource for cancer research. The tool also integrates functional enrichment analysis, which allows users to identify enriched gene ontology terms and biological pathways associated with differentially expressed genes. In conclusion, both GEPIA and OncoDB tools have become essential tools for cancer research, providing researchers with a comprehensive platforms to analyze gene expression patterns, identify potential therapeutic targets, and develop personalized cancer treatments. In this work, the GEPIA and OncoDB databases were used for the expression validation analysis of the hub genes across KIRC samples relative to controls. Moreover, GEPIA database was further utilized for survival analysis as well.
Subcellular localization and protein expression validation of hub genes
The Human Protein Atlas (HPA, https://www.proteinatlas.org/) database is an extensive and accessible repository of information on the human proteome [34]. It comprises a comprehensive collection of information on tissue and cell-specific protein expression patterns in humans, providing an invaluable resource for researchers and scientists around the world [34]. The database leverages a wide range of cutting-edge technologies, including high-throughput proteome analysis, antibody-based profiling, and molecular-imaging techniques. Additionally, the HPA offers interactive tools and resources for data visualization and analysis, enabling users to explore and interpret the data in various ways. In this work, the HPA database was used to identify the subcellular localization of the proteins encoded by the hub genes in KIRC cells as well as to validate hub gene expression at the protein level across KIRC samples relative to controls based on the immunohistochemical images.
Methylation analysis of hub genes
Mutational analysis, mutation based survival analysis, and co-express gene analysis of hub genes
cBioPortal (https://www.cbioportal.org/) is a user-friendly, open-access platform that enables researchers and clinicians to perform complex genomic analysis on cancer data [35]. It offers visualization and analysis tools for studying complex cancer genomics data, allowing users to explore and analyze large-scale clinical datasets with ease. The database includes genomic data from over 270 cancer studies, including more than 30,000 patients across 20 different types of cancer [35]. The platform incorporates multiple tools, such as data visualization, differential expression analysis, and gene set enrichment analysis, that enable researchers to query the database in various ways. The ability to stratify data by patient characteristics and cancer subtype greatly enhances the potential for personalized cancer treatment. cBioPortal provides an invaluable resource of genomic data for researchers and clinicians, helping revolutionize cancer diagnosis and treatment. In this work, this tool was used for the mutational analysis, mutation-based survival analysis, and co-express gene analysis of hub genes in TCGA KIRC samples.
Functional enrichment analysis
The Gene Ontology (GO) analysis provides functional annotation of the gene(s) of interest. While KEGG (Kyoto Encyclopedia of Genes and Genomes) provides an interpretation of the user-defined genes in biological pathways [36]. The GO and KEGG analyses of the hub and hub gene-enriched genes were performed using the GSEA program [37].
Immune cell infiltration analysis
The TIMER database (http://timer.cistrome.org/) is a comprehensive tool that allows for the estimation of immune cells in a variety of tumor types [38]. The database was developed by researchers at the University of Texas MD Anderson Cancer Center and provides a wealth of information on the immune status of tumors, including the abundance of various immune cells, such as T cells, B cells, natural killer cells, and dendritic cells [38]. This information can be used to better understand the immune response to tumors and develop more effective cancer treatments. In this research, levels of immune cell infiltration in KIRC were plotted against hub gene expression.
miRNA network analysis
The ENCORI platform (https://rna.sysu.edu.cn/encori/) was developed for analyzing the miRNA-mRNA, RNA-RNA, and lncRNA-miRNA networks across various types of cancer [39]. In this investigation, the ENCORI database was used to create the miRNA network of the identified hub genes [39].
Hub genes’ drug prediction analysis
We used the DrugBank database (https://go.drugbank.com/) to uncover a variety of drugs associated with the identified hub genes because we believe that the identified hub genes could be promising therapeutic targets. This database compiles information on chemotherapeutic drugs that target hub genes from a variety of reliable sources [40].
In vitro validation of cisplatin resistance-inducing genes via targeted bisulfide-seq and RT-qPCR analyses
Cell lines: Human RCC cell lines (786-O and A-498) and a normal renal tubular epithelial cell line (HK-2) were purchased from the American Type Culture Collection (ATCC, USA) and cultivated in accordance with the manufacturer’s instructions.
Construction of cisplatin resistance cell strains: Cisplatin resistance of cells was induced using a previously established pulse therapy technique [41] combined with continuous stepwise exposure to the drug. Finally, the cisplatin resistance cells were acquired at a cisplatin concentration of 0.5 μg/mL.
Total RNA and DNA extraction: Total RNA extraction from both KIRC and normal cell lines was done by isopycnic centrifugation as described previously [42]. The extracted RNA was then processed for the DNA digestion step of incubation with RNase-free DNase I (Roche, Germany) at 37°C for 15 min. The DNA extraction was done following the organic method [43]. The quality of the extracted RNA and DNA was checked by a 2100 Bioanalyzer (Agilent Technologies, Germany).
Targeted bisulfite-seq and RT-qPCR analyses: DNA samples were sent to Beijing Genomics Institute (BGI) company for targeted bisulfite-seq analysis. Following targeted bisulfite-seq analyses, methylation values were normalized as beta values. The obtained beta values against hub genes in RCC and a normal control cell line were compared to identify differences in the methylation levels.
RT-qPCR validation analysis: The specific protocols are as follows: First, the PrimeScript™ RT reagent kit (Takara, Japan) was used for reverse transcription of the extracted RNA from cisplatin-resistant RCC cell lines (786-O and A-498) and a normal control cell line (HK-2) into complementary DNA. Then, the RT-qPCR was carried out on an ABI ViiA 7 Real Time PCR System (Thermo Fisher, USA) with a SuperReal SYBR Green Premix Plus (Tiangen Biotech, China) as a fluorescent dye. GAPDH was chosen as the internal reference in the present study. All the experiments were done in triplicate independently. All the primers for each hub gene are shown as following:
GAPDHF 5-ACCCACTCCTCCACCTTTGAC-3′,
GAPDHR 5′-CTGTTGCTGTAGCCAAATTCG-3′ [44].
MCL1F 5′-GGACATCAAAAACGAAGACG-3′,
MCL1R 5′- GCAGCTTTCTTGGTTTATGG-3′ [45].
IGF1RF 5′-GGCACAATTACTGCTCCAAAGAC-3′,
IGF1RR 5′-CAAGGCCCTTTCTCCCCAC-3′ [46].
CCND1F 5′-CGATGCTGGAGGTCTGCGA-3′,
CCND1R 5′-AGAGGCCACGAACATGCAAG-3′ [47].
PTENF 5′-TGGCATACACCAAATATAAGAGC-3′,
PTENR 5′-TCCCTTATCAGATACATGACTTTCA-3′ [48].
Statistics details
For enrichment analysis, we used Fisher’s Exact test for computing statistical difference [49]. Correlational analyses were carried out using the Pearson method. The 2−ΔΔCt method was employed to evaluate the relative expression of each hub gene through RT-qPCR [50]. Log-rank test was used to compare survival curves between normal control and KIRC patients. Cox proportional hazards model was utilized to study the relationship between survival time and the predictive variable [51]. For comparisons, a student t-test was adopted in the current study. All the analyses were carried out in R version 3.6.3 software.
Results
Extraction of cisplatin resistance-associated genes, PPI network, module analysis, and hub genes recognition
After a brief literature search, a total of 124 genes were found associated with cisplatin resistance across KIRC in Homo sapiens. The obtained 124 genes were evaluated for constructing the PPI network through the STRING database. As highlighted in Fig. 1A, the constructed PPI consisted of 124 nodes, where nodes and edges represent proteins and protein–protein associations, respectively. Further, the constructed PPI was analyzed using the MCODE plugin application of the Cytoscape software to identify a significant module consisting of 17 nodes (Fig. 1B). Then, the identified significant module was subjected to Cytohubba analysis for the identification of hub genes based on the degree method. Cytohubba analysis revealed MCL1 (myeloid cell leukemia-1), Insulin-like growth factor 2 receptor (IGF1R), CCND1 (Cyclin D1), and Phosphatase and TENsin homolog deleted on chromosome 10 (PTEN) genes were the hub genes having the highest degree scores (Fig. 1C).
Expression and correlation analysis of MCL1, IGF1R, CCND1, and PTEN with different clinical variables using the TCGA dataset
The TCGA dataset was used in conjunction with the UALCAN tool to analyze the expressions of MCL1, IGF1R, CCND1, and PTEN. The results showed a significant up-regulation of MCL1 and CCND1, while IGF1R and PTEN showed a significant down-regulation at both mRNA and protein levels in KIRC samples compared to controls (Figs. 2A–2C). This difference was statistically significant (p < 0.05).
Furthermore, we also investigated the potential implications of dysregulation of these genes in KIRC patients with different clinicopathological parameters. Information on cancer stage, race, gender, and age was retrieved from UALCAN, and these parameters were analyzed in relation to mRNA expression of the hub genes in the KIRC cohort. Results revealed a significant up-regulation of MCL1 and CCND1, and a significant down-regulation of IGF1R and PTEN in KIRC patients with varying cancer stages, races, genders, and age groups compared to normal controls (Fig. 3).
Expression verification and survival analyses of MCL1, IGF1R, CCND1, and PTEN
In order to enhance the validity of our findings, we incorporated two additional databases, namely GEPIA and OncoDB, to perform expression verification analysis. For this purpose, we evaluated the expression levels of the MCL1, IGF1R, CCND1, and PTEN, as well as their impact on survival in KIRC and control tissues (Fig. 4). Our results indicated that MCL1 and CCND1 exhibited significantly higher mRNA expression levels (p < 0.05) while IGF1R and PTEN mRNA levels were significantly lower (p < 0.05) in KIRC samples compared to healthy samples (Figs. 4A and 4B). Furthermore, we employed the GEPIA database to investigate the hub genes’ impact on overall survival (OS), and the analysis revealed a significant (p < 0.05) correlation between dysregulated MCL1, IGF1R, CCND1, PTEN, and poor OS for KIRC patients (Fig. 4C).
Subcellular localization and protein expression validation analyses of MCL1, IGF1R, CCND1, and PTEN
The location of MCL1, IGF1R, CCND1, and PTEN within KIRC cells was determined using HPA. MCL1 was found to be primarily located in mitochondria (Fig. 5A), IGF1R was localized in the endoplasmic reticulum (Fig. 5A), CCND1 was enriched in the nucleoplasm (Fig. 5A), and PTEN was present in both the cytosol and nucleoplasm (Fig. 5A). The expressions of these proteins were confirmed in KIRC samples through immunohistochemistry (IHC) analysis using data obtained from HPA. The IHC results revealed that MCL1 and CCND1 had higher expressions (staining = medium and high) in KIRC samples (Fig. 5B) compared to control samples (staining = not detected and low) (Fig. 5B). Conversely, the expressions of IGF1R and PTEN were lower (staining = medium and not detected) in KIRC samples (Fig. 5B) compared to control samples (staining = high and medium) (Fig. 5B).
Effect of DNA mutations and promoter methylation on MCL1, IGF1R, CCND1, and PTEN dysregulation and KIRC patients’ survival
The dysregulation of MCL1, IGF1R, CCND1, and PTEN expressions was found to be correlated with various clinical parameters of KIRC and the worst OS. Therefore, exploring the potential regulatory mechanisms involved in the overexpression of these hub genes could have clinical significance. Firstly, we used the cBioPortal database to identify genetic mutations in the hub genes in the TCGA KIRC cohort. Low genetic alteration frequencies were detected in MCL1 (0.7%), IGF1R (0%), and CCND1 (2.5%), indicating their lesser involvement in hub gene expression regulation (Fig. 6A). PTEN had a higher genetic alteration frequency of 5%, which was linked to its down-regulation in KIRC patients (Fig. 6A). The top co-expressed genes were identified through cBioPortal, which included MCL1-WEE1, IGF1R-ARPC3, CCND1-AVPRB1, and PTEN-PTENP1 (Fig. 6B). Genetic alterations in the hub genes were also associated with the worst OS of KIRC patients (Fig. 6C). Finally, promoter methylation levels of the hub genes in normal control and KIRC samples were examined using the UALCAN and OncoDB databases. We found that the promoter methylation levels of MCL1 and CCNB1 were lower, while the promoter methylation levels of IGF1R and PTEN were higher in KIRC samples than in normal tissues (Fig. 7). Thus, abnormal promoter DNA methylation levels may lead to aberrant expression of MCL1, IGF1R, CCND1, and PTEN in KIRC.
GO and biological pathways analysis
The GSEA program was used to conduct GO and KEGG analyses of the identified hub genes. We obtained GO functional enrichments of the MCL1, IGF1R, CCND1, and PTEN genes with a p-value < 0.05. GO and KEGG pathways in KIRC. In the CC, “Insulin receptor complex, Alphav-beta integrin-IGF-1-IGF1R complex, Bcl-2 family protein complex, and Arp2/3 protein complex”, etc., were significantly associated with the analyzed genes (Fig. 8A). Concerning MF, the “Re-entry into mitotic cell cycle, Reg. of systemic arterial blood pressure by vasopressin, and Pos, Reg. phospholipase A2 activity”, etc., were closely associated the hub genes (Fig. 8B). In BP, some vital functions including “Insulin-activated receptor activity, Insulin-like growth factor-activated receptor activity, Insulin Binding, and BH3 binding”, etc., were significantly associated with hub genes (Fig. 8C). Hub genes and their co-expressed genes associated KEGG pathways include “Melanoma, Thyroid cancer, Glioma, Bladder cancer, and Prostate cancer”, etc. (Fig. 8D).
Immune cell analysis of MCL1, IGF1R, CCND1, and PTEN
Next, TIMER was utilized to uncover associations between the expression of MCL1, IGF1R, CCND1, and PTEN genes, and immune cell influx (CD8+ T, CD4+ T, and macrophages). A significant positive correlation (p < 0.05) was observed between CD8+ T and CD4+ T cell infiltration and the expression levels of MCL1, IGF1R, CCND1, and PTEN genes (Fig. 9), while a negative correlation was observed with macrophages (Fig. 9).
miRNA network of the MCL1, IGF1R, CCND1, and PTEN
Additionally, we established a miRNA-mRNA co-regulatory network of the aforementioned genes through ENCORI and Cytoscape. The constructed miRNA-mRNA co-regulatory network had a total of 454 miRNAs and 4 mRNAs. Based on the degree method, Cytohubba analysis further identified a miRNA (has-mir-17-5p) as a potential inducer of KIRC, as it targets all hub genes simultaneously. These findings suggest that the identified axis of has-mir-27a-3p and MCL1, IGF1R, CCND1, and PTEN genes may have a role in drug resistance development across KIRC (Fig. 10).
Drug prediction analysis of the MCL1, IGF1R, CCND1, and PTEN
For KIRC patients suffering from cisplatin drug resistance, the use of alternative drugs for medical treatment is the first option. Therefore, a selection of appropriate alternate candidate drugs is required. In the current study, via the DrugBank database, we explored some potential drugs, that can reverse the gene expressions of MCL1, IGF1R, CCND1, and PTEN hub genes. Alvocidib, Estradiol, Tretinoin, Capsaicin, Dronabinol, Metribolone, Calcitriol, Acetaminophen, Acitretin, Cyclosporine, Azacitidine, Genistein, and Resveratrol drugs (Table 1) could be useful to target MCL1, IGF1R, CCND1, and PTEN hub genes to counteract cisplatin resistance once clinical significance is established via preclinical/clinical studies.
| Sr. No. | Hub gene | Drug name | Effect | Reference | Group |
|---|---|---|---|---|---|
| 1 | MCL1 | Alvocidib | Decrease expression of MCL1 mRNA | A20630 | Approved |
| Estradiol | A21152 | ||||
| Tretinoin | A20405 | ||||
| 2 | IGF1R | Capsaicin | Increase expression of IGF1R mRNA | A21513 | Approved |
| Dronabinol | A22085 | ||||
| Metribolone | A21310 | ||||
| Calcitriol | A22301 | ||||
| 3 | CCND1 | Acetaminophen | Decrease expression of CCND1 mRNA | A20420 | Approved |
| Acitretin | A20453 | ||||
| Cyclosporine | A20661 | ||||
| 4 | PTEN | Estradiol | Increase expression of PTEN mRNA | A21179 | Approved |
| Azacitidine | A20985 | ||||
| Genistein | A22790 | ||||
| Resveratrol | A22790 | ||||
| Tretinoin | A24453 |
In vitro validation of cisplatin resistance-inducing genes via targeted bisulfide-seq
In the current study, by performing targeted bisulfite-seq analyses of 2 cisplatin-resistant RCC cell lines, including 786-O and A-498, and the normal renal tubular epithelial cell line HK-2, methylation levels of the MCL1, IGF1R, CCND1, and PTEN were validated. Methylation levels of these genes were validated using beta values. As shown in Fig. 11A, it was noticed that beta values of MCL1 and CCND1 were lower, while beta values of IGF1R and PTEN were higher in the cisplatin-resistant RCC cell lines (786-O and A-498) as compared to the normal control cell line (HK-2).
RT-qPCR validation analysis of MCL1, IGF1R, CCND1, and PTEN
To further validate the aforementioned bioinformatics analysis, the mRNA expression levels of these hub genes (MCL1, IGF1R, CCND1, and PTEN) were obtained by RT-qPCR experiment in cisplatin resistant RCC cell lines (786-O and A-498) and a normal control cell line (HK-2). As suggested in Fig. 11B, there were notable differences between the cisplatin resistant 786-O and A-498 cell lines and the HK-2 cell line in all four hub genes, including MCL1, IGF1R, CCND1, and PTEN (Fig. 11B). Interestingly, MCL1 and CCND1 genes had up-regulation, while IGF1R and PTEN had notable down-regulation in the cisplatin-resistant RCC cell line relative to the normal control cell line (HK-2) (Fig. 11B) as previously predicted by the TCGA dataset analysis.
Discussion
Cisplatin resistance refers to the ability of cancer cells to become resistant to the effects of cisplatin, to treat a variety of cancers, including KIRC [52]. Cisplatin resistance can develop due to a variety of mechanisms, including dysregulation of important genes, and changes in cellular signaling pathways [53]. According to scientific research, proteins with homologous similarities from phylogenetic evolution tend to conserve their protein domains, making them essentially alike [54]. This means that domain-conserved proteins may interact with each other when influenced by external factors such as xenobiotics [55]. In drug resistance scenarios, targeting a single protein may activate or inhibit the function of homologous proteins to compensate for the loss of the targeted protein’s function [56]. However, there is a need for more clarification and declaration of the protein-protein interactions among the trigger genes linked to cisplatin resistance in KIRC. To this end, the current study sets out to investigate the protein-protein interaction among the genes and identify some important genes linked to cisplatin resistance in KIRC. The identified hub genes will help in understanding the molecular mechanisms underlying cisplatin resistance and developing new therapeutic strategies for KIRC patients. Therefore, this study was initiated for identifying cisplatin resistance-associated biomarker hub genes across KIRC through in silico and in vitro methodologies.
To do so, initially, 124 cisplatin resistance-linked genes were identified from the published literature. Then, PPI of the identified genes was constructed and screened out via Cytoscape software to identify hub genes. Based on the degree method, the MCL1, IGF1R, CCND1, and PTEN genes were identified as the most prominent drug resistance-associated genes in KIRC.
The MCL1 gene encodes for a protein called myeloid cell leukemia 1, which plays a crucial role in regulating programmed cell death, or apoptosis [57,58]. In normal cells, MCL1 helps maintain a balance between cell growth and cell death [59]. However, in cancer cells, MCL1 expression is often up-regulated, leading to increased cell survival and resistance to chemotherapy drugs. This makes MCL1 an attractive target for cancer therapy development [60]. Several studies have shown that targeting MCL1 can induce apoptosis in cancer cells and sensitize them to chemotherapy [61,62]. Moreover, MCL1 has been implicated in resistance to a range of cancer therapies, including those targeting the BCL-2 family of proteins, which are also involved in apoptosis regulation [63,64]. Therefore, there is significant interest in developing drugs that can inhibit MCL1 expression or function as a potential cancer treatment strategy. In summary, the MCL1 gene plays a critical role in cancer by promoting cell survival and resistance to chemotherapy. Targeting MCL1 represents a promising therapeutic approach for overcoming resistance to cancer treatment and improving patient outcomes.
The IGF1R (Insulin-like Growth Factor 1 Receptor) gene encodes for a transmembrane receptor protein that binds to insulin-like growth factors (IGFs) and plays a key role in cell growth and survival [65,66]. In normal cells, the IGF1R signaling pathway is tightly regulated, but in cancer cells, it can become dysregulated, leading to uncontrolled cell growth and proliferation [67]. The overexpression of IGF1R has been observed in a variety of cancers, including breast, lung, prostate, and colon cancer [68], making it an attractive target for cancer therapy. Studies have shown that blocking the IGF1R pathway can inhibit cancer cell growth and promote apoptosis [69,70]. In preclinical studies, IGF1R inhibitors have demonstrated promising anti-tumor activity, particularly in combination with other cancer therapies [68]. Additionally, the IGF1R pathway has been linked to resistance to several cancer treatments, including chemotherapy and targeted therapies, highlighting its potential as a therapeutic target for overcoming drug resistance [68]. Despite initial enthusiasm, clinical trials of IGF1R inhibitors have shown mixed results, with some trials showing limited efficacy and others showing promising results in specific patient populations [71]. Overall, the IGF1R gene is an important player in cancer pathogenesis and a promising target for cancer therapy development.
The CCND1 gene, also known as cyclin D1, plays a crucial role in the cell cycle by regulating the G1/S transition [72,73]. However, overexpression of CCND1 has been linked to drug resistance in cancer cells [74]. Studies have shown that increased expression of CCND1 in cancer cells can lead to the activation of various drug-resistance pathways, such as the up-regulation of drug transporters and anti-apoptotic proteins [75]. Furthermore, inhibition of CCND1 has been found to sensitize cancer cells to chemotherapy and radiation therapy [76]. Therefore, targeting the CCND1 gene may offer a promising strategy to overcome drug resistance in cancer treatment. Overall, the CCND1 gene appears to be a key player in drug resistance in cancer cells, representing a valuable target for the development of novel cancer treatments.
The PTEN gene is a tumor suppressor gene that plays a critical role in regulating cell growth and proliferation [77]. PTEN is often mutated or deleted in many different cancers, leading to the loss of its function and uncontrolled cell growth [78]. Recent studies have also shown that the loss of PTEN function is associated with drug resistance in cancer cells [79,80]. In particular, PTEN loss has been linked to resistance to targeted therapies such as tyrosine kinase inhibitors and immune checkpoint inhibitors [81]. The mechanisms by which PTEN loss confers drug resistance are not yet fully understood but may involve altered signaling pathways, changes in the tumor microenvironment, and alterations in DNA repair mechanisms. Understanding the role of PTEN in drug resistance may lead to the development of novel therapeutic strategies to overcome this resistance in cancer patients.
We next analyzed and validated the mRNA and protein expression levels of MCL1, IGF1R, CCND1, and PTEN using UALCAN, GEPIA, OncoDB, and HPA databases on TCGA datasets. Our research took into account different factors such as cancer stages, races, genders, ages, and pathological data. Analysis of the UALCAN, GEPIA, and OncoDB databases revealed a significant up-regulation of MCL1 and CCND1 mRNA expressions, and a significant down-regulation of IGF1R and PTEN mRNA expressions in KIRC samples compared to controls. We also found that the protein expression of MCL1, IGF1R, CCND1, and PTEN in KIRC tissue samples was consistent with their corresponding mRNA expression levels using the HPA database and immunohistochemical staining. Therefore, we infer that the high expressions of MCL1 and CCND1 and low expressions of IGF1R and PTEN may play critical roles in developing cisplatin drug resistance across KIRC. The analysis of promoter methylation has revealed that an abnormal pattern of methylation significantly contributes to the dysregulation of MCL1, IGF1R, CCND1, and PTEN expression in KIRC patients compared to normal controls. MCL1 methylation has been associated with poor overall survival in KIRC patients [82], while IGF1R and CCND1 methylation have been linked to tumor aggressiveness [75,83]. PTEN methylation, on the other hand, has been associated with the loss of tumor suppressor function in KIRC [84]. Subsequently, we utilized the Kaplan–Meier plotter and discovered that a high expression of MCL1 and CCND1, along with low expressions of IGF1R and PTEN, could lead to a shorter survival time for KIRC patients. Therefore, the expression levels of MCL1, IGF1R, CCND1, and PTEN can serve as an independent risk factor for a poor prognosis of KIRC.
The TIMER analysis revealed a significant correlation between the expression of MCL1, IGF1R, CCND1, and PTEN and the presence of CD8+ T cells, CD4+ T cells, and macrophages. The CD8+ T and CD4+ T cells are the most abundant types of T lymphocytes in tumor microenvironments and capable of eliminating tumor cells, their differentiation into dysfunctional states may render them ineffective at responding to immunotherapy [85,86]. This suggests that the dysregulation of MCL1, IGF1R, CCND1, and PTEN expression could be responsible for KIRC aggressiveness and resistance to immunotherapy by impairing the functionality of CD8+ T and CD4+ T cells. Furthermore, the miRNA hsa-mir-17-5p was found to be a common modulator of MCL1, IGF1R, CCND1, and PTEN hub gene expression in KIRC patients through miRNA network analysis. The cisplatin resistance-associated role of hsa-mir-17-5p is well-established in various types of cancer, including colorectal cancer, head and neck cancer, and leukemia [87,88]. However, this study is the first to report the cisplatin resistance-associated role of hsa-mir-17-5p as a modulator of MCL1, IGF1R, CCND1, and PTEN hub genes. Lastly, KEGG analysis further validated the connection of MCL1, IGF1R, CCND1, and PTEN hub genes with different cisplatin resistance-associated pathways, such as Endocrine resistance, AMPK signaling pathway, and the Foxo signaling pathway [89–91].
Conclusion
This study revealed that MCL1, IGF1R, CCND1, and PTEN, are the central genes implicated in cisplatin resistance in KIRC. These genes demonstrated significant expression variability and had an impact on the OS of KIRC patients. These genes were also found modulating multiple cisplatin-associated pathways. Targeting MCL1, IGF1R, CCND1, and PTEN genes with different other drugs explored via the present study, such as palbociclib, methotrexate, bortezomib, fluorouracil, sorafenib, dasatinib, carboplatin, paclitaxel, gemcitabine, imatinib, doxorubicin, and vorinostat, may have the potential to overcome cisplatin drug resistance in KIRC.
Acknowledgments
None.
Funding Statement
The authors would like to extend their sincere appreciation to the Researchers Supporting Project Number (RSPD2023R986), King Saud University, Riyadh, Saudi Arabia.
Availability of Data and Materials
No datasets were generated or analyzed during the current study.
Ethics Approval
Not applicable.
Conflicts of Interest
The authors declare that they have no conflicts of interest to report regarding the present study.
References
Untitled section
References
- 1.Xiong, X., Chen, C., Yang, J., Ma, L., Wang, X.et al. (2022). Characterization of the basement membrane in kidney renal clear cell carcinoma to guide clinical therapy. Frontiers in Oncology , 12, 1024956.;
- 2.Rao, H., Li, X., Liu, M., Liu, J., Feng, W.et al. (2021). Multilevel regulation of β-catenin activity by SETD2 suppresses the transition from polycystic kidney disease to clear cell renal cell carcinoma SETD2 deficiency promotes the transition from PKD to ccRCC. Cancer Research , 81(13), 3554–3567;
- 3.Capitanio, U., Montorsi, F. (2016). Renal cancer. The Lancet , 387(10021), 894–906.
- 4.Xing, Q., Luan, J., Liu, S., Ma, L., Wang, Y. (2022). Six RNA binding proteins (RBPs) related prognostic model predicts overall survival for clear cell renal cell carcinoma and is associated with immune infiltration. Bosnian Journal of Basic Medical Sciences , 22(3), 435–452;
- 5.Huang, X., Shu, C., Chen, L., Yao, B. (2018). Impact of sex, body mass index and initial pathologic diagnosis age on the incidence and prognosis of different types of cancer. Oncology Reports , 40(3), 1359–1369;
- 6.Liu, Y., Wu, D., Chen, H., Yan, L., Xiang, Q.et al. (2023). Construction and verification of a novel prognostic risk model for kidney renal clear cell carcinoma based on immunity-related genes. Frontiers in Genetics , 14, 1107294;
- 7.Tsui, K. H., Shvarts, O., Smith, R. B., Figlin, R., de Kernion, J. B.et al. (2000). Renal cell carcinoma: Prognostic significance of incidentally detected tumors. The Journal of Urology , 163(2), 426–430;
- 8.Sorbellini, M., Kattan, M. W., Snyder, M. E., Reuter, V., Motzer, R.et al. (2005). A postoperative prognostic nomogram predicting recurrence for patients with conventional clear cell renal cell carcinoma. The Journal of Urology , 173(1), 48–51;
- 9.Dhote, R., Pellicer-Coeuret, M., Thiounn, N., Debre, B., Vidal-Trecan, G. (2000). Risk factors for adult renal cell carcinoma: A systematic review and implications for prevention. BJU International , 86(1), 20–27;
- 10.Sial, N., Ahmad, M., Hussain, M., Iqbal, M. J., Hameed, Y.et al. (2021). CTHRC1 expression is a novel shared diagnostic and prognostic biomarker of survival in six different human cancer subtypes. Scientific Reports , 11(1), 19873;
- 11.Smyth, M. J., Ngiow, S. F., Ribas, A., Teng, M. W. (2016). Combination cancer immunotherapies tailored to the tumour microenvironment. Nature Reviews Clinical Oncology , 13(3), 143–158;
- 12.Sial, N., Saeed, S., Ahmad, M., Hameed, Y., Rehman, A.et al. (2021). Multi-omics analysis identified TMED2 as a shared potential biomarker in six subtypes of human cancer. International Journal of General Medicine , 14, 7025–7042;
- 13.Zhang, G., Chen, X., Fang, J., Tai, P., Chen, A.et al. (2022). Cuproptosis status affects treatment options about immunotherapy and targeted therapy for patients with kidney renal clear cell carcinoma. Frontiers in Immunology , 13, 954440;
- 14.Krishnan, S., Chadha, A. S., Suh, Y., Chen, H. C., Rao, A.et al. (2016). Focal radiation therapy dose escalation improves overall survival in locally advanced pancreatic cancer patients receiving induction chemotherapy and consolidative chemoradiation. International Journal of Radiation Oncology, Biology, Physics , 94(4), 755–765;
- 15.Brown, A., Kumar, S., Tchounwou, P. B. (2019). Cisplatin-based chemotherapy of human cancers. Journal of Cancer Science and Clinical Therapeutics , 11 (4), 97.
- 16.Ranasinghe, R., Mathai, M. L., Zulli, A. (2022). Cisplatin for cancer therapy and overcoming chemoresistance. Heliyon , 8(9), e10608;
- 17.Haddad, R., O’Neill, A., Rabinowits, G., Tishler, R., Khuri, F.et al. (2013). Induction chemotherapy followed by concurrent chemoradiotherapy (sequential chemoradiotherapy) versus concurrent chemoradiotherapy alone in locally advanced head and neck cancer (PARADIGM): A randomised phase 3 trial. The Lancet Oncology , 14(3), 257–264;
- 18.Negi, P., Kingsley, P. A., Srivastava, H., Sharma, S. K. (2016). Three weekly versus weekly cisplatin as radiosensitizer in head and neck cancer: A decision dilemma. Asian Pacific Journal of Cancer Prevention , 17(4), 1617–1623;
- 19.Peiretti, M., Zapardiel, I., Zanagnolo, V., Landoni, F., Morrow, C.et al. (2012). Management of recurrent cervical cancer: A review of the literature. Surgical Oncology , 21(2), e59–e66;
- 20.Taniguchi, K., Wada, M., Kohno, K., Nakamura, T., Kawabe, T.et al. (1996). A human canalicular multispecific organic anion transporter (cMOAT) gene is overexpressed in cisplatin-resistant human cancer cell lines with decreased drug accumulation. Cancer Research , 56(18), 4124–4129;
- 21.Liu, L., Li, Y., Liu, S., Duan, Q., Chen, L.et al. (2017). Downregulation of miR-193a-3p inhibits cell growth and migration in renal cell carcinoma by targeting PTEN. Tumour Biology , 39(6). 10.1177/1010428317711951;
- 22.Sui, X., Chen, R., Wang, Z., Huang, Z., Kong, N.et al. (2013). Autophagy and chemotherapy resistance: A promising therapeutic target for cancer treatment. Cell Death & Disease , 4(10), e838.
- 23.Yang, Q. K., Chen, T., Wang, S. Q., Zhang, X. J., Yao, Z. X. (2020). Apatinib as targeted therapy for advanced bone and soft tissue sarcoma: A dilemma of reversing multidrug resistance while suffering drug resistance itself. Angiogenesis , 23(3), 279–298;
- 24.Morris, R., Black, K. A., Stollar, E. J. (2022). Uncovering protein function: From classification to complexes. Essays in Biochemistry , 66(3), 255–285;
- 25.Milroy, L. G., Grossmann, T. N., Hennig, S., Brunsveld, L., Ottmann, C. (2014). Modulators of protein-protein interactions. Chemical Reviews , 114(9), 4695–4748;
- 26.Niu, Y., Otasek, D., Jurisica, I. (2010). Evaluation of linguistic features useful in extraction of interactions from PubMed; application to annotating known, high-throughput and predicted interactions in I2D. Bioinformatics , 26(1), 111–119;
- 27.Wheeler, D. L., Barrett, T., Benson, D. A., Bryant, S. H., Canese, K.et al. (2007). Database resources of the national center for biotechnology information. Nucleic Acids Research , 35(suppl_1), D5–D12;
- 28.Mering, C. V., Huynen, M., Jaeggi, D., Schmidt, S., Bork, P.et al. (2003). STRING: A database of predicted functional associations between proteins. Nucleic Acids Research , 31(1), 258–261.
- 29.Demchak, B., Hull, T., Reich, M., Liefeld, T., Smoot, M.et al. (2014). Cytoscape: The network visualization tool for GenomeSpace workflows. F1000Research , 3, 151;
- 30.Chin, C. H., Chen, S. H., Wu, H. H., Ho, C. W., Ko, M. T.et al. (2014). cytoHubba: Identifying hub objects and sub-networks from complex interactome. BMC Systems Biology , 8(4), 1–7.
- 31.Chandrashekar, D. S., Bashel, B., Balasubramanya, S. A. H., Creighton, C. J., Ponce-Rodriguez, I.et al. (2017). UALCAN: A portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia , 19(8), 649–658;
- 32.Tang, Z., Li, C., Kang, B., Gao, G., Li, C.et al. (2017). GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Research , 45(W1), W98–W102;
- 33.Tang, G., Cho, M., Wang, X. (2022). OncoDB: An interactive online database for analysis of gene expression and viral infection in cancer. Nucleic Acid Research , 50(D1), D1334–D1339.
- 34.Thul, P. J., Lindskog, C. (2018). The human protein atlas: A spatial map of the human proteome. Protein Science , 27(1), 233–244. 10.1002/pro.3307;
- 35.Gao, J., Aksoy, B. A., Dogrusoz, U., Dresdner, G., Gross, B.et al. (2013). Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Science , 6(269), pl1.
- 36.Zhang, J., Xing, Z., Ma, M., Wang, N., Cai, Y. D.et al. (2014). Gene ontology and KEGG enrichment analyses of genes related to age-related macular degeneration. Biomed Research International , 2015, 450386.
- 37.Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L.et al. (2005). Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America , 102(43), 15545–15550;
- 38.Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q.et al. (2017). TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Research , 77(21), e108–e110;
- 39.Huang, D. P., Zeng, Y. H., Yuan, W. Q., Huang, X. F., Chen, S. Q.et al. (2021). Bioinformatics analyses of potential miRNA-mRNA regulatory axis in HBV-related hepatocellular carcinoma. International Journal of Medical Sciences , 18(2), 335–346;
- 40.Freshour, S. L., Kiwala, S., Cotto, K. C., Coffman, A. C., McMichael, J. F.et al. (2021). Integration of the drug-gene interaction database (DGIdb 4.0) with open crowdsource efforts. Nucleic Acids Research , 49(D1), D1144–D1151;
- 41.Liu, H., Zhang, Y. (2020). Effects of miR-222 on cisplatin resistance of renal cancer cell strains and related mechanisms. Indian Journal of Biochemistry and Biophysics (IJBB) , 57(4), 382–388.
- 42.Feng, Z., Li, Y., McKnight, K. L., Hensley, L., Lanford, R. E.et al. (2015). Human pDCs preferentially sense enveloped hepatitis A virions. The Journal of Clinical Investigation , 125(1), 169–176;
- 43.Ghatak, S., Muthukumaran, R. B., Nachimuthu, S. K. (2013). A simple method of genomic DNA extraction from human samples for PCR-RFLP analysis. Journal of Biomolecular Techniques , 24(4), 224–231;
- 44.Jafri, H., Mushtaq, S., Baig, S. (2021). Detection of kras gene in colorectal cancer patients through liquid biopsy: A cost-effective method. Journal of College of Physicians and Surgeons Pakistan , 31(10), 1174–1178.
- 45.Cui, J., Placzek, W. J. (2016). PTBP1 modulation of MCL1 expression regulates cellular apoptosis induced by antitubulin chemotherapeutics. Cell Death & Differentiation , 23(10), 1681–1690.
- 46.Flannery, C. A., Rowzee, A. M., Choe, G. H., Saleh, F. L., Radford, C. C.et al. (2016). Development of a quantitative PCR assay for detection of human insulin-like growth factor receptor and insulin receptor isoforms. Endocrinology , 157(4), 1702–1708;
- 47.Bièche, I., Olivi, M., Noguès, C., Vidaud, M., Lidereau, R. (2002). Prognostic value of CCND1 gene status in sporadic breast tumours, as determined by real-time quantitative PCR assays. British Journal of Cancer , 86(4), 580–586.
- 48.Ioffe, Y. J., Chiappinelli, K. B., Mutch, D. G., Zighelboim, I., Goodfellow, P. J. (2012). Phosphatase and tensin homolog (PTEN) pseudogene expression in endometrial cancer: A conserved regulatory mechanism important in tumorigenesis? Gynecologic Oncology , 124(2), 340–346;
- 49.Kim, H. Y. (2017). Statistical notes for clinical researchers: Chi-squared test and Fisher’s exact test. Restorative Dentistry & Endodontics , 42(2), 152–155.
- 50.Livak, K. J., Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods , 25(4), 402–408;
- 51.Bland, J. M., Altman, D. G. (2004). The logrank test. British Medical Journal , 328(7447), 1073;
- 52.Dasari, S., Tchounwou, P. B. (2014). Cisplatin in cancer therapy: Molecular mechanisms of action. European Journal of Pharmacology , 740(Suppl. A), 364–378;
- 53.Lugones, Y., Loren, P. (2022). Cisplatin resistance: Genetic and epigenetic factors involved. Biomolecules , 12(10), 12101365.
- 54.Marsh, J. A., Teichmann, S. A. (2014). Parallel dynamics and evolution: Protein conformational fluctuations and assembly reflect evolutionary changes in sequence and structure. BioEssays , 36(2), 209–218;
- 55.Mordvinov, V., Pakharukova, M. (2022). Xenobiotic-metabolizing enzymes in trematodes. Biomedicines , 10(12), 10123039.
- 56.Zhong, Q., Xiao, X., Qiu, Y., Xu, Z., Chen, C.et al. (2023). Protein posttranslational modifications in health and diseases: Functions, regulatory mechanisms, and therapeutic implications. MedComm , 4(3), e261;
- 57.Reed, J. C. (1994). Bcl-2 and the regulation of programmed cell death. Journal of Cell Biology , 124(1–2), 1–6;
- 58.Zhao, T., He, Q., Xie, S., Zhan, H., Jiang, C.et al. (2023). A novel Mcl-1 inhibitor synergizes with venetoclax to induce apoptosis in cancer cells. Molecular Medicine , 29(1), 00565-7.
- 59.Perciavalle, R. M., Opferman, J. T. (2013). Delving deeper: MCL-1’s contributions to normal and cancer biology. Trends in Cell Biology , 23(1), 22–29;
- 60.Wang, Z., Fukushima, H., Gao, D., Inuzuka, H., Wan, L.et al. (2011). The two faces of FBW7 in cancer drug resistance. BioEssays , 33(11), 851–859;
- 61.Brotin, E., Meryet-Figuière, M., Simonin, K., Duval, R. E., Villedieu, M.et al. (2010). Bcl-XL and MCL-1 constitute pertinent targets in ovarian carcinoma and their concomitant inhibition is sufficient to induce apoptosis. International Journal of Cancer , 126(4), 885–895;
- 62.Kang, M. H., Reynolds, C. P. (2009). Bcl-2 inhibitors: Targeting mitochondrial apoptotic pathways in cancer therapy. Clinical Cancer Research , 15(4), 1126–1132;
- 63.Knight, T., Luedtke, D., Edwards, H., Taub, J. W., Ge, Y. (2019). A delicate balance-The BCL-2 family and its role in apoptosis, oncogenesis, and cancer therapeutics. Biochemical Pharmacology , 162, 250–261;
- 64.Czabotar, P. E., Lessene, G., Strasser, A., Adams, J. M. (2014). Control of apoptosis by the BCL-2 protein family: Implications for physiology and therapy. Nature Reviews Molecular Cell Biology , 15(1), 49–63;
- 65.Li, R., Pourpak, A., Morris, S. W. (2009). Inhibition of the insulin-like growth factor-1 receptor (IGF1R) tyrosine kinase as a novel cancer therapy approach. Journal of Medicinal Chemistry , 52(16), 4981–5004;
- 66.Gligorijević, N., Dobrijević, Z., Šunderić, M., Robajac, D., Četić, D.et al. (2022). The insulin-like growth factor system and colorectal cancer. Life , 12(8), 1274;
- 67.Amirkhah, R., Schmitz, U., Linnebacher, M., Wolkenhauer, O., Farazmand, A. (2015). MicroRNA-mRNA interactions in colorectal cancer and their role in tumor progression. Genes, Chromosomes and Cancer , 54(3), 129–141;
- 68.Singh, P., Alex, J. M., Bast, F. (2014). Insulin receptor (IR) and insulin-like growth factor receptor 1 (IGF-1R) signaling systems: Novel treatment strategies for cancer. Medical Oncology , 31(1), 1–14.
- 69.van der Meel, R., Oliveira, S., Altintas, I., Heukers, R., Pieters, EH.et al. (2013). Inhibition of tumor growth by targeted anti-EGFR/IGF-1R nanobullets depends on efficient blocking of cell survival pathways. Molecular Pharmaceutics , 10(10), 3717–3727;
- 70.Riedemann, J., Macaulay, V. (2006). IGF1R signalling and its inhibition. Endocrine Related Cancer , 13(Supplement_1), S33;
- 71.Gombos, A., Metzger-Filho, O., Dal Lago, L., Awada-Hussein, A. (2012). Clinical development of insulin-like growth factor receptor–1 (IGF-1R) inhibitors: At the crossroad? Investigational New Drugs , 30(6), 2433–2442;
- 72.Ramos-Garcia, P., Gil-Montoya, J., Scully, C., Ayén, A., González-Ruiz, L.et al. (2017). An update on the implications of cyclin D1 in oral carcinogenesis. Oral Diseases , 23(7), 897–912;
- 73.Wang, J., Su, W., Zhang, T., Zhang, S., Lei, H.et al. (2023). Aberrant Cyclin D1 splicing in cancer: From molecular mechanism to therapeutic modulation. Cell Death & Disease , 14(4), 244.
- 74.Sewify, E. M., Afifi, O. A., Mosad, E., Zaki, A. H., El Gammal, S. A. (2014). Cyclin D1 amplification in multiple myeloma is associated with multidrug resistance expression. Clinical Lymphoma Myeloma and Leukemia , 14(3), 215–222.
- 75.Pan, G., Liu, Y., Shang, L., Zhou, F., Yang, S. (2021). EMT-associated microRNAs and their roles in cancer stemness and drug resistance. Cancer Communications , 41(3), 199–217;
- 76.Sun, H., Ding, C., Zhang, H., Gao, J. (2016). Let‐7 miRNAs sensitize breast cancer stem cells to radiation‐induced repression through inhibition of the cyclin D1/Akt1/Wnt1 signaling pathway. Molecular Medicine Reports , 14(4), 3285–3292;
- 77.Liu, A., Zhu, Y., Chen, W., Merlino, G., Yu, Y. (2022). PTEN dual lipid and protein-phosphatase function in tumor progression. Cancers , 14, 14153666.
- 78.Steelman, L., Abrams, S., Whelan, J., Bertrand, F., Ludwig, D.et al. (2008). Contributions of the Raf/MEK/ERK, PI3K/PTEN/Akt/mTOR and Jak/STAT pathways to leukemia. Leukemia , 22(4), 686–707;
- 79.Wang, F., Li, T., Zhang, B., Li, H., Wu, Q.et al. (2013). MicroRNA-19a/b regulates multidrug resistance in human gastric cancer cells by targeting PTEN. Biochemical and Biophysical Research Communications , 434(3), 688–694;
- 80.Housman, G., Byler, S., Heerboth, S., Lapinska, K., Longacre, M.et al. (2014). Drug resistance in cancer: An overview. Cancers , 6(3), 1769–1792;
- 81.Lim, S. Y., Menzies, A. M., Rizos, H. (2017). Mechanisms and strategies to overcome resistance to molecularly targeted therapy for melanoma. Cancer , 123(S11), 2118–2129;
- 82.Thai, S. F., Jones, C. P., Robinette, B. L., Ren, H., Vallant, B.et al. (2021). Effects of copper nanoparticles on mRNA and small RNA expression in human hepatocellular carcinoma (HepG2) cells. Journal of Nanoscience and Nanotechnology , 21(10), 5083–5098;
- 83.Xu, J., Chang, W. H., Fong, L. W. R., Weiss, R. H., Yu, S. L.et al. (2019). Targeting the insulin-like growth factor-1 receptor in MTAP-deficient renal cell carcinoma. Signal Transduction and Targeted Therapy , 4(1), 2;
- 84.Bazzichetto, C., Conciatori, F., Pallocca, M., Falcone, I., Fanciulli, M.et al. (2019). PTEN as a prognostic/predictive biomarker in cancer: An unfulfilled promise? Cancers , 11(4), 435;
- 85.Berghuis, D., Santos, S. J., Baelde, H. J., Taminiau, A. H., Maarten Egeler, R.et al. (2011). Pro-inflammatory chemokine-chemokine receptor interactions within the Ewing sarcoma microenvironment determine CD8+ T-lymphocyte infiltration and affect tumour progression. The Journal of Pathology , 223(3), 347–357;
- 86.Maimela, N. R., Liu, S., Zhang, Y. (2019). Fates of CD8+ T cells in tumor microenvironment. Computational and Structural Biotechnology Journal , 17 1–13;
- 87.Sur, D., Balacescu, L., Cainap, S. S., Visan, S., Pop, L.et al. (2021). Predictive efficacy of miR-125b-5p, miR-17-5p, and miR-185-5p in liver metastasis and chemotherapy response among advanced stage colorectal cancer patients. Frontiers in Oncology , 11, 651380;
- 88.Cuzziol, C. I., Marzochi, L. L., Possebon, V. S., Kawasaki-Oyama, R. S., Mattos, M. F.et al. (2022). Regulation of VEGFA, KRAS, and NFE2L2 Oncogenes by microRNAs in head and neck cancer. International Journal of Molecular Sciences , 23(13), 7483;
- 89.Tian, J., Cheng, L., Kong, E., Gu, W., Jiang, Y.et al. (2022). modulates cisplatin resistance and angiogenesis through AKT1/GSK3β/VEGFA pathway in cervical cancer. Reproductive Biology and Endocrinology , 20, 1–12. 10.1186/s12958-022-00995-2;
- 90.Rashtchizadeh, N., Argani, H., Ghorbanihaghjo, A., Sanajou, D., Hosseini, V.et al. (2019). AMPK: A promising molecular target for combating cisplatin toxicities. Biochemical Pharmacology , 163(1), 94–100;
- 91.Fernández de Mattos S., Villalonga P., Clardy J., Lam E. W. (2008). FOXO3a mediates the cytotoxic effects of cisplatin in colon cancer cells. Molecular Cancer Therapeutics , 7(10), 3237–3246;
Associated Data
Data Availability Statement
No datasets were generated or analyzed during the current study.