Identification of POLQ as a key gene in cervical cancer progression using integrated bioinformatics analysis and experimental validation

  • Authors:
    • Yuqin Zang
    • Ruqian Zhao
    • Tao Wang
    • Yueqian Gao
    • Lingli Chen
    • Shiqi Liu
    • Yingmei Wang
    • Fengxia Xue
  • View Affiliations

  • Published online on: April 26, 2023     https://doi.org/10.3892/mmr.2023.13002
  • Article Number: 115
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

As the most common gynecologic malignancy worldwide, cervical cancer (CC) is a serious hazard to health. Therefore, the present study aimed to identify the key genes in CC progression using integrated bioinformatics analysis and experimental validation. The mRNA microarray GSE63514 and microRNA (miRNA) microarray GSE86100 were obtained from the Gene Expression Omnibus database, and the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) in the progression of CC were identified. Thereafter, GO and KEGG functional enrichment analysis, protein‑protein interaction (PPI) network and significant subnetworks construction, and miRNA‑target regulatory network construction were performed. Based on the results of integrated bioinformatics analysis, the DEGs structural maintenance of chromosomes 4 (SMC4), ATPase family, AAA domain‑containing 2 (ATAD2) and DNA polymerase θ (POLQ) were identified as hub genes in the PPI network and were involved in the first significant subnetwork. In addition, these DEGs were predicted to be regulated by miR‑106B, miR‑17‑5P, miR‑20A and miR‑20B, which were identified as DEMs. Of note, SMC4 and ATAD2 are tumor‑promotors in CC. In the present study, small interfering (si)RNAs were used to knock down POLQ expression. Cell Counting Kit‑8, Transwell, cell cycle and apoptosis analyses revealed that the downregulation of POLQ restrained cell proliferation, migration and invasion, and promoted apoptosis and the arrest of the cell cycle in the G2 phase. In conclusion, POLQ, which may have a close interaction with SMC4 and ATAD2, may serve a vital role in the progression of CC.

Introduction

Cervical cancer (CC) is a serious hazard to health worldwide; it has the highest incidence and mortality rates among gynecological cancers (1). In 2020, >6 million people were diagnosed with CC and >3.4 million people with CC succumbed globally (1). The majority of cases of CC are caused by persistent infection with high-risk human papillomavirus (HPV), a virus comprising a double-stranded circular DNA that can incorporate into the host genome (2). Following the incorporation of HPV DNA, the E2 protein is deactivated, causing the dysregulated expression of the E6 and E7 oncoproteins, thereby driving the development of the cervical epithelium from normal to cervical intraepithelial neoplasia (CIN) and finally to CC (2,3). However, fully understanding the molecular mechanisms involved in CC pathogenesis is a distant achievement, and new biomarkers and therapeutic targets need to be explored.

