Skip to main content
  • Research article
  • Open access
  • Published:

Correlates of the molecular vaginal microbiota composition of African women



Sociodemographic, behavioral and clinical correlates of the vaginal microbiome (VMB) as characterized by molecular methods have not been adequately studied. VMB dominated by bacteria other than lactobacilli may cause inflammation, which may facilitate HIV acquisition and other adverse reproductive health outcomes.


We characterized the VMB of women in Kenya, Rwanda, South Africa and Tanzania (KRST) using a 16S rDNA phylogenetic microarray. Cytokines were quantified in cervicovaginal lavages. Potential sociodemographic, behavioral, and clinical correlates were also evaluated.


Three hundred thirteen samples from 230 women were available for analysis. Five VMB clusters were identified: one cluster each dominated by Lactobacillus crispatus (KRST-I) and L. iners (KRST-II), and three clusters not dominated by a single species but containing multiple (facultative) anaerobes (KRST-III/IV/V). Women in clusters KRST-I and II had lower mean concentrations of interleukin (IL)-1α (p < 0.001) and Granulocyte Colony Stimulating Factor (G-CSF) (p = 0.01), but higher concentrations of interferon-γ-induced protein (IP-10) (p < 0.01) than women in clusters KRST-III/IV/V. A lower proportion of women in cluster KRST-I tested positive for bacterial sexually transmitted infections (STIs; ptrend = 0.07) and urinary tract infection (UTI; p = 0.06), and a higher proportion of women in clusters KRST-I and II had vaginal candidiasis (ptrend = 0.09), but these associations did not reach statistical significance. Women who reported unusual vaginal discharge were more likely to belong to clusters KRST-III/IV/V (p = 0.05).


Vaginal dysbiosis in African women was significantly associated with vaginal inflammation; the associations with increased prevalence of STIs and UTI, and decreased prevalence of vaginal candidiasis, should be confirmed in larger studies.

Peer Review reports


Lactobacilli-dominated vaginal microbiota (VMB) have traditionally been considered to promote reproductive health of women and their fetuses by maintaining a low vaginal pH (<4.5), which restricts the growth of other bacteria and yeasts [1]. The clinical conditions caused by an imbalanced VMB include bacterial vaginosis (BV) and bacterial vaginitis [2]. In addition, VMB-associated bacterial communities may be influenced by other micro-organisms in the vagina, such as Candida species, Trichomonas vaginalis, and other sexually transmitted pathogens [3-6].

BV has traditionally been characterized as a reduction of vaginal lactobacilli and an overgrowth of other (facultative) anaerobic bacteria. In clinical settings, BV is typically diagnosed using Amsel criteria (three of the following 4 criteria should be present: 1) clue cells on wet mount microscopy; 2) a ‘fishy’ odour after adding 10% KOH to vaginal secretions; 3) vaginal pH > 4.5; and 4) thin, homogenous vaginal discharge [2]. In 1991, Nugent and colleagues developed a method that could be repeated by scoring a Gram stained slide based on microscopic visualization of three bacterial morphotypes (a Nugent score of 0–3 is considered normal vaginal microbiota, 4–6 intermediate microbiota, and 7–10 BV-positive) [7]. Nugent scoring is considered the gold standard for BV diagnosis and is typically used in research settings.

Since 2002, an increasing number of VMB studies have used molecular characterization of vaginal bacterial communities, such as next generation sequencing, quantitative PCR or microarray analysis of bacterial 16S rRNA genes. A recent review of 63 molecular VMB studies conducted between 2008 and 2013 concluded that lactobacilli-dominated VMB are indeed associated with a healthy vaginal micro-environment (but that Lactobacillus crispatus is more beneficial than L. iners) and that BV is best described as a polybacterial dysbiosis [6]. In most studies, the extent of dysbiosis correlated well with Nugent score and vaginal pH but not with the other Amsel criteria [3,6]. Some studies reported systematic VMB differences across ethnic groups, with Black and Hispanic American women being less likely to have a VMB dominated by L. crispatus, and having a higher average vaginal pH, than White and Asian American women [8-10]. However, data on VMB associations with genital immune responses and other potential sociodemographic, behavioral and clinical correlates are scarce and inconsistent [6].

The vaginal micro-environment is also important in the context of vaginal product development, such as vaginal microbicides for HIV prevention. Candidate products should not disturb, and would ideally promote, a lactobacilli-dominated VMB, and should not induce vaginal inflammation [11,12]. In vaginal microbicide safety trials, the VMB has traditionally been assessed by Nugent scoring, and vaginal inflammation by visual inspection during pelvic examination with or without colposcopy and quantification of a select number of cytokines in cervicovaginal lavages (CVLs) [12-15]. However, these data have always been difficult to interpret due to insufficient understanding of the normal background variation in women of different ages, behaviors and clinical conditions.

Here, we describe the bacterial composition of the VMB of different groups of women in four African countries (Kenya, Rwanda, South Africa and Tanzania) using traditional VMB characterization methods (Amsel criteria and Nugent scoring) as well as molecular methods (16S rDNA phylogenetic microarray). Potential VMB correlates, including cytokine concentrations in CVLs, and sociodemographic, behavioral, and clinical characteristics were also evaluated.


Study design and ethical approvals

