Analysis of the correlation between cervical HPV infection, cervical lesions and vaginal microecology

Background: Vaginal microbiota is involved in human papillomavirus (HPV) infection and cervical cancer (CC) progression, and the specific changes in vaginal microbial composition during this process remains uncertain.

Objective: This study aimed to observe the changes in the specific composition of vaginal microorganisms in different cervical lesions and identify biomarkers at different stages of lesions.

Methods: In this study we used the illumina high-throughput gene sequencing technology to determine the V4 region of 16SrRNA and observed the vaginal microbial composition in different cervical lesions.

Results: The vaginal microbiota of patients with high-risk HPV infection and cervical lesions is significantly different from that of the normal population, but there is no significant difference in the richness of vaginal microbes. The diversity of vaginal species in CC patients is higher than that in high-risk HPV infection or CIN patients. The main manifestation is an increase in the diversity of vaginal microbes, a decrease in the relative abundance of cyanobacteria and Lactobacillus, and an increase in the relative abundance of dialister, peptonephila and other miscellaneous bacteria. There are characteristic vaginal biomarker in normal women, high risk HPV patients and CC patients. In detail, the biomarker in the normal group was varibaculum, the biomarker in the high-risk HPV group was saccharopolyspora, the biomarker of the CC group was the Proteobacteria, Corynebacterium, Coprococcus, Peptococcus and Ruminococcus.

Conclusions: The study indicated that the compositions of vaginal microbes in different cervical lesions is different. The vaginal microbial composition has a certain diagnostic effect on healthy women, patients with high-risk HPV infection and cervical lesions. These microbes may serve as potential biomarkers for CC. It also provided an effective way for the treatment of HPV infections and cervical lesions.

1 Background

Cervical cancer (CC) is one of the most common malignancies in women which associated with HPV infection in different regions and populations (Crow, 2012). HPV infection is a common female reproductive tract disease, 70% of women may be infected with this disease in their lifetime (Xia et al., 2023), but so far there is still a lack of effective therapy. How to effectively control the replication of HPV virus has become a challenge for clinicians.

The human microbiota is considered to be the “second human genome” (Integrative HMP (iHMP) Research Network Consortium, 2019). As an important part of microbiota, vaginal microbiome plays a vital role in preventing various infectious diseases of the reproductive tract. Recent studies have shown that there is a relatively relationship between vaginal microbiome, HPV infection and cervical lesions (Chao et al., 2019).

The vaginal microecology is a complex and dynamic system, mainly composed of related anatomical structures, human immunity, internal and external environments, and vaginal microbiome. The normal vaginal micro-ecosystem is parasitized by a variety of bacteria such as lactobacillus, streptococcus, and a various microorganisms such as mycoplasma and viruses, and jointly resist the invasion of external pathogenic bacteria (Liang et al., 2019). When the vaginal microecological balance is broken, it will cause a decrease in the number of lactobacilli, an increase in the number of other miscellaneous bacteria, anaerobic bacteria, and changes in vaginal microbial diversity (Borgdorff et al., 2016), which may lead to a series of reproductive tract diseases. A meta-analysis based on related articles published from 2003 to 2017 showed that among sexually active women, those with vaginal microecological imbalance are more likely to develop persistent HPV infections and increase the risk of cervical intraepithelial neoplasia (CIN) (Brusselaers et al., 2019). Conversely, patients with HPV infection and cervical lesions are also more likely to develop vaginal infectious diseases. The essence of HPV infection is an infectious disease of the reproductive tract. In addition to symptomatic treatment of pathogenic bacteria, attention should be paid to the microecological balance. Therefore, restoring the balance of micro-ecology may be another important way to prevent and treat CC in the future. V4 region is one of the 9 variable regions of 16SrRNA which has a moderate length and high coverage of bacteria. It has been widely used in microbial sequencing and analysis in recent years (Caporaso et al., 2012). In this study we used the illumina high-throughput gene sequencing technology to determine the V4 region of 16SrRNA and observed the vaginal microbial composition in different cervical lesions. The purpose is to explore the causes of HPV infection and cervical lesions, which will provide relevant theoretical basis for the development of vaginal microecological therapy and microecological preparations, and provide new methods for the diagnosis and treatment of HPV infections and cervical lesions.

2 Materials and methods

2.1 Patients

This study was approved by the Ethics Committee of Suzhou Sci-tech City Hospital. According to the Declaration of Helsinki, written informed consent was obtained from all patients.

All patients were screened strictly according to the following criteria: (1) Having a history of sexual life; (2) Age 20-65 years old; (3) Race is a Han nationality; (4) No history of sexual life, vaginal treatment and flushing 3 days before sampling; (5) No history of antibiotic use within 2 weeks before sampling; (6) No history of radiotherapy and chemotherapy; (7) No history of cervical and vaginal surgery; (8) No sexually transmitted diseases; (9) No abnormal vaginal bleeding or vaginal prolapse and other vaginal diseases; (10) No long-term use of sex hormones, immunosuppressants and steroids; (11) No immune diseases, other malignant tumors and mental diseases; (12) Not during menstruation, pregnancy or lactation. If any criterion is not met, it cannot be included in our study.

