The mechanism(s) by which bacterial communities impact susceptibility to infectious diseases, such as HIV, and maintain female genital tract (FGT) health are poorly understood. Evaluation of FGT bacteria has predominantly been limited to studies of species abundance, but not bacterial function. We therefore sought to examine the relationship of bacterial community composition and function with mucosal epithelial barrier health in the context of bacterial vaginosis (BV) using metaproteomic, metagenomic, and in vitro approaches. We found highly diverse bacterial communities dominated by Gardnerella vaginalis associated with host epithelial barrier disruption and enhanced immune activation, and low diversity communities dominated by Lactobacillus species that associated with lower Nugent scores, reduced pH, and expression of host mucosal proteins important for maintaining epithelial integrity. Importantly, proteomic signatures of disrupted epithelial integrity associated with G. vaginalis-dominated communities in the absence of clinical BV diagnosis. Because traditional clinical assessments did not capture this, it likely represents a larger underrepresented phenomenon in populations with high prevalence of G. vaginalis. We finally demonstrated that soluble products derived from G. vaginalis inhibited wound healing, while those derived from L. iners did not, providing insight into functional mechanisms by which FGT bacterial communities affect epithelial barrier integrity.

The female genital tract (FGT) is a key mucosal surface in the context of HIV transmission. Lactobacillus species are beneficial to the FGT, while Garderella vaginalis and other anaerobic bacteria are detrimental. Bacterial vaginosis (BV) is an inflammatory condition characterized by an outgrowth of G. vaginalis and other anaerobes, which is linked to increased HIV acquisition rates. However, the mechanism behind this remains unknown. Here, we used a novel proteomic approach to simultaneously evaluate host and bacterial functions in the FGT. We found that women with G. vaginalis-dominated FGT bacterial communities always displayed markers of decreased epithelial barrier integrity, and decreased wound healing capacity. We also demonstrated that the abundance of proteins from G. vaginalis associated with these signatures of disrupted epithelial integrity. Finally, we showed that products derived from G. vaginalis prevented healing of wounded cell monolayers while products derived from L. iners maintained the ability of the cell monolayers to close wounds. This study provides novel mechanistic insight into the link between BV and increased HIV acquisition rates.

Funding: This work was supported by the Canadian Institutes for Health Research (ADB: TMI-138658), the National Institutes of Health (NIH) (JL: MEHP Al 111738; RMN: NIAID U19 AI078172-01, NICHD HD40539), the Public Health Agency of Canada (GW, ADB), and start-up funds from the University of Washington Department of Pharmaceutics to NRK. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

We sought to better understand the relationship between mucosal health and bacterial diversity using a combination of metaproteomics and metagenomics, which to our knowledge represents the first attempt at integrating these approaches to study the FGT. Indeed the functional diversity of the bacterial proteome, and how this relates to FGT health and inflammation has not been assessed comprehensively, and has largely been limited to 16S rRNA gene sequencing. Thus, we hypothesized that bacterial protein factors can influence FGT mucosal health and affect disease susceptibility. Here we characterized FGT bacterial communities in two distinct human cohorts, longitudinally and cross-sectionally, in asymptomatic and symptomatic women with BV, uncovering bacterial-host interactions leading to wound healing impairment.

Maintenance of the mucosal barrier is critical for preventing invading microorganisms, including HIV, from penetrating into tissues and entering circulation [ 12 ]. Bacterial diversity in the FGT has been strongly associated with negative consequences for FGT mucosa. Highly diverse communities dominated by G. vaginalis and Prevotella are associated with upregulated expression of Toll-like Receptor (TLR) and NFkB pathways, leading to increased pro-inflammatory cytokine concentrations and activation of immune cells [ 13 ]. While it is widely appreciated that BV is associated with inflammation, the mechanism that elicits this inflammation or the bacterial proteins associated with inflammation remain unresolved [ 14 ], which may partly explain the limited effectiveness of antimicrobial treatment for BV [ 15 – 17 ]. Bacterial metabolites including hydrogen peroxide, antimicrobial peptides, and acids that reduce the FGT pH have been proposed to have an important impact in sustaining mucosal health [ 3 ]. Furthermore, the integrity of mucosal epithelial surfaces has been shown to depend on bacterial community composition in other diseases [ 18 ], and has been proposed to be important in the FGT during bacterial dysbiosis [ 19 ], but this has not been extensively studied. Each of these factors likely impact disease susceptibility independently, and a Lactobacillus-dominant microbiota likely contributes to many of these factors to maintain the function of the healthy FGT and inhibit infections. Taken together these studies suggest that host-microbe interactions are key to understanding negative consequences on vaginal health, yet this interaction remains poorly defined in human cohorts [ 20 ].