Samples for phylogenetic microarray testing as well as data on potential VMB correlates were used from two studies [16,17]. The study contributing most samples and data was a multi-country prospective observational cohort study aimed at characterizing novel safety biomarkers for vaginal HIV microbicide development in East and South Africa (referred to as the Vaginal Biomarkers Study). The study was conducted in 2011–2012 with a cohort of 430 women from three African countries (Kenya, Rwanda and South Africa; Additional file 1). The participants included 109 HIV-negative adult women each in Kenya and South Africa, 30 HIV-negative adolescent women each in Kenya and South Africa, 30 HIV-negative pregnant women each in Kenya and South Africa, 31 HIV-negative women using traditional vaginal practice in South Africa, 30 HIV-negative adult women at high risk for HIV (mostly female sex workers) in Rwanda, and 30 HIV-positive women in Rwanda.

The second study was an intensive longitudinal study cohort of women conducted in 2009 for evaluating the impact of traditional vaginal practices on the vaginal micro-environment in Northwest Tanzania (referred to as the Tanzania Study; Additional file 1). Study participants were 100 women working in bars, guest-houses and other food and recreational facilities located in three towns adjacent to large gold or diamond mines [17,18].

Both studies were approved by all relevant institutional and national ethics committees (Additional file 1). All participants (or their guardians in the case of minors) provided written informed consent and received a modest reimbursement for each study visit (Additional file 1).

Study procedures

In the Vaginal Biomarkers Study, women were screened, and eligible consenting women were enrolled within four days of the last day of their menstrual period (visit 1). Most study groups described above included healthy, non-pregnant, HIV-negative women between 16 and 35 years of age. The pregnant women group included women who were at most 14 weeks into gestation as determined by abdominal ultrasound, and the HIV-positive women group consisted of women on antiretroviral treatment (ART) for at least six months, currently asymptomatic and with a CD4 count above 350 cells/μl. Women were excluded from all groups if they had a history of hysterectomy or other genital tract surgery in the three months before the screening visit; never had had penetrative vaginal intercourse; were enrolled in an HIV prevention study involving investigational products; had internal and/or external genital warts at screening and/or enrollment; or were breastfeeding and less than six months post-partum at the time of enrollment. Once enrolled, women returned for six follow-up visits (visits 2 to 5 at biweekly intervals over two menstrual cycles, and visits 6 and 7 at three and six months after visit 5), but this paper focuses on the screening and enrollment (visit 1) visits. The median time between these visits was 25 days (interquartile range (IQR) 14 – 39 days). At screening, women underwent face-to-face interviews, blood and urine sample collection with real-time HIV, pregnancy and urinary tract infection (UTI) testing, a pelvic examination with sample collection, and a general physical exam. Samples were subsequently tested for several sexually transmitted infections (STIs; see below), UTIs, BV by Amsel and Nugent criteria, and vaginal candidiasis. At enrollment, additional samples were collected for microarray testing and immunological assessments (see below).

In the Tanzania Study, participants were enrolled at any time during their menstrual cycle and followed every two to three days for four weeks (12 visits in total). For the microarray testing, two pairs of samples from 20 women (two to four weeks apart) were selected. At each visit women underwent a face-to-face interview, physical examination, pelvic examination, and sample collection for STI and BV by Nugent scoring. Samples for microarray analysis were collected either at enrollment (8 women) or visit 6 (12 women).

Participants in both studies received counseling and condoms free of charge. Women who tested positive for curable STIs, UTI, or symptomatic BV or vaginal candidiasis were treated by study clinicians using national treatment guidelines. All HIV-positive and pregnant women were linked to appropriate care in local public clinics.

Diagnostic and immunological testing

Cervicovaginal samples were collected in the following order: vaginal pH measurement, vaginal samples, CVLs, and endocervical specimens. All diagnostic tests were conducted at the study sites in Africa, and similar tests were used in the two studies, unless otherwise stated (see Additional file 1 for diagnostic details). HIV status was determined by locally approved rapid testing algorithm. Plasma samples were tested for herpes simplex type 2 (HSV-2) antibodies and for syphilis by Rapid Plasma Reagin test with confirmation by a Treponema pallidum-specific test. Endocervical swabs were tested for Neisseria gonorrhoeae and Chlamydia trachomatis by PCR. Vaginal swabs were used to prepare a wet mount (detection of trichomonads and clue cells, and after KOH addition yeasts and amine smell), a Gram stain for Nugent scoring (done centrally at the Institute of Tropical Medicine in Antwerp), and to inoculate Trichomonas vaginalis cultures. UTIs were diagnosed by the presence of white blood cells on a urine dipstick test and pregnancy by urine hCG test. Vaginal pH was measured using pH paper strips (pH range 3.6-6.1).

CVLs were obtained by irrigating the cervix and lateral vaginal walls with 10 ml of sterile normal saline (5 ml of saline in Tanzania), immediately stored at 4–8°C, and processed within two hours of collection. They were centrifuged for 10 minutes at 1,000 rpm (3,500 rpm in Tanzania), and the resulting supernatant and pellet were stored separately at −80°C until shipment. Soluble markers of inflammation in CVLs were quantified by Bio-Plex (Bio-Rad Laboratories NV-SA, Nazareth, Belgium) or ELISA at the ITM in Antwerp in the Vaginal Biomarkers Study [19], and by an in house multiplex bead immunoassay at St. Georges University in London in the Tanzania Study as described previously [20,21].

