<figure><img src="image_1.png" /></figure>

# Analyzing Cervical Microbiome Composition in HIV-Infected Women with Different HPV Infection Profiles: A Pilot Study in Thailand

Kanya Preechasuth <fnref n="1" />,* ID, Lionel Brazier ²ID, Woottichai Khamduang <fnref n="1" />,³ID, Sayamon Hongjaisee <fnref n="3" />,<fnref n="4" /> ID,
Nantawan Wangsaeng <fnref n="5" /> and Nicole Ngo-Giang-Huong <fnref n="2" />,<fnref n="3" /> ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Preechasuth, K.; Brazier, L.; Khamduang, W.; Hongjaisee, S.; Wangsaeng, N.; Ngo-Giang-Huong, N. Analyzing Cervical Microbiome Composition in HIV-Infected Women with Different HPV Infection Profiles: A Pilot Study in Thailand.
*Microorganisms* **2024**, *12*, 1298.
https://doi.org/10.3390/microorganisms12071298

Academic Editors: Moriya Tsuji
and Francesco Broccolo

Received: 6 May 2024

Revised: 14 June 2024

Accepted: 20 June 2024

Published: 26 June 2024

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** We conducted a pilot study to analyze the microbiome in cervical samples of women living with HIV with various profiles of HPV infections. The participants had an average age of 41.5 years. Sequence analysis of 16S rRNA V3 gene amplicons was performed using next-generation sequencing technology (Ion Torrent PGM™). The bioinformatics pipeline was analyzed using the Find, Rapidly, OTUs with Galaxy Solution system (FROGS). Common genera were determined to identify Community State Types (CSTs). The cervical microbiome profiles showed a dominance of lactobacilli in 56% (five out of nine) of samples. All three women with normal cervical cells and high-risk HPV infection were classified as CST IV, characterized by anaerobic bacteria associated with bacterial vaginitis, such as *Gardnerella*, *Prevotella*, *Atopobium*, and *Sneathia*. Among the two women with abnormal cervical cells and high-risk HPV infection, one was classified as CST III, and the other had an unclassified profile dominated by *L. helveticus*. Four women with normal cervical cells and no HPV infection exhibited various CSTs. Our study demonstrated the feasibility of the protocol in analyzing the cervical microbiome. However, further analysis with a larger number of longitudinal samples is necessary to determine the role of cervical microbiota in HPV persistence, clearance, or the development of precancerous lesions.

**Keywords:** vaginal microbiome; human papillomavirus; NGS; community state types

## 1. Introduction

The cervicovaginal microbiome composition is thought to play a major role in the genital tract health of women [1,2]. The cervicovaginal microbiome comprises numerous bacterial species, of which lactobacilli are the most abundant species. Next-generation sequencing (NGS) approaches have been commonly used since they allow for the in-depth analysis of bacterial 16S ribosomal RNA (rRNA) and provide a detailed microbial community structure. Based on the relative abundance of rRNA, cervical bacterial communities have been categorized into five groups named Community State Types (CSTs). The dominant bacteria species in CST I, II, III, and V are *Lactobacillus crispatus*, *L. gasseri*, *L. iners*, and *L. jensenii*, respectively, while CST IV has a larger diversity and includes anaerobic bacteria such as members of genus *Atopobium*, *Prevotella*, *Sneathia*, or *Gardnerella* [3].

<footnote n="1">¹ Division of Clinical Microbiology, Faculty of Associated Medical Sciences, Chiang Mai University, Chiang Mai 50200, Thailand; woottichai.k@cmu.ac.th</footnote>

<footnote n="2">² Maladies Infectieuses et Vecteurs: Écologie, Génétique, Évolution et Contrôle (MIVEGEC), Centre National de la Recherche Scientifique (CNRS UMR5290), Institut de Recherche Pour le Développement (IRD224), Université of Montpellier, 34394 Montpellier CEDEX 5, France; lionel.brazier@ird.fr (L.B.); nicole.ngo-giang-huong2@ird.fr (N.N.-G.-H.)</footnote>

