Skip to main content

Observation of the cervical microbiome in the progression of cervical intraepithelial neoplasia

Abstract

Objective

Cervical microbial community in the cervical intraepithelial neoplasia and cervical cancer patients was analysed to study its composition, diversity and signalling pathways by high-throughput 16S rDNA sequencing,and the candidate genes associated with occurrence and progression of cervical intraepithelial neoplasia were screened out and the model was established to predict the evolution of cervical intraepithelial neoplasia malignant transformation from the cervical microbial genes aspect.

Methods

Cervical tissues of normal, cervical intraepithelial neoplasia and cervical cancer patients without receiving any treatment were collected. The correlation between candidate genes and cervical intraepithelial neoplasia progression was initially determined by analyzing the microbial flora. Real-time fluorescence quantitative PCR was used to detect the expression of candidate genes in different cervical tissues, ROC curve and logistic regression was used to analyse and predict the risk factors related to the occurrence and progression of cervical intraepithelial neoplasia. Finally, the early warning model of cervical intraepithelial neoplasia occurrence and progression is established.

Results

Cervical tissues from normal, cervical intraepithelial neoplasia and cervical cancer patients were collected for microbial community high-throughput 16S rDNA sequencing. The analysis revealed five different pathways related to cervical intraepithelial neoplasia. 10 candidate genes were selected by further bioinformatics analysis and preliminary screening. Real time PCR, ROC curve and Logistic regression analysis showed that human papillomavirus infection, TCT severity, ABCG2, TDG, PCNA were independent risk factors for cervical intraepithelial neoplasia. We used these indicators to establish a random forest model. Seven models were built through different combinations. The model 4 (ABCG2 + PCNA + TDG) was the best early warning model for the occurrence and progression of CIN.

Conclusions

A total of 5 differential pathways and 10 candidate genes related to occurrence and progression of cervical intraepithelial neoplasia were found in cervical microbial community. This study firstly identified the genes from cervical microbial community that play an important role in the occurrence and progression of cervical intraepithelial neoplasia. At the same time, the early warning model including ABCG2 + PCNA+TDG genes provided a new idea and target for clinical prediction and blocking the evolution of cervical intraepithelial neoplasia malignant transformation from the aspect of cervical microbiological related genes.

Peer Review reports

Background

Cervical intraepithelial neoplasia (CIN) is the precancerous lesion of cervical cancer. A lot of studies have demonstrated that human papilloma virus (HPV) is the main carcinogen responsible for CIN and cervical cancer. However, some other studies found that not all patients infected with HPV must become CIN or cervical cancer [1]. Although it is now believed that adjuvant factors other than HPV play a key role in the development of cancer, most of the potential mechanisms of this carcinogenic effect are still unknown [2]. It has been shown that the cause of human diseases is not only a single pathogen, but also involves the overall changes in the human microbiology group [3]. In recent years, with the application of metagenomic principles and the development of high-throughput sequencing analysis, research on the relationship between microorganisms and human diseases has been initiated. At present, the most common sequencing methods to identify the microbiome are pyrosequencing and 16S rRNA sequencing. MitraA [4] recently performed 16S rRNA gene amplification of the vaginal wall microorganisms of 52 cases of LSIL, 92 cases of HSIL, 5 cases of ICC, and 20 normal controls. The results indicated that vaginal microbial diversity is associated with the severity of CIN disease and that microbes can participate in regulating the persistence of viral infections and disease progression. The role of Cervical microbial community in the progression of CIN has gradually been recognized, and its synergy with HPV in CIN and cervical cancer is expected to become a hot spot in cervical disease research. Through random forest model, vaginal microbiome-derived bacterial markers can be used as a predictive model to predict the CIN malignant transformation, indicates that vaginal microbiome may play a role as biomarker. Current research on the correlation between CIN and vaginal microbes is mostly about the composition of the vaginal microbial flora and the pathogenic mechanism of the bacterial species, the severity of the bacterial flora and CIN, the relationship between the changes of the CIN flora and cervical cancer, and there are no studies involving genes and pathways. And no researchers have focused on the role of cervical tissue microecology and its related genes in the progression of CIN. Therefore, we analysed the cervical microbial community by high-throughput 16S rDNA sequencing, bioinformatics analysis and real time PCR to study its composition, diversity and signalling pathways in patients with CIN and cervical cancer. Finally, we screened out the candidate genes associated with occurrence and progression of CIN and established the best model to predict the evolution of CIN malignant transformation. Through this study, we proposed the important role of cervical microbial community and its related genes in the process of cervical cells carcinogenesis which is never discovered before.

Materials and Methods

Selection of study cases

Thirty-eight cases of CIN tissues (9, 11 and 18 cases were CIN1, 2 and 3 respectively) were randomly selected in Affiliated Cancer Hospital of Guangxi Medical University from May 2015 to July 2015; 14 cases of normal cervical tissue (taken from patients with uterine fibroids requiring hysterectomy) and 10 cervical cancer tissues were subjected to high-throughput 16S rDNA sequencing of microbial communities, and a series of analyses were performed. Fifty-two cases of CIN tissue, 38 cases of normal cervical tissue and 30 cases of cervical cancer tissue from January 2017 to December 2018 were again selected for realtimefluorescence quantitative PCR detection. All cases were discovered for the first time and confirmed by histopathology. No treatment was performed before the operation, and human papillomavirus detection and cervical cytology were performed before treatment. No trichomonas, Candida infection or bacterial vaginosis was detected in the vaginal secretions within 3 days before sampling. Subjects were required to abstain from sexual intercourse 3 days before sampling. No drugs affecting the vaginal flora were used before sampling.

Cervical sample collection and method

Samples were placed in a sterile tube for cryopreservation and immediately stored in liquid nitrogen. After the samples were collected, they were transferred to the laboratory for storage at − 80 °C. For the isolation, extraction and purification of total bacterial DNA, mechanical (magnetic bead repeated beating method, Fast Prep FP120) combined with enzymatic methods (QIAam DNA Mini Kit, QIAGEN, Valencia, CA, USA) were used to efficiently extract relevant microbial DNA. DNA samples are detected by fluorescence quantification and agarose gel electrophoresis. We collected 1 μL for fluorescence quantitative detection (instrument: QubitFluorometer,manufacturer: Thermo Fisher), and 5 μL DNA for electrophoresis detection (agarose gel electrophoresis utilizes 1% agarose gel electrophoresis at 150 V for approximately 40 min) of integrity and presence of RNA or protein and secondary metabolite contamination.

16SrDNA V4 region target fragment library construction

The total 30 ng of DNA was used as the template, and the V4 region of the bacterial 16S rDNA was used as the target. The universal primers fused with the Miseq platform sequencing platform were used for primer design and synthesis, and New England Biolabs’ Phusion High-Fidelity PCR Master Mix with GC Buffer was used. High-efficiency and high-fidelity enzymes were employed for PCR: bacterial 16S rDNAprimers: 515F(5′-GTGCCAGCMGCCGCGGTAA-3′)/806R(5′-GGACTACHVGGGTWTCTAAT-3′); fusion primer design: forward primer for fusion V4 region: adapter+bracode+GTGCCAGCMGCCGCGGTAA F; reverse primer for fusion V4 region: adapter+bracode+GGACTACHVGGGTWTCTAAT R. The PCR reaction system (50 μL) consisted of DNA (30 ng) PCR Primer Cocktail*(4 μL) PCR Master Mix (25 μL)H2O(as needed). The PCR amplification reaction conditions were as follows: 94 °C pre-denaturation for 3 min; 98 °C for 45 s, 55 °C for 45 s, and 72 °C for 45 s for 30 cycles; 72 °C extension for 7 min. Library fragment recovery, purification and detection were performed using magnetic beads to screen the target Amplicon fragments; an Agilent 2100 Bioanalyzer (reagent: Agilent DNA 1000 Kit, cat No.5067–1504, manufacturer: Agilent) was used to detect the range of insert fragments in the library. An ABI StepOnePlus Real-Time PCR System (TaqMan Probe) was used to quantify the concentration of the library. According to the results of the library test, the samples are mixed and loaded on the machine, and the sample amount of each library was calculated. All libraries were mixed at a ratio of 1:1, and after mixing, the samples are fully shaken and sequenced using the Miseq sequencing platform. The sequencing was commissioned by Huada Gene Corporation.

