Abstract

Aim. The aim of this study is to demonstrate the expression and clinicopathological significance of complement C1q B chain (C1QB) in cervical cancer. Methods. In total, 120 cervical cancer tissues, as well as 20 samples each of high-grade squamous intraepithelial lesions (HSILs), low-grade squamous intraepithelial lesions (LSILs), and benign cervical tissue, were collected to evaluate the expression of C1QB protein via immunohistochemical staining. We conducted an integrated analysis of C1QB mRNA expression in cervical cancer using public microarrays and RNA-seq data sets by calculating standard mean differences (SMDs). Simultaneously, we explored the relations of C1QB with clinicopathological parameters and the expression of P16, Ki-67, and P53. Results. The expression of C1QB protein was higher in cervical cancer samples than that in benign cervical tissue, LSIL, and HSIL samples (). A combined SMD of 0.65 (95% CI: [0.52, 0.79], ) revealed upregulation of C1QB mRNA in cervical cancer. C1QB expression may also be related to the depth of infiltration, lymphovascular invasion, and perineural invasion in cervical cancer (). We also found that C1QB protein expression was positively correlated with P16 and Ki-67 expression in cervical cancer (). The gene set enrichment analysis showed that C1QB may participate in apoptosis and autophagy. A relationship was predicted between C1QB expression and drug sensitivity to cisplatin, paclitaxel, and docetaxel. Conclusion. We confirmed the overexpression of C1QB in cervical cancer at both mRNA and protein levels for the first time. C1QB may serve as an oncogene in the tumorigenesis of cervical cancer, but this possibility requires further study.

1. Introduction

Cervical cancer originates from embryonic paramesonephric ducts and currently ranks as the most frequent cancer of the female reproductive tract. Some estimates from 2020 indicate that approximately 604,000 cases of cervical cancer occurred and caused 342,000 deaths [1]. Human papillomavirus (HPV) infection, especially high-risk HPV subtype 16/18 infection, plays a crucial role in the tumorigenesis of cervical cancer [24], and although HPV vaccination can prevent cervical intraepithelial neoplasia and cervical cancer to some extent, some patients still suffer from cervical cancer unrelated to HPV.

The tumorigenesis and development of cervical cancer are complex, and it is affected by multiple factors [1, 46]. Recently, single nucleotide polymorphisms and post-transcriptional regulation have been associated with tumorigenesis in cervical cancer [7, 8]. In this context, the expression of complement C1q B chain (C1QB), located in 1p36.12, may be relevant. Complement C1QB encodes the B-chain of serum complement protein C1QB and acts as the first recognition subunit in the classical complement pathway. C1QB expression has been reported in the tumor microenvironment of multiple cancers, such as ovarian cancer, prostate cancer, glioma, and osteosarcoma [912]. Overexpression of the C1QB gene at the mRNA level has also been shown in gastric cancer [13]. However, dysregulation of C1QB has not yet been identified in cervical cancer.

The aim of the present study was to use immunohistochemical (IHC) staining to examine the expression status of the C1QB protein in cervical cancer tissues. We also explored the possible correlation between C1QB expression and the expression of Ki-67, P16, and P53. In addition, we used multicenter high-throughput data sets to determine the expression level of C1QB in cervical cancer tissues at the mRNA level. The overall goal was to clarify the molecular mechanism of cervical cancer and to identify potential therapeutic targets.

2. Materials and Methods

2.1. Collection of Experimental Samples

We collected 120 cervical cancer tissues, as well as 20 samples each of high-grade squamous intraepithelial lesions (HSILs), low-grade squamous intraepithelial lesions (LSILs), and benign cervical epithelial tissues, from Guangxi Medical University Cancer Hospital. All the clinical cases were patients who had undergone surgery in our hospital from January 1, 2018, to October 31, 2021. This study was approved by the Ethics Committee of Guangxi Medical University Cancer Hospital (No. 2022007), and all patients provided signed informed consent. This study conformed to the standards set by the Declaration of Helsinki.

2.2. IHC Staining

