The interplay between HPV, other Sexually Transmissible Infections and genital microbiome on cervical microenvironment (MicroCervixHPV study)

Background The importance of Cervicovaginal Microbiota in protecting against infections (such as HPV) is already well established, namely through Lactobacillus spp., as well as the mechanism through which HPV leads to Cervical Neoplasia. However, it is not possible to classify HPV as a complete carcinogen. Thus, the importance of exploring Cervicovaginal dysbiosis with the intention of deciphering this interaction with HPV, takes on greater relevance. The main objectives of this study were: 1) Comparison of the MCV composition of women with or without HPV and women with ASCUS or LSIL; 2) Characterization of cytokines present in the vaginal microenvironment; 3) Evaluation of the blood count ratios as prognostic systemic inflammatory biomarkers; 4) Correlation between MCV, HPV serotypes and cytokines. Methods This was a retrospective, observational, multicenter, cross-sectional study. CVM analysis was performed by isolation RNA and sequencing on a NGS platform. Cytokine concentrations of CVM were obtained through Multiplex platform. Statistical analysis was performed in SPSS v 26.0. An α of 0.05 was considered statistically significant. Results Highlighting the core of the study, CVM types of CST I and CST IV were found to influence the emergence of cervical lesions. Neutrophil-to-Lymphocyte ratio was found to impact the prognosis of ASCUS. Within CVM, Lactobacillus prevent the growth of other CST IV species, while the latter express symbiotic relationships with each other and show affinity for specific HPV serotypes. At last, RANTES chemokine is significantly elevated in cervicovaginal infections. Conclusion The importance of using vaginal cytokine profiles and CVM is highlighted in the hypothesis of prevention of Cervical Neoplasia development, as well as in its use as a prognostic biomarker. Taken together, these insights are one step closer to personalized medicine.


