Correlates of the molecular vaginal microbiota composition of African women
- Raju Gautam1,
- Hanneke Borgdorff2,
- Vicky Jespers3,
- Suzanna C Francis4,
- Rita Verhelst5,
- Mary Mwaura6,
- Sinead Delany-Moretlwe7,
- Gilles Ndayisaba8,
- Jordan K Kyongo3,
- Liselotte Hardy3,
- Joris Menten3,
- Tania Crucitti3,
- Evgeni Tsivtsivadze9,
- Frank Schuren9,
- Janneke HHM van de Wijgert1, 8, 10Email author and
- for the Vaginal Biomarkers Study Group
© Gautam et al.; licensee BioMed Central. 2015
Received: 19 August 2014
Accepted: 11 February 2015
Published: 21 February 2015
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.
KeywordsBacterial vaginosis Vaginal microbiome Vaginal microbiota Lactobacillus Candidiasis Sexually transmitted infections Urinary tract infections Women HIV Africa
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 . The clinical conditions caused by an imbalanced VMB include bacterial vaginosis (BV) and bacterial vaginitis . 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 . 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) . 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 . 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 .
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).
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 , 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.
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) . 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 . We used normalized S/B ratios to estimate bacterial loads, referred to from here onwards as ‘abundance’.
Clustering analysis was performed using Python 2.7.2 version (Python Software Foundation, https://www.python.org/download/releases/2.7.2/). 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 , and the Shannon diversity index . 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 . 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.
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 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
Sociodemographic, behavioral, and reproductive history correlates of VMB clusters
KRST-I n (%) N = 19
KRST-II n (%) N = 136
KRST-pIII-V n (%) N = 53
Sociodemographic characteristics at screening
Adult women, KE
Adolescent women, KE
Pregnant women, KE
Adult women, RSA
IVP users, RSA + TZ
HIV-positive and high risk women, RW
16 -17 years
18 -24 years
25 -29 years
30 -35 years
Socio-economic status (composite)a:
Sexual risk taking (composite)b:
Number of lifetime sexual partnersc:
Two to three
More than three
Number of sexual partners last 3 monthsd:
Vaginal sex previous morning/eveninge:
Frequency of sex last 3 monthsf:
11 - 30 times
Condom use last sexual contactf:
Freq of unprotected sex (last 3 months)g:
Put something in vagina to dry/tighten vagina before sex (enrollment)h
Clean vagina after sex (enrollment):
Vaginal practice for washing:
Yes, but not the evening before visit
Yes, including evening before visit
Products used to wash/clean/dry vagina
Water or fingers only
Water + soap
Reproductive history at screening
Median gravidity (IQR)
2 (1, 3)
2 (1, 3)
2 (1, 3)
Median parity (IQR)
2 (1, 2.5)
1.5 (0, 2)
1 (0, 2)
Median months since last deliveryi† (IQR)
41 (28, 54)
33 (19, 48)
56 (28, 82)
Current contraceptive use:
None, not pregnant
Combined oral contraceptives
Condoms + IUD
Associations between VMB clusters, Nugent scores and Amsel criteria
Clinical correlates of VMB clusters.
KRST-I n (%) N = 19
KRST-II n (%) N = 136
KRST-pIII-V n (%) N = 53
Self reported symptoms at screening
Vaginal discharge reported:
Clinician observed findings at screening
Cervical and/or vaginal epithelial abnormalities:
Any colposcopy finding?:
Laboratory-confirmed reproductive and urinary tract infections
BV by Nugent score
BV by Amsel criteria
Candidiasis on wet mountd
Urinary tract infection by diptstick test
Treatment in the 14 days prior to enrollment
Any systemic antibioticsd:
Any vaginal antibiotics: No
Bacterial vaginosis requiring treatment:
Candidiasis requiring treatmente:
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
Immunological correlates of VMB clusters
KRST-I (n = 19)
KRST-II (n = 136)
KRST-pIII-V (n = 53)
Mean (95% CI)
OR/aOR (95% CI)
Mean (95% CI)
OR (95% CI)
aOR (95% CI)
Mean (95% CI)
OR (95% CI)
aOR (95% CI)
1.08 (0.86, 1.3)
1.22 (1.14, 1.31)
1.61 (0.69, 3.77)
1.58 (0.58, 4.33)
1.64 (1.48, 1.81)
7.4 (2.63, 20.82)
11.5† (3.29, 40.22)
0.65 (0.33, 0.97)
0.83 (0.71, 0.95)
1.4 (0.72, 2.74)
1.34 (1.11, 1.57)
3.33 (1.58, 7.03)
4.65 (4.36, 4.94)
4.83 (4.74, 4.91)
1.9 (0.76, 4.76)
4.87 (4.74, 5)
2.28 (0.79, 6.6)
0.64 (0.33, 0.95)
0.9 (0.79, 1)
1.68 (0.85, 3.3)
0.9 (0.68, 1.11)
1.68 (0.8, 3.54)
2.06 (1.83, 2.28)
2.19 (2.09, 2.28)
1.41 (0.62, 3.2)
2.25 (2.04, 2.45)
1.65 (0.68, 3.97)
0.07 (−0.18, 0.32)
0.13 (0.05, 0.21)
1.31 (0.45, 3.81)
0.24 (0.12, 0.35)
2.11 (0.66, 6.7)
2.54 (2.36, 2.73)
2.62 (2.51, 2.73)
1.22 (0.57, 2.6)
0.98 (0.39, 2.43)
2.15 (1.93, 2.38)
0.47 (0.21, 1.06)
0.16† (0.05, 0.47)
0.26 (0.05, 0.47)
0.32 (0.25, 0.39)
1.38 (0.47, 4.03)
0.13 (−0.01, 0.27)
0.58 (0.19, 1.78)
1.84 (1.57, 2.11)
1.91 (1.79, 2.02)
1.16 (0.57, 2.38)
0.71 (0.26, 1.92)
2.08 (1.93, 2.23)
1.83 (0.81, 4.15)
3.6†† (1.02, 12.73)
0.69 (0.23, 1.15)
0.94 (0.83, 1.05)
1.57 (0.84, 2.94)
1.73 (0.79, 3.82)
0.7 (0.48, 0.93)
1.02 (0.52, 1.98)
0.61 (0.25, 1.52)
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 , and it has been suggested that pregnancy is associated with a pro-inflammatory state . 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)  and inflammation in the genital tract results in the attraction of CD4+ target cells for HIV  as well as shedding of HIV in HIV-positive women . 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 .
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 . 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 . 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 ).
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 . 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.
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: email@example.com.
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Fiscella K, Klebanoff MA. Are racial differences in vaginal pH explained by vaginal flora? Am J Obstet Gynecol. 2004;191:747–50.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- van de Wijgert J, Shattock R. Vaginal microbicides: moving ahead after an unexpected setback. Editorial review. AIDS. 2007;21:2369–76.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Fredricks DN, Fiedler TL, Marrazzo JM. Molecular identification of bacteria associated with bacterial vaginosis. N Engl J Med. 2005;353:1899–911.View ArticlePubMedGoogle Scholar
- Quackenbush J. Microarray data normalization and transformation. Nat Genet. 2002;32:496–501.View ArticlePubMedGoogle Scholar
- 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
- 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.View ArticlePubMedGoogle Scholar
- Shannon CE. A mathematical theory of communication. Bell Syst Tech J. 1948;27:379–423.View ArticleGoogle Scholar
- 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.Google Scholar
- Vasquez A, Jakobsson T, Ahrne S, Forsum U, Molin G. Vaginal lactobacillus flora of healthy Swedish women. J Clin Microbiol. 2002;40:2746–9.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.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/4.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.