Formalin-fixed and paraffin-embedded tissue slides were deparaffinized with xylene and rehydrated with ethanol, followed by antigen retrieval with ethylenediaminetetraacetic acid buffer (). The IHC staining was performed with an HRP-Polymer anti-Mouse/Rabbit IHC Kit (MaxVision™) and anti-C1QB polycolonal antibody (Abcam, EPR2981, dilution 1 : 70), following the manufacturers’ instructions. Two pathologists evaluated the results of the staining by microscopy. The positive cells in the visual field were scored with the following criteria: 0–5% (0 points), 6–25% (1 point), 26–50% (2 points), 51–75% (3 points), and >75% (4 points). The staining intensity was scored as follows: no staining (0 points), weak staining (1 point), moderate staining (2 points), and strong staining (3 points). We calculated the final IHC staining score for C1QB by multiplying the intensity score and positive cell score. We then divided the expression level of C1QB into three levels: (−), (+), and (++).

We also assessed the expression levels of P16, Ki-67, and P53 by IHC staining. The IHC staining results for P16 (Maxim, mouse-anti-human monoclonal antibody, MAB-0673) were rated as follows: continuous diffuse tan staining appearing in the cell nucleus and/or cytoplasm was marked as positive staining, whereas no staining, sporadic staining, or focal staining was marked as negative staining. Cells showing tan staining nuclei after IHC for Ki-67 (Maxim, rabbit-anti-human monoclonal antibody, RMA-0731) were rated as positive cells. A proportion of positive cells ≥10% was deemed positive staining for Ki-67; otherwise, the results were deemed negative staining. For P53 staining (Maxim, mouse-anti-human monoclonal antibody, MAB-0674), a proportion of positive cells >0% and <80% was defined as wild-type P53 (P53WT); otherwise, the staining result was defined as mutant-type P53 (P53MT). Simultaneously, hosphate buffer saline (PBS) was utilized to replace primary antibodies as negative controls.

2.3. Retrieval of High-Throughput Data Sets Related to Cervical Cancer

We searched The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium databases for tertiary RNA sequencing (RNA-seq) data sets of cervical cancer. Normal uterine cervical samples were acquired from the Genotype-Tissue Expression (GTEx) database. The sample size was further expanded using the Sequence Read Archive, Oncomine, ArrayExpress, and Gene Expression Omnibus databases to retrieve microarrays of cervical cancer. The following criteria were used for inclusion of a data set: (1) samples were from cervical cancer patients or cervical cancer cell lines; (2) samples and clinical cases did not receive intervention by drugs, radiation, or siRNA; (3) both cervical cancer samples and benign controls included more than 3 cases; (4) the probe matched with an official gene symbol; and (5) the expression profiles contained C1QB. Figure 1 demonstrates the technological process of retrieving datasets. As of July 31, 2021, 19 high-throughput data sets were obtained for our study. We combined the microarrays from the same platform and removed batch effects using the sva package of R [14].

2.4. Gene Set Enrichment Analysis

The underlying mechanisms of C1QB in cervical cancer were explored by gene set enrichment analysis (GSEA) with the TCGA-CESC cohort using the GSEA software 4.0.8. The reference gene set files “c2.cp.kegg.v7.2.symbols.gmt” and “c2.cp.reactome.v7.2.symbols.gmt” were downloaded from MSigDB. The software was used to calculate the normalized enrichment score (NSE), -value, and adjusted (FDR-q value).

2.5. Prediction of the Relation between C1QB Expression and Drug Sensitivity

We predicted the half-maximal inhibitory concentration (IC50) of every drug in the Genomics of Drug Sensitivity in Cancer database for each sample in the TCGA-CESC cohort based on ridge regression using the pRRophetic package in R software [15]. The results were visualized as box plots.

2.6. Statistical Analysis

Enumeration data were shown as proportions. The differences were analyzed with the χ2 test using the SPSS software 24.0. Spearman’s rank correlation was used to explore the correlation between C1QB expression and P16, Ki-67, and P53 expression. Student’s t-test was used to compare the expression levels of C1QB mRNA between cervical cancer and non-tumor control tissues using the GraphPad Prism 8 software. An integrated study was then conducted to calculate the standard mean difference (SMD) and to examine heterogeneity using the Stata v15.1 software (College Station, TX, USA). The integrated SMD was visualized by a forest plot, and the summary receiver operating characteristic (sROC) curve was drawn. The ROC and sROC curves were used to assess the ability of C1QB to distinguish cervical cancer from normal uterine cervical samples. The area under the curve (AUC) ≥0.7 displayed a moderate discriminatory capacity. Wilcoxon’s test was used to compare the differences in IC50 between the high C1QB and low C1QB expression groups. A two-tailed indicates a statistically significant difference.