Background
Sexually Transmissible Infections (STI) are a heterogeneous group of diseases, inherent to human condition, associated with the core function of a species survivalreproduction.One of the most common STI is Human Papillomavirus (HPV) persistent infection.HPV is responsible, on one hand by an acute infection, that can be cleared out, while on the other hand, can lead to a chronic infection, with possible progression to cervical dysplasia and subsequent cervical neoplasia (Alizon et al., 2017).
Cervical neoplasia (CN) remains the fourth more common cancer among women, right behind breast cancer.It is accountable for around 570 000 cases and 311 000 deaths per year, according to Global Cancer Observatory 2018 data.This type of cancer affects mainly middle-aged women, particularly in countries with less resources, and it represents the major source of cancer-related deaths in women living in almost all the African continent.(Arbyn et al., 2020).The main risk factors linked to CN are smoking, oral contraceptive use and multiple gestations.
HPV per se can be separated in low-risk (lr-HPV) and high-risk HPV (hr-HPV), based on carcinogenic capability, with lr-HPV being correlated with anogenital warts, whereas hr-HPV is linked to cervical intraepithelial neoplasia (CIN) and CN, maNLRy due to HPV-16 and HPV-18, which represents nearly 70% of all CN.(Clifford et al., 2003;Clifford et al., 2005;Braaten and Laufer, 2008;Xiong et al., 2020) Since the onset of HPV infection, and if clearance of the virus is absent, a period of time of roughly 15 years, happens until the development of CN, walking through defined stages of CIN, corresponding to a continuous oncogenic process (Moscicki et al., 2012;Kyrgiou et al., 2017).
In the cervical microenvironment microbiome arises from vertical transmission, at the time of vaginal delivery and from birth until adulthood, microbiota is unstable and suffer from considerable modifications, being highly susceptible to external influences, such as the mode of delivery (vaginal vs. C-section), ethnicity, mode of lactation (breastfeeding vs. bottle-feeding), household exposures (e.g.: siblings, pets, among others), antibiotic usage, chronic diseases, inter alia.When the microbiome gets more stable, it that can only be marginally modified by diet, level of physical activity, habits (tobacco, alcohol consumption) and pathologies (Stewart et al., 2018;Senn et al., 2020;Alpuim Costa et al., 2021).
Focusing on cervicovaginal microbiota (CVM), after the advent of 16s rRNA High-throughput sequencing technologies (16s-HTS) Ravel et al. was able to identify and divide CVM into 5 groups, baptized as Community State Types (CST).(Ravel et al., 2011) This subdivision is attributable to the low diversity of CVM (in a homeostatic state) and considers the more predominant bacteria in each composition.Thus, the main bacteria in CST I is Lactobacillus crispatus, in CST II is L. gasseri, in CST III is L. iners and in CST V is L. jensenii.In CST IV, the microbiota is slightly different, with a resemblance to bacterial vaginosis, presenting high diversity, with most anaerobic species, such as Gardenerella, Megasphera, Sneathia, Atopobium and Prevotella, combined with lower levels of Lactobacillus (CST I, II, III and IV).This composition is consistent with a dysbiotic state (Olmsted et al., 2003;Kyrgiou et al., 2017;Ilhan et al., 2019;Curty et al., 2020;Norenhag et al., 2020).
The "healthy" CVM comprises low diversity with widely domination of Lactobacillus species, as stated above.The main reason for this scenario is supported by the ability of these bacteria to produce lactic acid and hydrogen peroxide, which exert a growth control on other bacterial species and viruses (Kyrgiou et al., 2017;Tachedjian et al., 2017;Curty et al., 2020;van de Wijgert and Verwijs, 2020) CVM is modulated by hormonal changeability, associated with the reproductive cycle, acquiring a more unstable composition, during menstruation (Ravel et al., 2011;Gajer et al., 2012) and contrarywise greater stability and less diversity during pregnancy (MacIntyre et al., 2015;Kyrgiou et al., 2017).Other factors (e.g.: sexual activity, lactation, use of oral contraceptives, vaginal douche, and stress) might influence the hormonal microenvironment of cervix and vagina.(Hickey et al., 2013;Amabebe and Anumba, 2018;Curty et al., 2020) Finally, the CST-IV bacteria can enhance the concentration of pro-inflammatory cytokines , increasing CCR5+ CD4+ T cells in vaginal mucosa, which are the targets of HIV (Gosmann et al., 2017;Torcia, 2019); as well as, increments SCFA production that, in turn, can regulate histone acetylation, which holds the possibility of reactivation of latent HIV-1 proviruses (Imai et al., 2012;Curty et al., 2020).
The primary outcome of our study was to understand the interplay between HPV infection, CVM composition and CIN stages, by comparing cervicovaginal microbiotas of healthy women without HPV, HPV infected women and women in several CIN stages.Furthermore, there was the intention to characterize the vaginal cytokine microenvironment, to evaluate the usage of blood count ratios as a prognostic biomarker of systemic inflammation and, finally, to correlate between CVM, HPV genotypes and cytokines.

Study design and participants
This project is a retrospective, observational, multicentric and cross-sectional study, performed at the Genetics Laboratory of ISAMB, Faculdade de Medicina da Universidade de Lisboa, Dr. Joaquim Chaves, clinical analysis laboratory, and Health and Technology Research Center (H&TRC) -Escola Superior de Saúde de Lisboa, from Instituto Politećnico de Lisboa.The microbiome data was already stored in Joaquim Chaves Saude Lab.
By assuming a 10% of significancy level, a statistical power of 90% is achieved.According to Human Papillomavirus and Related Diseases Report (Bruni et al., 2021), 865 women are diagnosed with cervical neoplasia and, by assuming a standard deviation of 0,5, it is possible to extrapolate a need for a minimum sample size number of 68 patients with, at least, 17 patients per group.Therefore, we obtained the minimum quantity of patients, adding 57 more, and so, there was a total sample size of 125 patients.
The study was divided into 4 groups: Control Group 1: participant non-infected with persistent HPV and without cervical lesions; Control Group 2: participant infected with persistent HPV and without cervical lesions; Experimental Group 1: participant infected with persistent HPV and with ASCUS; Experimental Group 2: participant infected with persistent HPV and with LSIL; This study only includes patients with the following characteristics: Age equal or higher than 18 years old and the informed consent from each participant.
The main exclusion criterion was the presence of any known abnormal clinical or laboratory change that, as per the investigator, might interfere with the safety or efficacy assessment or any study procedure.Patient received an informed consent and a demographic questionnaire, via email, that was fulfilled with the aim to adjust the analysis to potential confounding factors.