Sequencing data processing

Data quality control was performed with QIIME software [5] to filter low-quality data, and the linker and primer sequence, polybase N, poly A/T tail and low-quality bases at the end of the sequence were removed offline to obtain high-quality data. Data splitting was performed using the barcode sequence to split the data into different sample data. The allowed number of mismatches between the barcode sequence and sequencing reads was 0 bp; the barcode identification sequence and PCR amplification primer sequence were cut off, and the number of sequences in each sample was statistically analysed. Tag splicing was performed with FLASH software [6] to splice a pair of overlapping sequences in each sample. The tags were truncated from the first low-quality base site where the number of consecutive low-quality values (default setting <=10) reached the set length (default of 5). Tags less than 70% of the average length of tags in the collection were filtered out after interception.

Species classification and abundance and analysis of pathways related to CIN

Species classification analysis

The software USEARCH (v7.0.1090) [7] was used to cluster the spliced tags into OTUs, and the sequences were clustered into operational taxonomic units (OTUs) based on 97% similarity. With a similarity of 97%, the number of OTUs for each sample was obtained, and the Venn diagram was used to show the number of common and unique OTUs for multiple samples and to visually show the overlap of OTUs between samples. The species represented by OTUs were combined to identify the core microorganisms in different environments. According to the abundance file of each sample OTU in each sample, the OTU of each sample or group was calculated to generate a Venn diagram using the VennDiagram package in R (v3.0.3) language. The obtained OTU representative sequence was annotated to obtain the species classification.

Species abundance and pathway analysis related to CIN

PICRUSt software uses marker gene data and reference genome databases to predict the functional composition of metagenomics. This method was first published in Nature Biotechnology in 2013 by Langille et al. [8]. We used the PICRUSt software in the KEGG database: Kyoto Encyclopaedia of Genes and Genomes database (59.3) (http://www.genome.jp/kegg/). After normalization according to the different 16 s rDNA copy numbers contained in each species, the possible abundance of each species in each sample was obtained, and the abundance value of each function was multiplied by the abundance value of the species and the abundance value of its potential function. Next, we combined the respective abundance values of all functions in the same sample to convert them into relative abundance values. A non-parametric rank sum test on the KEGG classification data was performed, using P < 0.05 and a false discovery rate (FDR) < 0.3 to denote a difference. COREMINE was used for biological annotation and text mining of differential pathways to further illustrate the relationship between differential signalling pathways and CIN.

Bioinformatics analysis of differential pathways, preliminary screening of candidate genes and biological annotation of their correlation with the progress of bacteria and CIN

The KEGG (https://www.kegg.jp/kegg/) database was again applied to query the genes contained in the five differential pathways. The search tool STRING online database was employed to construct a protein-protein interaction network diagram (PPI) for the genes of each pathway, Cytoscape was used to visualize the analytic results, and then MOCDE was applied to cluster the results according to the given network diagram to identify areas of dense connections. The selection criteria were as follows: MOCDE scores> 5 points, degree cut-off = 2, node score cut-off = 0.2, Max depth = 100 and k-score = 2. Combined with the PubMed literature search, 10 candidate differential genes were initially screened out. The correlation between candidate differential genes and the progress of bacteria and CIN was biologically annotated using the online software COREMINE (http://www.coremine.com/medical/).

Expression verification and analysis of candidate genes

Tissue total RNA was extracted using the TRIzol method, and the total RNA concentration and OD value (260/280) were determined using the SYnergy multifunctional microplate fluorescence analyser. The total RNA reverse transcription kit and real-time fluorescent quantitative PCR kit SYBR Premix Ex Taq™ II (TliRNaseH Plus) were provided by Takara, Japan. The synthesis of primers was completed by Takara Corporation of Japan, using β-Actin (ACTB) as the internal reference gene. The primer sequences are shown in Table 1. The fluorescence quantitative PCR reaction conditions were as follows: 95 °C predenaturation for 5 min, 95 °C for 30 s, 1 cycle; 95 °C for5 s, 60 °C for 30 s, 40 cycles; and 60 °C for 30 s, 40 cycles. ROC curve analysis was used to predict the value of candidate genes for early warning signs of CIN occurrence and progression, and logistic regression was used to analyse the risk factors related to CIN occurrence and progression.

Table 1 Fluorescence quantitative PCR primer sequence

Results

Analysis of clinical data of samples

Among the 62 samples, the oldest was 62 years old and the youngest was 22 years old. The average age was 33 ± 5.12 years. There was no significant difference in age between groups (P > 0.05), and there were no significant differences in parity or contraceptive methods (Table 2) (P > 0.05). The infection rate of high-risk HPV (hrHPV) also gradually increased with the aggravation of cervical lesions (P < 0.05).

Table 2 Analysis of clinical data of samples

Sequencing data and OTU construction

A total of 60 samples in 5 groups of specimens were analysed by high-throughput sequencing (2 samples were unqualified): Group 1 (normal group) 14 cases; Group 2 (CIN1 group) 9 cases; Group 3 (CIN2 group) 11 cases; Group 4 (CIN3 group) 17 cases; Group 5 (cervical cancer group) 9 cases. The 16S rDNA gene fragments were extracted from the samples, with the exception of low-quality sequences. The tags were spliced into tags through the overlap relationship between reads. A total of 1,873,442 tags were obtained for all samples, with an average of 31,224 tags per sample; the average tag length was 253 bp, and the SD value was 7 bp. The statistical results for the tags from each sample were as follows: A1-A14: normal group; B1-B9: CIN1 group; C1-C11: CIN2 group; D1-D17: CIN3 group; E1-E9: cervical cancer group. A total of 1896 OTUs were generated from 60 samples. For each sample OTU, the shared and unique OTUs between samples or groups and the Venn diagram between samples or groups are shown in Fig. 1.

Fig. 1
figure 1

OTU venn diagram (Group 1, Normal group; Group 2, CIN1 group; Group 3, CIN2 group; Group 4, CIN3 group; Group 5, Cervical cancer group)

Species classification analysis results

From the class level analysis, β-Proteobacteria, γ-Proteobacteria, Bacillus, α-Proteobacteria, and Clostridium were the dominant strains in the normal group, CIN group and cervical cancer group, and each group had different dominant strains, namely, the normal group consisted of δ-Proteobacteria, the CIN1 and CIN2 groups were Sphingomyces, and the CIN3 and cervical cancer groups were Bacteroides (Table 3).

Table 3 Class classification level strains

At the subordinate level, the dominant strains of Lactobacillus accounted for a large proportion in each group, but some samples in the normal group still did not show the dominant Lactobacillus. The other main microorganisms in CIN1/2 were Stenotrophom, Chitinophaga, and Acinetobacter in CIN3 and in cervical cancer were Halomonas, Shewanella, and Acinetobacter (Table 4).

Table 4 Genus classification level strains

Species abundance and pathway analysis using the CIN group flora

The OTU is classified into 23 classes, 41 orders, 56 families, and 73 genera in several classification levels of class, order, family, and genus. The species abundance of each sample is shown in Fig. 2 and Fig. 3. Using PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) software for species abundance analysis, combined with KEGG data, the results show that: the genes contained in the CIN group are mainly involved in the following pathways: Transporters, ABC transporters, DNA repair and recombination proteins, General function prediction only, Purine metabolism, Two-component system (Fig. 4).

Fig. 2
figure 2

Analysis of pathways in which the flora contains genes mainly involved

Fig. 3
figure 3

Linear relationship diagram between differential pathway and CIN

Fig. 4
figure 4

Interaction analysis between gens obtained by differential signaling pathways and CIN progression-related proteein genes

COREMINE biological annotation of differential pathways

The biological annotations of the screened differential pathways were carried out using COREMINE to further explain the relationship between the differential pathways and CIN progression. As shown in the figure, using the differential pathway name and CIN (cervical intraepithelial neoplasia) as keywords to conduct a co-occurrence analysis of the literature, a total of 5 differential pathways were found in the annotation in the network: ABC transporters, Base excision repair, Energy metabolism, Lipid biosynthesis proteins, Protein kinases. It is suggested that these differential pathways are related to CIN (Fig. 3).

Interaction analysis between genes obtained by differential signalling pathways and CIN progression-related protein genes

A total of 79 human genes in the above 5 differential signalling pathways were queried, such as APEX2, POLB, and POLL, among others. Using the String tool (http://string-db.org/), a comprehensive analysis of the protein interaction between 79 genes and 24 known genes related to CIN progression was performed to explain the overall relationship between this gene and CIN progression. The genes that interacted with 24 known genes to form a protein interaction network were ABCC1, ABCC3, ABCC5, ABCC6, ABCG2, APEX2, ARAF, ATM, FEN1, GBA, HMGB1, LIG1, LIG3, MCM3AP, MDK, OGG1, PCNA, POLB, SMUG1, TDG, UNG, and XRCC1 (Fig. 4), indicating that these genes might be functionally related and suggesting that the selected genes might also be related to CIN. In particular, ATM and 11 genes related to CIN progression (TP53, TP63, EGFR, ERBB2, PIK3CA, CDKN2A, CTNNB1, FHIT, MYC, STK11, TP53BP1) interacted directly; PCNA and 10 genes related to CIN progression Genes (CDKN2A, TP53, EGFR, MYC, CCL2, ERBB2, MKI67, CTNNB1, FHIT, TP53BP1) interacted directly; ABCG2 and 5 genes related to CIN progression (CTNNB1, EGFR, ERBB2, MYC, TP53) interacted directly; MDK and 5 genes related to CIN progression (MAPK1, TP53, ERBB2, PIK3CA, EGFR, MYC, PTGS2) interacted directly.

Protein-protein interaction network of 5 differential pathway genes

The STRING online database was used to construct a PPI network diagram for the genes of each pathway, and Cytoscape and MOCDE were applied to visually analyse and cluster the PPI and then obtain 5 protein interaction maps. Among them, the larger the node, the higher was the credibility that the gene was related to other genes in the network. The details are shown in Figs. 5, 6, 7, 8, and 9.

Fig. 5
figure 5

Interaction diagram of ABC transporters pathway genes

Fig. 6
figure 6

Interaction diagram of base excision repair pathway genes

Fig. 7
figure 7

Interaction relationship diagram of energy metabolism pathway

Fig. 8
figure 8

Interaction diagram of lipid biosynthesis protein pathway genes

Fig. 9
figure 9

Interaction diagram of protein kinases pathway genes

Literature search and genetic screening results

Through PUBMED, a literature search was conducted with keywords such as gene name, cervical cancer, cervical intraepithelial neoplasia, CIN, cancer, etc., to gain a preliminary understanding of the relationship between genes and tumours. Through protein interaction mapping and a literature search, we initially selected ATM, ABCG2, PCNA, XRCC1, HMGB1, OGG1, LIG1, SMUG1, FEN1, and TDG as candidate differential genes. Among them, ATM and ABCG2 were present in the ABC transporters pathway, and PCNA, XRCC1, HMGB1, OGG1, LIG1, SMUG1, FEN1, and TDG were present in the Base excision repair pathway.

The relationship between candidate genes and the progress of bacteria and CIN

Biological annotation and text mining were carried out through the COREMINE database to further illustrate the relationship between the 10 candidate differential genes and the progress of bacteria and CIN. As shown in Fig. 10, the co-occurrence analysis of the literature with the key words of gene name, CIN (Cervical intraepithelial neoplasia), bacterial name (Stenotrophom, Chitinophaga, Acinetobacter, Halomonas, Shewanella) and disease progression (Disease Progression) resulted in 10 genes all in the annotation network. Among them, PCNA, HMGB1, OGG1, SMUG1, and XRCC1 had a direct network relationship with the progress of CIN, while ATM, ABCG2, TDG, LIG1, and FEN1 had an indirect relationship. Excluding Chitinophaga, information on which could not be found, the remaining four bacterial genera were researched. The figure shows a complex network relationship between candidate genes and the progress of bacteria and CIN.

Fig. 10
figure 10

The linear relationship between candidate differential genes and the progress of bacteria and CIN

Expression verification of candidate genes

Patient clinical characteristics

Among the 38 normal, 52 CIN and 30 cervical cancer patients, average age was 43.55 ± 9.54, 42.12 ± 9.67, and 49.10 ± 10.49 years. The average age of the cervical cancer group was higher than the normal and CIN groups (P<0.05), and the number menopausal patients was greater in the cervical cancer group than the normal and CIN groups (P<0.05). The HPV infection rate in the normal group, CIN group, and cervical cancer group was 39.47, 90.38, and 93.33%, respectively. The HPV infection rate gradually increased with the degree of cervical lesions (P < 0.05). TCT also increased with the severity of cervical lesions (P < 0.05). However, there were no significant differences in the pregnancy and parity numbers among the included cases (P > 0.05). Age, menopause, HPV infection, and the severity of TCT might affect the occurrence and development of CIN (Table 5).

Table 5 Clinical characteristics of patients in each group

Comparison of target gene mRNA expression in different cervical tissues

The fluorescence quantitative PCR results showed that ATM, ABCG2, PCNA, OGG1, TDG, LIG1, and HMGB1 expression levels were statistically significant in three different cervical tissues (P < 0.05); There were no differences in the expression levels of XRCC1, SMUG1, and FEN1 (P > 0.05). Among them, The expression of ATM was higher in the cervical cancer group than the CIN group (P < 0.05). There was no difference between the normal and CIN group and the cervical cancer group (P > 0.05). The expression levels of PCNA, TDG, and LIG1 were elevated in cervical cancer group compared with the normal and CIN groups (P < 0.05). There were no differences between the normal group and CIN group (P > 0.05). The expression of ABCG2 was lower in the cervical cancer group than the normal and CIN groups (P < 0.05). There were no differences between the normal and CIN groups (P > 0.05). The expression level of OGG1 was lower in the CIN group than the normal group (P < 0.05) but higher than in cervical cancer group (P < 0.05). Compared with the cervical cancer group, the normal group showed no differences (P > 0.05). The expression of HMGB1 gradually increased with the aggravation of cervical lesions (P < 0.05). See Table 6 and Fig. 11 for details.

Table 6 Real-time fluorescence quantitative PCR detection of target gene mRNA expression in each tissue (copies/μl RNA x ± S)
Fig. 11
figure 11

ATM, ABCG2, PCNA, OGG1, TDG, LIG1, and HMGB1 relative expression in different cervical tissues

ROC curve analysis of differential genes with respect to the occurrence and progression of early warning CIN

ROC curve analysis showed the highest value for genetic OGG1 for early warning signs of CIN, with an AUC of 0.68 (P < 0.05). The HMGB1 gene had an AUC of 0.63 (P < 0.05). The sensitivity of OGG1 and HMGB1 as warning signs of CIN was 57.7 and 75.0%, and the specificity was 78.9 and 57.9%, respectively. The values for the residual genes TDG, ATM, ABCG2, LIG1, and PCNA revealed no association with CIN. Early warning of the progress of CIN was observed for HMGB1, LIG1, ABCG2, TDG, and PCNA with an AUC of 0.74, 0.68, 0.67, 0.65, and 0.65 (P < 0.05). The sensitivity was 86.7, 83.3, 96.7, 50.0, and 76.7%, and the specificity was 57.9, 52.6,39.5, 81.6 and 57.9%, respectively. The values for the residual genes ATMand OGG1 showed no association with early warning of CIN progress. Genes with value for the occurrence and progression of CIN early warning signs were OGG1 and HMGB1 and HMGB1, LIG1, ABCG2, TDG, and PCNA, respectively. HMGB1genes were valuable for early warning of the occurrence and progression of CIN (see Tables 7 and 8 for details).

Table 7 The AUC value, sensitivity, and specificity of differentially expressed genes for early warning of CIN occurrence
Table 8 The AUC value, sensitivity and specificity of differentially expressed genes for early warning of CIN progression

Logistic regression analysis of risk factors for CIN occurrence and progression

Logistic regression analysis that promotes risk factors for CIN

Age, menopause, HPV infection, TCT and ROC analysis of meaningful genes were subjected to single-factor logistic regression analysis, and meaningful factors were then analysed using multi-factor logistic regression analysis. The results revealed that during the progression of the normal cervix to CIN, HPV positivity and TCT severity were risk factors for CIN (see Table 9 for details).

Table 9 Logistic regression analysis of risk factors that promote CIN

Logistic regression analysis of risk factors for CIN progression

Age, menopause, HPV infection, TCT and ROC analysis of meaningful genes were subjected to single-factor logistic regression analysis, and meaningful factors were then subjected to multi-factor logistic regression analysis. The results showed that the low expression of ABCG2 and high expression of PCNA and TDG were risk factors for the progression of CIN (see Table 10 for details).

Table 10 Logistic regression analysis of risk factors for CIN progression

Establishment andselection of Random Forest Models

According to the results of the multivariate logistic regression analysis, the risk factors promoting the occurrence and development of CIN, such as HPV positivity, TCT severity, ABCG2, PCNA and TDG genes, were used to establish random forest models with different combinations of these indicators. A total of 7 models were established, as shown in Table 11. According to the accuracy, AUC value and OOB error value of each model, the best model was selected. Among the various models, the model 4 accuracy rate and AUC value were relatively large, and the OOB error value was the smallest. Therefore, model 4 (ABCG2 + PCNA+TDG) was selected as the best early warning model for the occurrence and progress of CIN.

Table 11 The occurrence and progression of CIN evaluated by the SCC random forest model

Discussion

The relationship between microbiology and CIN progression

Persistent high-risk human papillomavirus infection is undoubtedly the main carcinogen leading to CIN and cervical cancer. However, studies have found that not all patients with cervical HPV infection experience development into cancer [1]. Although it is now believed that auxiliary factors other than HPV play a key role in the development of tumours, most of the potential mechanisms of this carcinogenic effect are still unknown [2]. These carcinogenic factors include genetic factors, environmental factors, immune function, cellular defense, and body-specific genes and cytokines [9]. Evidence shows that the cause of human diseases is not only a single pathogen, but also involves the overall changes in the human microbiology group [3]. However, few studies have examined the progress of the microbiome in CIN at this stage. The available small numbers of cervical and vaginal microbiological studies based on high-throughput sequencing have mainly focused on other populations. The structural characteristics of the cervical flora in the normal and disease states in the Chinese population are not clear. Therefore, in this study, normal cervical tissues, CIN and cervical cancer patients in southern China served as the research objectives. Using high-throughput sequencing analysis technology based on bacterial 16S rDNA, the composition of the cervical flora was analysed and compared in the above population using bioinformatics methods. To construct an atlas of the composition of the significant specific flora of the disease, we selected specific flora or species that were significantly related to the disease. This project had the following main research goals: 1. preliminary understanding of the structure and diversity of the bacterial community in normal cervical tissue; 2. preliminary understanding of the structure of the bacterial community in patients with CIN and cervical cancer without treatment by pathology and diversity characteristics; 3. Attempt to find and identify functional pathways of genes present in related flora or species related to CIN progression.

The difference and clinical significance of the main bacterial groups among CIN, cervical cancer and normal cervix

In this study, Illumina’sMesiq high-throughput sequencing platform was used to conduct high-throughput analysis of cervical tissue-related bacterial communities, and the data obtained were subjected to quality control such as removal of single reads; qualified sequences were annotated by OTUs and evaluated by inter-sample diversity analysis such as dilution curves, among others. We confirmed that the amount of sequencing in this study was sufficient to cover cervical bacterial species, and the sample richness and uniformity between groups were good. In addition, there were no significant differences in age, parity, or contraceptive method, among others, in the selected cases, and some interfering factors were excluded. In general, cervical tissues were identified from 23 bacterial strains and 73 genera. Regarding the level distribution of vaginal-related bacterial communities, studies [10] have indicated that the phyla Firmicutes, Actinobacteria, Bacteroidetes, and Proteobacteria, and non-wall bacteria Tenericutes, Fusobacteria and TM7, are the main microorganisms. It has also been reported that HSIL vaginal microorganisms are characterized by higher levels of Sneathiasanguinegens, Anaerococcustetradius and Peptostreptococcusanaerobius, and lower levels of LSIL in Lactobacillus jannaschii. We found that the main microorganisms in CIN1/2 were Stenotrophom, Chitinophaga, and Acinetobacter, and the main microorganisms in CIN3 and cervical cancer were Halomonas, Shewanella, Acinetobacter. By comparison with the results of several existing high-throughput sequencing analyses of vaginal flora [4, 10, 11], the relative abundance of cervical bacterial communities at the phylum and genus levels and other species revealed certain differences, considering the possible relationships with factors such as disease state, race, and region.

Studies have shown that Lactobacillus is the most dominant bacteria in the vagina [12]. Symbiotic vaginal Lactobacillus is believed to produce species-specific metabolites and bacteriocins by maintaining a constant pH value, thereby destroying pathogen biofilms [13]. However, most current studies on the microbiota focus on the vaginal microbiota, and there are few reports on the microbiome in cervical tissue. At the subordinate level of this study, Lactobacillus was still dominant in most of the normal group and CIN group, but the proportion was significantly reduced in the CIN and cervical cancer groups. This result further supports the ability of Lactobacillus to maintain a normal vagina and cervix and provide protection from the environment. However, a small number of samples in the normal group still lacked a dominant Lactobacillus genus, potentially due to the differences between cervical tissue and vaginal microbes. The small sample size was not ruled out, which will be expanded in future detailed research.

The role of the microbiome in the progression of CIN is gradually becoming recognized. MitraA [4] recently performed 16S rRNA gene amplification of the vaginal wall microbes of 52 cases of LSIL, 92 cases of HSIL, 5 cases of ICC and 20 cases of normal controls. The results showed that the vaginal microbial diversity was associated with the severity of CIN disease. We can participate in regulating the persistence of viral infections and disease progression. However, no reports have examined the mechanisms and gene functions involved in the progression of CIN related to the cervical microbiome. The 16S rRNA is an important tool for microbial community research, but it does not provide direct evidence of community function. In 2013, Langille et al. [8] described a software that uses marker gene data and reference genome databases to predict the functional composition of metagenomics, namely PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States). PICRUSt uses the extended ancestral state reconstruction algorithm to predict the presence of specific gene families, and then combines the gene families to estimate the compound metagenomics. Using 16S information, PICRUSt can retrieve key results from the human microbiome and accurately predict the abundance of gene families associated with the host and environmental communities. In 2016, Tseng and other scholars [14] used 16S rRNA sequencing to identify the microbiome changes in gastric cancer tissues before and after surgery, and PICRUSt software to predict the classification of gene functional groups related to the microbiome. In this study, PICRUSt software was used to analyse the selected pathways, and the obtained differential pathway information was biologically annotated by COREMINE. Five differential pathways were found: ABC transporters, Base excision repair, Energy metabolism, Lipid biosynthesis proteins, Protein kinases Exist, suggesting that these 5 differential pathways were related to CIN. Using the String tool, among these differential signalling pathways, the human genes ATM and PCNA showed the closest relationship with the progression of CIN. We will further examine the relationship between these pathways and related genes and the progression of CIN to inhibit the progression of this disease.

Although HPV infection is common in sexually active women, most infections are temporary. Only some women infected with the HPV wasprogressed to precancerous lesions and invasive cervical cancer. A large number of research reports suggest that persistent high-risk human papillomavirus (hrHPV) infection is the greatest risk factor for the occurrence and development of CIN and cervical cancer. The mechanism by which HPV infection continues to cause cancer is not well understood. Studies have suggested that the hrHPV subtype itself produces two proto-oncoproteins, E6 and E7, which are key to its pathogenesis [15]. E6 and E7 proto-oncoproteins can change the terminal differentiation of host keratinocytes and combine with host cell cycle regulatory proteins to destroy the negative regulation of the cell cycle, leading to abnormal cell cycle regulation [16]. Additionally, normal tissues eventually develop into CIN and cervical cancer. Recently, Hu Z et al. [17] used whole-genome sequencing and capture sequencing to determine for the first time the impact of HPV-infected human specific integration sites and key fragile integration sites on tumour cells. Whole-genome analysis of HPV integration in CIN and cervical cancer is expected clarify the carcinogenic mechanism and block the occurrence and development of CIN.

The relationship between candidate genes screened by differential pathways and CIN and cervical cancer

In this study, five differential pathways related to the progression of CIN in the cervical microbial community were analysed. The KEGG database was used to search for the pathway genes, and approximately 1442 genes were queried. The STRING online database was used to construct a PPI network for the genes in each pathway. Cytoscape was used to visualize the network diagram, and then MOCDE was applied to further screen the genes. Combined with the literature search, 10 candidate differential genesrelated to CIN and cervical cancerwere initially screened: ATM, ABCG2, PCNA, XRCC1 HMGB1, OGG1, LIG1, SMUG1, FEN1, and TDG.

A literature search revealed that ATM is highly expressed in cervical cancer, and ATM expression inhibition can enhance the sensitivity of cervical cancer to radiotherapy and chemotherapy [18]. In addition to being related to cervical cancer, ATM is also related to other malignant tumours. For example, Santarpia et al. [19] found that low ATM expression in breast cancer was related to a high DNA mutation rate, tumour progression and increased angiogenesis, and ATM expression was related to patient prognosis. Kim et al. [20] usedimmunohistochemical analysis of 321 gastric cancer patients and found low ATM levels in 205 gastric cancer patients and high expression in 116 gastric cancer patients. The 5-year survival rate and overall survival rate of patients in the ATM low expression group were compared. The patients in the ATM high expression group showed worse outcomes (P = 0.017, P = 0.027).

ABCG2 is a member of the G subfamily of the ATP-binding cassette (ABC) transporter superfamily. It is a multidrug resistance gene, also known as breast cancer resistance protein (BCRP), and it is found in a variety of malignant tumours. High expression was found in, e.g., liver cancer and pancreatic cancer [21, 22], but low expression in colorectal and cervical cancer [23].

PCNA is the core protein in DNA synthesis, replication, and damage repair and a key factor in cell cycle regulation. It can be used as an indicator of tumour cell proliferation and tumour prognosis. Studies have shown that the expression of PCNA in cervical tissues increases with the increase in CIN and cervical cancer grade [24]. A systematic review showed that the expression of PCNA is closely related to the staging and prognosis of cervical cancer. Cervical cancer with high PCNA expression has a lower survival rate and a worse prognosis [25]. In addition to its high expression in cervical cancer, PCNA is highly expressed in rapidly dividing tumour cells, which in most cases is related to poor prognosis and is an effective target for tumour treatment [26].

HMGB1 is a member of the HMG superfamily. By participating in the functions of cell DNA recombination, replication, repair, and gene transcription, it is widely involved in tumorigenesis, growth, invasion and metastasis. The literature also indicates that HMGB1 is highly expressed in cervical cancer tissues [27]. In ovarian and gastric cancer, HMGB1 knockdown can inhibit cancer cell proliferation, migration and invasion [28, 29]. Wu et al. [30] conducted a meta-analysis of the relationship between HMGB1 and various tumours and found that overexpression of HMGB1 was associated with a poor prognosis of various tumours. HMGB1 is a potential marker for the prognosis of various tumours.

OGG1 is an enzyme that specifically recognizes 8-dihydro-8-oxoguanine (8-oxoG) in the body and excises and repairs it. The presence of this enzyme reduces the risk of DNA mutations and tumorigenesis [31]. Kubo et al. [32] detected low OGG1 expression in oesophageal cancer and its association with lymphatic metastasis and tumour staging. This study showed that the DNA repair pathway involved in OGG1 was closely related to oesophageal cancer, but its expression in CIN and cervical cancer has not yet been reported.

At present, most studies have shown that the XRCC1 gene polymorphism is related to the risk of cervical cancer, and few studies have examined its expression in cervical cancer [33]. Abdel-Fatah et al. [34] found that XRCC1 is highly expressed in ovarian cancer and is closely related to tumour staging, platinum resistance, and death outcomes. Patients with positive XRCC1 expression are more likely to have adverse clinicopathological and survival outcomes, which is considered a new predictive marker for ovarian cancer.

FEN1 is a multifunctional nuclease with a special structure. Its gene mutation may cause autoimmune diseases, chronic inflammation, and tumour susceptibility. He et al. [35] found that FEN1 was up-regulated in cervical cancer tissues. The use of FEN1 inhibitors in combination with paclitaxel can significantly improve the efficacy of paclitaxel in cervical cancer. In addition to its high expression in cervical cancer, FEN1 is overexpressed in many tumours such as breast and lung cancer [36, 37].

Few studies have examined LIG1, SMUG1, and TDG in CIN and cervical cancer. At present, their specific expression in cervical cancer has not been described, but they are closely related to other malignant tumours. It is known that abnormal DNA methylation is an important feature of tumorigenesis. Increased demethylation may be involved in tumorigenesis, and TDG plays an important role in DNA demethylation [38]. Yang et al. [39] found low TDG in early breast cancer and its relationship with a poor prognosis of early breast cancer. Zhong et al. [40] used the Whole Transcriptome Association Study (GWAS) to identify a relationship of SMUG1 with the risk of pancreatic cancer, making it a new candidate gene for predicting pancreatic cancer. At present, the literature shows that LIG1 is a DNA ligase that is usually overexpressed in tumours. LIG1 is considered to be necessary for ligating Okazaki fragments during DNA replication. Therefore, it is also necessary for cell survival and the focus of current clinical efforts. Research on DNA ligase inhibitors to inhibit cancer progression have been reported [41].

We used the COREMINE database to further conduct biological annotation and text mining regarding the relationship between these 10 candidate genes and the progress of bacteria and CIN, and we found that the 10 genes were all present in the annotation network. Among them, PCNA, HMGB1, OGG1, SMUG1, and XRCC1 showed a direct network relationship with the progress of CIN, and ATM, ABCG2, TDG, LIG1, FEN1 and CIN progression showed an indirect network relationship. Chitinophaga could not be found, but information on the other four bacterial genera was available. Among them, ABCG2 had a direct network relationship with the other four bacterial genera, ATM, PCNA and Acinetobacter. There was a direct network relationship, and Acinetobacter (Acinetobacter) was related to HPV infection and CIN progression [42]. OGG1 showed a direct network relationship with Stenotrophom. As can be seen from the network diagram in which these genes showed a complex network relationship with the progress of bacteria and CIN. However, this result has not yet clarified the specific mechanism of genes and bacteria and the progress of CIN, nor the specific role of bacteria in the progression of CIN. Maybe these genes are the cancer promoting genes groupin the cervical microbial community. In follow-up studies, we will further verify the expression of the selected candidate differential genes in cervical tissues. We hope to find specific genes and conduct detailed studies to block CIN progression.

Expression verification and analysis of candidate genes in CIN

As mentioned above, persistent hrHPV infection is the main reason for the occurrence and development of CIN and cervical cancer. Some studies have suggested that the hrHPV subtype itself produces E6 and E7, two key proto-oncoproteins in its pathogenesis. If HPV oncogenes such as E6 and E7 are integrated into the high-risk area of the patient’s chromosome, this process will interfere with tumour suppressor genes such as P53 or PRb, among others, which may cause CIN and even cervical cancer [43]. In our study, analysis of the patients’ clinical data showed that age, menopause, TCT test results, and HPV infection differed among the three groups of patients, and the HPV infection rate gradually increased with the severity of cervical lesions. Additionally, the cervical cancer group was obviously higher than the normal group. The data once again confirmed that the occurrence of cervical lesions was closely related to HPV, consistent with domestic and foreign reports. TCT examination played an important role in early warning of the occurrence of CIN, and TCT examination combined with HPV detection showed important value in the early diagnosis of CIN. Studies have shown that the residual rate, recurrence rate and the incidence of aggressive lesions after CIN conization are related to age and menopause [44], but whether they can promote the occurrence and progression of CIN remains to be examined.

This study verified the expression of 10 candidate differential genes screened from cervical microbes related to CIN progression pathways in different cervical tissues through fluorescence quantitative PCR: ATM, ABCG2, PCNA, OGG1, TDG, LIG1, and HMGB1. The expression in three different cervical tissues was statistically significant (P < 0.05). Other reports have also indicated that ATM, PCNA, and HMGB1 are highly expressed in cervical cancer tissues [24, 27, 45], and ABCG2 is expressed at low levels in cervical cancer tissues [23]. The results of this experiment are consistent with those reports. However, expression levels of TDG, LIG1 and OGG1 have not been examined in cervical cancer tissues. In this experiment, we found that TDG and LIG1 were expressed at high levels in cervical cancer, and OGG1 expression was low in cervical cancer (P < 0.05).

Through ROC curve analysis, we found that the genes that were valuable for early warning of CIN occurrence and progress were OGG1 and HMGB1 and HMGB1, LIG1, ABCG2, TDG, and PCNA, and the gene with value for early warning of CIN occurrence and progress was HMGB1. This result is further confirmed that these genes are the key carcinogenic genes in the cervical microbial community. Additionally, the expression of HMGB1 increased with the severity of cervical lesions, suggesting that HMGB1 was closely related to the occurrence and progression of CIN, and HMGB1 might be a new marker for early warning of the occurrence and progression of CIN. We further searched the literature to understand the expression of HMGB1 in various tumours, and we found that HMGB1 was expressed at high levels in other tumours such as ovarian cancer and gastric cancer; overexpression of HMGB1 was related to the poor prognosis of various tumours. HMGB1 is associated with and may be a potential marker in the prognosis of various tumours [28,29,30]. Studies have shown that HMGB1 is a chromatin component that is ubiquitous in mammalian cells. It has a variety of biological functions and plays an important role in cell migration, inflammation, cell differentiation and tumour metastasis. It has a role as a signal regulator of various biological functions such as the mediation of infection, injury and inflammation, promotion of autophagy, induction of cell death, and activation of natural immunity, playing a key role in the occurrence and development of tumors [46].

Logistic regression analysis showed that HPV infection and the severity of TCT were independent risk factors that could promote the occurrence of CIN, and low expression of ABCG2 and high expression of TDG and PCNA were risk factors that could promote the progression of CIN. ABCG2 is a member of the G subfamily of the ATP-binding cassette transporter superfamily. It is a multidrug resistance gene, also known as breast cancer resistance protein (BCRP), which has anti-tumour and anti-toxic effects [47]. Gupta et al. [23] used fluorescent quantitative PCR and immunohistochemistry to compare the expression of ABCG2 in normal and cervical cancer tissues, and they found low expression of ABCG2 in cervical cancer tissues, consistent with our results. Its down-regulation in tumours may play a role in tumour development by increasing tissue exposure to oncogene toxins and excessive production of nitric oxide. PCNA is the core protein in DNA synthesis, replication, and damage repair and a key factor in cell cycle regulation. It can be used as an indicator of tumour cell proliferation and tumour prognosis. Studies have shown that the expression of PCNA in cervical tissues increases with the increase in CIN and cervical cancer grade [24].TDG plays an important role in DNA demethylation [38]. But few studies have examined TDG in CIN and cervical cancer, so its role in the progression of CIN requires further analysis. It may be related to DNA methylation during the malignant transformation of cervical cells.

Construction and screening of the early prediction model for evolution of CIN malignant transformation

According to the results of multivariate logistic regression analysis, risk factors that promote the occurrence and development of CIN, such as HPV positivity, TCT severity, ABCG2, PCNA and TDG genes, were used to establish random forest models with different combinations of these indicators. The results showed that model 4 (ABCG2 + PCNA+TDG) had the greatest accuracy and AUC value and the smallest OOB error value. Therefore, thecombination of ABCG2 + PCNA+TDG genes was selected as the best early prediction model for evolution of CIN malignant transformation [48].

Summary

To our knowledge, this is the first study to use 16S rDNA sequencing to detect vaginal microbes at different levels of CIN and normal and cervical cancer tissues. The analysis suggests that with the severity of disease, the microbes show diversity, combined with high-risk HPV infection and disease levels. The relationship between vaginal microbial diversity may be related to the severity of high-risk HPV infection and CIN disease. This study is expected to clarify the microbiome changes related to the progression of CIN. Subsequently, we preliminarily screened out 10 candidate differential genes related to the progression of CIN in the cervical microbial community through bioinformatics analysis. Biological annotation and text mining of the relationships between genes, bacteria and CIN progress through the COREMINE database revealed a complex network relationship between candidate differential genes and the progress of bacteria and CIN. These 10 genes may be important candidate genes for community microbes to regulate CIN progression. Research based on these genes may provide new targets for blocking the progression of CIN. In the verification of these genes, we used the fluorescent quantitative PCR method to verify the expression of the candidate genes and found that ATM, ABCG2, PCNA, OGG1, TDG, LIG1, and HMGB1 were significantly expressed in three different cervical tissues. They might cooperate with HPV to participate in the occurrence and progression of CIN. Through ROC and logistic regression analyses, the factors associated with the occurrence and progression of CIN were initially screened out as HPV infection, TCT severity, HMGB1, ABCG2, and TDG. We used these indicators to establish a random forest model. Seven models were built through different combinations. The model 4 (ABCG2 + PCNA+TDG) had the highest accuracy and AUC value and the smallest OOB error value. Therefore, the combination of ABCG2 + PCNA+TDG genes was selected for the best early prediction model for evolution of CIN malignant transformation.

Current studies have found that the vaginal microecology is closely related to the occurrence and progression of CIN, but there has not been a study on the gene pathways contained in the flora or strains. This study firstly identified the genes from cervical microbial community that play an important role in the occurrence and progression of CIN. At the same time, the early warning model including ABCG2 + PCNA+TDG genes provided a new idea and target for clinical predictionand blocking the evolution of CIN malignant transformation from the aspect of cervical microbiological related genes. We comprehensively analyzed the risk factors for the occurrence and progression of CIN from clinical factors, microbial factors, and genetic levels. For the first time, we combined clinical factors and differential genes, and used the random forest model to select the best early warning model for the occurrence and progression of CIN. The purpose of this study is to find out the functional pathways of genes contained in the relevant flora and species related to the progress of CIN, to find the best early warning model for the occurrence and progression of CIN, and may provide an accurate clinical screening of CIN.

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the NCBI repository, http://www.ncbi.nlm.nih.gov/bioproject/815961. Other datasets are available from the corresponding author on reasonable request.

Abbreviations

CIN:

Cervical intraepithelial neoplasia

PCR:

Polymerase chain reaction

ROC:

Receiver operating characteristic curve

HPV:

Human papilloma virus

TCT:

Thinprep cytology test

LSIL:

Low-grade squamous intraepithelial lesion

HSIL:

High-grade squamous intraepithelial lesion

ICC:

Invasive cervical cancer

DNA:

Deoxyribo Nucleic Acid

OTUs:

Operational taxonomic units

AUC:

Area under the curve of ROC

KEGG:

Kyoto Encyclopaedia of Genes and Genomes database

FDR:

False discovery rate

OD value:

Optical density value

GWAS:

Whole Transcriptome Association Study

BCRP:

Breast cancer resistance protein

References

  1. Misra JS, Das V, Srivastava AN, Singh U. Chhavi: Role of different etiological factors in progression of cervical intraepithelial neoplasia. Diagn Cytopathol. 2006;34(10):682–5.

    Article  PubMed  Google Scholar 

  2. Syrjanen K. New concepts on risk factors of HPV and novel screening strategies for cervical cancer precursors. Eur J Gynaecol Oncol. 2008;29(3):205–21.

    CAS  PubMed  Google Scholar 

  3. Smith MI, Yatsunenko T, Manary MJ, Trehan I, Mkakosya R, Cheng J, et al. Gut microbiomes of Malawian twin pairs discordant for kwashiorkor. Science. 2013;339(6119):548–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Mitra A, MacIntyre DA, Lee YS, Smith A, Marchesi JR, Lehne B, et al. Cervical intraepithelial neoplasia disease progression is associated with increased vaginal microbiome diversity. Sci Rep. 2015;5:16865.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.

    Article  CAS  PubMed  Google Scholar 

  8. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Yildiz IZ, Usubutun A, Firat P, Ayhan A, Kucukali T. Efficiency of immunohistochemical p16 expression and HPV typing in cervical squamous intraepithelial lesion grading and review of the p16 literature. Pathol Res Pract. 2007;203(6):445–9.

    Article  CAS  PubMed  Google Scholar 

  10. Oh HY, Kim BS, Seo SS, Kong JS, Lee JK, Park SY, et al. The association of uterine cervical microbiota with an increased risk for cervical intraepithelial neoplasia in Korea. Clin Microbiol Infect. 2015;21(7):674 e671–9.

    Article  Google Scholar 

  11. Mitchell CM, Haick A, Nkwopara E, Garcia R, Rendi M, Agnew K, et al. Colonization of the upper genital tract by vaginal bacterial species in nonpregnant women. Am J Obstet Gynecol. 2015;212(5):611 e611–9.

    Article  Google Scholar 

  12. Witkin SS, Linhares IM, Giraldo P. Bacterial flora of the female genital tract: function and immune regulation. Best Pract Res Clin Obstet Gynaecol. 2007;21(3):347–54.

    Article  PubMed  Google Scholar 

  13. McMillan A, Dell M, Zellar MP, Cribby S, Martz S, Hong E, et al. Disruption of urogenital biofilms by lactobacilli. Colloids Surf B Biointerfaces. 2011;86(1):58–64.

    Article  CAS  PubMed  Google Scholar 

  14. Tseng CH, Lin JT, Ho HJ, Lai ZL, Wang CB, Tang SL, et al. Gastric microbiota and predicted gene functions are altered after subtotal gastrectomy in patients with gastric cancer. Sci Rep. 2016;6:20701.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Zhang J, Li S, Yan Q, Chen X, Yang Y, Liu X, et al. Interferon-beta induced microRNA-129-5p down-regulates HPV-18 E6 and E7 viral gene expression by targeting SP1 in cervical cancer cells. Plos One. 2013;8(12):e81366.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Doorbar J, Griffin H. Intrabody strategies for the treatment of human papillomavirus-associated disease. Expert Opin Biol Ther. 2007;7(5):677–89.

    Article  CAS  PubMed  Google Scholar 

  17. Hu Z, Zhu D, Wang W, Li W, Jia W, Zeng X, et al. Genome-wide profiling of HPV integration in cervical cancer identifies clustered genomic hot spots and a potential microhomology-mediated integration mechanism. Nat Genet. 2015;47(2):158–63.

    Article  CAS  PubMed  Google Scholar 

  18. Teng PN, Bateman NW, Darcy KM, Hamilton CA, Maxwell GL, Bakkenist CJ, et al. Pharmacologic inhibition of ATR and ATM offers clinically important distinctions to enhancing platinum or radiation response in ovarian, endometrial, and cervical cancer cells. Gynecol Oncol. 2015;136(3):554–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Santarpia L, Iwamoto T, Di Leo A, Hayashi N, Bottai G, Stampfer M, et al. DNA repair gene patterns as prognostic and predictive factors in molecular breast cancer subtypes. Oncologist. 2013;18(10):1063–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Bazarsad S, Kim JY, Zhang X, Kim KY, Lee DY, Ryu MH, et al. Ataxia-Telangiectasia-Mutated Protein Expression as a Prognostic Marker in Adenoid Cystic Carcinoma of the Salivary Glands. Yonsei Med J. 2018;59(6):717–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Sukowati CH, Rosso N, Pascut D, Anfuso B, Torre G, Francalanci P, et al. Gene and functional up-regulation of the BCRP/ABCG2 transporter in hepatocellular carcinoma. BMC Gastroenterol. 2012;12:160.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Wang F, Xue X, Wei J, An Y, Yao J, Cai H, et al. hsa-miR-520h downregulates ABCG2 in pancreatic cancer cells to inhibit migration, invasion, and side populations. Br J Cancer. 2010;103(4):567–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Gupta N, Martin PM, Miyauchi S, Ananth S, Herdman AV, Martindale RG, et al. Down-regulation of BCRP/ABCG2 in colorectal and cervical cancer. Biochem Biophys Res Commun. 2006;343(2):571–7.

    Article  CAS  PubMed  Google Scholar 

  24. Luo Y, Wang Q, Tian P, Jia Y. Highly expressed CHAF1A and PCNA are positively associated with malignancy of cervical squamous cell carcinoma. Xi Bao Yu Fen Zi Mian Yi Xue Za Zhi. 2017;33(12):1696–701.

    PubMed  Google Scholar 

  25. Lv Q, Zhang J, Yi Y, Huang Y, Wang Y, Wang Y, et al. Proliferating Cell Nuclear Antigen Has an Association with Prognosis and Risks Factors of Cancer Patients: a Systematic Review. Mol Neurobiol. 2016;53(9):6209–17.

    Article  CAS  PubMed  Google Scholar 

  26. Wang SC. PCNA: a silent housekeeper or a potential therapeutic target? Trends Pharmacol Sci. 2014;35(4):178–86.

    Article  PubMed  Google Scholar 

  27. Pang X, Zhang Y, Zhang S. High-mobility group box 1 is overexpressed in cervical carcinoma and promotes cell invasion and migration in vitro. Oncol Rep. 2017;37(2):831–40.

    Article  CAS  PubMed  Google Scholar 

  28. Chen J, Xi B, Zhao Y, Yu Y, Zhang J, Wang C. High-mobility group protein B1 (HMGB1) is a novel biomarker for human ovarian cancer. Gynecol Oncol. 2012;126(1):109–17.

    Article  CAS  PubMed  Google Scholar 

  29. Zhang J, Kou YB, Zhu JS, Chen WX, Li S. Knockdown of HMGB1 inhibits growth and invasion of gastric cancer cells through the NF-kappaB pathway in vitro and in vivo. Int J Oncol. 2014;44(4):1268–76.

    Article  CAS  PubMed  Google Scholar 

  30. Wu T, Zhang W, Yang G, Li H, Chen Q, Song R, et al. HMGB1 overexpression as a prognostic factor for survival in cancer: a meta-analysis and systematic review. Oncotarget. 2016;7(31):50417–27.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Vlahopoulos S, Adamaki M, Khoury N, Zoumpourlis V, Boldogh I. Roles of DNA repair enzyme OGG1 in innate immunity and its significance for lung cancer. Pharmacol Ther. 2019;194:59–72.

    Article  CAS  PubMed  Google Scholar 

  32. Kubo N, Morita M, Nakashima Y, Kitao H, Egashira A, Saeki H, et al. Oxidative DNA damage in human esophageal cancer: clinicopathological analysis of 8-hydroxydeoxyguanosine and its repair enzyme. Dis Esophagus. 2014;27(3):285–93.

    Article  CAS  PubMed  Google Scholar 

  33. Zeng X, Zhang Y, Yue T, Zhang T, Wang J, Xue Y, et al. Association between XRCC1 polymorphisms and the risk of cervical cancer: a meta-analysis based on 4895 subjects. Oncotarget. 2017;8(2):2249–60.

    Article  PubMed  Google Scholar 

  34. Abdel-Fatah T, Sultana R, Abbotts R, Hawkes C, Seedhouse C, Chan S, et al. Clinicopathological and functional significance of XRCC1 expression in ovarian cancer. Int J Cancer. 2013;132(12):2778–86.

    Article  CAS  PubMed  Google Scholar 

  35. He L, Yang H, Zhou S, Zhu H, Mao H, Ma Z, et al. Synergistic antitumor effect of combined paclitaxel with FEN1 inhibitor in cervical cancer cells. DNA Repair (Amst). 2018;63:1–9.

    Article  CAS  Google Scholar 

  36. He L, Zhang Y, Sun H, Jiang F, Yang H, Wu H, et al. Targeting DNA Flap Endonuclease 1 to Impede Breast Cancer Progression. EBioMedicine. 2016;14:32–43.

    Article  PubMed  PubMed Central  Google Scholar 

  37. He L, Luo L, Zhu H, Yang H, Zhang Y, Wu H, et al. FEN1 promotes tumor progression and confers cisplatin resistance in non-small-cell lung cancer. Mol Oncol. 2017;11(6):640–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Kohli RM, Zhang Y. TET enzymes, TDG and the dynamics of DNA demethylation. Nature. 2013;502(7472):472–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Yang L, Yu SJ, Hong Q, Yang Y, Shao ZM. Reduced Expression of TET1, TET2, TET3 and TDG mRNAs Are Associated with Poor Prognosis of Patients with Early Breast Cancer. Plos One. 2015;10(7):e0133896.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Zhong J, Jermusyk A, Wu L, Hoskins JW, Collins I, Mocci E, et al. A Transcriptome-Wide Association Study Identifies Novel Candidate Susceptibility Genes for Pancreatic Cancer. J Natl Cancer Inst. 2020;112(10):1003–12.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Tomkinson AE, Naila T, Khattri Bhandari S. Altered DNA ligase activity in human disease. Mutagenesis. 2020;35(1):51–60.

    Article  CAS  PubMed  Google Scholar 

  42. Chao XP, Sun TT, Wang S, Fan QB, Shi HH, Zhu L, et al. Correlation between the diversity of vaginal microbiota and the risk of high-risk human papillomavirus infection. Int J Gynecol Cancer. 2019;29(1):28–34.

    Article  PubMed  Google Scholar 

  43. Ruttkay-Nedecky B, Jimenez Jimenez AM, Nejdl L, Chudobova D, Gumulec J, Masarik M, et al. Relevance of infection with human papillomavirus: the role of the p53 tumor suppressor protein and E6/E7 zinc finger proteins (Review). Int J Oncol. 2013;43(6):1754–62.

    Article  CAS  PubMed  Google Scholar 

  44. Simoes RB, Campaner AB. Post-cervical conization outcomes in patients with high-grade intraepithelial lesions. APMIS. 2013;121(12):1153–61.

    Article  PubMed  Google Scholar 

  45. Roossink F, Wieringa HW, Noordhuis MG, ten Hoor KA, Kok M, Slagter-Menkema L, et al. The role of ATM and 53BP1 as predictive markers in cervical cancer. Int J Cancer. 2012;131(9):2056–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Yang H, Wang H, Chavan SS, Andersson U. High Mobility Group Box Protein 1 (HMGB1): The Prototypical Endogenous Danger Molecule. Mol Med. 2015;21(Suppl 1):S6–S12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Khunweeraphong N, Szollosi D, Stockner T, Kuchler K. The ABCG2 multidrug transporter is a pump gated by a valve and an extracellular lid. Nat Commun. 2019;10(1):5433.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Breiman L. Random forests. Mach Learn. 2001;45(1):5–32.

    Article  Google Scholar 

Download references

Funding

This study was supported by the Thousand-Hundred-Ten Talent Project of Guangxi Province.

(NO.2015226), the 139 Medical High-level Talents Training Plan of Guangxi Province (NO.G201903032) and the Medical and Health Appropriate Technology Development and Promotion Application Project of Guangxi Province (NO.S2018031) and the Special Budget for Key Laboratories of Guangxi Province (NO.GKZ201613) and the Key Laboratory of Early Prevention and Treatment of Regional High-incidence Tumors/Guangxi Key Laborator, Ministry of Education (NO.GXK201605).

Author information

Authors and Affiliations

Authors

Contributions

Li Li acted as the guarantor of this study and conceived the idea. He Wang collected the data and wrote this paper. Yanming Jiang and Yuejuan Liang performed the statistical analysis. Lingjia Wei and Wei Zhang performed the manuscript editing. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Li Li.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional committee and with the 1964 Helsinki Declaration and its later comparable ethical standards. This study was approved by the Ethics Committee of Guangxi Medical University Cancer Hospital. Written informed consent was obtained from the participants or their families.

Consent for publication

Informed consent was obtained from all individual participants included in the study.

Competing interests

No potential conflicts of interest relevant to this article were reported.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, H., Jiang, Y., Liang, Y. et al. Observation of the cervical microbiome in the progression of cervical intraepithelial neoplasia. BMC Cancer 22, 362 (2022). https://doi.org/10.1186/s12885-022-09452-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09452-0

Keywords