<footnote n="3">³ LUCENT International Collaboration, Faculty of Associated Medical Sciences, Chiang Mai University, Chiang Mai 50200, Thailand; sayamon.ho@cmu.ac.th</footnote>

<footnote n="4">⁴ Research Institute for Health Sciences, Chiang Mai University, Chiang Mai 50200, Thailand</footnote>

<footnote n="5">⁵ AMS-PHPT Research Collaboration, Faculty of Associated Medical Sciences, Chiang Mai University, Chiang Mai 50200, Thailand; nantawan.wangsaeng@phpt.org</footnote>

<footnote>* Correspondence: kanya.p@cmu.ac.th</footnote>
Cervical cancer is the fourth leading cause of cancer-related deaths in women world-
wide, with most of the cancers caused by HPV [4]. As a result of the advances in sequencing
technologies, numerous research studies have assessed the relationships between the cervi-
cal microbiome and infection/persistence/clearance of human papillomaviruses (HPV)
in women's cervix. A growing body of evidence suggests that the increased diversity of
vaginal microbiota combined with a reduction in the relative abundance of lactobacilli may
be involved in HPV acquisition, severity of cervical intra-epithelia neoplasia (CIN) and the
development of cervical cancer [1,2]. However, the cervical microbiome composition can
vary with ethnicity [3,5] and host immunity [6], which in turn may impact HPV infection
and persistence. Thus, analysis of the cervical microbiome should consider these factors.

Women living with HIV (WLWH) are at higher risk of having high-risk (HR)-HPV infection and an abnormal Pap smear compared to HIV-uninfected women [7–9]. Furthermore, it has been reported that the microbiomes of women living with HIV and co-infected with HPV have a higher frequency to be categorized into CST III and CST IV [10]. WLWH had higher bacterial richness than HIV-negative women. *Mycoplasma* *tales*, *Pseudomonadales*, and *Staphylococcus* were the most abundant operational taxonomic units (OTUs) associated with high-grade squamous intraepithelial lesions (HSIL) and it suggested that inflammation due to chronic cervical mycoplasma infection may contribute to HPV-dependent dysplasia [11].

Data on the cervical microbiome of 25 women in Thailand indicate that two micro-
biome communities are frequently found predominant of lactobacilli, with *L. iners* being
the most common species and non-lactobacilli-dominated with *Gardnerella* and *Atopobium*
the most frequently found [12,13]. However, there is still a lack of extensive profiles of the
cervical microbiome composition among HIV-positive women with different cervical cancer
screening results (including Pap smear outcomes and HPV genotypes) in Thailand. Further
studies are required to understand the role of cervical microbiota in HPV persistence,
clearance, and the development of precancerous lesions.

In this pilot study, we explored the microbiome communities in WLWH using NGS
technology with varying cervical cancer screening results, including Pap smear and HPV
genotypes. We demonstrated the feasibility of using NGS technology to longitudinally in-
vestigate the role of cervical microbiota in HPV persistence, clearance, and the development
of precancerous lesions.

**2. Materials and Methods**

2.1. *Subjects and Sample Collection*

Cervical samples used in this study were collected from ten women enrolled in the Papillo V study [14] and kept stored at -70 °C until use. The samples were randomly selected based on the results of Pap smear and HPV testing. In the “Normal” group, cervical swabs were collected from women with normal Pap smear (Normal) and no presence of HPV (Negative); Normal/Negative (NN). A second group included samples from women with high-grade squamous intraepithelial lesions (HSIL) on Pap smear and infection with high-risk (HR) HPV; HSIL/HR-HPV (HH). The third group included samples from women with normal Pap smear and infection with HR-HPV; Normal/HR-HPV (NH).

## 2.2. DNA Extraction