3. Results

3.1. The Positive Rate of C1QB in Different Groups

The positive rate of C1QB protein expression in the cervical cancer group was 75.8%, which was higher than that in the LSIL group and in the benign cervical tissue group (; Table 1). The positive rate of C1QB protein expression was higher in the HSIL group than that in the LSIL and benign cervical tissue groups (; Table 1). Figure 2 shows the IHC staining of C1QB protein in benign cervical tissues, LSIL, HSIL, and cervical cancer tissues.

3.2. Relationships between C1QB Protein Expression and Clinicopathological Parameters

Cervical cancer with deeper infiltration (depth ≥1/2 muscle layer, ), lymphovascular invasion (), and perineural invasion () tended to have positive C1QB expression, and the difference was statistically significant (Table 2). The expression of C1QB in adenocarcinoma, squamous cell carcinoma, and adenosquamous carcinoma was also statistically different (). By contrast, the differences in C1QB expression in persons of different nationalities or ages or in tumors of different sizes, clinical stages, lymph node metastasis, gross appearance, and HPV infection were not statistically significant ().

3.3. Correlations between the Expression of C1QB and the Expression of P16, Ki-67, and P53

IHC staining revealed a positive correlation between C1QB expression and the expression of Ki-67 (Spearman’s , ; Table 3) and P16 (Spearman’s , ; Table 3). However, no statistically significant correlation was found between C1QB and P53 expression (; Table 3). Figure 3 shows the IHC staining patterns for C1QB, P16, Ki-67, and P53 in cervical cancer tissues.

3.4. Microarray and RNA-seq Validation of C1QB mRNA Expression in Cervical Cancer

As of December 1, 2021, 18 external microarrays and one RNA-seq data set were collected (Table 4). Integration of the microarrays from the same platforms left us with ten cohorts. In seven data sets, the expression of C1QB mRNA was significantly higher in cervical cancer tissues than that in non-tumor cervical tissues, and the difference was statistically significant (TCGA_GTEx, ; GPL96, ; GPL570, ; GPL570, ; GPL6244, ; GSE138080-GPL4133, ; GSE39001-GPL201, ; Figures 4(a)–4(g)). However, the difference was not significant in the three other data sets (GSE55940-GPL16238, ; GSE46857-GPL7025, ; GSE7410-GPL1708, ; Figures 4(h)–4(j)). A comprehensive calculation of the SMD verified the upregulation of C1QB in cervical cancer tissues (SMD = 0.65, 95% CI [0.52, 0.79], ; Figure 5(a)), and the publication bias was not statistically significant (; Figure 5(b)).

Figures 6(a) and 6(b) illustrate the ROC curves for the ten datasets. Among the nine ROC curves, the largest AUC was 0.8430 (; Figure 6(c)). The AUC of the sROC was 0.79, indicating a moderate discriminatory ability of C1QB in cervical cancer (Figure 6(d)).

3.5. Gene Set Enrichment Analysis

Through GSEA, we identified that the group with high C1QB expression was mainly enriched in some immune-related gene sets, such as “KEGG_B_CELL_RECPTOR_SIGNALING_PATHWAY” (NSE = 2.764, , FDR-), “KEGG_T_CELL_RECPTOR_SIGNALING_PATHWAY” (NSE = 2.707, , FDR-), and “REACTOME_NEUTROPHIL_DEGRANULATION” (NSE = 2.980, , FDR-). Other gene sets, such as “KEGG_APOPTOSIS” (NSE = 2.543, , FDR-) and “REACTOME_TOLL_LIKE_RECEPTOR_CASCADES” (NSE = 2.927, , FDR-), were unrelated to immune function. The group with low C1QB expression was mainly enriched in metabolism-related gene sets (Figures 7(a) and 7(b)).

3.6. Relationship between C1QB Expression and Cervical Cancer Drug Sensitivity

Our prediction of the IC50 with ridge regression revealed that the group with high C1QB expression tended to have a high IC50 for chemotherapy drugs, such as cisplatin, paclitaxel, docetaxel, nilotinib, erlotinib, and mitomycin C (; Figures 8(a)–8(f)). This indicated that patients with high C1QB expression might have tumors with low drug sensitivity.

