- 1Department of Medical and Surgical Sciences (DIMEC), Alma Mater Studiorum - University of Bologna, Bologna, Italy
- 2Institute of Biomedical Technologies, National Research Council, Segrate, Italy
- 3National Biodiversity Future Center S.c.a.r.l., Palermo, Italy
- 4Section of Microbiology, Department of Medical and Surgical Sciences (DIMEC), Alma Mater Studiorum - University of Bologna, Bologna, Italy
- 5International PhD College, Collegio Superiore of Alma Mater Studiorum, University of Bologna, Bologna, Italy
- 6Department of Medical Sciences, University of Ferrara, Ferrara, Italy
- 7Microbiology Unit, IRCCS Azienda Ospedaliero-Universitaria di Bologna, Bologna, Italy
- 8Department of Animal Science, Food and Nutrition (DIANA), Università Cattolica Del Sacro Cuore, Piacenza, Italy
Introduction: The vaginal microbiota (VMB), predominantly shaped by Lactobacillus species, plays a crucial role in maintaining vaginal health and preventing infections. However, its delicate balance is increasingly challenged by the widespread use of antibiotics and the consequent rise in antibiotic resistance genes (ARGs). While the VMB has been recognized as a potential reservoir for ARGs, limited research has explored how microbial composition, antibiotic exposure, and individual behavioral factors converge to shape the vaginal resistome.
Materials and methods: In this cross-sectional study, vaginal swabs were collected from 105 reproductive-age Caucasian women. The VMB composition was characterized and classified into Community State Types (CSTs), by sequencing the hypervariable V3-V4 regions of the bacterial 16S rRNA gene. In order to highlight common patterns of abundance among taxa, a co-abundance groups (CAGs) analysis was performed. We assessed the distribution of 14 ARGs conferring resistance to macrolides, tetracyclines, beta-lactams, and quinolones along with two associated transposons by means of PCR. An overall composite ARGs score was also calculated. Participants completed detailed questionnaires assessing demographics and behavioral factors, with a particular focus on both health- and antibiotic-related behaviors. Statistical analyses examined associations between ARG prevalence, vaginal microbiome composition and relevant exposures.
Results: CSTs I and III were the most prevalent, with the most frequently detected ARGs being erm(F), tet(M), erm(B), erm(A), and tet(W), each present in over 65% of participants. ARG presence was positively correlated with a higher vaginal microbiome alpha-diversity. Moreover, BV-associated bacterial taxa showed strong associations with ARGs, while Lactobacillus species generally exhibited negative correlations. Smoking, a higher body mass index (BMI), presence of Candida spp. and a history of antibiotic use were significantly associated with increased ARG prevalence, whereas oral contraceptive use and higher diet quality scores were negatively associated. Correlating together the relative abundances of the microbial CAGs and the presence/absence of specific ARGs, we found a positive association between several genes related to macrolide and tetracycline resistance and the Gardnerella-Prevotella CAG. Additionally, the Gardnerella-Prevotella, and the Streptococcus CAGs were positively correlated to the total ARG score, whereas the L. crispatus/jenesenii/gasseri CAG was negatively correlated.
Conclusions: These findings underscore the role of the VMB as a dynamic reservoir of ARGs and highlight the influence of individual lifestyle and antibiotic-related behaviors on ARG dissemination in the vaginal niche. This supports the need for integrated public health strategies that combine antibiotic stewardship with targeted lifestyle and behavioral interventions, as well as the development of individualized therapeutic approaches.
1 Introduction
The vaginal microbiota (VMB) is typically characterized by a low bacterial diversity, in stark contrast to other mucosal sites such as the gut, where microbial diversity is high. In healthy, reproductive-aged women, the VMB comprises a variety of aerobic and anaerobic bacterial genera and species, with the genus Lactobacillus generally dominating. A Lactobacillus-dominated VMB supports vaginal health by preventing the colonization and overgrowth of pathogenic microorganisms, maintaining homeostasis and eubiosis, and has been associated with a lower risk of infections (Wang et al., 2024; Greenbaum et al., 2019; Ceccarani et al., 2019).
Next-generation sequencing of the hypervariable regions of the 16S rRNA gene has enabled detailed analyses of the composition of the VMB, leading to its classification into five main Community State Types (CSTs). Four of these are dominated by individual Lactobacillus species (CST I - L. crispatus, CST II - L. gasseri, CST III - L. iners, CST V - L. jensenii), while CST IV is characterized by a heterogeneous assemblage of strict and facultative anaerobes, including Gardnerella, Atopobium, Mobiluncus, and Prevotella (Borgogna et al., 2020; France et al., 2022).
Over a woman’s lifespan, the VMB, which plays a pivotal role in women’s health, can undergo significant changes in response to both local and systemic factors, potentially leading to dysbiosis. In literature, a long but non-exhaustive list of risk factors associated with vaginal dysbiosis has been reported. These include age, body mass index (BMI), hormonal status, diet, ethnicity, smoking, lifestyle, antibiotic use, and the presence of urogenital infections. Although essential for treating bacterial infections, antibiotics can inadvertently disrupt the delicate vaginal microbial ecosystem, resulting in negative health consequences (Abou Chacra et al., 2022; Ahrens et al., 2020; Morsli et al., 2024; Das et al., 2023).
The global increase in antibiotic use, often beyond clinical indications, has sparked concern regarding its long-term impact on human health. In particular, the overuse and misuse of antibiotics has contributed to the emergence of antibiotic-resistant bacteria and the spread of antibiotic resistance (AR), posing a significant threat to public health. Tetracyclines are the most common antibiotics in agricultural soil and food production, due to their widespread usage and strong persistence (Wang et al., 2024). Moreover, tetracyclines together with other antibiotic classes, such as macrolides and quinolones, which disrupt DNA and protein synthesis, have displayed in vitro activity against different pathogens, thus becoming the recommended drugs for clinical treatment in different settings (Georgakopoulou et al., 2024). Selective pressure exerted by environmental antibiotics facilitates the spread of AR through horizontal gene transfer, predominantly via mobile genetic elements such as plasmids and transposons (Sirichoat et al., 2020; Flórez et al., 2016).
In particular, the link between antibiotic use and the emergence of resistance in the VMB is of special concern, since the VMB has indeed been identified as a significant reservoir of AR determinants (Jeters et al., 2009; Severgnini et al., 2021). Moreover, the increasing prevalence of antibiotic resistance genes (ARGs) in commensal bacteria is problematic, as these organisms can serve as reservoirs, transferring resistance determinants to pathogenic species (Melkumyan et al., 2015).
Individuality and social aspects play a role in antibiotic adherence behavior. About the first, both the fear of illness and the lack of knowledge about how antibiotics function (Nortey et al., 2023) lead to the misuse of antibiotics. In terms of sociality, the family also represents an environment capable of affecting the health behaviors of individuals and situations of unhealthy practices (e.g., drug consumption) (Bagley et al., 2022).
A recent study highlighted the impact of family experiences and individual attitudes on antibiotic-related behaviors (Castellano et al., 2024), underscoring the need to explore how these psychosocial elements influence the presence of ARGs in the vaginal environment.
In this context, it’s important to investigate the distribution of resistance genes in the VMB, the relationship between antibiotic consumption and the presence of ARGs in the vaginal milieu, and how individual practices, awareness, and knowledge regarding antibiotic use may influence the VMB and the prevalence of ARGs.
In a previous paper (Severgnini et al., 2021), we assessed the presence of only 4 resistance markers (i.e., erm(B) and erm(F) conferring resistance to macrolides; tet(W) and tet(M) to tetracyclines) in the vaginal environment of women at different gestational ages.
With the aim of expanding and further exploring the dynamics of the vaginal environment, in this cross-sectional study, we assessed the distribution of a wider panel of ARGs (n=14) conferring resistance to macrolides, tetracyclines, beta-lactams and quinolones, along with two associated transposons, across vaginal CSTs in non-pregnant reproductive-age women. In particular, the presence of the following determinants of resistance was evaluated: erm(A), erm(B), erm(F), tet(M), tet(M)-Tn916, tet(O), tet(W), tet(Q), tet(Q)-rteA, blaOXA-2, blaTEM, blaZ, blaSHV, blaCTX-M, qnrA and qepA. As an innovative aspect, we also examined potential factors (exposures) associated with the presence of these genes and the overall composite ARGs score, with a particular focus on individual knowledge, attitudes, and habits regarding antibiotic use, employing a multivariate approach. The detection of ARGs was subsequently correlated with the bacterial composition of the vaginal microbiome.
2 Materials and methods
2.1 Study setting and population
From November 2023 to April 2024, 123 Caucasian women of reproductive age were enrolled in the study. The participants were volunteer students attending degree courses at the University of Bologna, Italy. At enrolment, the exclusion criteria were: (i) antibiotic use in the month prior to sampling; (ii) use of vaginal douches or topical agents in the last two weeks; (iii) age < 18 years; (iv) pregnancy; (v) menstruating at the time of sampling; (vi) HIV infection; (vii) presence of chronic conditions (e.g., diabetes, autoimmune diseases, malignancies); (viii) drug addiction or heavy smoking (>15 cigarettes/day). Moreover, women with urogenital infections due to sexually transmitted pathogens (i.e., Chlamydia trachomatis, Neisseria gonorrhoeae, Trichomonas vaginalis, Mycoplasma genitalium) were further excluded when diagnostic test results were available.
Several demographic (e.g., BMI, educational level, marital status), health-related (e.g., presence of Candida spp., smoking habits) and antibiotic-related behavioral (e.g., individual practice and awareness about antibiotics) data were collected (see specific paragraphs on exposures).
After obtaining written informed consent from all participants, they underwent two self-collected vaginal samplings at the same time point. The first swab (E-swab, Copan, Brescia, Italy) was used for diagnostic tests to exclude the presence of sexually-transmitted infections (STIs) by NAAT (Nucleic Acid Amplification Test, Alinity m STI Assay, Abbott Molecular Inc, Des Plaines, IL, USA). The second was collected with a sterile cotton bud, re-suspended in 1 mL of sterile saline, and stored at -80 °C until use.
The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board (IRB) of the University of Bologna, protocol number 0122421.
2.2 Amplification of the various ARGs
Nucleic acids were extracted from vaginal swabs by means of the DNeasy Blood & Tissue Kit (QIAGEN GmbH, Hilden, Germany) according to the manufacturer’s instructions. From the DNA eluate, each sample was screened for the presence of selected resistance genes and specific transposons that confer resistance to different classes of antibiotics by means of end-point PCR assays. The following genes were investigated: erm(A), erm(B) and erm(F) conferring resistance to macrolides; tet(M), tet(O), tet(W), and tet(Q), conferring resistance to tetracyclines; blaOXA-2, blaTEM, blaZ, blaSHV and blaCTX-M associated with β-lactam resistance; and qnrA and qepA linked to quinolone resistance. Additionally, the presence of two transposons associated with AR, namely tet(M)-Tn916 and tet(Q)-rteA, was also assessed. Detailed information on primers used and PCR conditions are displayed in the Supplementary Materials (Supplementary Table S1). These results were dichotomized as 0=absence of antibiotic resistance and 1=presence of antibiotic resistance.
An overall composite score for the ARGs was developed by summing the presence/absence of all the genes detected in at least one sample. The total possible score ranged from 0 to 10, with higher scores indicating greater antibiotic resistance. The score was also categorized into tertiles based on its distribution: low (≤3), medium (4-5), and high (≥6).
2.3 Vaginal microbiome molecular profiling
The hypervariable V3-V4 regions of the bacterial 16S rRNA gene were amplified from genomic DNA extracted from vaginal swabs. The PCR conditions and primer sequences were obtained from the Illumina 16S Sample Preparation Guide (https://support.illumina.com/documents/documentation/chemistry_documentation/16s/16s-metagenomic-library-prep-guide-15044223-b.pdf) (Illumina, San Diego, CA, USA), with the primers originally described in Klindworth et al., 2013. Final indexed libraries were prepared by equimolar (4 nmol/L) pooling, denaturation, and dilution to 6 pmol/L before loading on a MiSeq flow cell (Illumina) for a 2 × 300 bp paired-end run. 16S rRNA sequences were processed according to the same methods as described in Severgnini et al., 2021, using the latest release of the SILVA database (v. 138, Quast et al., 2013) for classification.
Depth of sequencing was set to the lowest sequenced sample (n=6,627 reads), in order to compensate for the sequencing unevenness of the samples and to provide a consistent minimum amount for the downstream analysis, carried out through the QIIME pipeline (version 1.9.0, Caporaso et al., 2010). Alpha-diversity was estimated according to Chao1, Shannon Index, Observed Species, and the Faith’s phylogenetic tree diversity metrics (“PD whole tree”), whereas beta-diversity was evaluated using unweighted and weighted UniFrac distances (Lozupone et al., 2011). Since all alpha-diversity indices performed similarly, we focused on the PD whole tree for the univariate and multivariate analyses (see below).
Lactobacillus species-level characterization was performed as in Severgnini et al., 2022 by BLAST-aligning all reads belonging to the Lactobacillaceae family to a custom reference database made up of all available reference sequences in the NIH-NCBI database (https://ftp.ncbi.nlm.nih.gov/genomes/GENOME_REPORTS/prokaryotes.txt) of 17 species commonly found in the vaginal environment, including a total of 3,392 genomes. In case of multiple matches with the same confidence, the taxonomy was reset to “Unclassified Lactobacillus” at the genus level.
Community-state types (CST) of the vaginal microbial communities were determined from the taxonomic profiles using VALENCIA, a nearest centroid-based tool that classifies samples into 5 major CST according to the similarity to a set of about 13,000 reference microbial profiles (France et al, 2020).
In order to highlight common patterns of abundance among taxa, a co-abundance groups (CAGs) analysis was performed as previously described (Claesson et al, 2012). Briefly, only the bacterial taxa present at >0.5% of abundance in at least 10% of the samples (n≥11) were selected, in order to exclude minor and transient contributors of the gut microbiota. This resulted in a subset of 21 taxa on the whole dataset. The co-abundance between each pair of taxa was evaluated by calculating Spearman’s correlation coefficient and displayed in a heatmap, hierarchically clustered using Pearson’s correlation metric and average linkage. Only associations having a Benjamini-Hochberg adjusted p<0.05 were used to build the hierarchical clustering. Permutational multivariate analysis of variance (P-MANOVA) was used to assess that the computed CAGs were significantly different from each other, using p<0.05 as a threshold.
2.4 Exposures
For the present analysis, exposures were selected based on theoretical knowledge because they are crucial for contextualizing the results, as they could potentially influence both antibiotic use and health outcomes (Morsli et al., 2024). These factors included the following demographic, health- and antibiotic-related behavioral data of the participants: age (continuous), educational level (categorized as high school vs university degree or higher), marital status (married/cohabiting vs single/separated), body mass index (BMI, calculated as weight in kg divided by height in cm squared), contraceptive use, smoking habits (non-smokers vs smokers), presence/absence of Candida spp., adherence to a healthy diet such as the Mediterranean diet (measured calculating the MEDILITE score, Sofi et al., 2017). In addition, we investigated individual practice and awareness about antibiotics and AR, and participants’ family of origin approach and behavior toward antibiotics, through ad hoc questions, to identify possible misuse of these drugs. An adapted version of the “Antibiotic Resistance: Multi-Country Public Awareness Survey,” developed by the World Health Organization (WHO, 2015), was employed. This revised survey consisted of 25 items organized into several sections. The first section contained two questions regarding the amount of antibiotics taken in the last year and month, along with four questions focused on best practices for antibiotic use where participants’ individual scores were derived from the sum of their answers. The second section included nine statements aimed at assessing participants’ understanding of antibiotic resistance, with responses recorded as True or False. Individual scores were calculated by summing the responses. The final section comprised ten statements evaluating participants’ beliefs and attitudes towards antibiotic use and resistance, utilizing a 5-point Likert scale (1 = strongly disagree to 5 = strongly agree). The total scores from these items were used as individual scores to gauge “awareness,” which refers to the level of self-awareness and understanding regarding strategies to combat AR. Eight ad hoc questions were included to investigate participants’ family background and attitudes towards antibiotics, with the aim to identify family models which used antibiotics appropriately (e.g., avoiding self-prescribing and adhering to prescribed dosages). In order to better explore the variable related to family antibiotic improper use, we categorized the total score in tertiles based on the distribution as low (≤10), medium (11–17) and high (≥18).
2.5 Statistical analysis
Sample characteristics were summarized using mean and standard deviation (SD) for continuous variables and counts and percentages for categorical variables.
Statistical evaluation of alpha-diversity indices for the microbiota profiles was performed by non-parametric Monte Carlo-based tests. Beta-diversity differences were assessed by a permutation test with pseudo F-ratios using the “adonis” function from R package “vegan” (version 2.0-10, https://cran.r-project.org/package=vegan). Relative abundances were analyzed first setting all “0” to a value of 10–6 in order to compensate for missing data due to the normalized ASV table (where the “0” can be due to insufficient depth of coverage for assessing the ASV presence) and, then, performing a centered log-ratio (CLR) transformation in order to deal with the compositional nature of the data; finally, a Kruskal-Wallis test, followed by a Dunn post hoc evaluation was performed. Correlation between presence of specific ARGs and bacterial taxa abundances was performed by point-biserial correlation, whereas correlation between the total ARG score and bacterial taxa abundances was based on the Spearman’s correlation coefficient (with a p-value ≤ 0.05 considered as significant).
Two regression models were performed to identify factors associated with antibiotic-resistance genes. Binary logistic regression models were constructed to calculate odds ratios (ORs) and 95% confidence intervals (CIs) to estimate the association between exposures of interest and the presence of each of the antibiotic resistance genes.
Linear regression models were used to calculate the exponentiated β-coefficients with 95% CIs to examine the associations between exposures and the continuous antibiotic-resistance gene score.
We first performed univariate analyses including all the following variables in the model separately: age, educational level, marital status, smoking habits, BMI, presence of Candida, family history of inappropriate antibiotic use, antibiotics taken in the past year, past compliance with antibiotic use, awareness of antibiotic resistance, adherence to MEDILITE, and biodiversity index (PD whole tree). Multivariate models were performed, including factors associated with a p-value ≤ 0.10 in univariate analyses. All analyses were performed with STATA (version 15.0, StataCorp LP, College Station, Texas, USA) and Matlab (v 2008b, Natick, MA, USA) and the functions from the Fathom Toolbox (Jones, 2017). Two-sided p-values of 0.05 and 0.10 were considered statistically significant or borderline significant, respectively.
2.6 Data availability
Raw sequencing data for this project are available in NCBI Short-Read Archive (SRA) under accession number PRJNA1188525 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1188525).
3 Results
3.1 Study population and exposures
A total of 123 participants were initially recruited. After excluding individuals with missing vaginal microbiota data (n = 7) and/or incomplete antibiotic behavior questionnaires (n = 14), 105 women were included in the final analysis. The mean age of participants was 21.6 ± 2.6 years (range: 19–30), and the mean BMI was 21.5 ± 2.8 kg/m² (range: 16.2–30.1).
Vaginal samples were categorized into the five main CSTs according to their microbial composition: CST I (41.9%, n = 44), CST II (5.7%, n = 6), CST III (30.5%, n = 32), CST IV (17.1%, n = 18), and CST V (4.8%, n = 5).
Many participants reported some level of inappropriate antibiotic use within their family of origin, with 32.4% indicating a high tendency toward misuse. Regarding individual antibiotic consumption, 42.9% of participants reported not having taken antibiotics in the previous year.
Overall, participants demonstrated high awareness of antibiotic resistance (mean score: 42 ± 4), along with strong adherence to recommended antibiotic use practices (mean score: 9 ± 2). Baseline characteristics of the sample, including behavioral factors and antibiotic-related variables, are presented in Table 1.

Table 1. Participants’ baseline characteristics for sociodemographic, behavioral, and antibiotic-related variables (n=105).
3.2 Detection of ARGs
Regarding ARGs distribution, erm(F) was the most detected one, with 72.2% of positivity (80 cases), followed by: tet(M) (74.3%, 78 cases), erm(B) (68.6%,72 cases), erm(A) (66.7%, 70 cases), tet(W) (65.7%, 69 cases), blaZ (36.2%, 38 cases), tet(Q) (35.2%, 37 cases), and blaTEM (16.2%, 17 cases). Furthermore, the less frequently detected genes were blaOXA-2 and tet(O), both with 2.9% of positivity (3 cases), qnrA with 1.9% (2 cases) and qepA with only 1% (1 case). No samples tested positive for blaSHV or blaCTX-M.
Positivity for specific ARGs was frequently associated with the presence of the corresponding mobile genetic elements. Specifically, tet(Q)-rteA was identified in 100% of the tet(Q)-positive samples (37/37), and tet(M)-Tn916 was detected in 96.2% of the tet(M)-positive samples (75/78).
3.3 Factors associated with ARGs
We initially performed univariate analyses to explore potential associations between selected factors and the presence of ARGs (data not shown). Variables with a p-value ≤ 0.10 were subsequently included in multivariate models to assess their independent association with the outcomes of interest.
Several associations remained statistically significant in the multivariate analyses (Table 2). In detail, we found that the microbial alpha-diversity, estimated by the Faith’s phylogenetic distance metric (PD whole tree), was positively associated with tet(M), tet(O), tet(Q), erm(A), and blaZ, with odds ratios ranging from 1.57 to 2.22. Faith’s PD was also strongly associated with the total ARG score (Exp(β-coefficient) 1.82, 95%CI 1.53-2.17), whereas borderline significance was observed with blaTEM (p-value ≤ 0.10).
Contraceptive use reduced the odds of tet(M) by 70% compared to non-use. Smoking increased the odds of erm(B) by 3.47 times and the Exp(β-coefficient) of the total score increased by 2.09 times. A higher BMI was positively associated with erm(A) and a higher total ARG score, with a borderline positive association with tet(O). Diet quality, as assessed by the MEDILITE score, was inversely associated with tet(O) detection, suggesting a potential protective effect of adherence to a Mediterranean diet. The presence of Candida spp. was significantly associated with a higher total ARG score (OR = 3.85).
Regarding antibiotic behavior, a higher degree of inappropriate antibiotic use in the family was associated with a slight increase in the odds of erm(A) and a borderline association with blaOXA-2. Antibiotic use in the past year was positively associated with tet(M). Past compliance of antibiotic usage mainly reduced the odds of tet(M) and blaZ. Moreover, a borderline inverse association with erm(B) and the total score was also noticed (with odds ratios ranging from 0.47 to 0.80). Awareness about antibiotic resistance decreased the odds of tet(Q) by 14%.
3.4 Analysis of the vaginal microbiota associated with ARGs
The multivariate analysis indicated that the alpha-diversity of the vaginal microbiota (according to the PD whole tree metric) was positively associated with the presence of several ARGs. At the same time, the microbial profiles between samples positive or negative to several ARGs were found to be statistically different. In particular, we confirmed the different microbial composition of samples positive/negative for erm(A), erm(B), tet(M), tet(O), tet(Q), and blaZ (unweighted UniFrac distance, p ≤ 0.008, adonis test of pseudo F-ratios). Moreover, also other ARGs (i.e.: erm(F) and tet(W)) had different microbial profiles (p=0.001 and p=0.002, respectively, unweighted UniFrac), whereas blaTEM, borderline associated to biodiversity, showed only a trend towards separation (p=0.191 unweighted UniFrac). Finally, microbial profiles for blaOXA-2, qnrA and qepA were not significantly different, due to the very low fraction of samples positives for these ARGs (n=3, n=2 and n=1, respectively) (Supplementary Figure S1 in the Supplementary Materials).
Interestingly, significant changes in both alpha- and beta-diversity were found when analyzing the vaginal microbiota and the number of ARGs per sample. In fact, a clear trend towards an increase in alpha-diversity with increasing ARGs was evident, both considering the total ARG score and its categorization into tertiles. Moreover, the beta-diversity analysis showed a significant separation among the samples according to the ARG score tertiles (p<0.008 and p<0.024 for all pairwise comparisons, unweighted and weighted UniFrac distances, respectively) and also a pattern along principal coordinate 1 moving from the right to the left of the x-axis with increasing total ARG scores (Figure 1).

Figure 1. Boxplots of the alpha-diversity estimation according to Faith’s phylogenetic diversity metric with the samples grouped according to the total ARG score (A) or the total ARG score as categorized in tertiles (B). Blue dashed lines represent the trend of the distribution medians; boxplots report the average (blue solid line) and median (red solid line) and individual samples as dots superimposed to the box-and-whiskers plots. Principal coordinate analysis (PCoA) plots based on the unweighted UniFrac distance among samples. Points represent single samples, colored, in each plot, according to the total ARG score as categorized in tertiles (C) or the total ARG score (D); marker shapes for the single samples are represented according to the sample CST; centroids represent the average coordinate of the samples within the same category and ellipses are the SEM-based confidence intervals; for each plot, the first and second coordinate are represented.
Considering the bacterial taxa involved in the vaginal microbiota, the abundance of some genera such as Prevotella, Porphyromonas, Dialister, Finegoldia, and Peptoniphilus were significantly (p<0.05, Kruskal-Wallis and Dunn’s post hoc pairwise tests) increased in both the “medium” and the “high” categorized ARG score tertiles as compared to “low”, whereas Gardnerella, Streptococcus, Megasphaera, Haemophilus, Anaerococcus, and Fusobacterium increased in “high” score samples only (Supplementary Figure S2 in the Supplementary Materials). Species from the Lactobacillus genus did not show a significant decrease with increasing total ARG scores but, as a whole, they accounted for about 90% of the relative abundance in “low” samples, while shifting to about 70% in “medium” and to about 55% in “high” ones (Figure 2).

Figure 2. Barplots representing the average relative abundance for the main bacterial taxa (average rel. ab >0.5% in at least one experimental group) over the samples grouped by the categorized total ARG score. The score was categorized into tertiles based on its distribution: low (≤3 ARGs), medium (4–5 ARGs), and high (≥6 ARGs).
3.5 Correlation between bacterial taxa and ARGs
Co-abundance analysis highlighted the presence of 4 groups (CAGs) with a similar pattern of abundance of their members over all the samples (Supplementary Figure S3 in the Supplementary Materials), which were labelled according to the most representative taxon/taxa (Supplementary Table S2 in the Supplementary Materials). Notably, while L. crispatus, L. jensenii and L. gasseri all belonged to the same CAG, L. iners stood apart together with Ureaplasma; the two other CAGs were those of taxa more frequently associated to vaginal dysbiosis, comprising Gardnerella, Prevotella, Atopobium, Dialister and Megasphaera in one and Streptococcus, Porphyromonas, Anaerococcus, and Finegoldia in the other.
Correlating together the relative abundances of the microbial CAGs and the presence/absence of specific ARGs over all the samples (Figure 3) highlighted some interesting features, such as the positive association between several genes related to resistance against macrolides (i.e.: erm(A), erm(F)) and tetracycline (i.e.: tet(M), tet(Q)) and the Gardnerella-Prevotella CAG, whereas the Streptococcus CAG showed some weak positive correlations, in particular, to tet(W) and blaZ. On the other hand, the L. crispatus/jenesenii/gasseri CAG was negatively associated with all ARGs, particularly to tet(M) and tet(Q), whereas the L. iners CAG was negatively correlated to blaZ resistance. Notably, a positive correlation to quite all the ARGs was observed for the “Other” CAG, which groups together all the low-abundance taxa (i.e.: relative abundance >0.5% in <10% of the samples). At the same time, the Gardnerella-Prevotella, the Streptococcus and the “Other” CAGs were positively correlated to the total ARG score, whereas the L. crispatus/jenesenii/gasseri CAG was negatively correlated.

Figure 3. Heatmap showing the correlation coefficients between the relative abundances of the bacterial co-abundance groups (left) and taxa (right) and the total ARGs score (top) and the presence/absence of specific ARGs (bottom). The point-biserial and the Spearman’s correlation coefficients were used for single ARGs and total ARG score, respectively. ARGs and taxa are organized clustering together patterns of similar correlation (as determined by hierarchically clustering according to Pearson’s correlation and average linkage). Color intensity is proportional to the absolute magnitude of correlation. Positive correlations are reported in red, whereas negative ones are in blue. White stars in the heatmap indicate the statistically significant correlations (p ≤ 0.05).
When evaluating single taxa, we found that, among the Gardnerella-Prevotella CAG, Prevotella and Dialister showed the higher number of ARGs (erm(A), erm(B), erm(F), tet(M), tet(Q)), whereas Atopobium, Gardnerella, and Megasphaera all showed a correlation to tet(Q); all these genera were positively correlated to the total ARG score. At the same time, in the Streptococcus CAG, several positive correlations were found for Finegoldia (i.e.: erm(A), erm(F), tet(M), tet(Q)), Porphyromonas and Anaerococcus (i.e.: erm(A), erm(B), erm(F), tet(M), tet(Q), tet(W), tet(O), blaZ); again, all of them were directly correlated to the total ARG score. Conversely, in the L. crispatus/jenesenii/gasseri CAG, L. crispatus was negatively associated with all ARGs, particularly to tet(M) and tet(Q), while the other lactobacilli in this CAG showed a similar behavior, except with a positive association between L. gasseri and erm(A), and between unclassified members of the Lactobacillus genus and qnrA. Finally, among the bacteria comprised in the “Other” CAG, in particular, Sneathia showed a positive correlation to qepA and tet(Q), whereas Escherichia-Shigella and Haemophilus were correlated to blaTEM and blaZ genes.
4 Discussion
Antibiotic resistance represents one of the greatest challenges facing modern healthcare, driven by the widespread use of antibiotics and the consequent dissemination of ARGs. These genes enable bacteria to survive antibiotic treatments, potentially compromising the effectiveness of antibiotic therapies. Investigating ARGs in the vaginal microbiota is therefore of significant relevance for both women’s and public health.
Bacteria in the vaginal flora can develop and transmit drug resistance determinants, particularly against tetracyclines, macrolides and beta-lactams, which are the most commonly used classes of antibiotics. Additionally, resistance to quinolones has recently emerged as a notable and potentially concerning trend.
In this study, we assessed the distribution of selected ARGs conferring resistance to tetracyclines, macrolides, beta-lactams and quinolones within CSTs of 105 Caucasian women in their reproductive age. We investigated the main factors (exposures) potentially linked with the presence of these resistance genes and with an overall composite ARGs score, with a particular focus on individual knowledge and attitudes toward antibiotics, antibiotic consumption and practices, and the characteristics of the vaginal microbiome, particularly its bacterial composition.
In particular, we assessed the presence of 14 ARGs and two related transposons: erm(A), erm(B), erm(F), tet(M), tet(M)-Tn916, tet(O), tet(W), tet(Q), tet(Q)-rteA, blaOXA-2, blaTEM, blaZ, blaSHV, blaCTX-M, qnrA and qepA. The exposures considered included the following demographic, and health-related behavioral data of the participants: age, educational, marital status, BMI, contraceptive use, smoking habits, presence/absence of Candida spp., and adherence to a healthy diet such as the Mediterranean diet, individual practice and awareness about antibiotics and antibiotic resistance.
Although a broader panel of ARGs has previously been investigated by other authors using next-generation sequencing (Bostwick et al., 2016), this is the first study to associate the presence of selected ARGs with individual lifestyle variables and antibiotic-related behaviors, in order to highlight their influence on ARG dissemination in the vaginal environment.
Our 16S rRNA sequencing analysis revealed that CSTs I and III were the most prevalent, followed by CST IV, with CSTs II and V being less common. These findings align with previous studies examining the vaginal microbiota in Caucasian reproductive-aged women (Ravel et al., 2011; Tamarelle et al., 2024).
Among the ARGs researched, the most prevalent ones were erm(F), tet(M), erm(B), erm(A) and tet(W), while blaZ, tet(Q) and blaTEM were less frequent. Notably, nearly all tet(M) and tet(Q)-positive samples also harbored their respective mobile elements Tn916 and rteA. These results are in line with previous studies on the gut and cervicovaginal microbiome, which revealed that the most abundant resistance determinants are related to tetracyclines and macrolides (Milanović et al., 2017; Roachford et al., 2021; Fadeyi et al., 2023). Similarly, the co-presence of tet(M) and Tn916 has been previously observed in VMB of pregnant women (Severgnini et al., 2021), suggesting the potential for horizontal gene transfer of tetracycline resistance genes between commensal and pathogenic vaginal bacteria (Chopra and Roberts, 2001).
When exploring the potential associations between the presences of ARGs, both as single gene and as composite ARGs score, and the considered exposures, interesting data emerged. For instance, we found that the use of contraceptives was negatively associated with the presence of tet(M). Research on the impact of oral contraceptives pills (OCP) has predominantly shown an increase in the relative abundance of beneficial Lactobacillus species and a reduction in bacterial vaginosis (BV)-associated taxa, thus resulting in a protective effect (Vodstrcil et al., 2013; Brooks et al., 2017). Conversely, copper intrauterine devices (Cu-IUDs) were associated with an increased BV risk (Brown et al., 2023), whereas depot medroxyprogesterone acetate (DMPA) and vaginal rings (VR) showed mixed results (Molatlhegi et al., 2020; Crucitti et al., 2018). All contraceptive users in our cohort reported using OCPs, supporting previous findings of their beneficial effect on vaginal health.
Smoking was positively associated with the presence of erm(B) and total ARG score. Bradshaw and colleagues showed that there was a direct and positive correlation between smoking habits and BV, whereby the greatest risk of BV was present in smokers of ≥30 cigarettes per week (Bradshaw et al., 2014). Smoking has been proven to have a negative impact on vaginal microbiota (VMB) associated with an increased risk for multiple and varied bacterial infections (Bagaitkar et al., 2008; Brotman et al., 2014). To our best knowledge, our work is the first to report an independent association between smoking habit and antibiotic resistance in the vaginal environment, although an increased AR for smokers has been previously reported for urinary tract infections treatments (Lorenzo-Gómez et al., 2020).
Similar data were observed with BMI. A higher BMI was positively associated with erm(A), tet(O) and total ARG score. Although the cohort had similar BMI values in the normal range, this association reflects previous findings linking higher BMI with increased ARG prevalence in the first trimester of pregnancy (Severgnini et al., 2021). Additionally, a higher MEDILITE score was negatively associated with tet(O). Our findings suggest that maintaining an appropriate body weight and adhering to a balanced diet are beneficial factors in preserving a healthy vaginal microbiome. Thoma and colleagues’ study on a cohort of 1735 nonpregnant women, found that those who consumed healthier diets were less likely to develop BV (Thoma et al., 2011). Moreover, a recent study that examined the relationship between dietary intake and VMB depicted how higher consumption of low-fat dairy, fruit, vitamin D, yogurt, and fiber were associated with a L. crispatus-dominated CST (Rosen et al., 2022).
Antibiotic-related behaviors emerged as influential factors: indeed, inappropriate use and low compliance were associated with increased ARG prevalence, whereas past compliance or higher awareness of antibiotic resistance were negatively associated with selected ARGs. A direct correlation between an increase in antibiotic use, mainly driven by inappropriate use and an increase in ARGs has been already reported (Cantón et al., 2013; Kunhikannan et al., 2021). These findings reinforce the importance of appropriate antibiotic use in limiting ARG dissemination. Besides, the presence of Candida spp. was significantly positively associated with the total ARG score, supporting existing evidence that antibiotics can promote Candida spp. overgrowth by disrupting the bacterial microbiota (Spinillo et al., 1999; Armstrong et al., 2024).
Alpha-diversity increased in tandem with the number of ARGs per sample, and higher ARG burden was linked to the enrichment of BV-associated bacterial taxa. In particular, Anaerococcus, Prevotella and Dialister were positively associated with resistance genes against macrolides and tetracyclines. It is worth mentioning that bacteria from different sites appear to be exchanging genetic material (Salyers et al., 2004). Also, some of the bacteria associated with the presence of ARGs in this study are commensal bacteria present in the gut microbiota and it is well-established that the human gut microbiota is an important reservoir for AR (Francino, 2016; Salyers et al., 2004). Thus, we can speculate that a bacterial translocation from the gastrointestinal-tract to the vagina through the gut-vagina axis has occurred (Amabebe and Anumba, 2020; Takada et al., 2023).
On the other hand, the genus Lactobacillus was negatively associated with all ARGs, with a few exceptions: positive association between L. gasseri and erm(A), L. iners and tet(Q) and, between unclassified members of the Lactobacillus genus and qnrA. This suggests that even commensal and health-promoting bacteria may acquire and harbor ARGs (Fadeyi et al., 2023).
This aspect has been highlighted by several studies suggesting that lactobacilli adapt to their environment by acquiring resistance genes from other bacteria through horizontal gene transfer. This phenomenon is exacerbated by several factors, such as the increasingly selective pressure exerted by antimicrobial treatments on lactobacilli colonizing the human gastrointestinal and vaginal tracts, as well as the widespread use of antibiotics in the food chain, where lactobacilli are often intentionally added as starter cultures (Campedelli et al., 2018; Colautti et al., 2022; Huang et al., 2024).
In this context, numerous studies have investigated how lactobacilli may act as reservoirs of antibiotic resistance genes that could potentially be transferred to pathogenic species (Campedelli et al., 2018; Colautti et al., 2022; Fadeyi et al., 2023). The most commonly studied resistances are those to tetracyclines and macrolides, associated with tet and erm genes, respectively (Colautti et al., 2022). Further studies are needed to better understand the role of Lactobacillus species as ARG reservoirs in relation to different CSTs (i.e., L. crispatus, L. gasseri, L. iners, and L. jensenii).
Other interesting results emerged when correlating together the relative abundances of microbial CAGs and the presence of ARGs. In particular, the L. crispatus/jenesenii/gasseri CAG was negatively correlated to the total ARG score, thus representing a potential ‘protective’ vaginal Lactobacillus community (Djusse et al., 2025). On the other hand, L. iners CAG exhibited neither positive nor negative correlation, strengthening the idea that the role of this Lactobacillus species in vaginal health remains unclear (Petrova et al., 2017). Indeed, L. iners has been considered a ‘transitional’ poorly protective species, typically associated with dysbiotic conditions. At the same time, it can be detected in normal conditions in a large subset of women, being its presence associated with young age and unprotected sexual practices (Novak et al., 2022).
Overall, our research delves into the complexity of the dynamics between the distribution of a panel of ARGs in the VMB of women in their reproductive age, with a particular emphasis on resistance to tetracyclines, macrolides, beta-lactams and quinolones. The high prevalence of ARGs such as erm(F), tet(M), erm(B) along with mobile genetic elements (transposons), underscores the potential for horizontal gene transfer within the vaginal ecosystem. Importantly, we identified several health-related behavioral, clinical, and demographic factors like contraceptive use, smoking, BMI, diet, presence of Candida spp. and antibiotic behavior that were significantly associated with the presence or absence of specific ARGs and to the total number of ARGs. These correlations reinforce the impact of lifestyle, health and antibiotic consumption on dynamics of antibiotic resistance. Notably, our findings suggest a protective role for oral contraceptives and healthy dietary patterns, while highlighting smoking and inappropriate antibiotic use as risk factors.
Taxonomic analysis further revealed a clear relationship between the presence of ARGs and an increase in bacterial diversity, particularly with the enrichment of BV-associated bacteria, while Lactobacillus spp. was negatively associated with ARG abundance. Interestingly, a few ARGs were also associated with commensal bacteria, suggesting that ARG dissemination may occur even in the absence of obvious risk factors.
Together, these findings emphasize the need for integrated public health strategies that combine antibiotic stewardship with broader lifestyle and behavioral interventions. Individualized therapeutic approaches, based on an individual’s resistome, may also represent a promising avenue for future treatment strategies.
We are fully aware of some limitations of this study. At first, even though we accounted for key factors potentially related to the distribution of ARGs in the VMB such as antibiotic behavior, contraceptive methods, smoking habits and diet, residual confounding cannot be excluded, especially regarding variables such as sexual behavior and microbiota exchanges with partners. Also, a shotgun metagenomic sequencing approach could provide a better comprehension of the resistance profiles actually present in some specific bacterial species. Additionally, the cross-sectional design limits causal inference, as both exposures and outcomes are measured simultaneously, thus it is difficult to establish a clear temporal sequence. Lastly, there is a sampling limitation, as the sample is not representative of the general Italian population. It consists of young women of reproductive age who voluntarily participated in the study and have a high level of education. This intrinsic selection bias substantially limits the generalizability of the findings, which should, therefore, be interpreted with caution. In particular, the results cannot be extrapolated to women of different age groups or educational backgrounds.
These findings lay the groundwork for longitudinal and functional studies to: (i) clarify causal relationships, (ii) explore ARG transmission via the gut-vagina axis, (iii) enable strain-level characterization of ARG-harboring bacteria through targeted bacterial isolation, (iv) analyze the functional activity/expression of detected ARGs, (v) assess a wider panel of genes, including those conferring resistance to 5-nitroimidazole, an antimicrobial commonly used in the vaginal setting (Williams et al., 2025).
Data availability statement
Raw sequencing data for this project are available in NCBI Short-Read Archive (SRA) under accession number PRJNA1188525 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1188525).
Ethics statement
The studies involving humans were approved by Institutional Review Board (IRB) of the University of Bologna, protocol number 0122421. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
PC: Investigation, Writing – original draft, Data curation. CCe: Writing – original draft, Investigation, Data curation. MED: Investigation, Writing – original draft, Data curation. MM: Writing – original draft, Conceptualization. SM: Data curation, Investigation, Writing – original draft. TC: Data curation, Writing – original draft, Investigation. SC: Investigation, Data curation, Writing – original draft. FP: Methodology, Data curation, Investigation, Writing – original draft. MS: Investigation, Data curation, Writing – review & editing, Methodology. CF: Writing – original draft, Conceptualization. MD: Writing – original draft, Investigation. CC: Data curation, Supervision, Methodology, Investigation, Writing – original draft. AM: Conceptualization, Writing – original draft.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This research received no specific grant from any funding agency in the public, commercial, or not for-profit sectors. Paola Castellano was supported by a doctoral fellowship of the National Operation Program (PON) Research and Innovation 2014-2020 (CCI. 2014IT16M2OP005), ESF REACT-EU resources, Action IV.4 “Doctorates and research contracts on themes of innovation” and Action IV.5 “Doctorates on Green themes.” CUP code J35F21003350006. Marielle Ezekielle Djusse was supported by the European Union-funded grant -NextGenerationEU from the National Recovery and Resilience Plan (NRRP) Mission 4, Component 1, Investment 3.4 (DM 118/2023) -Digital and Environmental Transitions -CUP: J33C23002390002. 14 The research group acknowledges the National Recovery and Resilience Plan (NRRP), Mission 4 Component 2 Investment 1.4 -Call for tender No. 3138 of 16 December 2021, rectified by Decree n.3175 of 18 December 2021 of Italian Ministry of University and Research funded by the European Union -NextGenerationEU; Award Number: Project code CN_00000033, Concession Decree No. 1034 of 17 June 2022 adopted by the Italian Ministry of University and Research, CUP: B83C22002930006, Project title “National Biodiversity Future Center -NBFC”, for supporting Camilla Ceccarani’s research grant. Silvia Conti was supported by the European Union-funded grant -NextGenerationEU -Mission 4, Component 2, Investiment 1.1, 2023-PRIN-VS_001 -PRIN 2022 -CUP F53D23007190006; and by the INVAT project, Fondo Ordinario dell’Ente 2022, National Research Council, CUP B53C22010140001.
Acknowledgments
We would like to thank the Clinical Laboratory Scientists of the Microbiology Unit of IRCCS Policlinico Sant’Orsola in Bologna for their assistance during the study.
Conflict of interest
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2025.1640992/full#supplementary-material
References
Abou Chacra, L., Fenollar, F., and Diop, K. (2022). Bacterial vaginosis: what do we currently know? Front. Cell Infect. Microbiol. 11. doi: 10.3389/fcimb.2021.672429
Ahrens, P., Andersen, L. O., Lilje, B., Johannesen, T. B., Dahl, E. G., Baig, S., et al. (2020). Changes in the vaginal microbiota following antibiotic treatment for Mycoplasma genitalium, Chlamydia trachomatis and bacterial vaginosis. PLoS One 15, e0236036. doi: 10.1371/journal.pone.0236036
Amabebe, E. and Anumba, D. O. C. (2020). Female gut and genital tract microbiota-induced crosstalk and differential effects of short-chain fatty acids on immune sequelae. Front. Immunol. 11. doi: 10.3389/fimmu.2020.02184
Armstrong, E., Hemmerling, A., Miller, S., Huibner, S., Kulikova, M., Liu, R., et al. (2024). Vaginal fungi are associated with treatment-induced shifts in the vaginal microbiota and with a distinct genital immune profile. Microbiol. Spectr. 12, e0350123. doi: 10.1128/spectrum.03501-23
Bagaitkar, J., Demuth, D. R., and Scott, D. A. (2008). Tobacco use increases susceptibility to bacterial infection. Tob Induc Dis. 4, 12. doi: 10.1186/1617-9625-4-12
Bagley, S. M., Schoenberger, S. F., dellaBitta, V., Lunze, K., Barron, K., Hadland, S. E., et al. (2022). An exploration of young adults with opioid use disorder and how their perceptions of family members’ Beliefs affects medication treatment. J. Addict. Med. 16, 689–694. doi: 10.1097/ADM.0000000000001001
Borgogna, J. C., Shardell, M. D., Santori, E. K., Nelson, T. M., Rath, J. M., Glover, E. D., et al. (2020). The vaginal metabolome and microbiota of cervical HPV-positive and HPV-negative women: a cross-sectional analysis. BJOG 127, 182–192. doi: 10.1111/1471-0528.15981
Bostwick, D. G., Woody, J., Hunt, C., and Budd, W. (2016). Antimicrobial resistance genes and modelling of treatment failure in bacterial vaginosis: clinical study of 289 symptomatic women. J. Med. Microbiol. 65, 377–386. doi: 10.1099/jmm.0.000236
Bradshaw, C. S., Walker, S. M., Vodstrcil, L. A., Bilardi, J. E., Law, M., Hocking, J. S., et al. (2014). The influence of behaviors and relationships on the vaginal microbiota of women and their female partners: the WOW Health Study. J. Infect. Dis. 209, 1562–1572. doi: 10.1093/infdis/jit664
Brooks, J. P., Edwards, D. J., Blithe, D. L., Fettweis, J. M., Serrano, M. G., Sheth, N. U., et al. (2017). Effects of combined oral contraceptives, depot medroxyprogesterone acetate and the levonorgestrel-releasing intrauterine system on the vaginal microbiome. Contraception 95, 405–413. doi: 10.1016/j.contraception.2016.11.006
Brotman, R. M., He, X., Gajer, P., Fadrosh, D., Sharma, E., Mongodin, E. F., et al. (2014). Association between cigarette smoking and the vaginal microbiota: a pilot study. BMC Infect. Dis. 14, 471. doi: 10.1186/1471-2334-14-471
Brown, B. P., Feng, C., Tanko, R. F., Jaumdally, S. Z., Bunjun, R., Dabee, S., et al. (2023). Copper intrauterine device increases vaginal concentrations of inflammatory anaerobes and depletes lactobacilli compared to hormonal options in a randomized trial. Nat. Commun. 14, 499. doi: 10.1038/s41467-023-36002-4
Campedelli, I., Mathur, H., Salvetti, E., Clarke, S., Rea, M. C., Torriani, S., et al. (2018). Genus-wide assessment of antibiotic resistance in Lactobacillus spp. Appl. Environ. Microbiol. 85, e01738–e01718. doi: 10.1128/AEM.01738-18
Cantón, R., Horcajada, J. P., Oliver, A., Garbajosa, P. R., and Vila, J. (2013). Inappropriate use of antibiotics in hospitals: the complex relationship between antibiotic use and antimicrobial resistance. Enferm Infecc Microbiol. Clin. 31 Suppl 4, 3–11. doi: 10.1016/S0213-005X(13)70126-5
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Castellano, P., Russo, P. M., and Mazzetti, M. (2024). Exploring the influence of family attitudes and individual psychological factors on antibiotic utilization: A pilot study. Healthcare (Basel). 12, 1213. doi: 10.3390/healthcare12121213
Ceccarani, C., Foschi, C., Parolin, C., D’Antuono, A., Gaspari, V., Consolandi, C., et al. (2019). Diversity of vaginal microbiome and metabolome during genital infections. Sci. Rep. 9, 14095. doi: 10.1038/s41598-019-50410-x
Chopra, I. and Roberts, M. (2001). Tetracycline antibiotics: mode of action, applications, molecular biology, and epidemiology of bacterial resistance. Microbiol. Mol. Biol. Rev. 65, 232–260. doi: 10.1128/MMBR.65.2.232-260.2001
Claesson, M. J., Jeffery, I. B., Conde, S., Power, S. E., O’Connor, E. M., Cusack, S., et al. (2012). Gut microbiota composition correlates with diet and health in the elderly. Nature 488, 178–184. doi: 10.1038/nature11319
Colautti, A., Arnoldi, M., Comi, G., and Iacumin, L. (2022). Antibiotic resistance and virulence factors in lactobacilli: something to carefully consider. Food Microbiol. 103, 103934. doi: 10.1016/j.fm.2021.103934
Crucitti, T., Hardy, L., van de Wijgert, J., Agaba, S., Buyze, J., Kestelyn, E., et al. (2018). Contraceptive rings promote vaginal lactobacilli in a high bacterial vaginosis prevalence population: A randomised, open-label longitudinal study in Rwandan women. PLoS One 13, e0201003. doi: 10.1371/journal.pone.0201003
Das, S., Bhattacharjee, M. J., Mukherjee, A. K., and Khan, M. R. (2023). Recent advances in understanding of multifaceted changes in the vaginal microenvironment: implications in vaginal health and therapeutics. Crit. Rev. Microbiol. 49, 256–282. doi: 10.1080/1040841X.2022.2049696
Djusse, M. E., Prinelli, F., Camboni, T., Ceccarani, C., Consolandi, C., Conti, S., et al. (2025). Dietary habits and vaginal environment: can a beneficial impact be expected? Front. Cell Infect. Microbiol. 15. doi: 10.3389/fcimb.2025.1582283
Fadeyi, T. E., Oyedemi, O. T., Awe, O. O., and Ayeni, F. (2023). Antibiotic use in infants within the first year of life is associated with the appearance of antibiotic-resistant genes in their feces. PeerJ 11, e15015. doi: 10.7717/peerj.15015
Flórez, A. B., Campedelli, I., Delgado, S., Alegría, Á., Salvetti, E., Felis, G. E., et al. (2016). Antibiotic susceptibility profiles of dairy leuconostoc, analysis of the genetic basis of atypical resistances and transfer of genes in vitro and in a food matrix. PloS One 11, e0145203. doi: 10.1371/journal.pone.0145203
France, M., Alizadeh, M., Brown, S., Ma, B., and Ravel, J. (2022). Towards a deeper understanding of the vaginal microbiota. Nat. Microbiol. 7, 367–378. doi: 10.1038/s41564-022-01083-2
France, M. T., Ma, B., Gajer, P., Brown, S., Humphrys, M. S., Holm, J. B., et al. (2020). VALENCIA: a nearest centroid classification method for vaginal microbial communities based on composition. Microbiome 8, 166. doi: 10.1186/s40168-020-00934-6
Francino, M. P. (2016). Antibiotics and the human gut microbiome: dysbioses and accumulation of resistances. Front. Microbiol. 6. doi: 10.3389/fmicb.2015.01543
Georgakopoulou, V. E., Lempesis, I. G., Tarantinos, K., Sklapani, P., Trakas, N., and Spandidos, D. A. (2024). Atypical pneumonia (Review). Exp. Ther. Med. 28, 424. doi: 10.3892/etm.2024.12713
Greenbaum, S., Greenbaum, G., Moran-Gilad, J., and Weintraub, A. Y. (2019). Ecological dynamics of the vaginal microbiome in relation to health and disease. Am. J. Obstet Gynecol. 220, 324–335. doi: 10.1016/j.ajog.2018.11.1089
Huang, L., Guo, R., Li, S., Wu, X., Zhang, Y., Guo, S., et al. (2024). A multi-kingdom collection of 33,804 reference genomes for the human vaginal microbiome. Nat. Microbiol. 9, 2185–2200. doi: 10.1038/s41564-024-01751-5
Jeters, R. T., Rivera, A. J., Boucek, L. M., Stumpf, R. M., Leigh, S. R., and Salyers, A. A. (2009). Antibiotic resistance genes in the vaginal microbiota of primates not normally exposed to antibiotics. Microb. Drug Resist. 15, 309–315. doi: 10.1089/mdr.2009.0052
Jones, D. L. (2017). Fathom Toolbox for MATLAB: software for multivariate ecological and oceanographic data analysis (St. Petersburg, FL, USA: College of Marine Science, University of South Florida). Available online at: https://www.usf.edu/marine-science/research/matlab-resources/fathom-toolbox-for-matlab.aspx.
Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., et al. (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41, e1. doi: 10.1093/nar/gks808
Kunhikannan, S., Thomas, C. J., Franks, A. E., Mahadevaiah, S., Kumar, S., and Petrovski, S. (2021). Environmental hotspots for antibiotic resistance genes. Microbiologyopen 10, e1197. doi: 10.1002/mbo3.1197
Lorenzo-Gómez, M. F., Santos-Antunes, M. T., Nieto-Huertos, A., Lorenzo-Gómez, A., Marquez-Sanchez, M. T., Flores-Fraile, M. C., et al. (2020). The influence of smoking on bacterial resistance after vaccine or antibiotic prophylaxis against recurrent urinary tract infections. Actas Urol Esp (Engl Ed). 44, 497–504. doi: 10.1016/j.acuro.2020.04.002
Lozupone, C., Lladser, M. E., Knights, D., Stombaugh, J., and Knight, R. (2011). UniFrac: an effective distance metric for microbial community comparison. ISME J. 5, 169–172. doi: 10.1038/ismej.2010.133
Melkumyan, A. R., Priputnevich, T. V., Ankirskaya, A. S., Murav’eva, V. V., and Lubasovskaya, L. A. (2015). Effects of antibiotic treatment on the lactobacillus composition of vaginal microbiota. Bull. Exp. Biol. Med. 158, 766–768. doi: 10.1007/s10517-015-2857-1
Milanović, V., Osimani, A., Aquilanti, L., Tavoletti, S., Garofalo, C., Polverigiani, S., et al. (2017). Occurrence of antibiotic resistance genes in the fecal DNA of healthy omnivores, ovo-lacto vegetarians and vegans. Mol. Nutr. Food Res. 61. doi: 10.1002/mnfr.201601098
Molatlhegi, R. P., Liebenberg, L. J., Leslie, A., Noel-Romas, L., Mabhula, A., Mchunu, N., et al. (2020). Plasma concentration of injectable contraceptive correlates with reduced cervicovaginal growth factor expression in South African women. Mucosal Immunol. 13, 449–459. doi: 10.1038/s41385-019-0249-y
Morsli, M., Gimenez., E., Magnan, C., Salipante, F., Huberlant, S., Letouzey, V., et al. (2024). The association between lifestyle factors and the composition of the vaginal microbiota: a review. Eur. J. Clin. Microbiol. Infect. Dis. 43, 1869–1881. doi: 10.1007/s10096-024-04915-7
Nortey, R. A., Kretchy, I. A., Koduah, A., and Buabeng, K. O. (2023). Biopsychosocial analysis of antibiotic use for the prevention or management of COVID-19 infections: A scoping review. Res. Soc. Adm Pharm. 19, 573–581. doi: 10.1016/j.sapharm.2022.11.011
Novak, J., Ravel, J., Ma, B., Ferreira, C. S. T., Tristão, A. D. R., Silva, M. G., et al. (2022). Characteristics associated with Lactobacillus iners-dominated vaginal microbiota. Sex Transm Infect. 98, 353–359. doi: 10.1136/sextrans-2020-054824
Petrova, M. I., Reid, G., Vaneechoutte, M., and Lebeer, S. (2017). Lactobacillus iners: friend or foe? Trends Microbiol. 25, 182–191. doi: 10.1016/j.tim.2016.11.007
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Ravel, J., Gajer, P., Abdo, Z., Schneider, G. M., Koenig, S. S., McCulle, S. L., et al. (2011). Vaginal microbiome of reproductive-age women. Proc. Natl. Acad. Sci. U.S.A. 108 Suppl 1, 4680–4687. doi: 10.1073/pnas.1002611107
Roachford, O. S. E., Alleyne, A. T., Kuelbs, C., Torralba, M. G., and Nelson, K. E. (2021). The cervicovaginal microbiome and its resistome in a random selection of Afro-Caribbean women. Hum. Microbiome J. 20, 100079. doi: 10.1016/j.humic.2021.100079
Rosen, E. M., Martin, C. L., Siega-Riz, A. M., Dole, N., Basta, P. V., Serrano, M., et al. (2022). Is prenatal diet associated with the composition of the vaginal microbiome? Paediatr. Perinat Epidemiol. 36, 243–253. doi: 10.1111/ppe.12830
Salyers, A. A., Gupta, A., and Wang, Y. (2004). Human intestinal bacteria as reservoirs for antibiotic resistance genes. Trends Microbiol. 12, 412–416. doi: 10.1016/j.tim.2004.07.004
Severgnini, M., Camboni, T., Ceccarani, C., Morselli, S., Cantiani, A., Zagonari, S., et al. (2021). Distribution of ermB, ermF, tet(W), and tet(M) Resistance Genes in the Vaginal Ecosystem of Women during Pregnancy and Puerperium. Pathogens 10, 1546. doi: 10.3390/pathogens10121546
Severgnini, M., Morselli, S., Camboni, T., Ceccarani, C., Laghi, L., Zagonari, S., et al. (2022). A deep look at the vaginal environment during pregnancy and puerperium. Front. Cell Infect. Microbiol. 12. doi: 10.3389/fcimb.2022.838405
Sirichoat, A., Flórez, A. B., Vázquez, L., Buppasiri, P., Panya, M., Lulitanond, V., et al. (2020). Antibiotic susceptibility profiles of lactic acid bacteria from the human vagina and genetic basis of acquired resistances. Int. J. Mol. Sci. 21, 2594. doi: 10.3390/ijms21072594
Sofi, F., Dinu, M., Pagliai, G., Marcucci, R., and Casini, A. (2017). Validation of a literature-based adherence score to Mediterranean diet: the MEDI-LITE score. Int. J. Food Sci. Nutr. 68, 757–762. doi: 10.1080/09637486.2017.1287884
Spinillo, A., Capuzzo, E., Acciano, S., De Santolo, A., and Zara, F. (1999). Effect of antibiotic use on the prevalence of symptomatic vulvovaginal candidiasis. Am. J. Obstet Gynecol. 180, 14–17. doi: 10.1016/s0002-9378(99)70141-9
Takada, K., Melnikov, V. G., Kobayashi, R., Komine-Aizawa, S., Tsuji, N. M., and Hayakawa, S. (2023). Female reproductive tract-organ axes. Front. Immunol. 14. doi: 10.3389/fimmu.2023.1110001
Tamarelle, J., Thiébaut, A. C. M., de Barbeyrac, B., Bébéar, C., Bourret, A., Fauconnier, A., et al. (2024). Vaginal microbiota stability over 18 months in young student women in France. Eur. J. Clin. Microbiol. Infect. Dis. 43, 2277–2292. doi: 10.1007/s10096-024-04943-3
Thoma, M. E., Klebanoff, M. A., Rovner, A. J., Nansel, T. R., Neggers, Y., Andrews, W. W., et al. (2011). Bacterial vaginosis is associated with variation in dietary indices. J. Nutr. 141, 1698–1704. doi: 10.3945/jn.111.140541
Vodstrcil, L. A., Hocking, J. S., Law, M., Walker, S., Tabrizi, S. N., Fairley, C. K., et al. (2013). Hormonal contraception is associated with a reduced risk of bacterial vaginosis: a systematic review and meta-analysis. PLoS One 8, e73055. doi: 10.1371/journal.pone.0073055
Wang, N., Chen, L., Yi, K., Zhang, B., Li, C., and Zhou, X. (2024). The effects of microbiota on reproductive health: A review. Crit. Rev. Food Sci. Nutr. 64, 1486–1507. doi: 10.1080/10408398.2022.2117784
Wang, Z., Li, Y., Wang, J., and Li, S. (2024). Tetracycline antibiotics in agricultural soil: Dissipation kinetics, transformation pathways, and structure-related toxicity. Sci. Total Environ. 949, 175126. doi: 10.1016/j.scitotenv.2024.175126
Williams, A., Ravel, J., Armstrong, E., Huibner, S., Rutt, L., Kaul, R., et al. (2025). Temporal dynamics of the vaginal microbiome and host immune markers before, during, and after metronidazole treatment for bacterial vaginosis. mSystems 10, e0038025. doi: 10.1128/msystems.00380-25
Keywords: vaginal microbiota, antimicrobial resistance, ARG, resistome, women’s health
Citation: Castellano P, Ceccarani C, Djusse ME, Mazzetti M, Morselli S, Camboni T, Conti S, Prinelli F, Severgnini M, Foschi C, Dall’Asta M, Consolandi C and Marangoni A (2025) Linking antibiotic resistance genes in the vaginal microbiota to health-related behaviors and antibiotic awareness in reproductive-age women: a cross-sectional study. Front. Cell. Infect. Microbiol. 15:1640992. doi: 10.3389/fcimb.2025.1640992
Received: 04 June 2025; Accepted: 03 September 2025;
Published: 18 September 2025.
Edited by:
Rodolfo García-Contreras, National Autonomous University of Mexico, MexicoReviewed by:
Valentina Margarita, University of Sassari, ItalyDonghyeok Seol, Seoul National University Bundang Hospital, Republic of Korea
Copyright © 2025 Castellano, Ceccarani, Djusse, Mazzetti, Morselli, Camboni, Conti, Prinelli, Severgnini, Foschi, Dall’Asta, Consolandi and Marangoni. 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.
*Correspondence: Claudio Foschi, Y2xhdWRpby5mb3NjaGkyQHVuaWJvLml0
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share last authorship