- Research
- Open access
- Published:
ESCO2’s oncogenic role in human tumors: a pan-cancer analysis and experimental validation
BMC Cancer volume 24, Article number: 452 (2024)
Abstract
Purpose
Establishment of sister chromatid cohesion N-acetyltransferase 2 (ESCO2) is involved in the mitotic S-phase adhesins acetylation and is responsible for bridging two sister chromatids. However, present ESCO2 cancer research is limited to a few cancers. No systematic pan-cancer analysis has been conducted to investigate its role in diagnosis, prognosis, and effector function.
Methods
We thoroughly examined the ESCO2 carcinogenesis in pan-cancer by combining public databases such as The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression Project (GTEx), UALCAN and Tumor Immune Single-cell Hub (TISCH). The analysis includes differential expression analysis, survival analysis, cellular effector function, gene mutation, single cell analysis, and tumor immune cell infiltration. Furthermore, we confirmed ESCO2’s impacts on clear cell renal cell carcinoma (ccRCC) cells’ proliferative and invasive capacities in vitro.
Results
In our study, 30 of 33 cancer types exhibited considerably greater levels of ESCO2 expression in tumor tissue using TCGA and GTEx databases, whereas acute myeloid leukemia (LAML) exhibited significantly lower levels. Kaplan-Meier survival analyses in adrenocortical carcinoma (ACC), kidney chromophobe (KICH), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), brain lower grade glioma (LGG), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), mesothelioma (MESO), and pancreatic adenocarcinoma (PAAD) demonstrated that tumor patients with high ESCO2 expression have short survival periods. However, in thymoma (THYM), colon adenocarcinoma (COAD) and rectum adenocarcinoma (READ), ESCO2 was a favorable prognostic factor. Moreover, ESCO2 expression positively correlates with tumor stage and tumor size in several cancers, including LIHC, KIRC, KIRP and LUAD. Function analysis revealed that ESCO2 participates in mitosis, cell cycle, DNA damage repair, and other processes. CDK1 was identified as a downstream gene regulated by ESCO2. Furthermore, ESCO2 might also be implicated in immune cell infiltration. Finally, ESCO2’S knockdown significantly inhibited the A498 and T24 cells’ proliferation, invasion, and migration.
Conclusions
In conclusion, ESCO2 is a possible pan-cancer biomarker and oncogene that can reliably predict the prognosis of cancer patients. ESCO2 was also implicated in the cell cycle and proliferation regulation. In a nutshell, ESCO2 is a therapeutically viable and dependable target.
Introduction
Cancer in most nations gradually displaces cardiovascular disease as the leading cause of premature mortality [1]. According to statistics, there will be 19.3 million new cancer cases and almost 10 million cancer-related deaths in 2020. The global cancer burden is expected to reach 28.4 million patients by 2040 [2]. Cancer treatment has steadily transitioned into the era of precision therapy because of the advanced high-throughput sequencing technology. Advanced cancer patients frequently receive targeted therapy and immunotherapy in the clinic. For instance, palbociclib (a CDK4/6 inhibitor) in conjunction with letrozole (24.8 months) resulted in considerably longer progression-free survival than letrozole alone (14.5 months) in patients with previously untreated ER-positive, HER2-negative advanced breast cancer [3]. Due to the identification of molecular targets and the growing understanding of their cellular effects, small molecule inhibitors have been created as a primary therapeutic approach for cancer treatment. The discovery of new molecular targets may contribute to advancing combined therapies, eradicating drug resistance, and increasing therapeutic efficacy. Small molecule inhibitors have a bright future; thus, scientists should investigate novel biomarkers to find potential therapeutic targets.
The cell cycle is a tightly controlled process that supports genetic material replication and cell proliferation. One of the distinguishing characteristics of malignancies is abnormal cell proliferation brought on by cell cycle dysregulation [4]. Dysregulation of the cancer cell cycle facilitates cell proliferation, which is driven by excessive mitotic signaling, inhibitory checkpoint failure, or both. Scientists have demonstrated that it is effective to target proteins associated with the cell cycle in order to limit tumor growth. The most notable successes in targeting cell cycle mechanisms have been inhibitors of CDK4 and CDK6 [5,6,7]. The introduction of these CDK4/6 inhibitors for clinical use constitutes a milestone in the treatment of breast cancer and may have broad ramifications for the management of other tumor types [8]. Despite the current success of CDK4/6 inhibitors, cancer therapy that targets cell cycle proteins is still in its infancy. An in-depth exploration of cell cycle regulatory mechanisms and their role in cancer will guide current cancer treatments and identify new therapeutic targets.
Establishment of sister chromatid cohesion N-acetyltransferase 2 (ESCO2) was identified as an effective target for cancer therapy, which is a pivotal protein in the cell division process [9]. ESCO2, the human homolog of yeast ECO1, works on proliferating cell nuclear antigen (PCNA) to stimulate sister chromatid cohesion [10, 11]. Acetylation of the SMC3 subunit of the adhesive protein via ESCO2 acetyltransferase facilitates sister chromatid cohesion, inhibiting cohesin release from chromatin [12]. The cohesin is a multiprotein complex whose typical role is to bind sister chromatids from S-phase to anaphase to prevent premature segregation of sister chromatids and to ensure equal segregation of chromosomes [13, 14]. Hence, the presence of ESCO2 ensures correct chromosomal segregation and makes recombinant DNA repair possible. Recent research has demonstrated that ESCO2 is strongly linked to the formation of several cancer types [15,16,17,18,19]. ESCO2 knockdown in human gastric cancer cell lines in vitro significantly inhibited cell proliferation and induced apoptosis by regulating P53 [17]. Meanwhile, ESCO2 can promote LUAD cell proliferation and metastasis by promoting aerobic glycolysis [16]. Similarly, Fu et al. identified that 53BP1-MDS ring-like structure disruption caused by ESCO2 depletion in colorectal cancer cells, which reduced the effectiveness of non-homologous end joining (NHEJ) repair and made cancer cells more susceptible to chemotherapy [20]. However, Guo et al. discovered that ESCO2 overexpression in colon cancer reduced MMP2’s transcriptional activity to limit tumor metastasis [19]. ESCO2’s involvement in cancer is therefore controversial. Investigating ESCO2’s biological function in other malignancies is necessary. A pan-cancer analysis can investigate the oncogene in several malignancies, allowing oncogenes’ effector function and differential expression analyses.
We identified ESCO2’s abnormal expression value in tumor and normal samples by pan-cancer analysis and then validated its protein value in clinical tissues. Subsequently, we confirmed the value of ESCO2 in assessing prognosis and performed gene enrichment analysis, single cell analysis and immuno-infiltration analysis for ESCO2.Moreover, we conducted in vitro experiments to corroborate our bioinformatics findings. ESCO2 is a unique and useful biomarker and a prospective molecular target.
Materials and methods
Data preparation and differential expression analysis
The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/), which has 11,315 total samples, was employed to download the bulk mRNA sequencing profiles of 33 distinct cancer samples and their corresponding normal tissues. All gene expression data were standardized using the Transcripts Per Kilobase of Exon Model per Million Mapped Reads (TPM) for all samples. In the meanwhile, we collected the relevant clinical information on cancer patients. The Progression Free Interval (PFI), Disease Free Interval (DFI), and Disease-Specific Survival (DSS) timings were retrieved from UCSC Xena(https://xena.ucsc.edu/), and the Overall Survival (OS) time was acquired from TCGA.
Then, using the TIMER [21] and Xena Shiny databases (https://shiny.hiplot.cn/ucsc-xena-shiny/)-powerful online platforms that incorporated mRNA sequencing data from the TCGA and GTEx datasets [22]-we analyzed the ESCO2’s mRNA expression value between tumor and normal samples for 33 types of cancer. Furthermore, the ESCO2’s protein level was determined using the UALCAN online platform (https://ualcan.path.uab.edu/index.html). Finally, three pairs of paraffin-embedded kidney renal clear cell carcinoma (KIRC), bladder urothelial carcinoma (BLCA), and adjacent samples were collected from the Second Hospital of Tianjin Medical University. The pathologic type of these individuals was clear cell renal cell carcinoma (ccRCC, i.e., KIRC) or BLCA, and none of them received preoperative chemotherapy. More crucially, all patients have approved the use of the surgical material for academic research and publications. All methods were approved by The Institutional Review Committee and the Medical Ethics Committee of the Second Hospital of Tianjin Medical University. Written informed consents have been obtained from all subjects.
Survival analysis and relationship with clinical stage
We filtered the cancer patients to assess the ESCO2’s impacts on the cancer cases’ prognosis. (a) The average value of several samples from the same patient was thought to represent the ESCO2 expression level. (b) cancer patients with a follow-up time of zero were excluded. After that, the survival R package was applied to run the survival analysis. The survfit function was utilized to model the Kaplan-Meier survival curve. The best-cutoff point was discovered using the “surv cutpoint” function in a survminer R package, which can repeatedly check all viable cutting points to get the highest rank statistic. We next divided each cancer type’s patients into two groups using the appropriate cutoff point for gene expression. The Kaplan-Meier survival curves were compared based on the two-sided long-rank test. Moreover, we performed Receiver operating characteristic curve (ROC) testing to assess the ESCO2’s ability to predict cancer across all types. Then, we investigated the relationship between ESCO2’s expression level and clinical parameters, including tumor size, stage, and grade.
Functional enrichment analysis of ESCO2
We extracted the RNA sequencing profiles of cancer types where ESCO2 has prognostic value to investigate ESCO2’s carcinogenesis in malignancies. Hence, 12 cancer types, including thymoma (THYM), colon adenocarcinoma (COAD), rectum adenocarcinoma (READ), adrenocortical carcinoma (ACC), kidney chromophobe (KICH), KIRC, kidney renal papillary cell carcinoma (KIRP), brain lower grade glioma (LGG), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), mesothelioma (MESO), and pancreatic adenocarcinoma (PAAD), were rolled into the next analysis. For each cancer type, individuals were grouped into ESCO2-high and ESCO2-low subgroups relied on the median value of ESCO2, respectively. After that, gene set enriched in two groups were detected by gene set enrichment analysis (GSEA). Gene sets, the 50 hall markers, were obtained from MSigDB [23] (https://www.gsea-msigdb.org/gsea/msigdb).
The hunt for ESCO2-related regulatory genes
The above data shows that ESCO2 is a negative oncogene in most cancers. Therefore, we extracted the RNA-seq matrix from KICH, KIRC, LGG, LIHC, LUAD, and PAAD to study ESCO2-related regulatory genes. ESCO2 had a negative impact on these cancer types, which also accounted for more than 150 cases. Similarly, we split patients into ESCO2-high and ESCO2-low groups based on the median value of ESCO2 in each cancer type. Following that, differentially expressed genes (DEGs) between ESCO2-high and the ESCO2-low group were performed using limma R packge [24](p < 0.05, Log FC > 1). After identifying the intersecting genes of DEGs in each cancer, we calculated the association between ESCO2 and these DEGs by spearman’s correlation analysis for each cancer type, respectively (p < 0.05, Cor > 0.4). Intersecting genes for these highly correlated genes were considered ESCO2-related regulatory genes. Then, using the clusterProfiler R [25] package, the Kyoto Encyclopedia of Genes and Genomes (KEGG) [26,27,28] and Gene Ontology (GO) pathway enrichment analyses were performed on ESCO2-related regulatory genes. Incorporating ESCO2-related regulatory genes into the String database (https://cn.string-db.org/) to form a PPI network enables the identification of core genes.
Single-cell analysis of ESCO2 in cancers
Tumor Immune Single-cell Hub (TISCH, http://tisch.comp-genomics.org/home/) is a robust database created for multiple single-cell analyses containing nearly 190 single-cell databases. We estimated the ESCO2 expression level across multiple cancer types in each cell type. Our study focused on the distribution of ESCO2 in KIRC cell subpopulations. The expression levels of ESCO2 in each cell type were quantified and visualized by a heatmap, scatter diagrams, and violin plots. In addition to the single-cell dataset of KIRC contained in the TISCH database, we also downloaded additional single-cell sequencing data from GSE156632 [29] in the Gene Expression Omnibus (GEO) database to validate ESCO2 expression level in each cell types.
Cell culture and siRNA transfection
The Cell Resource Center Affiliated with the Chinese Academy of Medical Sciences provided the human ccRCC cell lines A498 and human BLCA cell lines T24. In a 5% CO2 humidified cell incubator at 37 °C, A498 and T24 cells were grown in Dulbecco’s modified Eagle medium (DMEM) (Gibco BRL Life Technologies Inc., USA) with 10% FBS (Gibco BRL Life Technologies Inc., USA) and 1% penicillin-streptomycin (Hy-clone, CA, USA).
Small interfering RNAs (siRNAs) were ordered from GenePharma. The si-NC and si-ESCO2 siRNA were transfected into A498 and T24 cells with a Transfection reagent. When cells occupied 70–80% of a 50 mm Petri dish, transfection experiments were conducted; the whole medium was aspirated, and 2 mL of opti-MEM was applied to each well. The siRNA (GenePharma, Shanghai, China) and Transfection reagent (Cat#:L3000075, Thermo Fisher Scientific, United States) were added to 500 µL opti-MEM, mixed, and incubated for 20 min at room temperature. After that, the siRNA-Transfection reagent mixture was added to each well. Each well was switched to a complete medium 6 h after the transfection assay. The sequences for ESCO2 siRNA are as follows: si-ESCO2-1 (forward: GCAAAUCAAGGCUCACCAUTT; reverse: AUGGUGAGCCUUGAUUUGCTT) and si-ESCO2-2 (forward: CUCUUAGACCAGGAUUAUCTT; reverse: GAUAAUCCUGGUCUAAGAGUG).
qPCR and western blotting analysis
After the cells were transfected for 48 and 72 h, we harvested cells for RNA and protein extraction, respectively. Following the protocol, the total RNAs of cells were extracted using TRIzol reagents (Invitrogen). We then used cDNA reverse transcription kits to synthesize cDNA, and the mRNA expression level was calculated through qPCR assay in q225. We used the 2−ΔΔCt method to calculate GAPDH as an internal reference for mRNA. The ESCO2’s sequences of primers are as follows: F, CACTGGGACGCACCCAAAA, R, CACTTGCCTTGTCGCAAAAG (Sangon biotech). The following are the primer sequences for the GAPDH gene.: F, CAGGAGGCATTGCTGATGAT, R, GAAGGCTGGGGCTCATTT (Sangon biotech).
The tissues or cells were collected and lysed in RIPA lysis buffer (Cat#:89,901, Thermo Fisher Scientific, USA) with supplementation of a protease inhibitor cocktail (Cat#:539,131, Roche, Switzerland) in the ice for 30 min. The total protein extracts were measured by the ABC method (Cat#:23,227, Thermo Fisher Scientific, USA), separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and transferred to the nitrocellulose membrane. The membrane was blocked using 5% non-fat milk or 5% BSA (Solarbio, Beijing, China), then incubated with primary antibodies anti- (ESCO2,Cat#:23525-1-AP, diluted 1:2000, Proteintech, USA; CDK1, Cat#:23525-1-AP, diluted 1:2000, Proteintech, USA). The proper secondary antibody is chosen based on the primary antibody’s type, including goat anti-mouse HRP-conjugated IgG (Cat#:SA00001-1, 1:2000 dilutions, Proteintech, USA) and goat anti-rabbit HRP-conjugated IgG (Cat#:SA00001-2, 1:2000 dilutions, Proteintech, USA). The target band signals were captured with the help of a BIO-RAD ChemiDoc XRS chemiluminescence system (Bio-Rad Inc., CA, USA). The densitometry of all Western blots was analysed using a Gel-Pro analyzer (Media Cybernetics, USA). Imagej software was used to perform quantitative analysis with β-actin as a control (1.52a, National Institutes of Health, Bethesda, MD, USA).
Cell proliferation and colony formation assay
Briefly, 2000 cells per well were planted onto 96-well plates for the cell proliferation assay, and transfection assay was conducted after cell fixation. Cells were cultured in the incubator for 12, 24, 48, and 72 h. After that, each well received 10 µL of Cell Counting Kit-8 reaction solution (CCK8, Cat#:CA1210, Solarbio Science & Technology Co., Ltd., Beijing, China) and was incubated for 1 h. The light absorption values (OD) for si-NC and si-ESCO2 groups were captured at 450 nm on a Varioskan Flash Multimode plate reader (Thermo Fisher Scientific Inc., USA).
Cancer cells per well were seeded onto a 12-well plate, and the medium was changed every two days. After being stained with 0.5% crystal violet (Cat#: C0121, Beyotime) for 20 min and fixed with 4% paraformaldehyde at room temperature for 10 min, colonies were counted 14 days later.
Ethynyl-2-Deoxyuridine (EdU) incorporation assay
Cancer cells were seeded in 3.5 cm confocal dishes (Thermo Fisher Scientific Inc., USA) with 50,000 cells per dish. An appropriate amount of 50 µM EdU (Cat#:K1077, ApexBio Technology LLC) medium was prepared by diluting the EdU solution with cell culture medium at a ratio of 1000:1.1 ml of 50 µM EdU medium was added after the entire medium had been aspirated, and it was then incubated for 2 h. After discarding the medium, it was washed with PBS, fixed with 4% paraformaldehyde, and then permeabilized with 0.5% TritonX-100 (Cat#:X100, Sigma-Aldrich, USA). 200 μm of 1x Apollo (Cat#:K1077, ApexBio Technology LLC) staining reaction solution was added to each dish, incubated in the dark for 30 min. Nuclei were stained with DAPI (Cat#:MBD0015, Sigma-Aldrich, USA) for 10 min. Using a Nikon ECLIPSE 80i fluorescence microscope, images were captured, and ImageJ was used to count the cell number (1.52a, National Institutes of Health, Bethesda, MD, USA).
Transwell migration and wound healing assay
Matrigel (BD Biosciences, Franklin Lakes, NJ, USA) was evenly spread in transwell chambers (Cat#: 3422, Corning Inc., Corning, NY, USA), and cells starved for 24 h were seeded in a serum-free medium in the upper chamber. Medium containing 10% FBS was added to the lower chamber as a catalyst. The upper chamber was removed and separated after roughly 24 h, then fixed with methanol for 10 min, rinsed with PBS, and stained with 0.1% crystal violet (Cat#:G1064, Solarbio Science & Technology Co., Ltd., Beijing, China). Under a microscope, cells that had been pierced were viewed and counted (CKX41, Olympus, Tokyo, Japan).
Cells were seeded uniformly in micro-insert 4-wells in a µ-dish35 mm, high ibiTreat (ibidi GmbH, Germany). The culture inserts were removed when the cell density reached 90%. The 500 mm wide scratching gaps were washed with phosphate-buffered solution (PBS, Cat#:D1020, Solarbio Science & Technology Co., Ltd., Beijing, China) and were replaced with fresh medium. The dynamic changes of wells were recorded at 0, 12, 24 and 48 h on a microscope (Olympus, Japan).
Statistical analysis
Using the Wilcox test, comparisons between the two groups were made. For analyzing three or more groups, the one-way ANOVA test was performed. For all statistical calculations, R studio and GraphPad software were used. The threshold for significance was defined at P < 0.05.
Results
ESCO2 expression level in various cancer
Typically, tumor tissue exhibits higher or lower oncogene expression levels than normal tissue. Using the TCGA database, our work found that ESCO2 is substantially expressed in the tumor tissue of most cancer types. However, the ESCO2 expression value did not significantly change in KICH, PAAD, pheochromocytoma and paraganglioma (PCPG), or READ (Fig. 1A). Considering only a few normal samples available in the TCGA database, we combined the TCGA and GTEx databases. 30 of 33 cancer types exhibited considerably greater levels of ESCO2 expression in tumor tissue, whereas acute myeloid leukemia (LAML) exhibited significantly lower levels. In the ACC and PCPG, the ESCO2 expression levels of the tumor and normal tissues were similar (Fig. 1B). ESCO2 expression was raised in most cancers, supporting its oncogene function.
Furthermore, we evaluated the ESCO2 protein levels. We only found significantly higher protein levels of ESCO2 in tumor tissues of uterine corpus endometrial carcinoma (UCEC), breast invasive carcinoma (BRCA), and LUAD due to the UALCAN database’s limited inclusion of cancer types (Fig. 1C). To further explore the alterations in ESCO2’s protein levels in carcinogenesis, we collected samples of KIRC and BLCA. In KIRC and BLCA, the protein level of tumor samples was noticeably greater than para-cancerous ones (Fig. 1D).
ESCO2’s prognostic value across cancers
Using the TCGA database and focusing on OS, DFI, DSS, and PFI, we conducted a Kaplan-Meier survival analysis in each cancer to examine the association between ESCO2 and the cancer patient’s prognosis. Kaplan-Meier survival analyses in ACC, KICH, KIRC, KIRP, LGG, LIHC, LUAD, MESO, and PAAD demonstrated that tumor patients with high ESCO2 expression have short survival periods (Fig. 2D, E, F, G, H, I, J, K and L). While ESCO2’s expression was linked to better outcomes in patients with THYM, COAD, and READ (Fig. 2A, B and C). ROC analysis revealed that the ESCO2 performed well in predictive accuracy, with AUC values in ACC, KICH, KIRP, LGG, MESO, LIHC, and LUAD performing well (Figure S1A). High ESCO2 expression was linked to a worse prognosis in cancer patients with KIRP, LIHC, and PAAD, according to Kaplan-Meier survival analysis for DFI (Figure S1B). Similarly, as shown in Figures S2 and S3, cancer patients with increased ESCO2 expression had worse DSS and PFI across various malignancies. Our investigation revealed that ESCO2 might enhance tumor initiation and progression.
ESCO2 was positively related to tumor progression
To confirm the ESCO2 role in cancer progression, a correlation analysis was performed between the ESCO2’s mRNA value and clinical features. We discovered that ESCO2 expression was significantly increased in advanced tumor stage in several malignancies, including KICH, KIRC, KIRP, LUAD, and LIHC. On the contrary, ESCO2 expression was higher in lower tumor stages than in higher tumor stages in COAD (Fig. 3A).
In LIHC, KIRC, LGG, and PAAD, ESCO2 expression was positively correlated with tumor grade (Fig. 3C). Similarly, ESCO2 may also affect the tumor size: Tumor size and ESCO2 expression were positively linked in LIHC, KIRC, KIRP, and LUAD. (Figs. 3D). Furthermore, ESCO2 was crucial for the metastasis of COAD, KIRC, and LUAD (Fig. 3B). Tumor stage, size, and metastasis are the most widely used clinical indicators for evaluating the progression of cancer. Our investigation verified the connection between ESCO2 and the aforementioned indicators, indicating that ESCO2 is involved in tumor progression.
Analysis of ESCO2-related regulatory pathways
We gathered RNA sequencing profiles of cancer types where ESCO2 could significantly affect patients’ overall survival times to examine the ESCO2 oncogenic role in malignancies. For each cancer type described in the methods, cancer cases were grouped into ESCO2-high and -low subgroups according to the median value of ESCO2, respectively. Gene set enriched in two groups were then detected by gene set enrichment analysis. For visualization, we chose the top-10 NSE-ranked enriched pathways. Intriguingly, GSEA results were remarkably consistent across the 12 cancer types, indicating the validity of our findings (Fig. 4A, B, C, D, E, F, G, H, I, J, K and L). Genes upregulated in the ESCO2-high group showed the E2F targets enrichment, G2M checkpoint, mitotic spindle and Mtorc1 signaling, which were associated with cell cycle and proliferation. Meanwhile, these ESCO2-high patients exhibited highly enriched DNA repair. In addition to cell proliferation-related pathways, glycolysis and fatty acid metabolism pathways were enriched in the ESCO2-high group in several cancer types. Furthermore, ESCO2 seemed to participate in the modulation of cancer inflammation, such as IL2 STAT5 signaling.
To investigate ESCO2-related regulatory genes, we divided cancer patients into ESCO2-high and ESCO2-low groups based on the median ESCO2 value in six cancer types specified in the materials and methods. We then obtained 2419 intersecting genes, which mainly enriched on histone binding, ATP hydrolysis activity and helicase activity (Fig. 5A and B). The results of the GO analysis also confirmed that ESCO2 was primarily responsible for controlling cell growth. After spearman’s correlation analysis, we further identified the ESCO2-related regulatory genes. 249 ESCO2-related regulatory genes were collected in our study, which was highly correlated with ESCO2 across cancer types (Fig. 5C). These ESCO2-related regulatory genes were primarily involved in cell-division-related processes, like the mitotic cell cycle and DNA metabolic processes (Fig. 5F). We successfully identified the core genes by building the PPI network of ESCO2-related regulatory genes. The regulatory genes network’s nucleus is CDK1 (Fig. 5D). Cell cycle kinase family member Cyclin-dependent kinase 1 (CDK1) significantly impacts cell cycle progression [30, 31]. After the ESCO2 knockdown, we discovered that the CDK1 expression was significantly decreased (Fig. 5E). Therefore, CDK1 may be an important signaling pathway regulated by ESCO2 in KIRC.
ESCO2 analysis at the single-cell level
Recent studies have shown that cell cycle proteins play an additional role in tumor development, affecting not only tumor cells but also their microenvironment and modulating anti-tumor immune responses [32]. For instance, inhibition of CDK4/6 enhances tumor cell immunogenicity through multiple mechanisms [7]. Meanwhile, CDK4/6 inhibitors substantially reduce the proliferation of regulatory T cells and encourage cytotoxic T cells to eliminate tumor cells [33]. Therefore, we next explored the ESCO2’s role in the anti-tumor immune response. By scRNA-seq, we can investigate the ESCO2 at the single-cell level. Using single-cell sequencing of COAD, KIRC, KICH, LIHC and PAAD, we observed that ESCO2 was expressed not only in malignant cells but also in endothelial and immune cells, including T cells, macrophages, NK cells and B cells (Fig. 6A). Significantly, ESCO2 was strongly expressed in proliferating T cells. We found that ESCO2 was almost exclusively expressed in proliferating T cells in KIRC (Fig. 6B). To further validate the results observed in the TICSH database, we subjected all quality single cells of KIRC(GSE156632) to single-cell processing procedures. Similarly, ESCO2 was significantly expressed in proliferating T cells, but we observed a higher expression of ESCO2 in tumor cells and macrophages as well (Fig. 6C, D and E).
ESCO2 Promotes Human ccRCC and BLCA cell proliferation and invasion in vitro
To confirm ESCO2’s biological effector role in carcinogenesis, we first knocked it down in human ccRCC cells. qt-PCR and Western blotting results showed that ESCO2 was successfully knocked down (Fig. 7A). Results from the CCK8 assay, colony formation assay, and EdU staining showed that the A498 cells’ ability to proliferate was significantly decreased by ESCO2 suppression than the control group (Fig. 7B, C and F). This is consistent with our previous bioinformatic analysis that ESCO2 significantly regulates cell cycle progression. To further identify the potential role of ESCO2 in the invasive and migration of A498 cells. According to the findings of Transwell and Wound Healing, silencing ESCO2 inhibits invasion and migration (Fig. 7D and E). To confirm ESCO2 as a potential pan-cancer biomarker, we further validated the carcinogenesis of ESCO2 in bladder cancer cell lines. According to the results of the experiments in vitro, the proliferative capacity of T24 cells was significantly inhibited after ESCO2’s knockdown (Fig. 8A, B, C and F). Furthermore, the invasive ability of the T24 cells after knockdown of ESCO2 significantly decreased compared with the control group (Fig. 8D and E). We confirmed, therefore, that ESCO2 is required for the proliferation and invasion of human ccRCC and BLCA cells.
Discussion
Cancer treatment has evolved significantly, with the introduction of immunotherapy and targeted medicines increasing patient survival rates for those with advanced or metastatic cancer, but overall clinical outcomes remain disappointing. The challenge of treating advanced cancers has inspired researchers to explore the underlying mechanisms that contribute to cancer growth, which will help identify potentially effective therapeutic targets. Unquestionably, one of the most promising therapeutic targets is cell cycle related protein. Palbociclib, the first CDK4/6-specific inhibitor, was created in 2004 and has proven effective against a variety of human cancer cell lines [34]. All three CDK4/6 inhibitors (palbociclib, ribociclib, and abemaciclib) are currently approved by the United States Food and Drug Administration (FDA) for the treatment of breast cancer. In addition, numerous inhibitors of cell-cycle proteins, including CDK9, CDK2, and CDK5, are in clinical trials [35, 36]. Cell cycle regulating proteins such as CDCA4, CDCA8, and KIF2C have been linked to tumor growth and progression, affecting the proliferation, migration, invasion and metastasis of cancer cell lines [37,38,39]. ESCO2 was first reported in Roberts syndrome, whose inactivation mutation causes Roberts Syndrome [11]. Subsequent researchers have identified that ESCO2 is crucial in controlling cell mitosis and preserving genomic stability because it bridges two sister chromatids and is involved in the mitotic S-phase adhesins acetylation [40, 41]. As described in the context, current researchers consider ESCO2 an excellent therapeutic target as an oncogene that promotes cancer development. However, current cancer research on ESCO2 is limited to a few cancers, such as lung, stomach, and colon cancers. No studies are focusing on ESCO2 in multiple cancers that can shed light on the similarities of ESCO2 in cancer. In this study, we identified that ESCO2 was a reliable biomarker for cancer patients and could accurately predict the cancer patient’s prognosis. Furthermore, ESCO2 participated in mitosis, the cell cycle, DNA damage repair, other processes, and tumor immune infiltration. Finally, we confirmed that ESCO2 is essential for the proliferation and invasion of human ccRCC and BLCA cells in vitro.
First, we evaluated the ESCO2’s mRNA expression value in 33 cancers by TCGA and GTEx databases. The analysis showed that ESCO2 expression was upregulated in most cancer tissues except LAML, ACC and PCPG. Due to post-transcriptional processing, protein and mRNA gene expression levels may differ. Thus, we validated ESCO2 protein levels in KIRC and BLCA clinical samples. ESCO2 protein values were considerably greater in malignant tissues than in para-cancerous ones, consistent with the mRNA analysis. These results confirmed that ESCO2 was upregulated in various cancers, suggesting a promising future for ESCO2 in cancer diagnosis. Additionally, OS, DSS, DFI, and PFI analyses all revealed that ESCO2 was significantly associated with cancer patient’s prognosis. Evaluated ESCO2 resulted in poorer prognosis in ACC, BLCA, BRCA, CESC, COAD, GBM, HNSC, KIRC, LGG, LIHC, LUAD, LUSC, and PAAD. However, in THYM, COAD and READ, ESCO2 was a favorable prognostic factor. ESCO2 has been shown to prevent cancer metastasis in a recent trial on colon cancer, which is consistent with the results of our bioinformatics research [19]. According to the findings above, ESCO2 is crucial in determining a cancer patient’s prognosis and can serve as a reliable prognostic biomarker. To some extent, tumor grading and clinical staging can show how tumors progress. Our results show that ESCO2 expression was positively associated with tumor stage and size in LIHC, KIRC, KIRP and LUAD. Furthermore, ESCO2 was crucial for the metastasis of COAD, KIRC, and LUAD. Our research identified that ESCO2 is a proto-oncogene that is linked to the development and progression of tumors.
Emerging research has confirmed that ESCO2 participated in apoptosis and cell proliferation-related pathways, such as the P53 and mTOR pathway [15, 17]. Furthermore, ESCO2 could stimulate aerobic glycolysis in LUAD cells by upregulating PKM2 and downregulating PKM1 expression [16]. The GSEA results in our study were unexpectedly uniform among the 12 different cancers, indicating the validity of our findings. The cell cycle and proliferation-related E2F targets, G2M checkpoint, mitotic spindle, and Mtorc1 signaling were concentrated in the genes upregulated in the ESCO2-high group. These ESCO2-high patients also displayed greatly enhanced DNA repair. In sum up, our bioinformatics analysis confirmed ESCO2’s contribution to the control of cell proliferation and cycle. Therefore, we conducted experimental validation to verify the results of the bioinformatics analysis. We knocked down ESCO2 in A498 and T24 cells, and ESCO2’S knockdown significantly inhibited the proliferation, invasion, and migration of A498 and T24 cells. This result demonstrated that ESCO2’s role in cancer might be unified, further indicating the possibility that it could be an effective therapeutic target. Furthermore, we investigated its potential regulatory mechanisms. We found the core of the ESCO2-CDK1 regulatory network. When ESCO2 was inhibited, CDK1 expression decreased, demonstrating that CDK1 is a downstream gene controlled by ESCO2. Cell cycle progression is the central event in all proliferating cells and is primarily regulated by cell cycle-dependent kinases (CDK) [42].
CDK1 is the only CDK required for the G2-M transition and regulates G1 progression and G1-S transition [43]. In recent years, CDK1 has been suggested as a therapeutic target for cancer. CDK1 overexpression has been found in many cancers, including gastric cancer, ovarian cancer, oral squamous cell carcinoma, liver cancer, and breast cancer [44]. CDK1 inhibitors could be a potential small-molecule drug. CDK1 inhibitor RO3306 could improve the efficacy of sorafenib treatment by targeting cancer stem cells in a preclinical model of hepatocellular carcinoma [45]. In addition, Several CDK1 inhibitors, including Rigosertib (phase II/III) and Zotaraciclib, have begun phase I clinical trials for treating pancreatic cancer and glioma [46, 47]. As a CDK1’s upstream regulatory gene, ESCO2 took part in the control of the cell cycle. Hence, ESCO2 is likely to be a promising therapeutic target in the future. Finally, single-cell analysis observed that ESCO2 was expressed not only in malignant cells but also in endothelial and immune cells, especially proliferating T cells. The GSEA results also indicate that ESCO2 is involved in cancer inflammation control, including IL2 STAT5 signaling, which raises the possibility that ESCO2 also controls the tumor immune microenvironment. Studies have shown that inhibition of cell cycle protein not only induce tumor cell cycle arrest, but also promote anti-tumor immunity. Inhibition of CDK4 /6, for example, increases levels of PD-L1 protein. More excitingly, combining CDK4/6 inhibitor with anti-PD-1 immunotherapy promotes tumor regression and significantly increases overall survival in mouse tumor models [48]. Therefore, research into ESCO2’s function in the tumor immune microenvironment is essential.
This study has the following limitations, though. First, this work does not apply any unique clinical cohorts to evaluate the diagnostic and prognostic significance of ESCO2 in cancer, despite its extensive use of sequencing data from public databases. Second, this study only performed ESCO2’s functional experiment in ccRCC cells; the upstream and downstream pathways of ESCO2 were not comprehensively examined, and the precise molecular mechanism of ESCO2 regulation is still unknown. Third, we did not conduct experiments in vivo, which is an important issue. We should attempt to resolve this issue in the future.
In conclusion, ESCO2 is a potential biomarker and oncogene for pan-cancer that can accurately predict the cancer patient’s prognosis in ACC, KICH, KIRC, KIRP, LGG, LIHC, LUAD, MESO, PAAD, THYM, COAD, and READ. Furthermore, our bioinformatics results discovered that ESCO2 is involved in cell division and cell cycle regulation and verified that ESCO2 is essential for the proliferation and invasion of human ccRCC and BLCA cells in vitro. In a nutshell, ESCO2 is a potential and reliable therapeutic target.
Data availability
The data analyzed in this study can be downloaded from the GEO (GSE156632, https://www.ncbi.nlm.nih.gov/geo/) and TCGA (https://portal.gdc.cancer.gov/).
Abbreviations
- DEG:
-
differentially expressed genes
- ORR:
-
objective response rate
- QC:
-
quality control
- tSNE:
-
t-Distributed Stochastic Neighbor Embedding
- ROC:
-
Receiver operating characteristic curve
- GSEA:
-
gene set enrichment analysis
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- TME:
-
tumor microenvironment
References
Bray F, Laversanne M, Weiderpass E, Soerjomataram I. The ever-increasing importance of cancer as a leading cause of premature death worldwide. Cancer. 2021;127(16):3029–30.
Sung H, Ferlay J, Siegel RL, et al. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and Mortality Worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.
Finn RS, Martin M, Rugo HS, et al. Palbociclib and Letrozole in Advanced breast Cancer. N Engl J Med. 2016;375(20):1925–36.
Matthews HK, Bertoli C, de Bruin R. Cell cycle control in cancer. Nat Rev Mol Cell Biol. 2022;23(1):74–88.
Weiss JM, Csoszi T, Maglakelidze M, et al. Myelopreservation with the CDK4/6 inhibitor trilaciclib in patients with small-cell lung cancer receiving first-line chemotherapy: a phase Ib/randomized phase II trial. Ann Oncol. 2019;30(10):1613–21.
Suski JM, Braun M, Strmiska V, Sicinski P. Targeting cell-cycle machinery in cancer. Cancer Cell. 2021;39(6):759–78.
Goel S, Bergholz JS, Zhao JJ. Targeting CDK4 and CDK6 in cancer. Nat Rev Cancer. 2022;22(6):356–72.
Wang R, Xu K, Gao F, Huang J, Guan X. Clinical considerations of CDK4/6 inhibitors in triple-negative breast cancer. Biochim Biophys Acta Rev Cancer. 2021;1876(2):188590.
Liu Z, Cheng X, Pang B, et al. Effects of ESCO2 or its methylation on the prognosis, clinical characteristics, immune microenvironment, and pathogenesis of low-grade glioma. Int Immunopharmacol. 2022;104:108399.
Bender D, Da Silva E, Chen J, et al. Multivalent interaction of ESCO2 with the replication machinery is required for sister chromatid cohesion in vertebrates. Proc Natl Acad Sci U S A. 2020;117(2):1081–9.
Vega H, Waisfisz Q, Gordillo M, et al. Roberts syndrome is caused by mutations in ESCO2, a human homolog of yeast ECO1 that is essential for the establishment of sister chromatid cohesion. Nat Genet. 2005;37(5):468–70.
Jevitt AM, Rankin BD, Chen J, Rankin S. The cohesin modifier ESCO2 is stable during DNA replication. Chromosome Res. 2023;31(1):6.
Michaelis C, Ciosk R, Nasmyth K. Cohesins: chromosomal proteins that prevent premature separation of sister chromatids. Cell. 1997;91(1):35–45.
Guacci V, Koshland D, Strunnikov A. A direct link between sister chromatid cohesion and chromosome condensation revealed through the analysis of MCD1 in S. Cerevisiae. Cell. 1997;91(1):47–57.
Wang QL, Liu L. Establishment of cohesion 1 homolog 2 facilitates cell aggressive behaviors and induces poor prognosis in renal cell carcinoma. J Clin Lab Anal. 2020;34(5):e23163.
Zhu HE, Li T, Shi S, Chen DX, Chen W, Chen H. ESCO2 promotes lung adenocarcinoma progression by regulating hnRNPA1 acetylation. J Exp Clin Cancer Res. 2021;40(1):64.
Chen H, Zhang L, He W, et al. ESCO2 knockdown inhibits cell proliferation and induces apoptosis in human gastric cancer cells. Biochem Biophys Res Commun. 2018;496(2):475–81.
Zhang T, Wang B, Gu B et al. Genetic and Molecular Characterization Revealed the Prognosis Efficiency of Histone Acetylation in Pan-Digestive Cancers. J Oncol. 2022. 2022: 3938652.
Guo XB, Huang B, Pan YH, Su SG, Li Y. ESCO2 inhibits tumor metastasis via transcriptionally repressing MMP2 in colorectal cancer. Cancer Manag Res. 2018;10:6157–66.
Fu J, Zhou S, Xu H, et al. ATM-ESCO2-SMC3 axis promotes 53BP1 recruitment in response to DNA damage and safeguards genome integrity by stabilizing cohesin complex. Nucleic Acids Res. 2023;51(14):7376–91.
Li T, Fu J, Zeng Z, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509–14.
The GTEx. Consortium atlas of genetic regulatory effects across human tissues. Science. 2020;369(6509):1318–30.
Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1(6):417–25.
Ritchie ME, Phipson B, Wu D, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.
Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–92.
Zhang M, Zhai W, Miao J, et al. Single cell analysis reveals intra-tumour heterogeneity, microenvironment and potential diagnosis markers for clear cell renal cell carcinoma. Clin Transl Med. 2022;12(5):e713.
Izadi S, Nikkhoo A, Hojjat-Farsangi M, et al. CDK1 in breast Cancer: implications for theranostic potential. Anticancer Agents Med Chem. 2020;20(7):758–67.
Ren L, Yang Y, Li W, et al. CDK1 serves as a therapeutic target of adrenocortical carcinoma via regulating epithelial-mesenchymal transition, G2/M phase transition, and PANoptosis. J Transl Med. 2022;20(1):444.
Goel S, DeCristo MJ, Watt AC, et al. CDK4/6 inhibition triggers anti-tumour immunity. Nature. 2017;548(7668):471–5.
Deng J, Wang ES, Jenkins RW, et al. CDK4/6 inhibition augments Antitumor Immunity by enhancing T-cell activation. Cancer Discov. 2018;8(2):216–33.
Fry DW, Harvey PJ, Keller PR, et al. Specific inhibition of cyclin-dependent kinase 4/6 by PD 0332991 and associated antitumor activity in human tumor xenografts. Mol Cancer Ther. 2004;3(11):1427–38.
Diamond JR, Boni V, Lim E, et al. First-in-human dose-escalation study of cyclin-dependent kinase 9 inhibitor VIP152 in patients with Advanced malignancies shows early signs of clinical efficacy. Clin Cancer Res. 2022;28(7):1285–93.
Kumar SK, LaPlant B, Chng WJ, et al. Dinaciclib, a novel CDK inhibitor, demonstrates encouraging single-agent activity in patients with relapsed multiple myeloma. Blood. 2015;125(3):443–8.
Fang H, Sheng S, Chen B, et al. A Pan-cancer analysis of the oncogenic role of Cell Division Cycle-Associated protein 4 (CDCA4) in human tumors. Front Immunol. 2022;13:826337.
Xiang C, Sun WH, Ke Y, Yu X, Wang Y. CDCA8 contributes to the development and progression of thyroid Cancer through regulating CDK1. J Cancer. 2022;13(7):2322–35.
Wei S, Dai M, Zhang C, et al. KIF2C: a novel link between Wnt/β-catenin and mTORC1 signaling in the pathogenesis of hepatocellular carcinoma. Protein Cell. 2021;12(10):788–809.
Alomer RM, da Silva E, Chen J, et al. Esco1 and Esco2 regulate distinct cohesin functions during cell cycle progression. Proc Natl Acad Sci U S A. 2017;114(37):9906–11.
Sun H, Zhang J, Xin S, et al. Cul4-Ddb1 ubiquitin ligases facilitate DNA replication-coupled sister chromatid cohesion through regulation of cohesin acetyltransferase Esco2. PLoS Genet. 2019;15(2):e1007685.
Evan GI, Vousden KH. Proliferation, cell cycle and apoptosis in cancer. Nature. 2001;411(6835):342–8.
Santamaría D, Barrière C, Cerqueira A, et al. Cdk1 is sufficient to drive the mammalian cell cycle. Nature. 2007;448(7155):811–5.
Sofi S, Mehraj U, Qayoom H, et al. Targeting cyclin-dependent kinase 1 (CDK1) in cancer: molecular docking and dynamic simulations of potential CDK1 inhibitors. Med Oncol. 2022;39(9):133.
Wu CX, Wang XQ, Chok SH, et al. Blocking CDK1/PDK1/β-Catenin signaling by CDK1 inhibitor RO3306 increased the efficacy of sorafenib treatment by targeting cancer stem cells in a preclinical model of hepatocellular carcinoma. Theranostics. 2018;8(14):3737–50.
Wu J, Yuan Y, Long Priel DA, et al. Phase I study of Zotiraciclib in Combination with Temozolomide for patients with recurrent high-grade Astrocytomas. Clin Cancer Res. 2021;27(12):3298–306.
O’Neil BH, Scott AJ, Ma WW, et al. A phase II/III randomized study to compare the efficacy and safety of rigosertib plus gemcitabine versus gemcitabine alone in patients with previously untreated metastatic pancreatic cancer. Ann Oncol. 2015;26(9):1923–9.
Zhang J, Bu X, Wang H, et al. Cyclin D-CDK4 kinase destabilizes PD-L1 via cullin 3-SPOP to control cancer immune surveillance. Nature. 2018;553(7686):91–5.
Acknowledgements
We thank Yubo Feng for assistance with the statistical methods.
Funding
This work was supported by the National Natural Science Foundation of China (grant numbers: 22176142 and 22076138), Young Elite Scientists Sponsorship Program by Tianjin (No. TJSQNTJ-2020-07).
Author information
Authors and Affiliations
Contributions
Conceptualization: YZ, ZZ, YB, and JZ. Methodology: YH, DC, ZZ, QF, and BY. Software: DC, ZZ, QF, and YJ. Validation: YH, DC, YB, and YJ. Writing-original draft: YH and DC. Writing-review and editing: YZ, ZZ, YB, and JZ. Funding acquisition: ZZ, JZ. Supervision: YZ and ZZ. All authors reviewed and revised the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The authors confirm that all methods were conducted according to the principles of the Declaration of Helsinki and were approved by The Institutional Review Committee and the Medical Ethics Committee of the Second Hospital of Tianjin Medical University. Written informed consents have been obtained from all subjects.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Huang, Y., Chen, D., Bai, Y. et al. ESCO2’s oncogenic role in human tumors: a pan-cancer analysis and experimental validation. BMC Cancer 24, 452 (2024). https://doi.org/10.1186/s12885-024-12213-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12885-024-12213-w