ABSTRACT
Bacterial vaginosis (BV) is a common vaginal disorder characterized by the decrease of lactobacilli and overgrowth of Gardnerella vaginalis and resident anaerobic vaginal bacteria. In the present work, the effects of rifaximin vaginal tablets on vaginal microbiota and metabolome of women affected by BV were investigated by combining quantitative PCR and a metabolomic approach based on 1H nuclear magnetic resonance. To highlight the general trends of the bacterial communities and metabolomic profiles in response to the antibiotic/placebo therapy, a multivariate statistical strategy was set up based on the trajectories traced by vaginal samples in a principal component analysis space. Our data demonstrated the efficacy of rifaximin in restoring a health-like condition in terms of both bacterial communities and metabolomic features. In particular, rifaximin treatment was significantly associated with an increase in the lactobacillus/BV-related bacteria ratio, as well as with an increase in lactic acid concentration and a decrease of a pool of metabolites typically produced by BV-related bacteria (acetic acid, succinate, short-chain fatty acids, and biogenic amines). Among the tested dosages of rifaximin (100 and 25 mg for 5 days and 100 mg for 2 days), 25 mg for 5 days was found to be the most effective.
INTRODUCTION
Bacterial vaginosis (BV) is a complex polymicrobial vaginal disorder associated with an increase in the taxonomic richness and diversity of the vaginal microbiota. BV is microbiologically characterized by replacement of the lactobacillus-predominant vaginal microbiota by potential pathogenic anaerobic bacteria (1, 2). The diagnosis of BV is based on Amsel's criteria (1) and Nugent score determinations (3). The current recommended treatment of BV includes metronidazole (oral or vaginal) or clindamycin (vaginal) (4); however, the short-term (30 days) cure rate is often poor, and recurrences are common (5, 6).
Rifaximin is a new candidate for the local treatment of BV thanks to its antibacterial activity, which covers Gardnerella vaginalis and other pathogens responsible for urogenital infections (7, 8). Rifaximin is a semisynthetic rifamycin derivative characterized by low systemic absorption and good antibacterial activity. In common with the structural analogue rifampin and other members of the rifamycin class, it acts on the β subunit of the bacterial RNA polymerase to inhibit RNA synthesis (9, 10). Recently, the efficacy of rifaximin vaginal tablets in the treatment of BV was demonstrated by evaluating the rate of clinical remission (11) and the restoration of normal vaginal communities (12) and proteome profiles (13).
BV, as well as antibiotics used for the treatment of this disturbance, cause perturbations of the vaginal ecosystem, which are reflected in the overall profile of the metabolites produced by the host-bacterium metaorganism. The study of the comprehensive modifications occurring in the metabolic profiles of living systems actually represents the main field of metabolomics. Metabolomics is, by definition, the study of the effects of an external perturbation on a biological system by looking for changes in its metabolome, that is, the complete set of small-molecule metabolites within it (14). Some of these molecules represent nodes of a high number of metabolic cycles, so that little variations in the balance between the microorganisms of an environment can lead to amplified and more easily studied metabolome modifications (15, 16).
High-resolution 1H nuclear magnetic resonance (HR-NMR) spectroscopy, when applied to metabolomic studies, has provided significant information on a wide range of pathologies, like cancer (17), meningitis (18), and coronary heart disease (19), as well as a variety of gastrointestinal diseases (20). This is because an NMR profile contains qualitative and quantitative information on the hundreds of different small molecules present in a sample at >1 μM concentration (21). Moreover, NMR-based metabolomics makes no assumption on the identity of the metabolites that are relevant for the selected pathology, because information on the significant metabolic pattern alterations is directly obtained through statistical analysis of the NMR profiles (22, 23).
The present work describes a combined approach based on quantitative PCR (qPCR) of bacterial 16S rRNA genes and metabolomics to characterize the vaginal environment of women affected by BV as well as to search for modifications induced by the administration of different dosages of rifaximin vaginal tablets. To the best of our knowledge, this is the first work describing the vaginal environment through a metabolomic approach based on NMR, as well as the first describing, by means of this technique, the evolution of a vaginal syndrome following the administration of a drug.
MATERIALS AND METHODS
Study design and sample preparation.Vaginal samples from 92 premenopausal women affected by BV were analyzed in the present study. Demographic characteristics and sexual habits of these patients are reported in Table 1. The patients had previously been enrolled in a multicenter, double-blind, randomized, placebo-controlled study (EudraCT; 2009-011826-32) performed to compare the efficacy of rifaximin vaginal tablets to that of the placebo for the treatment of BV (11, 12). Diagnosis of BV was made using both Amsel's criteria and Nugent score at the screening visit (V1). Patients with a Nugent score of >3 and positive for ≥3 of Amsel's criteria were considered positive for BV. At the randomization visit (V2), the patients were distributed into 4 treatment groups: group A received a 100-mg rifaximin vaginal tablet once daily for 5 days (n = 23), group B received a 25-mg rifaximin vaginal tablet once daily for 5 days (n = 23), group C received a 100-mg rifaximin vaginal tablet once daily for the first 2 days and a placebo vaginal tablet for the remaining 3 days (n = 23), and group D received a placebo vaginal tablet once daily for 5 days (n = 23). The medication was dispensed by the investigator in accordance with a computer-generated random allocation sequence (SAS release 9.2; SAS Institute, Cary, NC, USA) prepared by a statistician who was not blind to the treatment regimen. Seven days after the end of the therapy, the first follow-up visit (V3) was performed. Patients showing remission according to Amsel's criteria (≤2) and Gram stain Nugent score (≤3) attended the second follow-up visit (V4) 28 days after the end of the treatment. Participating women were required to abstain from intercourse during the 5-day treatment period and for 3 days before the follow-up visits, to avoid the use of intravaginal products, including douches, sprays, tampons, spermicides, gels, foams, and diaphragms, and to adopt an adequate contraceptive method during the study. Women were excluded from the study if they were pregnant or breast-feeding, anticipated menses at screening or follow-up visits, or had received systemic or vaginal antimicrobial therapy in the 2 weeks before the study.
Demographic characteristics and sexual habits of the patients analyzed
A control group consisted of 10 healthy subjects who had no signs of vaginal tract infection and had never had BV. These patients were enrolled in the Heilig Hart Hospital of Tienen (Belgium) and signed informed consent in accordance with the approval from the local ethics committee.
Rifaximin and placebo vaginal tablets were purchased from Alfa Wassermann S.p.A., Bologna, Italy.
Standardized vaginal rinsings with 2 ml of saline were collected at V1, V3, and V4 by flushing and reaspirating the fluid through a 22-gauge needle in the left, central, and right upper vaginal vaults (24) and stored at −80°C until use.
One ml of each vaginal rinsing was centrifuged at 9,500 × g for 15 min to separate the pellet, which was processed for bacterial DNA isolation, from the supernatant used for metabolomic analysis.
DNA isolation was carried out as previously described (25). The DNA amount was quantified using a NanoDrop ND-1000 (NanoDrop Technologies, Wilmington, DE).
qPCR.Quantitative PCR (qPCR) was performed on DNA samples extracted from vaginal fluids using a LightCycler instrument (Roche, Mannheim, Germany) and SYBR green I as the reporter fluorophore. Genus- or species-specific primer sets targeted to the 16S rRNA gene or 16S-23S rRNA spacer region were used to amplify bacteria belonging to Lactobacillus (Bact-0011f/Lab-0677r [26]), Gardnerella vaginalis (F-GV1/R-GV3 [27]), Atopobium (c-Atopo-f/c-Atopo-r [28]), Prevotella (g-Prevo-f/g-Prevo-r [29]), Veillonella (VeilloF/VeilloR [30]), Mycoplasma hominis (MycF/MycR [31]), and Mobiluncus (Mob-s/Mob-as [32]). The reaction mixture contained each primer at 0.5 μM, 4 μl of LightCycler-FastStart DNA master SYBR green I (Roche), and 2 μl of template (final volume of 20 μl). The thermal cycling conditions were optimized as an initial denaturation step at 95°C for 10 min followed by 30 (Lactobacillus, Atopobium, G. vaginalis, Veillonella, and Mobiluncus), 35 (Prevotella), or 40 (M. hominis) cycles of denaturation at 95°C for 15 s; primer annealing at 63°C (Lactobacillus and Mobiluncus), 62°C (Veillonella and M. hominis), or 60°C (Atopobium, Prevotella, and G. vaginalis) for 20 s; extension at 72°C for 45 s (Lactobacillus, Atopobium, Prevotella, G. vaginalis, and Veillonella) or 30 s (M. hominis and Mobiluncus); and a fluorescence acquisition step at 82°C (M. hominis), 85°C (Lactobacillus, Atopobium, G. vaginalis, and Veillonella), 87°C (Prevotella), or 88°C (Mobiluncus) for 5 s. Amplifications were performed in triplicate for each primer set. Data were expressed as ng of DNA of the targeted genus or species per μg of total DNA extracted from the vaginal sample.
1H NMR measurements.One ml of thawed sample was centrifuged at 14,000 × g for 5 min and then added to 160 μl of a D2O solution of 6.25 mM 3-(trimethylsilyl)-propionic-2,2,3,3–d4 acid sodium salt (TSP). The pH was adjusted to 7.00 with the addition of 0.5 M solution of HCl or NaOH. 1H-NMR spectra were recorded at 298 K with an Avance spectrometer (Bruker, Milan, Italy) operating at a frequency of 600.13 MHz and equipped with an autosampler with 60 holders. The HOD residual signal was suppressed by applying the first increment of the nuclear Overhauser effect spectroscopy (NOESY) pulse sequence and a spoil gradient (33, 34). This was done by employing the NOESYGPPR1D sequence, part of the standard pulse sequence library. Each spectrum was acquired using 32,000 data points over a 7,211.54 Hz spectral width and adding 256 transients. A recycle delay of 5 s and a 90° pulse of 11.4 μs were set up. Acquisition time (2.27 s) and recycle delay were adjusted to be 5 times longer than the longitudinal relaxation time of the protons under investigation, which was considered to be not longer than 1.4 s.
The peaks were assigned by comparing their chemical shift and multiplicity to data in the literature (35) and the Amix software data bank (ver. 3.9.7; Bruker, Italy).
Data analysis.Statistical analyses were performed using the R computational language (36). Each NMR spectrum was processed by means of scripts developed in-house as follows. The spectrum baseline was adjusted by employing the peak identification algorithm described by Coombes et al. (37) and implemented by Liland et al. (38) under the name “baseline.peakDetection”. Spectra were normalized by means of the probabilistic quotient normalization method (39). Similarities among the analyzed samples, as well as trends in their microbial composition or metabolome profile, were investigated by means of principal component analysis (PCA). For this purpose, the data were means centered and subjected to singular value decomposition (40). PCA calculates linear combinations of the calculated variables, points of NMR spectra, or bacterial concentrations to build new axes, called principal components (PC), oriented along the maximum variance directions. The results of PC analysis were presented in the text as biplots that are superimpositions of scoreplots and loadingplots. The scoreplot is the representation of the samples in the PC space, while the loadingplot is a graphical representation of the variables which contribute most to the spreading of the samples. The scoreplot is one of the best options to highlight the trends characterizing the samples, while the loadingplot highlights the molecules determining the trends so identified. Differences in the amounts of bacterial genera and species, in metabolite concentrations, and in PCA data were analyzed using Wilcoxon's signed-rank test. Unless differently specified, a P value below 0.01 was considered significant in all statistical tests.
RESULTS
Clinical results.Remission (R) or nonremission (N) was evaluated at the first (V3) and second (V4) follow-up visits by the combined assessment of Amsel's criteria (clinical cure) and the Gram stain Nugent scoring system (Fig. 1). The highest remission rate at V3 was found for group B (12/23), which also showed the best remission maintenance at V4 (7/12). Notably, the lowest rate of remission at V3 was found for the placebo group (4/19), which recorded a good maintenance of spontaneous remission at V4 (2/4). Similarly, in the work of Donders et al. (11), treatment group B showed the best therapeutic cure and maintenance versus placebo. Demographic parameters and sexual behavior (Table 1) did not influence clinical remission, as previously reported (11).
Remission (R) or nonremission (N) evaluated at the follow-up visits V3 and V4 for the treatment groups A (100 mg/day for 5 days), B (25 mg/day for 5 days), C (100 mg/day for 2 days), and D (placebo for 5 days) and for the totality of the women (TOT).
Short-term impact of rifaximin/placebo on vaginal microbiota.Lactobacillus, G. vaginalis, Atopobium, Prevotella, Veillonella, M. hominis, and Mobiluncus were quantified by qPCR (see Table S1 in the supplemental material) in order to investigate the effects of rifaximin on the principal bacterial genera and species which are unbalanced in the presence of BV (41–44). Table 2 lists the bacterial groups that significantly changed between V1 and V3 in women going into remission at V3. A significant decrease was observed for G. vaginalis, Atopobium, and Prevotella, confirming that these bacteria are specific markers for BV condition (42). To obtain an overview of the variations in bacterial populations associated with remission, a PCA was performed on the qPCR data related to BV-affected women at V1 and V3 and to healthy controls (Fig. 2A). The first two PCs accounted for 44.5% of the whole variance of the investigated samples. Each group of women was spread considerably on the built space, showing that a great source of variability was represented by the individual differences existing among the vaginal microbial communities of different subjects. The general impact exerted by rifaximin/placebo treatment on the vaginal microbiota composition was highlighted by plotting the medians from the groups of women. The medians from the samples collected at V1 and V3 from BV-affected women going into remission (V1→R-V3) were significantly different, showing that a part of the sample variability was due to the evolution of BV syndrome. For these women, remission caused an increase of the lactobacillus/BV-related bacteria ratio and, consequently, a trajectory approaching that of healthy controls in the scoreplot. In this respect, Atopobium, Prevotella, Veillonella, and M. hominis underwent the greatest changes, as revealed by the loadingplot superimposed on the scoreplot. Interestingly, the lactobacillus/BV-related bacteria ratio also significantly increased in the women where BV persisted at V3 (V1→N-V3). In detail, the distance covered by the medians of such samples was 74% of the one covered by women going into remission. Notably, on average, women going into remission showed at the baseline a more favorable microbiological environment than the others.
Concentration of bacterial groups which significantly varied (P < 0.05) between V1 and V3 in women going into remission at V3 (n = 29)
(A) Biplot of a PCA performed on the autoscaled qPCR data related to BV-affected women at visits V1 and V3 and healthy control women. The squares and circles represent BV-affected women at V1 and V3, respectively, with the open and filled symbols indicating women in remission (R) or not in remission (N), respectively, after rifaximin or placebo treatment. The gray triangles represent healthy women (H). Median values of the sample groups are indicated as V1, R-V3, N-V3, and H. (B) Biplot of a PCA performed on qPCR data related to BV-affected women at visit V1 and to healthy women. R and N indicate the median values of women grouped according to the response to rifaximin (RRifa [n = 25] and NRifa [n = 44], respectively) or placebo (RPbo [n = 4] and NPbo [n = 19], respectively) at V3, while H indicates the median values of healthy women (H; n = 10). Both trajectories (V1→R-V3 and V1→N-V3) indicate significant differences. Abbreviations: Lact, Lactobacillus; Gard, G. vaginalis; Atop, Atopobium; Prev, Prevotella; Veil, Veillonella; Mhom, M. hominis; Mobi, Mobiluncus; expl.var, explained variance.
Correlation between baseline microbiota and remission.Samples collected from BV-affected women at V1 were evaluated by PCA in a relationship between remission condition at V3 (N or R) and the type of treatment applied (rifaximin or placebo) and were compared to healthy controls (H) (Fig. 2B). Although no significant difference was found among the microbiota associated with the 5 analyzed groups, women treated with placebo showed an interesting trend, which suggested a cause/effect relationship between high lactobacillus/BV-related bacteria ratios at the baseline and the possibility of going into spontaneous remission. This trend also suggests that women who went into spontaneous remission did not have true cases of BV.
Long-term impact of placebo on vaginal microbiota.Figure 3 illustrates the trajectories covered by women who underwent the second follow-up visit, V4, in the space built by means of PCA on the qPCR data. In the case of placebo treatment (Fig. 3A), women who still had BV at V3 and women who were in remission at both V3 and V4 were characterized by the lowest and highest ratios between lactobacilli and BV-related bacteria registered at V1, respectively. This ratio was intermediate in women who went into remission at V3 but presenting relapse at V4. After the first visit, neither group of women showed any significant modification of the microbiota or correlations with BV evolution.
Biplot of a PCA calculated on the autoscaled qPCR data related to healthy women (triangle) and BV-affected women in remission at both V3 and V4 (circles), in remission at V3 but not at V4 (squares), and not in remission at V3 (no symbol). To limit the sample superimposition, data are represented separately for women subjected to placebo (A) and rifaximin (B) treatment. Median values of the samples collected from healthy women (H) and BV-affected women at V1, V3, and V4 are shown. The number of women (n) for each sample group is reported in Fig. 1. Continuous lines indicate significant differences between V1–V3 and V3–V4 trajectories, while dashed lines indicate nonsignificant differences. Full V4 symbols indicate significant differences between V1 and V4. expl.var, explained variance.
Long-term impact of rifaximin on vaginal microbiota.The trajectories covered by women treated with rifaximin (Fig. 3B) were clearly different from those of women treated with placebo. V3 and V4 samples of women in remission at both follow-up visits were localized in the same position, characterized by a favorable balance between lactobacilli and BV-related bacteria. This position was very close to that of healthy controls and significantly distant from V1. In contrast to those of women in remission at both follow-up visits, the trajectories of V4 samples of women who were not in remission at the second follow-up visit moved away significantly from the V3 trajectory in the direction of V1. For these women, V1 and V4 samples were characterized by a less favorable balance in the microbiota than V3 samples. In women who still had BV after rifaximin treatment, V3 samples were closer to V1 samples and were characterized by the highest concentrations of BV-related bacteria. Nevertheless, the V3–V1 distance was statistically significant even for this group of patients.
Impact of rifaximin on metabolome.An example of an HR-NMR spectrum, describing the metabolomic profile of vaginal fluid collected from a woman affected by BV, is shown in Fig. 4. The spectrum was dominated in intensity and variance by the peaks ascribable to lactate, whose area ranged between 7% and 72% of the entire spectrum. Molecules associated with the growth of BV-related bacteria could be identified, such as short-chain fatty acids (acetate, propionate, butyrate, isobutyrate, and isovaleric acid) (45, 46), organic acids (succinate, formiate, and fumarate) (47), and biogenic amines (putrescine, spermidine, and tyramine) (6, 48). Peaks from amino acids and molecules from purine catabolism (hypoxanthine and xanthine) were also identified.
1H-NMR spectrum of a vaginal fluid sample collected from a BV-affected woman at the baseline.
Table 3 lists the molecules whose concentrations significantly changed between V1 and V3 in women who were in remission at V3. Most of these molecules are known to be modulated by the balance between lactobacilli and BV-related bacteria, such as lactate, which increased with remission, and biogenic amines, short-chain fatty acids, acetate, and succinate, which decreased in parallel (6, 45–48).
Concentration of metabolites which significantly varied (P < 0.01) between V1 and V3 in women going into remission at V3 (n = 29)
A complete overview of all samples considered in the present investigation could be obtained by means of PCA, as shown in Fig. 5A. The first PC, representing 99.1% of the total variance among the samples, was dominated by lactate, while the second, representing 0.5% of the total variance, brought information mainly on acetate and succinate, molecules produced by BV-related bacteria. PCA shows 3 types of trends with significant differences registered on PC2: (i) when the remission condition persisted at V4, the medians from V3 and V4 samples were superimposed and were very close to those of healthy samples and different from those of V1 samples; (ii) the average metabolome of women going into remission at V3 but presenting relapse at V4 showed a high variation at V3 with a return to the baseline at V4; and (iii) the metabolome of women who were not in remission at V3 underwent minor modifications compared to women going into remission, but these variations were still significant along PC2 (P < 0.05).
(A) Biplot of a PCA calculated on the entire set of NMR spectra. Median values of the samples collected from healthy women (triangle) and BV-affected women in remission at both V3 and V4 (circles), in remission at V3 but not at V4 (squares), and not in remission at V3 (no symbols) are shown. (B) Biplot of a PCA calculated on NMR data related to BV-affected women at visit V1 and to healthy women. R and N indicate the median values of women grouped according to the response to rifaximin (RRifa, n = 25; NRifa, n = 44) or placebo (RPbo, n = 4; NPbo, n = 19) at V3, while H indicates the median values of healthy women (H; n = 10). Continuous lines in panel A indicate significant differences between V1–V3 and V3–V4 trajectories, while dashed lines indicate nonsignificant differences. Full V4 symbols indicate significant differences between V1 and V4. expl.var, explained variance.
Correlation between baseline metabolome and remission.Figure 5B illustrates the correlation between remission at V3 and vaginal fluid metabolome at V1. Treatment with placebo showed a positive correlation between lactate concentration and the possibility to go into remission, while the opposite was observed when rifaximin was employed. These data, although not statistically significant, suggest that a high lactate concentration is beneficial for spontaneous remission, while the success of rifaximin treatment is independent of lactate presence. In analogy with the considerations on the microbiota, the high lactate concentration at baseline of women going into spontaneous remission suggests that Amsel and Nugent criteria could be not exhaustive for a correct diagnosis of BV.
Correlation between placebo treatment and BV-related metabolites.To study in detail the relationship between BV and short-chain fatty acids, organic acids, and biogenic amine concentration, a PC model was built on the concentrations of these molecules, as shown in Fig. 6. A picture coherent with clinical and microbiological data seemed to emerge. In the case of placebo treatment (Fig. 6A), the highest concentrations of BV-related molecules were found in V1 samples of women still affected by BV at V3, while the lowest concentrations of these metabolites were detected in V1 samples of women who maintained remission for the long term. Women in remission at V3 but still affected by BV at V4 showed an intermediate situation. This means that a connection between the position of the samples along PC1 at V1 and the possibility to go into remission could be noticed. However, the treatment did not cause significant modifications to the metabolome, so the samples collected at V3 and V4 appeared along PC1 very close to the samples collected at V1.
Biplot of a PCA calculated on the autoscaled concentrations of the molecules associated with the presence of BV-related bacteria measured in healthy women (triangle) and BV-affected women in remission at both V3 and V4 (circles), in remission at V3 but not at V4 (squares), and not in remission at V3 (no symbol). To limit the sample superimposition, data are represented separately for women subjected to placebo (A) or rifaximin (B) treatment. Median values of the samples collected from healthy women (H) and BV-affected women at V1, V3, and V4 are shown. The number of women (n) for each sample group is reported in Fig. 1. Continuous lines indicate significant differences between V1–V3 and V3–V4 trajectories, while dashed lines indicate nonsignificant differences. Full V4 symbols indicate significant differences between V1 and V4. expl.var, explained variance.
Correlation between rifaximin treatment and BV-related metabolites.In women treated with rifaximin (Fig. 6B), clinical remission and microbial restoration were associated with a significant reduction of the concentration of the molecules produced by BV-related bacteria (acetate, succinate, butyrate, putrescine, and tyramine). V4 samples collected from women who maintained remission were very close to V3 samples, while a relapse resulted in a return of V4 samples very close to V1 samples. A slight but still significant reduction in concentration of molecules produced by BV-related bacteria also could be observed in women failing to go into remission at V3.
Effects of different dosages of rifaximin on the vaginal microbiota.Figure 7A shows the scoreplot of a PCA calculated on the differences existing for each woman between the vaginal microbiota composition at V1 and V3. The two PCs accounted for 46.4% of the entire variability. The proximity of results for placebo (RD) and those for women who had still BV at V3 (N) indicates that spontaneous remission is not associated with modifications induced in the microbial communities by placebo, as changes are minor and similar to those determined by an ineffective treatment. This finding further confirms the association between spontaneous remission and favorable microbiota at the baseline. Differences in the microbiota were more marked for rifaximin treatment groups A, B, and C, with the highest variations registered for group B.
Biplots of PCA models built on the autoscaled differences for V3–V1 between qPCR (A) and metabolic (B) data. The medians of the samples from women in remission or not in remission at V3 are indicated as R and N, respectively, while rifaximin dosing regimens are indicated as subscripts. Abbreviations: Lact, Lactobacillus; Gard, G. vaginalis; Atop, Atopobium; Prev, Prevotella; Veil, Veillonella; Mhom, M. hominis; Mobi, Mobiluncus; expl.var, explained variance.
Effects of different dosages of rifaximin on the vaginal metabolome.The effects of the different dosages of rifaximin on the metabolome are shown in Fig. 7B. Very similar modifications were found between V1 and V3 in women in spontaneous remission at V3 induced by placebo (RD) and in women who still had BV at V3 (N). Major differences from ineffective treatment were observed for treatment B, followed by treatments C and A. The loadingplot superimposed on the scoreplot shows that treatment B induced the highest increase of lactate and proline, which was associated with the highest decreases of organic acids and biogenic amines.
DISCUSSION
An understanding of the composition and metabolic activities of the vaginal microbial ecosystem in relation to a vaginal disease is essential for comprehensively understanding the etiology of that disease and for rationally designing successful strategies of prevention and therapy. Although BV has been studied in detail regarding the implications of the vaginal microbiota composition (42–44), the vaginal metabolomic framework remains largely unexplored: to date, only one gas chromatography-mass spectrometry-based study has been focused on the vaginal metabolome in BV (49), but no data are available to correlate the metabolomics features with the potential remission induced by antibiotic therapy.
To study BV from a holistic point of view, a combination of two elements are necessary: (i) analytical techniques able to account for the complexity of the vaginal ecosystem, and (ii) multivariate statistical techniques able to bring out the general features of the vaginal samples from the confounding sources of variability that always characterize complex microbial ecosystems.
In the present work, the effects of rifaximin vaginal tablets on vaginal microbiota and metabolome of BV-affected women were investigated by combining qPCR, a well-established technique for the characterization of vaginal bacterial communities (25, 31, 50), and a metabolomic approach based on 1H-NMR, which had never been used to study vaginal ecology. To highlight the general trends of the bacterial communities and metabolomic profiles in response to the antibiotic/placebo therapy, a multivariate statistical strategy was set up based on the trajectories traced by vaginal samples in a PCA space.
In general, our data demonstrated the efficacy of rifaximin in restoring a health-like condition in terms of both bacterial communities and metabolomics features, in agreement with the clinical observations previously reported (11). In particular, a significant increase of the lactobacillus/BV-related bacteria ratio was found in women affected by BV and treated with rifaximin. This increase was observed not only in women in remission but also in women not in remission, even if the trajectory traced in the PC space built on the qPCR data was shorter for the latter.
Metabolomic data showed that the restoration of a Lactobacillus-dominated flora, which occurred following rifaximin treatment, was strongly associated with the increase of lactic acid concentration and the decrease of the pool of molecules associated with the metabolism of BV-related bacteria, such as acetic acid, succinate, short-chain fatty acids, and biogenic amines (6, 45–48). Placebo treatment did not cause any significant variation in the metabolome in cases of spontaneous remission and BV persistence. In addition, the PCA approach developed in the present work allowed us to visualize a BV relapse occurring at the second follow-up visit as a return of the vaginal samples toward the baseline situation in the PC spaces built on both qPCR and metabolomic data.
An interesting finding of the present study is the hypothesis of a spontaneous remission from BV on the basis of the vaginal microbial and metabolomic features present at the baseline. In fact, women going into remission after placebo treatment showed at the baseline the highest ratios of lactobacilli to BV-related bacteria and lactate to BV-related metabolites. These women also could be considered not to be true cases of BV; thus, antibiotic treatment would not be appropriate. On the contrary, the use of rifaximin treatment appeared to be useful to induce remission from BV in women who presented both unfavorable microbiomes and metabolomes. This observation suggests the use of molecular or metabolomics tools to predict the need for a specific antibiotic treatment to induce remission from BV.
The impact of different dosages of rifaximin on the vaginal microbiota and metabolome was also investigated in the present study. The major differences from ineffective treatment were induced by 25 mg of rifaximin for 5 days, in agreement with the findings of Cruciani et al. (12), who described this dosage as the best treatment regimen in terms of restoration of a healthy-like microbiota.
The metabolome of a biological system, being downstream of the genome, transcriptome, and proteome, is particularly sensitive to external agents, even to those determining little effects on the different metabolic fluxes. NMR is ideally tailored for metabolomic investigations due to its high reproducibility, easing the application of pattern recognition techniques and multivariate analyses on minimally optimized spectra (51). A strategy based on the combination of microbiological and metabolomic observations allowed a comprehensive study of the effects of the antibiotic rifaximin on the vaginal environment of women affected by BV. Relationships were found between the doses of the antibiotic agent and restoration of healthy conditions, as well as between initial conditions and evolution of the disease. This allows us to envisage possibilities to finely tailor treatment strategies against BV according to vaginal environmental characteristics.
ACKNOWLEDGMENTS
This work was supported by a research grant provided by Alfa Wassermann S.p.A. (no. 15/2012). G. Donders received advisory fees, lecture fees, and grant support from Alfa Wassermann S.p.A. F. Calanni is an employee of Alfa Wassermann S.p.A.
FOOTNOTES
- Received 5 February 2014.
- Returned for modification 19 March 2014.
- Accepted 28 March 2014.
- Accepted manuscript posted online 7 April 2014.
Supplemental material for this article may be found at http://dx.doi.org/10.1128/AAC.02469-14.
- Copyright © 2014, American Society for Microbiology. All Rights Reserved.