Bacterial genomic DNA was extracted using the NucleoSpin®Tissue (Macherey-Nagel, Düren, Germany) according to the manufacturer's instructions. The quality and quantity of DNA was measured by the Nanodrop as preliminary data. The DNA was stored at -70 °C until further processing.
## 2.3. DNA Libraries Preparation

A fragment of 181 bp of the V3 region of the 16S rRNA gene was amplified using primers Probio_Uni and Probio_Rev that included a 10-mer barcode and adaptor sequences needed for the Ion Torrent PGMTM sequencing technology. The libraries were prepared using the fusion method as previously reported with slight modifications [15,16]. Briefly, each DNA extract was amplified in four PCR reactions to avoid bias. For each PCR reaction, one microliter of DNA extract was added to 4 µL of 5× High-Fidelity PCR master mix, (Thermo Fisher, Paris, France), 0.4 µL of 10 µM Probio_Uni primer, 0.4 µL of 10 µM Probio_Rev primer, 0.4 µL of 5 mM dNTP, and 0.4 units Q5 DNA polymerase (New England Biolabs, Ipswich, MA, USA) in a final volume of 20 µL. PCR conditions were as follows: denaturation at 98 °C for 30 s, followed by 30 cycles at 98 °C for 30 s, at 50 °C for 30 s, at 72 °C 15 s, and final extension at 72 °C for 10 min. The four PCR products were pooled and purified twice using the AGENCOURT AMPure XP reagent (Beckman Coulter, Inc., Atlanta, GA, USA), then quantified with the Agilent 2100 bioanalyzer and high-sensitivity DNA kit (Aligent Technologies, Wald Bronn, Germany).

## 2.4. 16S RNA Gene Sequencing

The purified DNA libraries were pooled in equal concentrations (100 pM) then underwent emulsion PCR with the Ion PGM™ Hi-Q™ OT2 Kit in the Ion OneTouch™ 2 System (Life Technologies, Paris, France) according to the manufacturer's instructions. The bar-coded amplicon libraries were sequenced using the Ion PGM™ Hi-Q™ sequencing kit after loading on an Ion 318™ chips (Life Technologies, Paris, France). The raw sequences were pre-processed by deleting the adapter, then data files were converted to the Fastq format by the Ion Torrent Suite™ Software version 5.12.0 (Life Technologies, Paris, France). Acceptable sequence reads from Ion Torrent platform include chip loading efficiency above 60–70%, the number of reads meeting the expected output for the specific chip type are 5 million read for Ion 318 chips and reads with quality scores ≥ Q20.

## 2.5. Sequence Analysis

The Fastq data files were analyzed using the bioinformatics pipeline of Fast, Rapidly, OTUs with Galaxy Solution (FROGS) 3.2 [17]. Shortly, the raw sequences from Fastq files were pre-processed (FROGS pre-process) for deletion of sequences with unexpected lengths or with ambiguous bases, merging of reads, and dereplication. Sequences obtained were then clustered (FROGS clustering swarm) and PCR chimera in each sample were removed (FROGS remove chimera). Operational Taxonomic Unit (OTU) abundance was analyzed and OTUs with more than 10 were kept (FROGS OTU filters). Taxonomic identification of OTUs was performed using Ribosomal database project (RDP) classifier and the Silva138.1 16S database (FROGS Affiliation OTU). Finally, taxonomy distribution within each sample was obtained (FROGS affiliation stat). Genus-level identification was accepted at more than 97% sequencing identity. The person performing NGS and OTU classification was blinded to the women and samples' characteristics.

# 3. Results

## 3.1. Characteristics of Women and Samples