Mucosal surfaces exposed to the external environment contain distinct bacterial communities that exist in relationship with the host and can contribute to health and functioning. These bacterial communities have been linked to several human diseases and overall health [ 1 ], and can vary between individuals, but also over time within the same person [ 2 ]. In the female genital tract (FGT), colonization by Lactobacillus species and other lactate-producing bacteria helps to inhibit colonization by pathogenic bacteria [ 3 ]. However, colonization by more diverse communities of anaerobic bacteria, notably Gardnerella vaginalis, is common [ 4 ], and often associated with the development of bacterial vaginosis (BV) [ 5 ]. BV is highly prevalent, affecting 4–58% of women globally; some areas, such as sub-Saharan Africa have rates as high as 55% [ 6 ]. BV is associated with significant health consequences, including pre-term birth, post-partum endometriosis, pelvic inflammatory disease, upper reproductive tract infections, and increased susceptibility to sexually transmitted infections (STI’s) [ 7 , 8 ], with HIV being highly significant [ 9 , 10 ]. Indeed, a recent meta-analysis linked BV to a 60% increase in HIV acquisition rates [ 11 ]. However, while these relationships between microbial composition and vaginal health have been described epidemiologically, there is limited understanding about the mechanisms underlying the impact of bacterial dysbiosis on the vaginal mucosa.

The association of bacterial communities with barrier integrity proteins led us to hypothesize that wound-healing capacity may be supported or inhibited by specific bacterial species and/or their products. We thus performed a classical wound-healing assay wherein we cultivated relevant cervical cell line (HeLa CCL-2) in the presence of supernatants derived from cultures of L. iners or G. vaginalis. Prior to adding culture supernatants, a wound was induced by scratching HeLa cell monolayers. Incubation of scratched monolayers with L. iners culture supernatant did not alter wound healing compared to the control incubations. However, incubation with G. vaginalis culture supernatants significantly reduced wound healing after 24 hours compared to both the control and L. iners conditions ( Fig 5 ). These results confirm a relationship between soluble compounds produced by the major bacterial species of the G1 and G2 profiles and wound healing capacity. This implicates these species as important components or drivers of epithelial barrier repair, maintenance, and disruption in the FGT.

Due to the strong association of epithelial development pathways with different bacterial groups, we compared cornified envelope factors INVO and SPR1A to bacterial proteins. Nineteen bacterial proteins had strong associations in at least one comparison against either INVO or SPR1A after correcting for multiple comparisons. Proteins from L. iners that positively correlated with INVO and SPR1A were involved in Catabolism and Energy Metabolism pathways, including glycolysis and homolactic fermentation of sugars (Embden-Meyerhoff-Parnas (EMP) pathway) ( Fig 4A/4B ). These included a putative fructose 1,6-bisphosphate aldolase (PFBA), PK, GAP-DH, and LDH, as well as a ferritin-like protein (FLP), which is important for sequestering excess iron and preventing oxidative damage [ 26 ]. Bacterial proteins that negatively correlated with INVO and SPR1A belonged to alternate sugar metabolism pathways (phosphoketolase pathway), transport functions, and amino acid catabolism. The majority of these belonged to G. vaginalis ( Fig 4C/4D ), including D-xylulose 5-phosphate/D-fructose 6-phosphate phosphoketolase (XFBP), a putative sugar-binding secreted protein (P-SBSP), MAL-E ABC, an extracellular solute binding protein (ESBP), and glycine oxidase (GOx). A membrane protein from Prevotella sp., sharing sequence homology with SusD-like (starch-binding) protein, was also negatively associated ( S4 Fig ). Many associations with vaginal pH were also observed, including negative associations with enzymes from L. iners (PK, LDH, elongation factor tu, and FLP) ( S5A Fig ), and positive associations with enzymes from G. vaginalis, including GOx, which catalyzes the conversion of glycine into glyoxalate, ammonia, and hydrogen peroxide ( S5B Fig ). Therefore a clear relationship between metabolic function, epithelial barrier protein levels, and vaginal pH was observed, demonstrating that these microbial pathways may be an important component of mucosal barrier disruption and vaginal health.

