Characterization of the vaginal microbiota of healthy Canadian women through the menstrual cycle
© Chaban et al.; licensee BioMed Central Ltd. 2014
Received: 13 December 2013
Accepted: 31 March 2014
Published: 4 July 2014
The vaginal microbial community plays a vital role in maintaining women’s health. Understanding the precise bacterial composition is challenging because of the diverse and difficult-to-culture nature of many bacterial constituents, necessitating culture-independent methodology. During a natural menstrual cycle, physiological changes could have an impact on bacterial growth, colonization, and community structure. The objective of this study was to assess the stability of the vaginal microbiome of healthy Canadian women throughout a menstrual cycle by using cpn 60-based microbiota analysis. Vaginal swabs from 27 naturally cycling reproductive-age women were collected weekly through a single menstrual cycle. Polymerase chain reaction (PCR) was performed to amplify the universal target region of the cpn 60 gene and generate amplicons representative of the microbial community. Amplicons were pyrosequenced, assembled into operational taxonomic units, and analyzed. Samples were also assayed for total 16S rRNA gene content and Gardnerella vaginalis by quantitative PCR and screened for the presence of Mollicutes by using family and genus-specific PCR.
Overall, the vaginal microbiome of most women remained relatively stable throughout the menstrual cycle, with little variation in diversity and only modest fluctuations in species richness. Microbiomes between women were more different than were those collected consecutively from individual women. Clustering of microbial profiles revealed the expected groupings dominated by Lactobacillus crispatus, Lactobacillus iners, and Lactobacillus jensenii. Interestingly, two additional clusters were dominated by either Bifidobacterium breve or a heterogeneous mixture of nonlactobacilli. Direct G. vaginalis quantification correlated strongly with its pyrosequencing-read abundance, and Mollicutes, including Mycoplasma hominis, Ureaplasma parvum, and Ureaplasma urealyticum, were detected in most samples.
Our cpn 60-based investigation of the vaginal microbiome demonstrated that in healthy women most vaginal microbiomes remained stable through their menstrual cycle. Of interest in these findings was the presence of Bifidobacteriales beyond just Gardnerella species. Bifidobacteriales are frequently underrepresented in 16S rRNA gene-based studies, and their detection by cpn 60-based investigation suggests that their significance in the vaginal community may be underappreciated.
KeywordsVaginal microbiome cpn 60 Menstrual cycle Bifidobacteriales Bifidobacterium Lactobacillus Mollicutes Mycoplasma Ureaplasma Gardnerella Vaginal bacteria
It has long been recognized that the microbial community of the lower genital tract plays a vital role in maintaining the reproductive health of women. The vaginal microbiota of reproductive-aged women has traditionally been characterized by culture-based techniques as dominated by Lactobacillus species, which, among other roles, produce lactic acid, biosurfactants, hydrogen peroxide, and other factors that create an inhospitable environment for pathogenic bacteria[2–7]. Detailed community profiling with culture-independent techniques has demonstrated that “healthy” microbial communities are usually dominated by Lactobacillus crispatus, Lactobacillus iners, Lactobacillus gasseri, Lactobacillus jensenii, or a combination of these species, and in a small portion of women, “mixed” profiles are depleted of lactobacilli and consist of bacteria such as Gardnerella, Prevotella, Atopobium, Megasphaera, and Streptococcus[8–14]. The latter has challenged our traditional understanding of the “healthy vaginal microbiome,” and raises questions about the structure and function of this community and the host response to it.
The relative compositional stability of the vaginal microbiome is quite remarkable, given the variability in the host ecosystem associated with the menstrual cycle, sexual contact, and introduction of bacteria from the skin and external environment. In particular, the menstrual cycle creates an ever-changing vaginal environment, with ovulation, menses, and corresponding fluctuations of estrogen and progesterone levels affecting bacterial attachment to the vaginal epithelium, cervical mucus production, pH and redox potential, and glycogen levels. Culture-based studies in which the microbial community is characterized at several time points over the cycle have reported a range of findings, from no change in community composition, to some variation in the aerobic versus anaerobic communities over a menstrual cycle, and a greater proportion of non-Lactobacillus species present during menses[19–23]. This interindividual variability is mirrored in culture-independent studies, with some women maintaining a consistent microbial community throughout multicycle sampling, others having fluctuations timed with menses, and some having random fluctuations with no apparent cause[13, 24, 25]. Interestingly, in at least some cases, the overall functional characteristics of the community are predicted to be maintained, despite the fluctuations in the bacterial composition, because the shifts in relative dominance may be limited to different lactic acid-producing bacterial species.
Our understanding of the vaginal microbiota to date had been predominantly shaped by culture-based or 16S rRNA gene culture-independent studies. Although both approaches contribute a wealth of information on the microbiota composition, they have limitations. The move to culture-independent studies was spurred by the labor-intensive methodologies associated with bacteriologic culture (impractical for large studies), and limited ability to grow the diverse array of organisms present, as well as substantial difficulty in accurate determination of relative abundance. 16S rRNA gene-based methods have overcome some of the limits of culture, but are known to have amplification biases and offer limited resolution for some taxa. Alternative molecular targets, like the universal cpn 60 gene, have been used to gain a different perspective of microbial communities in a culture-independent fashion.
Similar to 16S rRNA gene-based studies, cpn 60 surveys of the vaginal microbiome reveal many largely lactobacilli-dominated profiles, and better taxonomic resolution of some nonlactobacilli groups like Gardnerella and Prevotella[28–31]. For future studies of the vaginal microbiome, it is very important to understand whether timing of sampling during the menstrual cycle would result in significant variability that must be accounted for in cross-sectional studies of the vaginal microbiome in populations.
The primary objective of this study was to characterize the vaginal microbiome over a natural menstrual cycle among a cohort of healthy, asymptomatic, Canadian women, by using the cpn 60 gene target. Secondarily, we aimed to probe for Mollicutes and Bifidobacteria not well detected by other high-throughput sequencing methods, and to better understand the temporal stability and/or variation of the vaginal microbiome within individual women.
Participants and study design
This longitudinal study was designed to collect vaginal-swab samples from a cohort of healthy women weekly over a single menstrual cycle. Women were eligible to participate if they demonstrated comprehension of the English language to a level necessary to provide informed consent, were at least 18 years of age, and had regular menstrual cycles. Individuals were excluded if they were pregnant or planning to become pregnant during the study period, had a chronic autoimmune or inflammatory condition, had an intrauterine device in situ, used hormonal contraceptives, or were currently taking or had taken antimicrobial medications (for example, antibiotic or antifungal therapy) within 4 weeks of enrollment. This study received ethics approval from the Institutional Review Board at the University of British Columbia (certificate no. H09-00860). The sample size was based on previous longitudinal studies of the vaginal microbiome in which a range of seven to 49 subjects provided sufficient numbers for microbiome investigation[13, 25, 32–34].
Data and sample collection
Healthy reproductive-aged women were recruited from two cohorts of participants in Phase-III clinical trials of HPV vaccines, a private Obstetrics and Gynecology practice, and through online and print advertisements placed in Vancouver, British Columbia, Canada. After obtaining informed consent, basic demographic and clinical data were collected by history and from clinical records. While conducting routine speculum examination (usually for pap smear screening), clinicians used a Dacron swab (Copan Diagnostics Inc., Murrieta, CA, USA) to sample the posterior fornix and lateral vaginal wall. This sample represented the first of the four samples collected through the menstrual cycle from charts. All participants were provided with diaries to log activities including sexual activity and any interval symptoms. They were also provided with self-collection kits that contained sterile flocked swabs designed with a break point on the handle (Puritan Medical Products Company LLC, Guilford, ME, USA), collection tubes containing 200 μl of DNAzol-Direct reagent (MRCGENE, Cincinnati, OH, USA) and detailed collection instructions. After hand washing, women were instructed to insert the swab into their vaginas to the half-way point on the handle, rotate the swab 3 times, remove the swab and place it into the provided collection tube, break off the handle, leaving the swab head in the reagent tube, and close and date the tube.
Samples were stored at ambient temperatures until all three self-collected samples were acquired at 1-week intervals. This self-sampling method was duplicative of previous study methods used by ourselves and colleagues in which we validated the high quality of self-sampling compared with clinician sampling[35, 36].
At the end of the collection period, participants returned samples and met with study staff to review their diaries and complete a final interview, in return for a modest honorarium ($20 CAD). All samples were de-identified, and by using information from participants assigned to a menstrual phase by using a calendar-based method: menstrual, day 1 (onset of menstruation) to cessation of bleeding (day 4 to 7); follicular, cessation of bleeding to day 12; periovulatory, day 13 to day 16; luteal, day 17 to days 26 to 32 (commencement of bleeding). If two samples were collected within the luteal phase, they were numbered sequentially as luteal-I and luteal-II. DNAzol-containing nucleic acid was separated from the swab head by centrifugation in the laboratory and used directly as template for PCR reactions.
Quantitative PCR (qPCR) and conventional PCR
Samples were evaluated for nucleic acid integrity by quantification of the human cytochrome C oxidase subunit 1 (cox 1) gene and bacterial 16S rRNA gene (V3 region) with SYBR Green assays, as described previously. Gardnerella vaginalis was quantified with primers and SYBR Green assay conditions determined previously. The presence of Mollicutes (Mycoplasma or Ureaplasma) was determined by targeting the 16S rRNA gene by using a conventional, semi-nested PCR, and Ureaplasma spp. were detected by specific PCR for the multiple-banded antigen gene.
cpn60 Universal Target (UT) PCR and pyrosequencing
PCR was carried out by using a cocktail of cpn 60 UT-specific primers consisting of a 1:3 molar ratio of primers H279/H280:H1612/H1613, as described previously[29, 42, 43]. Primer sets were modified at the 5' end with one of 24 unique decamer multiplexing identification (MID) sequences, as per the manufacturer’s recommendations (Roche, Brandford, CT, USA). In addition, the Mollicutes-specific 16S rRNA gene PCR product from 12 samples was pooled and pyrosequenced. Amplicons were pooled in equimolar concentrations to create libraries for sequencing on the GS FLX Titanium platform. Emulsion PCR and sequencing were performed at the National Research Council, Saskatoon, SK, Canada.
Analysis of operational taxonomic units (OTUs)
Pyrosequencing data were processed by using the default on-rig procedures from 454/Roche. MID-partitioned sequences were processed with the microbial Profiling Using Metagenomic Assembly (mPUMA) pipeline (http://mpuma.sourceforge.net) with default settings to generate operational taxonomic units (OTUs) with gsAssembler (Roche). OTUs were screened and filtered for chimeras with Chaban Chimera Checker (C3), followed by manual curation. OTU abundance was calculated based on mapping of sequence reads to OTU sequences by using Bowtie 2 in mPUMA. OTUs were identified by watered-Blast comparison to the cpn 60 reference database, cpnDB_nr (downloaded on March 21, 2013, from http://www.cpndb.ca,), and OTUs having the same best database reference were grouped together into nearest neighbor “species,” whereas OTUs having less than 55% identity to any reference sequence were removed from the dataset as non-cpn 60 sequence. Raw sequence data files were deposited to the NCBI Short Read Archive (BioProject PRJNA210319). Samples with fewer than 100 sequence reads after processing were removed, leaving a total of 76 samples from 27 women for analysis.
We examined the demographic and clinical parameters of the study cohort by using descriptive statistics.
Alpha and beta diversity comparisons
Quantitative Insights Into Microbial Ecology (QIIME) pipeline was used to calculate Shannon diversity index, Chao1 estimated number of species, and jackknifed beta diversity from Bray-Curtis distance matrices. For these analyses, the data grouped to the nearest neighbor “species” level were used, and all diversity measures were bootstrapped 100 times at 1,000 reads per sample or their sample maximum when less than 1,000. Rarefaction plots of Chao1 were generated to ensure that an adequate sampling depth for each sample was achieved.
Alpha diversity measures of the samples were averaged across the 100 bootstrapped datasets at 1,000 reads per sample, and compared among the four menstrual phases (menstrual, follicular, periovulatory, and luteal) by using a linear nested mixed-effects model for Shannon diversity and a Poisson nested mixed-effects model for Chao1 nested within subject. Post hoc Tukey comparisons were conducted when needed by using the “multcomp” package in R[47, 48].
Beta diversity (ecologic distance) was calculated as the average pair-wise distances from the 100 bootstrapped datasets at 1,000 reads per sample. We compared ecologic distances of samples within subjects with those between by comparing all within-subject pair-wise distances with all between-subject pair-wise distances by using ANOVA. In addition, we compared all pair-wise distances within menstrual phases pooled (that is, all within menstrual, and all within follicular, all within periovulatory, and all within luteal) with all between-phase distances (again pooled) by using ANOVA. These analyses do not take into account multiple samples per subject, as they are based on pair-wise distances among samples; they should therefore be considered exploratory.
As the average of the bootstrapped distance matrices did not include the samples with <1,000 reads, we generated another Bray-Curtis distance matrix on the proportion of reads of each nearest neighbor “species” per sample by using only taxa that represented at least 1% of the reads of at least one sample by using the vegan package in R. Hierarchic clustering was achieved by using this distance matrix and ward linkage. Bootstrap support for clusters was assessed by using average Jaccard similarities from the clusterboot function in the “fpc” package in R.
The results of the clustering are nearly identical when the bootstrapped averages for all samples with >1,000 reads are used instead. The additional samples contribute to some clusters, but do not change the structure of the clustering results.
In addition to the linear mixed-effects models on alpha diversity, and the comparisons of within versus between-woman ecologic distances, we used a graphic method to assess the “stability” of microbial communities across samples from the same subject. To explore how the relative abundance of taxa changed over time, we plotted proportional area plots of taxonomic composition for each woman by using taxa that represented at least 10% of at least one sample. Women with samples that had <25% change in the proportions of dominant taxa were considered stable. Stability by this measure, and stability as a measure of whether a woman remained in one of the previously defined clusters, or switched among clusters, were tested for associations with BMI (underweight, normal, overweight, obese), previous pregnancy, marital status (partnered versus single), alcohol use, vaginal sex, use of condoms, and use of unscented tampons during the study period by using Fisher Exact tests.
Correlation of sequencing-read abundance to qPCR for Gardnerella vaginalis
The proportion of total pyrosequencing reads obtained per sample corresponding to G. vaginalis (nearest neighbor “species”) were compared with the absolute amount of G. vaginalis DNA detected by qPCR in each sample by using a Spearman rank-order correlation test (SPSS Inc., Chicago, IL, USA).
Demographic and clinical data for study participants ( N = 27)
34.96 ± 4.17 (18-53)
23.47 ± 1.86 (16.8-36)
Current alcohol use
Sexual partners in past year
1.04 ± 0.28 (0-3)
Sexual partners in past 2 months
0.78 ± 0.23 (0-2)
Vaginal intercourse during study period
Condom use during study period
Cycle duration (days)
28.05 ± 0.68 (26-31)
Tampon use during study period
BV and STI history
Diagnosed with BV, yeast, or STI in lifetime
Diagnosed with BV in lifetime
BV episode in past year
Antimicrobial use in past 2 months a
Sample integrity evaluation
To confirm the presence of amplifiable nucleic acid in the DNAzol samples, samples were assayed for human cytochrome C oxidase gene and total bacterial 16S rRNA gene content. Human DNA was detectable in 72 (98.6%) of 73 samples tested, whereas total 16S rRNA gene content fell within the ranges of less than 104 copies/swab (33 samples), 104 to 105 copies/swab (31 samples) or 105 to 106 copies/swab (12 samples). All samples generated cpn 60-UT amplicons for pyrosequencing.
Microbiome profile generation
After sample collection and processing, an average of three vaginal samples per woman (range, 1 to 4; median, 3) were available for analysis. Microbial profiles were determined from an average of 7,585 sequence reads per sample (range from 125 to 37,419; median, 4,256), generating 567 unique OTUs. When combined at the nearest neighbor “species” level (OTU with same best reference database match), 73 bacterial “species” were identified with an average minimum percentage identity of 94.0% ± 7.1%, an average maximum percentage identity of 96.6% ± 5.2%, and an average overall percentage identify of 95.5% ± 5.5% to the reference sequence. A detailed summary of the nearest neighbor “species” identified in each sample is provided in an additional file (see Additional file1).
To evaluate whether the number of sequence reads per sample was a sufficient representation of the vaginal microbiota, rarefaction plots of Chao1-estimated number of species were calculated. Graphic depictions of the rarefaction plots for subsampled data are presented in an additional file (see Additional file2). In most samples, the Chao1 values remained relatively consistent throughout (flat line across the plot), indicating that the sample richness did not change significantly when 100 to 1,000 sequence reads per sample were examined. This consistency indicated that most of the community was captured at the sampling depth achieved.
Vaginal microbiome richness, diversity, and ecologic distance
The Shannon diversity index provides a quantitative measure of species diversity (richness and evenness), whereas the Chao1-estimated number of species provides a quantitative measure of species richness. No significant difference was found among the menstrual phases (linear mixed-effect model; likelihood ratio test, P = 1) in the average Shannon diversity of vaginal microbiome samples after accounting for multiple samples from each individual. However, a significant relationship between menstrual phase and the Chao1 estimated number of species was observed (Poisson mixed-effects model: likelihood ratio test, P = 0.03), with post hoc comparisons suggesting that the follicular phase had, on average, 1.3 more species than the luteal phase (adjusted P = 0.01), after multiple samples per individual were taken into account. None of the other menstrual-phase comparisons showed significant differences. Figures showing the average Shannon diversity index and Chao1 estimated number of species calculated for each menstrual phase are presented in an additional file (see Additional file3).
Vaginal microbiome composition
A relation did not appear to be present between clustering and menstrual phase, as most clusters contained samples from all phases (Figure 2). An inventory of the dominant species profile identified from the different menstrual-phase samples provided by each individual is provided in a supplementary file (see Additional file4).To gain a clearer picture of the stability of vaginal microbiome community composition over time, proportional area plots of the taxonomic composition of each woman’s microbiome were assembled (Figure 3). The plots grouped into four general types of vaginal microbiome profiles. The first group consisted of 12 women with relatively stable bacterial profiles over the sampling course (<25% variability in relative proportion of dominant bacterial species; Figure 3A). These women were also among those who did not switch between community clusters, as defined earlier by hierarchic clustering.The second group comprised nine women with microbiota profiles that consisted of the same bacterial organisms, but the relative proportions of each “species” fluctuated by >25% over time (Figure 3B). Five of these subjects also switched between clusters, as defined earlier, whereas the remaining four did not.The third group consisted of three women who showed dramatic changes in their microbiota over the sampling course (a complete change in the dominant members of the community; Figure 3C), whereas the fourth group contained two women with >25% variability in relative proportions of bacterial “species,” as well as an apparent introduction of new bacterial “species” (Figure 3D).
Differences in microbiome composition were reflected in analogous differences in the Shannon diversity index and Chao1-estimate values. Plots depicting temporal variation in the bootstrapped Shannon diversity index and Chao1 estimates for each individual are included in an additional file (see Additional file5).
No evidence was seen that “stability,” as defined either as <25% change in dominant taxa, or as cluster-switching, was related to BMI category, previous pregnancy, marital status, alcohol use, vaginal sex, use of condoms, or use of unscented tampons during the study period by Fisher Exact tests (all P > 0.05). Other clinical covariates were not investigated, as too few subjects existed per category (Table 1).
Evaluation of nearest-neighbor “species” abundance by G. vaginalis levels
In an effort to confirm independently whether the nearest-neighbor “species” sequence read abundance reflected actual bacterial species levels in the samples, G. vaginalis was quantified in the samples, and the quantitative PCR results compared with sequence read counts. G. vaginalis was chosen as a target because it was detected in 74% of samples at a wide range of sequence read abundances. The proportion of total reads obtained per sample corresponding to G. vaginalis was significantly correlated with the amount of G. vaginalis detected by qPCR in each sample (rs = 0.405; n = 76; P < 0.001).
Mycoplasma and Ureaplasma in the vaginal microbiome
A known limitation of cpn 60-based microbial profiling is that some species of Mollicutes do not contain this gene. To address this issue, all samples were screened for Mollicutes by using a conventional PCR-based test designed to detect Mycoplasma and/or Ureaplasma species. Mollicutes were detected in 59 (80.8%) of 73 samples, representing 23 (85.2%) of 27 women (Figure 2). When a Ureaplasma-specific PCR was applied to the same samples, 23 (31.5%) of 73 returned positive results from 10 (37.0%) of 27 women, comprising of 18 samples (24.7%) with U. parvum (eight women) and five samples (6.8%) with U. urealyticum (two women). To understand the composition of the Mollicutes-level PCR products further, PCR products from 12 samples (representing five women in either follicular, periovulatory, or luteal phases) were pooled and pyrosequenced to generate 54,926 reads. Mollicutes-like sequences accounted for 77.7% of reads, with clear identities (>95%) to Mycoplasma hominis and Ureaplasma spp. observed for 76.2% of the sequences (Ureaplasma cannot be speciated with this region of the 16S rRNA gene). The remaining 1.5% of Mollicutes-like reads grouped into seven distinct OTU that had only 85% to 90% identity to known species, suggesting that novel Mollicutes may be present in the microbiome.
Sequences matching Mycoplasma genitalium were not detected. Interestingly, this PCR also generated amplicons from several Staphylococcus species, including S. epidermidis (17.3% of reads), S. pasteuri (0.21% of reads), S. hyicus (0.04% of reads), and Staphylococcus spp. (2.9% of reads), as well as the Lactobacillus species L. iners (1.6% of reads), L. crispatus (0.04% of reads), and L. jensenii (0.07% of reads). This indicated that some of the PCR product generated was not Mollicutes-derived and suggests that this PCR assay has the potential to generate a false-positive result.
The results of our cpn 60-based study of the vaginal microbiota through a menstrual cycle demonstrate that in this cohort, only a few community states exist in healthy women, with L. crispatus, L. iners, and L. jensenii being the dominant members. This is consistent with previously reported studies[12, 13, 25]. Of note, Lactobacillus gasseri is prominent in other studies[12, 13], but it was not a dominant organism in any woman’s microbiota profile in this study.
In overall richness, diversity, and ecologic distance, we observed an element of a personalized compositional pattern in most women over time: profiles generated from samples collected at different menstrual phases from the same woman showed a higher degree of shared similarity in community structure versus profiles generated from different women but collected in the same menstrual phase[13, 51]. However, given the limited sample size, the extent of a “personalized” microbiota remains to be confirmed.
We did observe a statistically significant relation between menstrual phase and the Chao1 estimated number of species, in which follicular-phase samples had on average 1.3 more species than did luteal-phase samples. However, we do not believe that this observation is biologically significant, given that presence or absence of a single species can be affected by technical nuances and random sampling effects.
Genital Mycoplasma and Ureaplasma have been associated with female genital infections such as vaginitis, cervicitis, and pelvic inflammatory disease, and with a variety of other negative reproductive and neonatal outcomes[53–55]. M. hominis, M. genitalium, U. parvum, and U. urealyticum are regularly detected in vaginal samples by using culture or taxon-specific PCR methods, although the reported prevalence is wide-ranging (20% to 80% depending on the study population). These species are much more rarely reported in 16S rRNA gene-based microbiota analysis[12, 14, 51], likely because of universal primer bias. The lack of a cpn 60 gene in Mollicutes is not a universal phenomenon, because several species of Mycoplasma, including M. genitalium do have a cpn 60 gene. However, M. genitalium has not been detected in any cpn 60-based studies of the human vaginal microbiome to date[28–30]. Mycoplasma and Ureaplasma are commonly reported vaginal constituents[57, 58], further supporting the suggestion that these organisms are part of the vaginal microbiota of many clinically healthy women. As well, deep sequencing revealed seven novel Mollicutes-like sequences that could indicate uncharacterized species. Interestingly, deep pyrosequencing of the reportedly Mollicutes-specific PCR products from a pool of 12 PCR-positive samples showed that primer specificity was imperfect. The generation of amplicon from other organisms, particularly Staphylococcus and Lactobacillus species, that was the same size as the targeted Mollicutes product might result in some overreporting of Mollicutes detection in studies using these primers.
Our study diverges from several others in the number of vaginal profiles dominated by non-Gardnerella Bifidobacteriales. Bifidobacterium species (high G + C, Gram positive, Actinobacteria) are known members of the vaginal microbiota, and genera like Bifidobacterium and Alloscardovia have been isolated from vaginal samples in several culture-based investigations[19, 21, 59–62]. Culture-independent studies based on the 16S rRNA gene have also reported Bifidobacterium species in the vaginal microbiome[14, 63, 64], with one study finding it the dominant community member for two of 20 women in the study. More commonly, 16S rRNA gene-based investigations tend not to report Bifidobacterium in the vaginal microbiome[12–14, 51]. This phenomenon echoes previous findings in studies of the intestinal microbiome, in which culture-based investigations revealed a wealth of Bifidobacterium, whereas 16S rRNA gene-based culture-independent studies detected very few, leading to the discovery that many universal 16S rRNA gene PCR primers were a poor match for this genus[26, 66, 67].
Further hindering the identification of Bifidobacterium is its relatedness to G. vaginalis, a species that is prevalent in the vaginal microbiota, and belongs to the same taxonomic family. The cpn 60-based universal PCR protocol used in this study has been shown experimentally to represent bifidobacteria more accurately in intestinal microbiomes compared with 16S rRNA gene-based universal PCR, and the cpn 60 target sequence is clearly distinguishable between Gardnerella and Bifidobacterium (average sequence identity of only 75%, compared with 16S rRNA gene identities of >90%). The finding that five women in our study group had Bifidobacterium breve, Bifidobacterium longum, or Alloscardovia omnicolens as a dominant vaginal organism seems reasonable. Bifidobacteria are generally considered to be beneficial members of the intestinal microbiota, although their role in the vaginal microbiota has not yet been elucidated. It is conceivable that Bifidobacterium, lactic acid-producing bacteria, could have a protective or health-promoting effect in the vagina analogous to Lactobacillus. Bifidobacteria also appear to play an important role in early infant health and development, and their presence in the vaginal microbiota of healthy, reproductive-aged women could provide a means of transfer from the mother to the newborn during birth.
Measures were taken to evaluate sample quality, representativeness of the sequence read abundance, and known limitations of the technique used. Sample quality was assessed from the ability to amplify both human (via the cox 1 gene) and bacterial DNA (via the 16S rRNA and cpn 60 genes) from the samples. The fact that bacterial amplicon could be generated from all samples and human amplicon from all but one sample indicated that nucleic acid quality was maintained at a sufficient level for analysis. It was then determined empirically that cpn 60 amplicon generation and sequencing faithfully represented actual starting amounts of target in the example case of G. vaginalis. In addition, it was anticipated that Mollicutes species would be poorly characterized by cpn 60 analysis, so targeted investigation of this important group was separately undertaken. All of these measures added additional experimental data to the results generated from the cpn 60 analysis and helped to ensure that the descriptions of the vaginal microbiota under investigation were as faithful and representative as possible.
The results of this study suggest that the specific menstrual phase of a woman is not predictive of her vaginal microbiota at that time. For 18 (69%) women, all the samples collected throughout the menstrual cycle clustered together by composition. An additional six (23%) women transitioned between Lactobacillus- dominated clusters throughout the menstrual cycle, whereas only two (8%) women had sample compositions transition from the heterogeneous cluster to a Lactobacillus-dominated cluster. No apparent pattern of vaginal microbiome composition change was noted by menstrual phase, and no association with personal health practices or sexual activity. It is worth noting that the majority of the women in this study had microbiota dominated by Lactobacillus species. Although this is a commonly observed trend, we recognize that the observations made from this cohort may not apply to populations of women with a higher prevalence of non-Lactobacillus-dominated vaginal microbiota. Further studies involving more Canadian women over longer periods are needed to continue addressing important questions about the composition and stability of the vaginal microbiome. It is also important to ensure that these studies are done in different cohorts with variable ethnic and behavioral characteristics in multiple sites around the world.
The composition of the vaginal microbiome did not appear to be directly linked to the menstrual phases of women in our study. Some woman had vaginal microbial communities that were apparently stable in composition and abundant throughout the menstrual cycle, whereas others underwent moderate or dramatic shifts over time. Despite these shifts, the overall community composition tended toward one of only a few compositional structure types dominated by Lactobacillus species at any single time. The use of the protein-coding cpn 60 gene as a community barcode marker, unlike some other methods, revealed Bifidobacteriales-dominated vaginal profiles whose role in health warrants further investigation. The degree to which specific aberrations of a woman’s vaginal microbiome occur and lead to symptomatic disease requires urgent attention, given the prevalence of such illnesses and its adverse effects on quality of life and reproductive outcomes.
- 16S rRNA:
16-Svedberg ribosomal RNA
- cox 1 cytochrome C:
oxidase subunit 1 gene
- cpn 60:
- G. vaginalis :
- G + C:
guanine + cytosine
- L. crispatus :
- L. iners :
- L. jensenii :
- M. genitalium :
- M. hominis :
Microbial Profiling Using Metagenomic Assembly
National Center for Biotechnology Information
operational taxonomic unit
polymerase chain reaction
Quantitative Insights Into Microbial Ecology
quantitative polymerase chain reaction
- S. epidermidis :
- S. hyicus :
- S. pasteuri :
- U. parvum :
- U. urealyticum :
The authors thank the women who took the time and effort to participate in this study. Funding support was provided from a Canadian Institutes of Health Research (CIHR) Catalyst Award to DM.
- Ma B, Forney LJ, Ravel J: Vaginal microbiome: rethinking health and disease. Annu Rev Microbiol. 2012, 66: 371-389.View ArticlePubMedPubMed CentralGoogle Scholar
- Cadieux PA, Burton J, Devillard E, Reid G: Lactobacillus by-products inhibit the growth and virulence of uropathogenic Escherichia coli. J Physiol Pharmacol. 2009, 60: 13-18.PubMedGoogle Scholar
- Matu MN, Orinda GO, Njagi EN, Cohen CR, Bukusi EA: In vitro inhibitory activity of human vaginal lactobacilli against pathogenic bacteria associated with bacterial vaginosis in Kenyan women. Anaerobe. 2010, 16: 210-215.View ArticlePubMedGoogle Scholar
- Strus M, Brzychczy-Wloch M, Gosiewski T, Kochan P, Heczko PB: The in vitro effect of hydrogen peroxide on vaginal microbial communities. FEMS Immunol Med Microbiol. 2006, 48: 56-63.View ArticlePubMedGoogle Scholar
- O'Hanlon DE, Moench TR, Cone RA: In vaginal fluid, bacteria associated with bacterial vaginosis can be suppressed with lactic acid but not hydrogen peroxide. BMC Infect Dis. 2011, 11: 200-View ArticlePubMedPubMed CentralGoogle Scholar
- Reid G, McGroarty JA, Tomeczek L, Bruce AW: Identification and plasmid profiles of Lactobacillus species from the vagina of 100 healthy women. FEMS Immunol Med Microbiol. 1996, 15: 23-26.View ArticlePubMedGoogle Scholar
- Velraeds MMC, vander Mei HC, Reid G, Busscher HJ: Physicochemical and biochemical characterization of biosurfactants released by Lactobacillus strains. Colloids Surf B Biointerfaces. 1996, 8: 51-61.View ArticleGoogle Scholar
- Vitali B, Pugliese C, Biagi E, Candela M, Turroni S, Bellen G, Donders GG, Brigidi P: Dynamics of vaginal bacterial communities in women developing bacterial vaginosis, candidiasis, or no infection, analyzed by PCR-denaturing gradient gel electrophoresis and real-time PCR. Appl Environ Microbiol. 2007, 73: 5731-5741.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou X, Brown CJ, Abdo Z, Davis CC, Hansmann MA, Joyce P, Foster JA, Forney LJ: Differences in the composition of vaginal microbial communities found in healthy Caucasian and black women. ISME J. 2007, 1: 121-133.View ArticlePubMedGoogle Scholar
- Yamamoto T, Zhou X, Williams CJ, Hochwalt A, Forney LJ: Bacterial populations in the vaginas of healthy adolescent women. J Pediatr Adolesc Gynecol. 2009, 22: 11-18.View ArticlePubMedGoogle Scholar
- Zhou X, Hansmann MA, Davis CC, Suzuki H, Brown CJ, Schutte U, Pierson JD, Forney LJ: The vaginal bacterial communities of Japanese women resemble those of women in other racial groups. FEMS Immunol Med Microbiol. 2010, 58: 169-181.View ArticlePubMedPubMed CentralGoogle Scholar
- Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SS, McCulle SL, Karlebach S, Gorle R, Russell J, Tacket CO, Brotman RM, Davis CC, Ault K, Peralta L, Forney LJ: Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci U S A. 2011, 108: 4680-4687.View ArticlePubMedPubMed CentralGoogle Scholar
- Gajer P, Brotman RM, Bai G, Sakamoto J, Schütte UM, Zhong X, Koenig SS, Fu L, Ma ZS, Zhou X, Abdo Z, Forney LJ, Ravel J: Temporal dynamics of the human vaginal microbiota. Sci Transl Med. 2012, 4: 132-152.View ArticleGoogle Scholar
- Drell T, Lillsaar T, Tummeleht L, Simm J, Aaspõllu A, Väin E, Saarma I, Salumets A, Donders GG, Metsis M: Characterization of the vaginal micro- and mycobiome in asymptomatic reproductive-age Estonian women. PLoS One. 2013, 8: e54379-View ArticlePubMedPubMed CentralGoogle Scholar
- Chan RC, Bruce AW, Reid G: Adherence of cervical, vaginal and distal urethral normal microbial flora to human uroepithelial cells and the inhibition of adherence of gram-negative uropathogens by competitive exclusion. J Urol. 1984, 131: 596-601.PubMedGoogle Scholar
- Semmens JP, Wagner G: Estrogen deprivation and vaginal function in postmenopausal women. JAMA. 1982, 248: 445-448.View ArticlePubMedGoogle Scholar
- Wagner G, Ottesen B: Vaginal physiology during menstruation. Ann Intern Med. 1982, 96: 921-923.View ArticlePubMedGoogle Scholar
- Brenner RM, West NB: Hormonal regulation of the reproductive tract in female mammals. Annu Rev Physiol. 1975, 37: 273-302.View ArticlePubMedGoogle Scholar
- Bartlett JG, Onderdonk AB, Drude E, Goldstein C, Anderka M, Alpert S, McCormack WM: Quantitative bacteriology of the vaginal flora. J Infect Dis. 1977, 136: 271-277.View ArticlePubMedGoogle Scholar
- Sautter RL, Brown WJ: Sequential vaginal cultures from normal young women. J Clin Microbiol. 1980, 11: 479-484.PubMedPubMed CentralGoogle Scholar
- Wilks M, Tabaqchali S: Quantitative bacteriology of the vaginal flora during the menstrual cycle. J Med Microbiol. 1987, 24: 241-245.View ArticlePubMedGoogle Scholar
- Priestley CJ, Jones BM, Dhar J, Goodwin L: What is normal vaginal flora?. Genitourin Med. 1997, 73: 23-28.PubMedPubMed CentralGoogle Scholar
- Eschenbach DA, Thwin SS, Patton DL, Hooton TM, Stapleton AE, Agnew K, Winter C, Meier A, Stamm WE: Influence of the normal menstrual cycle on vaginal tissue, discharge, and microflora. Clin Infect Dis. 2000, 30: 901-907.View ArticlePubMedGoogle Scholar
- Srinivasan S, Liu C, Mitchell CM, Fiedler TL, Thomas KK, Agnew KJ, Marrazzo JM, Fredricks DN: Temporal variability of human vaginal bacteria and relationship with bacterial vaginosis. PLoS One. 2010, 5: e10197-View ArticlePubMedPubMed CentralGoogle Scholar
- Hickey RJ, Abdo Z, Zhou X, Nemeth K, Hansmann M, Osborn TW, Wang F, Forney LJ: Effects of tampons and menses on the composition and diversity of vaginal microbial communities over time. BJOG. 2013, 120: 695-704.View ArticlePubMedGoogle Scholar
- Hill JE, Fernando WM, Zello GA, Tyler RT, Dahl WJ, Van Kessel AG: Improvement of the representation of bifidobacteria in fecal microbiota metagenomic libraries by application of the cpn60 universal primer cocktail. Appl Environ Microbiol. 2010, 76: 4550-4552.View ArticlePubMedPubMed CentralGoogle Scholar
- Links MG, Dumonceaux TJ, Hemmingsen SM, Hill JE: The chaperonin-60 universal target is a barcode for bacteria that enables de novo assembly of metagenomic sequence data. PLoS One. 2012, 7: e49755-View ArticlePubMedPubMed CentralGoogle Scholar
- Hill JE, Goh SH, Money DM, Doyle M, Li A, Crosby WL, Links M, Leung A, Chan D, Hemmingsen SM: Characterization of vaginal microflora of healthy, nonpregnant women by chaperonin-60 sequence-based methods. Am J Obstet Gynecol. 2005, 193: 682-692.View ArticlePubMedGoogle Scholar
- Schellenberg J, Links MG, Hill JE, Dumonceaux TJ, Peters GA, Tyler S, Ball TB, Severini A, Plummer FA: Pyrosequencing of the chaperonin-60 universal target as a tool for determining microbial community composition. Appl Environ Microbiol. 2009, 75: 2889-2898.View ArticlePubMedPubMed CentralGoogle Scholar
- Schellenberg JJ, Links MG, Hill JE, Dumonceaux TJ, Kimani J, Jaoko W, Wachihi C, Mungai JN, Peters GA, Tyler S, Graham M, Severini A, Fowke KR, Ball TB, Plummer FA: Molecular definition of vaginal microbiota in East African commercial sex workers. Appl Environ Microbiol. 2011, 77: 4066-4074.View ArticlePubMedPubMed CentralGoogle Scholar
- Paramel Jayaprakash T, Schellenberg JJ, Hill JE: Resolution and characterization of distinct cpn60-based subgroups of Gardnerella vaginalis in the vaginal microbiota. PLoS One. 2012, 7: e43009-View ArticlePubMedPubMed CentralGoogle Scholar
- Keane FE, Ison CA, Taylor-Robinson D: A longitudinal study of the vaginal flora over a menstrual cycle. Int J STD AIDS. 1997, 8: 489-494.View ArticlePubMedGoogle Scholar
- Thoma ME, Gray RH, Kiwanuka N, Aluma S, Wang MC, Sewankambo N, Wawer MJ: Longitudinal changes in vaginal microbiota composition assessed by Gram stain among never sexually active pre- and postmenarcheal adolescents in Rakai, Uganda. J Pediatr Adolesc Gynecol. 2011, 24: 42-47.View ArticlePubMedPubMed CentralGoogle Scholar
- Brotman RM, Ravel J, Cone RA, Zenilman JM: Rapid fluctuation of the vaginal microbiota measured by Gram stain analysis. Sex Transm Infect. 2010, 86: 297-302.View ArticlePubMedPubMed CentralGoogle Scholar
- Petignat P, Hankins C, Walmsley S, Money D, Provencher D, Pourreaux K, Kornegay J, Rouah F, Coutlee F: Self-sampling is associated with increased detection of human papillomavirus DNA in the genital tract of HIV-seropositive women. Clin Infect Dis. 2005, 41: 527-534.View ArticlePubMedGoogle Scholar
- Ogilvie GS, Patrick DM, Schulzer M, Sellors JW, Petric M, Chambers K, White R, FitzGerald JM: Diagnostic accuracy of self-collected vaginal specimens for human papillomavirus compared to clinician collected human papillomavirus specimens: a meta-analysis. Sex Transm Infect. 2005, 81: 207-212.View ArticlePubMedPubMed CentralGoogle Scholar
- Chaban B, Albert A, Links MG, Gardy J, Tang P, Hill JE: Characterization of the upper respiratory tract microbiomes of patients with pandemic H1N1 influenza. PLoS One. 2013, 8: e69559-View ArticlePubMedPubMed CentralGoogle Scholar
- Menard JP, Fenollar F, Henry M, Bretelle F, Raoult D: Molecular quantification of Gardnerella vaginalis and Atopobium vaginae loads to predict bacterial vaginosis. Clin Infect Dis. 2008, 47: 33-43.View ArticlePubMedGoogle Scholar
- Dumonceaux TJ, Schellenberg J, Goleski V, Hill JE, Jaoko W, Kimani J, Money D, Ball TB, Plummer FA, Severini A: Multiplex detection of bacteria associated with normal microbiota and with bacterial vaginosis in vaginal swabs by use of oligonucleotide-coupled fluorescent microspheres. J Clin Microbiol. 2009, 47: 4067-4077.View ArticlePubMedPubMed CentralGoogle Scholar
- van Kuppeveld FJ, van der Logt JT, Angulo AF, van Zoest MJ, Quint WG, Niesters HG, Galama JM, Melchers WJ: Genus- and species-specific identification of mycoplasmas by 16S rRNA amplification. Appl Environ Microbiol. 1993, 59: 655-PubMedPubMed CentralGoogle Scholar
- Watson HL, Blalock DK, Cassell GH: Variable antigens of Ureaplasma urealyticum containing both serovar-specific and serovar-cross-reactive epitopes. Infect Immun. 1990, 58: 3679-3688.PubMedPubMed CentralGoogle Scholar
- Hill JE, Town JR, Hemmingsen SM: Improved template representation in cpn60 polymerase chain reaction (PCR) product libraries generated from complex templates by application of a specific mixture of PCR primers. Environ Microbiol. 2006, 8: 741-746.View ArticlePubMedGoogle Scholar
- Schellenberg J, Links MG, Hill JE, Hemmingsen SM, Peters GA, Dumonceaux TJ: Pyrosequencing of chaperonin-60 (cpn60) amplicons as a means of determining microbial community composition. Methods Mol Biol. 2011, 733: 143-158.View ArticlePubMedGoogle Scholar
- Links MG, Chaban B, Hemmingsen SM, Muirhead K, Hill JE: mPUMA: a computational approach to microbiota analysis by de novo assembly of operational taxonomic units based on protein-coding barcode sequences. Microbiome. 2013, 1: 23-View ArticlePubMedPubMed CentralGoogle Scholar
- Hill JE, Penny SL, Crowell KG, Goh SH, Hemmingsen SM: cpnDB: a chaperonin sequence database. Genome Res. 2004, 14: 1669-1675.View ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7: 335-336.View ArticlePubMedPubMed CentralGoogle Scholar
- Hothorn T, Bretz F, Westfall P: Simultaneous inference in general parametric models. Biomed J. 2008, 50: 346-363.Google Scholar
- Team RDC: R: A Language and Environment for Statistical Computing. 2012, Vienna, Austria: R Foundation for Statistical Computing, URL http://www.R-project.org/. ISBN 3-900051-07-0Google Scholar
- Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O'Hara RB, Simpson GL, Solymos P, Henry M, Stevens H, Wagner H: Vegan: Community Ecology Package. R package version 2.0-3. 2012,http:CRAN.R-project.org/package=vegan,Google Scholar
- Hennig C: fpc: Flexible Procedures For Clustering. R Package Version 2.1-6. 2013,http://CRAN.R-project.org/package=fpc,Google Scholar
- Hickey RJ, Zhou X, Pierson JD, Ravel J, Forney LJ: Understanding vaginal microbiome complexity from an ecological perspective. Transl Res. 2012, 160: 267-282.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou J, Jiang YH, Deng Y, Shi Z, Zhou BY, Xue K, Wu L, He Z, Yang Y: Random sampling process leads to overestimation of beta-diversity of microbial communities. MBio. 2013, 4: e00324-e00313.PubMedPubMed CentralGoogle Scholar
- Schlicht MJ, Lovrich SD, Sartin JS, Karpinsky P, Callister SM, Agger WA: High prevalence of genital mycoplasmas among sexually active young adults with urethritis or cervicitis symptoms in La Crosse, Wisconsin. J Clin Microbiol. 2004, 42: 4636-4640.View ArticlePubMedPubMed CentralGoogle Scholar
- Taylor-Robinson D: Infections due to species of Mycoplasma and Ureaplasma: an update. Clin Infect Dis. 1996, 23: 671-682.View ArticlePubMedGoogle Scholar
- Waites KB, Katz B, Schelonka RL: Mycoplasmas and ureaplasmas as neonatal pathogens. Clin Microbiol Rev. 2005, 18: 757-789.View ArticlePubMedPubMed CentralGoogle Scholar
- Hummelen R, Fernandes AD, Macklaim JM, Dickson RJ, Changalucha J, Gloor GB, Reid G: Deep sequencing of the vaginal microbiota of women with HIV. PLoS One. 2010, 5: e12078-View ArticlePubMedPubMed CentralGoogle Scholar
- Tibaldi C, Cappello N, Latino MA, Masuelli G, Marini S, Benedetto C: Vaginal and endocervical microorganisms in symptomatic and asymptomatic non-pregnant females: risk factors and rates of occurrence. Clin Microbiol Infect. 2009, 15: 670-679.View ArticlePubMedGoogle Scholar
- McKechnie ML, Hillman RJ, Jones R, Lowe PC, Couldwell DL, Davies SC, King F, Kong F, Gilbert GL: The prevalence of urogenital micro-organisms detected by a multiplex PCR-reverse line blot assay in women attending three sexual health clinics in Sydney, Australia. J Med Microbiol. 2011, 60: 1010-1016.View ArticlePubMedGoogle Scholar
- Hillier SL, Krohn MA, Rabe LK, Klebanoff SJ, Eschenbach DA: The normal vaginal flora, H2O2-producing lactobacilli, and bacterial vaginosis in pregnant women. Clin Infect Dis. 1993, 16: S273-S281.View ArticlePubMedGoogle Scholar
- Verhelst R, Verstraelen H, Claeys G, Verschraegen G, Delanghe J, Van Simaey L, De Ganck C, Temmerman M, Vaneechoutte M: Cloning of 16S rRNA genes amplified from normal and disturbed vaginal microflora suggests a strong association between Atopobium vaginae, Gardnerella vaginalis and bacterial vaginosis. BMC Microbiol. 2004, 4: 16-View ArticlePubMedPubMed CentralGoogle Scholar
- Santiago GL, Cools P, Verstraelen H, Trog M, Missine G, El Aila N, Verhelst R, Tency I, Claeys G, Temmerman M, Vaneechoutte M: Longitudinal study of the dynamics of vaginal microflora during two consecutive menstrual cycles. PLoS One. 2011, 6: e28180-View ArticlePubMedGoogle Scholar
- Korshunov VM, Gudieva ZA, Efimov BA, Pikina AP, Smeianov VV, Reid G, Korshunova OV, Tiutiunnik VL, Stepin II: [The vaginal Bifidobacterium flora in women of reproductive age]. Zh Mikrobiol Epidemiol Immunobiol. 1999, 4: 74-78.PubMedGoogle Scholar
- Sundquist A, Bigdeli S, Jalili R, Druzin ML, Waller S, Pullen KM, El-Sayed YY, Taslimi MM, Batzoglou S, Ronaghi M: Bacterial flora-typing with targeted, chip-based pyrosequencing. BMC Microbiol. 2007, 7: 108-View ArticlePubMedPubMed CentralGoogle Scholar
- Burton JP, Dixon JL, Reid G: Detection of Bifidobacterium species and Gardnerella vaginalis in the vagina using PCR and denaturing gradient gel electrophoresis (DGGE). Int J Gynaecol Obstet. 2003, 81: 61-63.View ArticlePubMedGoogle Scholar
- Hyman RW, Fukushima M, Diamond L, Kumm J, Giudice LC, Davis RW: Microbes on the human vaginal epithelium. Proc Natl Acad Sci U S A. 2005, 102: 7952-7957.View ArticlePubMedPubMed CentralGoogle Scholar
- Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO: Development of the human infant intestinal microbiota. PLoS Biol. 2007, 5: e177-View ArticlePubMedPubMed CentralGoogle Scholar
- Sim K, Cox MJ, Wopereis H, Martin R, Knol J, Li MS, Cookson WO, Moffatt MF, Kroll JS: Improved detection of bifidobacteria with optimised 16S rRNA-gene based pyrosequencing. PLoS One. 2012, 7: e32543-View ArticlePubMedPubMed CentralGoogle Scholar
- Oakley BB, Fiedler TL, Marrazzo JM, Fredricks DN: Diversity of human vaginal bacterial communities and associations with clinically defined bacterial vaginosis. Appl Environ Microbiol. 2008, 74: 4898-4909.View ArticlePubMedPubMed CentralGoogle Scholar
- Perdigón G, Locascio M, Medici M, Pesce de Ruiz Holgado A, Oliver G: Interaction of bifidobacteria with the gut and their influence in the immune function. Biocell. 2003, 27: 1-9.PubMedGoogle Scholar
- Azad MB, Konya T, Maughan H, Guttman DS, Field CJ, Chari RS, Sears MR, Becker AB, Scott JA, Kozyrskyj AL: Gut microbiota of healthy Canadian infants: profiles by mode of delivery and infant diet at 4 months. CMAJ. 2013, 185: 385-394.View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.