This article is mentioned in:

Introduction

Cervical cancer is the fourth leading cause of cancer-associated mortality worldwide, puts an increasing burden on low- and high-income countries alike, and was associated with >500,000 new cases and 300,000 deaths worldwide in 2018. Organized screening and human papillomavirus (HPV) vaccination can effectively reduce the incidence and mortality of cervical cancer (1). However, the mortality rate shows an increasing trend in underdeveloped and developing countries (2). Cervical cancer is commonly caused by high-risk human papillomavirus (HR-HPV) infection, and HPV DNA is found in >90% of malignant cervical lesions. HPV16 and 18, as HR-HPV, account for ~70% of cases, and HPV16 is the most common in female patients with cervical cancer (3). Therefore, it is essential to study the mechanism of HPV16 infection-induced cervical cancer.

Squamous cell carcinoma is the most common type of cervical cancer, accounting for ~75% of cases, and involves a series of precancerous lesion stages, namely, low-grade squamous intraepithelial lesion (LSIL) and high-risk squamous intraepithelial lesion (HSIL) (4). Persistent infection with high-risk HPV can lead to malignant transformation of the normal cervix, although this can take several years to decades (5–7). HPV infection can lead to a tissue inflammatory response, activating the Janus kinase (JAK)/STAT signaling pathway and promoting HR-HPV replication (8). Previous studies have mainly focused on the association between HPV16 and STAT3, whereas few studies have considered STAT1 (9–12).

JAK/STAT signaling is involved in numerous physiological cell processes, such as proliferation, differentiation, apoptosis and immune system regulation; however, aberrant STAT regulation may lead to pathological events, including malignant cell transformation and metastasis (13). STAT1 was the first STAT protein identified in the interferon (IFN) signal transduction pathways, and is a 91-kDa cytosolic protein that adopts a structure with six highly conserved domains (14). It serves important roles in cytokine-induced signaling pathways, and can act as an antiviral and antibacterial mediator, growth inhibitor and inducer of apoptosis (13). STAT1 may have different roles in the progression of cancer. For example, various studies have indicated that STAT1, as a tumor suppressor, restrains the proliferation of cancer cells, including lung cancer cells, esophageal squamous cell carcinoma cells and colorectal cancer cells (15–17). In addition, STAT1 knockout mice are more susceptible to experimentally induced tumors, including the development of mammary adenocarcinoma (18). By contrast, aberrant STAT1 expression has been identified in human cancer, including pleural mesothelioma, renal cell carcinoma and breast cancer (19–21). Patients with high levels of STAT1 and/or pSTAT1 in cancer tissues have poor clinical outcomes compared with patients with low levels (22,23). However, to the best of our knowledge, the effect of STAT1 on the progression of cervical cancer remains unknown.

The present study used human cervical tissues to investigate the associations among STAT1 expression, human cervical lesion progression and HPV16 viral load. The present study aimed to determine the association between STAT1 expression and the grade of cervical disease in HPV-infected patients.

Materials and methods

Study population

The present study included 141 cervical tissue specimens from patients treated at Shengjing Hospital, China Medical University (Shenyang, China) between January 2012 and December 2017. The specimens were divided into four groups: Normal, LSIL, HSIL and cancer. In the normal, LSIL and HSIL groups, cold knife conization was used to obtain the samples, while in the cancer group, samples were obtained during hysterectomy. The age of patients ranged between 20 and 70 years (median age, 43 years). Each HPV16-infected cervical tissue specimen underwent HPV flow-through hybridization detection and paraffin histopathological evaluation. The histopathological diagnosis was performed by two pathologists, the percentage of abnormal cells was recorded and the formalin-fixed, paraffin-embedded (FFPE) cervical specimens were graded histologically. All samples were collected with written informed consent from patients, and all related procedures were approved by the Ethics Committee of Shengjing Hospital, China Medical University (approval no. 2019PS476K).

HPV genotyping