(A) Heatmap of differentially abundant host proteins found in common between Cohort 1 (longitudinal) and Cohort 2 (cross-sectional) based on cervicovaginal microbiome profiles. (B) Functional annotation of differentially expressed proteins based on their gene ontology. Fisher’s exact test was used to compare annotated protein functions between G1 and G2 across both cohorts. (C) Lowered cornified envelope proteins (Involucrin and Small Proline-Rich Protein 1A) strongly correlate with G. vaginalis (G2) profiles, regardless of the presentation of clinical BV symptoms for both cohorts. Student’s t-test was used to compare the abundances of INVO and SPR1A between G1 and G2. BP: Biological process; CC: Cellular component.

Seventeen proteins were differentially abundant across both cohorts ( Fig 3A ). Host proteins more abundant in G2 profiles included apoptotic regulators (PRDX, NDKB, CADH1) and leukocyte migration factors (PLST), while G1 profiles showed increased keratinization, epidermis development, and cornified envelope (INVO, SPR1A) factors ( Fig 3B ). Of particular interest, the abundances of INVO and SPR1A were 14.7 and 7.2-fold lower in women G2 microbial profiles Cohorts 1 and 2, respectively. In Cohort 2, INVO and SPR1A were lower for women with G2 microbial profiles even if they had not been clinically diagnosed for BV ( Fig 3C ). These proteins are known to act as scaffolding for epidermal layers and are important for proper barrier function [ 25 ], and immunohistochemical analysis confirmed the presence of INVO and SPR1A in cervical and vaginal tissues, where they strongly associated with the squamous epithelium and stratum corneum in healthy FGT tissue ( S3 Fig ). Collectively these data show an association of heightened immune activation, apoptosis, and decreased epithelial barrier function in women with G. vaginalis-dominated bacterial profiles and that these effects are evident in G. vaginalis-dominated communities in the absence of clinical diagnosis.

Bacterial dysbiosis impacts HIV acquisition risk [ 10 , 11 , 24 ], reproductive health [ 7 ], and mucosal cellular activation [ 13 ], but the effect on the FGT is not well defined. Our analysis revealed that 69/434 (15.8%, 15 passing 5% FDR) and 64/434 (14.7%, 19 passing 5% FDR) host proteins were significantly differentially abundant between G1 and G2 profiles in Cohort 1 and 2, respectively. For Cohort 2, comparison based on bacterial groups rather than Nugent score criteria yielded greater host proteome differences, statistically (9.2% vs. 15.8%, P<0.05), and in magnitude (5 vs. 6 Log 2 Fold Change; S1G/S1H Fig ), suggesting that bacterial community composition, rather than clinical BV criteria, more accurately classifies mucosal inflammation. This comparison was not possible for Cohort 1, as all G2 profiles had clinically defined BV. Hierarchical cluster analysis revealed that longitudinal changes from G1 to G2 profiles in Cohort 1 were clearly distinguishable by two major branches of host proteins ( S2A/S2B Fig ). Proteins more abundant in G1 (Branch 1) associated with epidermis development and the cornified envelope, whereas G2 (Branch 2) showed increased factors involved in cytoskeletal-binding, threonine proteases involved in proteasome activity, as well as vesicular components and the melanosome. Many of these included S100 proteins and innate immune factors, important for antimicrobial defense based on gene ontology (DMBT1, CADH1, S10A7, EFHD2, S10AB, S10A6, TGM3, K2C1 S10A2). Similarly, in Cohort 2, hierarchical cluster analysis showed that proteins more abundant in G1 (Branch 1) associated with epidermis development, structural molecular activity, and the cornified envelope, while proteins elevated in G2 (Branch 2) also included ectoderm development and differentiation, although were related to cytoskeletal activity ( S2C/S2D Fig ). Many of these are important for leukocyte-mediated immunity and wounding responses based on their gene ontology (A1AT, IC1, GELS, CO3, PEBP1, PRDX1, PRDX2, CO4A, ANXA8).