Samples for microarray testing

Two sterile Copan vaginal swabs were collected per participant per visit. Copan vaginal swabs were shipped frozen to the ITM in Antwerp. Upon arrival, each swab tip was thawed at room temperature for 30 minutes, 1200 μl of diluted PBS (1 PBS: 9 saline, pH 7.4) was added, and the sample was vortexed for 15 seconds. An aliquot of 600 μl was sent on dry ice to TNO (Zeist, the Netherlands) for phylogenetic microarray analysis. We could not test all available samples by microarray due to funding constraints. In the Vaginal Biomarkers Study, only the samples from 216 women that were available in October 2011 were analyzed. These 216 women each contributed one enrollment sample, 34 women also contributed 61 follow-up samples, and two samples had missing clinical data. Tanzania Study participants contributed 20 enrollment and 20 follow-up samples. After excluding six poor quality samples, the total sample size was 313 samples from 232 women. All 313 samples were used for phylogenetic clustering and ecological analyses, but only one enrollment sample per woman with clinical data (N = 230) was used for all other analyses.

Microarray testing

The phylogenetic microarray (V-Chip, TNO, Zeist, The Netherlands) has been described previously [3,22]. Briefly, it contained 283 DNA hybridization probes that generated a consistent signal with a signal/background (S/B) ratio > 5, of which 74 16S probes were species-specific, 60 16S probes targeted multiple bacterial species within one genus, 42 16S probes were specific at family or order level, 85 targeted higher taxonomic levels, five were groEL probes, 14 were 18S probes, and three were viral probes. We focused our clustering analyses on these 283 probes, and all additional analyses on the 134 16S probes generating species or genus-specific signals. A probe targeting a bacterium classified by the Ribosomal Database Project as an uncultured bacterium in the Lachnospiraceae family matched perfectly with a bacterium recently named BV-associated bacterium 1 (BVAB1) in Genbank (Genbank entry AY724739.1) [23]. We refer to it as BVAB1 and included it in the 134 species/genus-specific probes.

Microarray sample preparation and labeling, amplification and hybridization were described elsewhere [9,22]. Imagene 5.6 software (BioDiscovery, Marina del Rey, USA) was used to read the scanned results and quantify the signal (S) and the background (B). Ratios for S and B were calculated and if S was not confidently above B (S > B + 2*standard deviation (SD) of B), the S/B ratio was set to 1. Slide normalization was performed by Lowess smoothing [24]. We used normalized S/B ratios to estimate bacterial loads, referred to from here onwards as ‘abundance’.

Statistical analysis

Clustering analysis was performed using Python 2.7.2 version (Python Software Foundation, Other statistical analyses were performed using R 3.0.2 version (R Development Core Team, 2013), MATLAB (R2012a, The Math Works, Natick, USA), STATA release 12 (StataCorp, TX, USA) and MS Excel (Microsoft, Washington, USA).

Neighborhood co-regularized multi-view spectral clustering of normalized S/B ratios was used to identify VMB clusters as described before [3,25]. These clusters were named KRST-I to KRST-V (with KRST denoting Kenya, Rwanda, South Africa, and Tanzania). For each sample, the probability of belonging to a particular cluster was calculated by probability decomposition of the co-occurrence matrix. A cut-off probability of 70% was used to assign samples to a cluster. For each cluster, the following microbial ecology parameters were computed: richness (median number of genera), evenness (expressed as a community organization value (Co-value) with 0 representing complete evenness and 100 complete unevenness [26], and the Shannon diversity index [27]. We focused the evenness calculations on the five most abundant bacteria in each cluster to reduce the influence of the long tail of minority species [26]. To compare cumulative Co-values per cluster, an average sample per cluster was generated by calculating median S/B ratios per genus across the samples in that cluster.

To assess associations between VMB clusters and potential correlates, we only included one sample per woman with clinical data (N = 230) and excluded women who could not be assigned to a cluster with at least 70% probability (N = 22). To improve statistical power, the three dysbiotic clusters (KRST-III, KRST-IV and KRST-V) were pooled, and the pooled cluster is referred to as KRST-pIII-V (p for pooled). In regression models, KRST-II and KRST-pIII-V were each compared to KRST-I. VMB correlates were grouped into three groups: 1) sociodemographics, sexual behavior and reproductive history; 2) self-reported symptoms, clinician-observed findings and antibiotic use; and 3) cervicovaginal immunology. Unadjusted associations were assessed by one way ANOVA for continuous variables and Fisher’s exact test for categorical variables; p values were adjusted for false discovery using the linear step up Benjamini-Hochberg procedure, assuming a false discovery rate q = 0.1 and a significance level α = 0.1. An adjusted p value of 0.01 (α = 0.1 * q = 0.1) was considered statistically significant. Adjusted associations were assessed by multinomial logistic regression models with stepwise backward elimination using p ≤ 0.2 as the cut-off in unadjusted models. The final model was selected based on the smallest Akaike Information Criteria (AIC). Microarray testing was done in two batches but we found no evidence for a batch effect in our analyses.


VMB clusters