125 patients were enrolled in our study who underwent HPV, Thinprep cytologic test (TCT), and colposcopy in the Second Hospital of Dalian Medical University and Suzhou Sci-tech City Hospital from June 2020 to December 2020. HPV typing was performed by PCR and reverse point hybridization which can detect 28 HPV types, including high-risk HPV: HPV16, 18, 31, 33, 35, 39, 45, 59, 52, 53, 56, 58, 59, and 66, common low-risk HPV: HPV68, CP8304, 6, 11, 43, 73, 82, 26, 40, 42, 44, 54, 61, and 83. According to the HPV and TCT results, they were divided into the normal group (n=27, N: HPV and TCT are negative) and the high-risk HPV group (n=40, H: high-risk HPV positive for at least 1 year, TCT negative). Based on postoperative pathological results, patients with shunt colposcopy, colposcopy biopsy or cervical conization were divided into CIN group (n=40, C: CIN2 or 3) and CC group (n=18, Ca: cervical cancer).

The types of HPV infection, age, menopausal history, pregnancy, birth history, and contraceptive methods in each group of patients were analyzed. Using 16SrRNA-V4 region gene amplification technology and Illlumina high-throughput sequencing technology to detect vaginal secretions retained in the GUHE Flora Storage buffer (Zhejiang Hangzhou Equipment Preparation 20190682, GUHE Laboratories, Hangzhou, China), the composition of the vaginal microbiota in each sample is verified. Using bioinformatics analysis methods to analyze the composition of the microbes, diversity of vaginal microbes, microbial markers and microbial predictive ability in each group.

2.2 Sampling method

The patients were maintained the lithotomy position, a sterile vaginal speculum without lubricant was used to gently and fully expose the cervix, rotate the sterile cotton swab on the upper 1/3 section of the vagina and the lower section of the vagina for 10-15 seconds, put the cotton swab with vaginal secretion into the storage solution (Zhejiang Hangzhou equipment preparation 20190682), stir and rinse for 30 seconds, discard the cotton swab and mix the storage solution, the specimen bottles were stored in the refrigerator at - 20°C.

2.3 DNA extraction

DNA samples from vaginal secretions were extracted using DNA isolation kit (GUHE Laboratories, Hangzhou, China). Agarose gel electrophoresis and Nano Drop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) were used to measure and detect DNA.

2.4 16SrRNA amplicon pyrosequencing

The forward primer 515F(5’-GTGCCAGCMGCCGCGGTAA-3’) and the reverse primer 806R(5’-GGACTACHVGGGTWTCTAAT-3’) were used for PCR amplification of the bacterial 16SrRNA gene V4 region. The PCR components contained Phusion High-Fidelity PCR Master Mix (25 μl), which include DNA Template (10 μl), Forward primer(3 μl), Reverse primer(3 μl), DMSO(3 μl), and ddH2O(6 μl). Thermal cycling includes initial denaturation for 30 s (98°C), 25 cycles of denaturation for 15 s (98°C), annealing for 15 s (58°C), extension for 15 s (72°C), final extension for 60 s (72°C). PCR amplicons were purified and quantified using AgencourtAMPure XP Beads (Beckman Coulter, Indianapolis, IN, USA) and the PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA). The Illumina NovoSeq6000 platform (GUHE Info technology Co., Ltd, Hangzhou, China) was used for paired-end 2×150 bp sequencing.

2.5 Sequence analysis

Sequencing data were processed by Quantitative Insights Into Microbial Ecology (QIIME, v1.9.1) (Caporaso et al., 2010). Low-quality sequences were filtered by specific criteria (Chen and Jiang, 2014).Operational taxonomic unit(OTU) was selected by Vsearch v2.4.4 (Rognes et al., 2016) and classified by VSEARCH to search representative sequences in the SILVA132 database (Quast et al., 2013). To prevent environmental contamination, all samples in the project were processed side by side using the same batches of reagents, at the same time, we have added a sample with a control check and a native sample of a standard strain. Decontam was used to remove contamination sequences.

2.6 Bioinformatics and statistical analysis

Sequence data were analyzed by R packages (v3.2.0) and QIIME. The OTU table was used for the calculation of OTU-level alpha diversity and generation of OTU-level ranked abundance curves. Taxa abundance at the phylum, class, order, family, genus and species levels were compared among samples or groups by Kruskal test from R stats package. Kruskal.test and Wilcox.testwere were used to analyze the alpha diversity of each group. Anosim test was used to determine whether the grouping was meaningful. Linear discriminant analysis was used for effect size (LEfSe) to detect significantly different taxa across groups (Segata et al., 2011).

2.7 Data analysis

SPSS 21.0 (SPSS, Chicago, Illinois, USA) was used for statistical analysis of clinical data. Continuous variables and categorical variables was analyzed by the t-test and χ2 test or Fisher’s test separately. P values <0.05 were considered to be statistically significant.

