Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2007, p. 5731-5741, Vol. 73, No. 18
0099-2240/07/$08.00+0 doi:10.1128/AEM.01251-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Pharmaceutical Sciences, CIRB-Center for Biotechnology, University of Bologna, Bologna, Italy,1 Femicare Clinical Research for Women, Tienen, Belgium,2 Departments of Obstetrics and Gynecology, General Hospital Heilig Hart, Tienen, and University Hospital Gasthuisberg, Catholic University of Leuven, Leuven, Belgium3
Received 6 June 2007/ Accepted 12 July 2007
|
|
|---|
|
|
|---|
The characterization of the composition and ecology of the vaginal microbial ecosystems in healthy women and in patients affected by infectious diseases, such as BV and CA, is important for an understanding of the etiology of these diseases and the development of efficient therapies. BV, widely believed to be the most common vaginal disorder affecting women of reproductive age, is defined as an alteration of the vaginal bacterial morphotypes characterized by an overgrowth of several anaerobic and microaerophilic bacteria, gram-positive cocci, and a Mycoplasma sp. and by a decreased prevalence of Lactobacillus species (17, 23). Incidence rates of BV range from 5 to 50% among both nonpregnant and pregnant women (22). CA is an infection by Candida spp., mostly Candida albicans. Yeast overgrowth can modify the normal vaginal flora. Up to 75% of women experience genital CA during their lifetime, and 5 to 8% have chronic, recurring CA, defined as four or more episodes in a 12-month period (19, 33).
In recent years, molecular techniques based on the analysis of rRNA gene sequences have been developed, providing powerful tools to reveal the phylogenetic diversity of the microorganisms found within complex ecosystems and understand community dynamics (12, 24, 41). PCR-denaturing gradient gel electrophoresis (PCR-DGGE) represents a rapid and reliable technique that has been used successfully to identify the bacterial compositions of different ecological niches, including the vaginal ecosystem (9, 10, 13, 14, 35). Sequencing of different 16S rRNA gene fragments allows the determination of which bacterial species are most common among the human specimens. Real-time PCR with genus- and species-specific primers can provide a quantification of selected bacteria present in complex microbial communities by measuring the amount of PCR products in each cycle as fluorescence (6, 11, 26, 27, 31).
In the present study, we utilized PCR-DGGE with universal primers for the bacterial 16S rrn operon to explore the composition of normal vaginal communities and assess changes related to BV and CA. In addition to using this qualitative approach, we quantified changes in Lactobacillus populations among women with different clinical profiles by using real-time PCR with 16S rRNA gene-targeted genus-specific primers.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Structure of vaginal microflora as analyzed by PCR-DGGE and real-time PCR
|
PCR-DGGE analysis.
Amplification of the V2-V3 region of the bacterial 16S rRNA gene was carried out using the universal eubacterial primers HDA1-GC (containing a GC clamp) and HDA2 (38), supplied by M-Medical, Milan, Italy. The amplification reactions were performed in a T Gradient thermal cycler (Biometra, Göttingen, Germany). Dynazyme II (Celbio, Milan, Italy) was used as a thermostable polymerase under conditions suggested by the supplier. The reaction mixture contained 0.5 µM of each primer, 200 µM of each deoxynucleoside triphosphate, 0.5 U of Dynazyme II, and 4 µl of the bacterial DNA template in a final volume of 50 µl. The thermal cycling program consisted of the following time and temperature profile: 95°C for 5 min; 30 cycles of 95°C for 30 s, 56°C for 30 s, and 72°C for 60 s; and 72°C for 8 min. Samples (5 µl) of the amplified products (200 bp) were subjected to gel electrophoresis in 2% agarose gels and visualized by ethidium bromide staining.
DGGE analysis was performed using the D-Code universal mutation detection system apparatus (Bio-Rad, Hercules, CA) with 20-cm by 20-cm by 0.75-mm gels. The sequence-specific separation of the PCR fragments was obtained in 8% (wt/vol) polyacrylamide (acrylamide-N,N'bisacrylamide; 40:3 [wt/vol]) gels in 0.5 x TAE buffer (20 mM Tris, 10 mM glacial acetic acid, 0.5 mM EDTA, pH 8). The denaturing gels contained a 30% to 50% gradient of urea and formamide increasing in the direction of electrophoresis. A 100% denaturing solution contained 40% (vol/vol) formamide and 7 M urea. A stacking gel containing 8% (wt/vol) polyacrylamide was applied onto the denaturing gel. A volume of 8 to 16 µl of PCR samples was loaded onto the stacking gel. Electrophoresis was conducted at a constant voltage of 130 V and a temperature of 60°C for approximately 6 to 7 h. Following electrophoresis, the gel was silver stained by using the protocol of Bassam et al. (7).
Sequencing of the V2-V3 region of the 16S rRNA gene.
DNA fragments of interest were excised from the denaturing gels with a sterile scalpel, washed once in 1 x PCR buffer, and incubated in 20 µl of the same buffer overnight at 4°C. Four microliters of the buffer solution was used as the template for PCR. Reamplification of the V2-V3 region was conducted as described above by employing the primers HDA1 (without the 5' GC clamp) and HDA2. The reamplified fragments were purified by using the QIAquick PCR purification kit (QIAGEN) and then subjected to automated sequence analysis of both DNA strands with HDA1 and HDA2 primers. BigDye terminators (ABI-PerkinElmer, Foster City, CA) were used with a PRISM 377 sequencer (ABI). The sequence identities were determined by comparison with rRNA gene sequences deposited in the GenBank database by using the BLAST algorithm (2).
Real-time quantitative PCR.
Quantitative PCR was performed with a LightCycler instrument (Roche, Mannheim, Germany), and SYBR green I fluorophore was used to correlate the amount of PCR product with the fluorescence signal. Total bacterial DNA extracted from vaginal fluids was amplified with the universal primer set HDA1/HDA2 (38) and the Lactobacillus genus-specific primer set Bact-0011f/Lab-0677r (21), which amplifies a 700-bp region inside the 16S rRNA genes. Three subsamples of each DNA extract were amplified in a final volume of 20 µl containing 4 mM of MgCl2, 0.5 µM of each primer, 2 µl of LightCycler-FastStart DNA Master SYBR green I (Roche), and 2 µl of either template or water (no-template control). The thermal cycling conditions used were as follows: an initial denaturation step at 95°C for 10 min followed by 40 cycles of denaturation at 95°C for 15 s, primer annealing at 56°C (HDA1/HDA2) or 63°C (Bact-0011f/Lab-0677r) for 25 s, extension at 72°C for 40 s, and an additional incubation step at 85°C for 5 s for fluorescence acquisition. For each step, the temperature transition rate was 20°C s–1. Melting curve analysis for PCR product identification was obtained by heating at 20°C s–1 to 95°C, cooling at 20°C s–1 to 60°C with a 15-s hold, and then slowly heating at 0.2°C s–1 to 99°C. Fluorescence readings were collected continuously during this heating to monitor the denaturation of PCR products.
Quantification of rrn operons of total eubacteria and lactobacilli was done by using standard curves made with known concentrations of Lactobacillus acidophilus NCFM genomic DNA, which has a size of 1.994 Mb and four rrn operons (1). Chromosomal DNA of L. acidophilus NCFM was extracted by using the DNeasy tissue kit (QIAGEN) and serially diluted from 105 to 102 molecules µl–1. Results obtained by PCR were converted to the average estimate of total eubacterial and Lactobacillus sp. rrn operons present in 1 ml of vaginal rinsing, and standard deviations (SD) were calculated. The ratio of Lactobacillus sp. rrn operons to total eubacterial 16S rRNA gene copies (relative abundance) was evaluated for each vaginal sample.
Statistical analysis.
Differences in the numbers of total eubacterial and lactobacillus rrn operons and relative abundances of lactobacilli among NI, BV, and CA groups were analyzed. Data were assessed for normal distribution by using the Kolmogorov-Smirnov test, and the Kruskal-Wallis test was performed to simultaneously compare the analysis groups. When significant differences (P < 0.05) were obtained, multiple-comparison Dunn's procedure was used to isolate the groups that differed from the others. Statistical analysis was carried out with SigmaStat software (Systat Sofware Inc., San Jose, CA).
Nucleotide sequence accession numbers.
The nucleotide sequences of the V2-V3 DGGE fragments have been deposited in the DDBJ nucleotide sequence database under accession numbers AB279893 to AB279983.
|
|
|---|
![]() View larger version (114K): [in a new window] |
FIG. 1. DGGE profiles of healthy women (NI group). A single DGGE profile was reported for each woman in the NI group. The bands correspond to bacterial taxa and GenBank accession numbers as follows: Lactobacillus iners, 3a (AB279893), 6b (AB279894), 8a (AB279895), and 10c (AB279896); Lactobacillus sp., 5a (AB279907) and 26a (AB279908); Lactobacillus acidophilus, 6a (AB279917), 10b (AB279918), 11a (AB279919), 22a (AB279920), 23a (AB279921), and 25b (AB279922); Lactobacillus gasseri, 10a (AB279926), 14a (AB279927), and 25a (AB279928); Lactobacillus vaginalis, 10d (AB279933) and 22b (AB279934); Gardnerella vaginalis, 5b (AB279938), 14c (AB279939), and 25c (AB279940); uncultured Gardnerella sp., 14b (AB279952).
|
![]() View larger version (96K): [in a new window] |
FIG. 4. DGGE profiles of patients developing both BV and CV (BV-CA). Three DGGE courses corresponding to the three monthly gynecologic visits were reported. The bands correspond to bacterial taxa and GenBank accession numbers as follows: Lactobacillus iners, 7e (AB279904), 9a (AB279905), and 15b (AB279906); Lactobacillus sp., 15a (AB279916); Gardnerella vaginalis, 7h (AB279949), 7i (AB279950), and 9b (AB279951); Leptotrichia amnionii, 7a (AB279958) and 7b (AB279959); Leptotrichia sanguinegens, 7c (AB279961); uncultured Prevotella sp., 7d (AB279966) and 7f (AB279967); uncultured Clostridium sp., 7g (AB279975).
|
L. iners, L. acidophilus, Lactobacillus plantarum, L. gasseri, and L. vaginalis remained the predominant species in the vaginas of the women in the BV subject group when the clinical evaluation registered an NI (Fig. 2). The PCR-DGGE patterns of women during visits during which BV was diagnosed (BV visits) showed a number of dominant fragments that was higher than that for women during NI visits, indicating that high bacterial diversity is associated with BV status, in accordance to what has been previously reported (9, 14, 36, 40). Generally, we observed a correlation between the BV condition and a relative lack of lactobacilli, accompanied by a profound increase in the quantity of other vaginal anaerobic bacteria related to Gardnerella vaginalis, Atopobium vaginae, Leptotrichia amnioni, Leptotrichia sanguinegens, Prevotella sp., Megasphaera sp., an uncultured Chloroflexi bacterium, an uncultured Clostridium sp., Streptococcus sp., Staphylococcus sp., and Veillonella montpellierensis. Subject 2, who had the clinical symptoms for BV at every visit, presented the most complex and diversified vaginal microbiota. Minor variations in DGGE profiles were observed between samples from BV and NI visits for subjects 12, 17, 20, and 24, for whom BV was diagnosed in one visit only. Notably, Gardnerella vaginalis was found in all vaginal samples corresponding to BV visits except for the third visit for subjects 20 and 12.
![]() View larger version (144K): [in a new window] |
FIG. 2. DGGE profiles of patients developing BV. Three DGGE courses corresponding to the three monthly gynecologic visits were reported. The bands correspond to bacterial taxa and GenBank accession numbers as follows: Lactobacillus iners, 17b (AB279897), 20e (AB279898), and 24b (AB279899); Lactobacillus acidophilus, 17a (AB279923) and 20d (AB279924); Lactobacillus plantarum, 20a (AB279937); Lactobacillus gasseri, 12a (AB279929) and 20b (AB279930); Lactobacillus vaginalis, 17c (AB279935); Atopobium vaginae, 2q (AB279953) and 12g (AB279954); Gardnerella vaginalis, 2m (AB279941), 2o (AB279942), 2p (AB279943), 17d (AB279944), 17e (AB279945), and 24c (AB279946); Leptotrichia amnionii, 2a (AB279955), 2b (AB279956), and 2d (AB279957); Leptotrichia sanguinegens, 24a (AB279960); uncultured Prevotella sp., 2c (AB279962), 2f (AB279963), 2g (AB279964), and 2h (AB279965); Prevotella sp., 12b (AB279968) and 12e (AB279969); uncultured Megasphaera sp., 2l (AB279970), 20f (AB279971), and 20g (AB279972); uncultured Chloroflexi bacterium, 2i (AB279973); uncultured Clostridium sp., 2n (AB279974); Streptococcus sp., 12c (AB279977); Staphylococcus sp., 2e (AB279978) and 20c (AB279979); Veillonella montpellierensis, 12d (AB279980) and 12f (AB279981).
|
![]() View larger version (98K): [in a new window] |
FIG. 3. DGGE profiles of patients developing CA. Three DGGE courses corresponding to the three monthly gynecologic visits were reported. The bands correspond to bacterial taxa and GenBank accession numbers as follows: Lactobacillus sp., 1a (AB279909), 1b (AB279910), 1c (AB279911), 1d (AB279912), 4b (AB279913), 13a (AB279914), and 19a (AB279915); Lactobacillus iners, 4c (AB279900), 16c (AB279901), 18a (AB279902), and 21b (AB279903); Lactobacillus gassseri, 16a (AB279931) and 21a (AB279932); Lactobacillus acidophilus, 16b (AB279925); Lactobacillus vaginalis, 4d (AB279936); Gardnerella vaginalis, 21c (AB279947) and 21d (AB279948); uncultured Ureaplasma sp., 4a (AB279982); uncultured Clostridium sp., 4e (AB279976); Shigella boydii, 18b (AB279983).
|
Multiple fragments of different sizes were sometimes visible in the DGGE gels (for example, bands 1a, 1b, 1c, and 1d of Fig. 3). After the minor DNA fragments were sequenced, they were found to be homologous to the species corresponding to the closest major fragment. Similar to findings reported by other authors (10, 38), the major fragment indicated the identity of the isolate, whereas the minor ones were probably PCR artifacts resulting from the highly folded (loops and stems) structure of V2-V3 region.
Figure 5 shows the frequency of the occurrence (in percentages) of the vaginal bacteria, detected by PCR-DGGE, in relation to the clinical status. Samples 7III and 9I were not included in this analysis and in the subsequent statistical evaluations because of the simultaneous presence of BV and CA conditions. The histograms indicate the increased complexity of vaginal ecology under BV conditions, characterized by the appearance of several anaerobic species and the decrease of L. vaginalis, L. gasseri, L. acidophilus, and Lactobacillus sp. The figure also shows that the composition of the vaginal bacterial flora did not undergo significant changes as a consequence of Candida infection, even if a diverse distribution of Lactobacillus species appeared. In particular, we observed increased frequencies of L. iners and Lactobacillus sp. and decreased frequencies of L. acidophilus, L. gasseri, and L. vaginalis. Notably, H2O2 is commonly produced by strains of L. acidophilus, L. gasseri, and L. vaginalis, while it is an uncommon metabolite among strains of L. iners (5, 34, 41).
![]() View larger version (33K): [in a new window] |
FIG. 5. Frequency of the occurrence of bacterial species in relation to clinical status of the subject. Frequencies were calculated as the ratios (percentage) of the number of visits characterized by the presence of each species to the total number of visits belonging to each visit group (NI, BV, and CA).
|
To test the significance of quantitative changes in Lactobacillus populations related to BV and CA infections, the numbers of 16S rRNA gene copies of total eubacteria and lactobacilli and the relative abundance of lactobacilli in samples from women during NI, BV, or CA visits (plotted in Fig. 6) were analyzed by a Kruskal-Wallis test due to nonnormal data distribution. The only significant difference among the groups was found for relative abundances of lactobacilli (P < 0.001). Dunn's procedure confirmed the existence of a significant difference in relative abundances of lactobacilli (P < 0.05) for each possible comparison (BV versus NI, CA versus NI, and CA versus BV). Figure 6C highlights the reduction of the relative abundance of lactobacilli in the vaginal flora of women during BV visits compared to those in the vaginal flora of women during NI and CA visits and the increase in the same parameter for women during CA visits compared to women during NI visits.
![]() View larger version (12K): [in a new window] |
FIG. 6. Real-time PCR evaluation of 16S rrn operons of total eubacteria (A) and lactobacilli (B) and relative abundance of lactobacilli (C) related to each group of visits characterized by a specific clinical status (NI, BV, and CA). The number of operons in 1 ml of vaginal rinsing is expressed as an absolute value. Lactobacillus relative abundance is expressed as the ratio of lactobacillus rrn operons to total eubacterial 16S rrn operons. The box for each group represents the interquartile range (25th to 75th percentile), and the line within this box is the median value. The bottom and top bars indicate the 10th and 90th percentiles, respectively. Outlier values are indicated (black circles).
|
![]() View larger version (13K): [in a new window] |
FIG. 7. Relative abundances of lactobacilli in the vaginal flora of women during NI visits, relative to the clinical history of the women. Lactobacillus relative abundance is expressed as the ratio of lactobacillus 16S rrn operons to total eubacterial 16S rrn operons. The box for each group represents the interquartile range (25th to 75th percentile), and the line within this box is the median value. The bottom and top bars indicate the 10th and 90th percentiles, respectively. Outlier values are indicated (black circles).
|
|
|
|---|
The present study provides basic information concerning the structure of bacterial communities in healthy vaginas and qualitative/quantitative changes in numerically dominant bacterial populations related to BV and CA diseases. Particular attention was paid to dynamics of lactobacilli. Two molecular techniques, PCR-DGGE and real-time PCR, were combined to achieve a detailed picture of vaginal ecology under normal and abnormal conditions. This integrated approach has been recently proposed for monitoring the effect of prebiotics, probiotics, and symbiotics on the fecal microbiota of healthy humans (37).
Consistent with data from previous studies (4, 10, 12, 24, 41), our DGGE data showed that the vaginal bacterial floras in healthy women were stable over time and dominated by a single species or group of closely related species of Lactobacillus (L. iners, L. acidophilus, L. gasseri, and L. vaginalis). The coexistence of multiple species of lactobacilli is rare because of the competitive exclusion of one species by another, preemptive colonization by a particular species or host factors that strongly influence which species are able to colonize the environment (41). However, in samples from a small number of healthy women, G. vaginalis was detected in conjunction with Lactobacillus species. This result suggests that the presence of lactobacilli may not exclude potential pathogens from the vagina; likewise, the presence of G. vaginalis is not strictly related to clinically manifested disturbances. In healthy vaginas, G. vaginalis could represent a "sentinel species," a term used to refer to an indigenous species that is particularly sensitive to changes in biological, physical, or chemical characteristics of the environment and responds through an increase in population size (12). The identification of potential sentinel species can help to predict whether there is an increased risk of an individual contracting an infection due to shifts in the composition of the microbial communities or changes in the abundance of specific populations. The sentinel species concept may be also used for prevention and early diagnosis of BV disturbance, which is characterized by the outgrowth of several potentially pathogenic microorganisms normally present but at much lower numbers (9, 14, 20, 40). Actually, the patients affected by BV enrolled in this clinical trial showed a relative scarcity of lactobacilli, together with an increased complexity of vaginal communities dominated by anaerobic bacteria related to Gardnerella vaginalis, Atopobium vaginae, Leptotrichia sp., Prevotella sp., Megasphaera sp., Chloroflexi, Clostridium sp., Streptococcus sp., Staphylococcus sp., and Veillonella sp. All these organisms have already been identified as components of the abnormal vaginal flora under conditions of BV disturbance (9, 20), except for Chloroflexi, Streptococcus sp., and Veillonella sp., which represent newly recognized populations not previously associated with BV. By contrast, we did not find an abnormal bacterial vaginal ecosystem in patients affected by CA, as the predominant presence of lactobacilli was observed with our DGGE-based qualitative analysis. However, a different frequency of particular Lactobacillus species was found in respect to the frequency detected under healthy conditions. In particular, we observed a decrease in H2O2-producing species (L. acidophilus, L. gasseri, and L. vaginalis) and an increase of non-H2O2-producing L. iners. This finding suggests some hypotheses regarding the role of lactobacilli in protecting women from Candida infection: (i) women with a vaginal flora dominated by H2O2-producing Lactobacillus species have a minor risk to contract CA; (ii) H2O2, rather than other metabolites produced by lactobacilli, such as lactic acid, could be responsible for the control of Candida overgrowth (33). In patients prone to developing CA as well as BV, the DGGE pattern of CA-associated lactobacilli was dominant, while the DGGE pattern of BV-associated bacteria was not detected in this group with mixed infections.
An important advantage of the DGGE technique is the ability to detect several bacterial populations not readily culturable. For example, L. iners does not grow on the major selective media for isolation of Lactobacillus, including MRS and Rogosa media (18). Likewise, the strict anaerobes Atopobium vaginae, Megasphaera, and Leptotrichia require specialized media and often grow slowly. The finding of these organisms as members of vaginal flora demonstrates how cultivation-based methods can be misleading.
Since DGGE can be considered a semiquantitative tool for monitoring bacterial populations, additional analysis with real-time PCR was carried out to obtain a quantitative estimation of changes in lactobacillus concentrations associated with BV and CA diseases. An analogous approach was used by Bartosch et al. (6) to monitor changes in intestinal bacterial community structure among different groups of elderly people. A statistically significant decrease in Lactobacillus relative abundance during BV visits was registered compared to the relative abundances during NI and CA visits, confirming the qualitative observations derived from DGGE analysis. Notably, the relative abundance of lactobacilli was significantly increased for women during CA visits, further supporting our hypothesis that the predisposition to Candida infection is not correlated with the predominance of lactobacilli but to the presence of particular Lactobacillus species. We also investigated the possibility of a link between the clinical history and the relative abundance of lactobacilli in vaginal fluids of women during NI visits, but no significant correlations were detected.
In conclusion, it has been demonstrated for the first time that the integrated use of DGGE and real-time PCR represents a successful strategy to monitor the qualitative and quantitative changes of vaginal bacterial communities correlated to different infectious diseases, such as BV and CA. In particular, the major findings obtained in this study are the following: (i) Lactobacillus species are the major constituents in vagina of healthy women and patients developing CA but not in patients developing BV; (ii) new components of vaginal microbiota in BV have been discovered; (iii) women prone to developing CA mainly host non-H2O2-producing Lactobacillus species; and (iv) in patients developing a mixed BV and CA infection, the colonization pattern of lactobacilli resembles the CA pattern.
Published ahead of print on 20 July 2007. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»