Differential gene expression analysis between cancer tissues and adjacent normal tissues was performed using The Cancer Genome Atlas (TCGA; https://cancergenome.nih.gov/) database and Gene Expression Profiling Interactive Analysis (GEPIA; http://gepia.cancer-pku.cn/) database. In the present study, the mRNA expression data of FGGY in both CRC and matched non-cancerous tissues were extracted from the GEO datasets (GSE20916 and GSE113513) (5,13) and TCGA (TCGA-COAD datasets) databases (14). Pearson correlation analysis was performed to evaluate the correlation between FGGY and p53 expression at the mRNA level in CRC using the GSE39582 GEO dataset (15) and TCGA-COAD datasets. Furthermore, the mRNA expression levels of FGGY in different types of digestive system cancer compared with in tissues from healthy controls were assessed using GEPIA (COAD datasets) (16).

The CRC cell lines, HCT-8 (cat. no. TCHu18), Caco2 (cat. no. TCHu146), SW480 (cat. no. TCHu172), HCT116 (cat. no. TCHu99) and RKO (cat. no. TCHu116), and the colorectal adenocarcinoma cell line HT-29 (cat. no. TCHu103) were obtained from the National Collection of Authenticated Cell Cultures. The wild-type HCT116 cell line (HCT116/p53^+/+) and p53-knockout HCT116 cell line (HCT116/p53^−/−) were provided by Dr Yao Lin (Fujian University of Traditional Chinese Medicine, Fuzhou, China), which were originally sourced from Dr. Bert Vogelstein (Johns Hopkins University, Baltimore, MD, USA). HT-29 cells were cultured in McCoy's 5A medium (Nanjing KeyGen Biotech Co., Ltd.), whereas RKO cells were grown in MEM-α medium (Thermo Fisher Scientific, Inc.), Caco2 cells were maintained in DMEM (Thermo Fisher Scientific, Inc.), and HCT-8, HCT116, HCT116/p53^+/+ and HCT116/p53^−/− cells were all maintained in RPMI1640 medium (Thermo Fisher Scientific, Inc.). Each cell line was incubated at 37°C in a humidified atmosphere containing 5% CO₂. The growth media for all cells were supplemented with 10% FBS (Thermo Fisher Scientific, Inc.), 100 U/ml penicillin and 100 μg/ml streptomycin (Hyclone; Cytiva). The HT-29 cell line was characterized by Suzhou Jianda Biotechnology Co., Ltd. using short tandem repeat profiling.

The total RNA was extracted from the cultured cells using RNAiso Plus reagent (Takara Biotechnology Co., Ltd.). RT into cDNA was performed using the PrimeScript RT reagent kit (Takara Biotechnology Co., Ltd.), strictly adhering to the guidelines provided by the manufacturer. In addition, a cDNA array of tissue samples, including 79 primary CRC and 15 adjacent noncancerous tissues (cat. no. HColA095Su01), was purchased from Shanghai Outdo Biotech Co., Ltd. The mRNA expression levels of FGGY and GAPDH were quantified using an ABI 7500 Fast Real-Time PCR System (Applied Biosystems; Thermo Fisher Scientific, Inc.) with SYBR Premix Ex Tag (Takara Biotechnology Co., Ltd.). The qPCR protocol consisted of an initial denaturation at 95°C for 10 min, followed by 40 cycles of denaturation at 95°C for 15 sec and annealing/extension at 60°C for 60 sec. GAPDH was used as the internal control. The primer sequences are listed in Table SI. The mRNA expression levels were determined using the following formula: The relative expression levels were calculated using the 2^−ΔΔCq method, where ΔΔCq=[Cq(target gene)-Cq(GAPDH)]sample-[Cq(target gene)-Cq(GAPDH)]control (17). The experimental procedures and protocols of this research were approved by the Human Ethics Committee of Shanghai Outdo Biotech Co., Ltd. (no. SHYJS-CP-1704010). The clinicopathological characteristics of the patients with CRC included in the cDNA array are presented in Table SII.

A TMA of CRC samples (including 101 CRC tissues and 79 adjacent noncancerous tissues; cat. no. HColA180Su15) and digestive system neoplasms [including colon adenocarcinoma (COAD), rectal adenocarcinoma (READ), pancreatic cancer (PAAD), esophageal cancer (ESCA), gastric cancer (STAD) and liver hepatocellular carcinoma (LIHC); including tumor tissues (COAD: 16; READ: 14; PAAD: 15; ESCA: 14; STAD: 14; LIHC: 16), normal tissues from healthy controls (colorectal tissues: 9; rectal tissues: 3; pancreatic tissues: 3; esophageal tissues: 5; gastric tissues: 7; liver tissues: 3) and adjacent noncancerous tissues (colorectal tissues: 16; rectal tissues: 14; pancreatic tissue: 15; esophageal tissues: 14; gastric tissues: 14; liver tissues: 16); cat. nos. HDgsC140PT01 and HorgC120PG04] were obtained from Shanghai Outdo Biotech Co., Ltd. IHC was performed to evaluate FGGY expression in clinical tumor samples, following previously described methods (5). Formalin-fixed, paraffin-embedded TMA sections were dehydrated using graded ethanol and antigen retrieval was performed by microwave heating the slides for 20 min in sodium citrate-hydrochloric acid buffer (Fuzhou Maixin Biotechnology Development Co., Ltd.). After permeabilization with 0.1% Triton X-100 (cat. no. T8200; Beijing Solarbio Science & Technology Co., Ltd.) for 10 min, the sections were then treated with 3% H₂O₂ for 10 min at room temperature, followed by incubation with a blocking reagent (cat. no. KIT-9710; Fuzhou Maixin Biotechnology Development Co., Ltd.) for 30 min at room temperature. Tissue sections were then incubated with an antibody against FGGY (1:900; cat. no. PA5-120316; Thermo Fisher Scientific, Inc.) overnight at 4°C, and with a secondary antibody (cat. no. KIT-9710; Fuzhou Maixin Biotechnology Development Co., Ltd.) for 1 h at room temperature. Immunoreactivity was then visualized using an enhanced DAB kit (cat. no. DAB-2031; Fuzhou Maixin Biotechnology Development Co., Ltd.). Images were captured using a Nano Zoomer 2.0 HT light microscope slide scanner (Hamamatsu Photonics K.K.) and processed using Nano Zoomer Digital Pathology View 1.6 software (Hamamatsu Photonics K.K.). To assess FGGY expression, a grading system was used; staining intensity was rated as follows: No staining, 0; weak, 1; moderate, 2; strong, 3; and the proportion of cells showing positive staining was rated as: 1-25%, 1; 26-50%, 2; 51-75%, 3; 76-100%, 4. The overall score for FGGY expression was determined by multiplying the intensity score with the percentage score.

For survival analysis, FGGY expression in CRC tissues was categorized into either low (final score, 0-6) or high (final score, 7-12). Kaplan-Meier survival curves were plotted for groups with high and low FGGY expression, and were analyzed using log-rank test. The experimental procedures and protocols of this research were approved by the Human Ethics Committee of Shanghai Outdo Biotech Co., Ltd. [nos. SHYJS-CP-1810005 (HColA180Su15), SHYJS-CP-1904003 (HDgsC140PT01), SHYJS-CP-1901007 (HorgC120PG04)]. The clinicopathological characteristics of patients with CRC are presented in Table SIII.

The specific short hairpin RNA (shRNA) targeting FGGY and a non-silencing control shRNA were designed and validated by Shanghai GeneChem Co., Ltd. The shRNA oligonucleotides were cloned into a GV248 lentiviral vector (Shanghai GeneChem Co., Ltd.) that included a green fluorescence protein (GFP) reporter gene. Briefly, a second-generation system was used and 293T cells were used as the interim cell line. For plasmid transfection, 293T cells were subcultured and seeded into 10-cm plates at a density of 5×10⁶ cells/15 ml. The DMEM was replaced with serum-free medium 2 h before transfection. A DNA solution containing 20 μg GV248 vector plasmid, 15 μg pHelper 1.0 carrier plasmid, 10 μg pHelper 2.0 vector plasmid and the corresponding volume of transfection reagent (cat. no. GRCT105; Shanghai GeneChem Co., Ltd.) was evenly mixed. The total volume was adjusted to 1 ml and incubated at room temperature for 15 min, then added dropwise to the 293T cells. After 6 h of transfection, the medium was replaced and 10 ml PBS was added. Subsequently, 20 ml fresh medium was added, and the cells were incubated at 37°C in a humidified atmosphere containing 5% CO₂ for an additional 48-72 h. The supernatant from the 293T cells was then collected, centrifuged at 4,000 × g for 10 min at 4°C to remove cell debris and filtered through a 0.45-μm filter. The supernatant was then centrifuged at 25,000 × g for 2 h at 4°C; after centrifugation, the supernatant was discarded, and viral storage solution was added to completely dissolve the lentiviral particles. For cell transduction, HCT116 and RKO cells were subcultured, seeded into 12-well plates at a density of 0.4×10⁵ cells/well and were cultured overnight. Lentivirus coding shRNA targeting FGGY (LV-FGGY-RNAi; hU6-MCS-Ubiquitin-EGFP-IRES-puromycin) or the control shRNA (sh-Ctrl; hU6-MCS-Ubiquitin-EGFP-IRES-puromycin) (18) was added at a multiplicity of infection (MOI) of 10, with medium and Hitrans P infection enhancement reagent (cat. no. REVG005; Shanghai GeneChem Co., Ltd.) reaching a total volume of 500 μl/well. The specific sequences of these shRNAs are detailed in Table SIV. The medium was changed after 12 h of transduction and fresh medium was added. The cells were then incubated at 37°C in a humidified atmosphere containing 5% CO₂ for an additional 72 h. After 72 h, the infection efficiency was verified by observing the number of GFP-positive cells under a fluorescence microscope. To confirm the knockdown of FGGY, the mRNA and protein expression levels of FGGY were assessed using RT-qPCR and western blot analysis, respectively. After verifying a significant decrease in FGGY expression at both the mRNA and protein levels, the cells were collected for subsequent experiments.

To investigate the effects of FGGY overexpression on cell proliferation, HT-29 cells were transduced with a lentiviral vector (GV492; Ubi-MCS-3FLAG-CBh-gcGFP-IRES-puromycin) (MOI, 10) that encoded the full-length human FGGY gene [NCBI Gene ID: 55277; OBiO Technology (Shanghai) Corp., Ltd.] or an empty vector control (GV492 without FGGY insert) at the same MOI. The protocol for the overexpression plasmid transduction was the same as that described for shRNA transduction. The medium was changed after 12 h of transduction and fresh medium was added. The cells were then incubated at 37°C in a humidified atmosphere containing 5% CO₂ for an additional 72 h. After 72 h, the infection efficiency was verified by observing the number of GFP-positive cells under a fluorescence microscope. Successfully transduced cells were selected using 2 μg/ml puromycin (cat. no. A1113803; Thermo Fisher Scientific, Inc.) for 1 week, and then maintained in medium containing 1 μg/ml puromycin. To confirm the overexpression of FGGY, the mRNA and protein expression levels of FGGY were assessed using RT-qPCR and western blot analysis, respectively. After verifying a significant decrease in FGGY expression at both the mRNA and protein levels, the cells were collected for subsequent experiments.

After lentiviral transduction and puromycin selection, cells were collected and lysed using RIPA lysis buffer (Beyotime Institute of Biotechnology) supplemented with 1 mM PMSF and protease inhibitors. The total protein concentration was determined using the Pierce BCA Protein Assay Kit (Thermo Fisher Scientific, Inc.). Protein lysates (50 μg) were separated by SDS-PAGE on 10% gels and were then transferred onto polyvinylidene difluoride membranes. These membranes were subsequently blocked with 5% nonfat milk in TBS-0.1% Tween for 2 h at room temperature, and then incubated with primary antibodies against FGGY (cat. no. PA5-120316; Thermo Fisher Scientific, Inc.), p53 (cat. no. 10442-1-AP; Proteintech Group, Inc.), p21 (cat. no. 2947S; Cell Signaling Technology, Inc.), PCNA (cat. no. 10205-2-AP; Proteintech Group, Inc.), Bcl-2 (cat. no. 3498; Cell Signaling Technology, Inc.), Bax (cat. no. 2772; Cell Signaling Technology, Inc.) or GAPDH (cat no. Abp57259; Abbkine Scientific Co., Ltd.) (all at a dilution of 1:1,000) at 4°C overnight. This was followed by incubation with goat anti-rabbit secondary antibodies conjugated with horseradish peroxidase (1:5,000; cat. no. 7074; Cell Signaling Technology, Inc.) for 2 h at room temperature. Protein bands were detected using an ECL imager (cat no. 32209; Thermo Fisher Scientific, Inc.) and their intensities were analyzed with ImageLab software (version 5.2.1; Bio-Rad Laboratories, Inc.). GAPDH expression was used as a loading control. Each experiment was performed in triplicate.

The cells were seeded in 96-well plates at a density 0.2×10⁵/ml and were cultured for 24, 48, 72, 96 or 120 h, after which 10 μl CCK-8 solution (cat. no. KTA1020; Abbkine Scientific Co., Ltd.) was added to each well. The plates were then incubated at 37°C for 2 h in an incubator with 5% CO₂ in the dark. Subsequently, the optical density of each well was measured at an absorbance of 450 nm using a microplate reader.

For the cell cycle analysis, transduced cells were harvested and fixed in 70% ethanol at 4°C overnight. These fixed cells were then centrifuged at 2,000 × g for 3 min at room temperature, washed and incubated with FxCycle™ PI/RNase staining solution (cat no. F10797; Thermo Fisher Scientific, Inc.) for 30 min at room temperature. Cell cycle progression was analyzed using a FACSCalibur flow cytometer (BD Biosciences), with Modfit LT version 3.0 software (Verity Software House, Inc.). For the apoptosis analysis, transfected cells were rinsed twice with ice-cold PBS and incubated with Annexin-V-AbFluor™ 647 Apoptosis Detection Kit (Abbkine Scientific Co., Ltd.) for 30 min at 37°C in the dark. The rate of apoptosis in these cells was also determined using FACSCalibur flow cytometer (BD Biosciences). Data were analyzed using CellQuest Pro software (version 6.0; BD Biosciences). Each assay was conducted in triplicate.

Transduced cells were seeded in 12-well plates at a density of 500 cells/well and were incubated at 37°C in a 5% CO₂ atmosphere for 10 days. At the end of this incubation period, the cells were fixed using 4% paraformaldehyde for 15 min at room temperature and were subsequently stained with 0.1% crystal violet for 15-20 min at room temperature. Colonies were defined as clusters containing ≥50 cells and were observed using a light microscope (Leica Microsystems GmbH) and images were captured with a camera. The colonies were manually counted and the relative change in colony formation was determined. Data were normalized to the control group. Each assay was conducted in triplicate.

A total of 22 male BALB/c nude mice (age, 6-8 weeks; weight, 20-22 g) were acquired from Shanghai SLAC Laboratory Animal Co. Ltd. The mice were housed in a controlled, pathogen-free environment, with a stable temperature (22°C), humidity (40-70%) and under a 12-h light/dark cycle. The mice were provided unrestricted access to food and water. All experimental procedures involving animals were approved by the Institutional Animal Care and Use Committee of the Fujian University of Traditional Chinese Medicine (approval no. FJTCM IACUC 2019052) and were conducted in strict accordance with international guidelines (19) for the ethical use of animals in research.

Transduced HCT116, RKO and HT-29 cells (1×10⁶) in 100 μl medium mixed with 50% Matrigel were subcutaneously injected into the flank of nude mice (n=5, sh-FGGY and sh-Ctrl groups; or n=6, FGGY overexpression and control groups). Tumor growth was monitored every other day using a vernier caliper, with the volume calculated as: 1/2 (larger diameter x smaller diameter²). At the end of the experimental time point (16 days post-inoculation for shRNA groups and 20 days post-inoculation for overexpression groups), the mice were anesthetized using 2% isoflurane inhalation for induction and 1.5% isoflurane inhalation for maintenance. Tumor imaging was captured using the IVIS Spectrum live-animal imaging system (PerkinElmer, Inc.). Signal intensity was quantified based on the photon count per second from the region of interest. After imaging, the mice were sacrificed by cervical dislocation, and the tissues were harvested for further analysis.

iTRAQ analysis was used to screen differentially expressed proteins (DEPs) after knockdown of FGGY in HCT116 cells (20,21). The protein samples were reduced with 10 mM DTT at 56°C for 1 h and alkylated with 55 mM IAM in the dark for 1 h at room temperature. The protein samples were then digested overnight at 37°C with sequencing-grade modified trypsin (Promega Corporation) at a protein-to-trypsin ratio of 20:1. After trypsin digestion, peptides were dried by vacuum centrifugation (10,000 × g, 2 h, 45°C) and reconstituted in 0.5 M TEAB. The liquid chromatography (LC)-MS/MS analysis was performed by CapitalBio Technology Co., Ltd. using a Q Exactive mass spectrometer (Thermo Fisher Scientific, Inc.). The peptide separation was achieved by reversed-phase chromatography on a DIONEX nano-UPLC system, using an Acclaim C18 PepMap100 nano-Trap column (75 μm × 2 cm) connected to an Acclaim PepMap RSLC C18 analytical column (75 μm × 25 cm, 2 μm particle size) (all from Thermo Fisher Scientific, Inc.). Before loading, the sample was dissolved in mobile phase A, which consisted of 2% acetonitrile and 0.1% formic acid. A linear gradient of mobile phase B (0.1% formic acid in 99.9% acetonitrile) was applied, increasing from 2 to 35% over 45 min, followed by a sharp increase to 80% mobile phase B in 1 min, at a flow rate of 300 nl/min. The nano-LC was coupled online with the Q Exactive mass spectrometer using a stainless-steel emitter attached to a nanospray ion source. Nitrogen gas temperature was 0.6 Mpa and nebulizer pressure was 50 psi. MS analysis was performed in a data-dependent manner, with full scans (350-1,600 m/z) acquired using an Orbitrap mass analyzer (Thermo Fisher Scientific, Inc.) at a mass resolution of 70,000 at 400 m/z in Q Exactive. The 20 most intense precursor ions from each survey scan were selected for MS/MS analysis, which was performed at a mass resolution of 17,500 at 400 m/z in the Orbitrap analyzer. All tandem mass spectra were generated using the higher-energy collision dissociation method, with a dynamic exclusion of 20 sec. Proteins with >1.2-fold change in expression between groups and with a P-value of <0.05 were identified as DEPs. These DEPs were identified using volcano plots and further analyzed using hierarchical clustering plots in R software (version 4.1.0; R Foundation for Statistical Computing; https://www.R-project.org/) using the ggplot2 package (version 3.3.5; https://cran.r-project.org/web/packages/ggplot2/) for volcano plots and pheatmap package (version 1.0.12; https://cran.r-project.org/web/packages/pheatmap/) for hierarchical clustering analysis. To ascertain the biological pathways affected by these DEPs, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, Reactome pathway and PANTHER pathway analyses were performed (22-24). The enriched pathways were identified from the DEPs using KEGG database (https://www.genome.jp/kegg/) with a significance cutoff of P<0.05 (25,26).

Transduced cells were lysed using SDT lysis buffer (Beyotime Institute of Biotechnology). Proteins were then digested with trypsin using filter-aided sample preparation. A total of 100 μg peptide mixture from each sample was labeled using iTRAQ reagent (SCIEX) adhering to the manufacturer's guidelines. Phosphopeptide enrichment was carried out, followed by LC-MS/MS analysis as per the manufacturer's instruction (L-3000 HPLC system; RIGOL Technologies) (27,28). Subsequently, identification and quantification of phosphorylated proteins were achieved using the MASCOT engine (version 2.2; Matrix Science), integrated into Proteome Discoverer 2.4.

Motif analysis was carried out using MeMe (http://meme-suite.org/). The amino acid sequences containing the modified site, and six upstream/downstream amino acids from the modified site (resulting in a total of 13 amino acid positions), were extracted and used for motif prediction (parameters used: Width, 13; occurrences, 20; background, species). Phosphopeptide sequences were analyzed with InterProScan software (version 5.52-86.0; EMBL-EBI)to identify protein domain signatures from the InterPro member database Pfam (22,23). Phosphopeptides showing ≥1.2-fold difference in expression in comparative analyses and with a P-value of <0.05 were considered differentially expressed. These phosphopeptides were visualized using volcano plots and further analyzed using hierarchical clustering plots. Proteins were categorized based on their KEGG annotations with P<0.05.

Cells with knockdown of FGGY or overexpression of FGGY were first washed with 1xPBS, and subsequently fixed with 4% paraformaldehyde for 10-15 min at room temperature. After fixation, the cells were washed and then incubated with β-gal staining solution (cat no. ab102534; Abcam) in a dry incubator at 37°C overnight. Images of the β-gal stained-cells were captured using a light microscope (Leica Microsystems GmbH) at a ×200 magnification. Each assay was conducted in triplicate.

Cells were initially rinsed with PBS after knockdown of FGGY, and were then fixed in 4% paraformaldehyde for 20 min at room temperature. The cells were subsequently permeabilized with 0.1% Triton X-100, blocked with 1% BSA (cat. no. SW3015; Beijing Solarbio Science & Technology Co., Ltd.) for 1 h at room temperature, and were incubated overnight at 4°C with primary antibodies against HP1γ (cat. no. 2619; Cell Signaling Technology, Inc.), or trimethylation of H3K9 (H3k9me3; cat. no. 49-1008; Thermo Fisher Scientific, Inc.), at a dilution of 1:200. Post-primary antibody incubation, the cells were washed with PBS and incubated with Goat anti-Rabbit IgG Alexa Fluor™ 647 (cat. no. A-21245; Thermo Fisher Scientific, Inc.) at a dilution of 1:200 for 1 h at room temperature, followed by Hoechst 33342 staining for 10 min at room temperature. All images were captured using a fluorescence microscope (Leica Microsystems GmbH) at a magnification of ×400. Each assay was conducted in triplicate.

Experiments were performed at least in triplicate and data are presented as the mean ± standard deviation. Statistical analysis was conducted using SPSS 26.0 software (IBM Corp). Comparisons between two datasets were made using a two-tailed unpaired Student's t-test for normally distributed data or Mann Whitney-U test for non-parametric data. Multiple group comparisons were analyzed using one-way ANOVA followed by Tukey's post hoc test or two-way ANOVA followed by Bonferroni's post hoc test. The associations between FGGY protein expression and clinicopathological characteristics were analyzed using χ² test or Fisher's exact test (when the expected counts were ≤5 in >20% of samples) for categorical variables. The correlation between genes was analyzed using Pearson correlation coefficient. P<0.05 was considered to indicate a statistically significant difference.

The mRNA and protein expression levels of FGGY in various types of digestive system neoplasms were determined by GEPIA database analysis and TMA analysis in COAD, READ, PAAD, ESCA, STAD and LIHC. There was an elevation in FGGY mRNA and protein expression in PAAD, COAD and READ tissues compared with in normal tissues from healthy controls (Fig. 1A-F). This increase was also observed in COAD and READ tissues compared with in adjacent noncancerous tissues (Fig. 1A-D). However, FGGY protein expression did not exhibit significant differences in ESCA, STAD and LIHC tissues compared with in normal tissues from healthy controls and adjacent noncancerous tissues (Fig. 1G-L). These findings suggested that FGGY could serve as a potential specific marker of CRC, with increased expression in CRC.

Furthermore, the analysis of TCGA and GEO datasets (GSE20916 and GSE13513) further confirmed the upregulation of FGGY mRNA expression in CRC tissues compared with in adjacent normal tissues (Fig. 2A-C). This finding was further verified at the mRNA level using a cDNA array by RT-qPCR analysis (Fig. 2D; n=15 pairs) and at the protein level using a TMA via IHC analysis (Fig. 2E; n=71 pairs) in CRC tissues compared with in adjacent normal tissues. Moreover, survival analysis revealed an association between high FGGY protein expression and shorter overall survival in patients with CRC (Fig. 2F; high expression, n=71; low expression, n=23; cutoff, 3). Analysis of the association of FGGY with clinicopathological characteristics revealed that higher FGGY expression was associated with N stage in patients with CRC (Table I).

Given that the endogenous expression levels of FGGY were relatively high in HCT116 and RKO cells, but low in HT-29 cells (Fig. S1A), the present study chose to knock down FGGY in HCT116 and RKO cells, and to overexpress FGGY in HT-29 cells. To investigate the role of FGGY in CRC cell viability, three distinct shRNAs specifically targeting FGGY were utilized, which were encoded into lentiviruses for cell transduction. Compared with in the sh-Ctrl group, transducing HCT116 cells with these shRNAs significantly reduced FGGY expression at both the mRNA and protein levels (Fig. S1B and C), and decreased cell viability (Fig. S1D). Since all three shRNAs exhibited a comparable effect on FGGY expression and cell viability, further experiments were conducted with a single construct, designated as sh-FGGY (specifically, sh-FGGY-1). RT-qPCR and western blot analysis confirmed a substantial decrease in FGGY expression at both the mRNA and protein levels in HCT116 and RKO cells after transduction with sh-FGGY (Fig. 3A and B). Additionally, FGGY knockdown significantly reduced the number of colonies formed and cell viability (Fig. 3C and D). Furthermore, FGGY knockdown resulted in an increased proportion of cells at the G₀/G₁ phase and a decreased proportion of cells at the S phase in both HCT116 and RKO cells (Fig. 3E). Moreover, Annexin V staining and subsequent flow cytometric analysis revealed that FGGY knockdown significantly elevated the percentage of early apoptotic cells (Fig. 3F). Additionally, FGGY knockdown increased the expression levels of the pro-apoptotic protein Bax and decreased the expression levels of anti-apoptotic Bcl-2 in HCT116 cells (Fig. S2). By contrast, HT-29 cells transduced with lentiviral vectors designed to overexpress FGGY exhibited elevated FGGY protein expression and increased cell viability (Fig. S3A and B).

Analysis of tumor growth in vivo using a xenograft mouse model revealed that, compared with in the sh-Ctrl group, FGGY knockdown significantly inhibited CRC tumor growth in nude mice, as indicated by a marked reduction in tumor volume (Fig. 4A), fluorescence intensity of tumor cells (Fig. 4B) and tumor weight (Fig. 4C). By contrast, overexpression of FGGY significantly promoted CRC tumor growth in nude mice, as evidenced by an increase in tumor volume (Fig. S3C), fluorescence intensity of tumor cells (Fig. S3D) and tumor weight (Fig. S3E) in nude mice.

To further explore the possible mechanisms by which FGGY knockdown inhibits tumor growth in CRC, iTRAQ and phosphoproteomics analyses were employed to identify DEPs in HCT116 cells after FGGY knockdown. As shown in Fig. 5, a total of 809 DEPs were identified, which included 236 proteins with increased expression and 573 with decreased expression (fold change >1.2, P<0.05) using iTRAQ methodology (Fig. 5A and B; Tables SV and SVI). Additionally, phosphoproteomics analysis identified 304 differentially expressed phosphopeptides, comprising 145 upregulated and 159 downregulated phosphopeptides (fold change >1.2, P<0.05) (Fig. 5D and E; Tables SVII and SVIII). The KEGG pathway enrichment analysis revealed that these proteins and phosphopeptides were enriched in various signaling pathways (Fig. 5C and F; Tables SIX and SX). KEGG analysis of iTRAQ revealed enrichment in both the 'Formation of senescence-associated heterochromatin foci (SAHF)' and the 'Senescence-associated secretory phenotype (SASP)' (Fig. 5C). These findings suggested that cellular senescence may serve a significant role in the effects of FGGY knockdown. Additionally, KEGG analysis of phosphoproteomics analysis revealed enrichment in both the 'p53 pathway' and 'p53 signaling pathway' (Fig. 5F). Given the well-established link between cellular senescence and the p53 pathway, the present study further investigated these two processes. Both cellular senescence and the p53 pathway are known to regulate cell death and tumor suppression, indicating that FGGY may modulate the p53 signaling pathway and influence cell senescence.

The present study aimed to elucidate how FGGY knockdown regulates cell senescence and the p53 pathway in CRC. Correlation analysis of microarray expression data from the GEO and TCGA databases revealed that a weak negative correlation existed between FGGY and p53 expression (Fig. 6A and B). Compared with in the control group, FGGY knockdown led to an increase in the protein expression levels of p53 and p21, while the expression levels of PCNA were decreased (Fig. 6C). Notably, when comparing the two cell lines, the inhibitory effects of FGGY knockdown on the viability of HCT116/p53^+/+ cells were attenuated in HCT116/p53^−/− cells (Fig. 6D and E). By contrast, HT-29 cells transduced with lentiviral vectors designed to overexpress FGGY exhibited decreased p53 and p21 protein expression compared with that in cells transduced with control vectors (Fig. S3F).

Cellular senescence is characterized by increased SA-β-gal activity, reduced cell proliferation, cell cycle arrest and the formation of SAHF (29-31). To further investigate the potential regulatory effects of FGGY on senescence, a SA-β-gal activity assay was conducted in both HCT116/p53^+/+ and HCT116/p53^−/− cells. The results revealed that β-gal activity was significantly increased after FGGY knockdown in HCT116/p53^+/+ cells, but β-gal activity remained unchanged after FGGY knockdown in HCT116/p53^−/− cells (Fig. 6F). Conversely, β-gal activity was decreased in HT-29 cells after FGGY overexpression (Fig. S3G). Moreover, FGGY knockdown significantly upregulated the expression of cell senescence markers HP1γ and H3k9me3m in HCT116/p53^+/+ cells, whereas in HCT116/p53^−/− cells, FGGY knockdown did not affect the protein expression of cell senescence markers HP1γ and H3K9me3 (Fig. 6G and H).