3 Results

3.1 Clinical data

97 of the 125 patients, were infected with HPV. In those 97 HPV patients, the most HPV subtype is type 16, accounting for 47.4% (n=46), followed by type 52 (22.7%) (n=22), type 58 (16.5%) (n=16) and type 18 (12.4%) (n=12) respectively. In the high-risk HPV group, HPV16 and/or 18 infections accounted for 52.5% (n=21). In the CIN group, HPV16 and/or 18 infections accounted for 50% (n=20), and 1 case without HPV infection, accounting for 2.5%. In the CC group, HPV16 and/or 18 infection accounted for 94.4% (n=17). The ages of the enrolled patients were 21 to 65 years old. There were no statistical differences in general clinical data among the groups including age, menopausal history, fertility status, contraceptive methods and BMI (Table 1).

Table 1

Table 1. Clinical data.

3.2 Microbial composition

The Relative abundance of phyla in four groups was shown in Figure 1. The Phylum Firmicutes occupied the highest proportion. The relative abundances of Firmicutes in the normal group, high-risk HPV group, CIN group and CC group were 60.1%, 68.7%, 68.6%, and 59.2%, respectively and there was no significant difference among the groups. The other top 10 phyla in vagina are Actinobacteria, Bacteroides, Proteobacteria, Fusobacteria, Tenericutes, Chlamydiae, TM7, Thermi, and Cyanobacteria (Figure 1A). The Wilcoxon test results showed that compared to the normal group, the relative abundance of Cyanobacteria and Proteobacteria in the CIN group was significantly decreased, and that of Thermi was significantly increased, while that of Cyanobacteria and Thermi in the CC group was significantly decreased. Compared with the CIN group, the relative abundance of Fusobacteria in the CC group was significantly decreased while that of Proteobacteria was significantly increased, In addition, the relative abundance of Cyanobacteria in the CIN group was significantly lower than that of high-risk HPV group (Figures 1B–E).

Figure 1

Figure 1. Comparative analysis of the relative abundance of phyla in four groups. (A) Bar chart of relative abundance of top 10 phyla in four groups. (B) Histogram of differential phylum between Group N and Group C. (C) Histogram of differential phylum between Group N and Group Ca. (D) Histogram of differential phylum between Group H and Group C. (E) Histogram of differential phylum between Group C and Group Ca.

The Relative abundance of genera in all groups was shown in Figure 2. The Lactobacillus occupied the highest proportion. The relative abundances in the normal group, high-risk HPV group, CIN group, and CC group were 48.9%, 57.4%, 53.8%, 33.6% and the remaining top 10 bacterial genera included Prevotella, Megasphaera, Streptococcus, Sneathia, Dialister, Anaerococcus, Ureaplasma, Peptoniphilus, and Mycoplasma (Figure 2A). The Wilcoxon test results showed that there were significant differences in bacterial genera among groups. For the top 10 bacterial genera, there are significant differences between the high-risk HPV group and the CC group or between the CIN group and the CC group. Compared with the high-risk HPV group, the relative abundance of Lactobacillus in the CC group was significantly reduced while that of Dialister was significantly increased. Compared with the CIN group, the relative abundance of Lactobacillus in the CC group was also significant reduced, while that of Dialister and Peptoniphiluswere significantly increased (Figures 2B–G).

Figure 2

Figure 2. Comparative analysis of the relative abundance of genera in four groups. (A) Bar chart of relative abundance of top 10 genera in four groups. (B) Histogram of differential genera between Group N and Group H. (C) Histogram of differential genera between Group N and Group C. (D) Histogram of differential genera between Group N and Group Ca. (E) Histogram of differential genera between Group H and Group C. (F) Histogram of differential genera between Group H and Group Ca. (G) Histogram of differential genera between Group C and Group Ca.

3.3 Microbial Alpha diversity

First, we described the Alpha diversity of the bacterial community by measuring Chao1 index and Shannon index. Second, we created the rarefaction curves to test the sequencing depth and the result showed that the amount of sampling is reasonable (Figure 3A). We found the Shannon index of the CC group was significantly higher than that of the HPV group (P=0.034) and the CIN group (P=0.019) (Figure 3B). But there was no significant difference of the Chao1 index between groups (P>0.05) (Figure 3C).

Figure 3

Figure 3. The microbial Alpha diversity analysis in different groups. (A) Shannon-Winner curve in different groups. (B) Shannon index box chart for each group. (C) Chao1 index box chart for each group. (*P<0.05).

3.4 Biomarkers