Cervical samples were collected from ten HIV-infected women on antiretroviral therapy [14]. Their characteristics and their Pap smear and HPV genotyping results are described in Table 1. Their mean age was 41.5 years, ranging from 31 to 61 years, and their mean CD4 count was 600 cells/mm³, ranging from 370 to 1253 cells/mm³.
<table><caption>**Table 1.** Samples' characteristics (N = 10).</caption><thead><tr><th>Results of Pap Smear/HPV Typing *</th><th>Sample Code</th><th>Ages</th><th>Level of CD4 (cells/mm³)</th><th>HPV Type</th></tr></thead><tbody><tr><td rowspan="3">HSIL/HR-HPV</td><td>HH 1</td><td>33</td><td>408</td><td>16</td></tr><tr><td>HH 2</td><td>42</td><td>637</td><td>52</td></tr><tr><td>HH 3</td><td>61</td><td>508</td><td>52</td></tr><tr><td rowspan="3">Normal/HR-HPV</td><td>NH 1</td><td>31</td><td>624</td><td>16, 33, 39, 58, 40</td></tr><tr><td>NH 2</td><td>40</td><td>537</td><td>52</td></tr><tr><td>NH 3</td><td>42</td><td>307</td><td>52</td></tr><tr><td rowspan="4">Normal/Negative</td><td>NN 1</td><td>31</td><td>464</td><td>Negative</td></tr><tr><td>NN 2</td><td>35</td><td>1253</td><td>Negative</td></tr><tr><td>NN 3</td><td>40</td><td>600</td><td>Negative</td></tr><tr><td>NN 4</td><td>60</td><td>601</td><td>Negative</td></tr></tbody></table>

## 3.2. Analysis of Cervical Microbiome

Amplification and sequencing were successfully performed on all samples. The summary of the run showed a loading density of 61%, usable reads at 71% with a total of 4,804,855 reads, and a median read length of 189 bp. Each sample contained an average of 448,695 reads (SD 97,128 reads) and sufficient OTU richness to represent species identity, shown in rarefaction curves reaching a plateau (Figure 1). The negative control had 1998 reads. One sample, HH 3 was excluded in the further analysis since the number obtained after OTU affiliation was 15,401, which could have occurred from contamination with exogeneous DNA.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Rarefaction curves calculated with OTUs richness of HH1 and HH2 (red), NH1-NH3 (green), and NN1-NN4 (blue).</figcaption></figure>

## 3.3. Composition of Bacterial and CST Typing

Raw FASTAQ file data were analyzed with FROGS [17] and showed that most common genus (*Lactobacillus*, *Gardnerella*, *Megashaera*, and *Atopobium*) were present in the vaginal microbiome. The relative abundance of the twelve most frequent genera in each sample is

<footnote>* HSIL: high-grade squamous intraepithelial lesions, HR-HPV: high-risk HPV, Normal: normal Pap smear result, Negative: no presence of HPV.</footnote>
shown in Figure 2. Three samples (HH1, NN2, NN3) belonged to CSTIII characterized by the predominance of *L. iners* (73–95%). Four samples (NH1, NH2, NH3, NN4) belonged to CST IV, characterized by diverse bacterial microbiota but predominantly of *G. vaginalis* (39–47%), Megasphaera, and Prevotella. One sample (NN1) was classified as CST III/IV as almost half was *L. iners* and it was diverse of anaerobic bacterial microbiota. One sample (HH2) including mostly *L. helviticus* was unclassified CST, and is shown in Table 2.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** The proportion of the twelve most common genera in each sample.</figcaption></figure>