However, unique functional signatures were observed between G1 and G2 FGT bacterial communities. Across both cohorts, the G1 group showed significant enrichment of proteins involved in transport and catabolism (6.9% higher), energy metabolism (5.6% higher), and folding, sorting, and degradation (4.9% higher), while G2 was highly significantly enriched in membrane transport functions (22% higher). Twelve bacterial proteins were significantly differentially abundant after multiple comparison correction between G1 and G2 in Cohort 1 ( S2 Table ). Proteins enriched in G1 mostly belonged to L. iners proteins and were involved in homolactic fermentation of carbohydrates including glyceraldehyde-3-phosphate dehydrogenase (GAP-DH), pyruvate kinase (PK), and lactate dehydrogenase (LDH). Proteins enriched in G2 were all G. vaginalis proteins and included a MalE-type ABC sugar transport system periplasmic component (MAL-E ABC) and an alpha-1,4 glucan phosphorylase, an enzyme that degrades starch and glycogen, suggesting that G. vaginalis directs its metabolism towards liberation and uptake of extracellular saccharides. Although many of these proteins were also differentially abundant between G1 and G2 in Cohort 2 they did not pass multiple comparison correction. Overall this shows that ‘core’ functional pathways necessary to host-associated bacterial life within the FGT include carbohydrate, amino acid, and translational machinery, while perturbations to membrane transport and carbohydrate catabolism are likely important for pathogenic states.

Bacterial proteome functions for both cohorts were determined by KEGG ontology. For both cohorts, G1 profiles (Lactobacillus dominant) showed increased abundance of proteins related to Transport and Catabolism and Energy Metabolism while G2 (Gardnerella dominant) showed increased abundances of proteins associated with Membrane Transport. Student’s t-tests were used to compare the abundances of each category between G1 and G2.

As bacterial diversity has been associated with other biological variables, such as concurrent STI’s [ 23 ] and hormonal contraceptive usage [ 10 ], we compared clinical characteristics between Lactobacillus and Gardnerella-dominant groups (Cohort 1- Table 1 ; Cohort 2- Table 2 ). With the exception of BV status, we found no differences between G1 and G2 with respect to age, contraceptive use, antimicrobial usage, last menstrual period, detectible STI’s, or sexual practices in either cohort. There were differences in Amsel’s criteria collected from Cohort 2 ( S1 Table ) between G1 and G2, where vaginal pH, clue and white blood cell presence was higher in women with G2 bacterial proteome profiles, in agreement with clinical BV status. Overall, there was no evidence to support that vaginal bacterial profiles were related to exogenous hormonal contraceptive use, the menstrual cycle, sexual behaviors, or concurrent STI’s in these cohorts.

Mucosal samples from a separate group of 31 women from North America (Cohort 2) were analyzed to further evaluate associations between FGT bacterial proteome diversity and BV. MS analysis of Cohort 2 samples showed similar trends to that of Cohort 1 ( Fig 1C , species-level taxonomy shown in S1B Fig ), including Lactobacillus spp.-dominant (G1) and G. vaginalis-dominant (G2) communities. A wider distribution of lactobacilli including L. iners, L. crispatus, and L. jensenii was observed in G1 in Cohort 2 than Cohort 1. Varying abundances of other BV-associated bacteria including Prevotella spp., Atopobium vaginae, Mobiliuncus mulieris, and Sneathia sp. were also observed. In Cohort 2, there was no difference in the species diversity of G1 compared to G2 ( S1D Fig ). While all 7 women with BV clustered into G2, 46% of participants demonstrated a G. vaginalis-dominated proteome despite a lack of clinical BV diagnosis, consistent with the observation that not all women have Lactobacillus-dominant FGT microflora despite low Nugent scores. Also similar to Cohort 1, G2 in Cohort 2 had higher overall microbial proteome burden than G1 (1.5-fold higher; S1F Fig ), indicating further changes in bacterial community function. This agrees with other studies showing Lactobacillus dominance varies between 37–90% of women, with greater diversity and variation in African women [ 5 , 13 , 22 ].

