Comprehensive analysis of the expression, prognostic significance, and regulation pathway of G2E3 in breast cancer

Background Loss of G2-specific E3-like (G2E3) protein sensitizes tumor cells to chemotherapy. However, the role of G2E3 in breast cancer development and patient’s prognosis is unclear. Here, we explored the expression, prognostic significance, and regulatory pathway of G2E3 in breast cancer. Methods TCGA and UALCAN database were utilized to explore G2E3 expression in breast cancer and normal tissues and its expression in breast cancer based on clinicopathological characteristics, respectively. The Kaplan–Meier plotter database was utilized to determine the effect of G2E3 on the prognosis of breast cancer patients. RT-PCR was utilized to validate the G2E3 expression in cancerous and normal breast tissues. Immunohistochemistry analysis was utilized to validate the prognostic effect of G2E3 expression in breast cancer patients and the relationship between G2E3 expression and lymphocyte infiltration levels. Receiver operating characteristic (ROC) curves were also generated to validate the diagnostic value of G2E3 expression in recurrence/distant organ metastasis and death. The STRING database, DAVID database, and Sanger-box tools were utilized to perform GO functional, KEGG pathway enrichment, and GSEA analysis. The TISIDB database was utilized to determine the relationship between G2E3 expression and tumor immunity. Finally, CTD database was utilized to screen for potential therapeutic compounds that could reduce the G2E3 mRNA expression. Results TCGA data presented that G2E3 expression was higher in breast cancer tissues than in normal breast tissues. This result was further validated by RT-PCR (P = 0.003). The Kaplan–Meier plotter database suggested that patients with high G2E3 mRNA expression had significantly shorter RFS and OS than patients with low G2E3 mRNA expression. Immunohistochemistry analysis of 156 breast cancer clinical specimens also validated patients with G2E3-positive expression had a significantly shorter DFS and OS than patients with G2E3-negative expression. Thus, G2E3 expression was an independent prognostic predictor of DFS and OS. The G2E3-positive expression also has a high diagnostic value for recurrence/distant organ metastasis and death. GSEA analysis revealed that G2E3 might be enriched in the E2F, PI3K/AKT/mTOR signaling, DNA repair pathways, and other cancer-related signaling pathways. The TISIDB database showed that G2E3 expression was significantly negatively associated with lymphocyte infiltration. This result was further validated in clinical breast cancer samples (P = 0.048; R = −0.158). Using the CTD database, we found that (+)-JQ1 compound, 1,2-dimethylhydrazine, and other compounds may decrease the G2E3 mRNA expression. These compounds could serve as potential therapeutic compounds for the clinical treatment of breast cancer. Conclusions G2E3 expression was higher in breast cancer tissues than in normal tissues. G2E3-positive expression was related to a worse survival outcome in patients with breast cancer. Genes co-expressed with G2E3 may be enriched in the breast cancer-related signaling pathways. The G2E3 expression was significantly negatively associated with lymphocyte infiltration. G2E3 may serve as a novel prognostic biomarker and therapeutic target for breast cancer. Supplementary Information The online version contains supplementary material available at 10.1186/s12957-022-02871-0.


Introduction
Breast cancer has been a malignant tumor with the highest morbidity in women. New breast cancer cases account for approximately 30% of the total number of new malignant tumors in women every year, which is a significant threat to their health [1]. Although breast cancer molecular typing and the emergence of comprehensive treatment methods have greatly improved the patient survival outcomes [2,3], local recurrence and distant metastasis remain the leading causes of breast cancer-related death [4,5]. Breast cancer was initially considered a local disease. However, they can migrate into lymph nodes and distant target organs, such as the bone, lung, and liver [6][7][8][9]. Current adjuvant therapy for breast cancer can also eliminate tumor cells that have spread to distant sites at the time of diagnosis and significantly improve the 10-year survival rate of women. However, about 40% of women who have already received adjuvant therapy would develop postoperative metastasis and eventually die of metastatic breast cancer [5,10]. Therefore, identifying the intrinsic driver genes of tumor metastasis, screening effective therapeutic targets, and developing targeted drugs may prevent tumor metastasis from the root cause [7,11].
G2-specific E3-like (G2E3) protein was first identified by Brooks et al. [12]. It is denoted as G2E3 as it is a G2-specific protein with a domain similar to many E3 proteins. Schmidt et al. found that G2E3 may be a potential target for chemo-sensitizing tumor cells and might participate in response to cisplatin [13]. G2E3 knockdown may promote apoptosis and inhibit cancer cell proliferation. G2E3 can be a molecular determinant of the DNA damage response (DDR) and cell survival, and its loss of expression can sensitize tumor cells to chemotherapy. However, the role of G2E3 in breast cancer development and patient's prognosis is still unclear. Therefore, we explored the prognostic significance of G2E3 expression in breast cancer patients and the effect of G2E3 on the malignant biological behavior of cancer cells.
This study primarily analyzed the G2E3 expression in cancerous and normal breast tissues, explored the effect of G2E3 high expression on patient' survival outcomes, and validated it using clinical breast cancer samples. Furthermore, we explored genes co-expressed with G2E3 and performed GO functional and KEGG pathway enrichment analyses. Following this, we analyzed the relationships between G2E3 expression, tumor immunity, and mutations of G2E3 in breast cancer. Finally, potential therapeutic compounds that can decrease the G2E3 mRNA expression in breast cancer were screened. Therefore, we identified an effective therapeutic target for breast cancer.