We identified five VMB clusters, which are visualized in a co-occurrence matrix (Additional file 1: Figure S1A-C). Thirty-five samples from 22 women had a probability of <70% belonging to one of the five clusters (Additional file 1: Figure S1A). These samples did not cluster together, and the sociodemographic, behavioral and clinical characteristics of the 22 women did not differ significantly from those of the 208 women who were assigned to a cluster (Additional file 1: Tables S1 and S2).

The VMB clusters were characterized using ecological parameters (below), Co-values (Figure 1), and a heatmap of S/B ratios of bacteria that were most abundant in this study or have been reported as important in previous studies [3] (Figure 2). Cluster KRST-I was dominated by L. crispatus and did not contain other bacterial taxa in high abundance. Only 19 women (9%) were assigned to this cluster. Cluster KRST-II was dominated by L. iners, but some samples also had high abundance of Gardnerella vaginalis, Atopobium vaginae, and Prevotella spp (Figure 2). Cluster KRST-II included the majority of women (n = 136, 65.4%). Clusters KRST-III, IV and V each contained multiple anaerobic species in high abundance (most notably G. vaginalis, A. vaginae, and Prevotella spp.) but in different proportions, and a lower abundance of lactobacilli than clusters KRST-I and II. About a quarter (25.5%) of the women were assigned to the combined cluster KRST-pIII-V. Cluster KRST-III had the highest abundance of Dialister, Megasphaera spp. and Mobiluncus spp. and the lowest abundance of L. iners. Cluster KRST-IV had the highest abundance of Prevotella spp. and the lowest abundance of BVAB1 and Megasphaera spp. Cluster KRST-V contained a higher abundance of L. iners and other Lactobacillus spp. than clusters KRST-III and IV.

Figure 1
figure 1

Median abundance and evenness per cluster. (a) For each median sample per cluster, the cumulative co-values for the five most abundant genera in descending order are plotted. In (b) the actual genera present in each cluster are shown; this clearly illustrates the dominance of lactobacilli in clusters KRST-I and II, and a more even distribution of the other 18 most abundant anaerobic genera in other clusters. Cluster KRST-II also contained Atopobium vaginae and Prevotella spp., but these are not visible in (b) because their abundance in most samples was low. They are visible in the heatmap in Figure 2.

Figure 2
figure 2

Bacterial composition of the cervicovaginal microbiome clusters. (a) shows the five clusters identified by Neighborhood Co-regularized Multi-view Spectral Clustering of microarray data. The white space in between clusters represents samples with less than 70% probability of belonging to a cluster. (b) Heatmap, showing normalized S/B ratios of 33 bacterial species/genera, including the most abundant species/genera per cluster and those traditionally known to be associated with BV and frequently reported in literature. 1Abbreviated probe name: additional targeted species in the same genus are provided in supplementary information of previously published work [3]. 2Abbreviated for Prevotella.

The bacterial diversity increased from cluster KRST-I (median Shannon index of 1.0) to cluster KRST-II (1.2) to clusters KRST-III, IV and V (2.1, 2.0, and 2.0, respectively; p <0.01). The evenness based on the five most abundant genera in each cluster was lower in clusters KRST-I and KRST-II compared to cluster KRST-pIII-V (Figure 1a). The richness in clusters KRST-I and KRST-II (a median of four and five genera per sample) was lower than in clusters KRST-III, IV and V (16, 14, and 14 genera per sample, respectively) (Figure 1b).

Sociodemographic, behavioral and reproductive history correlates of VMB clusters

VMB clustering was not associated with sociodemographic, behavioral or reproductive history characteristics, with the exception of time elapsed since last delivery (Table 1). With every additional month since last delivery, women were less likely to be assigned to cluster KRST-II than cluster KRST-pIII-V (p = 0.01). We did not conduct multivariable modeling with this group of variables because few of them were significant at p ≤ 0.01 in bivariable models.

Table 1 Sociodemographic, behavioral, and reproductive history correlates of VMB clusters

Associations between VMB clusters, Nugent scores and Amsel criteria

Diagnosis of BV by Nugent score and Amsel criteria were each strongly associated with VMB clustering (p < 0.0001) (Table 2, Figure 3). The majority (>88%) of women in clusters KRST-I and II were diagnosed as BV-negative by Nugent score and Amsel criteria. However, while the majority of women (89.6%) in cluster KRST-pIII-V was diagnosed as BV-positive by Nugent score, only 30% were diagnosed as BV-positive by Amsel criteria. Only 9.6% of women were diagnosed with intermediate microbiota by Nugent score and approximately equal proportions were found in each VMB cluster (Figure 3).

Table 2 Clinical correlates of VMB clusters.
Figure 3
figure 3

BV status distribution by Nugent score (a) and Amsel criteria (b) in the three VMB clusters. In Figure 3a, women were classified as BV-negative (BV-), BV intermediate (BVint), or BV-positive (BV+) based on Nugent score of 0–3, 4–6 and 7–10, respectively. In Figure 3b, women were classified as BV-positive if any three of the four Amsel criteria (i.e. clue cells > 20%, positive whiff test, vaginal pH > 4.5, and/or presence of unusual vaginal discharge) was positive. The error bars correspond to the 95% confidence interval of the proportions.

Clinical correlates of VMB clusters