Bacterial community composition for G1 and G2 was confirmed by 16S rRNA gene sequencing ( Fig 1B ). According to 16S rRNA gene sequencing, G1 communities were dominated by Lactobacillus spp. (26–99% of the total community), and G2 communities were dominated by Gardnerella, but at lower proportions than were detected by MS (17–66% of the total community). As with MS, bacterial genera detected in BV-positive individuals by 16S rRNA gene sequencing included Sneathia, Prevotella, Atopobium, Megasphera, and others. 16S rRNA gene sequencing also detected greater bacterial diversity in the G2 samples compared to G1 ( S1C Fig ). Several species detected by 16S and not by MS included Leptotrichia, Fastidiosipila, Shuttleworthia, and Aerococcus. Overall, this demonstrates significant heterogeneity in the structure of FGT bacterial communities between clinically defined BV and asymptomatic time points, that G. vaginalis and other anaerobes associate with BV, and that specific species dominate the bacterial proteome landscape in mucosal secretions.

(A) Bacterial genera distributions detected in Cohort 1 using MS. Bacterial abundances from the MS data were calculated by summing normalized total spectral counts for all proteins associated with a genus/species. Clustering of samples was performed using unsupervised hierarchical linkage with average Euclidean distances of the proportional bacterial abundance in each sample. (B) Bacterial genera distributions for Cohort 1 detected using 16S rRNA gene sequencing shows close concordance to MS, and clustered into two groups based on Lactobacillus (G1) or Gardnerella (G2) dominance. Bacterial genera detected in Cohort 2 using MS (C) showed similar grouping as Cohort 1. Nugent scores (top) and bacterial alpha diversity (Shannon’s H index, bottom) are also shown.

Cervicovaginal secretion samples from two cohorts of women were evaluated to understand the mucosal environment associated with bacterial dysbiosis. We first assessed mucosal changes in women at BV+ or BV- time points (Cohort 1, n = 10), through a combination of mass spectrometry (MS) and 16S rRNA gene sequencing. MS analysis identified 1123 unique proteins, including 434 human and 689 bacterial proteins from 64 species. To assess the diversity of the bacterial proteome, we quantified the relative proteome load of each bacterial genus in each sample by summing the total number of protein spectral counts assigned to each genus, an approach previously shown to directly correlate with colony-forming units [ 21 ]. We clustered the bacterial proteomes from the twenty samples using unsupervised hierarchical clustering. Two major bacterial proteomes were identified, dominated by either Lactobacillus iners (Group 1, or G1) or Gardnerella vaginalis (Group 2, or G2) ( Fig 1A , species-level taxonomy shown in S1A Fig ), which were used for downstream comparisons. In G1, L. iners proteins accounted for 87–100% of the total protein load while in G2, G. vaginalis proteins accounted for 48–96%. Compared to those in G1, the bacterial proteomes in G2 displayed significantly higher species diversity ( S1C Fig ). G2 profiles also had higher overall bacterial protein load when normalized to the total protein content (0.34 log 10 , +2.2 fold higher; S1E Fig ). L. iners dominated the FGT bacterial proteome of eight of the 10 patients from Cohort 1 at the time point without clinically diagnosed BV, but not the remaining two. Patient “10“, at the time point without BV, displayed a high abundance of G. vaginalis with a lower abundance of L. iners, and Patient “6”at the time point without BV had high levels of Lactococcus lactis and Streptococcus mitis. In contrast, all samples taken during episodes of BV had high abundances of proteins from G. vaginalis, Prevotella spp., Streptococcus mitis, Escherichia coli, and Atopobium vaginae, which have been previously identified as part of BV-associated bacterial communities [ 7 ].

