Abstract
The vaginal microbiome plays an influential role in several disease states in reproductive age women, including bacterial vaginosis (BV). While demographic characteristics are associated with differences in vaginal microbiome community structure, little is known about the influence of sexual and hygiene habits. Furthermore, associations between the vaginal microbiome and risk symptoms of bacterial vaginosis have not been fully elucidated. Using Bayesian network (BN) analysis of 16S rRNA gene sequence results, demographic and extensive questionnaire data, we describe both novel and previously documented associations between habits of women and their vaginal microbiome. The BN analysis approach shows promise in uncovering complex associations between disparate data types. Our findings based on this approach support published associations between specific microbiome members (e.g., Eggerthella, Gardnerella, Dialister, Sneathia and Ruminococcaceae), the Nugent score (a BV diagnostic) and vaginal pH (a risk symptom of BV). Additionally, we found that several microbiome members were directly connected to other risk symptoms of BV (such as vaginal discharge, odor, itch, irritation, and yeast infection) including L. jensenii, Corynebacteria, and Proteobacteria. No direct connections were found between the Nugent Score and risk symptoms of BV other than pH, indicating that the Nugent Score may not be the most useful criteria for assessment of clinical BV. We also found that demographics (i.e., age, ethnicity, previous pregnancy) were associated with the presence/absence of specific vaginal microbes. The resulting BN revealed several as-yet undocumented associations between birth control usage, menstrual hygiene practices and specific microbiome members. Many of these complex relationships were not identified using common analytical methods, i.e., ordination and PERMANOVA. While these associations require confirmatory follow-up study, our findings strongly suggest that future studies of the vaginal microbiome and vaginal pathologies should include detailed surveys of participants' sanitary, sexual and birth control habits, as these can act as confounders in the relationship between the microbiome and disease. Although the BN approach is powerful in revealing complex associations within multidimensional datasets, the need in some cases to discretize the data for use in BN analysis can result in loss of information. Future research is required to alleviate such limitations in constructing BN networks. Large sample sizes are also required in order to allow for the incorporation of a large number of variables (nodes) into the BN, particularly when studying associations between metadata and the microbiome. We believe that this approach is of great value, complementing other methods, to further our understanding of complex associations characteristic of microbiome research.