The proportions of women with laboratory-confirmed bacterial STIs (5.6% in KRST-I, 14.7% in KRST-II and 22.6% in KRST-pIII-V), HSV-2 (21.1%, 37.5%, and 37.7%) and HIV (0%, 5.2%, and 5.7%) increased from cluster KRST-I to KRST-pIII-V, indicating a trend although the evidence for an association was weak (ptrend = 0.07 for STIs, 0.34 for HSV2 and 0.44 for HIV) (Table 2, Additional file 1: Figure S2). No women in cluster KRST-I had an UTI compared to 55.6% in cluster KRST-II and 23.7% in cluster KRST-pIII-V (p = 0.06; Table 2). Women with abundant cervical mucus upon speculum examination were more likely to belong to cluster KRST-I, but the total number of women with abundant cervical mucus was small (p = 0.04; Table 2). Women who reported unusual vaginal discharge were more likely to belong to cluster KRST-pIII-V (p = 0.05; Table 2). There was no evidence for an association between VMB clustering and any other clinical characteristics. We did not conduct multivariable modeling with this group of variables due to the fact that few of them were significant at p ≤ 0.01 in bivariable models.

Immunological correlates of VMB clusters

The bivariable models showed strong evidence for an increase of the pro-inflammatory cytokines and growth factors IL-1α, IL-1β, G-CSF, and GM-CSF, and a decrease of IP-10, when moving from cluster KRST-I to KRST-pIII-V (Table 3). The MIP-1β concentration was highest in cluster KRST-II and lowest in cluster KRST-I. Six variables (IL-1α, IL-1β, G-CSF, GM-CSF, IP-10, and MIP-1β) qualified for inclusion in the multivariable model, but the final model (AIC = 296) included only four of them (IL-1α, G-CSF, IP-10, and MIP-1β). IL-1β, which was highly significant in the bivariable analysis, was eliminated due to its high correlation with IL-1α (Pearson’s correlation coefficient = 0.65). The associations in the final model were in the same direction as in the bivariable models.

Table 3 Immunological correlates of VMB clusters


We identified five vaginal microbiota clusters in a diverse group of women from four African countries. Two of these clusters were dominated by L. crispatus or L. iners, while the other three clusters consisted of various combinations of other (facultative) anaerobic bacteria in addition to L. iners. Studies in the USA, Europe and Asia have also reported clusters dominated by L. crispatus and L. iners, but L. crispatus clusters were typically more common and L. iners clusters less common than in our study [8,10,28,29]. Longitudinal studies have shown that L. crispatus protects women from vaginal dysbiosis more efficiently than L. iners (reviewed in 6). We did not identify any clusters dominated by L. gasseri or L. jensenii. Such clusters have been reported in studies in the USA, Europe and Asia but were less common than clusters dominated by L. crispatus or L. iners [8,10,28,29]. We found that Nugent scoring correlated well with molecular VMB clustering but the Amsel criteria did not, as has been found by others [reviewed in 6]. This may have clinical relevance since the Amsel criteria are often used in clinic settings to diagnose BV.

Women in the pooled dysbiotic cluster KRST-pIII-V had higher concentrations of several pro-inflammatory factors than women in clusters KRST-I and II. In the past, BV has often been described as a non-inflammatory syndrome, but more recent data consistently show that vaginal dysbiosis is associated with subclinical cervicovaginal immune activation [19,30]. The differences in inflammatory markers could also partially be explained by hormonal differences between the VMB clusters such as stage of the menstrual cycle, use of hormonal contraception and pregnancy. A recent systematic review showed that hormonal contraceptive use and pregnancy are negatively associated with vaginal dysbiosis [31], and it has been suggested that pregnancy is associated with a pro-inflammatory state [32]. However, in this study, hormonal contraceptive use and pregnancy were not associated with VMB clustering and a sensitivity analysis excluding the pregnant women and women from Tanzania (for whom we did not know the stage of the menstrual cycle at the time of sampling) did not change our results substantively (see footnotes Table 3). Inflammation related to vaginal dysbiosis is of concern because vaginal dysbiosis is very common (25.5% in this study) [33] and inflammation in the genital tract results in the attraction of CD4+ target cells for HIV [20] as well as shedding of HIV in HIV-positive women [34]. There was no evidence for an association between VMB clustering and HIV prevalence in this study (likely due to the small number of HIV infections) but we recently showed a strong association (using the same phylogenic microarray as in this study) in female sex workers in Kigali, Rwanda [3].

There was some evidence that a lower proportion of women in cluster KRST-I tested positive for STIs (ptrend = 0.07) and UTIs (p = 0.06) than in the other clusters, but the number of cases was small. These findings are, however, consistent with the findings of the above-mentioned study in female sex workers in Kigali that reported negative associations between a L. crispatus-dominated VMB and various STIs [3]. Other studies have also shown that lactobacilli-dominated VMB [reviewed in 6] or a Nugent score 0–3 [4,5,35] are negatively associated with both viral and bacterial STIs. The relationship between the VMB and UTIs has not been adequately studied, but recent studies have shown that the urine microbiome in women resembles their VMB [36]. In contrast, a higher proportion of women in clusters KRST-I and II had vaginal candidiasis than women in the other clusters (ptrend = 0.09), and this is consistent with findings of several other studies (reviewed in [6]).