In order to fully determine the role of vaginal microbiota and the differences in microbial between the groups, we performed the Anosim test. The results showed that there were significant differences between the groups(P=0.044) (Figure 4). Furthermore, the composition of vaginal microbes in each group was analyzed by LEfSe, and Figure 5 showed the clade diagram and bar graph of LDA value (only microbes with LDA>2).In this way, biomarkers (species with significant differences between groups) were screened between groups. The difference markers between groups selected by LEfSe are the dominant species, and these markers are distinguished by specific significant differences between this group and other groups. At the phylum level, the biomarker of the CC group was the Proteobacteria, and at the genus level, the biomarkers of the CC group were Corynebacterium (LDA=3.16), Coprococcus (LDA=2.27), Peptococcus (LDA=2.25) and Ruminococcus (LDA=2.08), the biomarker in the high-risk HPV group was saccharopolyspora (LDA=2.04), the biomarker in the normal group was varibaculum (LDA=2.34), and no biomarker have been found in the CIN group (Figure 5).

Figure 4

Figure 4. Anosim test box diagram of each group.

Figure 5

Figure 5. Biomarkers of each group. (A) The evolutionary branch diagram of microorganism with significant differences in each group. (B) LDA value histogram of microorganism with significant differences in each group.

3.5 Microbial prediction ability

In this study, 50 different bacterial communities were screened, a random forest graph was made (Figure 6), and an ROC curve was plotted (Figure 7). The results showed that the vaginal microbial composition had a predictive effect on each group. The prediction accuracy of HPV group (AUC=0.51) and CC group (AUC=0.68) was lower than that of the normal group (AUC=0.75) and CIN group (AUC=0.71).

Figure 6

Figure 6. Random-forest graph of top 50 different bacterial communities between each group.

Figure 7

Figure 7. ROC curve of each group.

4 Discussion

The female reproductive tract is a channel that communicates with the outside world which is easily invaded by external microorganisms, leading to HPV infections and other infectious diseases. The infection of high-risk HPV subtypes may lead to the CC (Yang et al., 2020). Consistent with the previous study, this study showed that more than 90% of patients with CIN or cervical cancer are infected with high-risk HPV. Among all HPV subtypes, types 16 and 18 are considered to have the highest carcinogenic effects (de Sanjosé et al., 2007) (Torre et al., 2015). However, most HPV infections subside within a short period of time, and only a few patients will gradually cause CC after prolonged infection. It is known that the vaginal microecology has played an important role in the process of HPV infections-induced cervical lesions (Kovachev, 2020). Vaginal microecology is affected by the factors such as age, menstruation, medication, vaginal lavage, contraceptive methods and other factors. It is a dynamic equilibrium system formed by the mutual restriction and coordination of microorganisms and the environment in the host (Mitra et al., 2016). If the vaginal microecological balance is broken and the immune system is compromised, foreign microorganisms will invade the reproductive system and reduce the immune system’s clearance rate of HPV and other pathogens (Chehoud et al., 2017).

Firmicutes is the dominant phyla and Lactobacillus is the dominant genera for normal vaginal microbiome. Similar to previous studies (Walther-António et al., 2016; Chao et al., 2019; Norenhag et al., 2020; So et al., 2020; Xie et al., 2020), our study found that Firmicutes did not differ significantly between the groups. Compared with the cervical lesions group, Lactobacillus were significantly reduced in the CC group, while the relative abundance of Dialister and Peptoniphilus was significantly increased. It has been suggested that pathogenic mechanism of Peptoniphilus may be related to its strong adhesion to the vaginal epithelium and the formation of biofilms (Patterson et al., 2010), while the pathogenic mechanism of Lactobacillus may vary with different strains (Ravel et al., 2011; Teka et al., 2023). The research on Proteobacteria is inconclusive (Chao et al., 2019; Cheng et al., 2020). Recent study has found that Proteobacteria in colon cancer are positively correlated with tumor burden (Lee et al., 2021). Our study also showed that Proteobacteria in CC patients has significantly increased compared with CIN patients. It indicated that Proteobacteria may be involved in tumor progression. Fusobacteria is another one of the common bacterial phyla in the vagina. Harrandahet al. found that Fusobacteria potentially enhanced the invasiveness and survival rate of cancer cells in the oral cavity (Harrandah et al., 2020), but our results demonstrated that the relative abundance of Fusobacteria has decreased significantly in the development of cervical lesions, The reasons for this difference may due to the selection of different cases, different geographic region and microenvironment. Therefore, it needs to be further explored.

There are few studies on the distribution and function of cyanobacteria in human microecology all over the world. Previous studies have shown that cyanobacteria are one of the common intestinal microorganisms of small intestinal adenocarcinoma (Pan et al., 2021). Cyanobacteria are also described in the breast. Some studies show that the relative abundance of cyanobacteria in Chinese women is significantly higher than that in Slovak women and it indicates that cyanobacteria are related to the geographical area. The relative abundance of cyanobacteria in primary breast cancer is lower than that of healthy women in China (Hadzega et al., 2021) and this result similar to ours. Although its mechanism is unclear, it suggests that cyanobacteria may play an inhibitory role in tumor progression. However, there are few studies on the physiological role and mechanism of cyanobacteria in the vaginal. Our findings may lay a foundation for future research. Previous studies have shown that thermi may be related to HPV16 (Huang et al., 2018). Other studies have shown that the relative abundance of thermi is higher in advanced lung cancer than in early lung cancer (Yu et al., 2016). Our study also suggests that the relative abundance of thermi in CIN group is high. It indicates that thermi may be related to tumors development. However, the specific mechanism remains to be further explored.