Microbiota analysis
The cervicovaginal microbiota samples was stored in the Joaquim Chaves Saúde laboratories, namely in ThinPrep ® medium and kept in refrigerator at 0°C.The Thinprep medium stabilizes the cells and microbiota, that were promptly extracted after harvesting.This analysis was performed through a PCR, with specific primers for the V1-V2 regions of the 16sRNA, due to the enhanced sensitivity and specificity (Loṕez-Aladid, 2023).The kit used was specific for 16S analysis, developed by Thermofisher.Afterwards, the samples were sequenced in the Ion GeneStudio S5 System, with Ion-Torrent Software, a system from Thermofisher Scientific.For control proposes, a sample of E. coli was used, in a concentration of 1 ng/mL.The rarefaction level was already part of the Ion-Torrent specifications.The identification of CST was performed by VALENCIA program (France and Ma, 2020), using Silva (www.arb-silva.de,s.d.) and Greengenes database (www.greengenes.secondgenome.com,s.d.).

Identification of other micro-organisms
The determination of other microorganisms was even performed on the same cervicovaginal exudate samples as the microbiota analysis and the HPV determination.
Mycoplasmas spp., essentially Ureaplasma urealyticum, Ureaplasma parvum and Mycoplasma hominis were evaluated by PCR with specific probes, according to the Amplisens/Florocenosis protocol for Mycoplasma.Chlamydia trachomatis and Neisseria gonorrhoeae were evaluated on a COBAS 4800 CT/NG platform from a company based in Switzerland.Finally, testing to identify other etiologic agents of sexually transmitted infections was performed using a CLART STIs clinical assay from Genomica of Spain.

Clinical data
In every of the three phases, clinical data, namely protein Creactive (PCR), neutrophil-lymphocyte ratio (NLR), plateletlymphocyte ratio (PLR) and lymphocyte-monocyte ratio (LMR), was collected from the patient file and was subsequently incorporated into the analysis.

Statistical analysis
Descriptive analysis was used for demographic characteristics and clinical conditions.Data was expressed as means and standard deviations for continuous variables and frequencies and percentages for categorical variables.
In the first phase and second phase, Shapiro-Wilk tests were performed, for continuous variables, with the aim to find out if normality was present.Also, Levene test was applied in continuous variables for checking homogeneity of variances.For categorical variables, Fisher test was implemented to evaluate if there is any statistically significant difference.The outcomes variables were compared between two groups (cases and controls) using statistical tests, such as Student t-test, if normality was present, or the alternative non-parametric Mann-Whitney U test (if the normality assumption is not satisfied).When the comparison is between more than 2 groups, ANOVA is performed, when normality is available, or Kruskal-Wallis test is applied, if normality is not fulfilled.Post-Hoc multiple comparisons, namely Scheffe and Bonferroni, were applied whenever ANOVA test resulted in statistically significant differences.Confounding variables were detected by Mantel-Haenszel statistics and adjusted with ANCOVA (Analysis of Covariance).Pearson and Spearman correlation coefficients and other association measures will be calculated to analyze relationships between continuous or ordinal variables, respectively.Also, linear relationship between variables will be analyzed using linear or logistic regression analysis (according to the nature of dependent variable).
For the third phase, when a comparison between two samples of the same individual was made, a t-student test for paired samples (in the assumption that normality is present) or a Wilcoxon test (in the assumption that normality was absent) was performed.Furthermore, when the different groups with the paired samples were compared, a one way ANOVA was performed, if normality was assumed, whereas a Friedman test was applied, when there was lacking normality.
All statistical analyses were performed using the Statistical Package for the Social Sciences (SPSS), version 24.0, for Windows.An alpha of 0.05 was considered statistically significant, and 95% confidence intervals were reported where appropriate.