We identified few additional correlates of VMB composition, perhaps due to limited statistical power. Women who had delivered a baby more recently were more likely to belong to clusters KRST-I and KRST-II (p = 0.01), whereas women who reported unusual vaginal discharge were more likely to belong to clusters KRST-III/IV/V (p = 0.05). The former could be explained by sex hormone levels but this seems unlikely since the median time period since the last delivery was at least 33 months in all three groups. Confounding by age or sexual behavior might be more likely. There was some evidence that women in the KRST-I cluster were older, had fewer sexual partners in the last three months, but more unprotected sex with a steady partner, than women in the other clusters.

Our study used samples and data from two studies. While we used the same VMB assessment methods in both studies, there were some differences in other assessments such as the way questions were asked, the test kits that were used for on-site diagnostic testing, and the platforms that were used in Antwerp and London for cytokine testing [37]. Other limitations include the cross-sectional nature of the study, the small sample sizes in some of the comparison groups (particularly the number of HIV and bacterial STI cases), and the imprecise timing of sexual behaviors (in the Vaginal Biomarkers Study), STI and UTI diagnoses (in the Vaginal Biomarkers Study) and stage of the menstrual cycle (in the Tanzania study) relative to microarray sampling.


Vaginal dysbiosis in African women was significantly associated with vaginal inflammation. The associations with increased prevalence of STIs (including HIV) and UTIs, and decreased prevalence of vaginal candidiasis, should be confirmed in larger studies.


  1. Redondo-Lopez V, Cook RL, Sobel JD. Emerging role of lactobacilli in the control and maintenance of the vaginal bacterial microflora. Rev Infect Dis. 1990;12:856–72.

    Article  CAS  PubMed  Google Scholar 

  2. Amsel R, Totten PA, Spiegel CA, Chen KC, Eschenbach D, Holmes KK. Nonspecific vaginitis. Diagnostic criteria and microbial and epidemiologic associations. Am J Med. 1983;74:14–22.

    Article  CAS  PubMed  Google Scholar 

  3. Borgdorff H, Tsivtsivadze E, Verhelst R, Marzorati M, Jurriaans S, Ndayisaba GF, et al. Lactobacillus-dominated cervicovaginal microbiota associated with reduced HIV/STI prevalence and genital HIV viral load in African women. ISME J. 2014; doi:10.1038/ismej.2014.26.

  4. Cherpes TL, Meyn LA, Krohn MA, Lurie JG, Hillier SL. Association between acquisition of herpes simplex virus type 2 in women and bacterial vaginosis. Clin Infect Dis. 2003;37:319–25.

    Article  PubMed  Google Scholar 

  5. Watts DH, Fazzari M, Minkoff H, Hillier SL, Sha B, Glesby M, et al. Effects of bacterial vaginosis and other genital infections on the natural history of human papillomavirus infection in HIV-1-infected and high-risk HIV-1-uninfected women. J Infect Dis. 2005;191:1129–39.

    Article  PubMed  Google Scholar 

  6. van de Wijgert JH, Borgdorff H, Verhelst R, Crucitti T, Francis S, Verstraelen H, et al. The vaginal microbiota: what have we learned after a decade of molecular characterization? PLoS One. 2014;9:e105998.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Nugent RP, Krohn MA, Hillier SL. Reliability of diagnosing bacterial vaginosis is improved by a standardized method of gram stain interpretation. J Clin Microbiol. 1991;29:297–301.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SS, McCulle SL, et al. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci U S A. 2011;108:4680–7.

    Article  CAS  PubMed  Google Scholar 

  9. Fiscella K, Klebanoff MA. Are racial differences in vaginal pH explained by vaginal flora? Am J Obstet Gynecol. 2004;191:747–50.

    Article  PubMed  Google Scholar 

  10. Zhou X, Hansmann MA, Davis CC, Suzuki H, Brown CJ, Schütte U, et al. The vaginal bacterial communities of Japanese women resemble those of women in other racial groups. FEMS Immunol Med Microbiol. 2010;58:169–81.

    Article  CAS  PubMed  Google Scholar 

  11. van de Wijgert J, Shattock R. Vaginal microbicides: moving ahead after an unexpected setback. Editorial review. AIDS. 2007;21:2369–76.

    Article  PubMed  Google Scholar 

  12. Jespers V, Millwood IY, Poynten IM, Van Damme L, Kaldor JM. The evolving design and methods for trials evaluating the safety and candidate vaginal microbicides: a systematic review. Sex Transm Dis. 2013;40:729–36.

    Article  CAS  PubMed  Google Scholar 

  13. Fichorova RN, Tucker LD, Anderson DJ. The molecular basis of nonoxynol-9-induced vaginal inflammation and its possible relevance to human immunodeficiency virus type 1 transmission. J Infect Dis. 2001;184:418–28.

    Article  CAS  PubMed  Google Scholar 

  14. Malonza IM, Mirembe F, Nakabiito C, Odusoga LO, Osinupebi OA, Hazari K, et al. Expanded phase I safety and acceptability study of 6% cellulose sulfate vaginal gel. AIDS. 2005;19:2157–63.

    Article  CAS  PubMed  Google Scholar 

  15. Doh AS, Ngoh N, Roddy R, Lai JJ, Linton K, Mauck C. Safety and acceptability of 6% cellulose sulfate vaginal gel applied four times per day for 14 days. Contraception. 2007;76:245–9.

    Article  CAS  PubMed  Google Scholar 

  16. Jespers V, Crucitti T, Menten J, Verhelst R, Mwaura M, Mandaliya K, et al. Prevalence and correlates of bacterial vaginosis in different subpopulations of women in sub-Saharan Africa: a cross-sectional study. PLoS One. 2014;9:e109670.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Francis SC, Baisley K, Hou Y, Herrera C, Andreasen A, Watson-Jones D, et al. Expression of soluble immune proteins in the female genital tract: Background variation among women at high risk for HIV in Tanzania. Abstract A-641-0263-09932, 20th International AIDS Conference, Melbourne, Australia.

  18. Kapiga SH, Ewings FM, Ao T, Chilongani J, Mongi A, Baisley K, et al. The epidemiology of HIV and HSV-2 infections among women participating in microbicide and vaccine feasibility studies in Northern Tanzania. PLoS One. 2013;8:e68825.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Kyongo JK, Crucitti T, Menten J, Hardy L, Cools P, Michiels J, et al. Genital tract immunological markers in Sub-Saharan African women with relevance to HIV risk and prevention. Clin Vaccine Immunol, in press.

  20. Biancotto A, Grivel JC, Iqlehart SJ, Vanpouille C, Lisco A, Sieg SE, et al. Abnormal activation and cytokine spectra in lymph nodes of people chronically infected with HIV-1. Blood. 2007;109:4272–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Fischetti L, Barry SM, Hope TJ, Shattock RJ. HIV-1 infection of human penile explant tissue and protection by candidate microbicides. AIDS. 2009;23:319–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Dols JA, Reid G, Kort R, Schuren FH, Tempelman H, Bontekoe TR, et al. PCR-based identification of eight Lactobacillus species and 18 hr-HPV genotypes in fixed cervical samples of South African women at risk of HIV and BV. Diagn Cytopathol. 2012;40:472–7.

    Article  PubMed  Google Scholar 

  23. Fredricks DN, Fiedler TL, Marrazzo JM. Molecular identification of bacteria associated with bacterial vaginosis. N Engl J Med. 2005;353:1899–911.

    Article  CAS  PubMed  Google Scholar 

  24. Quackenbush J. Microarray data normalization and transformation. Nat Genet. 2002;32:496–501.

    Article  CAS  PubMed  Google Scholar 

  25. Tsivtsivadze E, Borgdorff H, van de Wijgert JHHM, Schuren FHJ, Verhelst R, Heskes T. Neighborhood co-regularized multi-veiw spectral clustering of microbiome data. Nanjing: Proceedings of 2nd workshop on partially supervised learning; 2013.

    Google Scholar 

  26. Marzorati M, Maignien L, Verhelst A, Luta G, Sinnott R, Kerckhof FM, et al. Barcoded pyrosequencing analysis of the microbial community in a simulator of the human gastrointestinal tract showed a colon region-specific microbiota modulation for two plant-derived polysaccharide blends. Antonie Van Leeuwenhoek. 2013;103:409–20.

    Article  CAS  PubMed  Google Scholar 

  27. Shannon CE. A mathematical theory of communication. Bell Syst Tech J. 1948;27:379–423.

    Article  Google Scholar 

  28. Jespers V, Menten J, Smet H, Poradosu S, Abdellati S, Verhelst R, et al. Quantification of bacterial species of the vaginal microbiome in different groups of women, using nucleic acid amplification tests. BMC Microbiol. 2012;12. doi: 10.1186/1471-2180-12-83.

  29. Vasquez A, Jakobsson T, Ahrne S, Forsum U, Molin G. Vaginal lactobacillus flora of healthy Swedish women. J Clin Microbiol. 2002;40:2746–9.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Doerflinger SY, Throop AL, Herbst-Kralovetz MM. Bacteria in the vaginal microbiome alter the innate immune response and barrier properties of the human vaginal epithelia in a species-specific manner. J Infect Dis. 2014;209:1989–99.

    Article  CAS  PubMed  Google Scholar 

  31. van de Wijgert JH, Verwijs MC, Turner AN, Morrison CS. Hormonal contraception decreases bacterial vaginosis but oral contraception may increase candidiasis: implications for HIV transmission. AIDS. 2013;27:2141–53.

    Article  PubMed  Google Scholar 

  32. Anderson BL, Mendez-Figueroa H, Dahlke JD, Raker C, Hillier SL, Cu-Uvin S. Pregnancy-induced changes in immune protection of the genital tract: defining normal. Am J Obstet Gynecol. 2013;208:321–9.

    PubMed  PubMed Central  Google Scholar 

  33. Low N, Chersich MF, Schmidlin K, Egger M, Francis SC, van de Wijgert JH. Intravaginal practices, bacterial vaginosis, and HIV infection in women: individual participant data meta-analysis. PLoS Med. 2011;8:e1000416. doi:10.1371/journal.pmed.1000416.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Sha BE, Zariffard MR, Wang QJ, Chen HY, Bremer J, Cohen MH, et al. Female genital-tract HIV load correlates inversely with Lactobacillus species but positively with bacterial vaginosis and Mycoplasma hominis. J Infect Dis. 2005;191:25–32.

    Article  PubMed  Google Scholar 

  35. Spear GT, Sikaroodi M, Zariffard MR, Landay AL, French AL, Gillevet PM. Comparison of the diversity of the vaginal microbiota in HIV-infected and HIV-uninfected women with or without bacterial vaginosis. J Infect Dis. 2008;198:1131–40.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Siddiqui H, Nederbragt AJ, Lagesen K, Jeansson SL, Jakobsen KS. Assessing diversity of the female urine microbiota by high throughput sequencing of 16S rDNA amplicons. BMC Microbiol. 2011;11:244.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Fichorova RN, Richardson-Harman N, Alfano M, Belec L, Carbonneil C, Chen S, et al. Biological and technical variables affecting immunoassay recovery of cytokines from human serum and simulated vaginal fluid: a multicenter study. Anal Chem. 2008;80:4741–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


