Modeling the temporal dynamics of cervicovaginal microbiota identifies targets that may promote reproductive health

Cervicovaginal bacterial communities composed of diverse anaerobes with low Lactobacillus abundance are associated with poor reproductive outcomes such as preterm birth, infertility, cervicitis, and risk of sexually transmitted infections (STIs), including human immunodeficiency virus (HIV). Women in sub-Saharan Africa have a higher prevalence of these high-risk bacterial communities when compared to Western populations. However, the transition of cervicovaginal communities between high- and low-risk community states over time is not well described in African populations. We profiled the bacterial composition of 316 cervicovaginal swabs collected at 3-month intervals from 88 healthy young Black South African women with a median follow-up of 9 months per participant and developed a Markov-based model of transition dynamics that accurately predicted bacterial composition within a broader cross-sectional cohort. We found that Lactobacillus iners-dominant, but not Lactobacillus crispatus-dominant, communities have a high probability of transitioning to high-risk states. Simulating clinical interventions by manipulating the underlying transition probabilities, our model predicts that the population prevalence of low-risk microbial communities could most effectively be increased by manipulating the movement between L. iners- and L. crispatus-dominant communities. The Markov model we present here indicates that L. iners-dominant communities have a high probability of transitioning to higher-risk states. We additionally identify transitions to target to increase the prevalence of L. crispatus-dominant communities. These findings may help guide future intervention strategies targeted at reducing bacteria-associated adverse reproductive outcomes among women living in sub-Saharan Africa.