An GenoArray test kit (Hybribio Ltd.) was used to extract genomic DNA and perform HPV genotyping according to the manufacturer's protocols. The quality of the extracted DNA was ascertained by an internal control (β-globin primers), which was included in PCR mix. The amplification protocol was as follows: 9 min of denaturation at 95°C and 40 cycles of 20 sec of denaturation at 95°C, 30 sec of annealing at 55°C, and 30 sec of elongation at 72°C, followed by a final extension for 5 min at 72°C. This method was based on a flow-through hybridization principle and is able to classify 21 HPV genotypes, including 13 high-risk subtypes (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59 and 68), 5 low-risk subtypes (6, 11, 42, 43 and 44) and 3 subtypes that are common in China (53, 66 and CP8304) (24).

Immunohistochemical staining

STAT1 expression was examined using FFPE specimens. The tissue slices (3-µm thick) were deparaffinized in xylene and hydrated in an alcohol gradient series (100% for 10 min, 95% for 5 min, 85% for 5 min and 75% for 5 min). Specimens were boiled in citrate buffer (10 mM; pH 6.0) at 120°C for 5 min. Subsequently, 3% hydrogen peroxide was added dropwise and sections were incubated for 20 min at room temperature. After blocking in 5% BSA for 30 min at room temperature, STAT1 primary antibody (cat. no. AHO0832; dilution, 1:200; Thermo Fisher Scientific, Inc.) was used, and PBS was used as a negative control in place of the primary antibody. Subsequently, the sections were incubated in a refrigerator at 4°C overnight. The sections were washed with PBS and were incubated with horseradish peroxidase secondary antibody (cat. no. PV-6000; ZSGB-BIO, Inc.) at room temperature for 30 min. Subsequently, 3,3′diaminobenzidine was added for 30 sec for staining, and hematoxylin was used to counterstain the nuclei for 20 sec at room temperature. Finally, neutral gum mounting medium was used to seal the slides and images were captured with a DP70 digital camera-equipped light microscope (Olympus Corporation). The score for each tissue image was evaluated by two independent pathologists and the average score was used as the final score. A total of 10 high-power fields were randomly observed for each specimen, and the percentage of positively stained cells in each high-power field was calculated and scored. The interpretation criteria were the product of the staining intensity and the staining area score. The staining intensity was scored on a scale of 1–4: 1, no staining; 2, weak staining; 3, medium staining; 4, strong staining. The percentage of stained cells was scored on a scale of 1–4: 1, <5%; 2, 5–25%; 3, 25–50%; 4, ≥50%). The final product score was determined by multiplying these two values to obtain a score between 1 and 16: <6, low expression (−); ≥6, high expression (+).

Reverse transcription-quantitative (RT-q) PCR

Frozen tissue samples from 20 patients per group were selected for RT-qPCR. Total RNA was extracted from human tissues using the GeneJET RNA Purification kit (cat. no. K0731; Thermo Fisher Scientific, Inc.) according to the manufacturer's protocol. The cDNA was synthesized using a Prime Script RT-PCR kit (Takara Bio, Inc.). Samples were incubated with reverse transcriptase at 37°C for 15 min and 85°C for 5 sec, and the reverse program stopped at 4°C. qPCR assays were performed using a SYBR Premix ExTaq kit (Takara Bio, Inc.) and a 7,500-sequence detector (Applied Biosystems, Inc.). Amplification conditions were as follows: Initial step at 95°C for 30 sec, followed by 40 cycles of 95°C for 5 sec and 60°C for 34 sec. The sequences of the primers were as follows: STAT1, forward 5′-CTTACCCAGAATGCCCTGAT-3′, reverse 5′-CGAACTTGCTGCAGACTCTC-3′; and GAPDH, forward 5′-ACTGCCAACGTGTCAGTGGT-3′ and reverse 5′-GTGTCGCTGTTGAAGTCAGA-3′. GAPDH was used as the normalizing gene for determination of ΔCT values. Fold changes in gene expression were compared using 2−ΔΔCT values (25).

DNA isolation and HPV16 viral load quantification