Results
First sub-section: cervicovaginal microbiome analysis

Cohort characterization
As previously described, 125 women participated in this study and all their samples were collected during their medical appointments (without any modification of their normal medical routine) and analyzed by Dr. Joaquim Chaves Clinical Analysis Laboratory.Subsequently, according to their cytology results and HPV status, women were divided in 1) Control Group 1 (HPV -; normal cytology); Control Group 2 (HPV +; normal cytology); 3) Experimental Group 1 (ASCUS in cytology); Experimental Group 2 (LSIL in cytology).The whole demographic and clinical characteristics are available in Table 1.

Comparison of demographic and clinical variables between healthy women (with/without HPV) vs. ASCUS vs. LSIL
As observed in Table 1, in terms of demographic features, such as co-infection with other microorganisms and the type of CVM, there is no statistically significant difference between healthy women, ASCUS and LSIL, as well as, regarding clinical information, PCR, platelet-lymphocyte and lymphocyte-monocyte ratios.
Nevertheless, age appears to be a factor that is significantly dissimilar, since younger women had more infections (0,027), and the age differed significantly between healthy HPV negative women and ASCUS patients (Table 1).
Interestingly, neutrophils-lymphocyte ratio seems to have a statistically significant difference (p-value = 0,037) between healthy, ASCUS and LSIL women, maNLRy amid ASCUS and healthy HPV negative women, with a tendency to be increased, as well as among ASCUS and LSIL, with a tendency to be reduced (Table 1).
Since the two main variables that conditionate group formation are HPV status and cytology results, there was the intention to explore if age can influence either one of the variables.From what was appraised, age influences both the susceptibility to acquire HPV (p-value = 0,01) and the cytology result itself (p-value = 0,08).Moreover, the main difference is between LSIL and normal cytology, even behind Bonferroni correction.At last, co-infection with other microorganisms (pvalue = 0,650), especially with Ureaplasma spp.(p-value = 0,165), does not seem to influence cytology results in our population.

Comparison of cervico vaginal microbiota between healthy (with/without HPV) vs. ASCUS vs. LSIL women
Afterwards, we intend to verify if all the variables described above might have an impact in the cervicovaginal microbiota.
Furthermore, it is possible to comprehend that the abovementioned bacteria, essentially L. jensenii and Megaspharae type 1 are the ones that influence cytology result (p-value = 0,008), although there is an absence of impact in HPV status (p-value = 0,323), co-infection with other microorganisms (p-value = 0,328), neither with Ureaplasma spp.(p-value = 328).
Correlation between CVM and clinical and demographic data from healthy women (with/ without HPV) vs. ASCUS vs. LSIL Furthermore, it was assessed whether the CVM correlates in any way with the remaining clinical and demographic data.
Regarding clinical data, an intense and statistically significant positive correlation was observed between HPV status and cytology result (r = 0.631; P-value < 0.001), as expected.In addition, it was similarly found that the bacteria that constitute CVM had direct correlation with cytology results, namely L. iners (r = 0.211; P-value = 0.018), L. jensenii (r = 0.182; P-value = 0.042) and BVAB2 Megasphaera type 1 (r = 0.186; P-value = 0.038).Regarding HPV status, only the latter bacterial species were positively correlated (r = 0.183; P-value = 0.041).The possible interactions between the available laboratory values such as CRP, NLR, PLR and LMR were also evaluated.Therefore, there was an absence of relationships between the CVM and the blood count ratios, but the CRP was directly correlated with the concentration of the G. vaginalis bacterium (r = 0.345; p-value = 0.015).Additionally, a direct correlation was observed between the variable corresponding to the presence of other microorganisms and the concentration of L. jensenii (r = 0.182; P-value = 0.042).

