Patient enrollment
Cervical squamous cell carcinoma tissue, cervical adenocarcinoma tissue, and normal cervical control tissue (patients with hysterectomy due to uterine myoma and no cervical disease confirmed by pathology) specimens were collected by gynecological surgical resection at Xuzhou Central Hospital from 2021 to 2023. There were 44 cases of squamous cell carcinoma, 20 cases of adenocarcinoma, and 60 cases of normal cervical tissue. Additionally, Between January 2016 and July 2018, the Department of Gynecology at Xuzhou Central Hospital collected samples and data from 64 patients who were diagnosed with cervical cancer and had radical surgery. Fresh tissue was frozen in liquid nitrogen immediately after removal, formalin fixed and embedded. All samples used in this study were confirmed pathologically. All patients provided written informed consent prior to sample collection, and the Xuzhou Central Hospital's Ethical Committee approved the use of the clinical samples (XZXY-LJ-20210513-052).
Label-free relative quantification proteomics
Label-free MS proteomics quantitation (LC-MS/MS) was performed on three groups of samples: cervical adenocarcinoma tissues, cervical squamous carcinoma tissues and normal control tissues, to detect and analyze the significantly differentially expressed proteins in each group of samples. Maxquant Perseus software was used to analyze the differential expression of omics data. The data were imported to retrieve the output proteinGroups.txt file, and the corresponding LFQ intensity was selected. The data with Reverse, only identified by site and < s:1 > were removed, the LFQ intensity value was taken as the log base of 2, the disease and control groups were classified, and the effective data were screened. Significant differentially expressed genes (up 2-fold or down 0.5-fold Foldchange (FC), p < 0.05) were screened by t-test between the two groups. Proteins with pp-regulated or both down-regulated expression in adenocarcinoma (Adeno) and squamous carcinoma (Squa) tissues were analyzed using the InteractiVenn online analysis tool and presented in Venn diagrams. KEGG functional enrichment analysis was performed to identify differentially expressed proteins. The analyses of up-regulation and down-regulation were conducted separately and together. KEGG pathway enrichment analysis of differentially expressed genes was performed using DAVID 6.8 (or KOBAS) tool, and p < 0.05 was selected as the significant threshold. KEGG mapper was used for pathway map analysis of differentially expressed genes.
Cell culture
The Hela, SiHa, and C33A human cervical cancer cell lines were acquired from Shanghai Gaining Biological Technology Co., Ltd. (Shanghai, China). All cell lines underwent Short Tandem Repeat (STR) profiling authentication. The cells were cultivated in Dulbecco's Modified Eagle's Medium (GIBCO, Billings, MT, United States) supplemented with 10% fetal bovine serum (GIBCO) and 100 IU/mL penicillin-streptomycin combination (Sigma-Aldrich). Cultures of transfected cells were maintained in high-glucose Dulbecco's modified Eagle's medium (GIBCO, USA) supplemented with 10% FBS, 100 U/mL penicillin sodium, and 100 µg/mL streptomycin sulphate (Hyclone, USA). The cells were cultured in a humidified incubator with 5% CO2 at 37°C.
Immunohistochemistry (IHC)
RCN1 was investigated by IHC in cervical cancer tissues that were fixed in formalin and paraformaldehyde. PowerVisionTM Two-Step Histostaining Reagent (Zhongshan Golden Bridge Biotechnology, Beijing, China) was used for IHC staining. Cytoplasmic staining was independently evaluated by two pathologists in a double-blind fashion. Scores were assigned to the cytoplasm staining intensity (SI) as follows: 0 for negative, 1 for weak, 2 for moderate, or 3 for strong staining. Additionally, scores were assigned for the proportion (PP) of stained cells as follows: 0 (less than 5%), 1 (6–25%), 2 (26–50%), 3 (51–75%), and 4 (more than 75%). Finally, the histoscore (Q) was determined using the following formula: Q = SI×PP. For pathology evaluation, an immunoreactivity score (IRS) of 1–6 was considered low expression of RCN1, while IRS 8–12 was considered strong expression of RCN1.
Survival analysis
Cervical cancer samples were separated into high and low RCN1 groups using the optimal cutoff in the R survMisc package. Kaplan-Meier analysis was performed using log-rank test to investigate the relationship between the expression of RCN1 and survival.
Real-time quantitative PCR (RT-qPCR)
Total RNA was extracted from cervical cancer tissues or cell lines and reverse transcribed into cDNA. Then, RCN1 mRNA expression was assessed by RT-qPCR, as previously described. The primer sequences were as follows:
RCN1 FORWARD primer GGGAGGAGTTCACTGCCTTTCTG
RCN1 REVERSE primer TCATCCTGATCCACAAACCCATCC
GAPDH FORWARD primer CAGGAGGCATTGCTGATGAT
GAPDH REVERSE primer GAAGGCTGGGGCTCATTT
Western blot assay
As directed by the manufacturer, fresh tissues and cell extracts were generated using cell lysis reagent (Sigma-Aldrich, St. Louis, MO, USA), and the protein was measured using a BCA assay (P0010, Biyuntian Institute of Biotechnology, China). Equivalent amounts of protein were separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and then detected using a rabbit monoclonal anti-RCN1 antibody (P0023D) from the Biyuntian Institute of Biotechnology in China. GAPDH was evaluated as a loading control. Experiments were performed in triplicate.
Vectors construction and transfection
GeneChem Corporation (Shanghai, China) designed and constructed lentiviral particles harboring shRNA targeting RCN1 (sh-RCN1) and negative control shRNA (sh-NC). The shRNA targeting sequences are provided in Supplementary Table S1. One day before transfection, the C33A or Hela cells were plated in six-well plates. The cells were transfected with the lentiviruses for 12 hours, and then they were selected in puromycin for seven days. Western blotting and PCR were used to validate RCN-1 overexpression and knockdown.
Cell Counting Kit 8 (CCK-8) assay
To determine cell viability, the Cell Counting Kit 8 (TAOSHU, Japan) was utilized. Rapidly proliferating cells were seeded onto 96-well plates, and the optical density at 450 nm was measured using the Gen5 system (BioTek, Winooski, VT, USA) once every 24 hours for five days.
Scratch tests
Control and RCN1 overexpression cells were seeded at the same concentration (5×105 cells/well) in 6-well cell culture plates. After 24 h, a straight line was made perpendicularly in each well using a 200 µL lance tip. Serum-free medium was then added, and the plates were cultured an incubator with 5% CO2 at 37°C. The cell migration distance was measured under a light microscope at 0 and 24 hours.
Transwell assays
Transwell invasion assays were conducted using Transwell inserts (Corning, MA, USA) in 24-well plates, either with or without pre-coated Matrigel. Five hundred and forty-eight cells per well were seeded in the top chamber in 200 µL of 2% FBS and 500 µL of 20%. The medium in the lower chamber was supplemented with FBS. After a full day, the lower surface was fixed with methanol and stained with hematoxylin and eosin (H&E), while the upper surface was carefully removed using a cotton swab. Four fields of vision that were chosen at random were measured under a microscope.
IP-MS
C33A sh-RNA cell lines were detected and quantitatively analyzed by mass spectrometry, and the FLQ Intensity (relative quantitative signal values) was evaluated by T-testing. The FLQ Intensity Ratios in samples with p value < 0.05 were used as screening criteria to identify and screen highly possible interacting proteins.
Beads were washed three times by adding 500 µL TBS, vortexing, centrifuging for 30 s at 5000 g and removing the supernatant. Beads were digested in 100 µL of Buffer 1 (containing 2 M urea, 50 mM Tris-HCl pH 7.5 and 5 µg/mL Trypsin) for 30 m 27°C in a thermomixer, shaking at 800 rpm. After the initial digestion, the samples were centrifuged for 30 s at 7000 rpm, and the supernatant was removed and collected into fresh buffer. The beads were washed twice in 50 µL of buffer 2 (containing 2 M urea, 50 mM Tris-HCl pH 7.5 and 1 mM DTT), and the supernatants were pooled. The samples were left on the bench to continue digesting overnight at room temperature. Then, the samples were treated with 1 µL trifluoroacetic acid (TFA) to stop the digestion and were desalted in C18 stagetips. Tryptic peptides were desalted and centrifuged in a speedvac to dry. Then, Tryptic peptides were dissolved in 0.1% FA.
Liquid chromatography–tandem mass spectrometry (LC–MS/MS)
For LC-MS/MS analysis, the peptides were separated by a 60 min gradient elution at a flow rate of 0.22 µl/min with a Thermo Scientific EASY-nLC 1000 HPLC system, which was directly interfaced with a Thermo Scientific Orbitrap Exploris 480 mass spectrometer that was operated in the data-dependent acquisition mode using Xcalibur 3.0 software. This included a single full-scan mass spectrum in the orbitrap (300–2000 m/z, 70,000 resolution) followed by 20 data-dependent MS/MS scans at 27% normalized collision energy (HCD).
MS/MS spectra from raw data were imported to Proteome Discoverer (version 3.0, Thermo Scientific), and a database search was performed using MS CHIMERYS™ with a combination of human proteins (uniprot reference, version 2022-02-03, 20,509 entries) as well as common contaminants (PD_Contaminants, 344 entries). The intensity from PD was further processed in the Perseus computational platform (v 1.6.7.0). LFQ intensity values were log2 transformed, and the samples were grouped into experimental categories. Proteins not identified in 3 out of 3 replicates in at least one group were also removed. Missing values were imputed using normally distributed values with a 1.8 downshift (log2) and randomized 0.3 width (log2) for whole matrix values. Statistical analysis was performed to determine which proteins were significantly enriched in the C33A cell compared to the C93 cell samples or isotype IgG control sample (t-test with permutation-based False Discovery Rate (FDR) = 0.05 and S0 = 0.1).
In vivo assays
Two sets of BALB/c-nu mice were utilized to monitor tumor growth: oe-C33A-NC [mice implanted with control C33A cells (n = 8) that were infected with scrambled oeRNA]; and oe-C33A-RCN1 [mice inoculated with RCN1 oeRNA-infected C33A cells (n = 8)]. The mice were subcutaneously injected with 5×106 cells/100 µl. Every three days, the longest and shortest diameters of the developing tumors were measured using calipers, and the tumor volume (V) was calculated using the formula: V = (the longest diameter the shortest diameter2)/2. Cervical cancer metastasis was determined using a model of pulmonary metastasis. Live imaging was conducted once every 3 days, and the mice were euthanized 1 month after injection.
Statistical analysis
Categorical variables were quantified with n (%), whereas continuous variables were calculated using the meanSD of at least three independent experiments. Data analysis was done with SPSS 24.0 or Rstudio. Wilcoxon rank-sum tests were applied for two groups. Statistical significance was determined using the Student's two-tailed t-test, a parametric generalized linear model with random effects for the growth curve, and one-way ANOVA with LSD-t test for multiple groups. The Spearman's rank correlation test and the chi-square test were used in correlation analyses. In each two-sided statistical test, p value < 0.05 was considered statistically significant.