<table><caption>**Table 2.** Predominant microbiome in cervical samples by taxonomic distribution.</caption><thead><tr><th rowspan="2">Bacterial Microbiota</th><th colspan="8">Percentage of Bacteria Microbiota in Samples</th></tr><tr><th>HH1</th><th>HH2</th><th>NH1</th><th>NH2</th><th>NH3</th><th>NN1</th><th>NN2</th><th>NN3</th><th>NN4</th></tr></thead><tbody><tr><td>Lactobacillus iners</td><td>73.46</td><td></td><td>10.16</td><td></td><td></td><td>46.13</td><td>95.39</td><td>73.42</td><td></td></tr><tr><td>L. helveticus</td><td>25.40</td><td>89.15</td><td></td><td></td><td></td><td></td><td>3.66</td><td></td><td></td></tr><tr><td>Gardnerella vaginalis</td><td></td><td></td><td>38.73</td><td>46.94</td><td>12.87</td><td>21.98</td><td></td><td>19.46</td><td>71.68</td></tr><tr><td>Megasphaera</td><td></td><td></td><td>15.75</td><td>26.56</td><td>11.79</td><td>20.99</td><td></td><td></td><td></td></tr><tr><td>Sneathia</td><td></td><td></td><td>4.49</td><td></td><td>21.02</td><td></td><td></td><td></td><td>19.29</td></tr><tr><td>Atopobium vaginae</td><td></td><td></td><td>12.49</td><td>8.37</td><td>9.30</td><td>1.72</td><td></td><td>4.25</td><td></td></tr><tr><td>Prevotella</td><td></td><td></td><td></td><td>16.02</td><td>28.03</td><td></td><td></td><td></td><td></td></tr><tr><td>Fastidiosipila</td><td></td><td></td><td>11.88</td><td></td><td>2.57</td><td></td><td></td><td></td><td></td></tr><tr><td>Ureaplasma</td><td></td><td>0.14</td><td></td><td></td><td></td><td></td><td></td><td></td><td>2.32</td></tr><tr><td>Community State Type (CST)</td><td>III</td><td>Unclassified</td><td>IV</td><td>IV</td><td>IV</td><td>III/IV</td><td>III</td><td>III</td><td>IV</td></tr></tbody></table>

## 4. Discussion

Various regions of 16S rRNA can be used as sequencing targets, such as V3, V3-V4, V4 and V6-V7 regions [18,19]. Our study focused on the use of the V3 region of 16S rRNA to analyze the cervical microbiome, employing the Probio_Uni and Probio_Rev primers (181 bp) optimized for the Ion Torrent PGM™ system [15]. Despite the small size of the amplicons, our results successfully identified key bacterial phyla (Firmicutes, Actinobacteria, Bacteroidetes and Fusobacteria), as well as common genera and species. This demonstrates that the V3 region is effective for identifying CSTs in cervical microbiome studies. Our finding aligns with the study by Sirichoat et al., which indicated that the V3 region (approximately 215 bp) yielded the highest diversity of bacterial species compared to the V6-V7, V4, and V9 regions [19].
Our study was exploratory with a few samples in each group, limiting the depth of our conclusions. We observed that four women had CST III, which is characterized by the predominance of *L. iners* that has been associated with a higher risk of bacterial and viral infections and the development of precancerous lesions in Caucasian women [1,20]. However, *L. iners* is commonly found in Asian women and may not indicate a similar risk of viral infection or precancerous lesions in that context. In our study, only one of three women with the predominance of *L. iners* had a high-risk HPV infection.

We also found *L. helveticus* in one sample from a woman with high-risk HPV infection and altered cervical cells (HH2). This species was reported as one of vaginal microbiota in Indian women who have bacterial vaginitis [20,21]. Additionally, *L. helveticus* has been suggested to improve the vaginal health index in perimenopausal women with bacterial vaginosis, indicating a positive role in maintaining a healthy vaginal microbiome [22]. *L. crispatus* and *L. gasseri* which are common in White or European women, were absent in our study, consistent with their lesser prevalence in Asian, Black and Hispanic women [3]. The absence of these species suggests a potentially reduced ability to clear HPV in the population of Chinese women [23].

This pilot study aimed to assess the microbiome community in cervical samples of women living with HIV on antiretroviral therapy (ART). A previous meta-analysis reported that ART was associated with a lower prevalence of high-risk HPV compared to those not on ART [24]. Additionally, ART did not appear to disrupt microbial communities, as evidenced by the study by Liu et al., which showed that the vaginal microbiome remained stable post-ART initiation in women from sub-Saharan Africa [25]. Therefore, the microbiome results in our study are likely not affected by ART.

**5. Conclusions**

The protocol employed in this pilot study demonstrated the feasibility of identifying CSTs in most of the samples. However, due to the limited number of samples evaluated, no clear conclusion could be drawn regarding the relationship between CSTs and the risk of HPV infection or lesions. Future studies with larger and longitudinal sample sets are necessary to understand the role of cervical microbiota in HPV infection and persistence. This knowledge could potentially lead to treatments or prevention strategies for cervical precancerous lesions through the replacement with beneficial cervical microbiota.

**Author Contributions:** Conceptualization, K.P., W.K. and N.N.-G.-H.; methodology, K.P., L.B., N.W. and N.N.-G.-H.; validation, K.P. and L.B.; formal analysis, K.P. and L.B.; investigation, K.P. and L.B.; resources, L.B.; data curation, S.H. and N.N.-G.-H.; writing—original draft preparation, K.P. and N.N.-G.-H.; writing—review and editing, K.P., L.B., W.K., S.H., N.W. and N.N.-G.-H.; visualization, K.P. and S.H.; funding acquisition, K.P., W.K. and N.N.-G.-H. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Faculty of Associated Medical Sciences, Chiang Mai University, Thailand, the IRD Campagne de projection au Sud et mobilités Sud-Nord/Sud-Sud 2019, and Chiang Mai University, Thailand.

**Data Availability Statement:** The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

**Conflicts of Interest:** The authors declare no conflicts of interest.

**References**

1. Mitra, A.; MacIntyre, D.A.; Marchesi, J.R.; Lee, Y.S.; Bennett, P.R.; Kyrgiou, M. The Vaginal Microbiota, Human Papillomavirus Infection and Cervical Intraepithelial Neoplasia: What Do We Know and Where Are We Going Next? *Microbiome* **2016**, *4*, 58. [CrossRef] [PubMed]

2. Mitra, A.; MacIntyre, D.A.; Lee, Y.S.; Smith, A.; Marchesi, J.R.; Lehne, B.; Bhatia, R.; Lyons, D.; Paraskevaidis, E.; Li, J.V.; et al. Cervical Intraepithelial Neoplasia Disease Progression Is Associated with Increased Vaginal Microbiome Diversity. *Sci. Rep.* **2015**, *5*, 16865. [CrossRef] [PubMed]
3. Ravel, J.; Gajer, P.; Abdo, Z.; Schneider, G.M.; Koenig, S.S.K.; McCulle, S.L.; Karlebach, S.; Gorle, R.; Russell, J.; Tacket, C.O.; et al. Vaginal Microbiome of Reproductive-Age Women. *Proc. Natl. Acad. Sci. USA* **2011**, *108*, 4680–4687. [CrossRef] [PubMed]

4. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. *CA Cancer J. Clin.* **2021**, *71*, 209–249. [CrossRef] [PubMed]

5. Borgdorff, H.; Van Der Veer, C.; Van Houdt, R.; Alberts, C.J.; De Vries, H.J.; Bruisten, S.M.; Snijder, M.B.; Prins, M.; Geerlings, S.E.; Schim Van Der Loeff, M.F.; et al. The Association between Ethnicity and Vaginal Microbiota Composition in Amsterdam, The Netherlands. *PLoS ONE* **2017**, *12*, e0181135. [CrossRef] [PubMed]

6. Murall, C.L.; Rahmoun, M.; Selinger, C.; Baldellou, M.; Bernat, C.; Bonneau, M.; Boué, V.; Buisson, M.; Christophe, G.; D'Auria, G.; et al. Natural History, Dynamics, and Ecology of Human Papillomaviruses in Genital Infections of Young Women: Protocol of the PAPCLEAR Cohort Study. *BMJ Open* **2019**, *9*, e025129. [CrossRef] [PubMed]

7. Adler, D.H.; Wallace, M.; Bennie, T.; Mrubata, M.; Abar, B.; Meiring, T.L.; Williamson, A.L.; Bekker, L.G. Cervical Dysplasia and High-Risk Human Papillomavirus Infections among HIV-Infected and HIV-Uninfected Adolescent Females in South Africa. *Infect. Dis. Obstet. Gynecol.* **2014**, *14*, 498048. [CrossRef]

8. Hanisch, R.A.; Sow, P.S.; Toure, M.; Dem, A.; Dembele, B.; Toure, P.; Winer, R.L.; Hughes, J.P.; Gottlieb, G.S.; Feng, Q.; et al. Influence of HIV-1 and/or HIV-2 Infection and CD4 Count on Cervical HPV DNA Detection in Women from Senegal, West Africa. *J. Clin. Virol.* **2013**, *58*, 696–702. [CrossRef]

9. McDonald, A.C.; Tergas, A.I.; Kuhn, L.; Denny, L.; Wright, T.C. Distribution of Human Papillomavirus Genotypes among HIV-Positive and HIV-Negative Women in Cape Town, South Africa. *Front. Oncol.* **2014**, *4*, 48. [CrossRef]

10. Siqueira, J.D.; Curty, G.; Xutao, D.; Hofer, C.B.; Machado, E.S.; Seuánez, H.N.; Soares, M.A.; Delwart, E.; Soares, E.A. Composite Analysis of the Virome and Bacteriome of HIV/HPV Co-Infected Women Reveals Proxies for Immunodeficiency. *Viruses* **2019**, *11*, 422. [CrossRef]

11. Klein, C.; Gonzalez, D.; Samwel, K.; Kahesa, C.; Mwaiselage, J.; Aluthge, N.; Fernando, S.; West, J.T.; Wood, C.; Angeletti, P.C. Relationship between the Cervical Microbiome, HIV Status, and Precancerous Lesions. *MBio* **2019**, *10*, e02785-18. [CrossRef] [PubMed]

12. Sirichoat, A.; Buppasiri, P.; Engchanil, C.; Namwat, W.; Faksri, K.; Sankuntaw, N.; Pasomsub, E.; Chantratita, W.; Lulitanond, V. Characterization of Vaginal Microbiota in Thai Women. *PeerJ* **2018**, *2018*, e5977. [CrossRef] [PubMed]

13. Sasivimolrattana, T.; Chantratita, W.; Sensorn, I.; Chaiwongkot, A.; Oranratanaphan, S.; Bhattarakosol, P. Human Virome in Cervix Controlled by the Domination of Human Papillomavirus. *Viruses* **2022**, *14*, 2066. [CrossRef] [PubMed]

14. Delory, T.; Ngo-Giang-Huong, N.; Rangdaeng, S.; Chotivanich, N.; Limtrakul, A.; Putiyanun, C.; Suriyachai, P.; Matanasarawut, W.; Jarupanich, T.; Liampongsabuddhi, P.; et al. Human Papillomavirus Infection and Cervical Lesions in HIV Infected Women on Antiretroviral Treatment in Thailand. *J. Infect.* **2017**, *74*, 501–511. [CrossRef] [PubMed]

15. Milani, C.; Hevia, A.; Foroni, E.; Duranti, S.; Turroni, F.; Lugli, G.A.; Sanchez, B.; Martín, R.; Gueimonde, M.; van Sinderen, D.; et al. Assessing the Fecal Microbiota: An Optimized Ion Torrent 16S rRNA Gene-Based Analysis Protocol. *PLoS ONE* **2013**, *8*, e68739. [CrossRef] [PubMed]