Discussion

In this study, we demonstrated a novel metaproteomic approach to simultaneously assess bacterial diversity, abundance, and function, along with host barrier and inflammation processes, providing mechanistic insight relevant to women’s health. We described distinct vaginal bacterial proteome profiles that were dominated by Lactobacillus spp. (G1) or G. vaginalis (G2), where the latter associated with BV, increased community diversity, and significant divergence from normal metabolic function. We next demonstrated that bacterial functional profiles were significantly associated with cornified envelope factors in the FGT, and this was affected even in the absence of clinical diagnosis. Finally, we found that predominant species identified in this study, specifically G. vaginalis and L. iners, generate soluble products that disrupt or maintain the ability of cervical epithelial cells to repair and close wounds. Therefore, impaired wound healing is a potential mechanism by which key bacterial species may impact mucosal barrier function and therefore disease and/or HIV/STI infection risk. The association of vaginal inflammation and inflammatory vaginal bacteria with HIV susceptibility indicates that targeting this mechanism may lead to novel prevention strategies for HIV.

While increased diversity of bacterial communities has been linked to better mucosal functioning in the gut [27, 28], low-diversity bacterial communities are beneficial for the FGT [5], where increased diversity is strongly associated with BV [7]. Consistent with previous observations [5], many of the women in Cohort 2 with G. vaginalis-dominated communities were asymptomatic for BV (61%), further supporting the fact that Nugent score is underestimating the extent of non-Lactobacillus dominant communities. However, the effects on host epithelial pathways, including decreased integrity and increased inflammatory pathways were still evident in the absence of clinical diagnosis. Lactobacillus spp. and G. vaginalis proteins comprised more of the soluble proteome load than might be inferred from 16S rRNA gene sequencing, suggesting that these bacteria dominate the metabolic landscape of the FGT.

Metagenomic studies of the human microbiome have shown that core metabolic function is less variable than the community composition [2]. In agreement with this, we observed that the majority of assigned protein functions did not vary significantly, which likely represent core metabolic functions. However, some functions varied between G1 and G2, including increased carbohydrate metabolism, energy production, and folding/sorting functions in G1 to enhanced membrane transport and secretion of extracellular products in G2, with L. iners and G. vaginalis dominating these key functions. The increased abundances of enzymes important for sugar transport and starch and glycogen catabolism in G2 suggest that G. vaginalis may outcompete Lactobacillus spp. for the uptake of carbohydrate substrates. This agrees with a recent study showing that women with BV have significant metabolite alterations in cervicovaginal mucous, including lower levels of carbohydrates, amino acids, and lactate, accompanied by increased levels of amino acid catabolites and polyamines [29]. Overall, this demonstrates that increased bacterial diversity is associated with changes in key metabolic pathways, which allows for better understanding of dysbiosis in the FGT.

We found that G1 profiles from both cohorts strongly associated with cornified envelope factors, especially INVO and SPR1A, which are expressed in the upper layers of the vaginal and cervical epithelia, and aid in maintaining epithelial integrity. Our group has previously reported that increased levels of cervocovaginal CD4+ T cells associated with lower levels of cornified envelope factors [30], demonstrating the important link between vaginal epithelial integrity and HIV acquisition risk. G1 profiles were also associated with higher levels of antimicrobial peptides, such as dermcidin, which is important for host defense against microorganisms [31]. In comparison, the G2 bacterial profiles correlated with lower cornified envelope and epithelial barrier factors, increased cytoskeletal elements important for cell migration, and increased proteasome factors. This agrees with other studies showing that BV associates with activation of innate immune and inflammation pathways in the FGT, including increased complement [32], proteasome levels [33], and pro-inflammatory cytokines and activated CD4+ T-cells [13]. Importantly, G2 bacterial proteome profiles associated with decreased abundances of INVO and SPR1A regardless of clinical BV status. This finding demonstrates that current methods used to diagnose BV likely underestimate the true extent of bacterial dysbiosis on mucosal barrier function the FGT, as the Nugent Scores were poor predictors of BV, especially for Cohort 2. Thus, new methods to detect and treat G. vaginalis in the FGT could aid in reducing HIV acquisition risk by promoting mucosal and epithelial barrier integrity, and reduced inflammation.

Catabolic enzymes involved in homolactic fermentation of glucose from Lactobacillus, such as L-lactate dehydrogenase, correlated with higher epithelial barrier proteins, while membrane transporters, extracellular proteins, and alternate routes of carbohydrate metabolism (heterolactic fermentative or phosphoketolase pathways) from G. vaginalis were negatively correlated. In addition, GOx, was strongly correlated with increased vaginal pH, implicating a role of this enzyme in altered vaginal pH during dysbiosis. To our knowledge this is the first time these bacterial enzymes have been associated with epithelial disruption signatures and vaginal pH. Collectively, this shows a relationship between bacterial community structure, metabolic function, disruption of epithelial proteins important for barrier integrity, and overall vaginal health.

We also demonstrated that G. vaginalis culture supernatants inhibited healing of scratched HeLa cell monolayers while, L. iners culture supernatants maintained effective wound healing. Based on these data, G. vaginalis is likely an important component or a potential driver of subverting the wound healing process. While acknowledging that HeLa cell monolayers do not completely recapitulate the squamous epithelium or immune environment of the FGT, this nevertheless supports would healing as an underlying mechanism. Taken collectively, and considering the metaproteomic, metagenomic, and in vitro models, these data suggest that G. vaginalis releases a variety of extracellular products in the vaginal compartment that aid in uptake for nutrients, alter the vaginal microenvironment, contribute to innate immune activation, and prevent healing of the epithelial barrier. Future studies to identify exact protein pathways involved, how they may be altered, and more advance animal and engineered tissue models would help better decipher these host-microbiome interactions.

It is important to compare discuss the benefits and limitations of metaproteomics compared to 16S rRNA-based techniques to characterize microbial communities in the vaginal compartment. Both techniques are quantitative and spectral counts by MS have been shown to correlate directly to colony-forming units [21]. An advantage of 16S over MS is greater resolution of the overall community structure, and while we showed high sensitivity to identify species that were at 0.1% of the population by MS, 16S captured more overall bacterial species. It is likely that the larger dynamic range of the proteome over the genome is a large contributing factor to this observation. Both 16S and MS rely on curated databases to identify species and are subjected to this same limitation in availability and extensiveness of libraries. While databases for 16S rRNA genes are likely more comprehensive, proteomic libraries are growing and becoming more available. MS is advantageous in that it can provide direct species-level identification, which is not achievable through high-throughput 16S rRNA gene sequencing methods. Furthermore, metaproteomic analysis reveals bacterial functional and metabolic activity, which is not provided by 16S-based approaches. Prior studies have attempted to alleviate this using MS to correlate metabolite abundances with species abundances [34], through metagenomic studies [35], or by employing computational methods to estimate bacterial community functional capacity based on 16S rRNA gene signatures [36], but nevertheless represent indirect methods to evaluate bacterial community functionality. While 16S rRNA gene sequencing is a popular and well-validated method for studying microbial communities, the use of metaproteomic approaches provides complimentary and invaluable data on community structure, function, and host inflammation to better study host-bacterial relationships.

Our data provide novel mechanistic insight of how dysbiosis of vaginal bacterial communities may directly increase host susceptibility to infection through the disruption of epithelial barriers, inhibition of wound repair, and induction of inflammation. In the context of HIV transmission, inhibition of wound repair is under studied and may represent underlying mechanisms in other risk factors for HIV, including hormonal contraceptive usage, intravaginal practices, and other STI’s. These pathways may also impact the effectiveness or responsiveness to mucosa-targeted prevention technologies for other infections, such as microbicides or vaccines for HIV. In summary, this study delineated functional configurations of microbial communities that impact vaginal health during BV, providing new information on host-bacterial interactions, enabling future experiments to probe host-microbe relationships in the FGT that could have important implications for women’s health.