In recent years, the research on the diversity of the vaginal microbiome has gradually become a hot topic, and Alpha diversity is the main indicator to measure the diversity and richness of microbial communities. The diversity index is mainly used to measure the heterogeneity of the community, and the richness is mainly used to measure the number of species in a single sample. A study on the change of vaginal microecological under different cervical lesions showed that patients with normal cervix had a higher vaginal microbiome richness than patients with cervical lesions, but the difference was not statistically significant (Wu et al., 2020). Another study has also shown that there is no significant difference in the abundance and diversity of vaginal bacteria under different cervical lesions, but HPV infection significantly increases the abundance and diversity of vaginal bacteria (Chen et al., 2020). Mitra et al. found that following the degree of cervical lesions increases, the richness and diversity of vaginal bacteria also increase (Teka et al., 2023). In current study, the Shannon index was used to estimate the microbial diversity, and the Chao1 index was used to estimate the microbial richness. The results showed that there is no significant difference in microbial richness but the diversity of vaginal microbiome in patients with CC was significantly higher than that in patients with high-risk HPV infection and CIN patients. These different results in various studies on the diversity of vaginal microbiome may be related to different geographic regions, races, sample numbers, and sample collection techniques (Kovachev, 2020). A study involving almost all ethnic females showed that the risk of increased microbiota diversity and vaginal microecological imbalance in Caribbean African women is four times than that of European/Caucasian and African women (Mortaki et al., 2020), indicating that race is a key factor affecting the vaginal bacterial community. All the research objects of this study are from the Han nationality, it only represent the characteristics of vaginal microorganisms of the Han nationality. Based on the above studies, a worldwide multi-center study may be quite necessary.

It has been believed that there are characteristic biomarkers in patients with HPV infection and different degrees of cervical lesions, the biomarkers derived from the results of various studies are different. In recent years, studies found that Delftia and Gardnerella may be biomarkers of cervical precancerous lesions (Usyk et al., 2020; Wu et al., 2020) and bacteria such as Gardnerella vaginalis, Gardnerella, Peptostreptococcusanaerobius, Mobiluncuscurtisii, Prevotella genus and Fusobacterium nucleatum may be characteristic biomarkers of HPV16-positive patients (Yang et al., 2020). This study used LEfSe analysis to explore the characteristic biomarkers and showed that Proteobacteria, Corynebacterium, Coprococcus, Peptococcus and Ruminococcus may be used as biomarkers for CC patients. The biomarker for high-risk HPV patients is Saccharopolyspora, and the biomarker for normal patients is varibaculum. These characteristic bacterial genera are not only related to the degree of different cervical lesions, but may also be related to certain reproductive tract infections. Previous studies have shown that Proteobacteria (Wei et al., 2023), Corynebacterium (Wang et al., 2014) and Peptococcus (Lefèvre et al., 1985) are related to bacterial vaginosis. The relationship between CC and bacterial vaginosis is also worthy of further study. Varibaculum is a kind of actinomycete and previous studies have shown that its specific strain may exist in the intestine of premature infants (Brown et al., 2013). And a data searched in HMP showed that Varibaculum in healthy human microecology were mainly distributed in the vagina (Integrative HMP (iHMP) Research Network Consortium, 2019) and it is consistent with our results. Recently, a new strain of varibaculum has been extracted from the vagina of healthy Senegalese women (Fall et al., 2019). Therefore, the distribution and mechanism of varibaculum and its subordinate species may be a new research field of vaginal microecology. The composition of the vaginal microbiome is different between patients with CIN1 and CIN2, and this different composition has extremely high accuracy in predicting the severity of CIN (Lee et al., 2020). Our results also demonstrated that the composition of vaginal microbes may predict high-risk HPV infection and cervical lesions.

5 Conclusion

This study found that the vaginal microecology is different under various pathological conditions of the cervix, but no difference in the richness of vaginal microbes. The diversity of vaginal species in CC patients is mainly manifested as a reduction in the number of Cyanobacteria and Lactobacillus, accompanied by an increase in the number of Dialister, Peptoniphilus and other miscellaneous bacteria. Different cervical lesions have specific biomarkers, so the composition of the vaginal microecological flora has a diagnostic effect on cervical lesions. These results may provide new evidences for the diagnosis of HPV infection and cervical lesions. Due to the changes in the vaginal microecology under cervical lesions, restoring the balance of vaginal microecology may be an effective way for the treatment of HPV infection and cervical lesions.

Data availability statement

The raw data used to support the findings of this study are provided in the following links: https://www.jianguoyun.com/p/DdRG93UQ_svSDBjQidkFIAA.

Ethics statement

The studies involving humans were approved by The Ethical Committee of Suzhou Hospital, Affiliated Hospital of Medical School, Nanjing University approved the study protocol under the No. of IRB202110002RI. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.

Author contributions

ZF: Conceptualization, Investigation, Project administration, Supervision, Validation, Writing – review & editing, Writing – original draft. DH: Data curation, Investigation, Writing – review & editing. XF: Project administration, Resources, Supervision, Validation, Writing – review & editing. YZ: Data curation, Investigation, Visualization, Writing – original draft, Writing – review & editing. LZ: Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by Suzhou science and technology livelihood project (Grant number SYSD2020083).

Acknowledgments

The authors are grateful to the GUHE Info technology, Co., Ltd., (Hangzhou, China) for assistance with RNA sequencing and delivering the data to us.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Borgdorff, H., Armstrong, S. D., Tytgat, H. L., Ndayisaba, G. F., Wastling, J. M., van de Wijgert, J. H, et al. (2016). Unique Insights in the Cervicovaginal Lactobacillus iners and L. crispatus Proteomes and Their Associations with Microbiota Dysbiosis. PloS One 11, e0150767. doi: 10.1371/journal.pone.0150767

PubMed Abstract | Crossref Full Text | Google Scholar

Brown, C. T., Sharon, I., Thomas, B. C., Castelle, C. J., Morowitz, M. J., Banfield, J. F. (2013). Genome resolved analysis of a premature infant gut microbial community reveals a Varibaculum cambriense genome and a shift towards fermentation-based metabolism during the third week of life. Microbiome. 1, 30. doi: 10.1186/2049-2618-1-30

PubMed Abstract | Crossref Full Text | Google Scholar

Brusselaers, N., Shrestha, S., van de Wijgert, J., Verstraelen, H. (2019). Vaginal dysbiosis and the risk of human papillomavirus and cervical cancer: systematic review and meta-analysis. Am. J. Obstet Gynecol. 221, 9–18.e8. doi: 10.1016/j.ajog.2018.12.011

PubMed Abstract | Crossref Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | Crossref Full Text | Google Scholar

Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Huntley, J., Fierer, N., Owens, S. M., et al. (2012). Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624. doi: 10.1038/ismej.2012.8

PubMed Abstract | Crossref Full Text | Google Scholar

Chao, X. P., Sun, T. T., Wang, S., Fan, Q. B., Shi, H. H., Zhu, L., et al. (2019). Correlation between the diversity of vaginal microbiota and the risk of high-risk human papillomavirus infection. Int. J. Gynecol Cancer. 29, 28–34. doi: 10.1136/ijgc-2018-000032

PubMed Abstract | Crossref Full Text | Google Scholar

Chehoud, C., Stieh, D. J., Bailey, A. G., Laughlin, A. L., Allen, S. A., McCotter, K. L., et al. (2017). Associations of the vaginal microbiota with HIV infection, bacterial vaginosis, and demographic factors. AIDS. 31, 895–904. doi: 10.1097/QAD.0000000000001421

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, H., Jiang, W. (2014). Application of high-throughput sequencing in understanding human oral microbiome related with health and disease. Front. Microbiol. 5, 508. doi: 10.3389/fmicb.2014.00508

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, Y., Qiu, X., Wang, W., Li, D., Wu, A., Hong, Z., et al. (2020). Human papillomavirus infection and cervical intraepithelial neoplasia progression are associated with increased vaginal microbiome diversity in a Chinese cohort. BMC Infect. Dis. 20, 629. doi: 10.1186/s12879-020-05324-9

PubMed Abstract | Crossref Full Text | Google Scholar

Cheng, W., Xu, F., Gao, L., Liu, J. (2020). The correlation between the determination of vaginal micro-ecological composition and the outcome of HPV infection by high-throughput metagene sequencing information technology on the illumina platform. J. Infect. Public Health 13, 1961–1966. doi: 10.1016/j.jiph.2020.05.024

PubMed Abstract | Crossref Full Text | Google Scholar

Crow, J. M. (2012). HPV: The global burden. Nature. 488, S2–S3. doi: 10.1038/488S2a

PubMed Abstract | Crossref Full Text | Google Scholar

de Sanjosé, S., Diaz, M., Castellsagué, X., Clifford, G., Bruni, L., Muñoz, N., et al. (2007). Worldwide prevalence and genotype distribution of cervical human papillomavirus DNA in women with normal cytology: a meta-analysis. Lancet Infect. Dis. 7, 453–459. doi: 10.1016/S1473-3099(07)70158-5

PubMed Abstract | Crossref Full Text | Google Scholar

Fall, N. S., Lo, C. I., Fournier, P. E., Sokhna, C., Raoult, D., Fenollar, F., et al. (2019). Arcanobacterium ihumii sp. nov., Varibaculum vaginae sp. nov. and Tessaracoccus timonensis sp. nov., isolated from vaginal swabs from healthy Senegalese women. New Microbes New Infect. 31, 100585. doi: 10.1016/j.nmni.2019.100585

PubMed Abstract | Crossref Full Text | Google Scholar

Hadzega, D., Minarik, G., Karaba, M., Kalavska, K., Benca, J., Ciernikova, S., et al. (2021). Uncovering microbial composition in human breast cancer primary tumour tissue using transcriptomic RNA-seq. Int. J. Mol. Sci. 22, 9058. doi: 10.3390/ijms22169058

PubMed Abstract | Crossref Full Text | Google Scholar

Harrandah, A. M., Chukkapalli, S. S., Bhattacharyya, I., Progulske-Fox, A., Chan, E. (2020). Fusobacteria modulate oral carcinogenesis and promote cancer progression. J. Oral. Microbiol. 13, 1849493. doi: 10.1080/20002297.2020.1849493

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, X., Li, C., Li, F., Zhao, J., Wan, X., Wang, K. (2018). Cervicovaginal microbiota composition correlates with the acquisition of high-risk human papillomavirus types. Int. J. Cancer. 143, 621–634. doi: 10.1002/ijc.31342

PubMed Abstract | Crossref Full Text | Google Scholar

Integrative HMP (iHMP) Research Network Consortium (2019). The integrative human microbiome project. Nature 569, 641–648. doi: 10.1038/s41586-019-1238-8

PubMed Abstract | Crossref Full Text | Google Scholar

Kovachev, S. M. (2020). Cervical cancer and vaginal microbiota changes. Arch. Microbiol. 202, 323–327. doi: 10.1007/s00203-019-01747-4

PubMed Abstract | Crossref Full Text | Google Scholar

Lee, J. G., Lee, Y. R., Lee, A. R., Park, C. H., Han, D. S., Eun, C. S. (2021). Role of the global gut microbial community in the development of colitis-associated cancer in a murine model. BioMed. Pharmacother. 135, 111206. doi: 10.1016/j.biopha.2020.111206

PubMed Abstract | Crossref Full Text | Google Scholar

Lee, Y. H., Kang, G. U., Jeon, S. Y., Tagele, S. B., Pham, H. Q., Kim, M. S., et al. (2020). Vaginal microbiome-based bacterial signatures for predicting the severity of cervical intraepithelial neoplasia. Diagnostics (Basel). 10, 1013. doi: 10.3390/diagnostics10121013

PubMed Abstract | Crossref Full Text | Google Scholar

Lefèvre, J. C., Jean, M., Averous, S., Viraben, R., Blanc, C., Bauriaud, R., et al. (1985). [Etiology of bacterial vaginosis (non-specific vaginitis)]. J. Gynecol Obstet Biol. Reprod. (Paris). 14, 703–708.

PubMed Abstract | Google Scholar

Liang, Y., Chen, M., Qin, L., Wan, B., Wang, H. (2019). A meta-analysis of the relationship between vaginal microecology, human papillomavirus infection and cervical intraepithelial neoplasia. Infect. Agent Cancer. 14, 29. doi: 10.1186/s13027-019-0243-8

PubMed Abstract | Crossref Full Text | Google Scholar

Mitra, A., MacIntyre, D. A., Marchesi, J. R., Lee, Y. S., Bennett, P. R., Kyrgiou, M. (2016). The vaginal microbiota, human papillomavirus infection and cervical intraepithelial neoplasia: what do we know and where are we going next. Microbiome. 4, 58. doi: 10.1186/s40168-016-0203-0

PubMed Abstract | Crossref Full Text | Google Scholar

Mortaki, D., Gkegkes, I. D., Psomiadou, V., Blontzos, N., Prodromidou, A., Lefkopoulos, F., et al. (2020). Vaginal microbiota and human papillomavirus: a systematic review. J. Turk Ger Gynecol Assoc. 21, 193–200. doi: 10.4274/jtgga.

PubMed Abstract | Crossref Full Text | Google Scholar

Norenhag, J., Du, J., Olovsson, M., Verstraelen, H., Engstrand, L., Brusselaers, N. (2020). The vaginal microbiota, human papillomavirus and cervical dysplasia: a systematic review and network meta-analysis. BJOG. 127, 171–180. doi: 10.1111/1471-0528.15854

PubMed Abstract | Crossref Full Text | Google Scholar

Pan, H., Cheng, H., Wang, H., Ge, W., Yuan, M., Jiang, S., et al. (2021). Molecular profiling and identification of prognostic factors in Chinese patients with small bowel adenocarcinoma. Cancer Sci. 112 (11), 4758–4771. doi: 10.1111/cas.15119

PubMed Abstract | Crossref Full Text | Google Scholar

Patterson, J. L., Stull-Lane, A., Girerd, P. H., Jefferson, K. K. (2010). Analysis of adherence, biofilm formation and cytotoxicity suggests a greater virulence potential of Gardnerella vaginalis relative to other bacterial-vaginosis-associated anaerobes. Microbiol. (Reading). 156, 392–399. doi: 10.1099/mic.0.034280-0

Crossref Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | Crossref Full Text | Google Scholar

Ravel, J., Gajer, P., Abdo, Z., Schneider, G. M., Koenig, S. S., McCulle, S. L., et al. (2011). Vaginal microbiome of reproductive-age women. Proc. Natl. Acad. Sci. U S A. 108 Suppl 1, 4680–4687. doi: 10.1073/pnas.1002611107

PubMed Abstract | Crossref Full Text | Google Scholar

Rognes, T., Flouri, T., Nichols, B., Quince, C., Mahé, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ. 4, e2584. doi: 10.7717/peerj.2584

PubMed Abstract | Crossref Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12, R60. doi: 10.1186/gb-2011-12-6-r60

PubMed Abstract | Crossref Full Text | Google Scholar

So, K. A., Yang, E. J., Kim, N. R., Hong, S. R., Lee, J. H., Hwang, C. S., et al. (2020). Changes of vaginal microbiota during cervical carcinogenesis in women with human papillomavirus infection. PloS One 15, e0238705. doi: 10.1371/journal.pone.0238705

PubMed Abstract | Crossref Full Text | Google Scholar

Teka, B., Yoshida, K., Firdawoke, E., Chanyalew, Z., Gizaw, M., Addissie, A., et al. (2023). Cervicovaginal microbiota profiles in precancerous lesions and cervical cancer among Ethiopian women. Microorganisms. 11, 833. doi: 10.3390/microorganisms11040833

PubMed Abstract | Crossref Full Text | Google Scholar

Torre, L. A., Bray, F., Siegel, R. L., Ferlay, J., Lortet-Tieulent, J., Jemal, A. (2015). Global cancer statistics, 2012. CA Cancer J. Clin. 65, 87–108. doi: 10.3322/caac.21262

PubMed Abstract | Crossref Full Text | Google Scholar

Usyk, M., Zolnik, C. P., Castle, P. E., Porras, C., Herrero, R., Gradissimo, A., et al. (2020). Cervicovaginal microbiome and natural history of HPV in a longitudinal study. PloS Pathog. 16, e1008376. doi: 10.1371/journal.ppat.1008376

PubMed Abstract | Crossref Full Text | Google Scholar

Walther-António, M. R., Chen, J., Multinu, F., Hokenstad, A., Distad, T. J., Cheek, E. H., et al. (2016). Potential contribution of the uterine microbiome in the development of endometrial cancer. Genome Med. 8, 122. doi: 10.1186/s13073-016-0368-y

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, B., Xiao, B. B., Shang, C. G., Wang, K., Na, R. S., Nu, X. X., et al. (2014). Molecular analysis of the relationship between specific vaginal bacteria and bacterial vaginosis metronidazole therapy failure. Eur. J. Clin. Microbiol. Infect. Dis. 33, 1749–1756. doi: 10.1007/s10096-014-2128-5

PubMed Abstract | Crossref Full Text | Google Scholar

Wei, G., Liu, Q. Y., Wang, X. Y., Zhou, Z., Zhao, X., Zhou, W., et al. (2023). A probiotic nanozyme hydrogel regulates vaginal microenvironment for Candida vaginitis therapy. Sci. Adv. 9, eadg0949. doi: 10.1126/sciadv.adg0949

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, M., Gao, J., Wu, Y., Li, Y., Chen, Y., Zhao, F., et al. (2020). Characterization of vaginal microbiota in Chinese women with cervical squamous intra-epithelial neoplasia. Int. J. Gynecol Cancer. 30, 1500–1504. doi: 10.1136/ijgc-2020-001341

PubMed Abstract | Crossref Full Text | Google Scholar

Xia, S., Li, S., Li, H. (2023). HPV-infection status and urinary incontinence: a population-based analysis of the NHANES 2005-2016. World J. Urol. 41, 1597–1603. doi: 10.1007/s00345-023-04425-9

PubMed Abstract | Crossref Full Text | Google Scholar

Xie, Y., Feng, Y., Li, W., Zhan, F., Huang, G., Hu, H., et al. (2020). Revealing the disturbed vaginal micobiota caused by cervical cancer using high-throughput sequencing technology. Front. Cell Infect. Microbiol. 10, 538336. doi: 10.3389/fcimb.2020.538336

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, J., Wang, W., Wang, Z., Wang, Z., Wang, Y., Wang, J., et al. (2020). Prevalence, genotype distribution and risk factors of cervical HPV infection in Yangqu, China: a population-based survey of 10086 women. Hum. Vaccin Immunother. 16, 1645–1652. doi: 10.1080/21645515.2019.1689743

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, Q., Wang, Y., Wei, X., Zhu, J., Wang, X., Xie, X., et al. (2020). The alterations of vaginal microbiome in HPV16 infection as identified by shotgun metagenomic sequencing. Front. Cell Infect. Microbiol. 10, 286. doi: 10.3389/fcimb.2020.00286

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, G., Gail, M. H., Consonni, D., Carugno, M., Humphrys, M., Pesatori, A. C., et al. (2016). Characterizing human lung tissue microbiota and its relationship to epidemiological and clinical features. Genome Biol. 17, 163. doi: 10.1186/s13059-016-1021-1

PubMed Abstract | Crossref Full Text | Google Scholar