<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!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="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Cell. Infect. Microbiol.</journal-id>
<journal-title>Frontiers in Cellular and Infection Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Cell. Infect. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">2235-2988</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fcimb.2022.884272</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Cellular and Infection Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Altered Vaginal Microbiota Composition Correlates With Human Papillomavirus and Mucosal Immune Responses in Women With Symptomatic Cervical Ectopy</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>L&#xf3;pez-Filloy</surname><given-names>Mariana</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn004"><sup>&#x2021;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1696431"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Cortez</surname><given-names>Flor J.</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn004"><sup>&#x2021;</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Gheit</surname><given-names>Tarik</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/527499"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Cruz y Cruz</surname><given-names>Omar</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Cruz-Talonia</surname><given-names>Fernando</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ch&#xe1;vez-Torres</surname><given-names>Monserrat</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Arteaga-G&#xf3;mez</surname><given-names>Cristina</given-names>
</name>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Mancilla-Herrera</surname><given-names>Ismael</given-names>
</name>
<xref ref-type="aff" rid="aff6"><sup>6</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/531487"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Montesinos</surname><given-names>Juan J.</given-names>
</name>
<xref ref-type="aff" rid="aff7"><sup>7</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1017117"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Cort&#xe9;s-Morales</surname><given-names>V&#xed;ctor Adri&#xe1;n</given-names>
</name>
<xref ref-type="aff" rid="aff7"><sup>7</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/499545"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Aguilar</surname><given-names>Cecilia</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Tommasino</surname><given-names>Massimo</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="author-notes" rid="fn003"><sup>&#x2020;</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Pinto-Cardoso</surname><given-names>Sandra</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>*</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1765167"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Rocha-Zavaleta</surname><given-names>Leticia</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>*</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/854389"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Departamento de Biolog&#xed;a Molecular y Biotecnolog&#xed;a, Instituto de Investigaciones Biom&#xe9;dicas, Universidad Nacional Aut&#xf3;noma de M&#xe9;xico (UNAM)</institution>, <addr-line>Ciudad de M&#xe9;xico</addr-line>, <country>Mexico</country></aff>
<aff id="aff2"><sup>2</sup><institution>International Agency for Research on Cancer</institution>, <addr-line>Lyon</addr-line>, <country>France</country></aff>
<aff id="aff3"><sup>3</sup><institution>Cl&#xed;nica de Colposcopia Fundaci&#xf3;n &#x201c;Cruz Talonia&#x201d;</institution>, <addr-line>Ciudad de Mexico</addr-line>, <country>Mexico</country></aff>
<aff id="aff4"><sup>4</sup><institution>Centro de Investigaci&#xf3;n en Enfermedades Infecciosas, Instituto Nacional de Enfermedades Respiratorias Ismael Cos&#xed;o Villegas</institution>, <addr-line>Ciudad de M&#xe9;xico</addr-line>, <country>Mexico</country></aff>
<aff id="aff5"><sup>5</sup><institution>Deparatamento de Oncolog&#xed;a, Instituto Nacional de Perinatolog&#xed;a Isidro Espinosa de los Reyes</institution>, <addr-line>Ciudad de M&#xe9;xico</addr-line>, <country>Mexico</country></aff>
<aff id="aff6"><sup>6</sup><institution>Departamento de Infectolog&#xed;a e Inmunolog&#xed;a, Instituto Nacional de Perinatolog&#xed;a Isidro Espinosa de los Reyes</institution>, <addr-line>Ciudad de M&#xe9;xico</addr-line>, <country>Mexico</country></aff>
<aff id="aff7"><sup>7</sup><institution>Mesenchymal Stem Cells Laboratory, Oncology Research Unit, Oncology Hospital, National Medical Center, Instituto Mexicano del Seguro Social (IMSS)</institution>, <addr-line>Mexico City</addr-line>, <country>Mexico</country></aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Leticia Reyes, University of Wisconsin-Madison, United States</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Martin James Holland, University of London, United Kingdom; Jose Eleuterio Junior, Federal University of Ceara, Brazil</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Leticia Rocha-Zavaleta, <email xlink:href="mailto:lrochaz@iibiomedicas.unam.mx">lrochaz@iibiomedicas.unam.mx</email>; Sandra Pinto-Cardoso, <email xlink:href="mailto:sandra.pintocardoso.cieni@gmail.com">sandra.pintocardoso.cieni@gmail.com</email>
</p>
</fn>
<fn fn-type="present-address" id="fn003">
<p>&#x2020;Present address: Massimo Tommasino, Dipartimento di Farmacia - Scienze del Farmaco, University of Bari, Bari, Italy</p>
</fn>
<fn fn-type="equal" id="fn004">
<p>&#x2021;These authors have contributed equally to this work</p>
</fn>
<fn fn-type="other" id="fn002">
<p>This article was submitted to Microbiome in Health and Disease, a section of the journal Frontiers in Cellular and Infection Microbiology</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>17</day>
<month>05</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>12</volume>
<elocation-id>884272</elocation-id>
<history>
<date date-type="received">
<day>26</day>
<month>02</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>19</day>
<month>04</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 L&#xf3;pez-Filloy, Cortez, Gheit, Cruz y Cruz, Cruz-Talonia, Ch&#xe1;vez-Torres, Arteaga-G&#xf3;mez, Mancilla-Herrera, Montesinos, Cort&#xe9;s-Morales, Aguilar, Tommasino, Pinto-Cardoso and Rocha-Zavaleta</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>L&#xf3;pez-Filloy, Cortez, Gheit, Cruz y Cruz, Cruz-Talonia, Ch&#xe1;vez-Torres, Arteaga-G&#xf3;mez, Mancilla-Herrera, Montesinos, Cort&#xe9;s-Morales, Aguilar, Tommasino, Pinto-Cardoso and Rocha-Zavaleta</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>Cervical ectopy is a benign condition of the lower genital tract that is frequently detected in women of reproductive age. Although cervical ectopy is regarded as a physiological condition, some women experience symptoms such as leucorrhoea, persistent bleeding and recurrent vaginal infections that require medical intervention. Cervical ectopy has not been linked to cervical cancer, but it is thought to facilitate the acquisition of sexually transmitted diseases (STDs), like Human Papillomavirus (HPV) infection, as it provides a favorable microenvironment for virus infection and dissemination. We and others have described the presence of oncogenic HPV types in women with symptomatic cervical ectopy. The relevance of this finding and the impact of symptomatic cervical ectopy on the cervicovaginal microenvironment (vaginal microbiota, immune and inflammatory responses) are currently unknown. To shed some light into the interplay between HPV, the vaginal microbiota and mucosal immune and inflammatory responses in the context of this condition, we enrolled 156 women with symptomatic cervical ectopy and determined the presence of HPV using a type-specific multiplex genotyping assay. Overall, HPV was detected in 54.48% women, oncogenic HPV types were found in more than 90% of HPV-positive cases. The most prevalent HPV types were HPV16 (29.4%), HPV31 (21.17%) and HPV18 (15.29%). Next, we evaluated the vaginal microbial composition and diversity by 16S rDNA sequencing, and quantified levels of cytokines and chemokines by flow cytometry using bead-based multiplex assays in a sub-cohort of 63 women. IL-21 and CXCL9 were significantly upregulated in HPV-positive women (<italic>p</italic>=0.0002 and <italic>p</italic>=0.013, respectively). Women with symptomatic cervical ectopy and HPV infection had increased diversity (<italic>p</italic>&lt;0.001), and their vaginal microbiota was enriched in bacterial vaginosis-associated anaerobes (<italic>Sneathia</italic>, <italic>Shuttleworthia</italic>, <italic>Prevotella</italic>, and <italic>Atopobium</italic>) and depleted in <italic>Lactobacillus</italic> spp. Furthermore, the vaginal microbiota of women with symptomatic cervical ectopy and HPV infection correlated with vaginal inflammation (IL-1&#x3b2;, rho=0.56, <italic>p</italic>=0.0004) and increased mucosal homeostatic response (IL-22, rho=0.60, <italic>p</italic>=0.0001). Taken together, our results suggest that HPV infection and dysbiotic vaginal communities could favor a vaginal microenvironment that might delay the recovery of the cervical epithelium in women with symptomatic cervical ectopy and favor STDs acquisition.</p>
</abstract>
<kwd-group>
<kwd>human papillomavirus</kwd>
<kwd>vaginal microbiota</kwd>
<kwd>mucosal immune responses</kwd>
<kwd>genital inflammation</kwd>
<kwd>cervicovaginal mucus</kwd>
<kwd>cervical ectopy</kwd>
</kwd-group>
<contract-sponsor id="cn001">Instituto de Investigaciones Biom&#xe9;dicas, Universidad Nacional Aut&#xf3;noma de M&#xe9;xico<named-content content-type="fundref-id">10.13039/501100020704</named-content>
</contract-sponsor>
<counts>
<fig-count count="5"/>
<table-count count="3"/>
<equation-count count="0"/>
<ref-count count="61"/>
<page-count count="15"/>
<word-count count="7517"/>
</counts>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<title>Introduction</title>
<p>Over the last decades or so, female reproductive health has been at the forefront of research as strong evidence supports bidirectional relationships between the cervicovaginal microenvironment and cervical cancer, viral acquisition and persistence, gynecologic and obstetric diseases and other benign conditions like cervical ectopy (<xref ref-type="bibr" rid="B34">Laniewski et&#xa0;al., 2020</xref>). During adolescence, pregnancy, and in women who take hormonal contraceptives, eversion of the endocervix exposes the columnar epithelium to the vaginal environment, forming a cervical ectopy (<xref ref-type="bibr" rid="B27">Jacobson et&#xa0;al., 2000</xref>; <xref ref-type="bibr" rid="B51">Reich and Fritsch, 2014</xref>). Cervical ectopy tissue is eventually replaced by newly-formed squamous epithelium through an irreversible physiological process of squamous metaplasia (<xref ref-type="bibr" rid="B27">Jacobson et&#xa0;al., 2000</xref>). Cervical ectopy is generally considered a benign condition, however, there is evidence linking cervical ectopy to an increased risk of acquisition of sexually transmitted diseases (STDs), like <italic>Chlamydia trachomatis</italic> (<xref ref-type="bibr" rid="B35">Lee et&#xa0;al., 2006</xref>), Human Immunodeficiency Virus (HIV) (<xref ref-type="bibr" rid="B42">Myer et&#xa0;al., 2006</xref>), and Human Papillomavirus (HPV) (<xref ref-type="bibr" rid="B41">Monroy et&#xa0;al., 2010</xref>). In a previous work, we reported that a high proportion of women with cervical ectopy, who went on to receive treatment, were positive for the presence of high-risk HPV types (HR-HPV) (<xref ref-type="bibr" rid="B41">Monroy et&#xa0;al., 2010</xref>). Persistent infection with HR-HPV has been associated with cervical cancer (<xref ref-type="bibr" rid="B58">Walboomers et&#xa0;al., 1999</xref>). The clinical relevance of HPV infection in women with cervical ectopy remains unclear as cervical ectopy is not considered a malignant condition. Also, a number of women with cervical ectopy suffer from bothersome symptoms like abundant leucorrhoea, postcoital and intermenstrual bleeding, pelvic pain and recurrent cervicitis, that can impact their quality of life (<xref ref-type="bibr" rid="B21">Harry et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B23">Hua et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B14">Cekmez et&#xa0;al., 2016</xref>). The presence of chronic symptoms, along with the risk of further infections, warrants treatment of cervical ectopy, mainly cauterization of the area (<xref ref-type="bibr" rid="B37">Machado et&#xa0;al., 2008</xref>).</p>
<p>The natural history of cervical ectopy is influenced by the vaginal microenvironment and modified by the presence of reproductive hormones (estrogen), trauma (sexual activity), and the pH of the cervicovaginal mucus. Indeed, squamous metaplastic activity is stimulated by the acidic pH of the vagina (<xref ref-type="bibr" rid="B27">Jacobson et&#xa0;al., 2000</xref>), which is in turn regulated by the vaginal microbiota. Low vaginal pH (&lt; 4.5) is maintained by the activity of most species of <italic>Lactobacillus</italic> (<italic>L. crispatus</italic>, <italic>L. iners</italic>, <italic>L. jensenii</italic>, and <italic>L. gasseri</italic>), which colonize the cervicovaginal mucus (<xref ref-type="bibr" rid="B32">Lamont et&#xa0;al., 2011</xref>) and constitute what is considered an optimal vaginal microbiota. The interaction between epithelial cells, cervicovaginal mucus and vaginal microbiota, is important to preserve female reproductive health, by providing a biological barrier that prevents acquisition of STDs (<xref ref-type="bibr" rid="B4">Amabebe and Anumba, 2018</xref>), among other functions. Accordingly, recent studies have demonstrated that infection of the squamous epithelium with HR-HPV was associated with a vaginal microbiota dominated by bacterial vaginosis-associated bacteria (<xref ref-type="bibr" rid="B59">Wei et&#xa0;al., 2021</xref>) and depleted in <italic>Lactobacillus</italic> spp. (<xref ref-type="bibr" rid="B7">Andrade Pessoa Morales et&#xa0;al., 2021</xref>). Moreover, HPV-associated precancerous cervical lesions are also linked with changes in the vaginal microbial composition (<xref ref-type="bibr" rid="B45">Nieves-Ram&#xed;rez et&#xa0;al., 2021</xref>). Very little is known about the vaginal microbial ecology in women with symptomatic cervical ectopy, as a result, the potential interactions between the vaginal microbiota and HPV infection are also currently unknown. This is possibly related to the fact that cervical ectopy is considered a common condition that is often overlooked and understudied.</p>
<p>A general feature of women with symptomatic cervical ectopy is the abnormal discharge of cervicovaginal mucus, which is produced by the everted glandular cells. Cervicovaginal mucus contributes to mucosal immune responses, because it harbors both cytokines and antibodies which are mediators of innate and adaptive immunity. Genital inflammation, often determined by levels of pro-inflammatory cytokines in cervicovaginal mucus, has been extensively studied in women with HPV-associated squamous lesions and tumors (<xref ref-type="bibr" rid="B26">Iwata et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B16">Daniilidis et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B47">Otani et&#xa0;al., 2019</xref>), but lesser so in women with cervical ectopy (<xref ref-type="bibr" rid="B24">Hwang et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B31">Kyongo et&#xa0;al., 2012</xref>). Our group has shown the presence of mucosal IgA anti-HPV antibodies in a proportion of HPV-infected women with cervical ectopy (<xref ref-type="bibr" rid="B41">Monroy et&#xa0;al., 2010</xref>). Although cervical ectopy has been associated with increased levels of pro-inflammatory cytokines, the potential impact of HPV infection in the development of mucosal innate and inflammatory responses in women with symptomatic cervical ectopy has not been fully explored. Moreover, local immune responses are thought to be regulated by the vaginal microbiota (<xref ref-type="bibr" rid="B6">Anahtar et&#xa0;al., 2018</xref>). For instance, the presence of Gram-negative bacteria in a Lactobacilli-enriched genital tract induces a Th17 response and the production of IL-23 (<xref ref-type="bibr" rid="B60">Witkin et&#xa0;al., 2011</xref>). To our knowledge, the impact of the vaginal microbiota in regulating inflammatory and cellular immune responses to HPV in women with cervical ectopy is unknown.</p>
<p>Thus, in the present work we determined the prevalence of HPV genotypes in glandular tissue of women with symptomatic cervical ectopy. Mucosal immune responses were evaluated in cervicovaginal mucus of a selected group of HPV-positive and HPV-negative women, in order to determine their genital inflammation profile (cytokines and chemokines). Next, we also characterized the vaginal microbiota by 16S rRNA sequencing. Finally, associations between HPV infection, vaginal microbial communities and genital inflammation were explored. Our study aims at improving our understanding of the complex interactions of distinct elements within the cervicovaginal environment of women suffering from a common, but scarcely explored condition of the female genital tract.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="s2_1">
<title>Study Population</title>
<p>Female participants were selected from women attending the Colposcopy Clinic of the &#x201c;Cruz-Talonia&#x201d; Foundation, Mexico City, Mexico, a non-profit institution dedicated to preventing cervical cancer, between 2017 and 2018. Women who complained of abnormal vaginal discharge, intermenstrual and/or post-coital bleeding were examined and underwent colposcopy. Women with symptomatic cervical ectopy, defined as irregular patches of columnar epithelium with papillae and clefts, without ectocervical changes after acetic acid application, were invited to participate. Eligible women: older than 20 years of age, sexually active, not pregnant, without a previous history of HPV-associated pathologies, and HIV seronegative, were invited to participate. All participants agreed to participate in this study, signed an informed consent form and answered a validated questionnaire on sociodemographic characteristics, sexual behavior, and medical history. This study was approved by the Local Ethics Committee of the Colposcopy Clinic of the &#x201c;Cruz-Talonia&#x201d; Foundation (Reference Number: 17/02/2017/LRZ/CE05). All women consented to provide a tissue biopsy and a cervicovaginal mucus sample, and to give any remaining unwanted samples to research.</p>
</sec>
<sec id="s2_2">
<title>Sample Collection</title>
<p>All samples were obtained from participants at the seventh day of their menstrual cycle. Colposcopy-directed biopsies were obtained by trained specialists (colposcopists) experienced in the management of lower genital tract pathologies. Biopsies were transported on ice and processed immediately. Tissue biopsies were divided into fragments. One fragment was used for histopathology analysis to corroborate colposcopy-directed diagnosis, tissue biopsies were fixed in formalin and embedded in paraffin, 4 &#xb5;m sections were subsequently deparaffinized, rehydrated and stained with hematoxylin-eosin to be analyzed by experienced pathologists. The second fragment was used for DNA extraction. Cervicovaginal mucus was collected by washing the cervix with 1mL of sterile phosphate-buffered saline. Mucus samples were transported on ice and stored at -20&#xb0;C until used.</p>
</sec>
<sec id="s2_3">
<title>HPV Detection</title>
<p>DNA was isolated from cervical eversion biopsies using the QIAamp DNA kit (Qiagen GmbH, Hilden, Germany), according to the manufacturer&#xb4;s instructions. The presence and genotyping of HPV was performed using the type-specific multiplex genotyping (E7-MPG) assay, which combines multiplex polymerase chain reaction and bead-based Luminex technology (Luminex Co. Austin, TX, USA), at the International Agency for Research on Cancer (IARC), Lyon, France, as described elsewhere (<xref ref-type="bibr" rid="B15">Clifford et&#xa0;al., 2016</xref>). The assay was designed to simultaneously detect 21 mucosal HPV genotypes (HPV6, 11, 16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 70, 73, and 82), which have been classified as HR-HPV genotypes (HPV16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, and 59), putative HR-HPV types (pHR-HPV) (HPV26, 53, 66, 68, 70, 73, 82), and not carcinogenic or low risk types (LR-HPV) (HPV6, and 11) (<xref ref-type="bibr" rid="B25">IARC, 2007</xref>). The assay includes oligonucleotides for the amplification of &#x3b2;&#x2013;globin as a control for DNA integrity. Only samples that tested positive for the amplification of &#x3b2;&#x2013;globin were included in the study.</p>
</sec>
<sec id="s2_4">
<title>16S rRNA Gene Sequencing: DNA Extraction, Amplicon Generation and Sequencing</title>
<p>A total of 63 samples were available for 16S sequencing. DNA was extracted from cervicovaginal mucus using the PureLink&#x2122; Microbiome DNA Purification Kit (ThermoFisher Scientific, Waltham, MA, USA). Libraries were prepared as described in the MiSeq 16S rRNA gene Amplicon Sequencing protocol (Illumina, San Diego, CA, USA) with modifications. Primers for the V4 region were 515F (5&#xb4;- GTGCCAGCMGCCGCGGTAA-3&#xb4;and 806R (5&#xb4;- GGACTACHVGGGTWTCTAAT-3&#xb4;) (<xref ref-type="bibr" rid="B13">Caporaso et&#xa0;al., 2011</xref>). Briefly, the V4 region was amplified by PCR (triplicate reactions of 25 &#xb5;L per sample). Triplicate PCR amplifications were confirmed in a 2% agarose gel electrophoresis. Triplicate PCR reactions were pooled per sample and purified using AgenCourt AMPure XP beads (Beckman Dickson, Atlanta, GA, USA). Dual indices were attached using the Nextera XT Index Kit (Illumina). Next, indexed PCRs were purified twice using AgenCourt AMPure XP beads, and quantified using Qubit fluorometer (ThermoFisher Scientific, Waltham, MA, USA). Equimolar concentrations of indexed PCR libraries were pooled together at 4nM. Negative controls were also included to control for exogenous contaminants. Quality control and molarity of the final library was assessed on the Agilent&#xae; High Sensitivity DNA Kit/2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). Sequencing was performed on an Illumina MiSeq&#x2122; platform (Illumina) according to the manufacturer&#x2019;s specifications. The final library was sequenced at 8pM concentration with 25% PhiX to generate paired-end reads of 500 base-length (v2 2x250 cycles). A total of 13,115,362 sequences were obtained from a single run.</p>
</sec>
<sec id="s2_5">
<title>16S Data Analysis</title>
<p>Raw Illumina MiSeq sequences were processed using the Qiime2 (Quantitative Insights into Microbial Ecology, version 2019.4) (<xref ref-type="bibr" rid="B9">Bolyen et&#xa0;al., 2019</xref>). First, sequences that might arise from sequencing errors or low abundance were filtered out, chimeras were filtered out as well, then any PhiX reads were removed and construction of a feature table (selection of amplicon variant sequences (ASVs)) were performed using dada2 (<xref ref-type="bibr" rid="B12">Callahan et&#xa0;al., 2016</xref>), options were: &#x2013;p-trim-left-f 19, &#x2013;p-trim-left-r 20, &#x2013;p-trunc-len-r 180 and &#x2013;p-trunc-len-f 185. This method allows for distinct reads with as few as 1 nucleotide difference to be clustered into separate ASVs. A total number of 467 ASVs were identified, and a median of 177,565 (71,692-453,044) sequences were retained after dada2. Taxonomy was assigned to ASVs using the q2&#x2010;feature&#x2010;classifier using the classify&#x2010;sklearn na&#xef;ve Bayes taxonomy classifier against the GreenGenes (13_8) 99% reference sequences and trained to the region of the target sequences, in our case for the V4 region (515F-806R) (<xref ref-type="bibr" rid="B38">McDonald et&#xa0;al., 2012</xref>). Phylogenetic tree was constructed with fasttree2 (<italic>via</italic> q2&#x2010;phylogeny) after aligning all ASVs with mafft (<xref ref-type="bibr" rid="B28">Katoh et&#xa0;al., 2002</xref>; <xref ref-type="bibr" rid="B49">Price et&#xa0;al., 2010</xref>). Artifacts generated with Qiime2 were imported into Phyloseq for further manipulation, and graph visualization (R V.3.6.2, The R Foundation for Statistical Computing (<xref ref-type="bibr" rid="B39">McMurdie and Holmes, 2013</xref>). Before continuing with further manipulation, rarefaction was performed at a sampling depth of 64,522 sequences/sample (corresponding to 90% of the sample with the least sequences). Principal Coordinate analyses (PCoA) were used to visualize microbial communities and permutational multivariate analysis of variance (PERMANOVA, Adonis function) was used to assess differences in beta diversity (default parameters: 999 permutations). We used 2 metrics for beta diversity (Bray-Curtis dissimilarity index and weighted UniFrac) (<xref ref-type="bibr" rid="B36">Lozupone et&#xa0;al., 2011</xref>). Linear discriminant analysis (LDA) effect size (LEfSe) was used to detect significant differentially abundances features (taxa) between groups (default parameters were used, LDA threshold= 3, alpha value for factorial Kruskal-Wallis test 0.05, alpha value for pairwise Wilcoxon test 0.05 (<xref ref-type="bibr" rid="B53">Segata et&#xa0;al., 2011</xref>). Spearman correlations were performed using rcorr() function and visualized using corrplot() function.</p>
</sec>
<sec id="s2_6">
<title>Cytokines and Chemokines Determination</title>
<p>The following cytokines and chemokines were evaluated using bead-based multiplex immune-assays by flow cytometry according to the protocols of the manufacturer: Interleukin-1&#x3b2; (IL-1&#x3b2;), and Interleukin-12p70 (IL-12p70) (BD Cytometric Bead Array (CBA) Human Inflammatory Cytokines Kit. BD Biosciences, CA, USA); Interleukin-2 (IL-2), Interleukin-4 (IL-4), Interleukin-5 (IL-5), Interleukin-6 (IL-6), Interleukin-9 (IL-9), Interleukin-10 (IL-10), Interleukin-13 (IL-13), Interleukin-17A (IL-17A), Interleukin-17F (IL-17F), Interleukin-21 (IL-21), Interleukin-22 (IL-22), Interferon-&#x3b3; (IFN-&#x3b3;), and Tumor Necrosis Factor- &#x3b1; (TNF-&#x3b1;) (LEGENDplex&#x2122; Human Th Cytokine Panel (13-plex), BioLegend, CA, USA). MCP-1 (CCL2), RANTES (CCL5), IP-10 (CXCL10), Eotaxin (CCL11), TARC (CCL17), MIP-1&#x3b1; (CCL3), MIP-1&#x3b2; (CCL4), MIG (CXCL9), MIP-3&#x3b1; (CCL20), ENA-78 (CXCL5), GRO&#x3b1; (CXCL1), I-TAC (CXCL11) and IL-8 (CXCL8) (LEGENDplex&#x2122; Human Proinflammatory Chemokine Panel. BioLegend, CA, USA). Briefly, cervicovaginal mucus samples were thawed and centrifuged at 3000 x <italic>g</italic> for 15 min at 4&#xb0;C, supernatants were recovered and the total protein content was calculated using the DC Protein Assay (Bio-Rad Inc., CA, USA). Samples were analyzed in a FACS ARIA III flow cytometer (BD Biosciences, CA, USA). Cytokines and chemokines concentrations were calculated using the LEGENDplex&#x2122; Data Analysis Software v 7.0 (BioLegend, CA, USA). Cytokine and chemokine expression levels were normalized to total protein concentration. The lower limit of detection (LLOD) value, provided by the manufacturer, was assigned to those samples with cytokine or chemokine concentrations equal or below the reported limit. To control for inter-individual variation due to differences in the amount of protein contained in the samples, the final concentration of analytes was reported as pg/mg of total protein.</p>
</sec>
<sec id="s2_7">
<title>Statistical Analysis</title>
<p>Data was analyzed using the GraphPad Prism 7 software (GraphPad, USA). Differences in demographic and clinical factors between groups were analyzed using the Mann-Whitney test for continuous variables and Chi-square test for categorical values. Comparisons between two groups were performed using the Wilcoxon Rank Sum Test. Comparisons between three or more groups were performed using Kruskal-Wallis, <italic>p</italic> values were adjusted for multiple comparisons using Dunn&#x2019;s multiple comparisons test. When using this test, adjusted-<italic>p</italic> values were reported. Two-tailed <italic>p</italic> values (and adjusted <italic>p</italic> values for KW test) lesser than 0.05 were considered statistically significant.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<title>Results</title>
<sec id="s3_1">
<title>Prevalence of HPV in Women With Symptomatic Cervical Ectopy</title>
<p>A total of 156 samples from women with symptomatic cervical ectopy were available for HPV testing. A total of 19 HPV types were detected: HPV6, 11, 16, 18, 26, 31, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, and 82. The prevalence and distribution of HPV genotypes is shown in <xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>. Overall, 85 (54.48%) women tested positive for HPV infection. Both single and multiple HPV infections were detected. A total of 59 women had a single HPV type with 93.22% testing positive for HR-HPV or pHR-HPV genotypes, and 6.77% for LR-HPV genotypes. HPV infection with multiple types accounted for 30.58% of all HPV types (26/85), and all women in this category tested positive for HR-HPV types. The relative frequency of HPV genotypes detected as single and multiple infections is shown in <xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>. HPV16 was the most predominant, detected in 29.41% of the positive samples, followed by HPV31 (21.17%), and HPV18 (15.29%) (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref> and <xref ref-type="supplementary-material" rid="ST1"><bold>Supplementary Table&#xa0;1</bold></xref>). Women with HPV infection (HPV-pos) were younger than HPV negative (HPV-neg) women (<italic>p</italic>=0.0459, <xref ref-type="table" rid="T2"><bold>Table&#xa0;2</bold></xref> and <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure&#xa0;1</bold></xref>). Of interest, although all women reported to be sexually active and had had children, a total of 140 women reported never having had a Papanicolaou (PAP) smear test, with no difference being observed between HPV-neg and HPV-pos women. A total of 17 women were vaccinated against HPV with the 4-valent HPV vaccine targeting HPV-6, -11, -16, and -18 types (Gardasil<sup>&#xae;</sup>, MSD), again no difference was observed between HPV-neg (8/17) and HPV-pos (9/17) women. It is worth mentioning that 3 of the HPV-vaccinated women were positive for HPV-16, whereas 6 HPV-vaccinated women were infected with HPV types not included in the vaccine. Taken the latter in consideration, and the prevalence of HR-HPV types in our cohort, these women might be at increased risk for HPV-associated pathologies.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>HPV prevalence and genotype distribution in women with symptomatic cervical ectopy.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">HPV</th>
<th valign="top" align="center">Total (N=156)</th>
</tr>
<tr>
<th valign="top" align="left">Type</th>
<th valign="top" align="center">N (%)</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Negative</td>
<td valign="top" align="center">71 (45.51)</td>
</tr>
<tr>
<td valign="top" align="left">Any HPV</td>
<td valign="top" align="center">85 (54.48)</td>
</tr>
<tr>
<td valign="top" align="left">High Risk</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">16</td>
<td valign="top" align="center">25 (16.03)</td>
</tr>
<tr>
<td valign="top" align="left">18</td>
<td valign="top" align="center">13 (8.33)</td>
</tr>
<tr>
<td valign="top" align="left">31</td>
<td valign="top" align="center">18 (11.53)</td>
</tr>
<tr>
<td valign="top" align="left">35</td>
<td valign="top" align="center">1 (0.64)</td>
</tr>
<tr>
<td valign="top" align="left">39</td>
<td valign="top" align="center">10 (6.41)</td>
</tr>
<tr>
<td valign="top" align="left">45</td>
<td valign="top" align="center">4 (2.56)</td>
</tr>
<tr>
<td valign="top" align="left">51</td>
<td valign="top" align="center">10 (6.41)</td>
</tr>
<tr>
<td valign="top" align="left">52</td>
<td valign="top" align="center">3 (1.92)</td>
</tr>
<tr>
<td valign="top" align="left">56</td>
<td valign="top" align="center">3 (1.92)</td>
</tr>
<tr>
<td valign="top" align="left">58</td>
<td valign="top" align="center">8 (5.13)</td>
</tr>
<tr>
<td valign="top" align="left">59</td>
<td valign="top" align="center">7 (4.49)</td>
</tr>
<tr>
<td valign="top" align="left">Putative High Risk</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">26</td>
<td valign="top" align="center">8 (5.13)</td>
</tr>
<tr>
<td valign="top" align="left">53</td>
<td valign="top" align="center">9 (5.77)</td>
</tr>
<tr>
<td valign="top" align="left">66</td>
<td valign="top" align="center">2 (2.56)</td>
</tr>
<tr>
<td valign="top" align="left">68</td>
<td valign="top" align="center">2 (1.28)</td>
</tr>
<tr>
<td valign="top" align="left">73</td>
<td valign="top" align="center">1 (0.64)</td>
</tr>
<tr>
<td valign="top" align="left">82</td>
<td valign="top" align="center">1 (0.64)</td>
</tr>
<tr>
<td valign="top" align="left">Low Risk</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">6</td>
<td valign="top" align="center">5 (3.20)</td>
</tr>
<tr>
<td valign="top" align="left">11</td>
<td valign="top" align="center">8 (5.13)</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>HPV, Human Papillomavirus.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Histogram showing the distribution of human papillomavirus (HPV) genotypes detected as single and multiple infection, relative to the total number HPV-positive samples (N=85).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-12-884272-g001.tif"/>
</fig>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Clinical, demographic and risk factors associated with HPV infection.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Characteristic</th>
<th valign="top" align="center">HPV-pos (N=85)</th>
<th valign="top" align="center">HPV-neg (N=71)</th>
<th valign="top" align="center"><italic>p</italic> value </th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Age (years)</td>
<td valign="top" align="center">30.77</td>
<td valign="top" align="center">32.75</td>
<td valign="top" align="center">0.0459<sup>a*</sup>
</td>
</tr>
<tr>
<td valign="top" align="left">Parity</td>
<td valign="top" align="center">2.14</td>
<td valign="top" align="center">2.07</td>
<td valign="top" align="center">0.8949<xref ref-type="table-fn" rid="fnT2_1"><sup>a</sup></xref>
</td>
</tr>
<tr>
<td valign="top" align="left">Menarche (years)</td>
<td valign="top" align="center">12.64</td>
<td valign="top" align="center">12.31</td>
<td valign="top" align="center">0.1612<xref ref-type="table-fn" rid="fnT2_1"><sup>a</sup></xref>
</td>
</tr>
<tr>
<td valign="top" align="left">Age of sexual debut (years)</td>
<td valign="top" align="center">17.72</td>
<td valign="top" align="center">18.08</td>
<td valign="top" align="center">0.6256<xref ref-type="table-fn" rid="fnT2_1"><sup>a</sup></xref>
</td>
</tr>
<tr>
<td valign="top" align="left">Number of sexual partners</td>
<td valign="top" align="center">2.21</td>
<td valign="top" align="center">1.82</td>
<td valign="top" align="center">0.0550<xref ref-type="table-fn" rid="fnT2_1"><sup>a</sup></xref>
</td>
</tr>
<tr>
<td valign="top" align="left">Smoking status (%)</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.7196<xref ref-type="table-fn" rid="fnT2_2"><sup>b</sup></xref>
</td>
</tr>
<tr>
<td valign="top" align="left">Current smoker<break/>Non-smoker</td>
<td valign="top" align="center">16 (18.8)<break/>69 (81.2)</td>
<td valign="top" align="center">15 (21.1)<break/>56 (78.9)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">Contraception (%)</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.4949<xref ref-type="table-fn" rid="fnT2_2"><sup>b</sup></xref>
</td>
</tr>
<tr>
<td valign="top" align="left">None<break/>Physical<xref ref-type="table-fn" rid="fnT2_3"><sup>c</sup></xref>
<break/>Hormonal<xref ref-type="table-fn" rid="fnT2_4"><sup>d</sup></xref>
<break/>Chirurgical<xref ref-type="table-fn" rid="fnT2_5"><sup>e</sup></xref>
</td>
<td valign="top" align="center">24 (28.2)<break/>38 (44.7)<break/>9 (10.6)<break/>14 (16.5)</td>
<td valign="top" align="center">26 (36.6)<break/>27 (38.0)<break/>10 (14.1)<break/>8 (11.3)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">HPV vaccine (%)</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.8921<xref ref-type="table-fn" rid="fnT2_2"><sup>b</sup></xref>
</td>
</tr>
<tr>
<td valign="top" align="left">Yes<break/>No</td>
<td valign="top" align="center">9 (10.6)<break/>76 (89.4)</td>
<td valign="top" align="center">8 (11.7)<break/>63 (88.7)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">Prior pap smear (%)</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.7036<xref ref-type="table-fn" rid="fnT2_2"><sup>b</sup></xref>
</td>
</tr>
<tr>
<td valign="top" align="left">Yes<break/>No</td>
<td valign="top" align="center">8 (9.4)<break/>77 (90.6)</td>
<td valign="top" align="center">8 (11.7)<break/>63 (88.7)</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">Chronic diseases and conditions (%)</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.4781<xref ref-type="table-fn" rid="fnT2_2"><sup>b</sup></xref>
</td>
</tr>
<tr>
<td valign="top" align="left">Arthritis<break/>Asthma<break/>Diabetes<break/>Hypertension<break/>Other</td>
<td valign="top" align="center">1 (1.2)<break/>2 (2.4)<break/>1 (1.2)<break/>1 (1.2)<break/>3 (3.5)</td>
<td valign="top" align="center">0 (0.0)<break/>1 (1.4)<break/>3 (4.2)<break/>3 (4.2)<break/>1 (1.4)</td>
<td valign="top" align="center"/>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Values are given as mean and number (%).</p>
</fn>
<fn>
<p>HPV, Human Papillomavirus.</p>
</fn>
<fn>
<p>HPV-pos, HPV positive.</p>
</fn>
<fn>
<p>HPV-neg, HPV negative.</p>
</fn>
<fn id="fnT2_1">
<label>a</label>
<p>By Mann-Whitney test.</p>
</fn>
<fn id="fnT2_2">
<label>b</label>
<p>By Chi-square test.</p>
</fn>
<fn id="fnT2_3">
<label>c</label>
<p>Includes condoms, and intrauterine device (IUD).</p>
</fn>
<fn id="fnT2_4">
<label>d</label>
<p>Includes implants, injections, patches, and pills.</p>
</fn>
<fn id="fnT2_5">
<label>e</label>
<p>Tubal ligation surgery.</p>
</fn>
<fn>
<p>*p &lt; 0.05 was considered significant.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Subsequent analyses were performed in a sub-cohort of 63 women including 35 HPV-pos and 28 HPV-neg women. This sub-cohort was selected based on sample availability (characteristics and HPV prevalence and genotype distribution are summarized in <xref ref-type="supplementary-material" rid="ST2"><bold>Supplementary Tables&#xa0;2</bold></xref>, <xref ref-type="supplementary-material" rid="ST3"><bold>3</bold></xref> respectively).</p>
</sec>
<sec id="s3_2">
<title>Symptomatic Cervical Ectopy Is Associated With Generalized Inflammation While HPV Infection Is Associated With High Levels of IL-21 and CXCL9</title>
<p>Next, we explored the impact of HPV infection on mucosal inflammation as it is well documented that persistent HPV infection is associated with increased levels of local and systemic inflammatory cytokines (<xref ref-type="bibr" rid="B29">Kemp et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B52">Scott et&#xa0;al., 2013</xref>). Sixteen pro-inflammatory cytokines and 12 chemokines were determined in cervicovaginal mucus. Levels of IL-8, IL-1&#x3b2;, IL-22, CCL17, and CCL3 were comparable between HPV-pos and HPV-neg women (<xref ref-type="table" rid="T3"><bold>Table&#xa0;3</bold></xref>). These 3 cytokines and 2 chemokines were detected in more than 90% of women suggesting that most women with symptomatic cervical ectopy suffered from an inflammatory response independent of the presence of HPV. As shown in <xref ref-type="table" rid="T3"><bold>Table&#xa0;3</bold></xref>, levels of IL-21 (27.74 pg/mg in HPV-pos <italic>vs</italic> 4.89 pg/mg in HPV-neg women) and MIG (56.89 pg/mg in HPV-pos vs 15.34 pg/mg in HPV-neg women) were significantly higher in HPV-pos women (<italic>p</italic>=0.0002 and <italic>p</italic>=0.013, respectively), which is suggestive of an anti-viral cellular immune response.</p>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Cytokine and chemokine concentrations in cervicovaginal mucus from women with symptomatic cervical ectopy.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" rowspan="2" align="left">Analyte</th>
<th valign="top" colspan="2" align="center">HPV-pos</th>
<th valign="top" colspan="2" align="center">HPV-neg</th>
<th valign="top" rowspan="2" align="center"><italic>p</italic> value<xref ref-type="table-fn" rid="fnT3_1"><sup>a</sup></xref>
</th>
</tr>
<tr>
<th valign="top" colspan="2" align="center">N=35 (pg/mg)</th>
<th valign="top" colspan="2" align="center">N=28 (pg/mg)</th>
</tr>
<tr>
<th valign="top" align="left"/>
<th valign="top" align="center">Median</th>
<th valign="top" align="center">IQR</th>
<th valign="top" align="center">Median</th>
<th valign="top" align="center">IQR</th>
<th valign="top" align="center"/>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Cytokines</td>
<td valign="top" colspan="2" align="center"/>
<td valign="top" colspan="2" align="center"/>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-1&#x3b2;</td>
<td valign="top" align="center">693.10</td>
<td valign="top" align="center">(213.5-1574)</td>
<td valign="top" align="center">485.50</td>
<td valign="top" align="center">(313.9-1416)</td>
<td valign="top" align="center">0.833</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-2</td>
<td valign="top" align="center">4.03</td>
<td valign="top" align="center">(2.3-9.6)</td>
<td valign="top" align="center">1.79</td>
<td valign="top" align="center">(0.7-3.9)</td>
<td valign="top" align="center">0.114</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-4</td>
<td valign="top" align="center">5.44</td>
<td valign="top" align="center">(1.8-16.2)</td>
<td valign="top" align="center">1.95</td>
<td valign="top" align="center">(0.9-3.3)</td>
<td valign="top" align="center">0.145</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-6</td>
<td valign="top" align="center">13.08</td>
<td valign="top" align="center">(5.9-66.4)</td>
<td valign="top" align="center">13.89</td>
<td valign="top" align="center">(5.5-87)</td>
<td valign="top" align="center">0.922</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-8</td>
<td valign="top" align="center">6435.00</td>
<td valign="top" align="center">(2290-13743)</td>
<td valign="top" align="center">4591.00</td>
<td valign="top" align="center">(2957-8609)</td>
<td valign="top" align="center">0.226</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-9</td>
<td valign="top" align="center">3.39</td>
<td valign="top" align="center">(2.3-7.4)</td>
<td valign="top" align="center">2.19</td>
<td valign="top" align="center">(1.5-3.3)</td>
<td valign="top" align="center">0.080</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-12p</td>
<td valign="top" align="center">21.55</td>
<td valign="top" align="center">(16.9-26.2)</td>
<td valign="top" align="center">45.17</td>
<td valign="top" align="center">(35.2-83.6)</td>
<td valign="top" align="center">0.143</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-13</td>
<td valign="top" align="center">10.47</td>
<td valign="top" align="center">(5.2-25.7)</td>
<td valign="top" align="center">9.80</td>
<td valign="top" align="center">(3.6-15.9)</td>
<td valign="top" align="center">0.857</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-17B</td>
<td valign="top" align="center">4.24</td>
<td valign="top" align="center">(3.9-8.1)</td>
<td valign="top" align="center">1.60</td>
<td valign="top" align="center">(0.5-4.1)</td>
<td valign="top" align="center">0.143</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-21</td>
<td valign="top" align="center">16.24</td>
<td valign="top" align="center">(3.5-30.1)</td>
<td valign="top" align="center">2.02</td>
<td valign="top" align="center">(1.7-5.6)</td>
<td valign="top" align="center">0.0002*</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IL-22</td>
<td valign="top" align="center">983.60</td>
<td valign="top" align="center">(463.7-2402)</td>
<td valign="top" align="center">973.50</td>
<td valign="top" align="center">(348.6-1710)</td>
<td valign="top" align="center">0.518</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;TNF-&#x3b1;</td>
<td valign="top" align="center">26.20</td>
<td valign="top" align="center">(13.1-44)</td>
<td valign="top" align="center">17.45</td>
<td valign="top" align="center">(6.8-34.9)</td>
<td valign="top" align="center">0.413</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;IFN-&#x3b3;</td>
<td valign="top" align="center">20.18</td>
<td valign="top" align="center">(3.9-50.7)</td>
<td valign="top" align="center">8.64</td>
<td valign="top" align="center">(3.3-16.7)</td>
<td valign="top" align="center">0.445</td>
</tr>
<tr>
<td valign="top" align="left">Chemokines</td>
<td valign="top" colspan="2" align="center"/>
<td valign="top" colspan="2" align="center"/>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CXCL10</td>
<td valign="top" align="center">18.06</td>
<td valign="top" align="center">(10.6-35.2)</td>
<td valign="top" align="center">22.61</td>
<td valign="top" align="center">(7.4-41.5)</td>
<td valign="top" align="center">0.682</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CCL11</td>
<td valign="top" align="center">8.54</td>
<td valign="top" align="center">(6.2-11.8)</td>
<td valign="top" align="center">10.73</td>
<td valign="top" align="center">(4-26.5)</td>
<td valign="top" align="center">0.576</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CCL17</td>
<td valign="top" align="center">9.97</td>
<td valign="top" align="center">(7.4-15.5)</td>
<td valign="top" align="center">9.68</td>
<td valign="top" align="center">(3.9-24.1)</td>
<td valign="top" align="center">0.715</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CCL2</td>
<td valign="top" align="center">26.47</td>
<td valign="top" align="center">(9.7-85.6)</td>
<td valign="top" align="center">44.90</td>
<td valign="top" align="center">(20.4-85.1)</td>
<td valign="top" align="center">0.378</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CCL5</td>
<td valign="top" align="center">19.82</td>
<td valign="top" align="center">(11.5-40.4)</td>
<td valign="top" align="center">18.71</td>
<td valign="top" align="center">(10.2-61.4)</td>
<td valign="top" align="center">0.865</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CCL3</td>
<td valign="top" align="center">19.77</td>
<td valign="top" align="center">(10.2-28.3)</td>
<td valign="top" align="center">16.83</td>
<td valign="top" align="center">(7.7-39.7)</td>
<td valign="top" align="center">0.608</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CXCL9</td>
<td valign="top" align="center">22.26</td>
<td valign="top" align="center">(8.7-76.9)</td>
<td valign="top" align="center">11.90</td>
<td valign="top" align="center">(4.2-21.1)</td>
<td valign="top" align="center">0.013*</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CXCL5</td>
<td valign="top" align="center">53.24</td>
<td valign="top" align="center">(26.3-162.5)</td>
<td valign="top" align="center">26.46</td>
<td valign="top" align="center">(18-79.5)</td>
<td valign="top" align="center">0.069</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CCL20</td>
<td valign="top" align="center">20.84</td>
<td valign="top" align="center">(10.4-37.7)</td>
<td valign="top" align="center">21.49</td>
<td valign="top" align="center">(6.4-49.2)</td>
<td valign="top" align="center">0.814</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CXCL1</td>
<td valign="top" align="center">30.06</td>
<td valign="top" align="center">(18.5-77.1)</td>
<td valign="top" align="center">43.35</td>
<td valign="top" align="center">(12.6-75.4)</td>
<td valign="top" align="center">0.926</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CXCL11</td>
<td valign="top" align="center">7.06</td>
<td valign="top" align="center">(3.8-16.9)</td>
<td valign="top" align="center">11.22</td>
<td valign="top" align="center">(3.1-19.2)</td>
<td valign="top" align="center">0.771</td>
</tr>
<tr>
<td valign="top" align="left">&#x2003;CCL4</td>
<td valign="top" align="center">22.83</td>
<td valign="top" align="center">(12.6-36.8)</td>
<td valign="top" align="center">23.06</td>
<td valign="top" align="center">(6.4-52.2)</td>
<td valign="top" align="center">0.684</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Values are given as median (IQR 25-75).</p>
</fn>
<fn>
<p>HPV, Human Papillomavirus.</p>
</fn>
<fn>
<p>HPV-pos, HPV positive samples.</p>
</fn>
<fn>
<p>HPV-neg, HPV negative samples.</p>
</fn>
<fn id="fnT3_1">
<label>a</label>
<p>By Mann-Whitney test.</p>
</fn>
<fn>
<p>*p &lt; 0.05 was considered significant.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s3_3">
<title>Changes in Vaginal Microbiota Composition Are Associated With HPV Infection</title>
<p>Our results suggest that women with symptomatic cervical ectopy suffer from a state of generalized mucosal inflammation which is not entirely explained by the presence of HPV. Strong evidence indicates that inflammatory responses in the female genital tract are modulated by the local microbiota (<xref ref-type="bibr" rid="B5">Anahtar et&#xa0;al., 2015</xref>). Thus, we investigated the microbiota composition in cervicovaginal mucus of women with symptomatic cervical ectopy. Both richness and shannon diversity were increased in HPV-pos compared to HPV-neg women (<italic>p</italic>&lt;0.001, <xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2A</bold></xref>). Principal coordinates analysis (PCoA) was used to visualize the clustering of microbial communities stratified by HPV infection. HPV infection had an impact, albeit small, on the vaginal microbiota (explaining around 5% of variation, PERMANOVA <italic>p</italic>=0.019 for weighted UniFrac and <italic>p</italic>=0.013 for Bray-Curtis dissimilarity, <xref ref-type="fig" rid="f2"><bold>Figures&#xa0;2B, C</bold></xref>, respectively). Given that HPV infection had a small albeit measurable impact on the vaginal microbiota, we next stratified HPV-positive women into two groups, those with single HPV infection (representing the leading 3 most prevalent HR-HPV types, HPV16, HPV18 and HPV31) and those with multiple HPV infection. Women with single HPV infection had significantly increased microbial richness and diversity (Richness, adjusted p=0.0021 and Shannon, adjusted p=0.0021. <xref ref-type="supplementary-material" rid="SM2"><bold>Supplementary Figure&#xa0;2A</bold></xref>) compared to HPV-negative women. Women with multiple HPV infection had increased richness (Richness, adjusted p=0.041) but not diversity (Shannon, adjusted p&gt;0.05. <xref ref-type="supplementary-material" rid="SM2"><bold>Supplementary Figure&#xa0;2A</bold></xref>). PCoA was used to visualize the clustering of microbial communities of the 3 groups: HPV-negative, single HPV infection and multiple HPV infection (<xref ref-type="supplementary-material" rid="SM2"><bold>Supplementary Figures&#xa0;2B, C</bold></xref>). Single HPV infection had a significant impact on the vaginal microbiota (weighted UniFrac PERMANOVA <italic>p</italic>=0.014 and Bray-Curtis <italic>p</italic>=0.003). Multiple HPV infection had no impact on the vaginal microbiota (weighted UniFrac PERMANOVA <italic>p</italic>=0.149 and Bray-Curtis <italic>p</italic>=0.122). Results seem to indicate that differences in vaginal microbial composition between HPV-negative and HPV-positive women might be driven by single HPV infection, rather than multiple HPV infection. Therefore, next we restricted our analysis to the 3 predominant HR-HPV types found in our study. The prevalence of HPV16, HPV18 and HPV31 single infection was 8/35 (22.9%), 5/35 (14.3%) and 7/35 (20%), respectively.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>The vaginal microbiota (cervicovaginal mucus) was characterized in women with symptomatic cervical ectopy with or without HPV infection. <bold>(A)</bold> Two alpha diversity indexes were calculated, richness (observed species) and shannon. Boxplots depicting the median and interquartile range stratified by HPV infection. Groups were compared using the Wilcoxon Rank Sum test. Both weighted UniFrac <bold>(B)</bold> and Bray-Curtis dissimilarity <bold>(C)</bold> were used in conjunction with principal coordinate analysis (PCoA) to visualize microbial communities stratified by HPV infection. Differences in beta diversity were assessed using permutational multivariate analysis of variance (PERMANOVA). Each point was colored according to the HPV type for HPV-pos women and order by the number of HPV types found.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-12-884272-g002.tif"/>
</fig>
<p>HPV16, HPV31, and HPV18 were the most prevalent genotypes in our cohort. Thus, we further analyzed the microbial composition including only women with single HPV16, single HPV31, and single HPV18 infection. Bacterial richness and diversity were increased in HPV16-positive (Richness, adjusted <italic>p</italic>=0.059; Shannon, adjusted <italic>p</italic>=0.030) and HPV31-positive (Richness, adjusted <italic>p</italic>=0.049; Shannon, adjusted <italic>p</italic>=0.047) compared to HPV-neg women (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3A</bold></xref>). We found that microbiota diversity was not significantly affected by the presence of HPV18 (<italic>p</italic>&gt;0.1 for both alpha diversity indices) (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3A</bold></xref>). Beta diversity analyses showed bacterial community differences between HPV-neg women and HPV31-positive women (explaining 11.9% of variation, PERMANOVA <italic>p</italic>=0.004 for weighted UniFrac and 7% of variation, PERMANOVA <italic>p</italic>=0.012 for Bray-Curtis dissimilarity) (<xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3B, C</bold></xref>), and HPV18-positive women (PERMANOVA <italic>p</italic>=0.044 for weighted UniFrac, and PERMANOVA <italic>p</italic>=0.02 for Bray-Curtis dissimilarity) (<xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3B, C</bold></xref>). However, no differences in bacterial communities between HPV16-positive and HPV-neg women were observed (weighted UniFrac PERMANOVA <italic>p</italic>=0.114 and Bray-Curtis PERMANOVA <italic>p</italic>=0.061) (<xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3B, C</bold></xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>The impact of single HPV16, HPV31 and HPV18 on the vaginal microbiota diversity and community membership. <bold>(A)</bold> Two alpha diversity indexes were calculated, richness (observed species) and shannon. Boxplots depicting the median and interquartile range stratified by single HPV infection. Groups were compared using the Kruskal-Wallis test, <italic>p</italic> values were adjusted for multiple comparisons. Both weighted UniFrac <bold>(B)</bold> and Bray-Curtis dissimilarity <bold>(C)</bold> were used in conjunction with principal coordinate analysis (PCoA) to visualize microbial communities stratified by single HPV infection. Differences in beta diversity were assessed using permutational multivariate analysis of variance (PERMANOVA). Each point was colored according to the single HPV type. Confidence ellipses are shown for each group.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-12-884272-g003.tif"/>
</fig>
<p>The taxonomic profiling of the vaginal microbiota is shown in <xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4</bold></xref>. The top 30 genera stratified by HPV status, ordered by the relative abundance of <italic>Lactobacillus</italic> in each woman is shown in <xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4A</bold></xref>. The top three genera were <italic>Lactobacillus</italic>, <italic>Gardnerella</italic> and <italic>Atopobium</italic>. HPV31-positive women had increased <italic>Sneathia</italic>, <italic>Shuttleworthia</italic>, <italic>Pseudomonas</italic>, <italic>Prevotella</italic>, <italic>Megasphaera</italic>, <italic>Clostridium</italic>, <italic>Atopobium</italic>, and decreased <italic>Lactobacillus</italic> spp. and <italic>Gardnerella</italic> compared to HPV-neg women, whereas HPV-16 women showed increased <italic>Sneathia, Pseudomonas, Megasphaera, Atopobium</italic>, and decreased <italic>Lactobacillus</italic> spp. Although changes in microbiota diversity were not significant for HPV18-positive women, we observed a tendency for an increase in <italic>Sneathia, Prevotella</italic>, and <italic>Gardnerella</italic>, along with a decrease in <italic>Lactobacillus</italic> spp. (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4A</bold></xref>). LEfSe was used to detect significant differentially abundances features (taxa) between HPV-pos and HPV-neg women (LDA threshold of 3 or higher, other default parameters were used). A total of 17 features were identified. In particular, 3 genera were found to differentiate HPV-pos women (<italic>Atopobium</italic>, <italic>Megasphaera, Prevotella</italic>) (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4B</bold></xref>). In contrast, <italic>Lactobacillus</italic> was identified as a discriminant feature in HPV-neg women (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4B</bold></xref>). No discriminant features were found when analysing single vs multiple HPV infections or single HPV16, HPV31 and HPV18 infections.</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Taxonomic profile of the vaginal microbiota at genus level in women with symptomatic cervical ectopy. <bold>(A)</bold> Taxa barplots showing the taxonomic composition at genus level. The top 30 genera are represented, these account for 99.5% of all genera present. Genera are listed from the less abundant to the more abundant within each phylum, and ordered in ascending order (Tenericutes, Bacteroidetes, Fusobacteria, Proteobacteria, Actinobacteria, and Firmicutes). Color bars reflect the phyla they belong to and the intensity of the color their relative abundance, with light colors depicting the least abundant, and the darker colors the most abundant. <bold>(B)</bold> Histogram of the linear discriminant analysis scores reveals that 17 differential features were identified. <italic>Atopobium</italic>, <italic>Megasphaera</italic>, <italic>and Prevotella</italic> were found enriched in HPV-pos women while <italic>Lactobacillus</italic> was enriched in HPV-neg. Linear discriminant analysis Effect Size (LEfSe) analysis was performed with LDA values of 3.0 or higher.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-12-884272-g004.tif"/>
</fig>
</sec>
<sec id="s3_4">
<title>Associations Between Pro-Inflammatory Cytokines and Vaginal Microbiota in the Context of HPV Infection</title>
<p>Finally, we explored associations between genital inflammation and vaginal microbiota (richness, diversity and top 30 genera) in the context of HPV infection. We restricted our analysis to the cytokines and chemokines who were determined in 85% or more samples (IL-1&#x3b2;, IL-8, IL-22, CCL11, CCL17, CCL2, CCL5, CCL3, CXCL5, and CXCL1) (<xref ref-type="supplementary-material" rid="ST4"><bold>Supplementary Table&#xa0;4</bold></xref>). When the entire cohort was considered, irrespective of HPV status, we found strong inverse correlations (rho &gt;-0.60, <italic>p</italic>&lt;0.001) between <italic>Lactobacillus</italic> spp., bacterial-vaginosis associated anaerobes and vaginal richness and diversity, as would be expected (<xref ref-type="supplementary-material" rid="SM3"><bold>Supplementary Figure&#xa0;3</bold></xref>). At the same time, positive correlations (rho&gt;0.40, <italic>p</italic>&lt;0.001) were found between alpha diversity and most anaerobes. No associations were found between inflammation-associated cytokines and chemokines and the vaginal microbiota.</p>
<p>When stratifying for HPV infection, distinct correlations were found. In HPV-pos women, we found positive correlations between alpha diversity and IL-1&#x3b2; (rho= 0.56, <italic>p</italic>=0.0004) and IL-22 (rho= 0.60, p=0.0001) as shown in <xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>. Negative correlations (rho&gt;-0.60, <italic>p</italic>&lt;0.001) were found between <italic>Lactobacillus</italic> spp. and alpha diversity, and anaerobes like <italic>Atopobium</italic>, <italic>Clostridium</italic>, <italic>Megasphaera</italic>, <italic>Mycoplasma</italic>, <italic>Prevotella</italic>, and <italic>Sneathia</italic>. In contrast, no correlations were observed between richness, diversity and inflammation-associated cytokines and chemokines in HPV-neg women (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>). At the same time, positive correlations (rho&gt;0.40, <italic>p</italic>&lt;0.001) were found between alpha diversity indices and most anaerobes in both HPV-pos and HPV-neg women, however fewer correlations were found compared to HPV-pos women.</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Correlations between the vaginal microbiota and genital inflammation in women with symptomatic cervical ectopy with <bold>(A)</bold> or without <bold>(B)</bold> HPV infection. Correlations were computed using Spearman test, and visualized using corrplot. Only <italic>p</italic> values&lt; 0.01 are shown. Positive rho values are shown in blue, negative rho values shown in red, the intensity of the color is proportional to the strength of the rho value.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcimb-12-884272-g005.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>Discussion</title>
<p>Here, we present data on the prevalence of HPV infection and distribution of HR-HPV types, the presence and magnitude of genital inflammation in cervicovaginal mucus, the diversity and taxonomic composition of the vaginal microbiota, and associations in women with symptomatic cervical ectopy, a common gynaecological condition, associated with hormonal changes, that either resolves by itself or requires treatment when symptoms (bleeding and recurrent vaginal infections) persist.</p>
<p>The presence of cervical ectopy has been associated with an increased risk of sexually transmitted infections (<xref ref-type="bibr" rid="B35">Lee et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B41">Monroy et&#xa0;al., 2010</xref>). Here, we detected the presence of HPV in more than half of women with symptomatic cervical ectopy, in particular in younger women (20-30 years of age). Our results are comparable to other studies (<xref ref-type="bibr" rid="B44">Nielsen et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B10">Bruni et&#xa0;al., 2010</xref>). We found high prevalence of single and multiple HR-HPV and pHR-HPV types, including HPV16, HPV31 and HPV18. We and others have previously reported the presence of oncogenic HPV types in women with cervical ectopy (<xref ref-type="bibr" rid="B41">Monroy et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B24">Hwang et&#xa0;al., 2011</xref>), however, since endocervical cells are unable to sustain productive HPV life cycle, the clinical relevance and implications of HR-HPV infection in everted endocervical epithelium remain controversial. Interestingly, a refined model of cervical carcinogenesis, that considers the participation of endocervical cells as potential precursors of cervical cancer has been proposed (<xref ref-type="bibr" rid="B20">Doorbar and Griffin, 2019</xref>). In line with this thinking, cervical ectopy represents a challenge, because in this condition endocervical cells are exposed to the complex vaginal microenvironment containing chemical, immune and biological elements that may alter the natural evolution of columnar epithelium to metaplastic squamous tissue.</p>
<p>In women with cervical ectopy, endocervical cells protruding onto the ectocervix create an additional mucosal barrier with a natural vulnerability to physical trauma that may lead to increased genital inflammation in sexually active women (<xref ref-type="bibr" rid="B31">Kyongo et&#xa0;al., 2012</xref>). Accordingly, we found the co-expression of five important mediators of inflammation, IL-1&#x3b2;, IL-22, IL-8, CCL17, and CCL3 in the vast majority of women included in this study. Our results are consistent with previous reports showing increased levels of cervicovaginal pro-inflammatory cytokines in healthy women with cervical ectopy (<xref ref-type="bibr" rid="B24">Hwang et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B31">Kyongo et&#xa0;al., 2012</xref>). However, the impact of HPV infection on the inflammatory responses in women with symptomatic cervical ectopy has not been previously evaluated. We could not detect differences in pro-inflammatory cytokines between HPV-pos and HPV-neg women, suggesting that HPV infection does not increase inflammation further. Our results agree with Kriek and colleagues (<xref ref-type="bibr" rid="B30">Kriek et&#xa0;al., 2016</xref>), who did not find an association between pro-inflammatory cytokines and HPV in endocervical samples. In contrast, a pro-inflammatory microenvironment has been identified in squamous epithelium persistently infected by HPV, where E5, E6 and E7 oncoproteins were found to be the main activator of inflammation-associated pathways (<xref ref-type="bibr" rid="B22">Hemmat and Bannazadeh, 2019</xref>). Glandular cells are unable to sustain HPV replication cycle, and expression of oncoproteins. This might explain why the presence of HPV is not associated with increased genital inflammation in our cohort.</p>
<p>HPV is one of the most common pathogens affecting the cervical canal. Cellular immune responses are decisive for the resolution of HPV infections (<xref ref-type="bibr" rid="B56">Stanley, 2006</xref>). Here, we detected increased levels of IL-21 and CXCL9 in HPV-pos women. IL-21 is a crucial factor for the immunological control of viral infections. IL-21 is produced by CD4+ T lymphocytes and natural killer T (NKT) cells, and it modulates the anti-viral activities of cytotoxic T lymphocytes (CTL) and natural killer (NK) cells (<xref ref-type="bibr" rid="B55">Spolski and Leonard, 2008</xref>). IL-21 may contribute to both initial immune responses and control of HPV-associated genital warts (<xref ref-type="bibr" rid="B1">Abu El-Hamd et&#xa0;al., 2019</xref>). Thus, it is reasonable to assume that elevated concentrations of IL-21 detected in our cohort, might reflect a cellular immune response against HPV. As already stated, IL-21 is produced by NKT cells. Migration of NKT cells to sites of infection and inflammation is mediated by chemokines such as CXCL9 (<xref ref-type="bibr" rid="B17">Ding et&#xa0;al., 2016</xref>). We detected increased levels of CXCL9 in HPV-pos women, which may be stimulating the recruitment of IL-21-producing NKT cells. An association of HPV infection with elevated levels of cervicovaginal CXCL9 was also reported by Shannon and colleagues (<xref ref-type="bibr" rid="B54">Shannon et&#xa0;al., 2017</xref>), however, gynecological conditions of women included in that study were not reported. Interestingly, the expression of CXCL9 can also be induced by bacterial lipopolysaccharide (LPS) (<xref ref-type="bibr" rid="B8">Bandow et&#xa0;al., 2012</xref>), suggesting that some bacteria within the vaginal microbiota may also influence the expression of chemokines and as a consequence the recruitment of cells to this site.</p>
<p>As shown by others, bacterial communities in the vagina are able to influence host mucosal immune responses (<xref ref-type="bibr" rid="B3">Agostinis et&#xa0;al., 2019</xref>). A vaginal microbiota is considered optimal if dominated by <italic>Lactobacillus</italic> species and a low microbial diversity. We found that the vaginal microbiota of HPV-pos women was highly diverse and differed in community composition (beta diversity) with a decrease in <italic>Lactobacillus</italic> dominance compared to HPV-neg women. We also found a strong correlation between alpha diversity metrics and IL-1&#x3b2; and IL-22, both of which are produced at sites of inflammation. These results are in line with a previous study showing correlations between highly diverse bacterial communities and genital pro-inflammatory cytokines (<xref ref-type="bibr" rid="B5">Anahtar et&#xa0;al., 2015</xref>). Also, changes in vaginal microbial composition have been shown to associate with female genital tract conditions including pelvic inflammatory disease (<xref ref-type="bibr" rid="B43">Ness et&#xa0;al., 2005</xref>). To our knowledge, this is the first report showing associations between specific pro-inflammatory cytokines, and increased bacterial diversity in women with symptomatic cervical ectopy in the context of HPV infection. Our results show that high bacterial diversity correlates with both IL-1&#x3b2; and IL-22. Also, NKT-recruiting factor CXCL9 was increased in the cervical mucus of HPV-pos women. Both results are of great interest to us. IL-1&#x3b2; induces the production of IL-22 by NKT cells as a response to bacterial invasion (<xref ref-type="bibr" rid="B18">Doisne et&#xa0;al., 2011</xref>). Likewise, levels of IL-22 are increased in women with HPV-associated genital warts (<xref ref-type="bibr" rid="B50">Radwa et&#xa0;al., 2021</xref>), thus it is likely that both, the altered bacterial communities and the presence of HPV may contribute to the elevated levels of IL-22 detected in our cohort.</p>
<p>The link between vaginal dysbiosis and HPV infection, persistence, premalignant lesions (<xref ref-type="bibr" rid="B33">Laniewski et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B57">Usyk et&#xa0;al., 2020</xref>), and ultimately cervical cancer has been clearly established (<xref ref-type="bibr" rid="B11">Brusselaers et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B61">Xie et&#xa0;al., 2020</xref>). Given that cervical ectopy is not considered a malignant condition, the potential association between HPV and vaginal dysbiosis in women with this condition remains unclear. As stated previously, a <italic>Lactobacillus</italic>-dominant is considered an optimal vaginal microbiota, whereas the presence of strict anaerobes is considered a state of dysbiosis (<xref ref-type="bibr" rid="B48">Pramanick et&#xa0;al., 2022</xref>). In our cohort, HPV-neg women were mostly <italic>Lactobacillus</italic> dominant, whereas HPV-pos women showed a significant decreased in Lactobacilli and increase in anaerobes such as <italic>Sneathia, Pseudomonas, Megasphaera, Atopobium</italic>, <italic>Shuttleworthia</italic>, <italic>Prevotella</italic>, and <italic>Clostridium</italic>. The presence of anaerobes is associated with high rates of bacterial vaginosis recurrence (<xref ref-type="bibr" rid="B46">Onderdonk et&#xa0;al., 2016</xref>) and cervicovaginal inflammation (<xref ref-type="bibr" rid="B40">Mitchell and Marrazzo, 2014</xref>). Furthermore, byproducts of dysbiotic bacteria metabolism may cause degradation of the mucosal surfaces of the glandular epithelium (<xref ref-type="bibr" rid="B2">Africa et&#xa0;al., 2014</xref>), and increase the vaginal pH in a number of women (<xref ref-type="bibr" rid="B19">Donders et&#xa0;al., 2016</xref>). In contrast, most <italic>Lactobacillus</italic> spp. are known to produce lactic acid and favor a low pH in the vagina. Importantly, low pH (pH 4.5) is a critical factor for the natural replacement of ectopic tissue with metaplastic squamous epithelium. Thus, the loss or decrease in the relative abundance of Lactobacilli may delay the recovery of the external stratified epithelium in women with cervical ectopy, possibly increasing their risk of acquiring STDs, including HPV and HIV.</p>
<p>Our results should be interpreted with caution due to the limitations in statistical power inherent to small cohort sizes like ours. We also acknowledge that we did not adjust for confounding variables; when comparing clinical and demographic variables between these two groups, we found that the number of sexual partners was significantly higher in HPV-positive women as previously reported. All other variables were not significant, with the exception of age. There was a tendency for HPV-negative women to be older than HPV-positive women.</p>
<p>We are aware that our study has important limitations; first, the inclusion of a control group of women without cervical ectopy was not possible, so we cannot compare our results with those of age-matched women without cervical ectopy. However, previous reports show that cytokines, chemokines and vaginal microbial community profiles are different in women with cervical ectopy compared with women with a normal ectocervix (<xref ref-type="bibr" rid="B24">Hwang et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B31">Kyongo et&#xa0;al., 2012</xref>). Taking this latter in consideration, we hypothesize that our results might have been in line with those previously published. Another limitation is the conceivable presence of HPV-infected squamous cells, in the biopsy used for HPV DNA determination. Although we performed an inspection of the cervix after acetic acid application, and we separated women showing an acetowhite lesion, we cannot exclude the possibility of an ongoing exocervical HPV infection with no clinical manifestations in some women. We are aware that recovering pure glandular tissue from biopsies using microdissection would be the ideal methodology. The cross-sectional nature of our study prevents us from proposing any associated mechanisms or causative relationships. To this end, longitudinal studies are in process. Finally, we are aware of the difficulty of classifying to species level with this approach (16S), so we only reported taxonomic profiles at the genus level.</p>
<p>In conclusion, we showed the presence of oncogenic HPV types in a high proportion of women with symptomatic cervical ectopy. Notably, HPV was associated with increased levels of IL-21 and CXCL9 suggesting an active cellular immune response. We also found the co-expression of inflammation markers in these women irrespective of HPV status. We identified that changes in vaginal microbiota composition were associated with the presence of HPV. <italic>Lactobacillus</italic> dominance decreased in HPV-pos women, and bacterial vaginosis-associated anaerobes were increased. Correlations between inflammation-associated cytokines IL-1&#x3b2; and IL-22 and the vaginal microbiota were exclusively found in HPV-pos women. Taken together, our observations suggest that HR-HPV infection and vaginal dysbiosis coexist in women with symptomatic cervical ectopy. These women might be more susceptible to chronic inflammation and epithelial damage, and could benefit from a closer follow-up to improve their sexual and reproductive health. Our study provides a deeper understanding of the cervicovaginal microenvironment of women with symptomatic cervical ectopy, a benign gynecological condition.</p>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data Availability Statement</title>
<p>The raw 16S sequences have been deposited at the National Center for Biotechnology Information (NCBI)-Sequence Read Archive (SRA) under project PRJNA763215.</p>
</sec>
<sec id="s6" sec-type="ethics-statement">
<title>Ethics Statement</title>
<p>The studies involving human participants were reviewed and approved by Local Ethics Committee of the Colposcopy Clinic of the &#x201c;Cruz-Talonia&#x201d; Foundation (Reference Number: 17/02/2017/LRZ/CE05). All women gave written informed consent to participate in this study.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author Contributions</title>
<p>LR-Z conceived the project, designed experiments, coordinated study operation, supervised laboratory staff, analyzed and discussed data, and wrote the manuscript. SP-C conducted 16S rRNA gene sequencing analysis, designed figures, analyzed and discussed data, and wrote the manuscript. ML-F performed flow cytometry analysis, sequencing reactions, analyzed data, designed figures and wrote the manuscript. FC generated clinical data base, processed tissue samples, extracted and analyzed DNA, and designed tables. TG conducted the E7-MPG assay to determine HPV genotypes. OC and FC-T coordinated study operation ant the Colposcopy Clinic of the &#x201c;Cruz-Talonia&#x201d; Foundation, and obtained tissue samples. MC-T extracted DNA from cervicovaginal mucus, constructed bacterial libraries, and performed 16S rRNA sequencing. CA-G analyzed cases and contributed to diagnosis confirmation. IM-H performed flow cytometry for chemokine determination. JM and VC-M performed flow cytometry analysis of cytokines. CA participated in the selection of patients, collection of patient data and processing of patient samples. MT coordinated study operation at the International Agency for Research in Cancer, analyzed and discussed data. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>This work was supported by a grant from PAPIIT, UNAM (grant number IN204419). LR-Z received a fellowship from Programa de Apoyos para la Superaci&#xf3;n del Personal Acad&#xe9;mico de la UNAM (PASPA. Number 523.01/2299DFA/2017). Grants from CONACyT (grant numbers 620702, 464745) were awarded to FC. SP-C was supported by funds from the Mexican Government (Programa Presupuestal P016, Anexo 13 del Decreto del Presupuesto de Egresos de la Federaci&#xf3;n).</p>
</sec>
<sec id="s9">
<title>Author Disclaimer</title>
<p>Where authors are identified as personnel of the International Agency for Research on Cancer/World Health Organization, the authors alone are responsible for the views expressed in this article and they do not necessarily represent the decisions, policy or views of the International Agency for Research on Cancer/World Health Organization.</p>
</sec>
<sec id="s10" sec-type="COI-statement">
<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 id="s11" sec-type="disclaimer">
<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>
</body>
<back>
<ack>
<title>Acknowledgments</title>
<p>First, the authors would like to thank all women taking part in this study. The authors wish to thank Mrs. Sandrine McKay-Chopin for excellent technical assistance and Dr. Santiago &#xc1;vila-Rios, Centro de Investigaci&#xf3;n en Enfermedades Infecciosas, Instituto Nacional de Enfermedades Respiratorias, for scientific support.</p>
</ack>
<sec id="s12" sec-type="supplementary-material">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fcimb.2022.884272/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fcimb.2022.884272/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Image_1.pdf" id="SM1" mimetype="application/pdf"/>
<supplementary-material xlink:href="Image_2.pdf" id="SM2" mimetype="application/pdf"/>
<supplementary-material xlink:href="Image_3.pdf" id="SM3" mimetype="application/pdf"/>
<supplementary-material xlink:href="Table_1.docx" id="ST1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Table_2.docx" id="ST2" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Table_3.docx" id="ST3" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Table_4.docx" id="ST4" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abu El-Hamd</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Assaf</surname> <given-names>H. A. R.</given-names>
</name>
<name>
<surname>Sedky</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Mohammad</surname> <given-names>S. H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Possible Role of Interleukin 21 and Interleukin 33 in Patients With Genital Warts</article-title>. <source>Dermatol. Ther.</source> <volume>32</volume>, <elocation-id>e13063</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/dth.13063</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Africa</surname> <given-names>C. W. J.</given-names>
</name>
<name>
<surname>Nel</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Stemmet</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Anaerobes and Bacterial Vaginosis in Pregnancy: Virulence Factors Contributing to Vaginal Colonization</article-title>. <source>Int. J. Environ. Res. Public Health</source> <volume>11</volume>, <fpage>6979</fpage>&#x2013;<lpage>7000</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/ijerph110706979</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Agostinis</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Mangogna</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Bossi</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Ricci</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Kishore</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Bulla</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Uterine Immunity and Microbiota : A Shifting Paradigm</article-title>. <source>Front. Immunol.</source> <volume>10</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fimmu.2019.02387</pub-id>
</citation>
</ref>
<ref id="B4">
<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>D. O. C.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The Vaginal Microenvironment: The Physiologic Role of <italic>Lactobacilli</italic>
</article-title>. <source>Front. Med. (Lausanne)</source> <volume>5</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmed.2018.00181</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Anahtar</surname> <given-names>M. N.</given-names>
</name>
<name>
<surname>Byrne</surname> <given-names>E. H.</given-names>
</name>
<name>
<surname>Doherty</surname> <given-names>K. E.</given-names>
</name>
<name>
<surname>Bowman</surname> <given-names>B. A.</given-names>
</name>
<name>
<surname>Yamamoto</surname> <given-names>H. S.</given-names>
</name>
<name>
<surname>Soumillon</surname> <given-names>M.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Cervicovaginal Bacteria are Major Modulator of Host Inflammatory Responses in the Female Genital Tract</article-title>. <source>Immunity</source> <volume>42</volume>, <fpage>965</fpage>&#x2013;<lpage>976</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.immuni.2015.04.019</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Anahtar</surname> <given-names>M. N.</given-names>
</name>
<name>
<surname>Gootenberg</surname> <given-names>D. B.</given-names>
</name>
<name>
<surname>Mitchell</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>Kwon</surname> <given-names>D. S.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Cervicovaginal Microbiota and Reproductive Health: The Virtue of Simplicity</article-title>. <source>Cell Host Microbe</source> <volume>23</volume>, <fpage>159</fpage>&#x2013;<lpage>168</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.chom.2018.01.013</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Andrade Pessoa Morales</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Marconi</surname> <given-names>C.</given-names>
</name>
<name>
<surname>El-Zein</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ravel</surname> <given-names>J.</given-names>
</name>
<name>
<surname>da Silva Pinto</surname> <given-names>G. V.</given-names>
</name>
<name>
<surname>Silveira</surname> <given-names>R.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Vaginal Microbiome Components as Correlates of Cervical Human Papillomavirus Infection</article-title>. <source>J. Infect. Dis.</source>, <fpage>jiab547</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/infdis/jiab547</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bandow</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kusuyama</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Shamoto</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Kakimoto</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Ohnishi</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Matsuguchi</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>LPS-Induced Chemokine Expression in Both MyD88-Dependent and &#x2013;Independent Manners is Regulated by Cot/Tpl2-ERK Axis in Macrophages</article-title>. <source>FEBS Lett.</source> <volume>586</volume>, <fpage>1540</fpage>&#x2013;<lpage>1546</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.febslet.2012.04.018</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bolyen</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Rideout</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Dillon</surname> <given-names>M. R.</given-names>
</name>
<name>
<surname>Bokulich</surname> <given-names>N. A.</given-names>
</name>
<name>
<surname>Abnet</surname> <given-names>C. C.</given-names>
</name>
<name>
<surname>Al-Ghalith</surname> <given-names>G. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Reproducible, Interactive, Scalable and Extensible Microbiome Data Science Using QIIME 2</article-title>. <source>Nat. Biotechnol.</source> <volume>37</volume>, <fpage>852</fpage>&#x2013;<lpage>857</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41587-019-0209-9</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bruni</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Diaz</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Castellsagu&#xe9;</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Ferrer</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Bosch</surname> <given-names>F. X.</given-names>
</name>
<name>
<surname>de Sanjos&#xe9;</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Cervical Human Papillomavirus Prevalence in 5 Continents: Meta-Analysis of 1 Million Women With Normal Cytological Findings</article-title>. <source>J. Infect. Dis.</source> <volume>202</volume>, <fpage>1789</fpage>&#x2013;<lpage>1799</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1086/657321</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brusselaers</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Shrestha</surname> <given-names>S.</given-names>
</name>
<name>
<surname>van de Wijgert</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Verstraelen</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Vaginal Dysbiosis and the Risk of Human Papillomavirus and Cervical Cancer: Systematic Review and Meta-Analysis</article-title>. <source>Am. J. Obstet. Gynecol</source>. <volume>221</volume>, <fpage>9</fpage>&#x2013;<lpage>18</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ajog.2018.12.011</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Callahan</surname> <given-names>B. J.</given-names>
</name>
<name>
<surname>McMurdie</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Rosen</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>A. W.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>A. J.</given-names>
</name>
<name>
<surname>Holmes</surname> <given-names>S. P.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>DADA2: High-Resolution Sample Inference From Illumina Amplicon Data</article-title>. <source>Nat. Methods</source> <volume>13</volume>, <fpage>581</fpage>&#x2013;<lpage>583</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nmeth.3869</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Caporaso</surname> <given-names>J. G.</given-names>
</name>
<name>
<surname>Lauber</surname> <given-names>C. L.</given-names>
</name>
<name>
<surname>Walters</surname> <given-names>W. A.</given-names>
</name>
<name>
<surname>Berg-Lyons</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Lozupone</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Turnbaugh</surname> <given-names>P. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Global Patterns of 16S rRNA Diversity at a Depth of Millions of Sequences Per Sample</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>108 Suppl 1</volume>, <fpage>4516</fpage>&#x2013;<lpage>4522</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.1000080107</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cekmez</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Sanlikan</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Gocmen</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Vural</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Turkmen</surname> <given-names>S. B.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Is Cryotherapy Friend or Foe for Symptomatic Cervical Ectopy</article-title>? <source>Med. Princ. Pract.</source> <volume>25</volume>, <fpage>8</fpage>&#x2013;<lpage>11</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1159/000441433</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Clifford</surname> <given-names>G. M.</given-names>
</name>
<name>
<surname>Vaccarella</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Franceschi</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Tenet</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Umulisa</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Tshomo</surname> <given-names>U.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>Comparison of Two Widely Used Human Papillomavirus Detection and Genotyping Methods, GP5+/6+ -Based PCR Followed by Reverse Line Blot Hybridization and Multiplex Type-Specific E7-Based PCR</article-title>. <source>J. Clin. Microbiol.</source> <volume>54</volume>, <fpage>2031</fpage>&#x2013;<lpage>2038</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/JCM.00618-16</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Daniilidis</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Koutsos</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Oikonomou</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Nasioutziki</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Hatziparadisi</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Tantanasis</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Cytokines of Cervical Mucosa and Human Papillomavirus Infection of the Cervix: A Descriptive Study</article-title>. <source>Acta Cytol.</source> <volume>60</volume>, <fpage>58</fpage>&#x2013;<lpage>64</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1159/000445161</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ding</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Xia</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Ding</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Fan</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>CXCL9: Evidence and Contradictions for its Role in Tumor Progression</article-title>. <source>Cancer Med.</source> <volume>5</volume>, <fpage>3246</fpage>&#x2013;<lpage>3259</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/cam4.934</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Doisne</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Soulard</surname> <given-names>V.</given-names>
</name>
<name>
<surname>B&#xe9;court</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Amniai</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Henrot</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Havenar-Daughton</surname> <given-names>C.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Cutting Edge: Crucial Role of IL-1 and IL-23 in the Innate IL-17 Response of Peripheral Lymph Node NK1-1-Invariant NKT Cells to Bacteria</article-title>. <source>J. Immunol.</source> <volume>186</volume>, <fpage>662</fpage>&#x2013;<lpage>666</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4049/jimmunol.1002725</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Donders</surname> <given-names>G. G. G.</given-names>
</name>
<name>
<surname>Gonzaga</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Marconi</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Donders</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Michiels</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Eggermont</surname> <given-names>N.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>Increased Vaginal pH in Ugandan Women: What Does it Indicate</article-title>? <source>Eur. J. Clin. Microbiol. Infect. Dis.</source> <volume>35</volume>, <fpage>1297</fpage>&#x2013;<lpage>1303</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10096-016-2664-2</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Doorbar</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Griffin</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Refining Our Understanding of Cervical Neoplasia and its Cellular Origins</article-title>. <source>Papillomavirus Res.</source> <volume>7</volume>, <fpage>176</fpage>&#x2013;<lpage>179</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.pvr.2019.04.005</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Harry</surname> <given-names>V. N.</given-names>
</name>
<name>
<surname>Cruickshank</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Parkin</surname> <given-names>D. E.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Auditing the Use of Colposcopy Versus General Gynecology Clinics to Investigate Women With Postcoital or Intermenstrual Bleeding: A Case for a New Outpatient Service</article-title>. <source>J. Low. Genit. Tract Dis.</source> <volume>11</volume>, <fpage>108</fpage>&#x2013;<lpage>111</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1097/01.lgt.0000245049.67965.39</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hemmat</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Bannazadeh</surname> <given-names>B. H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Association of Human Papillomavirus Infection and Inflammation in Cervical Cancer</article-title>. <source>Pathog. Dis.</source> <volume>77</volume>, <elocation-id>ftz048</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/femspd/ftz048</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hua</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Zeng</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Diao</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Using Platelet-Rich Plasma for the Treatment of Symptomatic Cervical Ectopy</article-title>. <source>Int. J. Gynecol. Obstet.</source> <volume>119</volume>, <fpage>26</fpage>&#x2013;<lpage>29</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijgo.2012-05.029</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hwang</surname> <given-names>L. Y.</given-names>
</name>
<name>
<surname>Scott</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Moscicki</surname> <given-names>A. B.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Higher Levels of Cervicovaginal Inflammatory and Regulatory Cytokines and Chemokines in Healthy Young Women With Immature Cervical Epithelium</article-title>. <source>J. Reprod. Immunol.</source> <volume>88</volume>, <fpage>66</fpage>&#x2013;<lpage>71</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jri.2010.07.008</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>IARC</collab>
</person-group> (<year>2007</year>). &#x201c;<article-title>Human Papillomaviruses</article-title>,&#x201d; in <source>IARC Monogr. Eval. Carcinog. Risks Hum</source>. (<publisher-loc>Geneva, Switzerland</publisher-loc>: <publisher-name>WHO Press</publisher-name>), vol. <volume>90</volume>, <fpage>1</fpage>&#x2013;<lpage>636</lpage>.</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Iwata</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Fujii</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Morii</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Saito</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Sugiyama</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Nishio</surname> <given-names>H.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Cytokine Profile in Cervical Mucosa of Japanese Patients With Cervical Intraepithelial Neoplasia</article-title>. <source>Int. J. Clin. Oncol.</source> <volume>20</volume>, <fpage>126</fpage>&#x2013;<lpage>133</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10147-014-0680-8</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jacobson</surname> <given-names>D. L.</given-names>
</name>
<name>
<surname>Peralta</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Graham</surname> <given-names>N. M. H.</given-names>
</name>
<name>
<surname>Zenilman</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Histologic Development of Cervical Ectopy. Relationship to Reproductive Hormones</article-title>. <source>Sex Transm. Dis.</source> <volume>27</volume>, <fpage>252</fpage>&#x2013;<lpage>258</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1097/00007435-200005000-00003</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Katoh</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Misawa</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kuma</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Miyata</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>MAFFT: A Novel Method for Rapid Multiple Sequence Alignment Based on Fast Fourier Transform</article-title>. <source>Nucleic Acids Res.</source> <volume>30</volume>, <fpage>3059</fpage>&#x2013;<lpage>3066</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkf436</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kemp</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Hildesheim</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Garc&#xed;a-Pi&#xf1;eres</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Williams</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Shearer</surname> <given-names>G. M.</given-names>
</name>
<name>
<surname>Rodriguez</surname> <given-names>A. C.</given-names>
</name>
<etal/>
</person-group>. (<year>2010</year>). <article-title>Elevated Systemic Levels of Inflammatory Cytokines in Older Women With Persistent Cervical Human Papillomavirus Infection</article-title>. <source>Cancer Epidemiol. Biomarkers Prev.</source> <volume>19</volume>, <fpage>1954</fpage>&#x2013;<lpage>1959</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1158/1055-9965.EPI-10-0184</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kriek</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Jaumdally</surname> <given-names>S. Z.</given-names>
</name>
<name>
<surname>Masson</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Little</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Mbulawa</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Gumbi</surname> <given-names>P. P.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>Female Genital Tract Inflammation, HIV Co-Infection and Persistent Mucosal Human Papillomavirus (HPV) Infections</article-title>. <source>Virology</source> <volume>493</volume>, <fpage>247</fpage>&#x2013;<lpage>254</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.virol.2016.03.02</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kyongo</surname> <given-names>J. K.</given-names>
</name>
<name>
<surname>Jespers</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Goovaerts</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Michiels</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Menten</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Fichorova</surname> <given-names>R. N.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>Searching for Lower Female Genital Tract Soluble and Cellular Biomarkers: Defining Levels and Predictors in a Cohort of Healthy Caucasian Women</article-title>. <source>PloS One</source> <volume>7</volume>, <elocation-id>e43951</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0043951</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lamont</surname> <given-names>R. F.</given-names>
</name>
<name>
<surname>Sobel</surname> <given-names>J. D.</given-names>
</name>
<name>
<surname>Akins</surname> <given-names>R. A.</given-names>
</name>
<name>
<surname>Hassan</surname> <given-names>S. S.</given-names>
</name>
<name>
<surname>Chaiworapongsa</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Kusanovic</surname> <given-names>J. P.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>The Vaginal Microbiome: New Information About Genital Tract Flora Using Molecular Based Techniques</article-title>. <source>BJOG</source> <volume>118</volume>, <fpage>533</fpage>&#x2013;<lpage>549</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1471-0528.2010.02840.x</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Laniewski</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Barnes</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Goulder</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Cui</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Roe</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Chase</surname> <given-names>D. M.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Linking Cervicovaginal Immune Signatures, HPV and Microbiota Composition in Non-Hispanic and Hispanic Women</article-title>. <source>Sci. Rep.</source> <volume>8</volume>, <fpage>7593</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-018-25879-7</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Laniewski</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Ilhan</surname> <given-names>Z. E.</given-names>
</name>
<name>
<surname>Herbst-Kralovetz</surname> <given-names>M. M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>The Microbiome and Gynaecological Cancer Development, Prevention and Therapy</article-title>. <source>Nat. Rev. Urol.</source> <volume>17</volume>, <fpage>232</fpage>&#x2013;<lpage>250</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41585-020-0286-z</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Tobin</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Foley</surname> <given-names>E.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Relationship of Cervical Ectopy to Chlamydia Infection in Young Women</article-title>. <source>J. Fam. Plann. Reprod. Health Care</source> <volume>32</volume>, <fpage>104</fpage>&#x2013;<lpage>106</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1783/147118906776276440</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lozupone</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Lladser</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Knights</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Stombaugh</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Knight</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>UniFrac: An Effective Distance Metric for Microbial Community Comparison</article-title>. <source>ISME J.</source> <volume>5</volume>, <fpage>169</fpage>&#x2013;<lpage>172</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/ismej.2010.133</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Machado</surname> <given-names>J. L. C.</given-names>
</name>
<name>
<surname>Whitaker</surname> <given-names>D. A. S.</given-names>
</name>
<name>
<surname>Barbosa</surname> <given-names>C. H.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Evidence for Benefits From Treating Cervical Ectopy: Literature Review</article-title>. <source>Sao Paulo Med. J.</source> <volume>126</volume>, <fpage>132</fpage>&#x2013;<lpage>139</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1590/s1516-31802008000200014</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McDonald</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Price</surname> <given-names>M. N.</given-names>
</name>
<name>
<surname>Goodrich</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Nawrocki</surname> <given-names>E. P.</given-names>
</name>
<name>
<surname>DeSantis</surname> <given-names>T. Z.</given-names>
</name>
<name>
<surname>Probst</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>An Improved Greengenes Taxonomy With Explicit Ranks for Ecological and Evolutionary Analyses of Bacteria and Archaea</article-title>. <source>ISME J.</source> <volume>6</volume>, <fpage>610</fpage>&#x2013;<lpage>618</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/ismej.2011.139</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McMurdie</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Holmes</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data</article-title>. <source>PloS One</source> <volume>8</volume>, <elocation-id>e61217</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0061217</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mitchell</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Marrazzo</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Bacterial Vaginosis and the Cervicovaginal Immune Response</article-title>. <source>Am. J. Reprod. Immunol.</source> <volume>71</volume>, <fpage>555</fpage>&#x2013;<lpage>563</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/aji.12264</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Monroy</surname> <given-names>O. L.</given-names>
</name>
<name>
<surname>Aguilar</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Lizano</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Cruz-Talonia</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Cruz</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Rocha-Zavaleta</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Prevalence of Human Papillomavirus Genotypes, and Mucosal IgA Anti-Viral Responses in Women With Cervical Ectopy</article-title>. <source>J. Clin. Virol.</source> <volume>47</volume>, <fpage>43</fpage>&#x2013;<lpage>48</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jcv.2009.10.008</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Myer</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wright</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Denny</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Kuhn</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Nested Case-Control Study of Cervical Mucosal Lesions, Ectopy, and Incident HIV Infection Among Women in Cape Town, South Africa</article-title>. <source>Sex Transm. Dis.</source> <volume>33</volume>, <fpage>683</fpage>&#x2013;<lpage>687</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1097/01.olq.0000216026.67352.f9</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ness</surname> <given-names>R. B.</given-names>
</name>
<name>
<surname>Kip</surname> <given-names>K. E.</given-names>
</name>
<name>
<surname>Hillier</surname> <given-names>S. L.</given-names>
</name>
<name>
<surname>Soper</surname> <given-names>D. E.</given-names>
</name>
<name>
<surname>Stamm</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Sweet</surname> <given-names>R. L.</given-names>
</name>
<etal/>
</person-group>. (<year>2005</year>). <article-title>A Cluster Analysis of Bacterial Vaginosis-Associated Microflora and Pelvic Inflammatory Disease</article-title>. <source>J. Epidemiol</source>. <volume>162</volume>, <fpage>585</fpage>&#x2013;<lpage>590</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/aje/kwi243</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nielsen</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Kjaer</surname> <given-names>S. K.</given-names>
</name>
<name>
<surname>Munk</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Iftner</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Type-Specific HPV Infection and Multiple HPV Types: Prevalence and Risk Factors Profile in Nearly 12,000 Younger and Older Danish Women</article-title>. <source>Sex Transm. Dis.</source> <volume>35</volume>, <fpage>276</fpage>&#x2013;<lpage>282</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1097/OLQ.0b013e31815ac5c7</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nieves-Ram&#xed;rez</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Partida-Rodr&#xed;guez</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Moran</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Serrano-V&#xe1;zquez</surname> <given-names>A.</given-names>
</name>
<name>
<surname>P&#xe9;rez-Ju&#xe1;rez</surname> <given-names>H.</given-names>
</name>
<name>
<surname>P&#xe9;rez-Rodr&#xed;guez</surname> <given-names>M. E.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Cervical Squamous Intraepithelial Lesions are Associated With Differences in the Vaginal Microbiota of Mexican Women</article-title>. <source>Microbiol. Spectr.</source> <volume>9</volume>, <elocation-id>e0014321</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/Spectrum.00143-21</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Onderdonk</surname> <given-names>A. B.</given-names>
</name>
<name>
<surname>Delaney</surname> <given-names>M. L.</given-names>
</name>
<name>
<surname>Fichorova</surname> <given-names>R. N.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The Human Microbiome During Bacterial Vaginosis</article-title>. <source>Clin. Microbiol. Rev.</source> <volume>29</volume>, <fpage>223</fpage>&#x2013;<lpage>238</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/CMR.00075-15</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Otani</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Fujii</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Kukimoto</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Yamamoto</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Tsukamoto</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Ichikawa</surname> <given-names>R.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Cytokine Expression Profiles in Cervical Mucus From Patients With Cervical Cancer and its Precursor Lesions</article-title>. <source>Cytokine</source> <volume>120</volume>, <fpage>210</fpage>&#x2013;<lpage>219</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.cyto.2019.05.011</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pramanick</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Nathani</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Warke</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Mayadeo</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Aranha</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Vaginal Dysbiotic Microbiome in Women With No Symptoms of Genital Infections</article-title>. <source>Front. Cell Infect. Microbiol.</source> <volume>11</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fcimb.2021.760459</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Price</surname> <given-names>M. N.</given-names>
</name>
<name>
<surname>Dehal</surname> <given-names>P. S.</given-names>
</name>
<name>
<surname>Arkin</surname> <given-names>A. P.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>FastTree 2-Approximately Maximum-Likelihood Trees for Large Alignments</article-title>. <source>PloS One</source> <volume>5</volume>, <elocation-id>e9490</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal,pone.0009490</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Radwa</surname> <given-names>E. M. M.</given-names>
</name>
<name>
<surname>Aya</surname> <given-names>Q. E. M. A.</given-names>
</name>
<name>
<surname>Fadia</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Maha</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Amal</surname> <given-names>H. A. G.</given-names>
</name>
<name>
<surname>Lina</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Serum Level of Interleukin-22 in Patients With Cutaneous Warts: A Case Control Study</article-title>. <source>J. Cosmet. Dermatol.</source> <volume>20</volume>, <fpage>1782</fpage>&#x2013;<lpage>1787</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/jocd.13779</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Reich</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Fritsch</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>The Developmental Origin of Cervical and Vaginal Epithelium and Their Clinical Consequences: A Systemic Review</article-title>. <source>J. Low Genit. Tract. Dis.</source> <volume>18</volume>, <fpage>358</fpage>&#x2013;<lpage>360</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1097/LGT.0000000000000023</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Scott</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Shvetsov</surname> <given-names>Y. B.</given-names>
</name>
<name>
<surname>Thompson</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Hernandez</surname> <given-names>B. Y.</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Wilkens</surname> <given-names>L. R.</given-names>
</name>
<etal/>
</person-group>. (<year>2013</year>). <article-title>Cervical Cytokines and Clearance of Incident Human Papillomavirus Infection: Hawaii HPV Cohort Study</article-title>. <source>Int. J. Cancer</source> <volume>133</volume>, <fpage>1187</fpage>&#x2013;<lpage>1196</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/ijc.28119</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Segata</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Izard</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Waldron</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Gevers</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Miropolsky</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Garrett</surname> <given-names>W. S.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Metagenomic Biomarker Discovery and Explanation</article-title>. <source>Genome Biol.</source> <volume>12</volume>, <fpage>R60</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/gb-2011-12-6-r60</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shannon</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Yi</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Perusini</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Gajer</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Humphrys</surname> <given-names>M. S.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Association of HPV Infection and Clearance With Cervicovaginal Immunology and the Vaginal Microbiota</article-title>. <source>Mucosal Immunol.</source> <volume>10</volume>, <fpage>1310</fpage>&#x2013;<lpage>1319</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/mi.2016.129</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Spolski</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Leonard</surname> <given-names>W. J.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Interleukin-21: Basic Biology and Implications for Cancer and Autoimmunity</article-title>. <source>Annu. Rev. Immunol.</source> <volume>26</volume>, <fpage>57</fpage>&#x2013;<lpage>77</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev.immunol.26.021607.090316</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stanley</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Immune Responses to Human Papillomavirus</article-title>. <source>Vaccine</source> <volume>24 Suppl 1</volume>, <fpage>S16</fpage>&#x2013;<lpage>S22</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.vaccine.2005.09.002</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Usyk</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Zolnik</surname> <given-names>C. P.</given-names>
</name>
<name>
<surname>Castle</surname> <given-names>P. E.</given-names>
</name>
<name>
<surname>Porras</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Herrero</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Grandissimo</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Cervicovaginal Microbiome and Natural History of HPV in a Longitudinal Study</article-title>. <source>PloS Pathog.</source> <volume>16</volume>, <elocation-id>e1008376</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.ppat.1008376</pub-id>
</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Walboomers</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Jacobs</surname> <given-names>M. V.</given-names>
</name>
<name>
<surname>Manos</surname> <given-names>M. M.</given-names>
</name>
<name>
<surname>Bosch</surname> <given-names>F. X.</given-names>
</name>
<name>
<surname>Kummer</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Shah</surname> <given-names>K. V.</given-names>
</name>
<etal/>
</person-group>. (<year>1999</year>). <article-title>Human Papillomavirus is a Necessary Cause of Invasive Cervical Cancer Worldwide</article-title>. <source>J. Pathol.</source> <volume>189</volume>, <fpage>12</fpage>&#x2013;<lpage>19</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/(SICI)1096-9896(199909)189:1&lt;12::AID-PATH431&gt;3.0.CO;2-F</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname> <given-names>Z. T.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H. L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C. F.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>G. L.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>S. M.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>S. L.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Depiction of Vaginal Microbiota in Women With High-Risk Human Papillomavirus Infection</article-title>. <source>Front. Public Health</source> <volume>8</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpubh.2020.587298</pub-id>
</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Witkin</surname> <given-names>S. S.</given-names>
</name>
<name>
<surname>Alvi</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Bongiovanni</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Linhares</surname> <given-names>I. M.</given-names>
</name>
<name>
<surname>Ledger</surname> <given-names>W. J.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Lactic Acid Stimulates Interleukin-23 Production by Peripheral Blood Mononuclear Cells Exposed to Bacterial Lipopolysaccharide</article-title>. <source>FEMS Immunol. Med. Microbiol.</source> <volume>61</volume>, <fpage>153</fpage>&#x2013;<lpage>158</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1574-695X.2010.00757.x</pub-id>
</citation>
</ref>
<ref id="B61">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Feng</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Zhan</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>H.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Revealing the Disturbed Vaginal Microbiota Caused by Cervical Cancer Using High-Throughput Sequencing Technology</article-title>. <source>Front. Cell. Infect. Microbiol.</source> <volume>10</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fcimb.2020538336</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>