With the development of gene sequencing technology, the abundant genome sequencing data of human diseases have been uploaded to public databases, such as the Gene Expression Omnibus database (GEO; http://www.ncbi.nlm.nih.gov/geo), and bioinformatics analysis has become a common method for assessing the molecular change in these diseases. The expression level of RNAs, such as mRNA and microRNA (miRNA), between different samples can be compared conveniently, and further information can also be obtained. In the present study, integrated bioinformatics analysis methods were used to identify differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) in the process of CC progression and to analyze the function, interaction and signaling pathways. Thereafter, a key gene network that probably has a significant influence on the progression of CC was screened out and verified in vitro. Overall, the present study may provide new insights into the molecular mechanism of CC progression and the search for new therapeutic targets.

Materials and methods

Selection of microarrays

The microarrays were obtained from the GEO database. The mRNA microarray GSE63514 is comprised of 24 normal cervical epithelium samples, 76 CIN samples and 28 CC samples, all of which have HPV infection and >90% of which have high-risk HPV infection. The miRNA microarray GSE86100 contains six samples of normal mucosa tissues with no HPV infection and six samples of CC tissues with high-risk HPV infection.

Identification of DEGs and DEMs

The present study analyzed the gene expression profiles using the analysis tool GEO2R (ncbi.nlm.nih.gov/geo/geo2r) on the GEO website to obtain putative DEGs and DEMs between different cervical samples. Normal cervical epithelium vs. CIN samples, and CIN vs. CC samples were compared to determine the DEGs. The overlapping genes exhibiting similar expression trends from the two comparisons were regarded as DEGs that were differentially expressed in both processes from normal to precancerous, and from precancerous to cancerous. Data were compared between normal tissue samples and CC tissue samples to determine the DEMs. P<0.05 and |Log FC|≥1, where FC is fold change, were set as the cut-off values.

Functional enrichment analyses of DEGs

Gene Ontology (GO) functional term enrichment, comprising biological processes, cellular components and molecular functions, as well as Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed using the DAVID 6.8 (https://david.ncifcrf.gov) website.

Protein-protein interaction (PPI) network

A PPI network is a visual display of the molecular interaction among proteins encoded by DEGs. Interaction information was exported from the STRING database 10.0 (string-db.org); a combined score of >0.4 was set as the threshold value. Cytoscape 3.5.0 (http://www.cytoscape.org) software for visualizing complex networks, was used for network construction, and subnetwork analysis was performed using the plug-in MCODE (http://apps.cytoscape.org/apps/mcode). The node degrees were analyzed using the Cytoscape NetworkAnalyzer. Further functional enrichment analysis of genes in the extracted subnetworks was conducted using the DAVID 6.8 website.

miRNA-target regulatory network

The online tool WebGestalt (http://www.webgestalt.org) was used to conduct miRNA-target gene enrichment analysis. By uploading DEGs to the website, miRNA-target interaction information can be obtained immediately, and a miRNA-target regulatory network can be generated using Cytoscape.

Literature search

A literature search was performed using terms structural maintenance of chromosomes 4 (‘SMC4’), ATPase family, AAA domain-containing 2 (‘ATAD2’) and DNA polymerase θ (‘POLQ’), combined with ‘cervical cancer’ on PubMed (pubmed.ncbi.nlm.nih.gov/) up to 20th October 2022.

Cell culture and transfection

According to the data from mRNA microarray GSE63514, all the samples are infected by HPV, which meant the DEGs identified may be associated with HPV infection. Therefore, the present study selected HeLa cells, a type of human CC cell line with HPV infection, to conduct the experiment validation. HeLa cells were purchased from the American Type Culture Collection and provided by the Tianjin Medical University. Cells were cultured and maintained in DMEM supplemented with 10% fetal bovine serum (FBS), 100 U/ml penicillin and 100 U/ml streptomycin in a humidified incubator at 37°C with 5% CO2. Three small interfering (si)RNAs targeting DNA polymerase θ (POLQ; si-POLQ) and negative control siRNA (si-NC) were synthesized by Shanghai GenePharma Co., Ltd.; the siRNA sequences are listed in Table I. Following the manufacturer's instruction, all siRNAs were transfected into the HeLa cells (density of 50%) at 16.7 nM using Lipofectamine® 3000 (Thermo Fisher Scientific, Inc.) at 37°C for 24 or 48 h. After transfection, the efficiency of transfection was evaluated by measuring mRNA levels at 24 h and protein levels at 48 h. Transfected cells were collected for subsequent experiments 48 h after transfection.

Table I.

siRNA sequences used in transfection experiments.

Table I.

siRNA sequences used in transfection experiments.

siRNASequence (sense)Sequence (antisense)
si-POLQ#1 5′-GCAAAGGCCUACUUCCCAUTT-3′ 5′-AUGGGAAGUAGGCCUUUGCTT-3′
si-POLQ#2 5′-GGAUAUAUUUCCUGUCCAATT-3′ 5′-UUGGACAGGAAAUAUAUCCTT-3′
si-POLQ#3 5′-CCAGAUACACAGGGAUUAATT-3′ 5′-UUAAUCCCUGUGUAUCUGGTT-3′
si-NC 5′-UUCUCCGAACGUGUCACGUTT-3′ 5′-ACGUGACACGUUCGGAGAATT-3′

[i] NC, negative control; POLQ, DNA polymerase θ; si, small interfering RNA.

Reverse transcription-quantitative PCR (RT-qPCR)

Following transfection with siRNA for 24 h, the total RNA of HeLa cells (with a density of 70%) was obtained using TRIzol® reagent (Thermo Fisher Scientific, Inc.) and subsequently reverse transcribed using the GoScript Reverse Transcription System (Promega Corporation) according to the manufacturer's protocols. qPCR was performed using the GoTaq qPCR Master Mix (Promega Corporation) and a 7500 Real-Time PCR Instrument (Thermo Fisher Scientific, Inc.) according to the manufacturer's instructions. The PCR cycling conditions were as follows: an initial denaturation step at 95°C for 2 min, followed by 40 cycles of 95°C for 15 sec and an annealing and extension step at 60°C for 1 min. The 2−ΔΔCq method was used for data analysis (4); gene expression levels were normalized to that of the housekeeping gene, GAPDH. PCR primers were synthesized by Tianyi Huiyuan Biotechnology Co. Ltd. and the sequences were: POLQ forward, 5′-TCCTACACCCATTCCAACATCTG-3′ and reverse, 5′-GTCTTTGAACCCATTTCTACTCCC-3′; GAPDH forward, 5′-GGTGGTCTCCTCTGACTTCAACA-3′ and reverse, 5′-GTTGCTGTAGCCAAATTCGTTGT-3′. These experiments were replicated three times.

Western blotting

Following transfection with siRNA for 48 h, total proteins were extracted from HeLa cells using RIPA lysis buffer and determined with the BCA Protein Assay Kit (Solarbio). 20 µg of protein were subjected to SDS-PAGE separation (10% separating gel and 5% stacking gel) and transferred onto PVDF membranes. Following blocking with 5% nonfat milk at room temperature for 2 h, the membranes were incubated with primary antibodies against POLQ (cat. No. orb48495; 1:1,000 dilution; Biorbyt, Ltd.) and GAPDH (cat. No. 97166; 1:1,000 dilution; CST) overnight at 4°C. The membranes were then incubated with the appropriate secondary antibody (cat. No. ZB-5301; and cat. No. ZB-5305; both 1:5,000; both ZSGB-BIO) for 1 h. Protein bands were visualized using the Immobilon Western HRP Substrate (MilliporeSigma) and the results were analyzed using the software Image J Version 2.0.0-rc-54 (imagej.net/software/fiji/).

Cell Counting Kit (CCK)-8 assay

CCK−8 (Beyotime Institute of Biotechnology) was used to assess cell viability and proliferation. HeLa cells were seeded in 96-well plates and transfected with si-NC, si-POLQ#2 or si-POLQ#3. After 1, 2, 3, 4 and 5 days, 10 µl of CCK-8 solution was added to each well. After 1.5 h of culture, the OD value at 450 nm was determined.

Transwell assay

The cell migration assay was conducted using 24-well Transwell chambers (Corning, Inc.). The transfected cells were suspended in serum-free medium, and 1×105 cells were seeded into the upper chamber; medium supplemented with 10% FBS was added to the lower compartment. After 24 h of culture, the migrated cells were fixed with 4% paraformaldehyde for 15 min at room temperature, stained with crystal violet ammonium oxalate solution for 10 sec at room temperature. A total of six fields for each chamber were examined with a light microscope at 100× magnification for cell counting. For the invasion assay, prior to the cell suspension being prepared, 70 µl of Matrigel was added to each upper Transwell chamber at 37°C for 1 h. The subsequent procedures were the same as those used for the migration assay.

Apoptosis and cell cycle analysis

HeLa cells transfected with siRNAs were collected with 25% trypsin without EDTA and washed with PBS. The annexin V-FITC/PI Apoptosis kit (BD Biosciences) was used for apoptosis detection. The cells were resuspended in 1X binding buffer and adjusted to 1×106 cells/ml. Thereafter, 5 µl of annexin V-FITC and 5 µl propidium iodide (PI) were added to 100 µl of cell suspension. After incubation for 15 min in the dark, 400 µl of the 1X binding buffer was added to the cell suspension. Apoptosis (early + late apoptotic cells) was detected within 1 h using a flow cytometer (Beckman Coulter, Inc.). For cell cycle analysis, HeLa cells were fixed with 75% cold ethanol, resuspended in 100 µl RNase A and incubated at 37°C for 30 min. Subsequently, cells were stained with 400 µl PI buffer and incubated in the dark at 4°C for 30 min. The DNA content was determined using flow cytometry.

Statistical analysis

The data were analyzed using one-way ANOVA followed by Dunnett's multiple comparisons test with GraphPad Prism 8.0 (Dotmatics). P<0.05 was considered to indicate a statistically significant difference.

Results

DEGs and DEMs

Based on the GSE63514 dataset, a total of 1,318 DEGs were screened upon comparison of the CIN sample data with those of normal cervical epithelium samples. Furthermore, 1,575 DEGs were screened between CC and CIN tissues. The expression heat maps of the top 100 up- and downregulated DEGs are shown in Fig. 1A and B. Among them, 250 genes, including 118 upregulated and 132 downregulated genes (data not shown), were common to both comparisons and presented the same expression trends (Fig. 1C). Based on the miRNA dataset, GSE86100, 166 DEMs were screened, including 74 upregulated and 92 downregulated DEMs, and the top 50 are shown in Fig. 1D.

GO and KEGG enrichment analyses

Using the online tool, DAVID, the GO term enrichments in terms of biological processes, cellular components and molecular functions, and KEGG pathway enrichment were analyzed; the top five significant enrichment terms for each category are presented in Fig. 2. For the biological processes, the upregulated DEGs were mainly enriched in ‘nuclear division’ and ‘mismatch repair’ pathway (Fig. 2A), whereas the downregulated DEGs were mainly enriched in the process of ‘icosanoid metabolic process’ and ‘arachidonic acid metabolism’ pathway (Fig. 2B).

PPI network

As depicted in Fig. 3A, the PPI network consisted of 123 nodes (including 67 nodes of upregulated DEGs and 56 nodes of downregulated DEGs) and 283 edges. Based on the NetworkAnalyzer tool, the top 15 hub proteins were RFC4, SMC4, STAT1, KNTC1, ATAD2, FBXO5, TRIP13, ESR1, CD44, CENPN, ANLN, CDK2, ECT2, KIF14 and POLQ (Table II). Using the plug-in MCODE, three significant subnetworks with scores >4 were extracted (Fig. 3B-D). For the first significant subnetwork, it contains 72 edges and 13 nodes, while 11 of these nodes belong to the top 15 hub proteins of the PPI network. For the second significant subnetwork, it contains 33 edges and 9 nodes, including IFIT3, USP18, IFI44L, STAT1, GBP5, IFI44, EPSTI1, RSAD2, GBP1. And for the third significant subnetwork, it is composed of 20 edges and 9 nodes, and contains the only two downregulated DEGs (ESR1 and AR) in these three significant subnetworks. The top five biological processes in the GO analysis of each subnetwork are presented in Fig. 3E-G.

Table II.

Top 15 hub genes identified in the protein-protein interaction network.

Table II.

Top 15 hub genes identified in the protein-protein interaction network.

Hub geneGene nameNode degree
RFC4Replication factor C subunit 419
SMC4Structural maintenance of chromosomes 418
STAT1Signal transducer and activator of transcription 117
KNTC1Kinetochore associated 117
ATAD2ATPase family, AAA domain-containing 216
FBXO5F-box protein 516
TRIP13Thyroid hormone receptor interactor 1315
ESR1Estrogen receptor 114
CD44CD44 molecule (Indian blood group)13
CENPNCentromere protein N13
ANLNAnillin actin-binding protein13
CDK2Cyclin dependent kinase 212
ECT2Epithelial cell transforming 212
KIF14Kinesin family member 1412
POLQDNA polymerase θ12
miRNA-target regulatory network

The miRNA-target regulatory network is presented in Fig. 4A. This network comprises 66 nodes, including 30 DEGs and 36 miRNAs, and 137 edges. The majority of DEGs in this regulatory network were upregulated and were included in the PPI network. It was found that approximately one-third of the miRNAs in the regulatory network were also the screened DEMs, aforementioned. Collectively, SMC4, ATAD2 and POLQ in the miRNA-target regulatory network were not only identified as hub genes in the PPI network, but are also the nodes of the first significant subnetwork (Table II and Fig. 3B). Furthermore, these genes were modulated by miR-106b, miR-17-5p, miR-20a and miR-20b, which were identified as DEMs in the CC samples (Figs. 4B and 1D). Therefore, it was hypothesized that the three DEGs (SMC4, ATAD2 and POLQ) may interact closely with each other in CC progression.

Literature search of SMC4 and ATAD2 in CC

Using ‘SMC4’, ‘ATAD2’ and ‘POLQ’, combined with ‘cervical cancer’, a total of six articles on the roles of SMC4 and ATAD2 in CC were obtained from PubMed (510). As in Table III, SMC4 expression was upregulated in CC and recognized as a hub gene of miRNA-transcription factor-mRNA network of HPV-positive CC (57). Moreover, SMC4 can negatively regulate CTCF which mediates the ribosomal RNA gene transcription in CC, and promote proliferation, migration, invasion, epithelial-mesenchymal transition, spheroid formation of CC cells (7,8). ATAD2 is also identified as one hub gene of the PPI network and involved in the circRNA-miRNA-mRNA network in CC (9). Besides, ATAD2 is overexpressed in CC tissues and related to adverse clinicopathologic features of CC patients, and the knockdown of ATAD2 can reduce the growth, invasion, migration and clonogenic potential of CC cells (10). To the best of the authors’ knowledge, the involvement of POLQ in CC has not yet been reported. As a result, the present study proceeded to validate the effect of POLQ on the biological behavior of HeLa cells.

Table III.

Studies of SMC4 and ATAD2 in CC.

Table III.

Studies of SMC4 and ATAD2 in CC.

Author, yearGeneResult(Refs.)
He, 2021SMC4SMC4 expression levels are elevated in CC tissues compared with adjacent normal tissues.(7)
SMC4 promotesproliferation, migration, invasion, epithelial-mesenchymal transition and spheroid formation of CC cells, as well as the expression of stem cell markers.
S, 2021 SMC4 is one of top 10 dysregulated hub proteins in CC.(6)
Yuan, 2020 SMC4 expression is upregulated and is as a hub gene in the miRNA-transcription factor- mRNA network of HPV-positive CC.(5)
Huang, 2013 SMC4 negatively regulates CTCF, which mediates ribosomal RNA gene transcription in CC.(8)
Yi, 2019ATAD2ATAD2 is identified as one of the PPI network hub genes and is involved in the circRNA-miRNA-mRNA network in the pathology of CC.(9)
Zheng, 2015 ATAD2 expression is highly elevated in CC tissue and related to cancer stage, lymph metastasis and worse patient survival.(10)
ATAD2 knockdown reduces the proliferation, invasion, migration and clonogenic potential of CC cells.

[i] ATAD, ATPase family, AAA domain-containing 2; CC, cervical cancer; circRNA, circular RNA; CTCF, CCCTC-binding factor; HPV, human papillomavirus; miRNA, microRNA; PPI, protein-protein interaction; SMC4, structural maintenance of chromosomes 4.

Effects of POLQ knockdown on CC cells

siRNAs were used to knock down the expression of POLQ in HeLa cells. The efficiency of three POLQ siRNAs was assessed using western blotting (Fig. 5A and B) and RT-qPCR (Fig. 5C). si-POLQ#2 and si-POLQ#3 exhibited greater efficiency and were selected for subsequent experimental validation. Cell proliferative ability was determined using the CCK-8 test, and the data showed that HeLa cells transfected with si-POLQ#2 and si-POLQ#3 had lower absorbance values compared with those transfected with si-NC (Fig. 5D. Thus, POLQ may be important for the proliferation of HeLa cells. To determine the effect of POLQ knockdown in migration and invasion, Transwell assays were performed and showed that the knockdown of POLQ inhibited the migration and invasion of HeLa cells (Fig. 5E-G). Apoptosis and cell cycle progression were assessed using flow cytometry. The percentage of apoptotic cells increased when POLQ was downregulated compared with that in the si-NC group (Fig. 5H and I). Furthermore, the percentage of cells arresting in the G2 phase were increased in the si-POLQ#2 group and si-POLQ#3 group than in the si-NC group (Fig. 5J and K), indicating that POLQ knockdown may impede HeLa cell cycle progression.

Discussion

CC is the most common gynecological malignancy and causes the most mortality associated with gynecological malignancy worldwide (1). To explore the molecular mechanism of CC progression, bioinformatics analysis among normal cervical epithelium samples, CIN samples and CC samples was performed based on the microarray datasets of GSE63514 and GSE86100. A total of 118 upregulated and 132 downregulated DEGs were identified in the process of progression from normal cervical epithelium to CIN and from CIN to CC. In addition, 166 DEMs were identified based on an analysis of the normal cervical epithelium and CC data. Based on the results of GO term and KEGG pathway enrichment analyses, PPI network analysis and miRNA-target regulatory network construction, the upregulated DEGs, SMC4, ATAD2 and POLQ, were identified as hub genes in the PPI network and as genes involved in the first significant subnetwork. Notably, in the miRNA-target regulatory network, these three DEGs were all regulated by miR-106B, miR-17-5P, miR-20A and miR-20B, which were identified as DEMs in the present study and which were shown to be modulated by the HPV proteins E6/E7 in a previous study (1114). Therefore, it was hypothesized that SMC4, ATAD2 and POLQ may interact closely with each other and serve vital roles in CC progression.

SMC4 serves a role in chromosome condensation and segregation (1517). As SMC4 contains an ATPase domain, its inactivation can cause cell death (15). The overexpression and tumor-promoting function of SMC4 in certain types of cancer, such as pancreatic ductal adenocarcinoma (16), endometrial cancer (17), lung adenocarcinoma (18) and glioma (19), have been reported in a number of studies. SMC4 has been identified as a crucial gene in CC development in other studies (5,6). Through in vitro experiments, knockdown of SMC4 has been shown to inhibit cell proliferation, migration and invasion, as well as colony formation and EMT in CC (7). In addition, the expression of stem cell markers, spheroid formation and activation of NF-κB pathways were reduced (7). In addition, SMC4 can negatively regulate CCCTC-binding factor (CTCF)-enhanced rRNA gene transcription (8).

As a member of the ATPase family, ATAD2 participates in chromatin dynamics, DNA replication and gene transcription (20). The overexpression of ATAD2 is detected in numerous advanced types of human cancer, including pancreatic cancer (21), gastric cancer (22), ovarian cancer (23), colorectal cancer (24) and esophageal carcinoma (25). Furthermore, ATAD2 is regarded as an epigenetic reader, transcriptional co-regulator, proliferation driver, apoptosis negative regulator and cell cycle progression inhibitor in various types of cancer as it regulates signaling pathways, such as MAPK, AKT, HIF1α, hedgehog and Rb/E2F-cMyc (20,26,27). The potential and prospects of ATAD2 as a drug target for cancer have been identified (20,26). Consistent with the present study, ATAD2 has been identified as a hub gene and regulated by miR-106B in CC in a study by Yi et al (9). Furthermore, the regulatory effect of miR-106B on ATAD2 was previously verified in papillary thyroid cancer by Sun et al (28). In addition, according to Zheng et al (10), ATAD2 may be important in the promotion of growth and metastasis in CC.

POLQ exhibits ATPase activity similar to ATAD2 and SMC4 (2933). POLQ has been revealed to serve a vital role in DNA replication and the repair of DNA double-strand breaks induced by different factors, such as endonucleases and Cas9 nickase (29). Furthermore, POLQ expression is upregulated in various types of cancer, including hepatocellular carcinoma (30), lung adenocarcinoma (31), breast cancer (32) and colorectal cancer (33), and it is associated with the poor prognosis (3234). POLQ is necessary for cell proliferation and migration, and inhibits apoptosis (30), and it is recognized as a promising drug target for the treatment of cancer (35). In addition, POLQ can be downregulated by miR-20B in XP-V tumor cells, which accelerates mismatch in DNA replication repairing (36). However, the role of POLQ in CC remains to be elucidated. Therefore, upon identification of POLQ as one of the key genes in the progression of CC, the present study proceeded to perform experimental validation in vitro, using si-POLQ transected HeLa cells. Results from the present study showed that proliferation, migration and invasion of HeLa cells were inhibited by POLQ knock down. In addition, apoptosis was increased and cell cycle was impeded at the G2 phase when POLQ was downregulated. These results suggested that POLQ may be important in the progression of CC.

In conclusion, based on the results of integrated bioinformatics analyses, the gene network comprising SMC4, ATAD2 and POLQ, which are all genes associated with ATPase activity and related to the genome regulation, may serve important roles in the progression of CC. Furthermore, the effect of these genes may be modulated by the HPV proteins E6/E7 and their regulatory miRNAs (miR-106B, miR-17-5P, miR-20A and miR-20B). To the best of our knowledge, this is the first study to elucidate the potential role of this gene network and to verify the promoting effect of POLQ in CC progression through in vitro experiments. Such findings may provide novel insights into the progression mechanism and novel therapeutic targets of CC. However, the present study has some limitations. The sample size of the miRNA microarray GSE86100 dataset was not large enough, and differences in terms of HPV infection were found between samples in the two microarrays, which may affect the accuracy of the results. The relationships between HPV infection and these DEGs were not well demonstrated and further studies are needed.

Acknowledgements

Not applicable.

Funding

The present study was supported by The National Natural Science Foundation of China (grant no. 81972448) and The Science and Technology Project of Tianjin, China (grant no. 20JCZDJC00330).

Availability of data and materials

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

Authors' contributions

YZ was responsible for conceiving the study, methodology, data analysis and writing the manuscript. RZ was responsible for the methodology and writing the manuscript. TW was responsible for data analysis and manuscript revision. YG performed experiments. LC and SL were responsible for data curation. YW and FX conceived the study, reviewed and edited the final manuscript, and were responsible for funding acquisition. YZ and RZ confirm the authenticity of all the raw data. All authors read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.

Patient consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

References

1 

Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A and Bray F: Global cancer statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 71:209–249. 2021. View Article : Google Scholar : PubMed/NCBI

2 

Wang T, Gong M, Cao Y, Zhao C, Lu Y, Zhou Y, Yao S, Chen J, Zhao C and Ju R: Persistent ferroptosis promotes cervical squamous intraepithelial lesion development and oncogenesis by regulating KRAS expression in patients with high risk-HPV infection. Cell Death Discov. 8:2012022. View Article : Google Scholar : PubMed/NCBI

3 

Sun J, Yue Y, Li R, Sun Q, Hu C, Ge X and Guan Q: Detection of HPV E6/E7 mRNA in the diagnosis of cervical cancer and precancerous lesions after kidney transplantation. Am J Transl Res. 13:7312–7317. 2021.PubMed/NCBI

4 

Schmittgen TD and Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 3:1101–1108. 2008. View Article : Google Scholar : PubMed/NCBI

5 

Yuan Y, Shi X, Li B, Peng M, Zhu T, Lv G, Liu L, Jin H, Li L and Qin D: Integrated analysis of key microRNAs /TFs /mRNAs/ in HPV-positive cervical cancer based on microRNA sequencing and bioinformatics analysis. Pathol Res Pract. 216:1529522020. View Article : Google Scholar : PubMed/NCBI

6 

S D A, . Pasumarthi D, Pasha A, Doneti R B S, Botlagunta M, B VL and Pawar SC: Identification of differentially expressed genes in cervical cancer patients by comparative transcriptome analysis. Biomed Res Int. 2021:88100742021.PubMed/NCBI

7 

He H, Zheng C and Tang Y: Overexpression of SMC4 predicts a poor prognosis in cervical cancer and facilitates cancer cell malignancy phenotype by activating NF-κB pathway. Hum Cell. 34:1888–1898. 2021. View Article : Google Scholar : PubMed/NCBI

8 

Huang K, Jia J, Wu C, Yao M, Li M, Jin J, Jiang C, Cai Y, Pei D, Pan G and Yao H: Ribosomal RNA gene transcription mediated by the master genome regulator protein CCCTC-binding factor (CTCF) is negatively regulated by the condensin complex. J Biol Chem. 288:26067–26077. 2013. View Article : Google Scholar : PubMed/NCBI

9 

Yi Y, Liu Y, Wu W, Wu K and Zhang W: Reconstruction and analysis of circRNA-miRNA-mRNA network in the pathology of cervical cancer. Oncol Rep. 41:2209–2225. 2019.PubMed/NCBI

10 

Zheng L, Li T, Zhang Y, Guo Y, Yao J, Dou L and Guo K: Oncogene ATAD2 promotes cell proliferation, invasion and migration in cervical cancer. Oncol Rep. 33:2337–2344. 2015. View Article : Google Scholar : PubMed/NCBI

11 

Liu X: Up-regulation of miR-20a by HPV16 E6 exerts growth-promoting effects by targeting PDCD6 in cervical carcinoma cells. Biomed Pharmacother. 102:996–1002. 2018. View Article : Google Scholar : PubMed/NCBI

12 

Zheng ZM and Wang X: Regulation of cellular miRNA expression by human papillomaviruses. Biochim Biophys Acta. 1809:668–677. 2011. View Article : Google Scholar : PubMed/NCBI

13 

Honegger A, Schilling D, Bastian S, Sponagel J, Kuryshev V, Sültmann H, Scheffner M, Hoppe-Seyler K and Hoppe-Seyler F: Dependence of intracellular and exosomal microRNAs on viral E6/E7 oncogene expression in HPV-positive tumor cells. PLoS Pathog. 11:e10047122015. View Article : Google Scholar : PubMed/NCBI

14 

Cheng Y, Geng L, Zhao L, Zuo P and Wang J: Human papillomavirus E6-regulated microRNA-20b promotes invasion in cervical cancer by targeting tissue inhibitor of metalloproteinase 2. Mol Med Rep. 16:5464–5470. 2017. View Article : Google Scholar : PubMed/NCBI

15 

Palou R, Dhanaraman T, Marrakchi R, Pascariu M, Tyers M and D'Amours D: Condensin ATPase motifs contribute differentially to the maintenance of chromosome morphology and genome stability. PLoS Biol. 16:e20039802018. View Article : Google Scholar : PubMed/NCBI

16 

Fukuhisa H, Seki N, Idichi T, Kurahara H, Yamada Y, Toda H, Kita Y, Kawasaki Y, Tanoue K, Mataki Y, et al: Gene regulation by antitumor miR-130b-5p in pancreatic ductal adenocarcinoma: The clinical significance of oncogenic EPS8. J Hum Genet. 64:521–534. 2019. View Article : Google Scholar : PubMed/NCBI

17 

Yan Y, Liu C, Zhang J, Li W, Yin X, Dong L, Pang S and Li X: SMC4 knockdown inhibits malignant biological behaviors of endometrial cancer cells by regulation of FoxO1 activity. Arch Biochem Biophys. 712:1090262021. View Article : Google Scholar : PubMed/NCBI

18 

Zhang C, Kuang M, Li M, Feng L, Zhang K and Cheng S: SMC4, which is essentially involved in lung development, is associated with lung adenocarcinoma progression. Sci Rep. 6:345082016. View Article : Google Scholar : PubMed/NCBI

19 

You A, Rao G, Wang J, Li J, Zhang Y, Gu J, Ge X, Zhang K, Gao X, Wu X, et al: MiR-433-3p restrains the proliferation, migration and invasion of glioma cells via targeting SMC4. Brain Res. 1767:1475632021. View Article : Google Scholar : PubMed/NCBI

20 

Nayak A, Dutta M and Roychowdhury A: Emerging oncogene ATAD2: Signaling cascades and therapeutic initiatives. Life Sci. 276:1193222021. View Article : Google Scholar : PubMed/NCBI

21 

Dutta M, Das B, Mohapatra D, Behera P, Senapati S and Roychowdhury A: MicroRNA-217 modulates pancreatic cancer progression via targeting ATAD2. Life Sci. 301:1205922022. View Article : Google Scholar : PubMed/NCBI

22 

Zhang MJ, Zhang CZ, Du WJ, Yang XZ and Chen ZP: ATAD2 is overexpressed in gastric cancer and serves as an independent poor prognostic biomarker. Clin Transl Oncol. 18:776–781. 2016. View Article : Google Scholar : PubMed/NCBI

23 

Liu Q, Liu H, Li L, Dong X, Ru X, Fan X, Wen T and Liu J: ATAD2 predicts poor outcomes in patients with ovarian cancer and is a marker of proliferation. Int J Oncol. 56:219–231. 2020.PubMed/NCBI

24 

Hou M, Huang R, Song Y, Feng D, Jiang Y and Liu M: ATAD2 overexpression is associated with progression and prognosis in colorectal cancer. Jpn J Clin Oncol. 46:222–227. 2016. View Article : Google Scholar : PubMed/NCBI

25 

Cao LJ, Zhang YJ, Dong SQ, Li XZ, Tong XT, Chen D, Wu ZY, Zheng XH, Xue WQ, Jia WH and Zhang JB: ATAD2 interacts with C/EBPβ to promote esophageal squamous cell carcinoma metastasis via TGF-β1/Smad3 signaling. J Exp Clin Cancer Res. 40:1092021. View Article : Google Scholar : PubMed/NCBI

26 

Hussain M, Zhou Y, Song Y, Hameed HMA, Jiang H, Tu Y and Zhang J: ATAD2 in cancer: A pharmacologically challenging but tractable target. Expert Opin Ther Targets. 22:85–96. 2018. View Article : Google Scholar : PubMed/NCBI

27 

Haitani T, Kobayashi M, Koyasu S, Akamatsu S, Suwa T, Onodera Y, Nam JM, Nguyen PTL, Menju T, Date H, et al: Proteolysis of a histone acetyl reader, ATAD2, induces chemoresistance of cancer cells under severe hypoxia by inhibiting cell cycle progression in S phase. Cancer Lett. 528:76–84. 2022. View Article : Google Scholar : PubMed/NCBI

28 

Sun W, Lan X, Zhang H, Wang Z, Dong W, He L, Zhang T, Zhang P, Liu J and Qin Y: NEAT1_2 functions as a competing endogenous RNA to regulate ATAD2 expression by sponging microRNA-106b-5p in papillary thyroid cancer. Cell Death Dis. 9:3802018. View Article : Google Scholar : PubMed/NCBI

29 

Wang Z, Song Y, Li S, Kurian S, Xiang R, Chiba T and Wu X: DNA polymerase θ (POLQ) is important for repair of DNA double-strand breaks caused by fork collapse. J Biol Chem. 294:3909–3919. 2019. View Article : Google Scholar : PubMed/NCBI

30 

Pan Q, Wang L, Liu Y, Li M, Zhang Y, Peng W, Deng T, Peng ML, Jiang JQ, Tang J, et al: Knockdown of POLQ interferes the development and progression of hepatocellular carcinoma through regulating cell proliferation, apoptosis and migration. Cancer Cell Int. 21:4822021. View Article : Google Scholar : PubMed/NCBI

31 

Shinmura K, Kato H, Kawanishi Y, Yoshimura K, Tsuchiya K, Takahara Y, Hosokawa S, Kawase A, Funai K and Sugimura H: POLQ overexpression is associated with an increased somatic mutation load and PLK4 overexpression in lung adenocarcinoma. Cancers (Basel). 11:7222019. View Article : Google Scholar : PubMed/NCBI

32 

Lemée F, Bergoglio V, Fernandez-Vidal A, Machado-Silva A, Pillaire MJ, Bieth A, Gentil C, Baker L, Martin AL, Leduc C, et al: DNA polymerase theta up-regulation is associated with poor survival in breast cancer, perturbs DNA replication, and promotes genetic instability. Proc Natl Acad Sci USA. 107:13390–13395. 2010. View Article : Google Scholar : PubMed/NCBI

33 

Pillaire MJ, Selves J, Gordien K, Gourraud PA, Gentil C, Danjoux M, Do C, Negre V, Bieth A, Guimbaud R, et al: A ‘DNA replication’ signature of progression and negative outcome in colorectal cancer. Oncogene. 29:876–887. 2010. View Article : Google Scholar : PubMed/NCBI

34 

Allera-Moreau C, Rouquette I, Lepage B, Oumouhou N, Walschaerts M, Leconte E, Schilling V, Gordien K, Brouchet L, Delisle MB, et al: DNA replication stress response involving PLK1, CDC6, POLQ, RAD51 and CLASPIN upregulation prognoses the outcome of early/mid-stage non-small cell lung cancer patients. Oncogenesis. 1:e302012. View Article : Google Scholar : PubMed/NCBI

35 

Yousefzadeh MJ and Wood RD: DNA polymerase POLQ and cellular defense against DNA damage. DNA Repair (Amst). 12:1–9. 2013. View Article : Google Scholar : PubMed/NCBI

36 

Guo J, Jiang Z, Li X, Wang X and Xiao Y: miR-20b downregulates polymerases κ and θ in XP-V tumor cells. Oncol Lett. 11:3790–3794. 2016. View Article : Google Scholar : PubMed/NCBI

Related Articles

Journal Cover

June-2023
Volume 27 Issue 6

Print ISSN: 1791-2997
Online ISSN:1791-3004

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Zang Y, Zhao R, Wang T, Gao Y, Chen L, Liu S, Wang Y and Xue F: Identification of <em>POLQ</em> as a key gene in cervical cancer progression using integrated bioinformatics analysis and experimental validation. Mol Med Rep 27: 115, 2023.
APA
Zang, Y., Zhao, R., Wang, T., Gao, Y., Chen, L., Liu, S. ... Xue, F. (2023). Identification of <em>POLQ</em> as a key gene in cervical cancer progression using integrated bioinformatics analysis and experimental validation. Molecular Medicine Reports, 27, 115. https://doi.org/10.3892/mmr.2023.13002
MLA
Zang, Y., Zhao, R., Wang, T., Gao, Y., Chen, L., Liu, S., Wang, Y., Xue, F."Identification of <em>POLQ</em> as a key gene in cervical cancer progression using integrated bioinformatics analysis and experimental validation". Molecular Medicine Reports 27.6 (2023): 115.
Chicago
Zang, Y., Zhao, R., Wang, T., Gao, Y., Chen, L., Liu, S., Wang, Y., Xue, F."Identification of <em>POLQ</em> as a key gene in cervical cancer progression using integrated bioinformatics analysis and experimental validation". Molecular Medicine Reports 27, no. 6 (2023): 115. https://doi.org/10.3892/mmr.2023.13002