The authors thank the study participants, the Vaginal Biomarkers Study Group, and Nynke van Berkum and TNO analysts for microarray testing. The Vaginal Biomarkers Study Group consists of: ICRH, Mombasa, Kenya: Kishor Mandaliya, Lou Dierick, Mary Mwaura, Walter Jaoko, Eunice Irungu, Christine Katingima, Mercy Maina, Jane Wanjiru Mazera, Josephine Gichuru, Grace Aketch Onuki, Mary Kiambi, Mary Thiong’o, Salome Wanjiku, Patricia Nduku, Carol Njeru, Bernard Mbogho, Sammy Wambua, Rachel Sidi Baya, Emmanuel Moffat Onduko, Patrick Katana Kombo, Simon Chengo Masha, Mary Ndinda John, Kevin Odeyo, Dora Ngala, Collins Odero. Wits Reproductive Health & HIV Institute, Johannesburg, South Africa: Sinead Delany-Moretlwe, Vinodh Aroon Edward, Krishnaveni Reddy, Nina Von Knorring, Ishania Mahabeer, Johannah Nkoleleng Mashilo, Ntombifuthi Mnyandu, Keneuoe Mokoatle, Siyabulela Nani, Gugu Tshabalala, Thembisile Hope Mngwevu, Noxolo Mtabane, Puseletso Maria Masalesa, Zodidi Kumase, Sefora Dipolelo Mohale, Mavis Mantshitseng Madi, Mandla Mlotshwa, Pholo Wilson Maenetje, Nishanee Arjun, Debra De Assis Rosa. Rinda Ubuzima, Kigali, Rwanda: Gilles F. Ndayisaba, Evelyne Kestelyn, Ammiel Gasarabwe, Servaas Van Eeckhoudt, Stephen Agaba, Rosette Busasa, Deogratias Nshimuyimana, Grace Umutoni, Vincent Karangwa, Claire Bukuru, Alice Fiat, Lambert Mwambarangwe, Viateur Musengamana, Jeanine Nyinawabega. ITM, Antwerp, Belgium: Vicky Jespers, Liselotte Hardy, Tania Crucitti, Joris Menten, Céline Schurmans, Harry van Loen, Anne Buvé, Jordan Kyongo, Kevin Ariën, Guido Vanham, Said Abdellati, Vicky Cuylaert, Wendy Thys, An Ielegems, Lieve Casier. ICRH University Ghent, Ghent, Belgium: Hans Verstraelen, Marleen Temmerman, Rita Verhelst, Mario Vaneechoutte, Piet Cools, Bart Saerens. AIGHD/AMC, Amsterdam, The Netherlands: Janneke van de Wijgert, Friso Janssen. MRC CTU, London, United Kingdom: Sheena McCormack, Sarah Joseph. LSHTM, London, United Kingdom: Suzanna Francis, Kathy Baisley, Richard Hayes. MITU, Mwanza, Tanzania: Saidi Kapiga, Aura Andreasen, John Changalucha, Kaballa Maganja, Clemens Masesa.