4. Discussion

This study confirmed the high expression of C1QB protein in clinicopathological specimens of cervical cancer. However, it also added the novel finding that this upregulation was also seen at the mRNA level in 1,191 samples. Our investigation of the relationships between C1QB protein expression and the clinicopathological parameters of cervical cancer patients also provided the first correlation between C1QB expression and the expression of P16, Ki-67, and P53. In addition, we identified the underlying function of C1QB in cervical cancer using GSEA.

A variable expression status of C1QB has been reported in some malignant tumors. For instance, upregulation of C1QB at the mRNA level was reported in gastric cancer and head and neck squamous cell carcinoma, but C1QB mRNA was downregulated in esophageal squamous cell carcinoma [13, 16]. The current study is the first to report a higher expression of C1QB protein in cervical cancer tissues than in non-tumor cervical tissues (based on the PubMed database, as of December 28, 2021). Another novel finding was that the positive rates of C1QB progressively increased from the benign cervical tissue (10%) to LSIL (15%), then to HSIL (65%), and finally to cervical cancer (75.8%). We also confirmed the upregulation of C1QB in cervical cancer at the mRNA level through a comprehensive analysis based on a large sample size (n of cervical cancer = 825 and n of non-tumor = 366). The collection of multicenter samples worldwide from China (n = 10), USA (n = 664), Mexico (n = 303), Netherlands (n = 80), UK (n = 77), India (n = 29), and Germany (n = 28) helped to reduce the influence of ethnicity on the results.

The dysregulation of C1QB in many malignant tumors makes it an attractive target for clinical research. One study has reported that an osteosarcoma patient with high C1QB expression tended to have a favorable outcome and that C1QB expression was correlated with the percent necrosis observed at definitive surgery [17]. Another study found that C1QB might be a protective factor in osteosarcoma patients [12]. However, a different study revealed that gastric cancer patients with high C1QB expression had a poor prognosis [13]. Similarly, a glioma-related study found that expression of C1QB mRNA was negatively related to the survival rate in patients with grade III glioma and glioblastoma [11]. In the present study, we found that positive C1QB expression was related to deeper invasion, lymphovascular invasion, and perineural invasion of cervical cancer, providing the first evidence for clinicopathological implications of C1QB expression in cervical cancer.

The combination of P16/Ki-67 detection has been widely applied in the auxiliary diagnosis of cervical cancer and cervical intraepithelial neoplasia [1822]. A prospective cohort study revealed that patients with negative dual staining for P16/Ki-67 had a low risk of precancerous lesions and cervical cancer [23]. Ki-67 has been recognized as a proliferation marker of malignant tumor cells [24, 25]. In our study, we demonstrated a positive correlation between C1QB expression and the expression of Ki-67 and P16 in cervical cancer, indicating that C1QB may act as a facilitating factor in cervical cancer tumorigenesis. However, this possibility still needs further validation.

The C1QB protein can bind to IgG molecules and activate the classical pathway of complement, and this is the acknowledged function of the C1QB complement [26]. However, in addition to its classical function, C1QB was demonstrated to accelerate primary hemostasis through an interaction with von Willebrand factor [27]. We explored the latent mechanisms of C1QB in cervical cancer by conducting GSEA, and we determined that C1QB may have a function related to apoptosis. Apoptosis is a programmed cell death that occurs regularly in a sequential order and ensures the balance between cell proliferation and cell death [28, 29]. Resistance to apoptosis is considered a fundamental capacity of tumor cells and plays a significant role in the tumorigenesis of numerous malignancies [2931]. In recent years, numerous studies have indicated that apoptosis and related proteins and factors, including the Bcl-2 superfamily, p53, and PI3K–Akt signaling, may be crucial therapeutic targets in the treatment of cancers and may also trigger multi-drug resistance of some cancers to chemotherapy [29, 30, 32]. One study reported that C1QB could upregulate the expression of Fas and TNF-α and induce apoptosis of ovarian cancer cells, either independently or through a caspase cascade reaction [9].

Another form of cell death is autophagy, which is conducted by lysosomes [28, 33]. The fundamental process of autophagy involves the transport of macroproteins or whole organelles into the lysosome and subsequent digestion by lysosomal enzymes. Deficient regulation of autophagy has been linked to tumor malignancy [34]. Some studies have revealed that HPV infection, a known cause of cervical cancer, may induce the dysregulation of p62 and Beclin-1 and inhibit autophagy, thereby promoting cervical cancer tumorigenesis [35]. Thus, autophagy remains an attractive area of research, and it may represent another potential therapeutic target for cancer treatment [3640]. In our study, we demonstrated that C1QB expression in cervical cancer might be related to both apoptosis and autophagy, suggesting the possible participation of C1QB in cervical cancer tumorigenesis. However, this possibility still needs further verification.

Chemotherapy is widely utilized as a pivotal adjuvant therapy for cervical cancer treatment. Cisplatin exerts an antitumor effect via interaction with the purine bases in DNA to generate DNA lesions, and it is broadly used for the treatment of cervical, ovarian, and lung cancers, among others [41, 42]. However, drug resistance has been a major challenge, limiting the use and efficacy of cisplatin. Paclitaxel is another commonly used drug for cervical cancer treatment [43, 44], and it is usually combined with cisplatin, carboplatin, or bevacizumab. One clinical trial reported that docetaxel combined with carboplatin may provide a favorable control rate for stage IV cervical cancer [45]. We predicted that the drug sensitivity of cervical cancer to cisplatin, paclitaxel, and docetaxel would depend on C1QB expression, and we found that patients with a low expression of C1QB tended to have a high drug sensitivity. This result may be helpful in future formulations of chemotherapy regimens, although more experiments are still needed for validation.

In brief, we have demonstrated the overexpression of C1QB mRNA and protein in cervical cancer in both clinical specimens and multicenter samples. The clinical significance and underlying mechanisms of C1QB in cervical cancer were explored through multiple approaches. Nevertheless, our study also had some limitations. The main limitation was that we were unable to address the prognostic significance of C1QB due to the lack of patient follow-up information. More experiments are also needed to explore the molecular mechanisms of C1QB, both in vitro and in vivo.

5. Conclusions

Collectively, the data presented here provide the first verification of the overexpression of C1QB in cervical cancer at both the mRNA and protein levels in 1,341 samples. This trend toward upregulation may be related to the depth of invasion, lymphovascular invasion, and perineural invasion in cervical cancer. Our GSEA data also indicate that C1QB may participate in apoptosis and autophagy processes. In fact, C1QB may be an oncogene in cervical cancer, but this needs further study.

Data Availability

The datasets analysed during the current study are available in the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/), ArrayExpress (https://www.ebi.ac.uk/arrayexpress/), The Genotype-Tissue Expression (https://gtexportal.org/home/), The Cancer Genome Atlas (https://portal.gdc.cancer.gov/), and Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra). The accessions of datasets were listed in Table 4 of our article.

Ethical Approval

This study was approved by the Ethics Committee of Guangxi Medical University Cancer Hospital (No. 2022007). This study conformed to the standards set by the Declaration of Helsinki.

Informed consent was received from all participants.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors participated in the interpretation of the studies and review of the manuscript. Y.W.D. and G.C.: designed the study. S.T.L. and Z.Q.L.: performed the statistical analysis, prepared the figures and tables, and wrote the original draft. X.Q.Y., J.Y.L., and J.H.C.: collected the data. S.T.L., Y.Y.P., and X.Y.C.: conducted the experiments.

Acknowledgments

The authors would like to thank Guangxi Key Laboratory of Medical Pathology for technical supports. The authors disclosed receipt of the following financial support for the research, authorship, and publication of this article: this work was supported by: Guangxi Zhuang Autonomous Region Health Commission Self-financed Scientific Research Project [grant number: Z20210882, Z20200124, Z20180979, Z20190529], Guangxi Medical High-level Key Talents Training “139” Program [grant number: 2020], Guangxi Higher Education Undergraduate Teaching Reform Project [grant number: 2021JGA142], Guangxi Educational Science Planning Key Project [grant number: 2021B167], the Fund of Future Academic Star of Guangxi Medical University [grant number: WLXSZX21117], and Guangxi Medical University 2021 Undergraduate Innovation and Entrepreneurship Training Program [grant number: 202110598124].