<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="brief-report" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Med.</journal-id>
<journal-title>Frontiers in Medicine</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Med.</abbrev-journal-title>
<issn pub-type="epub">2296-858X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmed.2025.1615571</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Medicine</subject>
<subj-group>
<subject>Brief Research Report</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Cervical microbiota diversity and functional shifts in the development of cervical high-grade squamous intraepithelial lesions</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" equal-contrib="yes" corresp="yes"><name><surname>Rosas Cancio-Su&#x00E1;rez</surname> <given-names>Marta</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="aff" rid="aff3"><sup>3</sup></xref><xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref><xref ref-type="author-notes" rid="fn0003"><sup>&#x2020;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2407920/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
</contrib>
<contrib contrib-type="author" corresp="yes" equal-contrib="yes"><name><surname>Moreno</surname> <given-names>Elena</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref><xref ref-type="author-notes" rid="fn0003"><sup>&#x2020;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1623681/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
</contrib>
<contrib contrib-type="author"><name><surname>del Valle Rubido</surname> <given-names>Cristina</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
</contrib>
<contrib contrib-type="author"><name><surname>Salvador</surname> <given-names>Marta</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
</contrib>
<contrib contrib-type="author"><name><surname>Moreno</surname> <given-names>Ana</given-names></name><xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
</contrib>
<contrib contrib-type="author"><name><surname>Luna</surname> <given-names>Laura</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
</contrib>
<contrib contrib-type="author"><name><surname>D&#x00ED;az-Garc&#x00ED;a</surname> <given-names>Claudio</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/3206880/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author"><name><surname>Tapia</surname> <given-names>Carlos</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
</contrib>
<contrib contrib-type="author"><name><surname>del Amo</surname> <given-names>Ana</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/3186717/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author"><name><surname>Moreno</surname> <given-names>Santiago</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1978981/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author"><name><surname>S&#x00E1;nchez-Conde</surname> <given-names>Matilde</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1566589/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
</contrib>
<contrib contrib-type="author"><name><surname>Serrano-Villar</surname> <given-names>Sergio</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="aff" rid="aff6"><sup>6</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/679558/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Infectious Diseases, IRYCIS, and Universidad de Alcal&#x00E1;, Hospital Universitario Ram&#x00F3;n y Cajal</institution>, <addr-line>Madrid</addr-line>, <country>Spain</country></aff>
<aff id="aff2"><sup>2</sup><institution>CIBERINFEC, Instituto de Salud Carlos III</institution>, <addr-line>Madrid</addr-line>, <country>Spain</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Medicine, University of Alcal&#x00E1; de Henares, Guadalajara Campus</institution>, <addr-line>Alcal&#x00E1; de Henares</addr-line>, <country>Spain</country></aff>
<aff id="aff4"><sup>4</sup><institution>Department of Ginecology, IRYCIS, Hospital Universitario Ram&#x00F3;n y Cajal</institution>, <addr-line>Madrid</addr-line>, <country>Spain</country></aff>
<aff id="aff5"><sup>5</sup><institution>Alpes Primary Care Center</institution>, <addr-line>Madrid</addr-line>, <country>Spain</country></aff>
<aff id="aff6"><sup>6</sup><institution>Department of Medicine, University Antonio de Nebrija, Life Sciences Campus</institution>, <addr-line>Madrid</addr-line>, <country>Spain</country></aff>
<author-notes>
<fn fn-type="edited-by" id="fn0001">
<p>Edited by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/275529/overview">Ant&#x00F3;nio Machado</ext-link>, University of the Azores, Portugal</p>
</fn>
<fn fn-type="edited-by" id="fn0002">
<p>Reviewed by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1871315/overview">Thanayod Sasivimolrattana</ext-link>, Mahidol University, Thailand</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3060602/overview">Bhagyashri Patil-Takbhate</ext-link>, Hospital and Research Centre, India</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3081651/overview">Irma Berenice Mulato Briones</ext-link>, National Polytechnic Institute, Mexico</p>
</fn>
<corresp id="c001">&#x002A;Correspondence: Elena Moreno, <email>emolmo@salud.madrid.org</email>; Marta Rosas Cancio-Su&#x00E1;rez, <email>marta.rosas@salud.madrid.org</email></corresp>
<fn fn-type="equal" id="fn0003"><p><sup>&#x2020;</sup>These authors have contributed equally to this work</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>24</day>
<month>10</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="collection">
<year>2025</year>
</pub-date>
<volume>12</volume>
<elocation-id>1615571</elocation-id>
<history>
<date date-type="received">
<day>26</day>
<month>04</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>08</day>
<month>10</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2025 Rosas Cancio-Su&#x00E1;rez, Moreno, del Valle Rubido, Salvador, Moreno, Luna, D&#x00ED;az-Garc&#x00ED;a, Tapia, del Amo, Moreno, S&#x00E1;nchez-Conde and Serrano-Villar.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Rosas Cancio-Su&#x00E1;rez, Moreno, del Valle Rubido, Salvador, Moreno, Luna, D&#x00ED;az-Garc&#x00ED;a, Tapia, del Amo, Moreno, S&#x00E1;nchez-Conde and Serrano-Villar</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Research on microbial changes in the cervix, where most human papillomavirus (HPV) complications arise, is limited. Here, we aimed to understand the specific role of the cervicovaginal microbiota in developing high-grade squamous intraepithelial lesions (HSIL) associated with HPV infection. Our results show higher diversity in the microbiota associated with HSIL, with the genera <italic>Parvimonas, Fastidiosipila</italic>, and <italic>Pseudomonas</italic> being the most abundant. Additionally, an imputed functional analysis revealed that pathways such as glycine, serine, threonine, and sulfur metabolism were enriched in cervical samples from women with HSIL. Identifying biomarkers that help prevent HSIL progression could benefit women at risk of developing HPV-related cancerous lesions.</p>
</abstract>
<kwd-group>
<kwd>HSIL</kwd>
<kwd>cervicovaginal microbiota</kwd>
<kwd>16S</kwd>
<kwd>HPV</kwd>
<kwd>bacterial functions</kwd>
</kwd-group>
<counts>
<fig-count count="2"/>
<table-count count="1"/>
<equation-count count="0"/>
<ref-count count="24"/>
<page-count count="7"/>
<word-count count="4924"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Obstetrics and Gynecology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="sec1">
<title>Background</title>
<p>The cervicovaginal microbiota has co-evolved with mucosal immune cells to facilitate fetal implantation and development during pregnancy and protect the female urogenital tract from infection (<xref ref-type="bibr" rid="ref1">1</xref>, <xref ref-type="bibr" rid="ref2">2</xref>). The correlation between cervicovaginal microbiota and its role in either fostering susceptibility or providing protection against HPV, the third leading cause of cervical cancer, is gaining attention (<xref ref-type="bibr" rid="ref3">3</xref>).</p>
<p>Research on bacterial communities has revealed that the optimal vaginal microbiota is predominantly composed of <italic>Lactobacillus crispatus</italic>, <italic>L. gasseri</italic>, and <italic>L. jenseni</italic>. These species secrete lactobiocins and biosurfactants that prevent other species from binding to the epithelium (<xref ref-type="bibr" rid="ref4">4</xref>). Interestingly, a vaginal microbiota dominated by <italic>L. iners</italic> has been associated with higher odds of HPV prevalence, including high-risk HPV and cervical cancer or dysplasia (<xref ref-type="bibr" rid="ref5">5</xref>). It has been established that the interaction between HPV, host immunity, and the vaginal microenvironment influences the risk of persistent HPV infection and the progression of cervical intraepithelial neoplasia (CIN) (<xref ref-type="bibr" rid="ref6">6</xref>, <xref ref-type="bibr" rid="ref7">7</xref>). However, while HPV-related complications primarily develop in the cervix, the cervical microbiota has been less studied. Comparative analyses not focused on HSIL further show cervix&#x2013;vagina differences, underscoring site-specific composition and detectability, and reporting variable Lactobacillus dominance (<xref ref-type="bibr" rid="ref8">8</xref>, <xref ref-type="bibr" rid="ref9">9</xref>). Furthermore, decreased levels of Lactobacillus have been reported in women with low HPV viral loads (<xref ref-type="bibr" rid="ref10">10</xref>).</p>
<p>Therefore, we conducted a study to explore the connections between HPV infection and squamous intraepithelial lesions, focusing on the composition and function of the cervical microbiota. We aimed to identify metabolites that can help identify risk markers for malignant or premalignant lesions, enhance current screening efforts, and identify potential therapeutic targets in the long term.</p>
</sec>
<sec sec-type="methods" id="sec2">
<title>Methods</title>
<sec id="sec3">
<title>Ethics</title>
<p>This study conformed to the principles of the Declaration of Helsinki and Good Clinical Practice Guidelines and was approved by the Independent Ethics Committee of the coordinating center Hospital Ram&#x00F3;n y Cajal (<email>ceic.hrc@salud.madrid.org</email>, approval 146-21l) and the participating centers. All participants provided informed consent before the initiation of the study procedures.</p>
</sec>
<sec id="sec4">
<title>Study population</title>
<p>We conducted an observational, prospective, single-center study. From July 2022 to July 2023, we systematically invited women undergoing routine cervical cytology consultations at their Primary Care Center or Gynecology department assigned to the project in Madrid, Spain. Participants were excluded if they were under 18&#x202F;years of age, declined to sign informed consent, had received systemic antibiotics in the previous month, or presented with systemic inflammatory conditions (such as chronic hepatitis B or C, intercurrent infections, or autoimmune diseases), in order to minimize major disturbances to the microbiota. We included 105 women classified into three groups: women who tested negative for HPV and had normal cytology (NEG), those who tested positive for HPV and had low-grade squamous intraepithelial lesions (LSIL), and those who tested positive and had high-grade squamous intraepithelial lesions (HSIL).</p>
<p>All the women provided written consent to participate in the study. Endocervical samples were collected using sterile endocervical brushes (Cervex-Brush&#x00AE; Combi, Rovers Medical Devices) by trained gynecologists or primary care physicians. The brush was gently inserted into the endocervical canal and rotated 360&#x00B0; to maximize epithelial contact while minimizing contamination from the vaginal introitus.</p>
<p>HPV detection was performed using the Cobas&#x00AE; 4,800 HPV Test, which individually detects HPV16 and HPV18, and reports the presence of 12 additional high-risk genotypes (HPV-31, &#x2212;33, &#x2212;35, &#x2212;39, &#x2212;45, &#x2212;51, &#x2212;52, &#x2212;56, &#x2212;58, &#x2212;59, &#x2212;66, and &#x2212;68) collectively as &#x201C;other HR-HPV.&#x201D; No specific genotyping was performed for these additional genotypes.</p>
</sec>
<sec id="sec5">
<title>16S rRNA sequencing analysis</title>
<p>Cryopreserved samples obtained in sterile iSWAB-Microbiome tubes were thawed and DNA was extracted using the QIAamp DNA Mini Kit. For each sample, the V3-V4 regions of the 16S rRNA gene were analyzed using Illumina MiSeq (Novogene Bioinformatics Technology Co., Ltd.). 50&#x202F;K reads with 250&#x202F;bp Pair-Ends inserts were analyzed for OTUs cluster and phylogenetic relationship construction. Different Alpha and Beta diversity indexes were calculated using Simpson, Shannon, and Principal Coordinate Analysis (PCoA) respectively. Differential abundance and function prediction analyses were performed by analyzing the compositions of microbiomes with bias correction (ANCOM-BC2) (<xref ref-type="bibr" rid="ref11">11</xref>) identifying significantly different orthologous genes using the Kyoto Encyclopedia of Genes and Genomes (KEGG/KO) and Enzyme Commission (EC) terms between groups in the dataset. Bioinformatic analysis was performed using the Novogene pipeline. We then selected these terms and used MicrobiomeAnalyst to identify the enriched metabolic pathways or biological processes associated with them.</p>
</sec>
<sec id="sec6">
<title>Statistical analysis</title>
<p>Study data were collected and managed using REDCap electronic data capture tools v13.1.27. We report qualitative variables as frequency distribution and quantitative variables as medians with their 25th-75th percentiles (P25-P75) or as means and standard deviation (SD), according to the data distribution. We performed comparisons between groups using the Chi-square test for categorical variables. Since the distribution of all the assessed variables deviated from normality, we used the Mann&#x2013;Whitney U test for between-group comparisons of continuous variables. We used Stata v18.0 (StataCorp LP, College Station, TX, USA) to perform the statistical tests. Beta diversity was primarily assessed using weighted UniFrac distances to account for differences in taxon abundance across groups. Differential abundance of bacterial taxa was assessed using ANCOM-BC2, with statistical significance defined as a false discovery rate (FDR)&#x2013;adjusted q-value &#x003C; 0.05. For functional enrichment analysis based on KEGG Orthologs (KO) and Enzyme Commission (EC) terms, we considered pathways with a raw <italic>p</italic>-value &#x003C; 0.01 as significant.</p>
</sec>
</sec>
<sec sec-type="results" id="sec7">
<title>Results</title>
<sec id="sec8">
<title>General characteristics</title>
<p>We included 105 women: 87, and 18 in the non-HSIL, and HSIL groups, respectively. The median age of the participants was similar across the groups: 41&#x202F;years in the non-HSIL group, 38.5&#x202F;years in the HSIL group. Most participants were Spanish (86.2, and 72.2% of the non-HSIL, and HSIL groups, respectively). A smaller proportion of participants were from Eastern Europe and South America.</p>
<p>We found differences among the groups in terms of health determinants such as smoking and HPV vaccination (see <xref ref-type="table" rid="tab1">Table 1</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S1</xref>). Women in the HSIL groups had a higher prevalence of current smoking than those in the non-HSIL group did and women in the HSIL group were more frequently in a stable relationship compared to those in the non-HSIL group (78% vs. 46%, <italic>p</italic>&#x202F;=&#x202F;0.031). Additionally, there were notable differences in HPV vaccination rates, with a higher percentage of vaccinated women in the non-HSIL group than in the HSIL group.</p>
<table-wrap position="float" id="tab1">
<label>Table 1</label>
<caption>
<p>General characteristics of the study participants dividing the groups in two categories (non-HSIL and HSIL).</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th colspan="2">Variable</th>
<th align="center" valign="top">non-HSIL (<italic>N</italic>&#x202F;=&#x202F;87)</th>
<th align="center" valign="top">HSIL (<italic>N</italic>&#x202F;=&#x202F;18)</th>
<th align="center" valign="top"><italic>p</italic>-value</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">Age, median (IQR)</td>
<td/>
<td align="center" valign="middle">41 (34&#x2013;48)</td>
<td align="center" valign="middle">38.5 (34.5&#x2013;45.5)</td>
<td align="center" valign="middle">0.9</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="3">Nationality, <italic>n</italic> (%)</td>
<td align="center" valign="top">Spain</td>
<td align="center" valign="top">75 (86.2%)</td>
<td align="center" valign="top">13 (72.2%)</td>
<td align="center" valign="middle" rowspan="3">0.136</td>
</tr>
<tr>
<td align="center" valign="top">Eastern Europe</td>
<td align="center" valign="top">4 (4.6%)</td>
<td align="center" valign="top">1 (5.5%)</td>
</tr>
<tr>
<td align="center" valign="top">South America</td>
<td align="center" valign="top">6 (6.9%)</td>
<td align="center" valign="top">4 (22.2%)</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="2">Current smoker, <italic>n</italic> (%)</td>
<td align="center" valign="top">No</td>
<td align="center" valign="top">45 (51.7%)</td>
<td align="center" valign="top">12 (66.7%)</td>
<td align="center" valign="middle" rowspan="2">0.874</td>
</tr>
<tr>
<td align="center" valign="top">Yes</td>
<td align="center" valign="top">14 (16.1%)</td>
<td align="center" valign="top">5 (27.8%)</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="4">Condom use, <italic>n</italic> (%)</td>
<td align="center" valign="top">Never</td>
<td align="center" valign="top">43 (49.4%)</td>
<td align="center" valign="top">13 (72.2%)</td>
<td align="center" valign="middle" rowspan="4">0.892</td>
</tr>
<tr>
<td align="center" valign="top">&#x003C;&#x202F;50%</td>
<td align="center" valign="top">2 (2.3%)</td>
<td align="center" valign="top">0 (0%)</td>
</tr>
<tr>
<td align="center" valign="top">&#x003E;&#x202F;50%</td>
<td align="center" valign="top">4 (4.6%)</td>
<td align="center" valign="top">1 (5.6%)</td>
</tr>
<tr>
<td align="center" valign="top">Always</td>
<td align="center" valign="top">10 (11.5%)</td>
<td align="center" valign="top">3 (16.7%)</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="3">Currently in a couple, <italic>n</italic> (%)</td>
<td align="center" valign="top">No:</td>
<td align="center" valign="top">19 (21.8%)</td>
<td align="center" valign="top">3 (16.7%)</td>
<td align="center" valign="middle" rowspan="3">&#x003C;0.05</td>
</tr>
<tr>
<td align="center" valign="top">Yes</td>
<td align="center" valign="top">40 (46.0%)</td>
<td align="center" valign="top">14 (77.8%)</td>
</tr>
<tr>
<td align="center" valign="top">N/A</td>
<td align="center" valign="top">28 (32.2%)</td>
<td align="center" valign="top">1 (5.6%)</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="2">STIs in the last year, <italic>n</italic> (%)</td>
<td align="center" valign="top">No</td>
<td align="center" valign="top">84 (96.5%)</td>
<td align="center" valign="top">17 (94.4%)</td>
<td align="center" valign="middle" rowspan="2">1.000</td>
</tr>
<tr>
<td align="center" valign="top">Yes</td>
<td align="center" valign="top">2 (2.3%)</td>
<td align="center" valign="top">1 (5.6%)</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="2">HPV vaccine, <italic>n</italic> (%)</td>
<td align="center" valign="top">Yes</td>
<td align="center" valign="top">18 (20.7%)</td>
<td align="center" valign="top">2 (11.1%)</td>
<td align="center" valign="middle" rowspan="2">0.540</td>
</tr>
<tr>
<td align="center" valign="top">No</td>
<td align="center" valign="top">69 (79.3%)</td>
<td align="center" valign="top">16 (88.9%)</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="3">HPV-AR, <italic>n</italic> (%)</td>
<td align="center" valign="top">HPV16/18</td>
<td align="center" valign="top">5 (5.7%)</td>
<td align="center" valign="top">5 (27.7%)</td>
<td align="center" valign="middle" rowspan="3">&#x003C;0.05</td>
</tr>
<tr>
<td align="center" valign="top">Other HR types</td>
<td align="center" valign="top">29 (33.3%)</td>
<td align="center" valign="top">11 (61.1%)</td>
</tr>
<tr>
<td align="center" valign="top">None</td>
<td align="center" valign="top">53 (60.9%)</td>
<td align="center" valign="top">2 (11.1%)</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>STI, Sexually transmitted infections; HPV, Human Papilloma Virus; HSIL, high-grade squamous intraepithelial lesions HR: High Risk.</p>
<p>Values not represented correspond to non-responses from patients. IQR, interquartile range.</p>
</table-wrap-foot>
</table-wrap>
<p>Among HPV-seropositive women, the distribution of HPV types was different between the non-HSIL and HSIL groups (<italic>p</italic>&#x202F;&#x003C;&#x202F;0.05). Most high-risk HPV genotypes were neither 16 nor 18. For these high-risk HPV genotypes, the specific type was not identified by polymerase chain reaction (PCR), which means that it could be any of the following: 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 59, 66, and 68. However, the strong enrichment of HPV16/18 in HSIL is consistent with their established role in cervical carcinogenesis, and supports the validity of our dichotomized grouping strategy for microbiota analyses. For more information, please refer to <xref ref-type="table" rid="tab1">Table 1</xref> and <xref ref-type="supplementary-material" rid="SM1">Supplementary Table S1</xref>.</p>
</sec>
<sec id="sec9">
<title>Microbiota diversity and specific taxa associations in the development of cervical HSIL</title>
<p>Changes in microbiota and metabolism have been documented across different stages of CC. We analyzed the differences between the different stages (Negative, LSIL and HSIL) (see <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure S1</xref>). However, we decided to focus on microbiota analyses by dividing the samples into two general groups: women who developed HSIL (HSIL) and those who did not develop HSIL (non-HSIL) to prioritize a clinically actionable and parsimonious grouping, since only HSIL requires treatment, and to focus on extreme phenotypes, which may help to identify more robust microbial markers.</p>
<p>To measure the richness and evenness of bacterial taxa within a community, we calculated alpha diversity indexes (<xref ref-type="fig" rid="fig1">Figures 1A</xref>,<xref ref-type="fig" rid="fig1">B</xref>). Alpha diversity, calculated using Shannon and Simpson indexes, varied significantly between the HSIL and non-HSIL groups. We also applied a dimensionality-reduction technique, Principal Coordinates Analysis (PCoA), to visualize the differences between microbial community configurations (beta diversity), and we found no clearly distinct clusters (<xref ref-type="fig" rid="fig1">Figure 1C</xref>). To identify specific taxa that differed significantly in abundance between groups, we used the ANCOM-BC2 method to compare the differences in absolute abundances. Women with HSIL showed increased abundances in the genus <italic>Parvimonas, Fastidiosipila,</italic> and <italic>Pseudomonas</italic>, while the genera <italic>Porphyromonas</italic>, <italic>Elizabethkingia</italic>, <italic>Bacteroides,</italic> and <italic>Akkermansia</italic> were depleted (<xref ref-type="fig" rid="fig1">Figure 1D</xref>).</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p>Diversity and abundance of cervicovaginal microbiota in the two compared groups (Non-HSIL and HSIL). <bold>(A)</bold> Boxplots based on Shannon index, showing the maximum, minimum, median, and abnormal values of the index from each group. Wilcoxon two-tailed test <italic>p</italic>-value: 0.036. <bold>(B)</bold> Boxplots based on the Simpson index, showing the maximum, minimum, median, and abnormal values of the index from each group. Wilcoxon two-tailed test <italic>p</italic>-value: 0.024. <bold>(C)</bold> PCoA based on the weighted Unifrac distance showing the two components explaining most of the variance (PC1 and PC2). The multiresponse permutation procedure (MRPP) test was performed to statistically evaluate the differences between the groups (MRPP significance: non-HSIL-HSIL&#x202F;=&#x202F;0.238). <bold>(D)</bold> The ANCOM-BC2 test was performed to analyze the abundance of genera in the two groups. The results were calculated as Log Fold Change (LFC). Red indicates genera with a positive LFC (higher abundance) in the HSIL group. Blue indicates genera with a negative LFC (lower abundance) in the HSIL group.</p>
</caption>
<graphic xlink:href="fmed-12-1615571-g001.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Four charts depicting microbial diversity and abundance. Chart A shows a boxplot of alpha diversity (Shannon index) comparing non-HSL and HSL groups, with HSL having higher diversity. Chart B presents a boxplot of alpha diversity (Simpson index) with a similar trend. Chart C displays a PCoA plot of beta diversity with red squares and blue circles representing two groups, surrounded by confidence ellipses. Chart D is a bar chart of abundance from ANCOM-BC2 analysis, showing log fold changes of specific bacteria, with red indicating positive change and blue negative.</alt-text>
</graphic>
</fig>
</sec>
<sec id="sec10">
<title>Functional analysis of the cervical microbiome reveals key pathways associated with the development of cervical HSIL</title>
<p>Next, we sought to understand the potential mechanisms by which the cervical microbiome influences the development of HSIL. We compared the bacterial imputed functions between the most extreme phenotypes, the non-HSIL and HSIL groups. We also applied the ANCOM-BC2 method, focusing on the KEGG and EC databases, to identify the most significant terms. We identified 37 KO and 17 EC terms that were significantly associated with HSIL (<xref ref-type="fig" rid="fig2">Figure 2</xref>; <xref ref-type="supplementary-material" rid="SM1">Supplementary Tables S2, S3</xref>).</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>Functional terms differentially found in the cervical HSIL group compared to the negative group. <bold>(A)</bold> KEGG terms obtained using ANCOM-BC2 and selected based on their LFC. <bold>(B)</bold> EC terms obtained using ANCOM-BC2 and selected based on their LFC. The genus in red had a positive LFC (higher abundance) in the group with HSIL. The genus in blue had a negative LFC (lower abundance) in the HSIL group. Genus subjected to pseudocount correction by ANCOM-BC2 is shown in green. <bold>(C)</bold> The main terms were obtained using KEGG mapper with the list of significant terms from the ANCOM-BC2 analysis. Numbers in parentheses denote the number of terms in each category. <bold>(D)</bold> Main terms obtained after performing enrichment analysis, as explained in the Methods section, for the KO terms (green) and EC terms (purple) obtained using ANCOM-BC2. The inverse of the <italic>p</italic>-value is represented on the x-axis for visualization purposes.</p>
</caption>
<graphic xlink:href="fmed-12-1615571-g002.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Bar charts and diagrams showing data comparisons between HSIL and non-HSIL groups. Chart A displays significant KO terms, with red bars indicating positive log fold changes and blue negative. Chart B shows significant EC terms, similarly color-coded. Diagram C outlines main KEGG pathways, divided into KO and EC terms, highlighting metabolic pathways and biosynthesis. Chart D visualizes enriched functions in HSIL, with bars for glycine, sulfur, starch, and lysine biosynthesis colored by term type.</alt-text>
</graphic>
</fig>
<p>This two-step approach allows us to first identify relevant genes through differential abundance.</p>
<p>Analysis and interpretation of the functional impact of these genes were performed by mapping them to pathways, reducing noise, and increasing the biological relevance of our results. In women with HSIL, several KO terms related to metabolic pathways, biosynthesis of secondary metabolites, and microbial metabolism were significantly enriched, including those associated with amino sugar and nucleotide metabolism, sulfur metabolism, and starch/sucrose metabolism. On the enzyme side, key EC terms related to lysine biosynthesis, carbohydrate degradation, and glyoxylate metabolism were more abundant, whereas pathways such as glycine/serine/threonine metabolism were depleted.</p>
</sec>
</sec>
<sec sec-type="discussion" id="sec11">
<title>Discussion</title>
<p>Our study revealed a more diverse vaginal microbiota in women with cervical dysplasia than that in healthy women, with differences in specific taxa and biological functions linked to the development of cervical HSIL. Increased microbial diversity in the cervical environment, as observed in women with HSIL, may reflect a shift away from a Lactobacillus-dominated state toward a polymicrobial community enriched in anaerobes and associated with local inflammation. Such transitions have been linked to epithelial barrier disruption, altered immune responses, and HPV persistence, all of which may contribute to lesion progression (<xref ref-type="bibr" rid="ref12">12</xref>, <xref ref-type="bibr" rid="ref13">13</xref>). The increased abundance of <italic>Parvimonas, Fastidiosipila</italic>, and <italic>Pseudomonas</italic> in women with HSIL suggests their potential role in cervical dysbiosis and HPV progression. Notably, <italic>Parvimonas</italic> spp. has been frequently associated with precancerous lesions and dysbiosis (<xref ref-type="bibr" rid="ref7">7</xref>), particularly in connection with sexually transmitted infections such as <italic>Chlamydia trachomatis</italic> and persistent HPV infection. Interestingly, <italic>Pseudomonas,</italic> while often considered part of the healthy endocervical microbiota, has been linked to various pathologies such as <italic>Chlamydia trachomatis</italic> infection in some <italic>in vitro</italic> studies (<xref ref-type="bibr" rid="ref14">14</xref>, <xref ref-type="bibr" rid="ref15">15</xref>), further supporting the idea that shifts in the microbiota may contribute to disease. The prominence of these genera, which are more commonly associated with the intestinal microbiota, suggests possible translocation or shared pathways between the gut and the cervicovaginal environment (<xref ref-type="bibr" rid="ref16">16</xref>, <xref ref-type="bibr" rid="ref17">17</xref>). Mechanisms such as bacterial ascent through the reproductive tract or hematogenous spread from gastrointestinal niches could explain this crossover. However, further studies are required to explore these possibilities.</p>
<p>The functional analysis further reinforces the association between microbial dysbiosis and HSIL. We found enriched gene pathways related to typical bacterial metabolic processes and enzyme pathways related to glycine, serine, and threonine metabolism, suggesting that the microbial community in HSIL cases exhibits heightened metabolic activity. The increase in glycosidases and a corresponding decrease in oxidoreductases also point to significant metabolic alterations, potentially contributing to a dysbiotic environment conducive to lesion progression. These findings are consistent with metabolic dysregulation linked to cancer development, including processes such as the Warburg effect and oxidative stress, which may promote DNA damage and promote the development of cancerous lesions (<xref ref-type="bibr" rid="ref18">18</xref>, <xref ref-type="bibr" rid="ref19">19</xref>).</p>
<p>Our study has some limitations. First, although we observed robust and consistent associations across diversity indices, taxonomic composition, and functional pathways, the relatively small number of HSIL cases (<italic>n</italic>&#x202F;=&#x202F;18) and the amplicon-based sequencing technique may have limited our ability to detect more subtle differences. Nonetheless, the effect sizes observed were sufficient to reach statistical significance despite this limitation. Second, certain behavioral and hormonal factors known to influence the cervicovaginal microbiota&#x2014;such as contraceptive use, recent antibiotic exposure, or sexual behavior&#x2014;were not systematically collected, which may represent sources of residual confounding. While our findings are representative of a well-defined Spanish population, we acknowledge the need for larger, multicenter, and multiethnic studies to confirm and extend these results. Importantly, the microbial shifts we observed are consistent with reports from other settings, supporting the biological plausibility and relevance of our findings. Although our study was not designed to assess the clinical utility of microbiota-based markers, the microbial and functional differences observed may help inform future efforts to develop adjunctive tools for risk stratification. We have previously shown the predictive value of microbiome and metabolomic signatures for anal HSIL progression in a prospective cohort (<xref ref-type="bibr" rid="ref20">20</xref>), and similar strategies are now being planned for cervical dysplasia.</p>
<p>Microbiome-based interventions, particularly intravaginal administration of <italic>Lactobacillus crispatus</italic> and long-term use of <italic>Lactobacillus rhamnosus</italic>, have demonstrated potential in improving cervical cytology and promoting high-risk HPV clearance in small randomized or controlled trials (<xref ref-type="bibr" rid="ref20">20</xref>, <xref ref-type="bibr" rid="ref21">21</xref>). While these findings are promising, they remain limited by small sample sizes and short follow-up durations. Larger, rigorously designed randomized controlled trials are urgently needed to confirm efficacy and guide clinical implementation. Additionally, most studies examining the relationship between microbiota and cervical dysplasia have focused on vaginal rather than cervical samples (<xref ref-type="bibr" rid="ref21">21</xref>) or have reported mixed results from different anatomical sites (<xref ref-type="bibr" rid="ref22">22</xref>). Future studies correlating both environments could provide a less invasive means of assessing women at a higher risk of progressing to HSIL by linking findings in the vaginal microbiota to those in the cervix.</p>
<p>Geographic location may also influence microbiota composition (<xref ref-type="bibr" rid="ref23">23</xref>). To the best of our knowledge, this study is the first to examine the link between the cervical microbiota and HSIL in women from Spain, and direct comparisons with the existing literature are challenging. This underscores the need for region-specific studies to better understand the microbial landscape and its relationship to cervical dysplasia. Lastly, although our functional analysis identified significant pathways, the presence of genes does not guarantee their activity, highlighting the importance of further experimental studies, particularly longitudinal studies, to establish causal relationships between cervical dysbiosis and CIN progression.</p>
<p>We acknowledge that cross-sectional studies such as ours cannot determine causality. Although key exclusion criteria were applied to reduce major confounders, we acknowledge that many other factors influencing microbiota composition&#x2014;such as diet, probiotic use, or menstrual cycle phase&#x2014;were not systematically collected and represent a limitation of this study. In our regrouped analysis (HSIL vs. non-HSIL), being currently in a couple was more frequent among women with HSIL, however, while this finding may reflect unmeasured aspects of sexual behavior or partner-related factors, it should be interpreted cautiously given the lack of detailed behavioral data. These limitations have been broadly discussed in our recent review (<xref ref-type="bibr" rid="ref24">24</xref>), where we emphasized the need for complementary experimental models &#x2014;such as experimental models, multi-omics approaches, and longitudinal cohort designs&#x2014;to better elucidate the causal relationships between the microbiota and viral pathogenesis. Although we observed differences in variables such as smoking, age, or nationality across groups, our study was not designed or powered to disentangle their individual effects on microbiota composition. Future hypothesis-driven studies with larger sample sizes and multivariable modeling will be required to evaluate the independent contribution of these host factors.</p>
<p>Despite these limitations, our study lays the groundwork for future research to explore the role of the cervicovaginal microbiome in HSIL development. The microbial taxa and metabolic pathways identified here could inform studies focused on identifying markers for cervical dysplasia progression, offering avenues for targeted screening and therapeutic interventions.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="sec12">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/<xref ref-type="supplementary-material" rid="SM1">Supplementary material</xref>.</p>
</sec>
<sec sec-type="ethics-statement" id="sec13">
<title>Ethics statement</title>
<p>The studies involving humans were approved by Independent Ethics Committee of the coordinating center Hospital Ram&#x00F3;n y Cajal (<email>ceic.hrc@salud.madrid.org</email>, approval 146-21l). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.</p>
</sec>
<sec sec-type="author-contributions" id="sec14">
<title>Author contributions</title>
<p>MR: Investigation, Writing &#x2013; review &#x0026; editing, Resources, Writing &#x2013; original draft, Methodology. EM: Writing &#x2013; review &#x0026; editing, Project administration. CV: Supervision, Methodology, Writing &#x2013; review &#x0026; editing, Resources. MS: Writing &#x2013; review &#x0026; editing, Methodology, Resources. AM: Resources, Writing &#x2013; review &#x0026; editing, Validation. LL: Project administration, Writing &#x2013; review &#x0026; editing, Visualization. CD-G: Software, Supervision, Writing &#x2013; review &#x0026; editing. CT: Investigation, Writing &#x2013; review &#x0026; editing, Software, Formal analysis. AA: Data curation, Writing &#x2013; review &#x0026; editing. SM: Formal analysis, Writing &#x2013; review &#x0026; editing. MS-C: Writing &#x2013; review &#x0026; editing, Visualization, Validation. SS-V: Formal analysis, Visualization, Writing &#x2013; review &#x0026; editing, Supervision, Funding acquisition.</p>
</sec>
<sec sec-type="funding-information" id="sec15">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research and/or publication of this article. This work was funded by the Instituto de Salud Carlos III and Fondos FEDER, Acci&#x00F3;n Estrat&#x00E9;gica en Salud (ICI20/00058, PI21/00041, PI24/00078, FORT23/00046), co-financed by the European Development Regional Fund &#x2018;A way to achieve Europe&#x2019; (ERDF), a Gilead Fellowship (GLD23/0059), and CIBERINFEC (IM23/INFEC/4), and the European Comission (EDCTP program, RIA2018D-2511).</p>
</sec>
<ack>
<p>We thank all the patients and healthcare workers who participated in the study.</p>
</ack>
<sec sec-type="COI-statement" id="sec16">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec sec-type="ai-statement" id="sec17">
<title>Generative AI statement</title>
<p>The authors declare that no Gen AI was used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p>
</sec>
<sec sec-type="disclaimer" id="sec18">
<title>Publisher&#x2019;s note</title>
<p>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.</p>
</sec>
<sec sec-type="supplementary-material" id="sec19">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fmed.2025.1615571/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fmed.2025.1615571/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Supplementary_file_1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1"><label>1.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Onyango</surname> <given-names>S</given-names></name> <name><surname>Mi</surname> <given-names>JD</given-names></name> <name><surname>Koech</surname> <given-names>A</given-names></name> <name><surname>Okiro</surname> <given-names>P</given-names></name> <name><surname>Temmerman</surname> <given-names>M</given-names></name> <name><surname>Von Dadelszen</surname> <given-names>P</given-names></name> <etal/></person-group>. <article-title>Microbiota dynamics, metabolic and immune interactions in the cervicovaginal environment and their role in spontaneous preterm birth</article-title>. <source>Front Immunol</source>. (<year>2023</year>) <volume>14</volume>:<fpage>1306473</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2023.1306473</pub-id>, PMID: <pub-id pub-id-type="pmid">38196946</pub-id></citation></ref>
<ref id="ref2"><label>2.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Campisciano</surname> <given-names>G</given-names></name> <name><surname>Iebba</surname> <given-names>V</given-names></name> <name><surname>Zito</surname> <given-names>G</given-names></name> <name><surname>Luppi</surname> <given-names>S</given-names></name> <name><surname>Martinelli</surname> <given-names>M</given-names></name> <name><surname>Fischer</surname> <given-names>L</given-names></name> <etal/></person-group>. <article-title><italic>Lactobacillus iners</italic> and <italic>gasseri</italic>, <italic>Prevotella bivia</italic> and HPV belong to the microbiological signature negatively affecting human reproduction</article-title>. <source>Microorganisms</source>. (<year>2020</year>) <volume>9</volume>:<fpage>39</fpage>. doi: <pub-id pub-id-type="doi">10.3390/microorganisms9010039</pub-id>, PMID: <pub-id pub-id-type="pmid">33375526</pub-id></citation></ref>
<ref id="ref3"><label>3.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hu</surname> <given-names>M</given-names></name> <name><surname>Yang</surname> <given-names>W</given-names></name> <name><surname>Yan</surname> <given-names>R</given-names></name> <name><surname>Chi</surname> <given-names>J</given-names></name> <name><surname>Xia</surname> <given-names>Q</given-names></name> <name><surname>Yang</surname> <given-names>Y</given-names></name> <etal/></person-group>. <article-title>Co-evolution of vaginal microbiome and cervical cancer</article-title>. <source>J Transl Med</source>. (<year>2024</year>) <volume>22</volume>:<fpage>559</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12967-024-05265-w</pub-id>, PMID: <pub-id pub-id-type="pmid">38863033</pub-id></citation></ref>
<ref id="ref4"><label>4.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Amabebe</surname> <given-names>E</given-names></name> <name><surname>Anumba</surname> <given-names>DOC</given-names></name></person-group>. <article-title>The vaginal microenvironment: the physiologic role of lactobacilli</article-title>. <source>Front Med</source>. (<year>2018</year>) <volume>5</volume>:<fpage>181</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmed.2018.00181</pub-id>, PMID: <pub-id pub-id-type="pmid">29951482</pub-id></citation></ref>
<ref id="ref5"><label>5.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Santella</surname> <given-names>B</given-names></name> <name><surname>Schettino</surname> <given-names>MT</given-names></name> <name><surname>Franci</surname> <given-names>G</given-names></name> <name><surname>De Franciscis</surname> <given-names>P</given-names></name> <name><surname>Colacurci</surname> <given-names>N</given-names></name> <name><surname>Schiattarella</surname> <given-names>A</given-names></name> <etal/></person-group>. <article-title>Microbiota and HPV: the role of viral infection on vaginal microbiota</article-title>. <source>J Med Virol</source>. (<year>2022</year>) <volume>94</volume>:<fpage>4478</fpage>&#x2013;<lpage>84</lpage>. doi: <pub-id pub-id-type="doi">10.1002/jmv.27837</pub-id>, PMID: <pub-id pub-id-type="pmid">35527233</pub-id></citation></ref>
<ref id="ref6"><label>6.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>Y</given-names></name> <name><surname>Wang</surname> <given-names>S</given-names></name> <name><surname>Liu</surname> <given-names>J</given-names></name> <name><surname>Su</surname> <given-names>M</given-names></name> <name><surname>Diao</surname> <given-names>X</given-names></name> <name><surname>Liang</surname> <given-names>X</given-names></name> <etal/></person-group>. <article-title>Characteristics of vaginal microbiota in various cervical intraepithelial neoplasia: a cross-sectional study</article-title>. <source>J Transl Med</source>. (<year>2023</year>) <volume>21</volume>:<fpage>816</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12967-023-04676-5</pub-id>, PMID: <pub-id pub-id-type="pmid">37974192</pub-id></citation></ref>
<ref id="ref7"><label>7.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peremykina</surname> <given-names>A</given-names></name> <name><surname>Cheranev</surname> <given-names>V</given-names></name> <name><surname>Krivoy</surname> <given-names>A</given-names></name> <name><surname>Andreev</surname> <given-names>AO</given-names></name> <name><surname>Repinskaia</surname> <given-names>Z</given-names></name> <name><surname>Asaturova</surname> <given-names>AV</given-names></name> <etal/></person-group>. <article-title>Microbiome markers in HPV-positive and HPV-negative women of reproductive age with ASCUS and SIL determined by V4 region of 16S rRNA gene sequencing</article-title>. <source>Front Microbiol</source>. (<year>2024</year>) <volume>15</volume>:<fpage>1334502</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2024.1334502</pub-id>, PMID: <pub-id pub-id-type="pmid">38550866</pub-id></citation></ref>
<ref id="ref8"><label>8.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sasivimolrattana</surname> <given-names>T</given-names></name> <name><surname>Chantratita</surname> <given-names>W</given-names></name> <name><surname>Sensorn</surname> <given-names>I</given-names></name> <name><surname>Chaiwongkot</surname> <given-names>A</given-names></name> <name><surname>Oranratanaphan</surname> <given-names>S</given-names></name> <name><surname>Bhattarakosol</surname> <given-names>P</given-names></name> <etal/></person-group>. <article-title>Cervical microbiome in women infected with HPV16 and high-risk HPVs</article-title>. <source>Int J Environ Res Public Health</source>. (<year>2022</year>) <volume>19</volume>:<fpage>14716</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ijerph192214716</pub-id>, PMID: <pub-id pub-id-type="pmid">36429432</pub-id></citation></ref>
<ref id="ref9"><label>9.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ata</surname> <given-names>B</given-names></name> <name><surname>Yildiz</surname> <given-names>S</given-names></name> <name><surname>Turkgeldi</surname> <given-names>E</given-names></name> <name><surname>Brocal</surname> <given-names>VP</given-names></name> <name><surname>Dinleyici</surname> <given-names>EC</given-names></name> <name><surname>Moya</surname> <given-names>A</given-names></name> <etal/></person-group>. <article-title>The Endobiota study: comparison of vaginal, cervical and gut microbiota between women with stage 3/4 endometriosis and healthy controls</article-title>. <source>Sci Rep</source>. (<year>2019</year>) <volume>9</volume>:<fpage>2204</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-019-39700-6</pub-id>, PMID: <pub-id pub-id-type="pmid">30778155</pub-id></citation></ref>
<ref id="ref10"><label>10.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Camargo</surname> <given-names>M</given-names></name> <name><surname>Vega</surname> <given-names>L</given-names></name> <name><surname>Mu&#x00F1;oz</surname> <given-names>M</given-names></name> <name><surname>S&#x00E1;nchez</surname> <given-names>R</given-names></name> <name><surname>Patarroyo</surname> <given-names>ME</given-names></name> <name><surname>Ram&#x00ED;rez</surname> <given-names>JD</given-names></name> <etal/></person-group>. <article-title>Changes in the cervical microbiota of women with different high-risk human papillomavirus loads</article-title>. <source>Viruses</source>. (<year>2022</year>) <volume>14</volume>:<fpage>2674</fpage>. doi: <pub-id pub-id-type="doi">10.3390/v14122674</pub-id>, PMID: <pub-id pub-id-type="pmid">36560678</pub-id></citation></ref>
<ref id="ref11"><label>11.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lin</surname> <given-names>H</given-names></name> <name><surname>Peddada</surname> <given-names>SD</given-names></name></person-group>. <article-title>Analysis of compositions of microbiomes with bias correction</article-title>. <source>Nat Commun</source>. (<year>2020</year>) <volume>11</volume>:<fpage>3514</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41467-020-17041-7</pub-id>, PMID: <pub-id pub-id-type="pmid">32665548</pub-id></citation></ref>
<ref id="ref12"><label>12.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Harris</surname> <given-names>PA</given-names></name> <name><surname>Taylor</surname> <given-names>R</given-names></name> <name><surname>Thielke</surname> <given-names>R</given-names></name> <name><surname>Payne</surname> <given-names>J</given-names></name> <name><surname>Gonzalez</surname> <given-names>N</given-names></name> <name><surname>Conde</surname> <given-names>JG</given-names></name></person-group>. <article-title>Research electronic data capture (REDCap)&#x2014;a metadata-driven methodology and workflow process for providing translational research informatics support</article-title>. <source>J Biomed Inform</source>. (<year>2009</year>) <volume>42</volume>:<fpage>377</fpage>&#x2013;<lpage>81</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jbi.2008.08.010</pub-id>, PMID: <pub-id pub-id-type="pmid">18929686</pub-id></citation></ref>
<ref id="ref13"><label>13.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>M</given-names></name> <name><surname>Yang</surname> <given-names>J</given-names></name> <name><surname>Zhou</surname> <given-names>L</given-names></name> <name><surname>Zhang</surname> <given-names>J</given-names></name> <name><surname>Li</surname> <given-names>Y</given-names></name> <name><surname>Chen</surname> <given-names>J</given-names></name> <etal/></person-group>. <article-title>Efficacy of a novel affitoxin targeting major outer membrane protein against <italic>Chlamydia trachomatis</italic> in vitro and in vivo</article-title>. <source>J Infect Dis</source>. (<year>2024</year>) <volume>230</volume>:<fpage>1476</fpage>&#x2013;<lpage>87</lpage>. doi: <pub-id pub-id-type="doi">10.1093/infdis/jiae257</pub-id>, PMID: <pub-id pub-id-type="pmid">38723186</pub-id></citation></ref>
<ref id="ref14"><label>14.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Erridge</surname> <given-names>C</given-names></name> <name><surname>Pridmore</surname> <given-names>A</given-names></name> <name><surname>Eley</surname> <given-names>A</given-names></name> <name><surname>Stewart</surname> <given-names>J</given-names></name> <name><surname>Poxton</surname> <given-names>IR</given-names></name></person-group>. <article-title>Lipopolysaccharides of <italic>Bacteroides fragilis</italic>, <italic>Chlamydia trachomatis</italic> and <italic>Pseudomonas aeruginosa</italic> signal via toll-like receptor 2</article-title>. <source>J Med Microbiol</source>. (<year>2004</year>) <volume>53</volume>:<fpage>735</fpage>&#x2013;<lpage>40</lpage>. doi: <pub-id pub-id-type="doi">10.1099/jmm.0.45598-0</pub-id>, PMID: <pub-id pub-id-type="pmid">15272059</pub-id></citation></ref>
<ref id="ref15"><label>15.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Howe</surname> <given-names>SE</given-names></name> <name><surname>Shillova</surname> <given-names>N</given-names></name> <name><surname>Konjufca</surname> <given-names>V</given-names></name></person-group>. <article-title>Dissemination of Chlamydia from the reproductive tract to the gastro-intestinal tract occurs in stages and relies on Chlamydia transport by host cells</article-title>. <source>PLoS Pathog</source>. (<year>2019</year>) <volume>15</volume>:<fpage>e1008207</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.ppat.1008207</pub-id>, PMID: <pub-id pub-id-type="pmid">31790512</pub-id></citation></ref>
<ref id="ref16"><label>16.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Campisciano</surname> <given-names>G</given-names></name> <name><surname>Zanotta</surname> <given-names>N</given-names></name> <name><surname>Quadrifoglio</surname> <given-names>M</given-names></name> <name><surname>Careri</surname> <given-names>A</given-names></name> <name><surname>Torresani</surname> <given-names>A</given-names></name> <name><surname>Cason</surname> <given-names>C</given-names></name> <etal/></person-group>. <article-title>The bacterial DNA profiling of chorionic villi and amniotic fluids reveals overlaps with maternal Oral, vaginal, and gut microbiomes</article-title>. <source>Int J Mol Sci</source>. (<year>2023</year>) <volume>24</volume>:<fpage>2873</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ijms24032873</pub-id>, PMID: <pub-id pub-id-type="pmid">36769194</pub-id></citation></ref>
<ref id="ref17"><label>17.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Re&#x017E;en</surname> <given-names>T</given-names></name> <name><surname>Rozman</surname> <given-names>D</given-names></name> <name><surname>Kov&#x00E1;cs</surname> <given-names>T</given-names></name> <name><surname>Kov&#x00E1;cs</surname> <given-names>P</given-names></name> <name><surname>Sipos</surname> <given-names>A</given-names></name> <name><surname>Bai</surname> <given-names>P</given-names></name> <etal/></person-group>. <article-title>The role of bile acids in carcinogenesis</article-title>. <source>Cell Mol Life Sci</source>. (<year>2022</year>) <volume>79</volume>:<fpage>243</fpage>. doi: <pub-id pub-id-type="doi">10.1007/s00018-022-04278-2</pub-id>, PMID: <pub-id pub-id-type="pmid">35429253</pub-id></citation></ref>
<ref id="ref18"><label>18.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Manos</surname> <given-names>J</given-names></name></person-group>. <article-title>The human microbiome in disease and pathology</article-title>. <source>APMIS</source>. (<year>2022</year>) <volume>130</volume>:<fpage>690</fpage>&#x2013;<lpage>705</lpage>. doi: <pub-id pub-id-type="doi">10.1111/apm.13225</pub-id>, PMID: <pub-id pub-id-type="pmid">35393656</pub-id></citation></ref>
<ref id="ref19"><label>19.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Serrano-Villar</surname> <given-names>S</given-names></name> <name><surname>Tincati</surname> <given-names>C</given-names></name> <name><surname>Raju</surname> <given-names>SC</given-names></name> <name><surname>S&#x00E1;enz</surname> <given-names>JS</given-names></name> <name><surname>Moreno</surname> <given-names>E</given-names></name> <name><surname>Bargiela</surname> <given-names>R</given-names></name> <etal/></person-group>. <article-title>Microbiome-derived cobalamin and succinyl-CoA as biomarkers for improved screening of anal cancer</article-title>. <source>Nat Med</source>. (<year>2023</year>) <volume>29</volume>:<fpage>1738</fpage>&#x2013;<lpage>49</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41591-023-02407-3</pub-id>, PMID: <pub-id pub-id-type="pmid">37464040</pub-id></citation></ref>
<ref id="ref20"><label>20.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ou</surname> <given-names>Y-C</given-names></name> <name><surname>Fu</surname> <given-names>H-C</given-names></name> <name><surname>Tseng</surname> <given-names>C-W</given-names></name> <name><surname>Wu</surname> <given-names>C-H</given-names></name> <name><surname>Tsai</surname> <given-names>C-C</given-names></name> <name><surname>Lin</surname> <given-names>H</given-names></name></person-group>. <article-title>The influence of probiotics on genital high-risk human papilloma virus clearance and quality of cervical smear: a randomized placebo-controlled trial</article-title>. <source>BMC Womens Health</source>. (<year>2019</year>) <volume>19</volume>:<fpage>103</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12905-019-0798-y</pub-id>, PMID: <pub-id pub-id-type="pmid">31340789</pub-id></citation></ref>
<ref id="ref21"><label>21.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>McKee</surname> <given-names>KS</given-names></name> <name><surname>Carter</surname> <given-names>KA</given-names></name> <name><surname>Bassis</surname> <given-names>C</given-names></name> <name><surname>Young</surname> <given-names>VB</given-names></name> <name><surname>Reed</surname> <given-names>B</given-names></name> <name><surname>Harper</surname> <given-names>DM</given-names></name> <etal/></person-group>. <article-title>The vaginal microbiota, high-risk human papillomavirus infection, and cervical cytology: results from a population-based study</article-title>. <source>Gynecol Pelvic Med</source>. (<year>2020</year>) <volume>3</volume>:<fpage>18</fpage>. doi: <pub-id pub-id-type="doi">10.21037/gpm-20-10</pub-id>, PMID: <pub-id pub-id-type="pmid">35252846</pub-id></citation></ref>
<ref id="ref22"><label>22.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Alizhan</surname> <given-names>D</given-names></name> <name><surname>Ukybassova</surname> <given-names>T</given-names></name> <name><surname>Bapayeva</surname> <given-names>G</given-names></name> <name><surname>Aimagambetova</surname> <given-names>G</given-names></name> <name><surname>Kongrtay</surname> <given-names>K</given-names></name> <name><surname>Kamzayeva</surname> <given-names>N</given-names></name> <etal/></person-group>. <article-title>Cervicovaginal microbiome: physiology, age-related changes, and protective role against human papillomavirus infection</article-title>. <source>J Clin Med</source>. (<year>2025</year>) <volume>14</volume>:<fpage>1521</fpage>. doi: <pub-id pub-id-type="doi">10.3390/jcm14051521</pub-id>, PMID: <pub-id pub-id-type="pmid">40094958</pub-id></citation></ref>
<ref id="ref23"><label>23.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yao</surname> <given-names>T</given-names></name> <name><surname>Wang</surname> <given-names>Z</given-names></name> <name><surname>Liang</surname> <given-names>X</given-names></name> <name><surname>Liu</surname> <given-names>C</given-names></name> <name><surname>Yu</surname> <given-names>Z</given-names></name> <name><surname>Han</surname> <given-names>X</given-names></name> <etal/></person-group>. <article-title>Signatures of vaginal microbiota by 16S rRNA gene: potential bio-geographical application in Chinese Han from three regions of China</article-title>. <source>Int J Legal Med</source>. (<year>2021</year>) <volume>135</volume>:<fpage>1213</fpage>&#x2013;<lpage>24</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00414-021-02525-y</pub-id>, PMID: <pub-id pub-id-type="pmid">33594458</pub-id></citation></ref>
<ref id="ref24"><label>24.</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Moreno</surname> <given-names>E</given-names></name> <name><surname>Ron</surname> <given-names>R</given-names></name> <name><surname>Serrano-Villar</surname> <given-names>S</given-names></name></person-group>. <article-title>The microbiota as a modulator of mucosal inflammation and HIV/HPV pathogenesis: from association to causation</article-title>. <source>Front Immunol</source>. (<year>2023</year>) <volume>14</volume>:<fpage>1072655</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2023.1072655</pub-id>, PMID: <pub-id pub-id-type="pmid">36756132</pub-id></citation></ref>
</ref-list>
</back>
</article>