16. Brazier, L.; Koumavor, C.K.; Renaud, N.; Prugnolle, F.; Thomas, F.; Ategbo, S.; Engoba, M.; Obengui, Leroy, E.M.; Durand, P.; et al. Evolution in Fecal Bacterial/Viral Composition in Infants of Two Central African Countries (Gabon and Republic of the Congo) during Their First Month of Life. *PLoS ONE* **2017**, *12*, e185569. [CrossRef] [PubMed]

17. Escudié, F.; Auer, L.; Bernard, M.; Mariadassou, M.; Cauquil, L.; Vidal, K.; Maman, S.; Hernandez-Raquet, G.; Combes, S.; Pascal, G. FROGS: Find, Rapidly, OTUs with Galaxy Solution. *Bioinformatics* **2018**, *34*, 1287–1294. [CrossRef] [PubMed]

18. Onywera, H.; Meiring, T.L. Comparative Analyses of Ion Torrent V4 and Illumina V3-V4 16S rRNA Gene Metabarcoding Methods for Characterization of Cervical Microbiota: Taxonomic and Functional Profiling. *Sci. Afr.* **2020**, *7*, e00278. [CrossRef]

19. Sirichoat, A.; Sankuntaw, N.; Engchanil, C.; Buppasiri, P.; Faksri, K.; Namwat, W.; Chantratita, W.; Lulitanond, V. Comparison of Different Hypervariable Regions of 16S rRNA for Taxonomic Profiling of Vaginal Microbiota Using Next-Generation Sequencing. *Arch. Microbiol.* **2021**, *203*, 1159–1166. [CrossRef]

20. Ceccarani, C.; Foschi, C.; Parolin, C.; D'Antuono, A.; Gaspari, V.; Consolandi, C.; Laghi, L.; Camboni, T.; Vitali, B.; Severgnini, M.; et al. Diversity of Vaginal Microbiome and Metabolome during Genital Infections. *Sci. Rep.* **2019**, *9*, 14095. [CrossRef]

21. Pramanick, R.; Nathani, N.; Warke, H.; Mayadeo, N.; Aranha, C. Vaginal Dysbiotic Microbiome in Women With No Symptoms of Genital Infections. *Front. Cell. Infect. Microbiol.* **2022**, *11*, 1369. [CrossRef] [PubMed]

22. Vivekanandan, V.; Khan, Z.H.; Venugopal, G.; Musunuru, B.; Mishra, P.; Srivastava, S.; Ramadass, B.; Subhadra, B. VagiBIOM Lactobacillus Suppository Improves Vaginal Health Index in Perimenopausal Women with Bacterial Vaginosis: A Randomized Control Trial. *Sci. Rep.* **2024**, *14*, 3317. [CrossRef] [PubMed]

23. Zhang, Y.; Xu, X.; Yu, L.; Shi, X.; Min, M.; Xiong, L.; Pan, J.; Zhang, Y.; Liu, P.; Wu, G.; et al. Vaginal Microbiota Changes Caused by HPV Infection in Chinese Women. *Front. Cell. Infect. Microbiol.* **2022**, *12*, 814668. [CrossRef] [PubMed]
24. Kelly, H.; Weiss, H.A.; Benavente, Y.; de Sanjose, S.; Mayaud, P. Association of Antiretroviral Therapy with High-Risk Human Papillomavirus, Cervical Intraepithelial Neoplasia, and Invasive Cervical Cancer in Women Living with HIV: A Systematic Review and Meta-Analysis. *lancet*. *HIV 2018*, *5*, e45-e58. [CrossRef]

25. Liu, C.M.; Packman, Z.R.; Abraham, A.G.; Serwadda, D.M.; Nalugoda, F.; Aziz, M.; Prodger, J.L.; Kaul, R.; Kalibbala, S.; Gray, R.H.; et al. The Effect of Antiretroviral Therapy Initiation on the Vaginal Microbiome in HIV-Infected Women. *Open Forum Infect. Dis.* **2019**, *6*, ofz328. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.