Introduction {#S0001}
============

Cervical carcinoma is one of the most common genital-system cancers. Among various malignant tumors, cervical cancer has shown the largest increase in incidence among women, especially in China, and continues an increasing trend of mortality as one of the 10 most prevalent cancers worldwide.[@CIT0001] Although many therapeutic options for cervical cancer are available (surgical treatments, chemotherapy, and radiotherapy), selection of an optimal treatment strategy for patients at different stages is yet to be confirmed by further studies. For example, radical hysterectomy and pelvic external-beam radiation therapy play major roles in the treatment of stages IB1 and IIA1 of cervical squamous cell carcinoma (CSCC). Nonetheless, which of the two treatments is a better option is still unclear.

Recent evidence suggests that CXC chemokine receptor type 7 (CXCR7), a novel receptor of stromal cell-derived factor 1 (CXCL12), is strongly expressed in different tumor types and tumor-associated vasculature but is only weakly expressed in many normal tissues.[@CIT0002]--[@CIT0005] Furthermore, CXCR7 is involved in the initiation, progression, invasiveness, and distant metastasis of cancers of the skin, esophagus, lung, and colon. It performs crucial functions in the biological behavior of tumors. Nevertheless, studies on the functions and mechanisms of action of CXCR7 in CSCC are limited.

Some reports have shown that CXCR7, along with a signaling receptor of epidermal growth factor (EGFR), triggers intracellular signaling.[@CIT0006] EGFR and CXC chemokine receptor type 4 (CXCR4, a receptor of CXCL12) have been reported to be associated with distant metastases and disease-free survival (DFS) in cervical cancer. Accordingly, whether CXCR7 has a similar effect on CXCR4 should be determined. The present study was aimed at evaluating the protein expression of CXCR7 and of its putative coreceptor EGFR by immunohistochemical (IHC) analysis as well as at determining their relation with patient characteristics and survival in CSCC. Additionally, we evaluated the potential usefulness of CXCR7 as a prognostic factor and a screening index for treatment options at CSCC stages IB1 and IIA1. These data have not been reported previously.

Materials and Methods {#S0002}
=====================

Patients {#S0002-S2001}
--------

Between March 2005 and March 2011, data on 106 cases with normal cervical epithelia (NCE), 120 cases of cervical intraepithelial neoplasia (CIN), and 165 cases of CSCC were retrieved from the Pathology Department in The First Affiliated Hospital of Guangxi Medical University and were analyzed retrospectively. All the tissue samples were fixed in formalin and embedded in paraffin. The patients with CSCC ranged in age from 29 to 77 years, with a mean age of 48 years. Among the 165 cases of CSCC, according to the International Federation of Gynecology and Obstetrics 2009 system (FIGO) criteria, 67 cases were FIGO stage IB1--IIA1 (40.6%), 40 cases were FIGO stage IIB--III (24.2%), and 58 cases were FIGO stage IV (35.2%). No subject underwent chemotherapy, radiotherapy, or other adjuvant therapeutic modalities before the tissues were sampled during the surgical operation. After surgical operation, all CSCC patients with high risk factors received radiotherapy.

A total of 106 NCE samples were obtained from patients with uterine myoma who underwent routine hysterectomy. Among the 120 cases of CIN and CIN I, the patients were treated with the surface damage method, whereas CIN II and III patients were treated with focus excision. Of 67 patients at FIGO stages IB1 and IIA1, 33 underwent radical hysterectomy, pelvic lymph node dissection, and para-aortic lymph node sampling (RH group), whereas 34 patients received pelvic external-beam radiation therapy with brachytherapy (EBRT group, total point A dose: 80--85 Gy).

Collection of Clinical Data {#S0002-S2002}
---------------------------

All the data including age, histological grade, disease stage, resection margins, vasoinvasion status, parametrial infiltration depth, lymph node metastasis status, and tumor size were retrieved from a combination of clinical and histopathological data, outpatient-clinic medical records, and communication with the patients. According to multidetector-row computed tomography, round or irregular nodes of \>10 mm in the major axis in the abdomen and pelvis were clinically diagnosed as lymph node metastasis. The tumor depth was determined as the vertical distance from the epidermis (granular layer) to the deepest lower margin of the SCC. The parametrial infiltration depth was categorized as \<15 mm or ≥15 mm because some studies have shown that cervical cancer more than 15 mm deep is likely to metastasize.[@CIT0007]

IHC Staining {#S0002-S2003}
------------

Expression levels of CXCR7 and EGFR in the resected and biopsy specimens from 391 patients were analyzed by IHC staining. In this procedure, tissue blocks were sectioned (at 4 µm) and deparaffinized with xylene prior to rehydration in a gradient of ethanol solutions. Antigens were retrieved with ethylenediamine tetraacetic acid via heating of the tissue samples in a thermostat-controlled water bath at 95 °C for 15 min, followed by incubation with 3% H~2~O~2~ for 10 min to block endogenous peroxidase. After incubation with anti-CXCR7 (monoclonal mouse IgG1 clone No. 11G8, 1:100, R&D Systems, Minneapolis, MN, USA), anti-EGFR (monoclonal rabbit IgG1 clone No. 1902--1, 1:100, Epitomics, Burlingame, CA, USA), anti-phosphorylated-ERK (anti-p-ERK; monoclonal mouse IgG1 clone No. sc-136521, 1:100, Santa Cruz Biotechnology, Dallas, TX, USA), anti-ERK (monoclonal mouse IgG2a clone No. sc-514302, 1:100, Santa Cruz Biotechnology), anti-p-AKT (monoclonal mouse IgG2b clone No. sc-514032, 1:100, Santa Cruz Biotechnology), and anti-AKT (monoclonal rabbit IgG clone No. sc-8312, 1:100, Santa Cruz Biotechnology) antibodies at 4 °C overnight, every section was probed with an appropriate secondary antibody for 50 min at room temperature. The tissue sections were stained with 3,3-diaminobenzidine tetrahydrochloride for 5 min. A negative control involved PBS instead of the primary antibodies, and a positive control included breast cancer tissue in lieu of cervical tissue.

The sections were analyzed by two pathologists without knowledge of the patients' clinical data. Tumor cells with a brown cytoplasm or membrane were considered positively stained. The staining intensity was scored on a four-point scale: 0, none; 1, weak brown; 2, moderate brown; and 3, strong brown. The percentage of positive tumor cells was classified into four categories: 0, \<10%; 1, 10%--49%; 2, 50%--74%; and 3, \>75%. A staining index (SI), which was calculated by multiplying the staining intensity by the percentage of positive cells, was assigned to all cases. The SI of 0--3 was classified as negative staining, whereas an SI of more than 3 was classified as positive staining.

Survival Indices {#S0002-S2004}
----------------

DFS and overall survival (OS) were evaluated to determine the time point of disease recurrence and death from any cause, respectively. DFS was assessed from the date of treatment until the date of either regional recurrence or distant metastasis. The end point for OS was either the date of death or last follow-up examination. We obtained follow-up information in March 2018 through regular outpatient reviews and telephone interviews with the patients or their relatives.

Cell Culture and Transfection {#S0002-S2005}
-----------------------------

Two human CSCC cell lines (HeLa and SiHa) were purchased from the Shanghai Institute of Biochemistry and Cell Biology (Shanghai, China), and were grown in Dulbecco's modified Eagle's medium (DMEM) that was supplemented with 10% heat-inactivated fetal bovine serum (FBS) and 1% penicillin/streptomycin mixture (all from Gibco; Thermo Fisher Scientific, Inc., Waltham, MA, USA). All cells were kept at 37°C in a humidified incubator containing 5% CO~2~.

Small interfering (si)RNA targeting CXCR7 (si-CXCR7) and non-targeting siRNA (si-NC) were chemically synthesized by GenePharma Co., Ltd (Shanghai, China). For transfection, cells were seeded into 6-well plates, and were incubated at 37°C with 5% CO~2~ for about 12 h. Si-CXCR7 or si-NC was transfected into cells using Lipofectamine 2000™ reagent (Invitrogen; Thermo Fisher Scientific, Inc.) in accordance with the manufacturer's protocol.

Reverse-Transcription Quantitative PCR (RT-qPCR) {#S0002-S2006}
------------------------------------------------

The isolation of total RNA was performed from si-CXCR7 or si-NC-transfected HeLa and SiHa cells by means of a TRIzol reagent (Thermo Fisher Scientific, Waltham, MA). For the determination of CXCR7 mRNA expression, reverse transcription was conducted using a Prime Script RT reagent Kit (Takara, Tokyo, Japan). The synthesized cDNA was then subjected to amplification through a SYBR Premix Ex Taq (Takara). GAPDH acted as a endogenous control for normalizing the CXCR7 mRNA expression. Relative CXCR7 mRNA was analyzed using the 2^--ΔΔCq^ method.

Cell Counting Kit-8 Assay {#S0002-S2007}
-------------------------

Cellular proliferation was assessed through CCK-8 assay. In detail, si-CXCR7 or si-NC-transfected HeLa and SiHa cells were collected, resuspended and then inoculated into 96-well plates with a density of 2,000 cells/well. subsequent to different time incubation, each well was supplemented with 10 µL of the CCK-8 reagent (Dojindo Molecular Technologies, Inc., Kumamoto, Japan), after which cells were incubated at 37 °C in 5% CO2 incubator for additional 2 h. The absorbance of each well was read at 450 nm wavelength using a microplate reader (BioTek, Winooski, VT, USA).

Transwell Invasion Assay {#S0002-S2008}
------------------------

HeLa and SiHa cells transfected with si-CXCR7 or si-NC were harvested after 48 h culture, washed with DMFM and resuspended in serum-free DMEM. A total of 200 µL cell suspension containing 5×10^4^ cells was inoculated into the upper compartments of Transwell inserts that were precoated with Matrigel (BD Biosciences, Franklin Lakes, New Jersey, USA). The lower compartments were placed with 600 µL DMEM containing 10% FBS. After 24 h culture, the invasive cells were fixed with methanol and stained with 0.5% crystal violet, and photographed using an inverted light microscope (Olympus, Tokyo, Japan). Five visual fields of each insert insert were randomly selected and the average number of invasive cells was counted.

Statistical Analysis {#S0002-S2009}
--------------------

Continuous variables were evaluated by Student's *t* test, whereas enumeration data were analyzed using the *χ*^2^ test. One-way ANOVA followed by the Student--Newman--Keuls post hoc test was performed to analyze the differences among multiple groups ([Figure 1](#F0001){ref-type="fig"}). To investigate whether the expression of CXCR7 in CSCC is associated with EGFR expression, Spearman's rank correlation test was conducted. Univariate survival curves were constructed by the Kaplan--Meier method, and their differences were analyzed by the logrank test. Univariate and multivariate Cox regression models were employed to process the data by means of SPSS 18.0 software (Chicago, IL, USA). Statistical significance was set to P \< 0.05.Figure 1Immunohistochemical CXCR7 and EGFR staining of tissue samples of CSCC. (**A**) Representative tissue sections with positive CXCR7 staining. (**B**) Representative sections with negative CXCR7 staining. (**C**) Representative sections with positive EGFR staining. (**D**) Representative tissue sections with negative EGFR staining. (**E**) Representative sections with positive ERK staining. (**F**) Representative sections with negative ERK staining. (**G**) Representative sections with positive p-ERK staining. (**H**) Representative sections with negative p-ERK staining. (**I**) Representative sections with positive AKT staining. (**J**) Representative sections with negative AKT staining. (**K**) Representative sections with positive p-AKT staining. (**L**) Representative sections with negative p-AKT staining. (**M**) Statistical analysis of the IHC results on CXCR7 and EGFR. (**N**) The ratio of phosphorylated versus total protein for ERK and AKT. \*P \< 0.05 vs. NCE. Magnification: (**A**) ×400; (**B**) ×400; (**C)** ×400; (**D**) ×400; (**E**) ×400; (**F**) ×400; (**G**) ×400; (**H**) ×400; (**I**) ×400; (**J**) ×400; (**K**) ×400; (**L**) ×400.

Results {#S0003}
=======

Quantitation of CXCR7, EGFR, ERK, p-ERK, AKT, And p-AKT in NCE, CIN, and CSCC {#S0003-S2001}
-----------------------------------------------------------------------------

IHC staining of CXCR7 in CSCC tissue sections was detected in the cytoplasm ([Figure 1A](#F0001){ref-type="fig"}), and negative CXCR7 staining is shown in [Figure 1B](#F0001){ref-type="fig"}. The expression of CXCR7 was significantly higher in CSCC than in NCE and CIN (P \< 0.05; [Table 1](#T0001){ref-type="table"}). Staining for EGFR was the strongest in the membrane and cytoplasm in the CSCC group ([Figure 1C](#F0001){ref-type="fig"}), with [Figure 1D](#F0001){ref-type="fig"} illustrating negative EGFR staining. The positivity prevalence of EGFR in CSCC was 64.2% (106 of 165 cases). EGFR was expressed in 62 out of 120 CIN cases (51.7%). Of the 106 NCE tissue samples, only 39 (36.8%) yielded uniformly detectable EGFR immunostaining. The positivity prevalence rates of EGFR were significantly higher among the CSCC and CIN tissue samples than among NCE tissue samples (P \< 0.05). In addition, we determined the amounts of ERK, p-ERK, AKT, and p-AKT in the same groups of tissue samples. Similar to the expression status of EGFR, the amounts of ERK ([Figure 1E](#F0001){ref-type="fig"} and [F](#F0001){ref-type="fig"}; P \< 0.05), p-ERK ([Figure 1G](#F0001){ref-type="fig"} and [H](#F0001){ref-type="fig"}; P \< 0.05), AKT ([Figure 1I](#F0001){ref-type="fig"} and [J](#F0001){ref-type="fig"}; P \< 0.05), and p-AKT ([Figure 1K](#F0001){ref-type="fig"} and [L](#F0001){ref-type="fig"}; P \< 0.05) were all higher in CSCC and CIN tissue samples than in the NCE tissue samples. The statistical analysis of the IHC results on CXCR7 and EGFR is presented in [Figure 1M](#F0001){ref-type="fig"}. The ratio of phosphorylated to total protein for ERK and AKT was analyzed and is shown in [Figure 1N](#F0001){ref-type="fig"}.Table 1Correlation Between CXCR7 Expression and Clinical DataCharacteristicNCXCR7P Value−+Diagnosis0.016 NCE1066343 CIN1205565 CSCC1656996Mean age16548.147.50.612Histological grade0.457 G1311021 G2753243 G3592732Disease stage0.014 IB1--IIA1673730 IIB--III401228 IV582038Resection margins0.488 −1435885 +221111Vasoinvasion0.280 −1396178 +26818Parametrial infiltration depth0.751 \<15 mm984256 ≥15 mm672740Lymph node metastasis0.000 −71629 +94787Tumor size0.000 \<40 mm69627 ≥40 mm96789[^2]

Clinical Characteristics in Relation to the Expression of CXCR7 in the Tumors of Patients with CSCC {#S0003-S2002}
---------------------------------------------------------------------------------------------------

Clinical characteristics of the 165 patients with CSCC are summarized in [Table 1](#T0001){ref-type="table"}. The mean age at diagnosis was 48 years (range, 29--77 years). The expression of CXCR7 was significantly associated with the stage of CSCC, and strong staining for CXCR7 corresponded to more advanced stages (P \< 0.05; [Table 1](#T0001){ref-type="table"}). Moreover, CXCR7 status was associated with lymph node metastasis and tumor size ≥40 mm in cervical cancer (both P \< 0.05; [Table 1](#T0001){ref-type="table"}). CXCR7 expression significantly correlated with EGFR expression (r = 0.880, P \< 0.05). No other associations between CXCR7 and clinical characteristics were found.

DFS and OS in Relation to CXCR7 and EGFR Expression in CSCC {#S0003-S2003}
-----------------------------------------------------------

To determine whether the expression levels of CXCR7 and EGFR correlate with DFS or OS in CSCC, patients were classified into CXCR7-negative cases (n = 69), CXCR7-positive cases (n = 96), EGFR-negative cases (n = 59) without CXCR7 expression (n = 53), and EGFR-positive cases (n = 106) with CXCR7 expression (n = 96). Follow-up information revealed that the median survival time was 57 months (range, 5--156 months), and 52 patients (31.5%) died of CSCC. Kaplan--Meier analysis suggested that patients with the expression of CXCR7 and EGFR in the tumor had significantly shorter DFS and OS than those without this expression pattern (P \< 0.05; [Figure 2A](#F0002){ref-type="fig"}--[F](#F0002){ref-type="fig"}).Figure 2DFS among the patients with CSCC in the presence or absence of tumor CXCR7 expression (**A**), presence or absence of tumor EGFR expression (**C**), and presence or absence of tumor CXCR7 and EGFR expression (**E**). OS of patients with CSCC in the presence or absence of tumor CXCR7 expression (**B**), presence or absence of tumor EGFR expression (**D**), and presence or absence of tumor CXCR7 and EGFR expression (**F**).

Univariate Cox regression analysis of all the clinical characteristics indicated that CXCR7 (hazard ratio \[HR\] = 0.260, 95% confidence interval \[CI\] 0.117--0.376, P = 0.0002), EGFR (HR = 0.394, 95% CI 0.226--0.687, P = 0.001), tumor size (HR = 0.468, 95% CI 0.268--0.816, P = 0.007), and lymph node metastasis (HR = 0.484, 95% CI 0.277--0.845, P = 0.011) were independently associated with DFS ([Table 2](#T0002){ref-type="table"}). Univariate analysis revealed that the expression levels of CXCR7 and EGFR and tumor size were prognostic factors of OS (HR = 0.300 with P = 0.001; HR = 0.450 with P = 0.004; and HR = 0.533 with P = 0.024, respectively). Lymph node metastasis was also found to be a prognostic factor of OS (HR = 0.549; P = 0.032; [Table 2](#T0002){ref-type="table"}).Table 2Univariate and Multivariate Analyses of Clinical and IHC Factors Among the 165 Patients with CSCCFactorsStatusCutoffDFSOSHR95% CI*P* ValueHR95% CIP Value**Univariate analysis** CXCR7Negative vs positiveNegative0.2600.117--0.3760.0000.3000.281--0.4900.001 EGFRNegative vs positiveNegative0.3940.226--0.6870.0010.4500.260--0.7790.004 Tumor size\<40 mm vs ≥40 mm\<40 mm0.4680.268--0.8160.0070.5330.308--0.9220.024 Lymph node metastasisNegative vs positiveNegative0.4840.277--0.8450.0110.5490.317--0.9510.032**Multivariate analysis** CXCR7Negative vs positiveNegative0.3510.200--0.6140.0000.4050.233--0.7020.001[^3][^4]

Multivariate Cox regression analysis revealed that only CXCR7 (HR = 0.351, 95% CI 0.200--0.614, P = 0.0007) was an independent predictor of DFS after adjustment for EGFR, tumor size, and lymph node metastasis status. In multivariate Cox proportional hazards analysis, CXCR7 was a significant prognostic factor of OS in CSCC (HR = 0.405, 95% CI 0.233--0.702, P = 0.001; [Table 2](#T0002){ref-type="table"}).

Survival Analysis for the Two Therapeutic Options (RH and EBRT) at Different CXCR7 Expression Levels {#S0003-S2004}
----------------------------------------------------------------------------------------------------

To test whether CXCR7 expression is useful as a screening index for treatment options at stages IB1 and IIA1 of CSCC, groups were compiled of CXCR7-negative cases (n = 37), including the RH group (n = 16) and EBRT group (n = 21), and CXCR7-positive cases (n = 30), including the RH group (n = 17) and EBRT group (n = 13). The median follow-up duration among patients with stage IB1 or IIA1 CSCC was 45 months (95% CI 41.4--55.3). The median follow-up duration was 43 and 48.5 months, for the RH group (95% CI 34.5--55.5) and EBRT group (95% CI 42.1--61.2), respectively. Furthermore, 19 patients (28.4%) died of stage IB1 or IIA1 CSCC. Kaplan--Meier analysis showed that cases with CXCR7 expression manifested no significant difference in OS and DFS between groups RH and EBRT (P = 0.141). Kaplan--Meier analysis of the cases of negative CXCR7 expression yielded similar results (P = 0.299).

Survival Outcomes {#S0003-S2005}
-----------------

The details of failure patterns of the two therapeutic options at different expression levels of CXCR7 are listed in [Table 3](#T0003){ref-type="table"}. A total of 33 patients experienced locoregional failure, and the most common failure site was the vagina. Approximately 36 cases experienced distant metastasis, and the most common failure site was the lungs. Among the 37 CXCR7-negative cases, no significant difference was found in the failure pattern between the RH group and EBRT group (recurrence, P = 0.104; locoregional failure, P = 0.508; and distant metastasis, P = 0.743). Among the 30 CXCR7-positive cases, the prevalence rates of recurrence and distant metastasis were significantly lower in the EBRT group than in the RH group (P = 0.007 and P = 0.025, respectively; [Table 3](#T0003){ref-type="table"}). By contrast, no significant difference between the RH group and EBRT group was found in terms of locoregional failure (P = 0.484).Table 3Failure Patterns of the Two Treatment Options at Different Expression Levels of CXCR7Failure PatternNCXCR7 (−)P ValueNCXCR7 (+)P ValueRHEBRTRHEBRTRecurrence0.1040.007 −171071028 +2061420155Locoregional failure0.5080.484 −2010101495 +176111688Distant metastasis0.7430.025 −187111349 +1991017134[^5]

Roles of CXCR7 Knockdown in CSCC {#S0003-S2006}
--------------------------------

To test the functions of CXCR7 in CSCC progression, we silenced its expression in HeLa and SiHa cells via transfecting with si-CXCR7. The HeLa and SiHa cells that were treated with si-NC served as the control. RT-qPCR analysis confirmed the successful knockdown of CXCR7 in HeLa and SiHa cells following si-CXCR7 transfection ([Figure 3A](#F0003){ref-type="fig"}, P \< 0.05). CCK-8 and Transwell invasion assays were respectively applied to examine the impacts of CXCR7 knockdown on HeLa and SiHa cells. The proliferative ([Figure 3B](#F0003){ref-type="fig"}, P \< 0.05) and invasive ([Figure 3C](#F0003){ref-type="fig"}, P \< 0.05) capacities of HeLa and SiHa cells were turned out to be obviously hindered after CXCR7 silencing. These results collectively demonstrated that CXCR7 plays a cancer-promoting part in CSCC progression.Figure 3The impacts of CXCR7 knockodwn on HeLa and SiHa cells. (**A**) RT-qPCR analysis was utilized to detect CXCR7 mRNA expression in HeLa and SiHa cells after si-CXCR7 or si-NC transfection. (**B**) CCK-8 assay was performed to measure the proliferation of CXCR7 deficient-HeLa and SiHa cells. (**C**) The invasive ability of si-CXCR7 or si-NC-transfected HeLa and SiHa cells was examined via Transwell invasion assay. \*P \< 0.05 vs. si-NC.

Discussion {#S0004}
==========

Aside from CXCR4, CXCR7 was recently identified as a novel G protein--coupled seven-span transmembrane receptor of CXCL12, which is strongly expressed in many cancers, such as esophageal squamous cell carcinoma, bronchogenic carcinoma, pancreatic cancer, and bladder cancer.[@CIT0008]--[@CIT0010] In the present study, IHC staining proved that CXCR7 is frequently expressed in the cytoplasm or membrane of cervical epithelial cells, suggesting that the source of CXCR7 in the tumor is tumor cells and that CXCR7 may play a part in tumorigenesis.

Furthermore, this study addressed CXCR7 expression in NCE, CIN, and CSCC and revealed that the highest expression of CXCR7 among the three cervical epithelial tissue types is present in CSCC lesions. Similar results have been reported on bladder cancer.[@CIT0010] Besides, our study shows that CXCR7 manifests stage-dependent expression in CSCC. High CXCR7 expression indicated more lymph node metastases and greater tumor size in patients with CSCC. These results suggest that CXCR7 expression may be conducive to metastasis and aggressive tumor growth, for the following reasons. First, rapid growth of tumor tissue results in a microenvironment featuring local hypoxia. Hypoxic preconditions induce the expression of *CXCR7* mRNA in squamous epithelial cervical carcinoma cells, which resemble bone marrow--derived mesenchymal stem cells. Moreover, the high expression of CXCR7 in malignant tumors may be due to hypoxia; this finding has been confirmed by several studies.[@CIT0003] Second, CXCR7 might participate in the formation of (and angiogenesis in) CSCC because CXCR7 is strongly expressed in tumor-associated blood vessels in breast and prostate cancer tissues.[@CIT0011],[@CIT0012] Other researchers have reported that selective CXCR7 upregulation in (or application of CXCR7 to) hypoxic lungs has a significant effect on pulmonary vascular-cell proliferation and vascular remodeling. Third, CXCR7 can promote the migration and adhesion of squamous epithelial cervical carcinoma cells, and a decrease in CXCR7 expression reduces the number of tumor cells.[@CIT0013]

Our data also mean that CSCC tumors overexpress EGFR, which was more strongly expressed in cervical cancer tissues than in the other two types of tissue. Of note, CXCR7 was found to be statistically significantly coexpressed with EGFR in the tumors; both were associated with DFS and OS in CSCC. The high expression of CXCR7 was related to shorter survival, and survival was poor among the patients with CSCC and positive EGFR staining in the tumor, indicating that CXCR7 may form heterodimers with EGFR for optimal intracellular signaling.[@CIT0014],[@CIT0015] Although our present study supports this hypothesis, functional studies on cervical cancer cell lines should be conducted to confirm this conclusion. Accordingly, a series of functional experiments were performed to determine the functions of CXCR7 on CSCC progression. The results indicated that interference of CXCR7 expression resulted in an obvious decrease of CSCC cell proliferation and invasion. Hence, CXCR7 worked as an oncogene in the malignancy of CSCC, and might be a potential target for treating patients with this disease.

Additionally, we noted that CXCR7, EGFR, tumor size, and lymph node metastasis were prognostic risk factors of DFS and OS in our univariate analysis. Only CXCR7 was an independent prognostic factor in our multivariate analysis, which suggested that tumor CXCR7 expression may represent a subset of "high-risk" patients or could serve as a prognostic biomarker of CSCC. On the other hand, our results are inconsistent with some studies.[@CIT0008],[@CIT0016] These discrepancies might be due to different histological sources or different patient enrollment conditions leading to variation in prognoses. CSCC is a squamous epithelial malignant tumor, whereas hepatic cellular cancer (in the other studies) is an epithelioglandular tumor. Therefore, the secretion of CXCR7 may be affected by different histological sources, thereby leading to the discrepancy in results between this and other studies. In addition, factors such as age, medical history, and clinical stage possibly affect the DFS and OS of patients with malignant tumors.

A number of treatment options are available for the FIGO stages IB1 and IIA1 of CSCC. Nonetheless, selecting a feasible therapeutic option is an issue that clinicians must address. Our study shows for the first time that external-beam radiation therapy with brachytherapy is a feasible therapeutic regimen for the FIGO stages IB1 and IIA1 of CSCC with high expression of CXCR7 in the tumor.

This study revealed no survival advantage between groups RH and EBRT during either the presence or absence of CXCR7 expression. The 37 patients with negative CXCR7 staining exhibited no significant difference between the RH group and EBRT group in terms of failure patterns. Among the 30 patients with CXCR7 expression, the failure patterns of RH and EBRT for CSCC stages IB1 and IIA1 were compared, and the data showed that pelvic external-beam radiation therapy with brachytherapy significantly reduced the probability of recurrence and distant metastasis. Therefore, to clinicians, when high CXCR7 expression is present, pelvic external-beam radiation therapy with brachytherapy may be recommended to minimize recurrence and promote local control. Nevertheless, clinical trials are needed to analyze the postoperative adverse reactions and side effects of pelvic external-beam radiation therapy with brachytherapy; the resultant data should make the screening index based on CXCR7 highly clinically relevant and convincing.

In conclusion, the findings of this retrospective study suggest that CXCR7 may cause aggressive tumor growth and may be a valuable prognostic biomarker in CSCC. CXCR7 and EGFR turned out to be related to cervical carcinogenesis and associated with CSCC progression, which may involve the activation of the ERK or AKT pathway. These observations could offer new molecular targets for treating CSCC. Nevertheless, further studies are necessary to identify the exact mechanisms of CXCR7 and EGFR actions on CSCC progression. Moreover, CXCR7 may serve as a screening index for treatment options (RH and EBRT) among patients with stage IB1 or IIA1 CSCC. Randomized controlled trials are necessary to find the optimal therapeutic regimen in the case of CXCR7 overexpression owing to the limitations of the retrospective design, our small sample size, and a lack of patients receiving pelvic external-beam radiation therapy and brachytherapy with or without concurrent chemotherapy.

Abbreviations {#S0005}
=============

CI, confidence interval; CIN, cervical intraepithelial neoplasia; CSCC, cervical squamous cell carcinoma; DFS, disease-free survival; EBRT group, pelvic external-beam radiation therapy with brachytherapy; FIGO, International Federation of Gynecology and Obstetrics 2009 system; HR, hazard ratio; IHC, immunohistochemical; NCE, normal cervical epithelia; OS, overall survival; p-AKT, phospho-AKT; RH group, radical hysterectomy, pelvic lymph node dissection, and para-aortic lymph node sampling; SI, staining index.

Consent for Publication {#S0006}
=======================

All the patients approved the publication of the images associated with this study.

Ethics Approval and Consent to Participate {#S0007}
==========================================

The study protocol was approved by the Ethics Committee of The First Affiliated Hospital of Guangxi Medical University (Guangxi, China). Every patient provided verbal informed consent because of the absence of an existing standard for written consent in our department.

Data Sharing Statement {#S0008}
======================

The datasets used and/or analyzed during the present study are available from the corresponding authors on reasonable request.

Author Contributions {#S0009}
====================

All authors contributed to data analysis, drafting and revising the article, gave final approval of the version to be published, and agree to be accountable for all aspects of the work.

Disclosure {#S0010}
==========

The authors report no conflicts of interest in this work.

[^1]: These authors contributed equally to this work

[^2]: **Abbreviations:** CIN, cervical intraepithelial neoplasia; CSCC, cervical squamous cell carcinoma; CXCR, chemokine receptor; EGFR, epidermal growth factor receptor; NCE, normal cervical epithelia.

[^3]: **Note:** Method: Cox regression model.

[^4]: **Abbreviations:** CI, confidence interval; CXCR, CXC chemokine receptor; DFS, disease-free survival; EGFR, epidermal growth factor receptor; HR, hazard ratio; OS, overall survival.

[^5]: **Abbreviations:** CXCR, chemokine receptor; EBRT, pelvic external-beam radiation therapy with brachytherapy; RH, radical hysterectomy, pelvic lymph node dissection, and para-aortic lymph node sampling.
