Skip to main content

Changes in vaginal microbiome after focused ultrasound treatment of high-risk human papillomavirus infection-related low-grade cervical lesions

Abstract

Background

In this study, the changes of vaginal microbiome after focused ultrasound (FU) treatment were evaluated to explore the possible mechanism of FU in the treatment of high-risk human papillomavirus (HR-HPV) infection.

Methods

This study was nested in the FU arm of a prospective cohort study. A total of 37 patients diagnosed with HR-HPV infection-related cervical low-grade squamous intraepithelial lesion (LSIL) who met the inclusion criteria were enrolled in this study from October 2020 to November 2021, and these patients were treated with FU. We used 16S ribosomal RNA (16S rRNA) gene amplicon sequencing to profile the vaginal microbiota composition of patients before and 3 months after FU treatment.

Results

After FU treatment, HR-HPV was cleared in 24 patients, with a clearance rate of 75.0% (24/32). Lactobacillus iners was the predominant species among all samples. No significant difference was found in alpha-diversity index before and 3 months after FU treatment (P > 0.05), but the rarefaction curves showed that the vaginal microbial diversity before FU treatment was higher than that after FU treatment. Linear discriminant analysis (LDA) effect size (LEfSe) showed that Bifidobacterium contributed the most to the difference between the two groups at the genus level, and the abundance after FU treatment was significantly higher than that before treatment (P = 0.000).

Conclusions

The decrease of vaginal microbial diversity may be related to the clearance of HR-HPV infection, and FU treatment contributed to the decrease of vaginal microbial diversity. Increased Bifidobacterium abundance in the vaginal microbiome may be associated with clearance of HR-HPV infection, and FU treatment may contribute to the increase in Bifidobacterium abundance.

Trial registration number: This study was registered in the Chinese Clinical Trial Registry on 23/11/2020 (ChiCTR2000040162).

Peer Review reports

Background

Human papillomavirus (HPV) remains one of the most important oncogenic DNA viruses and has been identified as the cause of approximately 5% of cancers worldwide [1]. To date, more than 200 HPV genotypes have been identified, of which about 15 high-risk genotypes are considered to be necessary causes of cervical cancer (CC) [2]. In addition, socio-environmental factors, epidemiological factors, and genetic factors related to the host are considered to be co-factors that act together with HPV to influence the risk of CC development [3,4,5]. Over the past decade, emerging evidence suggests that the vaginal microbiome also plays a role in HPV infection, persistence, and cervical carcinogenesis [6,7,8].

In the lower female reproductive tract (vagina and cervix), the predominance of Lactobacillus is associated with vaginal health, and the depletion of these microorganisms can lead to a number of adverse conditions, such as increased risk of acquiring sexually transmitted diseases (STDs), spontaneous miscarriage, preterm birth or pelvic inflammatory disease [9]. Bacterial vaginosis (BV), characterized by depletion of Lactobacillus and overgrowth of anaerobic bacteria, has been reported to be significantly associated with increased risk of HPV infection and decreased HPV clearance [10]. Meanwhile, several studies have shown that women with cervical intraepithelial neoplasia (CIN) and CC have a vaginal microbial profile of community state type (CST) IV, characterized by depletion of Lactobacillus spp. and a substantial increase in vaginal microbiome diversity compared to healthy women [8, 11, 12].

As a novel and noninvasive treatment technology, focused ultrasound (FU) has achieved good curative effect in the treatment of cervical lesions, and has attracted much attention in China. In economically underdeveloped regions like southwestern China, some women diagnosed with high-risk human papillomavirus (HR-HPV) infection-related cervical low-grade squamous intraepithelial lesion (LSIL) lack previous screening data, have limited access to follow-up, and therefore urgently require treatment rather than observation. Our previous studies have confirmed that FU can effectively eliminate HR-HPV infection and reverse LSIL lesions [13, 14]. However, the response of vaginal microbiota to FU treatment remains unknown. Therefore, the purpose of this study was to explore the changes in the vaginal microbiome following FU treatment, and to provide clues for the mechanism of FU in clearing HR-HPV infection from the perspective of vaginal microbiota.

Methods

Study population and sample collection

From October 2020 to November 2021, women who visited the gynecology clinic of the Affiliated Hospital of North Sichuan Medical College, Nanchong city, Southwest China, were enrolled in this study if they met the following inclusion criteria: (1) women aged 18–55 years with a sexual history, confirmed HR-HPV infection by HPV-DNA test, and cervical biopsy revealed histological LSIL; (2) non-pregnant, non-lactating and non-menstruating women; (3) no antiviral treatment and corticosteroids and sex hormones application within 3 months before sampling; (4) no systemic or vaginal topical antibiotic therapy or antifungal therapy within 1 month before sampling; (5) no vaginal douching and administration within 1 week before sampling; and (6) no sexual activity within 72 h before sampling. Women with acute inflammation of the reproductive tract, sexually transmitted infections (STIs), immunodeficiency diseases, severe cardiac, hepatic and renal dysfunction, coagulopathy, and those who cannot be followed up were excluded.

HPV testing was performed using a High-Risk Human Papillomavirus Genotyping Real Time PCR Kit (Shanghai ZJ Bio-Tech Co., Ltd. China), which tests for 15 HR-HPV genotypes (16, 18,31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, 68, and 82). Cytological results were reported according to the 2014 Bethesda System [15], and histological diagnosis was obtained from two experienced pathologists using a 2-Tier Terminology as recommended by the LAST guidelines [16].

Samples for vaginal microbial analysis were collected from the upper third of the vaginal lateral wall using a sterile swab under direct visualization prior to FU treatment. The vaginal swab samples were immediately frozen and stored at − 80 °C until extraction. The sampling process at the 3-month follow-up after FU treatment was the same as above. The samples before and after FU treatment were collected during the follicular phase of the patients.

Focused ultrasound treatment

Women were treated with the Ultrasound Therapeutic Device (Model-CZF, Chongqing Haifu Medical Technology Co., Ltd. China). Women were placed in the lithotomy position and a sterile speculum was used to expose the cervix. The treatment probe was placed in close contact with the cervix, and a circular scanning was performed at the speed of 5–10 mm/s with the external cervical os as the center. The treatment typically lasted for 3–5 min until the lesion and the squamocolumnar junction shrunk with moderate inward depression of the external cervical os.

DNA extraction and MiSeq sequencing

Microbial DNA was extracted using the FastDNA® Spin Kit for Soil (MP Biomedicals, USA) according to the manufacturer’s instructions. All DNA samples were quality checked and the concentration was quantified by NanoDrop 2000 (Thermo Fisher Scientific, USA). Bacterial 16S rRNA gene fragments (V3–V4) were amplified from the extracted DNA using primers 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). Amplification was performed in 20-μL reactions with TransStart Fast Pfu DNA Polymerase (TransGen Biotech, China), 5 μM of each primer and 10 ng of template. The reactions were performed under the following thermal profile: 95 °C for 3 min, followed by 25 cycles of 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 45 s, and one cycle of 72 °C for 10 min and a 4 °C hold. PCR products were examined by 2% agarose gel electrophoresis and then purified using an AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, USA). Amplicons were quantified using the Quantus™ Fluorometer (Promega, USA). Amplicons were subjected to paired-end sequencing on the Illumina MiSeq sequencing platform at Majorbio Bio-Pharm Technology Company (Shanghai, China).

Sequence analysis

After demultiplexing, the resulting sequences were merged with FLASH (v1.2.11) [17] and quality filtered with fastp (0.19.6) [18]. Then the high-quality sequences were de-noised using DADA2 [19] plugin in the QIIME 2 [20] (version 2020.2) pipeline with recommended parameters, which obtained single nucleotide resolution based on error profiles within samples. DADA2 denoised sequences were usually called amplicon sequence variants (ASVs). Taxonomic assignment of ASVs was performed using the Naive baye consensus taxonomy classifier implemented in QIIME 2 and the Silva 16S rRNA database ((Release 138).

Alpha-diversity indices were calculated using mothur [21]. The richness and diversity of each sample were evaluated with the Sobs index and Shannon index, respectively. Wilcoxon rank-sum test was used to compare the alpha-diversity index between the pre-treatment group and the post-treatment group, and the Sob index and Shannon index were used to construct the rarefaction curves for comparison between the two groups. Hierarchical clustering analysis was performed on all samples based on Bray–Curtis distances and unweighted pair-group method with arithmetic means (UPGMA) algorithm. Species differences between the two groups were analyzed at the genus level using the Wilcoxon rank-sum test based on the microbial community abundance data. Linear discriminant analysis (LDA) effect size (LEfSe) algorithm was used to detect the microbial markers that contributed the most to the difference between the two groups. Q-value (P-value adjusted by false discovery rate (FDR)) and P-value < 0.05 were considered significant in this study.

Results

Clinical characteristics of the study population

A total of 37 patients were enrolled in this study. Thirty-seven pairs of vaginal microbial samples were obtained before and 3 months after FU treatment, with a total of 74 samples. After PCR amplification, the target bands of the PCR products were too weak or undetectable in 6 samples for subsequent experiments. Therefore, only 68 samples were finally sequenced, including 32 samples before treatment and 36 samples after treatment. The HR-HPV subtype with the highest infection rate in the study cohort was HPV16 (37.5%), followed by HPV52 (21.9%), and HPV51 (15.6%). At the 6-month follow-up, 24 patients were cleared of HR-HPV (defined as negative for all HR-HPV subtypes), with a clearance rate of 75.0% (24/32), while all 32 patients had cytology results of negative for intraepithelial lesion or malignancy (NILM). The clinical characteristics of patients are shown in Table 1.

Table 1 Clinical characteristics of patients (n = 32)

Sequencing results

A total of 11,957,834 optimized reads were obtained from 68 samples, with an average read lengths of 425 bp. For normalization, the reads in each sample were randomly subsampled to the lowest number of 108,824 in sample A19 (post-treatment group), and a total of 9784 ASVs were obtained.

Vaginal microbiota richness and diversity

At the ASV level, microbial richness and diversity were estimated using Sobs index and Shannon index, respectively. There was no significant difference in the Sobs index between the pre-treatment group and the post-treatment group (382.7 ± 1172.7 vs. 67.0 ± 50.6, P = 0.138), and no significant difference was found in the Shannon index between the two groups (1.2 ± 1.8 vs. 0.8 ± 0.7, P = 0.199) (Fig. 1).

Fig. 1
figure 1

Vaginal microbial richness and diversity in two groups. A Sobs index. B Shannon index. Wilcoxon rank-sum test was used to compare differences between the two groups

Rarefaction curves were constructed based on the Sobs index and Shannon index, respectively. The rarefaction curve based on the Shannon index reached a flat state, indicating that most of the microbial diversity had been generated, while the microbial diversity before treatment was higher than that after treatment. The rarefaction curve based on the Sobs index also indicated that the microbial richness before treatment was higher than that after treatment (Fig. 2).

Fig. 2
figure 2

Rarefaction curves based on Sobs index and Shannon index in two groups. A Rarefaction curve based on Sobs index. B Rarefaction curve based on Shannon index

Taxonomy of the vaginal microbiota in two groups

Overall, 35 bacterial phyla were recovered across all samples, and Firmicutes was the most abundant phylum, followed by Actinobacteria, Proteobacteria, Bacteroidetes, and Fusobacteria. Firmicutes accounted for 74.47% and 76.00% of the pre-treatment group and the post-treatment group, respectively, and there was no significant difference between the two groups (P = 0.969). Moreover, there was no significant difference between the two groups in the abundance of the top 10 bacteria at the phylum level (P > 0.05) (Fig. 3).

Fig. 3
figure 3

Comparison of the top 10 bacterial abundances between the two groups (phylum level)

At the genus level, a total of 802 taxa were found across all samples, with Lactobacillus being the most dominant genus overall, followed by Gardnerella, Streptococcus, Bifidobacterium, and Atopobium. Lactobacillus accounted for 66.22% and 60.20% of the pre-treatment group and the post-treatment group, respectively, and there was no significant difference between the two groups (P = 0.913). Among the top 10 bacteria in abundance, Bifidobacterium accounted for 0.14% and 5.14% of the pre-treatment group and the post-treatment group, respectively. The abundance of Bifidobacterium after treatment was significantly higher than that before treatment (P < 0.001) (Fig. 4).

Fig. 4
figure 4

Comparison of the top 10 bacterial abundances between the two groups (genus level). ***P < 0.001

Identification of vaginal microbial markers correlated with FU

We used the LEfSe model to identify microbial markers related to FU treatment (Fig. 5). The threshold for the logarithmic LDA model score for discriminative features in this study was 2.0 (P < 0.01). FU treatment was strongly associated with Bifidobacterium at the genus level (LDA = 4.40, P = 0.000).

Fig. 5
figure 5

The unique taxa and microbial markers for the two groups. Shown is a histogram of LDA scores computed for features differentially abundant in the two groups

Characteristics of vaginal microbial community types

At the species level, hierarchical clustering analysis revealed that all samples clustered into five major community types: dominated by Lactobacillus iners, dominated by Streptococcus anginosus, dominated by Gardnerella vaginalis, dominated by Lactobacillus crispatus, and dominated by Bifidobacterium dentiu (Fig. 6). However, the pre-treatment and post-treatment samples cannot be clearly separated (Fig. 7). The community type dominated by Lactobacillus iners accounted for the highest proportion in the pre-treatment group and post-treatment group (Fig. 8).

Fig. 6
figure 6

Hierarchical clustering tree of all samples at the species level. Based on Bray–Curtis distances and unweighted pair-group method with arithmetic means (UPGMA) algorithm

Fig. 7
figure 7

Typing analysis of all samples at the species level. Based on Jensen-Shannon distance

Fig. 8
figure 8

Typing histogram of the two groups at the species level

Discussion

In this study, we applied the 16S rRNA gene sequencing approach to a longitudinal cohort of patients before and after FU treatment to analyze the changes in the vaginal microbiome. A significant increase in Bifidobacterium abundance was observed in patients with HR-HPV infection-related LSIL at 3-month follow-up after FU treatment. Meanwhile, we found that the richness and diversity of vaginal microbiota decreased after FU treatment by rarefaction curves. Lactobacillus was still the most abundant genus in the two groups. However, at the species level, Lactobacillus iners was the predominant bacterium in the two groups.

The alpha-diversity index was widely used in microbial research. Several studies have suggested that HPV infection and progression of cervical lesions were associated with increased vaginal microbial diversity [22, 23]. Chen et al. [22] compared the Chao index and Shannon index in the HPV-negative group, HPV-positive group, LSIL group, high-grade squamous intraepithelial lesion (HSIL) group, and CC group. The results showed that the means of Chao and Shannon indices in the HPV-negative group was significantly lower than those in the HPV-positive group, while the values of Chao index and Shannon index in the CC group were the highest. However, other studies did not find a relationship between the alpha-diversity index and HPV infection and cervical lesions. Chalifour et al. [3] found no significant difference in the Shannon index between the HPV-negative group, HPV-positive group, and different stages of cervical lesions. In this study, we compared the Sobs index and Shannon index of the pre-treatment group and post-treatment group, and the results showed that the difference was not statistically significant (P > 0.05). However, from the rarefaction curves based on Sobs index and Shannon index, it could be found that the microbial richness and diversity before FU treatment were higher than those after FU treatment. Moreover, the HR-HPV clearance rate in the study cohort was 75.0% (24/32) at the 6-month follow-up. Pinder et al. [24] compared the efficacy of thermal ablation with cryotherapy or loop excision in a randomized controlled trial, showing that at the 6-month follow-up, the HR-HPV-negative rates were 40%, 42% and 47% in the cryotherapy, thermal ablation and LEEP treatment groups, respectively. In terms of the efficacy of FU, our results showed that FU was not inferior to other modalities. The microbial sequencing results suggested that the decrease in vaginal microbial diversity may be associated with the clearance of HR-HPV infection, and FU therapy contributed to the decrease of vaginal microbial diversity in patients.

Ravel et al. [25] classified vaginal microbiota into 5 major types, termed community state types (CSTs), based on vaginal bacterial species composition and abundance in reproductive-aged women. Lactobacillus crispatus, Lactobacillus gasseri, Lactobacillus iners, and Lactobacillus jensenii were the dominant bacteria in CST I, II, III, and V, respectively, while CST IV had a lower proportion of Lactobacillus and a higher proportion of strict and facultative anaerobes. We used hierarchical clustering and microbiota typing analysis to cluster the microbial communities of all samples into 5 types, and found that the dominant microbial community before and after FU treatment was dominated by Lactobacillus iners. Our results were consisted with other studies. A meta-analysis showed that vaginal microbiota dominated by Lactobacillus iners were significantly positively associated with HPV infection, cervical dysplasia and CC compared with Lactobacillus crispatus [12]. Many studies have found that the vaginal microbiota dominated by Lactobacillus crispatus (CST I) was associated with vaginal health status, whereas the vaginal microbiota dominated by Lactobacillus iners (CST III) was more common in patients with HPV infection or cervical lesions [26,27,28]. Lactobacillus iners can persist under the drastically changing vaginal environment due to its ability to respond and regulate its genomic functions [29], which may be the reason why FU therapy failed to alter the state of the vaginal microbial community.

We also found that Bifidobacterium at the genus level contributed the most to the difference between the pre-treatment group and post-treatment group, and the abundance of Bifidobacterium after FU treatment was significantly higher than that before FU treatment. Bifidobacterium is among the dominant bacterial populations in the human gastrointestinal tract (GIT), and also colonizes the vagina and oral cavity [30]. Studies have confirmed that Bifidobacterium in the GIT plays an important role in human physiology and nutrition [31]. However, there were few studies on the characteristics of Bifidobacterium in the vagina. Schellenberg et al. isolated vaginal Bifidobacterium with a ‘‘double-strong’’ phenotype of acid- and hydrogen peroxide (H2O2)- production [32], and the acid-producing property was thought to be responsible for maintaining vaginal microbial homeostasis. Freitas et al. found that lactic acid was produced by all vaginal Bifidobacterium isolates, but only 32% of isolates produced hydrogen peroxide [33]. These findings suggested that Bifidobacterium may be the dominant member of certain vaginal microbial communities, and may have protective or health-promoting effects in the vagina analogous to Lactobacillus by producing lactic acid [34].

Chao et al. [35] investigated the potential vaginal microbiome biomarkers for HSIL, and found that a paucity of Bifidobacterium, Dialister, unidentified Prevotellaceae, Faecalibacterium, and Bacteroides was related with HSIL. The relative abundance of Bifidobacterium being under 0.01% maybe a good predictor of HSIL in HPV16 and/or 18 infected individuals. Cha et al. [36] analyzed the antiviral activity of Bifidobacterium adolescentis SPM1005-A in the SiHa cervical cancer cell line expressing HPV 16, and found that the Bifidobacterium strain had antiviral activity through suppression E6 and E7 oncogene expression. Abdolalipour et al. [37] found that the intravenous or oral administration of Bifidobacterium bifidum can effectively induce antitumor immune responses and inhibit tumor growth in C57BL/6 mice bearing transplanted TC-1 cell of HPV-related tumor, expressing HPV 16 E6/E7 oncogenes. Our findings suggested that the increased abundance of Bifidobacterium in the vaginal microbiota may be associated with the clearance of HR-HPV infection, and FU treatment may contribute to the increased abundance of Bifidobacterium. However, studies on the role of Bifidobacterium in vaginal microbiota, HPV infection and HPV-related cervical lesions were still very insufficient.

The strength of this study is that it preliminarily explored the changes in the vaginal microbiome after FU treatment in patients with HR-HPV infection-related LSIL. We found that the diversity of the vaginal microbiome was decreased and the abundance of Bifidobacterium was increased after FU treatment in a prospective cohort, which provided clues to the mechanism of FU in clearing HR-HPV infection. We also confirmed that the vaginal microbial community of patients with HR-HPV infection-related LSIL was dominated by Lactobacillus iners. The limitations of this study were the small sample size and lack of a control group. Further studies with larger sample size and longer follow-up time are needed to assess the correlation between vaginal microbiome dynamics and FU treatment, HR-HPV infection clearance and LSIL lesions regression.

Conclusion

This work investigated the effect of FU treatment on the vaginal microbial composition. Bifidobacterium as a biomarker of vaginal microbial differences was associated with FU treatment, and increased Bifidobacterium abundance may be associated with HR-HPV infection clearance. This suggests that the characteristics and roles of Bifidobacterium in the vaginal microbial community deserve further study.

Availability of data and materials

All reads obtained were submitted to NCBI Sequence Read Archive (SRA) under the accession number PRJNA826816 (https://www.ncbi.nlm.nih.gov/sra/PRJNA826816).

Abbreviations

FU:

Focused ultrasound

HR-HPV:

High-risk human papillomavirus

HPV:

Human papillomavirus

LSIL:

Low-grade squamous intraepithelial lesion

CC:

Cervical cancer

STI:

Sexually transmitted infections

LEfSe:

Linear discriminant analysis (LDA) effect size

ASV:

Amplicon sequence variant

HSIL:

High-grade squamous intraepithelial lesion

CST:

Community state type

GIT:

Gastrointestinal tract

TCT:

ThinPrep cytology test

References

  1. Plummer M, de Martel C, Vignat J, Ferlay J, Bray F, Franceschi S. Global burden of cancers attributable to infections in 2012: a synthetic analysis. Lancet Glob Health. 2016;4(9):e609–16.

    Article  Google Scholar 

  2. Bhatla N, Aoki D, Sharma DN, Sankaranarayanan R. Cancer of the cervix uteri. Int J Gynaecol Obstet. 2018;143(Suppl 2):22–36.

    Article  Google Scholar 

  3. Audirac-Chalifour A, Torres-Poveda K, Bahena-Roman M, Tellez-Sosa J, Martinez-Barnetche J, Cortina-Ceballos B, Lopez-Estrada G, Delgado-Romero K, Burguete-Garcia AI, Cantu D, et al. Cervical microbiome and cytokine profile at various stages of cervical cancer: a pilot study. PLoS ONE. 2016;11(4): e0153274.

    Article  Google Scholar 

  4. Almonte M, Albero G, Molano M, Carcamo C, Garcia PJ, Perez G. Risk factors for human papillomavirus exposure and co-factors for cervical cancer in Latin America and the Caribbean. Vaccine. 2008;26(Suppl 11):L16-36.

    Article  Google Scholar 

  5. Zhang X, Zhang L, Tian C, Yang L, Wang Z. Genetic variants and risk of cervical cancer: epidemiological evidence, meta-analysis and research review. BJOG. 2014;121(6):664–74.

    Article  CAS  Google Scholar 

  6. Laniewski P, Ilhan ZE, Herbst-Kralovetz MM. The microbiome and gynaecological cancer development, prevention and therapy. Nat Rev Urol. 2020;17(4):232–50.

    Article  CAS  Google Scholar 

  7. Shannon B, Yi TJ, Perusini S, Gajer P, Ma B, Humphrys MS, Thomas-Pavanel J, Chieza L, Janakiram P, Saunders M, et al. Association of HPV infection and clearance with cervicovaginal immunology and the vaginal microbiota. Mucosal Immunol. 2017;10(5):1310–9.

    Article  CAS  Google Scholar 

  8. Laniewski P, Barnes D, Goulder A, Cui H, Roe DJ, Chase DM, Herbst-Kralovetz MM. Linking cervicovaginal immune signatures, HPV and microbiota composition in cervical carcinogenesis in non-Hispanic and Hispanic women. Sci Rep. 2018;8(1):7593.

    Article  Google Scholar 

  9. Martin DH, Marrazzo JM. The vaginal microbiome: current understanding and future directions. J Infect Dis. 2016;214(Suppl 1):S36-41.

    Article  Google Scholar 

  10. Kero K, Rautava J, Syrjanen K, Grenman S, Syrjanen S. Association of asymptomatic bacterial vaginosis with persistence of female genital human papillomavirus infection. Eur J Clin Microbiol Infect Dis. 2017;36(11):2215–9.

    Article  CAS  Google Scholar 

  11. Di Paola M, Sani C, Clemente AM, Iossa A, Perissi E, Castronovo G, Tanturli M, Rivero D, Cozzolino F, Cavalieri D, et al. Characterization of cervico-vaginal microbiota in women developing persistent high-risk Human Papillomavirus infection. Sci Rep. 2017;7(1):10200.

    Article  Google Scholar 

  12. Norenhag J, Olovsson M, Verstraelen H, Engstrand L, Brusselaersb CN. The vaginal microbiota, human papillomavirus and cervical dysplasia: a systematic review and network meta-analysis. BIOG. 2020;127(2):171–80.

    CAS  Google Scholar 

  13. Wang W, Liu Y, Pu Y, Li C, Zhou H, Wang Z. Effectiveness of focused ultrasound for high risk human papillomavirus infection-related cervical lesions. Int J Hyperthermia. 2021;38(2):96–102.

    Article  CAS  Google Scholar 

  14. Fu Z, Fan Y, Wu C, Yan P, Ye Y, Yang H, Li C. Clinical efficacy and mechanism for focused ultrasound (FUS) in the management of cervical intraepithelial neoplasia 1 (CIN1). Int J Hyperthermia. 2020;37(1):339–45.

    Article  CAS  Google Scholar 

  15. Nayar R, Wilbur DC. The Pap test and Bethesda 2014. Cancer Cytopathol. 2015;123(5):271–81.

    Article  Google Scholar 

  16. Darragh TM, Colgan TJ, Fau-Cox JT, Cox JT, Fau-Heller DS, Heller DS, Fau-Henry MR, Henry MR, Fau-Luff RD, Luff RD, Fau-McCalmont T, McCalmont T, Fau-Nayar R, Nayar R, Fau-Palefsky JM, Palefsky JM, Fau-Stoler MH, Stoler MH, Fau-Wilkinson EJ, et al. The lower anogenital squamous terminology standardization project for HPV-associated lesions: background and consensus recommendations from the College of American Pathologists and the American Society for Colposcopy and Cervical Pathology. Arch Pathol Lab Med. 2012;136(10):1266–97.

    Article  Google Scholar 

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

    Article  CAS  Google Scholar 

  18. Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90.

    Article  Google Scholar 

  19. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3.

    Article  CAS  Google Scholar 

  20. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, Alexander H, Alm EJ, Arumugam M, Asnicar F, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.

    Article  CAS  Google Scholar 

  21. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.

    Article  CAS  Google Scholar 

  22. Chen Y, Qiu X, Wang W, Li D, Wu A, Hong Z, Di W, Qiu L. Human papillomavirus infection and cervical intraepithelial neoplasia progression are associated with increased vaginal microbiome diversity in a Chinese cohort. BMC Infect Dis. 2020;20(1):629.

    Article  CAS  Google Scholar 

  23. Mitra A, MacIntyre DA, Lee YS, Smith A, Marchesi JR, Lehne B, Bhatia R, Lyons D, Paraskevaidis E, Li JV, et al. Cervical intraepithelial neoplasia disease progression is associated with increased vaginal microbiome diversity. Sci Rep. 2015;5:16865.

    Article  CAS  Google Scholar 

  24. Pinder LF, Parham GP, Basu P, Muwonge R, Lucas E, Nyambe N, Sauvaget C, Mwanahamuntu MH, Sankaranarayanan R, Prendiville W. Thermal ablation versus cryotherapy or loop excision to treat women positive for cervical precancer on visual inspection with acetic acid test: pilot phase of a randomised controlled trial. Lancet Oncol. 2020;21(1):175–84.

    Article  CAS  Google Scholar 

  25. Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SS, McCulle SL, Karlebach S, Gorle R, Russell J, Tacket CO, et al. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci USA. 2011;108(Suppl 1):4680–7.

    Article  CAS  Google Scholar 

  26. Brotman RM, Shardell MD, Gajer P, Tracy JK, Zenilman JM, Ravel J, Gravitt PE. Interplay between the temporal dynamics of the vaginal microbiota and human papillomavirus detection. J Infect Dis. 2014;210(11):1723–33.

    Article  CAS  Google Scholar 

  27. Oh HY, Kim BS, Seo SS, Kong JS, Lee JK, Park SY, Hong KM, Kim HK, Kim MK. The association of uterine cervical microbiota with an increased risk for cervical intraepithelial neoplasia in Korea. Clin Microbiol Infect. 2015;21(7):e1-9.

    Article  Google Scholar 

  28. Zheng N, Guo R, Wang J, Zhou W, Ling Z. Contribution of lactobacillus iners to vaginal health and diseases: a systematic review. Front Cell Infect Microbiol. 2021;11(792787):1–12.

    Google Scholar 

  29. Macklaim JM, Fernandes AD, Bella JMD, Hammond J-A, Reid G, Gloor GB. Comparative meta-RNA-seq of the vaginal microbiota and differential expression by Lactobacillus iners in health and dysbiosis. Microbiome. 2013;1(12):1–11.

    Google Scholar 

  30. Hidalgo-Cantabrana C, Delgado S, Ruiz L, Ruas-Madiedo P, Sanchez B, Margolles A. Bifidobacteria and their health-promoting effects. Microbiol Spectr. 2017;5(3):1–19.

    Article  CAS  Google Scholar 

  31. Turroni F, Milani C, Duranti S, Ferrario C, Lugli GA, Mancabelli L, van Sinderen D, Ventura M. Bifidobacteria and the infant gut: an example of co-evolution and natural selection. Cell Mol Life Sci. 2018;75(1):103–18.

    Article  CAS  Google Scholar 

  32. Schellenberg JJ, Dumonceaux TJ, Hill JE, Kimani J, Jaoko W, Wachihi C, Mungai JN, Lane M, Fowke KR, Ball TB, et al. Selection, phenotyping and identification of acid and hydrogen peroxide producing bacteria from vaginal samples of Canadian and East African Women. PLoS ONE. 2012;7(7):1–13.

    Article  Google Scholar 

  33. Freitas AC, Hill JE. Quantification, isolation and characterization of Bifidobacterium from the vaginal microbiomes of reproductive aged women. Anaerobe. 2017;10(47):145–56.

    Article  Google Scholar 

  34. Chaban B, Links MG, Jayaprakash TP, Wagner EC, Bourque DK, ZoeLohn AYA, Schalkwyk J, Reid G, Hemmingsen SM, Hill JE, Money DM. Characterization of the vaginal microbiota of healthy Canadian women through the menstrual cycle. Microbiome. 2014;2(23):1–12.

    Google Scholar 

  35. Chao X, Wang L, Wang S, Lang J, Tan X, Fan Q, Shi H. Research of the potential vaginal microbiome biomarkers for high-grade squamous intraepithelial lesion. Front Med (Lausanne). 2021;8: 565001.

    Article  Google Scholar 

  36. Cha M-K, Lee D-K, An H-M, Lee S-W, Shin S-H, Kwon J-H, Kim K-J, Ha N-J. Antiviral activity of Bifidobacterium adolescentis SPM1005-A on human papillomavirus type 16. BMC Med. 2012;10(72):1–6.

    Google Scholar 

  37. Abdolalipour E, Mahooti M, Salehzadeh A, Torabi A, Mohebbi SR, Gorji A, Ghaemi A. Evaluation of the antitumor immune responses of probiotic Bifidobacterium bifidum in human papillomavirus-induced tumor model. Microb Pathog. 2020;145: 104207.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the University-level Doctoral Foundation of North Sichuan Medical College (CBY22-QDA21).

Author information

Authors and Affiliations

Authors

Contributions

WW and HZ contributed to the conception and design of the study. WW, YL, and YY organized the study cohort. JR performed the FU treatment. WW performed the data analysis and wrote the first draft of manuscript. HZ reviewed and edited the manuscript. All authors contributed to manuscript revision and read and approved the submitted version. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Honggui Zhou.

Ethics declarations

Ethics approval and consent to participate

This prospective cohort study was approved by the ethics committees of the Affiliated Hospital of North Sichuan Medical College (Reference: 2020ER124-1), and was registered in the Chinese Clinical Trial Registry on 23/11/2020 (ChiCTR2000040162). All experiments were performed in accordance with relevant guidelines and regulations and all participants signed a written informed consent form.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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, W., Liu, Y., Yang, Y. et al. Changes in vaginal microbiome after focused ultrasound treatment of high-risk human papillomavirus infection-related low-grade cervical lesions. BMC Infect Dis 23, 3 (2023). https://doi.org/10.1186/s12879-022-07937-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-022-07937-8

Keywords