DNA isolations in the FFPE samples were performed using a Qiagen FFPE kit (56404; Qiagen GmbH) as previously described (26). The present study evaluated the expression levels of HPV16 L1 to represent HPV16. Using GAPDH as an internal reference, fluorescent qPCR (SYBR-Green PCR kit; Qiagen GmbH) was used to quantitate HPV16 DNA content in tissues. The amplification assay was performed on an I-cycler (Bio-Rad Laboratories, Inc.) with the following conditions: Pre-denaturation at 95°C for 2 min, followed by 45 amplification cycles of 95°C for 5 sec, 60°C for 1 min, and with a final extension step at 72°C for 10 sec. The nucleotide sequences of the HPV16 L1 forward and reverse primer set used for amplification were 5′-ACTATTTTGGAGGACTGG-3′ and 5′-TCTTTAGGTGCTGGAGGT-3′, respectively. The relative amount of DNA was calculated using the 2−ΔΔCT method.

Bioinformatics analysis

Gene Expression Profiling Interactive Analysis (GEPIA; http://gepia.cancer-pku.cn/) is a visualization website based on The Cancer Genome Atlas data (27). GEPIA was used to verify the differential expression of STAT1 in cervical squamous cell carcinoma (CESC) and normal tissues. Expression data were first transformed using the formula log2 (TPM + 1) for the differential analysis, and log2FC was defined as the difference between the median of the tumor and the median of the normal tissue. An online KM plotter database (http://kmplot.com/analysis/, 2020.05.25) was used to analyze the overall survival rate of patients with different expression levels of STAT1 based on the clinical data of 304 patients with CESC (28). STAT1 was entered in the database to acquire CESC KM survival plots. The patients were divided into two cohorts according to the median expression of the gene. A logarithmic rank test was used for survival analysis. Hypothesis testing was performed, and the hazard ratio (HR) and 95% CI were calculated.

Statistical analysis

Pearson's χ2 test was used to compare the characteristics of patients according to STAT1 expression and HPV viral load. All analyses were performed using the SPSS statistical package (version 25.0; IBM Corp.). Data were presented as the means ± standard deviations. One-way ANOVA followed by Dunnett's post hoc test was used to determine statistically significant differences. To analyze the relationship between STAT1 protein expression and HPV viral load, Spearman's correlation analysis was performed for protein expression and HPV16 DNA load in the same samples using GraphPad Prism 6.0 software (GraphPad Software, Inc.). Receiver operating characteristic (ROC) curves and area under the curve (AUC) values were also calculated for STAT1 expression. All tests were two-tailed, and P<0.05 was considered to indicate a statistically significant difference.

Results

Association between STAT1 expression and the extent of cervical lesions

To investigate the expression levels of STAT1, immunohistochemical staining of 141 cervical tissue samples was performed. The results revealed that STAT1 expression gradually increased with the progression of cervical lesions (Fig. 1A and B). STAT1 protein expression was observed in 10.71 (3/28), 41.18 (14/34), 53.06 (26/49) and 90.00% (27/30) of normal cervical squamous epithelium, LSIL, HSIL and cervical squamous cell carcinoma samples, respectively (Table I). As shown in Fig. 1B, compared with normal tissues, STAT1 expression was significantly increased in LSIL (P=0.007; χ2=7.160), HSIL (P<0.001; χ2=13.610) and cancer samples (P<0.001; χ2=36.460). The proportion of samples with positive STAT1 expression was higher in HSIL than in LSIL, although the difference was not statistically significant (P=0.287; χ2=1.136). Furthermore, the proportion of samples with positive STAT1 expression was higher in cervical cancer than in HSIL (P<0.001; χ2=11.499). Frozen tissue samples from 20 patients per group were selected for RT-qPCR to detect STAT1 mRNA expression. Compared with that in cancer samples, STAT1 mRNA expression was markedly decreased in normal, LSIL and HSIL samples (Fig. 1C). These findings demonstrated that STAT1 expression gradually increased with increasing severity of lesions in cervical lesion tissues.

Figure 1. STAT1 expression in cervical tissues. (A) Immunohistochemical staining of STAT1 protein expression in normal cervical squamous epithelium, LSIL, HSIL and cervical squamous cell carcinoma. Magnification, ×200. (B) Box plot of the immunohistochemical staining scores of STAT1 in normal tissues, LSIL, HSIL and cancer. Statistical analysis was performed using the χ2 test. (C) Total RNA was isolated to analyze STAT1 mRNA expression (20 cases per group). Data are presented as the mean ± SD. ***P≤0.001. HSIL, high-grade squamous intraepithelial lesion; LSIL, low-grade squamous intraepithelial lesion.

Table I. Analysis of STAT1 expression and clinical factors of patients with cervical lesions.

Table I.

Analysis of STAT1 expression and clinical factors of patients with cervical lesions.

		STAT1
Characteristics	Total, n	High expression, n	Low expression, n	P-value
Age, years				0.147
≤43	63	27	36
>43	78	43	35
Lesion grade				<0.001
Normal	28	3	25
LSIL	34	14	20
HSIL	49	26	23
Cancer	30	27	3

[i] HSIL, high-grade squamous intraepithelial lesion; LSIL, low-grade squamous intraepithelial lesion.

Associations among HPV16 DNA load, extent of cervical lesions and STAT1 expression

To determine the distribution of HPV16 in the study samples, HPV16 L1 was detected. All normal group tissues were confirmed to be HPV16-negative. PCR detection was performed for 34 cases in the LSIL group, 49 in the HSIL group and 30 in the cancer group. The viral loads gradually increased from LSIL to HSIL to cervical squamous cell carcinoma, with mean values of 26.98±12.34, 51.30±21.377 and 169±79.47, respectively. These data demonstrated that, as the disease progressed, HPV16 viral load gradually increased (Fig. 2A). Spearman's correlation analysis was used to analyze the correlation between HPV16 viral load and STAT1 expression. The results demonstrated that HPV16 viral load was positively correlated with STAT1 expression (Fig. 2B).

Figure 2. HPV16 viral load in cervical lesions. (A) HPV16 L1 reverse transcription-quantitative PCR using GAPDH as an internal control. HPV16 viral load was assessed using the 2−ΔΔCT method. One-way ANOVA followed by Dunnett's post hoc test was used for statistical analysis. (B) Correlation between HPV16 DNA load and STAT1 expression in different cervical lesions. Spearman's correlation analysis was performed. **P<0.01. HPV, human papillomavirus; HSIL, high-grade squamous intraepithelial lesion; LSIL, low-grade squamous intraepithelial lesion.

Associations among STAT1 expression, HPV16 viral load and the severity of cervical lesions in female patients

To elucidate the associations among STAT1 protein expression, HPV16 viral load (individually or in combination) and the severity of the cervical lesions in patients, ROC analysis was performed. The AUC values for distinguishing LSIL and HSIL samples from cervical cancer samples were 0.824 (95% CI, 0.719-0.929) for STAT1 protein expression and 0.736 (95% CI, 0.598-0.873) for HPV16 viral load (Table II). The AUC values for distinguishing LSIL samples from HSIL and cervical cancer samples were 0.711 (95% CI, 0.592-0.830) for STAT1 protein expression and 0.698 (95% CI, 0.577-0.819) for HPV16 viral load. The AUC values were improved by using a combination of the STAT1 protein expression and HPV16 viral load, with an effective distinction between LSIL and HSIL samples and cervical cancer samples (sensitivity, 76.0%; specificity, 85.2%) obtained using the mean change frequencies of STAT1 protein expression and HPV16 viral load (AUC, 0.849; 95% CI, 0.749-0.949) when a 29% change frequency was used as the critical threshold (Fig. 3A; Table II). LSIL samples were distinguishable from cervical cancer and HSIL samples with a sensitivity of 56.6% and a specificity of 84.6% using the mean change frequencies of STAT1 protein expression and HPV16 viral load (AUC, 0.739; 95% CI, 0.627-0.852) when a 69% change frequency was used as the critical threshold (Fig. 3B; Table II). Concordant trends were observed when distinguishing HSIL samples from cervical cancer samples using a 32% change frequency as the critical threshold, which yielded a sensitivity of 88.0% and a specificity of 67.9% (AUC, 0.833; 95% CI, 0.718-0.947; Fig. 3C; Table II). The results revealed that STAT1 protein expression and HPV16 viral load, individually or in combination, could be used to differentiate between samples from patients with LSIL/HSIL and samples from patients with cervical cancer.

Figure 3. Receiver operating characteristic curves for the distinction between tissue samples using STAT1 expression and HPV16 viral load individually or in combination. (A) LSIL and HSIL vs. cervical cancer. (B) LSIL vs. HSIL and cervical cancer. (C) HSIL vs. cervical cancer. HPV, human papillomavirus; HSIL, high-grade squamous intraepithelial lesion; LSIL, low-grade squamous intraepithelial lesion.

Table II. AUC values obtained from receiver operating characteristic curve analysis of individual or combined STAT1 expression and HPV16 viral load.

Table II.

AUC values obtained from receiver operating characteristic curve analysis of individual or combined STAT1 expression and HPV16 viral load.

Group	AUC (95% CI)	Sensitivity, %	Specificity, %	Threshold, %
LSIL and HSIL vs. cervical cancer
STAT1	0.824 (0.719-0.929)	72.0	85.2	37
HPV16	0.736 (0.598-0.873)	48.0	94.4	33
STAT1 and HPV16	0.849 (0.749-0.949)	76.0	85.2	29
LSIL vs. HSIL and cervical cancer
STAT1	0.711 (0.592-0.830)	43.4	88.5	76
HPV16	0.698 (0.577-0.819)	86.8	50.0	60
STAT1 and HPV16	0.739 (0.627-0.852)	56.6	84.6	69
HSIL vs. cervical cancer
STAT1	0.806 (0.680-0.931)	72.0	82.1	50
HPV16	0.713 (0.563-0.863)	48.0	96.4	51
STAT1 and HPV16	0.833 (0.718-0.947)	88.0	67.9	32

[i] AUC, area under the curve; HPV, human papillomavirus; HSIL, high-grade squamous intraepithelial lesion; LSIL, low-grade squamous intraepithelial lesion.

Bioinformatics analysis

The present study used the GEPIA database to analyze STAT1 expression in CESC tissues, and non-cancerous tissues, and identified that STAT1 expression was significantly higher in CESC tissues than in non-cancerous tissues (Fig. 4A). To clarify the effect of gene expression on survival time, Kaplan Meier-plotter, an online survival analysis tool, was used to analyze overall patient survival rates based on STAT1 expression using the clinical data of 304 patients with cervical cancer, of which there were 204 in the high expression group and 100 in the low expression group (HR, 0.54; P=0.015; Fig. 4B). The results revealed that patients with cervical cancer with high STAT1 expression had longer survival.

Figure 4. Bioinformatics analysis of STAT1 expression in cervical tissues. (A) Box plot of STAT1 expression in CESC tissues and non-cancerous tissues. (B) Survival curves of patients with cervical cancer grouped by STAT1 expression. *P<0.05. CESC, cervical squamous cell carcinoma; HR, hazard ratio.

Discussion

The present study revealed that increased STAT1 expression was associated with the progression of cervical lesions. STAT1 expression in cervical epithelial lesions was significantly higher than in normal cervical epithelium samples. In a study on breast cancer by Hix et al (29), STAT1 was also highly expressed in mouse breast cancer cells, enhancing the growth and invasiveness of cancer cells, while the growth and invasiveness were decreased in STAT1 knockout cancer cells. Therefore, high STAT1 expression could promote the initiation and progression of tumors. However, STAT1 is also widely recognized as a tumor suppressor, and some studies have found that its expression is increased in early cervical intraepithelial neoplasia (CIN) I/II cervical lesions, decreased in CIN III/cervical carcinoma in situ and significantly upregulated in invasive cervical cancer. These findings suggest that STAT1 may have dual effects on cervical tumorigenesis. In early stages, STAT1 has a protective effect, while it is carcinogenic in aggressive tumors (30,31). STAT1 activation and enhanced expression in tumor cells are induced by different stimuli, such as IFN-α, IFN-β, IFN-γ, the action of oncogenes or loss of tumor suppressors, interaction with the tumor stroma, or genotoxic agents (32). Phosphorylated STAT1 is the activated form of STAT1. The present study detected the total STAT1 protein levels because both phosphorylated and unphosphorylated proteins could affect the outcome of cervical cancer, as described previously for breast cancer (33) and ovarian cancer (34). Tumors with a high STAT1 levels express IFN-related DNA damage resistance signature, which is linked to therapy resistance. Translation, autophagy and apoptosis in tumor cells are also altered by STAT1 and may contribute to therapy resistance. Furthermore, STAT1 enhances the invasion and metastasis of tumor cells (32).

HPV is a circular double-stranded DNA virus and HR-HPV is closely associated with the occurrence of cervical cancer (35). As a product of the number of infected cells and the number of viruses per infected cell, HPV viral load is influenced by two main factors: The extent of an HPV infection on the cervical surface and the level of viral production in the infection area. Changes in viral load and the association of these changes with disease risk may indicate that there is a complex interaction between HPV and the human host, and thus such changes could potentially serve as an additional predictive marker for the outcome of infection (36). Among the HPV subtypes, the number of patients with cervical cancer caused by HPV16 is the greatest (37). In 1995, the World Health Organization defined HPV16 as a viral tumor promoter. Due to the carcinogenicity of HPV, detection technology for HPV has gradually improved, and accumulated evidence indicates that the sensitivity of HPV DNA detection (97.6%) is higher than that of cytology. Therefore, in a number of countries, HPV DNA testing has been recommended as a primary screening method for cervical cancer (38). Furthermore, HPV16 DNA viral load in paraffin-embedded human cervical tissue samples was determined. Through qPCR, the present study revealed that viral load increased with the progression of cervical lesions. This is consistent with the findings of Camus et al (39) regarding cervical lesions. The development of cervical cancer is usually associated with an infection with HPV, particularly with a high-risk genotype of HPV (1). HPV16 is one of the most common genotypes, found in >35% of patients with cervical cancer (40,41). Changes in viral load and the associations of these changes with disease risk may imply a complex interaction between HPV and the human host, and potentially serve as an additional predictive marker for the outcomes of infection. A previous study suggested that HPV16 is not only particularly common, but also more oncogenic than other HPV types (42).

It has been reported that STAT1 expression is reduced in HFK cells transfected with the HPV16/31 genome and E6/E7 oncogenes (43). E6 and E7 serve as two key oncogenic transcription factors in the early stages of the HPV16 viral life cycle. E6 serves an important role in p53 degradation and can eliminate cell growth arrest, and E7 binds to and inactivates retinoblastoma protein (pRb), thus interfering with cell cycle regulation (44). E6 and E7 can synergistically promote the proliferation of infected cells, which promotes the occurrence and development of cervical cancer (35). By degrading p53, pRb and other tumor suppressor genes, these transformed keratinocytes have important effects on the expression of numerous target genes in cells (45).

The suppression of STAT1 by HPV proteins is necessary for the stable maintenance of viral episomes and genome amplification (11). By contrast, the present study revealed that STAT1 expression was positively correlated with HPV16 viral load in human cervical tissue samples. The reason for this relationship between STAT1 expression and HPV16 viral load is currently unclear, but three possible reasons were suggested. First, the cervical tissue consists of a complex microenvironment, including squamous epithelial cells, keratinocytes, macrophages, dendritic cells (DCs), neutrophils and natural killer (NK) cells. HPV16 infection is an ongoing process. Macrophages and DCs can recognize and phagocytize HPV16 viruses, induce an inflammatory response, and secrete pro-inflammatory cytokines, including IFN-γ, interleukin-1β and TNF-α. IFN-γ can increase STAT1 expression, even in HPV-infected cells (11). Second, a previous study revealed that STAT3, a regulator of proliferation, apoptosis, inflammation and tumorigenesis, is positively associated with HPV16 and promotes cervical cancer progression (46). HPV16 E6 has also been demonstrated to increase and activate STAT3 (47). In addition, STAT3 has been reported to be able to cooperate with epidermal growth factor receptor and HER2 to regulate STAT1 transcriptional activation in breast cancer cells (48). Therefore, it was speculated that HPV16 mediates STAT1 transcription by activating STAT3 in cervical cancer. Third, in the case of HPV-positive oropharyngeal squamous cell carcinoma with non-destructive p53 mutations, the E6/E7 oncoprotein has been reported to be less likely to control intracellular target genes, including STAT1, by inhibiting p53 expression (49). It may be hypothesized that STAT1 expression may be associated with the type of p53 mutations inhibited by E6 in cervical cancer tissues.

The present bioinformatics analysis demonstrated that high STAT1 expression was associated with improved survival in patients with cervical cancer. However, the biological mechanisms behind this remain to be determined. Accumulating evidence has demonstrated that increased STAT1 expression promotes tumor progression in multiple types of cancer, including breast cancer (29), and STAT1 has been demonstrated to be associated with poor survival, as well as chemotherapy and radiotherapy resistance, in breast tumors (50). By contrast, previous studies have also demonstrated that STAT1, as a tumor suppressor (51), restrains angiogenesis and cell proliferation (52). Furthermore, increased STAT1 expression is associated with improved disease-free survival and overall survival in ovarian cancer (34), and could induce the expression of pro-survival genes and resistance to genotoxic stress in human squamous carcinoma cells (53). In addition, STAT1 serves important roles in determining the immune responses in the tumor microenvironment, and clinical findings have demonstrated the importance of immune responses in suppressing tumorigenesis (54). STAT1 can manipulate tumor immunoregulation by elevating the cytotoxicity of NK cells and CD8+ T cells and influencing the surface receptors on tumor cells, including programmed death-ligand 1 and cytotoxic T-lymphocyte associated protein 4, which leads to the elimination of tumor cells by the immune system (32).

STAT1 can regulate various cellular processes, such as antimicrobial activities, tumor suppression, cell proliferation and cell death (32). Previous studies have demonstrated that STAT1 could mediate chemo- and radio-resistance during anticancer therapy (55,56), while a recent cervical cancer study suggested that STAT1 may contribute to the enhanced radiosensitivity by regulating intracellular levels of reactive oxygen species (57). The present study revealed that the patients with high STAT1 expression had improved survival. STAT1 treatment may increase the radiosensitivity to improve the prognosis of patients with cervical cancer treated by radiotherapy, and improve the survival of patients with cervical cancer who had low STAT1 expression.

In summary, the present study demonstrated that STAT1 expression was positively associated with the progression of cervical lesions in HPV16-infected patients. The combined detection of STAT1 expression and HPV16 viral load may provide an improved method for distinguishing the grade of cervical lesions, and STAT1 expression was demonstrated to be associated with a favorable prognosis in cervical cancer. These findings will help to further understand the development of cervical lesions and may provide improved methods for tumor therapy.

Acknowledgements

Not applicable.

Funding

The present study was supported by the National Natural Science Foundation of China (grant no. 81171581), Natural Science Foundation of Liaoning Province (grant no. LFWK201710) and 345 Talent Project of Shengjing Hospital, China Medical University.

Availability of data and materials

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

Authors' contributions

SWu, YW and ZS were responsible for the conception and design of the study. SWu, YW and YZ collected the human samples. SWu, YW, YY, ML, YL, CC, SWa and MY performed experiments. SWu and YW analyzed results and interpreted data. ZS supervised the study. SWu and YL drafted the manuscript. YL, CC, MY, YZ and ZS revised the manuscript. All authors read and approved the final manuscript.

Ethics approval and consent to participate

The Ethics Committee of Shengjing Hospital, China Medical University approved the present study (approval no. 2019PS476K). Written informed consent was obtained from all patients.

Patient consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

References