Affinity between CVM bacteria and HPV genotype
Furthermore, the affinity between CVM and HPV genotype was evaluated.Something to highlight is the fact that L. crispatus, the dominant species of CST I, does not show any kind of affinity with HPV.
Finally, we found that other invasive microorganisms, such as Trichomonas vaginalis, presented affinities with a greater number of genotypes, namely HPV-6 (r = 0.

Second sub-section: analysis of the correlation between MCV and the cytokine microenvironment
Comparison of cytokines values between healthy women (with/without HPV) vs. ASCUS vs. LSIL Regarding cytokines, the mean values per group are described in Table 3.
As it can be seen in Table 3, there are no statistically significant changes in any of the cytokines evaluated.
It was then assessed whether several variables, namely age, HPV status, cytology result, presence of other microorganisms, presence of Ureaplasma spp.and type of MCV, significantly influenced cytokine concentrations.It was concluded that most variables did not significantly affect cytokine concentrations (p-Value > 005), apart from the presence of microorganisms which apparently specifically affects the RANTES chemokine (p-Value = 0.033).
Correlation between CVM and cytokines in healthy women (with or without HPV) vs. ASCUS vs. LSIL Furthermore, aimed to understand the relationship between the constituent bacteria of the various types of MCV.
Although most bacteria and cytokines did not show statistically significant correlations, it is possible to highlight the existing negative correlation between the bacteria L. iners and FGF-2 (r= -0.366; p-Value = 0.043), as well as the intense positive correlation between L. gasseri and SCF (r = 0.558; p-Value = 0.001).Additionally, one of the main bacteria belonging to the CST IV-A microbiotic Additionally, correlations were assessed between cytokines and other clinical data, namely cytology result, HPV status, presence of microorganisms, presence of Ureaplasma spp.and type of CVM.Therefore, there were no statistically significant correlations between most variables, except for the variable related to the presence of microorganisms, with the chemokine RANTES (r = 0.401; p-Value = 0.028) and the type of CVM, which correlated directly with the cytokine IL-6 (r = 0.429; p-Value = 0.014).

Discussion
Despite the unequivocal development of the last few years, regarding prognosis, prevention, and treatment of cervical neoplasia, it is still a pathology with a considerable prevalence and associated morbidity and mortality, especially in low-income countries.One of the major recent objectives of the Medical Community is the eradication of this disease, with an underlying symbolism associated with the cure for cancer.
Thus, this study aimed to identify new prognostic biomarkers, to assist the gynecologist physician in his/her clinical decisionmaking process while simultaneously aiming to understand the role of MCV both in the pathophysiology, as well as in its use for prevention and/or therapy in cervical neoplasia.
According to the results obtained, the most prevalent type of CVM in each of the groups is CST I, dominated by L. crispatus, with the following prevalences: 62%, 61.5%, 42.9% and 58.3%, corresponding to healthy women without and with HPV, women with ASCUS and women with LSIL, respectively.This type of CVM is more present in healthy individuals and belonging to European populations (Curty et al., 2020).Comparing with other geographical regions, we found different prevalences, namely a predominance of CST type III in Asia, while in Africa and Australia the dominance is clear for CST type IV (Ravel et al., 2011;Chico et al., 2012;Bradshaw et al., 2013;Kyrgiou et al., 2017;Curty et al., 2020;Onywera et al., 2021).
Furthermore, the impact that the age of patients in the various study groups might have on the interpretation of the results was assessed.We observed that healthy HPV-negative and/or HPVpositive women had a mean age of 44 years, while the remaining women had 38 and 35 years, respectively.When this dissimilarity was further investigated, it was noticed that it was mainly between HPV-negative healthy women and women with ASCUS.Moreover, it was found that this inequality, in terms of age, was significant not only in relation to the cytology result, but also in the susceptibility to acquire HPV.In this case, the interaction is more marked between healthy women and women with LSIL.This relationship agrees with another study, which reports a 50% and 75% risk of acquiring HPV at the mean ages of 20 and 30 years, respectively (Burger et al., 2017).Additionally, another author realized that, as women age, HPV-induced cervical modifications decrease, reaching only 38% above 46 years of age (Kainz et al., 1995), which coincides with the results that we obtained.In addition, it should be noted that many older women have not received the HPV vaccine, contrary to younger women, which may influence both the susceptibility to HPV and the cytology result itself because, if the woman is infected with HPV, it will probably be a genotype with less carcinogenic risk.
Blood count ratios were analyzed with the aim of identifying possible prognostic biomarkers.According to the literature, we can understand that the ratio showed greater prominence and more reliable results is the neutrophil-lymphocyte ratio.Indeed, several studies show that a higher NLR correlates with a worse prognosis, namely through a decrease in both overall survival and disease-free survival (Templeton et al., 2014;Lima et al., 2021).In addition, an increased risk of locoregional metastasis, lymphovascular invasion and tumor size has even been assessed (Wu et al., 2017;Ma et al., 2018).Moreover, it was found that even applying this ratio in CIN, namely regarding its post-excision recurrence, an elevation of this marker, would be related to a decrease in the disease-free interval, which means, an increment in the risk of recurrence (Farzaneh et al., 2019).Given that our study was based on early stages of CIN, it is not possible to substantiate the existing evidence.However, there as a summing of further evidence to this prognostic biomarker, particularly regarding the clinical uncertainty associated with an ASCUS cytology.In this study, it was possible to understand that the NLR's differed significantly between the various groups.When this dissimilarity was further investigated, the conclusion reached was that it essentially resided among women in the ASCUS group, healthy women with negative HPV and women with LSIL.Thus, we understand that an increase in the NLR in a woman with ASCUS relates to a progression to LSIL, while a decrease in the NLR, in the same woman, is associated with a regression to a normal cytology.It therefore helps the gynecologist to decide what the next step should be when faced with a diagnosis of ASCUS.
Furthermore, the associations between the constituent bacteria of CVM and the numerous available parameters were scrutinized.In this way, it was assessed the interaction between CVM and the cytology result itself.
It was possible to observe that both L. jensenii and Megasphaera type 1 were significantly different between the various study groups, essentially between HPV negative healthy women and women with LSIL.The same positive correlation was found between the above bacteria and the cytology result.These results are in line with the existing literature.One of the studies, which associated CVM with cytology result, showed that a CVM of CST type IV, where Megasphaera type 1 is included, increased the risk of developing LSIL by 2-fold, HSIL by 3-fold and CN by 4-fold, respectively (Mitra et al., 2015;Mitra et al., 2016;Kyrgiou et al., 2017).Other studies report similar results, where an CVM is intrinsically associated with a higher risk of developing more advanced CIN and, consequently, CN (Oh et al., 2015;Piyathilake et al., 2016;Łaniewski et al., 2018;Łaniewski et al., 2020;Curty et al., 2020;Norenhag et al., 2020).As for the issue of L. iners, the dominant species of CST III, despite being a Lactobacillus spp., it is classified as a transitional species for microbiotic dysbiosis, given that in some studies, it presents a risk of interaction with other species integrating CST IV and, therefore, increases the risk of developing LSIL and HSIL (Kinney et al., 1998;Kyrgiou et al., 2017).This can be explained by the smaller genome of L. iners, which leads to the production of L-lactic acid, in contrast to the other members of Lactobacillus spp.which produce mainly D-lactic acid.The underlying problem is the different biological actions inherent to L-lactic acid and D-lactic acid, with D-lactic acid having anti-inflammatory characteristics, as opposed to L-lactic acid, which has pro-inflammatory characteristics, with the capacity to reactivate latent viral infections.Complementarily, L. iners also presents the ability to produce inerolysin, similar to vaginolysin, both toxins that allow the formation of pores in the vaginal membrane, with subsequent facilitation of HPV infection and, ultimately, cytological alteration (Rampersaud et al., 2011;Kyrgiou et al., 2017;Petrova et al., 2017;Ragaliauskas et al., 2019;Curty et al., 2020;Pleckaityte, 2020).
As for L. jensenii, it is positively correlated with advanced stages of cervical lesions, as well as increased susceptibility of infection by other microorganisms.Spurbeck et al. demonstrated the antiinflammatory properties of L. jensenii regarding the inhibition of Neisseria gonorrheae infection (Spurbeck and Arvidson, 2010).Another study showed that the presence of L. jensenii would allow achieving cure, in patients with bacterial vaginosis, about 1.67 times compared to the control (Mitchell et al., 2012).Therefore, our results are dissonant with the literature.Thus, our database was explored with the purpose to find an explanation, which justifies this dissimilarity.The most plausible explanation is that only 1 of 125 women participating in the study had a CVM with identification of L. jensenii.Therefore, the statistical power is too low to draw conclusions.
Furthermore, a direct correlation was found between the bacteria L. crispatus and L. gasseri, which was already expected, since although L. crispatus is the dominant species in most microbiotas.It is possible to observe a symbiotic relationship with other Lactobacillus species in order to maintain the homeostasis of the cervicovaginal microenvironment.Furthermore, both bacteria possess the ability to inhibit cell proliferation and induce apoptosis, essentially in tumor cells, making them crucial for protection against both infectious and oncological pathologies (Motevaseli et al., 2013).However, L. crispatus show negative correlations with species belonging to CST IV, compatible with microbiotic dysbiosis, namely A. vaginae, G. vaginalis and BVAB2 Megasphaera type 1.These associations are plausible given the very functioning of L. crispatus which possesses several strategies to reduce colonization by other more deleterious species.Firstly, it has the capacity to decrease the pH of the medium, converting itself into a medium with low pH, through the production of D-lactic acid, by a mechanism whose substrate is glycogen, deposited in the vaginal walls, resulting from the action of oestrogen.It can then also produce hydrogen peroxide and bacteriocins.Taken together, they prevent the growth of other harmful bacterial species (Kyrgiou et al., 2017;Tachedjian et al., 2017;Curty et al., 2020;van de Wijgert and Verwijs, 2020).Therefore, our results are in agreement with the literature.Furthermore, L. gasseri also expresses the same negative correlations with A. vaginae, G. vaginalis and BVAB2 Megasphaera type 1, but with less intensity than L. crispatus.In turn, something not yet described in the literature, was the positive interaction between BVAB2 Megasphaera type 1 with A. vaginae.Apparently, the deleterious species present a symbiosis among themselves and, when there is a significant depletion of Lactobacillus spp.and allow their emergence and growth, these species assist the translocation of other species, with similar characteristics, thus increasing their dominance in the medium and, leading to the conversion to a CST type IV microbiota.Looking at other data, it may be possible to explain the resistance to probiotic colonization in the cervicovaginal region or the resistance provided to the treatment of bacterial vaginosis (Mitchell et al., 2012;Srinivan et al., 2015;Kyrgiou et al., 2017;Łaniewski et al., 2020;Curty et al., 2020).
Then, it was checked, in a more detailed way, if CVM had affinity for the various genotypes of HPV causing infection.Thus, it was shown that L. crispatus does not show any affinity for any HPV genotype, thus constituting the most resistant type of CST against HPV infection.However, a CVM dominated by L. gasseri had a higher risk of being infected by an HPV of genotype 51, whereas a CST type V MCV, dominated by L. jensenii, had a higher risk of being infected by HPV-56.In turn, L. iners, being a transitional species, shows affinity for a higher number of genotypes, namely HPV-18, HPV-53 and HPV-59.Therefore, it was possible to conclude that all Lactobacillus spp. that show some kind of affinity, only correlate with the most common genotypes in normal cytology and low-grade lesions (So et al., 2019).However, it was observed that the remaining species present in CST IV, the presence of infection by other microorganisms and even colonization by Ureaplasma spp.manifested a greater number of affinities with various HPV genotypes.However, they were correlated to genotypes more present in HSIL, increasing the risk of developing CN (So et al., 2019).Of note, this interaction between CVM and affinity for certain HPV genotypes had not been described in the literature.
Finally, the focus was on a topic little explored in the literature, namely the characterization of the cervicovaginal microenvironment regarding the cytokines present.Therefore, when analyzing the literature on the topic, we found that one study reported that a state of dysbiosis associated essentially with CST IV expressed an increase in pro-inflammatory cytokines, especially TNF-a, IL-6, IL-8, RANTES, IFG-g and GMC-SF, while simultaneously manifesting a decrease in the concentration of cytokines with anti-inflammatory characteristics, such as IP-10, compared to women with so-called "healthy" CST (Hedges et al., 2006;Novak et al., 2007;Anderson et al., 2011;Anahtar et al., 2015;Kyrgiou et al., 2017;Łaniewski et al., 2018).Another cross-sectional study evaluated women with CN or dysplasia and healthy women with or without HPV infection and assessed the cervicovaginal microenvironment that had multiple proinflammatory cytokines (IL-36g, IP-10, MIP-1b and RANTES), haematopoietic (FLT3 ligand), growth factors and angiogenic factors (FGF-2, SCF, TRAIL), hormones (prolactin), chemokines (MIF and TNF-a) and osteopontin.All this microenvironment was found to correlate positively with vaginal inflammation, but also with microbiotic dysbiosis (Łaniewski et al., 2019; Łaniewski et al., 2020).Unfortunately, given the tiny sample size where cytokines were assessed (n= 39), it was not possible to obtain the vast majority of the aforementioned correlations.However, we were able to perceive the interaction between infection by other microorganisms and the chemokine RANTES.This chemokine, also called CCL5, is produced by platelets, and has been shown to play an important role in cancer progression, namely in invasion, metastasis, angiogenesis, and immune cell infiltration.It appears to function as a crucial regulator, allowing two-way communication between the microenvironment, modulated by the microbiota, and tumor cells.Furthermore, it promotes macrophage recruitment, collagen deposition and tumor relapse (Walens et al., 2019).Consequently, it is present from the early stages of the tumor and therefore could be an important chemokine for understanding tumor development.

Conclusion
The main limitations are the small number of confounding factors that may influence the results, as well as the small number of valid samples for cytokine analysis, which greatly reduced the information that we could have obtained.Furthermore, all the women only performed cytology, since the majority of cervical lesions were in the early stages, without indication for cervical biopsy, therefore it is a limitation itself, since there is not any histological confirmation of the abovementioned lesions.At last, since this is a retrospective project, there was an absence of information regarding the clinical status of the woman (such as immunosuppression), limiting the extrapolation of knowledge, since it can influence cervicovaginal microbiota composition and cervical neoplasia development.
In short, this study might contribute for another step in the development of a more personalized medicine, with more effective and efficient results and with better care for patients.Additionally, it has brought us closer to one of the most pressing challenges, namely the eradication of cervical cancer.

TABLE 1
Demographic and clinical characteristics of first phase analysis.
1Application of Kruskal-Wallis test. 2 Application of Chi-square test.All the bold values, in the p-value column, represent statistically significant results with p-value < 0,05.NA, not applicable.

TABLE 2
Cervicovaginal microbiota of first phase analysis.
1Application of Kruskal-Wallis test.All the bold values, in the p-value column, represent statistically significant results with p-value < 0,05.

TABLE 3
Cytokines analysis of second phase.
1 Kruskal-Wallis test.All the values in bold are statistically significant results with a p-value < 0,05.