CCLE database
CCLE database (https:// sites. broad insti tute. org/ ccle) was utilized to determine the G2E3 expression in a series of cancer cell lines [15].

UALCAN database
UALCAN database (ualcan.path.uab.edu/index.html) was utilized to explore the G2E3 expression in breast cancer patients based on sample type, patient age, individual cancer stage, menopause status, nodal metastasis status, and breast cancer subclasses [16].

Kaplan-Meier plotter database
The Kaplan-Meier plotter database (http:// kmplot. com) was utilized to determine the effect of high G2E3 mRNA expression on relapse-free survival (RFS) and overall survival (OS) in breast cancer patients [17].

RT-qPCR
Total RNA from 30 pairs of fresh cancerous and normal breast tissues were isolated using TRIzol solution (Solarbio Company). The extracted RNA was reversetranscribed by a cDNA synthesis kit (TaKaRa). qPCR was performed by SYBR Green PCR Master Mix (TaKaRa) and primers binding to the G2E3 and GAPDH. These primers were designed by Shanghai Sangon Biotech Co., Ltd. (G2E3: forward, 5′-CAG CAC TAT GAG CGT TGT  GAT GTT C-3′ and reverse 5′-ACC GTA ATG AGG AGC  AGG CTA AAT G-3′; GAPDH: forward, 5′-CCT TCC GTG  TCC CCACT-3′ and reverse, 5′-GCC TGC TTC ACC ACC  TTC -3′). The cycling protocol was 95 °C for 30 s (initial denaturation), followed by 40 denaturation cycles at 95 °C for 3 s, and finally annealing and extension at 60 °C for 30 s. Relative mRNA levels were calculated using the 2 −ΔΔCt method [18]. These breast cancer specimens were obtained from affiliated hospital of China Medical University at the time of surgery, and the basic clinicopathological characteristics of the 30 patients were presented in Supplementary Table 1. This study was approved by China Medical University Institutional Review Board.

Patients and breast cancer specimens
This study also included 156 patients diagnosed with invasive ductal carcinoma between January 2007 and November 2011 at affiliated hospital of China Medical University. All patients underwent routine operation and postoperative treatment. The inclusion criteria were as follows: (1) complete clinicopathological patient information was required and (2) no organ metastasis at the time of surgery. The exclusion criteria were as follows: (1) patients who did not have complete clinicopathological information, (2) patients who did not receive surgery or routine postoperative treatment, and (3) patients with unknown survival status. All patients were followed up for at least 10 years. Disease-free survival (DFS) was defined from the date of the operation to the date of recurrence/ metastases. OS was defined as the time from operation to death. The survival status of patients was determined using interviews/telephone calls.

Immunohistochemistry analysis
The clinical specimens obtained from 156 patients were fixed and embedded in paraffin. Then, they were sliced into 4-μm sections and were deparaffinized and rehydrated. The sections were incubated with an antibody against G2E3 (1:200; Bioss, bs-17000R) at 4 °C overnight. Then, these sections were incubated with a secondary antibody (Zhong Shan Jin Qiao) on the 2nd day at room temperature.
Evaluation of the levels of tumor-infiltrating lymphocytes was performed as previously described [19].

STRING analysis and GO functional and KEGG pathway enrichment
The STRING database (https:// string-db. org/) was applied to identify genes co-expressed with G2E3 [20]. Protein-protein interactions (PPI) were also determined. The DAVID database (https:// david. ncifcrf. gov) was utilized for GO functional and KEGG pathway enrichment analyses of genes co-expressed with G2E3 [21].

Sanger-box tools
Sanger-box tools (http:// www. sange rbox. com/ tool) were used for the following analysis: Gene Set Enrichment Analysis (GSEA); the relationships between G2E3 expression and ImmuneScore, StromalScore, and Immu-neScore; the relationships between G2E3 expression and immune checkpoint gene expression; the relationships between G2E3 expression and the number of tumor neoantigens, tumor mutational burden (TMB), microsatellite instability (MSI), and the mutation pattern of G2E3 in breast cancer; and the correlations between the expression of G2E3 in each tumor sample and that of the DNA repair genes (MMRs) and methyltransferase.

Comparative Toxicogenomics Database (CTD)
CTD is an innovative online database that provides literature-based data on the interactions between oncogenes and chemotherapeutic compounds. This tool was used to screen potential therapeutic compounds that could reduce G2E3 mRNA expression [23].

Statistical analysis
Correlations between G2E3 expression and age, tumor (T) grade, lymph nodes metastases (N) grade, histological grade, subtypes, menopausal status, relapse/metastasis, and death were analyzed using the chi-square test. Cox regression analyses were performed to identify independent predictors of DFS and OS. Hazard ratios (HRs) and 95% confidence intervals (CIs) were also calculated. The correlations between G2E3 expression and tumorinfiltrating lymphocytes (TILs) were analyzed using Pearson correlation. Survival curves were generated using the Kaplan-Meier test in SPSS 25.0. ROC curves were generated using SPSS 25.0 software. The significance was set at P < 0.05. Figure 1A presents the G2E3 expression in 31 normal human tissues using the GTEx database. We analyzed the data for a series of cell lines obtained from CCLE database. According to the organization source, the data could be classified into 21 organizations. We analyzed G2E3 expression in these 21 tissues (Fig. 1B). We then obtained the difference in G2E3 expression between cancer tissues and normal tissues in each tumor sample from The Cancer Genome Atlas (TCGA) database (Fig. 1C). Because there were fewer normal breast samples in TCGA, we further integrated the data of normal breast tissues from the GTEx database and the data of TCGA tumor tissues to analyze the G2E3 expression differences among the 27 tumor types (Fig. 1D). Our interest lies in breast cancer research. From the analysis depicted in Fig. 1 A and B, we found that the G2E3 gene has moderate expression in breast tissues compared with other human tissues. The expression of G2E3 was significantly higher in breast cancer tissues than in normal tissues, shown in Fig. 1 C and D.

Expression and prognostic effect of G2E3 in breast cancer explored by UALCAN and Kaplan-Meier plotter databases
Using the UALCAN database, we also found the expression of G2E3 mRNA was significantly higher in breast cancer than in normal tissues (P = 1.96E−03; Fig. 2A). In the subgroup of patients, those aged 21-40 years had the highest G2E3 mRNA expression (P = 4.68E−03;   Fig. 2C). In the subgroup of patients with menopausal status, those who were premenopausal had the highest G2E3 mRNA expression (P = 5.79E−04; Fig. 2D). In the subgroup of patients with nodal metastasis, N2 patients had the highest G2E3 mRNA expression (P = 1.60E−04; Fig. 2E). Finally, in the subgroup of patient subclasses, patients with luminal cancer had the highest G2E3 mRNA expression (P = 1.24E−03; Fig. 2F).

Validation of the expression and prognostic effect of G2E3 in breast cancer clinical specimens
First, we validated the expression of G2E3 mRNA in 30 pairs breast cancer and normal tissues by RT-PCR. We found G2E3 mRNA expression was significantly higher in breast cancer tissues than in normal tissues (P = 0.003; Fig. 3A). We explored the G2E3 expression by immunohistochemistry in 156 breast cancer specimens, presenting as G2E3-positive and -negative cells (Fig. 3B). We found that patients with G2E3-positive expression had a significantly shorter DFS (P < 0.01; Fig. 3C) and OS (P = 0.02; Fig. 3D) than patients with G2E3-negative expression. Therefore, we can conclude that G2E3-positive expression has adverse prognostic effects on the prognosis of breast cancer patients. We also explored the diagnostic value of G2E3 expression in recurrence/distant organ metastasis and death. For recurrence/distant organ metastasis, G2E3 expression had a high diagnostic value, wherein the area under the curve (AUC) was 0.7186 (P < 0.001; Fig. 3E). G2E3 expression also had a high diagnostic value for death, with an AUC of 0.7468 (P < 0.001; Fig. 3F).
We further analyzed the relationships between G2E3 expression and age, T grade, N grade, histological grade, subtypes, menopausal status, relapse/metastasis, and death in these 156 patients with breast cancer. We observed that G2E3-positive expression was related to a higher N grade, although the difference was not significant (P = 0.055). Moreover, G2E3 expression was  Table 1.
To identify independent predictors of prognosis, we performed Cox regression analysis. These results were displayed in Tables 2 and 3. For DFS, N grade (P < 0.001) and G2E3 expression (P = 0.004) were considered independent predictors of DFS in these 156 patients with breast cancer. For OS, the N grade was an independent predictor of OS (P = 0.003). However, G2E3 expression was not a significant independent predictor of overall with OS (P = 0.096). These results suggest that G2E3 is an independent predictor of DFS, but not OS.

PPI network construction and GO functional and KEGG pathway enrichment analysis
Using the STRING database, we constructed a PPI network of co-expressed genes with G2E3 (Fig. 4A). These genes are IAH1, LRRC17, ARHGEF39, and PRR11. We then performed functional and pathway enrichment analyses based on these genes using the DAVID database. GO functional enrichment analysis consisted of biological process (BP), cellular component (CC), and molecular function (MF). For BP, these genes were mainly enriched in processes such as mitotic cells' metaphase/anaphase transition and regulation of the metaphase/anaphase transition of the cell cycle (Fig. 4B). For CC, these genes were mainly enriched in components such as meiotic spindles, condensed chromosome outer kinetochores, and intercellular bridges (Fig. 4C). For MF, these genes were mainly enriched in histone kinase activity and protein serine/threonine/tyrosine kinase activity (Fig. 4D). KEGG pathway analysis suggested that they were mainly enriched in cell cycle and oocyte meiosis (Fig. 4E).

Relationships between G2E3 and tumor immunity in breast cancer
Tumor immunity has been found to be involved in the development of breast cancer [24][25][26]. Therefore, we evaluated the correlations between G2E3 and immunity in breast cancer. Using the TISIDB database, we first explored the correlation between G2E3 expression and lymphocyte infiltration (Fig. 5A). The top four lymphocytes negatively correlated to G2E3 expression most significantly were CD56 bright, CD56dim, monocyte, and pDC (Fig. 5B). We further experimentally validated the relationship between G2E3 expression and TILs levels in breast cancer clinical samples. We found G2E3 expression was significantly negatively associated with the TILs levels (P = 0.048; R = −0.158; Table 5). Typical images of G2E3 expression and TILs were shown in Fig. 5 C and D.
We also explored the correlations between G2E3 expression and immunomodulators and chemokines by the TISIDB database. Supplementary Fig. 2A shows the correlation between G2E3 expression and immuneinhibitors. The top four immune-inhibitors most significantly negatively associated with G2E3 expression were LGALS9, PDCD1, TGFβ1, and TGFBR1 ( Supplementary  Fig. 2B). Supplementary Fig. 2C shows the correlations between G2E3 expression and immune-stimulators. The top four immune-stimulators most significantly negatively associated with G2E3 expression were C10orf54, TNFRSF14, TNFRSF25, and TNFSF4 ( Supplementary  Fig. 2D). Supplementary Fig. 2E shows relationships   2F). Supplementary Fig. 3A shows the correlation between G2E3 expression and chemokines. The top four chemokines most significantly negatively related to G2E3 expression were CCL14, CCL19, CCL21, and CX3CL1 ( Supplementary Fig. 3B). Supplementary  Fig. 3C showed the correlation between G2E3 expression and chemokine receptors. The top four chemokine receptors most significantly negatively associated with G2E3 expression were C-C motif chemokine receptor (CCR)7, CCR10, CXCR3, and CXCR5 ( Supplementary Fig. 3D). We also explored the expression of G2E3 in different  Fig. 3E), similar to the results obtained from the UALCAN database (Fig. 2F). The expression of G2E3 was similar in all immune subtypes of breast cancer (Supplementary Fig. 3F). An increasing number of reports have indicated that the tumor immune microenvironment plays a vital role in the development of cancer [27][28][29]. Sanger-box tools provided the analysis of the immune and stromal scores of every tumor sample to observe the correlations between gene expression and immune scores. These data included ImmuneScore, correlations between gene expression and matrix score such as StromalScore, and the correlations between gene expression and ESTIMATE immune score such as ESTIMATE score in 33 tumors. We found that in breast cancer, G2E3 expression was significantly negatively associated with ImmuneScore (Supp. Fig. 4) and ESTIMATE score (Supplementary Fig. 6) but was not significantly related to StromalScore (Supplementary Fig. 5).
Under normal conditions, the immune system can effectively recognize and eliminate tumor cells [30]. However, tumor cells can adopt many strategies to suppress the body's immune system and prevent tumor cells from being eliminated [31]. All stages of the immune response survived. Tumor immunotherapy is a treatment method that can restore a normal immune response in the body [32,33]. Figure 5E depicts the analysis from Sangerbox tools that collected nearly 40 immune checkpoint genes and explored the correlations between G2E3 and immune checkpoint gene expression, such as TNFRSF14, NRP1, and CD44.
Tumor neoantigens are abnormal proteins encoded by genetic point and deletion mutations in cancer cells [34,35]. Using the immune activity of tumor neoantigens, neoantigen vaccines can be designed and synthesized according to mutations in tumor cells, and patients can be immunized to achieve therapeutic effects. Here, Sanger-box tools count the number of neoantigens in every tumor sample and analyze the correlation between G2E3 gene expression and the number of antigens. Unfortunately, there was no significant relationships between G2E3 expression and the number of breast cancer neoantigens ( Supplementary  Fig. 7). Sanger-box tools also revealed the relationship between G2E3 gene expression and TMB and MSI in breast cancer. These results are shown in Fig. 5F and G.

Mutations of G2E3 gene in breast cancer
Using Sanger-box tools, mutations in G2E3 in breast cancer were also explored. The somatic mutation rate of G2E3 in breast cancer was 0.41% (Fig. 6A).  DNA repair genes are mismatched repair genes in cells. The loss of function of critical genes in this mechanism would cause irreparable DNA replication errors, leading to increased somatic mutations [36]. Sanger-box tools used the expression profile data to evaluate the correlation between five DNA repair genes: MLH1, MSH2, MSH6, PMS2, EPCAM mutation, and G2E3 expression. The G2E3 expression was significantly positively related to mutations of these five genes in breast cancer (Fig. 6B).
DNA methylation can cause changes in the chromatin structure, DNA conformation, and other factors, thereby controlling gene expression. DNA methylation can only be achieved using DNA methyltransferases [37]. Therefore, we analyzed the relationship between G2E3 expression and the expression of the four methyltransferases. In breast cancer, the G2E3 expression was significantly positively related to the expression of these four methyltransferases (Fig. 6C).

Screening potential therapeutic compounds which can decrease G2E3 mRNA expression for breast cancer
By applying the CTD database, we aimed to screen for potential therapeutic compounds that could decrease the mRNA expression of the G2E3 oncogene in breast cancer patients. We found that (+)-JQ1 compound, 1,2-dimethylhydrazine, and other compounds can decrease the G2E3 mRNA expression, which may be potential therapeutic compounds for breast cancer. The results are presented in Table 6. The 3D structures of the top ten compounds were obtained from the PubChem database (https:// pubch em. ncbi. nlm. nih. gov) and are shown in Fig. 7.

Discussion
Brooks et al. first identified the G2E3 protein with a HECT (homologous to E6-associated protein) domain located on the C-terminus, suggesting that it might serve as a ubiquitin ligase or E3 [12]. They further found that G2E3 functions as a nucleocytoplasmic shuttling protein.
Its export relies on sequences in the carboxy-terminal half of the protein and is not dependent on the regular CRM1/Exportin 1-mediate process. Following DNA damage, G2E3 rapidly delocalizes from the nucleoli into the cytoplasm. The cell cycle phase-specific expression and rapidly induced subcellular localization of G2E3 play an important role in regulating the cell cycle and the DDR. Owing to the rapid response to DNA damage, Schmidt et al. found that G2E3 participates in the process by which chemotherapy drugs kill tumor cells [13]. One of the main functions of chemotherapeutic drugs that kill tumor cells is DNA damage. They result in DNA doublestrain breaks in cancer cells, which are often associated with cell death [38][39][40]. Phosphorylation of histone variant H2AX at Ser139 (γH2AX) is a marker of DDR [41]. However, in the process of cisplatin treatment and the ubiquitin system, the decrease in γH2AX is a marker of DDR. Schmidt et al. found that in cisplatin treatment, G2E3 depletion may decrease the phosphorylation of H2AX. G2E3 was required for transmitting the DDR signal to H2AX when cells were treated with cisplatin. Loss of G2E3 can cause p53 accumulation, cell apoptosis, and inhibition of tumor cell proliferation. When cells are depleted of the nucleoside analogs gemcitabine and G2E3, DNA replication can be further hindered. Simultaneously, the G2E3 expression decreased when tumor cells were treated with chemotherapeutic drugs, such as cisplatin. Therefore, they concluded that G2E3 might serve as a novel modulator of DDR, and its loss of expression can sensitize tumor cells to DNA damage. These findings predict that G2E3 is an oncogene for breast cancer. These findings are consistent with our results.
This study found G2E3 expression was significantly higher in cancer tissues than in normal tissues, validated using an online database and PCR. By the Kaplan-Meier plotter database, we found G2E3 high expression was significantly related to shorter RFS and OS. We further validated this by breast cancer clinical specimens. We also found G2E3-positive expression had a high diagnostic value for DFS and OS by ROC curves and was an independent prognostic predictor for DFS in breast cancer patients. Pathway enrichment analysis of G2E3 revealed that it was enriched in "HALLMARK_E2F_TARGETS, " "HALL-MARK_PI3K_AKT_MTOR_SIGNALING, " and other cancer-related pathways. We found that the expression of G2E3 was significantly negatively correlated with lymphocytes infiltration, which was also validated in clinical breast cancer specimens. We also explored the relationships between G2E3 expression and the tumor immune microenvironment and tumor neoantigens. Finally, we explored the mutation of the G2E3 gene in breast cancer and screened for potential therapeutic compounds that can decrease the G2E3 mRNA expression. We hypothesized that G2E3 could serve as a novel therapeutic target for breast cancer. However, there are still some minor issues with the results of the online database. From the results of Fig. 2 C, E, G, and H, we can find that G2E3 expression had a significant adverse effect on RFS and OS, but either stage 3 or N2 patients had the highest G2E3 expression, not stage 4 nor N3 patients which probably correlated with worse prognosis than stage 3 nor N2 patients [42][43][44][45]. We think the results may owe to the few numbers of breast cancer patients included in Fig. 2 C and E. Figure 2C only included 20-stage 4 breast cancer patients, and Fig. 2E only included 77 N3 breast cancer patients. The worse prognostic effects of G2E3 mRNA high expression on RFS and OS presented in Fig. 2 G and H were gained based on 2032 and 943 breast cancer patients respectively. Therefore, if included more stage 4 or N3 breast cancer patients, we may get more reasonable results.
Studies have found that G2E3 is critical for early embryonic development, especially for maintaining germline stem cells in Drosophila [46]. Brooks et al. found that G2E3 is an unusual ubiquitin ligase essential for early embryonic [47]. Li et al. found G2E3 methylation may affect embryo diapause by regulating the cell cycle [48]. Powell et al. found that G2E3 knockout can also increase the bodyweight of mice [49]. This finding suggests that G2E3 may regulate obesity. Finally, Zhang et al. found that cathepsin D can enhance invasion and metastasis in breast cancer by promoting hepsin ubiquitin-proteasome degradation. However, this regulation is mediated by cathepsin D [50]. Cathepsin D overexpression can significantly increase G2E3 expression, and blocking G2E3 expression can significantly inhibit hepsin degradation induced by cathepsin D. Immunohistochemical analysis also found G2E3 expression was significantly negatively related to the expression of hepsin and positively related to cathepsin D expression in breast cancer. This result suggests that G2E3 may regulate the malignant behavior of breast cancer cells.
This study also had some limitations. First, we did not validate the function of G2E3 using in vitro or in vivo experiments. Second, we did not explore the exact regulatory mechanisms of G2E3 in breast cancer development. Finally, we did not determine a relationship between G2E3 expression and breast cancer immunity. These limitations are the future study directions for our experiments.

Conclusions
In conclusion, we found that G2E3 was more highly expressed in breast cancer tissues than in normal breast tissues. G2E3-positive expression was related to a worse prognosis in patients with breast cancer. Genes co-expressed with G2E3 may be enriched in the E2F and PI3K/AKT/mTOR signaling pathways and other cancer-related signaling pathways. The G2E3 expression was significantly negatively correlated with lymphocyte infiltration. Thus, G2E3 may serve as a novel prognostic biomarker and therapeutic target for breast cancer.