This work was funded by the European and Developing Countries Clinical Trials Partnership (EDCTP; grant numbers IP_2007_33070_001 and CT_ct_05_32070_002), the Research Foundation Flanders (FWO - 3GA13210), Aids Fonds Netherlands (project number 201102), the European Commission (Framework Programme 7, CHAARM consortium), the International Partnership for Microbicides, and the UK Medical Research Council (fellowship to S. Francis, G1002369). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

The underlying data that were used for the analyses presented in this paper can be obtained from Professor JHHM van de Wijgert via email:

Author information

Authors and Affiliations



Corresponding author

Correspondence to Janneke HHM van de Wijgert.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived of and designed the study: JvdW, VJ, RV, TC, FS. Coordinated clinical data and sample collection: SCF, MM, SMD, GN, LH. Coordinated diagnostic testing: TC. Coordinated microarray testing: FS. Conducted immunology testing: JKK. Conducted statistical and bioinformatics analyses: RG, HB, JM, ET. All authors provided input in the analyses and manuscript writing and approved the final manuscript.

Additional file

Additional file 1:

Supplementary methods, tables and figures.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gautam, R., Borgdorff, H., Jespers, V. et al. Correlates of the molecular vaginal microbiota composition of African women. BMC Infect Dis 15, 86 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: