Introduction

Adolescent females and young women are at high risk of sexually transmitted infections (STI), particularly in sub-Saharan Africa (SSA)1,2. Young women in SSA are also faced with a high risk of unintended pregnancies3,4, which can be associated with increased maternal and infant mortality and morbidity, especially in developing countries5. Hormonal contraceptives (HC) play a crucial role in preventing unintended pregnancy. However, a number of observational studies have suggested that the use of HCs may influence women’s susceptibility to HIV infection6,7, by altering vaginal microbiota8,9,10,11,12,13,14,15 or immunity16,17,18,19,20,21.

The composition of the vaginal microbiota may impact HIV susceptibility8,9. Highly diverse communities, as well as the relative and absolute abundance of specific organisms, including Parvimonas, Prevotella spp. (such as P. amnii, P. bivia, and P. melaninogenica), Mycoplasma, Gemella asaccharolytica, Sneathia, Veillonella montpellieriensis, and Megasphaera, have been associated with HIV risk9,10. In contrast, lactobacilli, particularly non-iners Lactobacillus spp., have been associated with lower risk of HIV infection9,10,11. High-diversity vaginal communities are also strongly associated with high levels of genital pro-inflammatory cytokine concentrations20,21, which increase the risk for subsequent HIV infection19. In observational cohorts, HCs have been associated with cervical and vaginal inflammation16,17. Combined oral contraceptive (COC) use has been associated with increases in relative abundance of lactobacilli12. There are conflicting data as to whether HCs are protective or permissive for bacterial vaginosis (BV)13,14,15,18.

Estrogen is associated with glycogen deposition in epithelial cells, aiding lactic acid metabolism22,23. The progestin-only injectable depomedroxyprogesterone acetate (DMPA) may cause a hypoestrogenic state, which in turn may impair colonization of lactic acid-producing lactobacilli24. In observational studies, DMPA has been associated with decreased risk of BV; however, there are few data on the effects of the injectable norethisterone enanthate (Net-En) or combined contraceptive vaginal rings (CCVR)13,15. We hypothesized that the long-acting injectable progestin Net-En may alter vaginal microbiota, causing increased genital inflammation, but that COC and CCVR, due to their estrogen component, would not. Therefore, we assessed temporal changes in cervicovaginal microbiota and cytokines in adolescents randomized to three HCs: the long-acting progestin-only injectable Net-En, combined oral contraceptives (COC) Triphasil® or Nordette® and NuvaRing®, a CCVR containing etonorgesterol/ethinyl estradiol. We found that adolescents assigned to COC not only had more lactobacilli compared to those assigned to Net-En, but also those assigned to CCVR. CCVR use was associated with significantly elevated genital tract inflammation compared with both Net-En and COC. Together, these data suggest that COC may benefit women at high risk for BV and HIV.

Results

Cohort characteristics

Of the 180 screened for the parent study, 130 adolescent females (median age 17 years; interquartile range (IQR) 16–18) were enrolled into a randomized study with crossover (Fig. 1)25. Of these, 45 were randomized to NuvaRing® (hereon referred to as CCVR), 45 to Net-En, and 40 to COC (Table 1 and Fig. 1). All particiants consented to enrol in this substudy, which intended to assess mucosal changes induced by contraceptives. A total of 108 adolescents reached the crossover visit at 16 weeks, while 94 adolescents completed the 32-week exit visit. One participant assigned to the Net-En arm seroconverted before the crossover visit, one participant withdrew consent for mucosal sampling at the exit visit, and there was one pregnancy in the CCVR arm before the crossover visit; therefore, these three participants were exited from the substudy. At baseline, adolescents assigned to each study arm were similar in demographics, medical, and reproductive history (including age, body mass index (BMI), vaginal insertion practices and antibiotics use, as well as BV, Candida, and STI prevalence) (Table 1). Reported sexual risk behavior was similar across groups, as was the distribution of baseline HC use.

Fig. 1: Study overview and randomization.
figure 1

CONSORT diagram of the number of participants that completed each study visit and provided genital samples, and the distribution of participants on each of the hormonal contraceptive methods analyzed: combined oral contraceptives (COC), the Net-En injection, and the vaginally inserted combined contraceptive ring (CCVR) at crossover and exit. Note that at crossover, some participants could choose their next contraceptive method if assigned to CCVR as their first method. LTFU loss-to-follow-up.

Table 1 Characteristics of participants at baseline according to randomization arms.

Vaginal microbiota in South African adolescents

Of the 329 vaginal lateral wall samples taken from 130 adolescents at all visits, 319 samples passed sequencing and quality control measures (≥5000 reads/sample, screening: n = 126, crossover: n = 104, and exit: n = 89). Using soft k-means clustering26 with weighted UniFrac as the distance measure, we identified three major microbial community-state types (CSTs) (corresponding to CST-I, CST-III, and CST-IV described by Ravel et al.27, Fig. 2a and Supplementary Fig. 1). The most common community type, CST-IV (n = 145, 45.5%), had high microbial diversity comprising a diverse group of anaerobic bacteria with the most abundant species being Gardnerella vaginalis followed by Lachnovaginosum genomospecies (BVAB1), Megasphaera, L. iners, Prevotella spp. (including P. amnii, P. timonensis, and P. bivia), Atopobium vaginae, Sneathia, and Aerococcus christensenii (Fig. 2a). The latter two communities, CST-I and CST-III, were low-diversity communities dominated by L. crispatus (n = 76, 23.8%) and L. iners (n = 98, 30.7%), respectively (Fig. 2a and Supplementary Fig. 1). The majority of adolescents with a CST-IV community were correspondingly BV positive (Nugent scores 7–10; n = 126, 86.9%), whereas adolescents with CST-I and CST-III communities were mostly BV negative (Nugent 0–3; n = 73, 96.1% and n = 84, 85.7%, respectively). In comparison to studies performed in North American cohorts27,28, the vaginal community types CST-II (dominated by L. gasseri) and CST-V (dominated by L. jensenii) were not present in this cohort, similar to what has been found in other African cohorts20,21. Baseline alpha diversity (i.e., within-sample diversity) was slightly higher in adolescents who were randomized to the Net-En arm compared to those assigned to COC and CCVR (median Shannon Index (SI) 1.78 (IQR 0.64–2.27) vs. 0.87 (0.43–1.99) and 1.46 (0.61–2.14), respectively), but not significantly so (Table 1). Importantly, CSTs according to the randomized arm at baseline were evenly distributed.

Fig. 2: Vaginal community clusters in a South African adolescent cohort.
figure 2

a Barplot depicting the relative abundance of the most abundant bacteria in the baseline samples identified by 16S rRNA microbiome profiling. Samples (n = 126) are grouped by community-state type (CST) established using soft k-means clustering with weighted UniFrac distances and ordered by the most abundant species in each CST (CST-I: L. crispatus, CST-III: L. iners, and CST-IV: G. vaginalis). Alpha diversity for each sample (Shannon Index) is depicted below the barplot. bd Change in CST from baseline to crossover within study arms. Alluvial plot showing the change in CST distribution from baseline to crossover for matched participants randomized to one of the three study arms: b COC (n = 36), c Net-En (n = 32), and d CCVR (n = 34) in an intention-to-treat (ITT) analysis. Each line represents one adolescent and CST changes over time. The color of the line is based on the CST assigned at baseline. eg Alpha diversity at baseline and crossover within study arms. Boxplot representing in alpha diversity (Shannon Index) of vaginal microbiota for matched participants at baseline and crossover randomized to one of the three study arms: e COC (n = 36), f Net-En (n = 32), and g CCVR (n = 34). hj Alpha diversity within participants changing between hormonal contraceptive methods. Boxplots showing the alpha diverisity (Shannon Index) in h 62 vaginal samples from 29 participants changing from COC to Net-En or vice versa, i 52 vaginal samples from 26 participants changing from COC to CCVR or vice versa, and j 130 vaginal samples from 50 participants changing from Net-En to CCVR or vice versa. P values were generated using two-sided paired Wilcoxon signed-rank tests. Y axes are log10-transformed. Significance codes: *P < 0.05. Bounds of boxes show interquartile range (IQR) with the lower and upper hinges corresponding to the 25th and 75th percentiles, respectively, lines in the middle of the box indicate median, and top and bottom whiskers demonstrate value ranges within 1.5 × IQR from the hinge. Points beyond that are plotted individually. CCVR: combined contraceptive vaginal ring, COC: combined oral contraceptives. Source data are provided as a Source Data file.

Impact of HC on the vaginal microbiota

The impact of HC after 16 weeks of contraceptive use was first assessed in an intention-to-treat (ITT) analysis. There were no differences in reported sexual risk behavior between arms at crossover. Furthermore, no differences in the presence of yeast (by microscopy), BV (by Nugent scoring), or in HSV-2 serostatus were observed between study arms at 16 weeks (Table 2).

Table 2 Characteristics of participants with microbiome data at crossover according to assigned study arm (intention-to-treat).

At crossover, after 16 weeks of randomized contraceptive use, a significant difference in the distribution of vaginal CSTs between randomized arms was observed (P = 0.029) (Table 2). This significant difference was driven by the difference between Net-En and COC (P = 0.007). There was some evidence of a difference between Net-En and CCVR (P = 0.080) but no evidence of a difference between COC and CCVR (P = 0.479). The CST-III community type was the most prevalent in the COC arm and significantly more prevalent in this arm compared to the Net-En study arm (adj. P = 0.019, Table 2). In both the Net-En and the CCVR arms, the diverse CST-IV community type was most common followed by the L. crispatus-dominant CST-I community type in the Net-En arm and the L. iners-dominant CST-III community type in the CCVR arm. This corresponded to a high proportion of participants in the COC arm transitioning from CST-I or CST-IV to CST-III (43% and 35%, respectively), while the majority of participants with CST-III at baseline remained stable (67%) (Fig. 2b–d; Supplementary Table 1). In the CCVR arm, 67% of participants with CST-IV at baseline stayed in CST-IV, while 30% and 33% of participants with CST-I or CST-III respectively transitioned to CST-IV. Participants in the Net-En arm remained stable with no transitions away from CST-I and with 72% of participants with CST-IV at baseline not transitioning. However, none of the participants in the Net-En arm with a CST-III community at baseline stayed as such: 57% transitioned to CST-I and 43% transitioned to CST-IV. Using omnibus symmetry exact tests, there was no strong statistical evidence for a lack of symmetry in the CSTs at baseline and at crossover apart from some evidence of a trend for the Net-En group (P = 0.067) where all those with CST-III at baseline transitioned to another state at crossover.

A decrease in alpha diversity was observed in the COC and Net-En arms from baseline to crossover, while an increase was observed in the CCVR arm; however, this did not reach statistical significance (Fig. 2e–g and Table 2). Cross-sectionally, the vaginal microbiota of adolescents randomized to the COC arm had significantly lower alpha diversity at crossover (median SI: 0.90 (IQR 0.35–1.54)) compared to the microbiota of adolescents assigned to either the Net-En (median SI: 1.64 (IQR 0.75–2.22)) or CCVR (median SI: 1.66 (IQR 0.79–1.95), adj. P = 0.009 and P = 0.021, respectively, Table 2). No difference in vaginal alpha diversity was evident at crossover between adolescents randomized to the CCVR and Net-En arms (adj. P = 0.941). Using a linear regression model with Tukey’s post hoc test, the alpha diversity of the vaginal microbiota of adolescents in the COC arm was significantly lower than that of adolescents assigned to the CCVR arm (P = 0.018) and the Net-En arm (P = 0.007). After adjusting for possible confounders (antibiotic use (collinear with baseline STIs)), these differences remained significant (COC vs. CCVR, P = 0.007 and COC vs. Net-En, P = 0.003). We assessed for effect of modification by baseline CST by including interaction terms for study arm and baseline CST in the model; however, these were not significant (P = 0.401), suggesting no evidence that the effect of study arm on vaginal alpha diversity at crossover differs by baseline CST. The contribution of potential mediators, such as menstrual bleeding patterns and sexual risk behavior, of the effects of study arm on the microbiota, was assessed using a mediation analysis of time since the last menstrual period and condom use, respectively (Supplementary Table 2). Only a very small and nonsignificant proportion of the effects of study arm on the microbiota could be explained by the indirect effects of menstrual cycle and sexual risk behavior, with estimates of the proportion mediated ranging from 0.0004 for the time since the last menstrual period in the COC versus CCVR arm to 0.0371 in the Net-En versus CCVR arm (Supplementary Table 2).

To assess differences in community composition cross-sectionally, we performed a permutational multivariate analysis of variance (PERMANOVA) on phylogenetic distances at crossover. We found significant differences in beta diversity between the three assigned study arms at crossover (P = 0.026, R2 = 0.040). Pairwise comparisons revealed that this difference was driven primarily by the differences between COC and Net-En (P = 0.008, R2 = 0.048) and to a lesser degree CCVR and COC (P = 0.055, R2 = 0.031), while there were no significant differences between the Net-En and CCVR arm (P = 0.383, R2 = 0.012). After adjusting for baseline CST, the results from the model remained similar, and the interaction of study arm and baseline CST was not significant (P = 0.443), indicating that there was no effect modification by baseline CST. The differences in beta diversity between the study arms at crossover remained significantly different after adjusting for antibiotic use (P = 0.027, R2 = 0.040) again driven by differences between the COC and Net-En arm (P = 0.011, R2 = 0.048) and the COC and CCVR arm (R2 = 0.031, P = 0.061, Net-En vs. CCVR, R2 = 0.012, P = 0.395).

At the crossover visit, all adolescents who were assigned to either COC or Net-En at baseline were crossed over to CCVR for the final 16 weeks, whereas adolescents assigned to CCVR at baseline were given the choice between COCs or Net-En. Some adolescents changed HC method before mucosal sampling visits due to side effects (bleeding or spotting) or issues with adherence or discomfort with study product (COC: 4/46 (8.7%); Net-En: 0/60 (0.0%); CCVR: 17/80 (21.3%), Fig. 1). Therefore, a per-protocol (PP) analysis was applied to adolescents at crossover (n = 107) and exit (n = 92). Importantly, adolescents using Net-En and CCVR reported less consistent condom use than adolescents on COC, albeit not significantly so (Supplementary Table 3). In this analysis, the vaginal microbiota of adolescents using Net-En or CCVR at the time of sampling was most commonly assigned to the diverse CST-IV community type (43.5% and 45.5%, respectively), whereas participants on COC were fairly equally distributed between CSTs (Supplementary Table 3). Using a linear mixed-effect regression model with Tukey’s post hoc test, no difference in alpha diversity was observed between CCVR and Net-En users (adj. P = 0.980), while the vaginal microbiota of adolescents using COC had a significantly lower alpha diversity compared to the microbiota of those using Net-En or CCVR (adj. P = 0.004 and 0.002, respectively) (Supplementary Table 3). After adjusting for antibiotic use and baseline CST, the alpha diversity remained significantly different between COC users and both Net-En (adj. P = 0.012) and CCVR users (adj. P = 0.009). Once again, the interaction between baseline CST and HC method was nonsignificant (P = 0.772). There was some evidence of differences in microbiota composition (beta diversity) between adolescents on the three different contraceptive methods in the PP analysis; however, these did not reach statistical significance (PERMANOVA P = 0.057, R2 = 0.034).

In a paired analysis of matched samples from 29 adolescents who either changed from COC use to Net-En use or vice versa, the alpha diversity of the vaginal microbiota of adolescents when on COC was significantly lower than when on Net-En (median SI: 0.89 (IQR 0.34–1.52) vs. 1.60 (0.55–2.18), P = 0.028, Fig. 2h). Vaginal alpha diversity of adolescents when on COC was moderately lower than when using CCVR but not significantly so (n = 26, median SI: 0.75 (IQR 0.26–1.44) vs. 1.10 [0.45–2.05], P = 0.162) (Fig. 2i). In contrast, alpha diversity of the vaginal microbiota was similar within participants when using Net-En and CCVR (n = 50, median SI: 1.57 (IQR 0.61–2.11) vs. 1.62 (0.87–2.01), P = 0.643) (Fig. 2j).

Differentially abundant taxa by study arm

We used DESeq2 to identify differentially abundant taxa between (cross-sectionally) and within (pre- vs. post longitudinally) randomized study arms. Fifteen taxa were differentially abundant between vaginal specimens from adolescents in the COC versus Net-En study arms at crossover after adjusting for antibiotic use and sexual risk behavior (condom use) (Table 3). Some Lactobacillus spp. had higher relative abundance in the COC study arm compared with the other two arms, whereas several Prevotella spp. and Mageeibacillus indolicus (BVAB3) were more abundant in the Net-En arm and Prevotella amnii was more abundant in the CCVR arm compared to COC. Prevotella disiens, Gemella asaccharolytica, Mycoplasma hominis, and L. jensenii were more abundant in the CCVR arm compared to the Net-En arm. When comparing changes in relative abundance longitudinally between baseline and crossover in matched samples from the same participant within each randomized arm using DESeq2, we observed a decrease in the abundance of two Clostridium spp. and an increase in L. iners abundances between baseline and crossover within adolescents who initiated COC (Table 3). For adolescents in the Net-En arm, an increase in the abundances of P. disiens, M. hominis, and Campylobacter was observed after Net-En initiation (Table 3), while no taxa differed from baseline to crossover in the CCVR arm.

Table 3 DESeq2 analysis of taxa differentially abundant between (A) randomized study arms at crossover and (B) between baseline and crossover (intention-to-treat analysis).

For bacteria that have previously been associated with high or low HIV risk9,10, we compared read counts standardized to the median read depth in a paired analysis of participants changing between contraceptive methods longitudinally (Supplementary Fig. 2). The read counts of L. iners were significantly higher in participants when using COC than when using Net-En (P = 0.031, adj. P = 0140), while Prevotella was higher in the adolescents when on Net-En compared to COC (P = 0.032, adj. P = 0140, Supplementary Fig. 2a). Mycoplasma (P = 0.019, adj. P = 0.086), Parvimonas (P = 0.017, adj. P = 0.086), and Sneathia (P = 0.020, adj. P = 0.086) were significantly more abundant in adolescents when on CCVR compared to when using COC (Supplementary Fig. 2b). The read counts of L. iners and P. bivia were significantly higher in adolescents when using CCVR than when using Net-En (Supplementary Fig. 2c).

Relationship of HC with bacterial taxa and cervicovaginal cytokines

Of the 15 cytokines measured in cervicovaginal secretions, 13 were reliably detectable, including inflammatory cytokines previously associated with HIV risk (interleukin (IL)-1β, IL-6, IFN-γ, and TNF-α), soluble CD40 ligand (sCD40L), and Th17 cytokines (since Th17 cells have been identified as key target cells for viral replication during vaginal transmission in nonhuman primates)29,30. As shown previously20,21, standardized read counts of taxa associated with increased HIV susceptibility (such as Parvimonas micra, Prevotella spp., Gemella asaccharolytica, Mycoplasma hominis, and Sneathia) correlated positively with the concentrations of several inflammatory cytokines (Supplementary Fig. 3a). In contrast, the abundances of L. crispatus and L. iners were negatively correlated with the concentrations of several cytokines. The canonical Th17 cytokines IL-17A and IL-17F correlated positively with abundance of bacteria not previously thought of as inflammatory in the female genital tract (FGT), including Streptococcus anginosus, Staphylococcus, and Corynebacterium striatum (Supplementary Fig. 3a).

In a longitudinal, paired ITT analysis from baseline to crossover, there were no significant changes in the levels of measured cytokines in the COC or Net-En arms after adjusting for multiple comparisons (Fig. 3a, b). However, the levels of IL-1β, IL-6, IL-21, IL-33, TNF-α, and IFN-γ significantly increased between baseline and crossover in adolescents randomized to CCVR (Fig. 3c and Supplementary Table 4a). Changes in IL-1β, IL-6, IL-21, and TNF-α remained significant after adjusting for multiple comparisons using the Benjamini–Hochberg (BH) method. In an ITT cross-sectional analysis, seven cytokines (IL-1β, IL-6, IL-17A, IL-21, IL-25, and TNF-α) were significantly elevated in adolescents using CCVR versus adolescents using COCs. TNF-α remained significant after adjusting for multiple comparisons (adj P < 0.05, Supplementary Fig. 3b and Supplementary Table 4B). In a linear mixed-effect model with Tukey’s post hoc test adjusted for condom use, baseline alpha diversity and antibiotic use, as well as multiple comparisons, there were no significant differences in cytokine concentrations at crossover between COC and CCVR or Net-En, but multiple cytokines were significantly higher in women assigned to CCVR versus Net-En (IL-21, -25, -31, -33, IFN-γ, and sCD40L).

Fig. 3: Cytokine concentrations from baseline to crossover according to hormonal contraceptive arm in an intent-to-treat analysis.
figure 3

Boxplot showing the change in cytokine concentrations from baseline to crossover for matched participants within each of the three contraceptive arms: a adolescents randomized to COC (n = 36), b adolescents randomized to Net-En (n = 35), and c adolescents randomized to CCVR (n = 34). P values were generated using two-sided paired Wilcoxon signed-rank tests adjusted for multiple comparisons using the Benjamini–Hochberg (BH) method. Y axis log10-transformed. Significance codes: *P < 0.05, **P < 0.01, ***P < 0.001. Bounds of boxes show interquartile range (IQR) with the lower and upper hinges corresponding to the 25th and 75th percentiles, respectively, lines in the middle of the box indicate median, and top and bottom whiskers demonstrate the largest and lowest value within 1.5 × IQR from the hinge. Points beyond that are plotted individually. CCVR: combined contraceptive vaginal ring, COC: combined oral contraceptives. Source data are provided as a Source Data file.

Integration of microbiota and cytokine data

Since both Net-En and CCVR users had higher relative abundance of inflammation-associated taxa than COC users, but inflammation was the highest in CCVR, we performed multivariate logistic regression on the ITT cohort. To do this, unsupervised clustering of cytokine concentrations was used to binarize the overall genital inflammation status into inflammation-high and inflammation-low clusters (Supplementary Fig. 4a, c). The inflammation-high group had significantly higher alpha diversity compared to the inflammation-low group (median SI (IQR): 1.84 (1.00–2.22) vs. 0.83 (0.39–1.55), P < 0.001) and was more likely to be CST-IV (66% vs. 19%, P < 0.001) (Supplementary Fig. 4b). We used HC as the response variable and inflammation (high vs. low) and CST (I vs. III and IV) as our predictor variables (Supplementary Table 5). When comparing adolescents in the COC arm versus those in the Net-En arm at crossover (ITT), CST-III was significantly less prevalent than CST-I in adolescents in the Net-En arm after adjustment (OR: 0.1, 95% CI: 0.02–0.5, P = 0.009), while CST-IV (0.5: 95% CI: 0.1–2.1, P = 0.343) and high inflammation (OR: 1.6, 95% CI: 0.4–6.5, P = 0.488) were not. When comparing adolescents in the COC arm versus the CCVR arm, there were no associations with CST, although CCVR was associated with a nonsignficant increase in inflammation (OR: 3.1, 95% CI: 1.0–10.9, P = 0.072). When comparing adolescents in the CCVR and Net-En arms at crossover, no significant associations were found with inflammation group or CST (Supplementary Table 5).

We next explored the combined effect of HC on bacterial taxa and cervicovaginal cytokines implicated in HIV risk using a sparse Partial Least Squares Discriminant Analysis (sPLSDA). The bacterial taxa included in the analysis were Parvimonas micra, Prevotella spp. (including P. amnii, P. bivia, and P. melaninogenica), G. asaccharolytica, Sneathia, V. montpellieriensis, M. hominis, Megasphaera, L. iners, and L. crispatus due to their prior association with HIV risk or protection9,10. Cytokines IL-1β, IL-6, IFN-γ, and TNF-α were included in the analysis due to their prior association with HIV acquisition19. CCVR inititation was associated with an increase in component 1 (Fig. 4a). P. melaninogenica, IFN-γ, and IL-1β were all positively loaded on this component, while P. amnii was negatively loaded (Fig. 4a, b). Participants in the COC arm had an increase on component 2 between baseline and crossover, and L. iners loaded positively on component 2. On the other hand, Net-En initiation was associated with a positive shift in both components 1 and 2, associated with positive loading of P. bivia and M. hominis on component 1 and negative loading of Sneathia (i.e., Net-En use was associated with an increase in these taxa). However, the shifts between baseline and crossover for Net-En were small, likely because 50% were on Net-En at baseline (Table 1). Cytokine changes were not associated with either COC or Net-En arms in the sPLSDA (Fig. 4b).

Fig. 4: Effect of hormonal contraception on bacterial taxa and cervicovaginal cytokines implicated in HIV risk.
figure 4

a A multilevel sparse partial least-squares discriminant analysis (sPLSDA) of cytokines and bacterial taxa associated with HIV grouped before (baseline) and after (crossover) use of COC, Net-En, or CCVR, respectively. b Barplot showing the loadings of components 1 and 2 from the sPLSDA analysis. The absolute value is an indication of the importance of the bacteria or cytokine, while the sign indicates positive/negative correlations between the variables. CCVR: combined contraceptive vaginal ring, COC: combined oral contraceptives. Source data are provided as a Source Data file.

Discussion

Young African women are at high risk for both HIV infections and unintended pregnancies, and are in need of safe and effective HC. Here, we compare the effects of three HCs on vaginal bacterial communities and inflammation in a randomized trial. Due to the randomized design, we were poised to overcome many of the challenges that have plagued prior observational studies, such as confounding due to reductions in condom use by women using more effective contraceptives. Furthermore, this study focused on adolescent females, who are at the highest risk for both unintended pregnancies and HIV acquisition. We hypothesized that the progestin-only HCs would cause depletion of Lactobacillus in the vaginal microbiota, but that the estrogen component in both CCVR and COC would prevent this from occuring. Further, we hypothesized that Lactobacillus-depleted vaginal microbiota would be associated with increased inflammatory cytokines in cervicovaginal fluid. We found that adolescents assigned to COC not only had significantly less diverse vaginal microbiota compared to those assigned to Net-En, but also those assigned to CCVR. These differences were evident even when adjusting for antibiotic use as well as baseline CST. We found no evidence suggesting that the effect of HC use on vaginal alpha diveristy differed by baseline CST. Although this study was randomized, it was not blinded, and participants may have altered their sexual risk behavior due to perceived risk for pregnancy with use of different and potentially unfamiliar HC methods. Furthermore, use of HC may alter a woman’s menstrual bleeding patterns contributing to the changes oberved in the vaginal microbiota with HC use. However, these potential mediators could only explain a very small and nonsignificant proportion of the effects of study arm on the microbiota. In a longitudinal, paired analysis, alpha diversity was lower when adolescents were on COC compared to both CCVR and Net-En. Community composition differed between assigned study arms at crossover, with COC users more likely belonging to the L. iners-dominated CST-III type, and Net-En and CCVR users more likely to belong to the more diverse CST-IV type. Net-En users displayed a tendency to shift from an L. iners-dominated community toward one of the other CSTs. Taken together, these data suggest that COC may have a beneficial effect on vaginal microbiota composition. Although women randomized to COC in our study had higher relative abundance of L. iners, a non-H2O2-producing lactobacillus whose role in vaginal health is unclear31,32, L. iners abundance was associated with lower vaginal inflammatory cytokine concentrations, and therefore in this African context, may be favorable.

While certain vaginal bacterial species have been associated with an optimal mucosal environment and protection against sexually acquired pathogens, others have been associated with increased risk of STIs9,10. We found that species associated with risk for HIV were more abundant in adolescents randomized to Net-En (such as Prevotella, Sneathia, and Parvimonas) or CCVR (such as Prevotella, Mycoplasma, and Parvimonas) compared to COC. In agreement with our data, a study in Rwandan female sex workers, COC users had lower semiquantitative abundance of Prevotella, Sneathia/Leptotrichia amnionii, and Mycoplasma species33. These results have implications for adolescents in SSA as COCs are rarely used, injectable contraceptives are preferred34,35, and STI prevalence is high36,37.

Most epidemiological data investigating the association between injectable HC and HIV susceptibility have focused on intramuscular DMPA due to the limited use and availability of Net-En6,7. The few studies evaluating Net-En suggest that it may have lower HIV risk implications than DMPA, possibly reflecting different synthetic progestins38. Different progestins could have varied effects on the FGT ecosystem, due to variations in HC dose, steroid receptor binding, peak serum progestin concentrations, and the level of estrogen suppression, among other factors39,40,41. They may therefore also have different effects on HIV susceptibility. The recently completed Evidence for Contraceptive Options and HIV Outcomes (ECHO) Study42, which compared DMPA, copper intrauterine device (IUD), and levonorgesterol (LNG) implant, found no increased risk of HIV with DMPA-IM use compared to Copper-IUD or LNG implant. However, ECHO did not include any of the HC we studied here. Furthermore, ECHO did not include a contraceptive method with an estrogen component and therefore did not address the potential effects of topical or systemic estrogen. The high HIV incidence in contracepting women in the ECHO trial highlighted the urgent need to expand contraceptive method options and integrate HIV and sexual and reproductive health products and services, and therefore evaluation of a broad range of contraceptive options is an urgent public health priority.

There are limited data on the effects of CCVRs on vaginal microbiota composition, with most studies using Nugent scoring and culture-based methods43,44,45. In a recent open-label study, the impact of CCVR use on vaginal microbiota composition was evaluated pre- and 12 weeks post CCVR initiation by qPCR of key vaginal bacteria (i.e., Lactobacillus species, Gardnerella vaginalis, and Atopobium vaginae) in 120 Rwandan women44. The presence and concentrations of Lactobacillus species in vaginal fluid increased significantly, whereas the concentration of G. vaginalis and presence of A. vaginae decreased after initiation of CCVR, suggesting that the CCVR favored healthy vaginal microbe colonization over BV-associated anaerobes44. In our present study, measuring relative abundance rather than absolute quantities of bacteria, we observed the opposite. It is likely that both the relative abundance and the absolute concentrations of key bacteria are important factors for optimal vaginal health. Of more concern was our finding that several cervicovaginal cytokines (including those that could put women at risk for later HIV seroconversion19) were significantly elevated in adolescents randomized to CCVR compared to either COC or Net-En. In a sPLSDA analysis of cytokines and bacterial taxa associated with HIV risk or protection, CCVR use was associated with increases in both cytokines and bacteria associated with risk, while Net-En and COC use was associated with HIV risk-modifying bacteria only (P. bivia, M. hominis, and Sneathia for Net-En and L. iners for COC), possibly explaining the impact of CCVR use on genital inflammation.

Since the vaginal microbiota of adolescents initiating both Net-En and CCVR included higher relative abundance of bacterial species associated with genital inflammation, as well as diverse communities, the lack of inflammation in women intiating Net-En was surprising. Although Net-En users had higher relative abundance of inflammatory taxa, they also had a high relative abundance of L. crispatus, which may dampen any inflammation. Elevated cervicovaginal cytokines in adolescents using CCVRs could reflect an inflammatory cytokine response to the presence of a foreign body in the FGT or to biofilm formation on the ring45. Considering both changes in genital cytokines and microbial diversity associated with HC arm, a multivariate logistic regression analysis suggested that the difference in the diversity of the vaginal microbiota was more profound between study arms than cytokine inflammation categories. The CCVR group also had a higher number of Candida infections compared to the other groups. In Rwandan women, the percentage of women with vaginal yeast increased from 5 to 22% after CCVR initiation46, and in vitro data corroborate this47, further supporting this possibility. Many different vaginal rings are currently in clinical trials, and our findings may be specific to NuvaRing®. Furthermore, our study did not test the effects of continuous-use CCVR, which may be associated with less bleeding and therefore increased quantities of Lactobacillus48.

One limitation of this study was the use of the V4 region for 16S sequencing. The V4 regions did not allow for distinguishing some species-level OTUs, including lactobacilli species such as L. johnsonii, L. gasseri, and L. taiwanensis. Primers spanning the V3/V4 region would identify more species-level resolution in the vaginal microbiota than V4 alone. Our study was also limited by the number of participants enrolled and the loss-to-follow-up rates; however, the longitudinal study design provided us with the opportunity to perform both within-subject assessments comparing the time periods before and after randomization to one of the three HCs in addition to the cross-sectional comparisons of the difference between the study interventions. Due to the study design, participants randomized to COC could have received one of two COC that differed in their components. However, if anything, this would have introduced bias toward to null. Adherence to the randomization arm tended to be poorer for those using the ring, as has been reported in other trials49. However, this study differs from many of these prior ring studies as the adolescents in this study were seeking effective contraception and knew that they had the risk of pregnancy with nonadherence. Additionally, since the aim of the parent study was to assess feasibility and acceptability of methods in adolescents, reporting of problems with adherence and method changes in response were encouraged, albeit reliant on self-report. The study did not include an HC-unexposed group, and the screening visit in most cases did not represent a true, contraceptive-naive baseline control for comparison with no washout periods. Although highly variable, NET can been detected in blood for 6 months after discontinuation of Net-En injection, although at lower levels than needed for contraceptive efficacy50. However, randomization ensured that there was equal distribution between arms of previously used methods of contraception at crossover. Nevertheless, it is not known what impact these low levels of NET have on microbiota and inflammation in the presence of other hormones. Future randomized trials with a larger number of contraceptive-naive adolescents would be helpful to further study the biological impact of HC on FGT microbiota composition. Regardless, insights gained from this study can potentially contribute to the development of expanded, safe contraceptive options for young women in HIV-prevalent settings.

Methods

Study cohort

Adolescents were recruited through a parent study, UChoose, an open-label, randomized crossover study designed to evaluate the feasibility of different HC options among adolescents25. The parent study was approved by the Division of AIDS and the University of Cape Town (UCT) Health Science Research Ethics Committee and was conducted in full compliance with South African Good Clinical Practice (SA-GCP), ICH76 GCP guidelines, and ICMJE guidelines, and registered in the public registry database of ClinicalTrials.gov (NCT02404038). Approval for this substudy was obtained from the Human Research Ethics Committee at the University of Cape Town (HREC 801/2014). Participants were screened, those 18 years or older provided informed consent for the substudy, while assent from the participant and informed consent from a parent or legal guardian were obtained for participants younger than 18 years old. Eligibility criteria for enrollment for the parent study included either HC naive or willingness to change the method, no symptomatic STIs within the prior 40 days, no known sensitivity to any of the study products, and no intentions of becoming pregnant throughout the study period. Furthermore, participants were asked to abstain from inserting any nonstudy products or objects into the vagina throughout the duration of study. The parent study enrolled 130 HIV-negative, nonpregnant adolescent girls aged 15–19 years between September 22, 2015 and June 30, 2017, all of whom consented to participate in this substudy. Participants returned for enrolment within 40 days of the screening visit. At enrolment, participants were randomly assigned in a 1:1:1 ratio to one of three study arms for a 16- week period: Arm 1: injectable hormonal contraception (Net-En, containing 200 mg of the progestogen norethisterone enantate) once every 8 weeks, arm 2: combined contraceptive intravaginal ring (CCVR, NuvaRing®; MSD Pty Ltd, containing etonogestrel/ethinyl estradiol, 0.120 mg/0.015 mg per day) to be inserted once every 28 days (and removed after 21 days of each 28- day insertion), or arm 3: combined oral contraceptive pills (COC) (either Triphasil® or Nordette®, both containing ethinyl estradiol/levonorgestrel (Triphasil® (triphasic regimen): six tablets containing 30 μg/50 μg, five tablets containing 40 μg/75 μg, and ten tablets containing 30 μg/125 μg; Nordette® (monophasic regimen): all tablets contain 30 μg/150 μg daily for 21 days each month with a placebo tablet for 7 days (days 22–28)). After 16 weeks, participants returned at which time they switched over to another HC option for the final 16 weeks of the study. Participants receiving either Net-En or COC for the first 16 weeks of the study (arms 1 and 3) switched over to the CCVR, while those receiving the CCVR as the first method (arm 2) were given the choice between either Net-En or the daily COCs as their second method. After a further 16 weeks, participants returned for a final visit at the clinic and exited the study. Randomization was performed using random number sequence in Stata and provided to the pharmacist in sealed envelopes. Safety was measured by monitoring laboratory and clinical adverse events using the DAIDS table for grading the severity of adult and pediatric adverse events, version 1.0, December 2004 (clarification August 2009).

Power calculations

The sample size of the substudy was limited to that of the parent study, where power calculations were based on the study’s primary outcome, the relative acceptability of CCVR versus other modalities based on the total score for the ORTHO BC SAT questionnaire51 at 4 months after randomization. Power for the substudy was calculated for the aim of evaluating CD4 + cervical T-cell populations at 4 months post contraceptive method initiation52 as prior data on the microbial endpoints were not available. Based on estimates from previous research53,54, n = 50 per group provides 80% power to detect absolute differences of 10% in the mean frequencies of CD4 + HLADR + CD38 + cervical T cells between any two contraceptive groups by one-way ANOVA at alpha = 0.05.

Sample collection

At all study visits, a rapid HIV and a pregnancy test were performed and if positive, the participant was counseled and referred for management, and no further mucosal samples were collected. A detailed interviewer-assisted questionnaire assessing medical history, sexual behavior, the last menstrual cycle, adherence to study the product, intravaginal practices, adverse experiences, and antibiotic use was completed. The following genital tract samples were collected at the screening, crossover, and exit-study visits: cervical secretions collected using a Softcup® menstrual cup inserted for 30 min; two vulvo-vaginal swabs for STI testing, Nugent scoring, Candida screening, and pH measurement, and a lateral wall swab for 16S rRNA gene sequencing. Upon arrival at the laboratory, vaginal swabs and menstrual cup secretions were stored at −80 °C until testing. No samples were collected during menstruation; instead the visit was rescheduled.

STI and BV testing and treatment

Molecular testing for the following STIs: Chlamydia trachomatis, Neisseria gonorrhoeae, Trichomonas vaginalis, and Mycoplasma genitalium by multiplex PCR was performed at each visit as described55. If any of these laboratory-based tests were positive, appropriate targeted therapy was prescribed and recorded. Blood was obtained for HIV rapid test and herpes simplex virus 2 (HSV-2) serology. A vulvo-vaginal swab was collected for BV testing (Gram staining and Nugent scoring; BV negative (Nugent 0–3), intermediate (Nugent 4–6), or positive (Nugent 7–10)), and microscopy for Candida hyphae and spores. Treatment for BV was only offered to participants presenting with clinical symptoms according to South African syndromic management guidelines. Vaginal pH was measured using color-fixed indicator strips (Macherey-Nagel, Düren, Germany).

Sequencing of the V4 region of the 16S rRNA gene

Vaginal lateral wall swabs were collected for microbiome analysis using 16S rRNA gene sequencing. Swabs from 130 participants from 329 sample visits were thawed and treated with an enzyme cocktail consisting of mutanolysin (25kU/ml, Sigma Aldrich), lysozyme (450 kU/ml, Sigma Aldrich), and lysostaphin (4 kU, Sigma Aldrich) for 1 h at 37 °C. Microbial DNA was extracted using the Quick-DNATM Fungal/Bacterial Miniprep kit (Zymo Research) following the manufacturer’s protocol. Mechanical disruption was performed in a Qiagen TissueLyser LT for 5 min at 50 oz. The V4 hypervariable region of the bacterial 16 S rRNA gene was amplified by PCR using modified universal primers56: 515 F (5′-TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG NNN NNG TGC CAG CMG CCG CGG TAA-3′) and 806 R (5′-GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACA GNN NNN GGA CTA CHV GGG TWT CTA AT-3′) (Supplementary Table 6). For each PCR reaction, 1 μl of 515 F primer (5 μM), 1 μl of 806 R primer (5 μM), 12.5 μl, KAPA Hotstart ReadyMix, template DNA (10–20 ng/μl), and molecular-grade water was used for a final volume of 25 μl. PCR was performed in a thermal cycler (Geneamp PCR system 9700, Applied Biosystems) using the following program: 95 °C for 3 min, 35 cycles of 95 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s, and 72 °C for 5 min. Samples were purified using Agencourt AMPure XP beads (Beckman Coulter, Brea, CA, USA) and quantified using the Qubit dsDNA HS Assay (Life Technologies) according to the manufacturer’s protocols. Illumina sequencing adapters and dual-index barcodes were added to the purified amplicon products using limited-cycle PCR and the Nextera XT Index Kit (Illumina) (Supplementary Table 6). For each PCR reaction, 5 μl of forward primer (1 μM), 5 μl of reverse primer (1 μM), 25 μl of KAPA Hotstart ReadyMix, 5 μl of template DNA (amplicon products from the first PCR), and 10 μl of molecular-grade water was used. Limited-cycle PCR was performed using the following program: 95 °C for 3 min, eight cycles of 95 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s, and 72 °C for 5 min. Amplicons from 96 samples and positive (bacterial mock community) and negative (DNA extraction and PCR water) controls were pooled in equimolar amounts, and the resultant libraries purified by gel extraction (Qiagen) and quantified using the Qubit dsDNA HS Assay Kit (Life Technologies). The libraries were sequenced on the Illumina MiSeq platform (300-bp paired end) with v3 chemistry.

Bioinformatics analysis of the 16S rRNA gene-sequencing data

The sequencing reads were analyzed using the following in-house analysis pipeline (https://github.com/uct-cbio/16S-rDNA-pipeline). Following demultiplexing, raw reads were preprocessed as follows: forward and reverse reads were merged using usearch757, allowing a maximum of three mismatches and the resultant merged reads were quality-filtered using usearch7. Primer sequences were removed using a custom python script and reads truncated at 250 bp. Sequences were then dereplicated and clustered de novo into operational taxonomic units (OTUs) at 97% similarity using usearch7. Chimeric sequences were detected (against the Gold database) using UCHIME58 and removed. Individual sequences were assigned to the specific identifiers using a 97% similarity threshold. Taxonomic assignment was performed in QIIME 1.8.059 using the RDP classifier (using the default confidence level of 0.5) against the GreenGenes 13.8 reference taxonomy. To increase species-level resolution, we used the usearch_global command implemented in VSEARCH60 to search the de novo-picked OTUs’ representative sequences against our own Custom Vaginal 16 S Reference Database previously described21,61. All hits with ≥97% identity were accepted. The remaining OTUs (<97% identity) were manually curated using BLAST on NCBI’s nucleotide database (excluding uncultured organisms). OTUs that mapped to more than one species (with the same identity score) were annotated as follows: if an OTU mapped to two or three species, the OTU would be named Genus speciesA_speciesB or Genus speciesA_speciesB_speciesC, respectively; if an OTU mapped to more than three species but one species was clearly associated with vaginal microbiota (based on prior knowledge), the OTU was named Genus species_cluster, where “species” was selected based on the majority of hits21. One OTU that mapped to both L_crispatus and L. acidophilus will be referred to as L. crispatus throughout this paper due to its known presence in the vagina, despite us not being able to distinguish between these two species using the V4 region of the 16 S rRNA gene. The results were imported as a phyloseq object in R62. Samples with ≥5000 reads were selected for downstream analyses. The OTU table was normalized (i.e., transformed to relative abundance * median sample read depth), and filtered so that each OTU had to have at least ten counts in at least 20% of samples or have a relative abundance of at least 0.001% unless otherwise specified.

Cytokine measurements

The concentrations of IL-1β, IL-4, IL-6, IL-10, IL-17A, IL-17F, IL-21, IL-22, IL-23, IL-25, IL-31, IL-33, IFN-γ, sCD40L, and TNF-α in Softcup® cervical secretions were measured by Luminex Milliplex assay using the Bio-Plex Pro Human Th17 Cytokine Panel. The assay plates were read using a Bio-Plex Suspension Array Reader (Bio-Rad Laboratories Inc., USA) and Bio-plex manager software version 4 (Bio-Rad Laboratories, Inc.). A 5-parameter logistic regression formula was used to calculate sample concentrations from the standard curves. For samples with cytokine values below the limit of detection, half of the lowest-detectable value for the specific cytokines was assigned. For samples with values above the detection limit, the highest-detectable value was applied. Specimens from six participants were run across all plates (interplate controls), and samples from six participants were duplicated on each set of plates (intraplate controls) for quality control measures. Spearman’s rank test was used to measure intra-assay and inter-assay correlation coefficients to determine assay reliability and reproducibility. For the individual cytokines, a cutoff of 55% detectable samples was set for inclusion in analyses. Unsupervised clustering of the cytokine data was performed by partitioning around medoids (PAM) clustering available in the R package “cluster”26 with Euclidean distances and optimal k (the k producing the maximum average silhouette width).

Statistical analyses

All downstream statistical analyses were performed in RStudio using the packages phyloseq63 for beta-diversity analyses with weighted UniFrac distances, vegan64 for alpha-diversity estimates (Shannon Index), ordinations, and redundancy analysis. Microbiota community clusters were established by Fuzzy clustering using the R package “cluster”26 with optimal k, a membership exponent of 1.25, and weighted UniFrac as the dissimilarity measure. Differences in study population characteristics according to study arm were tested using Pearson’s Chi-squared test or Fisher’s exact test (when the expected value was <5) for count data and unpaired Student’s t test for differences in mean (parametric data) and unpaired Mann–Whitney U test for differences in medians (nonparametric data) with post hoc testing. Linear regression models were used to adjust for covariates such as antibiotic use and baseline CST. Baseline CST was also included as an interaction term in order to evaluate the impact of HC on the microbiota according to entry microbiota. Mediation analysis was conducted using the R package “mediation” comparing two arms at the time65. Multivariate logistic regression analyses with Tukey’s post hoc testing were used to examine the association of study arm with alpha diversity and inflammation category. Linear mixed- effect regression models for continous data and conditional logit models of binary responses and multinomial counts were used for repeated measure analysis, and Wilcoxon signed-rank and Friedman tests were used for paired, longitudinal analyses. The overall difference in microbial composition between groups was determined by permutational multivariate analysis of variance (PERMANOVA) using distance matrices (weighted UniFrac) with 999 permutations using the adonis2 function; homogeneity of variance between groups was assessed using the betadisper function using vegan64. OTUs were first merged at the lowest available taxonomic level using a custom script21. Spearman’s rank test was applied to test for correlation between nonparametric data. The R function nominalSymmetryTest in the “rcompanion” package66 was used to conduct omnibus symmetry exact tests for the CST transitions between baseline and crossover in each study arm. The model input included the participant ID, CST, and follow-up time in weeks for each sample, and study arm was modeled as a covariate. The study was not powered to calculate the 95% confidence interval to compare the study arms. Differential abundance analysis was performed using DESeq267. OTUs with an adjusted P value of <0.5 and an estimated fold change of >1.5 or <1/1.5 were included. Sparse partial least-squares discriminant analysis (sPLSDA) was used for variable selection using the “mixOmics” package in R68. The input variables used for sPLSDA analysis included log ratio (CLR)-transformed bacterial count data and log-transformed cytokine levels before (baseline) and after (crossover) HC use.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.