Skip to main content

Lactic acid from vaginal microbiota enhances cervicovaginal epithelial barrier integrity by promoting tight junction protein expression

Abstract

Background

Women with a cervicovaginal microbiota dominated by Lactobacillus spp. are at reduced risk of acquiring sexually transmitted infections including HIV, but the biological mechanisms involved remain poorly defined. Here, we performed metaproteomics on vaginal swab samples from young South African women (n = 113) and transcriptomics analysis of cervicovaginal epithelial cell cultures to examine the ability of lactic acid, a metabolite produced by cervicovaginal lactobacilli, to modulate genital epithelial barrier function.

Results

Compared to women with Lactobacillus-depleted microbiota, women dominated by vaginal lactobacilli exhibit higher abundance of bacterial lactate dehydrogenase, a key enzyme responsible for lactic acid production, which is independently associated with an increased abundance of epithelial barrier proteins. Physiological concentrations of lactic acid enhance epithelial cell culture barrier integrity and increase intercellular junctional molecule expression.

Conclusions

These findings reveal a novel ability of vaginal lactic acid to enhance genital epithelial barrier integrity that may help prevent invasion by sexually transmitted pathogens.

Video abstract

Background

The cervicovaginal microbiota has a substantial influence on sexual and reproductive health in women including their susceptibility to sexually transmitted infections (STIs) [1]. The composition of the cervicovaginal microbiota varies between women of different age, ethnicity and geographical location and is influenced by menstrual cycle phase and behavioural factors [2,3,4]. A cervicovaginal microbiota dominated by noninflammatory Lactobacillus spp. is considered optimal [5] and is associated with favourable health outcomes including lower rates of preterm birth [6,7,8] and reduced susceptibility to STIs including human immunodeficiency virus (HIV) [1, 9,10,11]. Women with a cervicovaginal microbiota dominated by L. crispatus have a 4.4-fold reduced risk of HIV acquisition compared to women with a nonoptimal microbiota lacking lactobacilli [10]. In contrast, the most common form of a nonoptimal microbiota, bacterial vaginosis (BV), characterised by an increase in anaerobic bacteria and depletion of lactobacilli, is associated with adverse pregnancy outcomes including low birth weight [12], spontaneous preterm birth [7] and increased risk of acquiring and transmitting STIs including HIV [1, 13,14,15,16].

The mucosal surface of the female reproductive tract (FRT) is the primary site of entry for many STIs including HIV and, under optimal conditions, possesses substantial physical and immunological defences against pathogen invasion [17]. Interactions between FRT epithelial cells are mediated by intercellular junctional molecules comprising tight junctions, which control paracellular permeability, along with adherens junctions and desmosomal proteins that represent a physical barrier to the passage of microorganisms [18]. This barrier plays an integral role in precluding HIV transmission by preventing virions present in the cervicovaginal lumen from accessing target cells typically located in the subepithelial lamina propria [19, 20].

Invading microorganisms can penetrate the columnar (endocervix and endometrium) and the stratified (ectocervix and vagina) epithelial cell layers through different mechanisms. The FRT epithelial barrier can be compromised by physical damage including microabrasions and/or disruption of the intercellular junctional barrier where the pathogen infiltrates between the intercellular space by paracellular penetration [20, 21]. Alternatively, pathogens can penetrate the epithelium by transcytosis where the apically endocytosed pathogen is transferred via the basal epithelial surface by cell-to-cell spread to susceptible target cells, as reported for HIV in vitro [22]. Paracellular penetration is considered a major mechanism for invasion by pathogens including HIV and has been observed in human cervical explants and following in vivo vaginal exposure of non-human primates [20].

Many microorganisms have evolved strategies to target junctional molecules to aid penetration. Viruses including human papillomavirus [23] and HIV [24, 25] directly modulate the integrity of the FRT epithelial barrier to increase pathogen transmigration and facilitate infection. The increased risk of STIs associated with BV may also be due in part to effects on epithelial barrier integrity, as cervicovaginal fluid from women with BV disrupts the endocervical epithelial barrier and increases transmission of HIV-infected lymphocytes in vitro [26].

Whilst there is an established relationship between an optimal Lactobacillus-dominant microbiota and favourable sexual and reproductive health outcomes, the mechanisms through which this occurs are still being uncovered. In addition to precluding replication of STIs [27,28,29], lactobacilli and their products may act directly on the FRT mucosa to promote beneficial properties including enhanced wound healing [30], epithelial barrier integrity [31] and a noninflammatory state [32, 33]. Whether metabolites produced by vaginal lactobacilli directly enhance the FRT epithelial barrier is unknown.

Lactobacilli produce the metabolite lactic acid (LA) at levels that can vary depending on the species and strain [34, 35]. In women with a microbiota dominated by optimal Lactobacillus spp., LA concentrations are approximately 110 mM (1% w/v) which acidifies the vagina to a pH of ~3.5 [4, 36]. LA has direct antiviral and antimicrobial activity against BV-associated bacteria [37], HIV [38, 39] and other STIs including herpes simplex virus 2 [40] and Neisseria gonorrhoeae [27]. We have shown that physiologically relevant concentrations of LA have an anti-inflammatory effect on FRT epithelial cells and inhibit pro-inflammatory cytokine and chemokine production elicited by viral and bacterial pathogen-associated molecular patterns [41, 42]. LA can exist in both L- and D-isomers, with different Lactobacillus spp. preferentially producing each isoform [34, 43]. The protonated form of LA predominates at a pH below its acid dissociation constant (pKa 3.86) and possesses virucidal and anti-inflammatory activity, in contrast to the lactate anion that predominates at pH > 3.86 [38, 41, 42]. These data, and others [29, 44, 45], indicate that metabolites produced by optimal FRT microbiota species are likely key mediators of the health benefits associated with Lactobacillus-dominant microbiota.

There is limited evidence that cervicovaginal bacteria and their metabolites can enhance FRT epithelial integrity. Zevin et al. identified an association between Lactobacillus-dominant FRT microbiota and expression of host proteins relevant for maintaining epithelial barrier function, observing a positive association observed between levels of bacterial L-lactate dehydrogenase (L-LDH, a key enzyme involved in LA production by lactobacilli) and epithelial barrier proteins [46]. We considered whether these observations indicated a role for vaginal LA in directly enhancing the FRT epithelial barrier and whether this may be mediated by modulation of intercellular junctional molecules. In this study, we explored the in vivo relationship between LA and epithelial barrier integrity in a large cohort of young women in South Africa with diverse FRT microbiota and used functional and transcriptomic in vitro analyses to demonstrate the ability of LA to directly enhance FRT epithelial barrier integrity.

Methods

Metaproteomics analysis of cervicovaginal swab samples

Analysis of the vaginal microbiome and secretome was performed on cervicovaginal fluid samples collected from young women in the Women’s Initiative in Sexual Health (WISH) study conducted in South Africa. A full description of the study design and cohort have been detailed previously [47] but included sexually active, HIV-negative, nonpregnant women aged 16–22 years. Lateral vaginal wall swabs were collected by the study nurse by rotating Dacron swabs 360° on the lateral vaginal wall before being placed in 1 ml phosphate-buffered saline (PBS), transported to the laboratory at 4° C and stored at −80 °C. Samples were computational randomised and the eluted material subjected to metaproteomic analysis using liquid chromatography-tandem mass spectrometry (LC-MS/MS) on a Q-Exactive quadrupole-Orbitrap MS (Thermo Fisher Scientific, MA, USA) coupled with a Dionex UltiMate 3000 nano-UPLC system (120 min per sample) [45]. Proteomic data were used to infer vaginal microbiome composition and host secretome as described previously [45]. Briefly, proteins were identified using two different databases: (i) UniProt database restricted to human and microbial entries (73,910,451, release August 2017) and filtered using the MetaNovo pipeline [48] and (ii) human proteins combined with a vaginal metagenome-based database obtained from the study of Afiuni-Zadeh et al. [49]. Taxonomy was assigned using UniProt, and relative abundance of each taxon was determined by aggregating the intensity-based absolute quantification (iBAQ) values of all proteins identified for each taxon. The taxonomic analysis of the proteins identified using both of the databases was similar; however, the first database resulted in the identification of a greater number of proteins compared to the vaginal metagenome database. The first database also showed a high degree of similarity compared to 16S rRNA gene sequence data from a subset of the same women and was thus used for downstream analysis [45]. As the majority of taxa identified had ≤ 2 proteins detected, a more stringent cut-off was applied to include only taxa with ≥ 3 detected proteins or 2 proteins detected in multiple samples. For LDH analysis, log2-transformed iBAQ values for bacterial L-LDH and D-LDH were aggregated and the participants categorized into low (< median) or high (≥ median) groups. The presence of BV was assessed by Nugent score (Nugent BV) [5], with a score ≥ 7 considered as BV, 4–6 as an intermediate and ≤ 3 as non-BV [50]. This study evaluated samples from 113 women recruited from the Cape Town study site who had sufficient sample and metadata available for analysis.

The limma R package [51] was then used to identify differentially abundant host proteins between these categories. Logistic regression was used to adjust for potential confounders including Lactobacillus spp. relative abundance, STIs, prostate-specific antigen and contraceptives.

Culture and LA treatment of cervicovaginal epithelial cells

The human ectocervical (Ect1/E6E7) and vaginal (VK2/E6E7) cell lines (purchased from the ATCC) were cultured and seeded for treatment on transwell supports as previously described [41]. Average LA concentration in cervicovaginal fluid is 1% (+/−0.2%) [36]; however, the multilayered structure of the epithelium and the presence of mucus means viable epithelial cells in basal layers are likely to be exposed to concentrations < 1%. Indeed, previous optimisations indicated 0.3% LA is the highest concentration tolerated by monolayers of cervicovaginal epithelial cells without eliciting cytotoxicity [41] and was thus utilised for these analyses. Cells were treated for 1 h with media containing 0.3% L-LA or D-LA (Sigma-Aldrich) at either pH 3.9 or 7 or low pH media control (pH 3.9, media acidified by hydrochloric acid, HCl) added to the apical compartment. HCl-containing media was replenished during the 1-h treatment to maintain a pH < 4.5 [41]. Following treatment, media was removed and cells washed with calcium and magnesium-free PBS (PBS−, Life Technologies, Carlsbad, CA). Fresh media was added before returning to culture for either a further 4 h (for RNA analysis) or 24 h (for all other experiments). Epithelial barrier integrity was assessed by transepithelial electrical resistance (TEER) using a Millicell-ERS Voltohmmeter and probe (Millipore, MA). The unit area resistance (Ω × cm2) was calculated by multiplying the sample resistance by the membrane area.

Impact of LA treatment on epithelial barrier gene expression in FRT epithelial cells

RNA was extracted using the RNeasy mini kit (Qiagen, Hilden, Germany) as previously described [41]. To determine gene expression in LA-treated cervicovaginal epithelial cells, next-generation RNA-Seq analysis was performed at the Australian Genome Research Facility (AGRF, Melbourne, Australia) using the HiSeq 2500 NGS platform (Illumina, San Diego, CA) and an Illumina bcl2fastq 2.20.0.422 pipeline with an output of 100 bp paired-end reads. Sequence read datasets were checked for quality using FastQC [52]. Reads were trimmed based on quality scores using the Trim Reads tool with default parameters within CLC Genomics Workbench 11.0 (CLC; Qiagen). Overlapping paired reads were merged, and then, all reads were mapped against the Homo sapiens reference genome (build version hg38) using the RNA-Seq tool in CLC. Differentially expressed gene (DEGs) analysis was performed in Degust (http://degust.erc.monash.edu/), using the voom/limma method [53, 54]. A false discovery rate (FDR) cut-off of < 0.05 and a Log2 fold change (FC) ≥ 0.5 were used to identify genes with significantly altered expression. Targeted interrogation of the expression of genes related to the tight junction barrier was performed using a gene list derived from the RT2 Profiler™ Human Tight Junctions PCR Array (Qiagen). To identify and visualise significantly overrepresented gene ontologies among the DEGs, we utilised the BiNGO application [55] within Cytoscape v.3.8.2 [56] to perform a hypergeometric test using the Benjamini and Hochberg FDR correction. Unranked lists of DEGs were tested against the July 02, 2021 Gene Ontology release (http://geneontology.org) [57].

Altered expression of junctional molecules was confirmed by qRT-PCR using RT2 SYBR Green Mastermix (Qiagen) and primers specific for the following tight junction factors and housekeeping genes; claudin-1 (CLDN1, Cat no. PPH02779A), claudin-4 (CLDN4, Cat no. PPH07330D), occludin (OCLN, Cat no. PPH02571B), junctional adhesion molecule A (F11R, Cat no. PPH02605A), zona occludens 2 (TJP2, Cat no. PPH09978B) and the housekeeping genes glyceraldehyde 3-phosphate dehydrogenase (GAPDH, Cat no. PPH00150F), hypoxanthine phosphoribosyltransferase 1 (HPRT1, Cat no. PPH01018C) and ribosomal protein lateral stalk subunit P0 (RPLP0, Cat no. PPH21138F; all from Qiagen). Amplification was performed with a thermocycle of 95 °C for 10 min, 40 cycles of 95 °C for 15 s and 40 °C for 60 s, followed by melt curve analysis. Primer specificity was verified by melt curve analysis, which indicated a single amplification product with a unique melting temperature for each gene target. Ct values were analysed to determine relative gene expression, which was standardised to the average expression of the housekeeping genes and calculated using the 2−ΔΔCT method [58].

Analysis of tight junction protein expression and localisation

Levels of tight junction proteins in ectocervical epithelial cells were determined by Western blot using primary antibodies against claudin-1, claudin-4, tight junction protein (TJP) 1 and TJP2 (Cat no. 51-9000, 36-4800, 33-9100 and 71-1400, respectively, all from Thermo Fisher Scientific), β-actin (Cat no. ab8224, Abcam) and secondary antibodies Alexa Fluor 680-goat anti-mouse IgG (Thermo Fisher Scientific) and Alexa Fluor 800-donkey anti-rabbit IgG (Li-COR, Lincoln, NE). Proteins were visualised with an Odyssey Infrared Imaging system and analysed with Image Studio Lite software (both from Li-COR). Protein levels were normalised to β-actin protein intensity.

Immunofluorescence and confocal microscopy

Claudin-4 protein location in Ect cells was visualised by immunofluorescence and confocal microscopy. Epithelial cells were fixed with ice-cold methanol, rehydrated with 1% foetal bovine serum (FBS) in PBS− and permeabilised with 0.2% Triton X-100 in PBS− for 10 min. After blocking with 3% BSA in PBS−/0.2% Triton X-100, cells were incubated with rabbit polyclonal anti-claudin-4 antibody (Invitrogen, Cat no. 36-4800, 1:25 dilution) and then Alexa Fluor 488-goat anti-rabbit IgG secondary antibody (Invitrogen, Cat no. A-11008, 1:200). Nuclei were visualised using Hoechst 33342 stain (Thermo Scientific). Transwell membranes and attached cells were excised and mounted on slides with ProLong Gold mounting medium (Invitrogen). Confocal Z-stacks were captured at 0.2 μm per section, 25 sections in total using a Nikon A1R confocal microscope (Nikon Corporation, Japan) at 60× magnification with an oil-immersion lens. Images were selected at random (at least 3 stacks per slide) and captured for quantitative analysis with ImageJ software (National Institute of Mental Health). Images were thresholded to eliminate background fluorescence. The sum of fluorescence intensity was calculated for the stack, and mean fluorescence intensity (MFI) was determined.

Transcriptomic analysis of vaginal epithelial cells exposed to bacterial culture supernatants

Bacterial culture supernatants were generated by inoculating 10 ml of culture media with 1 ml of 1 × 107 colony-forming units (CFU) of L. crispatus (ATCC 33197), L. jensenii (ATTC 25258), L. iners (ATTC 55195) or G. vaginalis (ATCC 14018). NYC-III media was used for propagation of L. crispatus, L. jensenii and L. iners, whilst Tryptic Soy Broth (TSB) was used for G. vaginalis. CFU to OD600 measures were pre-calculated for each strain. Cultures were grown anaerobically for 48 h to late stationary phase. L. crispatus, L. jensenii, L. iners and G. vaginalis reached 5.2 × 108 CFU/ml, 5.8 × 108 CFU/ml, 4.1 × 108 CFU/ml and 1.2 × 109 CFU/ml, respectively. The 10 ml of culture supernatant was clarified by centrifugation at 3000 × g for 10 min, sterile filtered (0.2 μm filter) and stored at −20 °C until use. Bacterial supernatants were diluted to 20% (v/v) in complete VK2 cell culture medium and added to VK2 cells for 13 h. VK2 cells were serum starved for 24 h in keratinocyte serum-free base media only, prior to addition of 20% bacterial culture supernatants. After the exposure, media was removed, cells were washed once with 1× PBS and 300 μl RNAlater (QIAGEN) was added to wells. Cells were mechanically detached and stored at −80 °C before total RNA extraction. Total RNA was extracted using the MasterPure™ Complete DNA and RNA purification kit (Lucigen, WI, USA) according to manufacturer’s instructions. Methodology for transcriptomic analyses by RNASeq including library preparation, RNA sequencing and read mapping has been detailed previously [29]. LA concentration within bacterial supernatant-containing media was quantified by using the D/L lactic acid assay kit (R-Biopharm AG, Darmstadt, Germany, Cat no. 11112821035). Concentration of the protonated form of LA within treatment media was calculated from the sample pH and LA concentration using the Henderson-Hasselbalch equation as previously described [36]. Total lactate and protonated LA were normalised to the final CFU/ml for L. crispatus, L. jensenii, L. iners and G. vaginalis achieved at the stationary phase.

Results

Bacterial lactate-dehydrogenase relative abundance is increased in cervicovaginal secretions from women with Lactobacillus-dominated microbiota

A previous proteomic analysis of cervicovaginal secretions from 41 women identified an increased relative abundance of cornified envelope factors (INVO and SPR1A), which are involved in maintaining epithelial barrier integrity, in women with a Lactobacillus-dominated as compared to a G. vaginalis-dominated microbiota [46]. Furthermore, levels of lactate dehydrogenase (LDH, a key enzyme in LA synthesis responsible for catalysing the conversion of pyruvate to lactate) from L. iners were positively associated with these factors. We therefore explored the relationship between LA isomers produced from different Lactobacillus spp. and proteins specifically involved in epithelial barrier integrity in vivo. We conducted a proteomic analysis utilising cervicovaginal fluid samples collected from 113 HIV-uninfected, nonpregnant women enrolled in the WISH study in Cape Town, South Africa (Table 1). Half of these women were Nugent-BV positive (50%, Table 1), whilst the remainder typically exhibited a microbiota dominated by L. iners or L. crispatus (Fig. 1A). LC-MS/MS analysis detected 3186 different peptides, 55% of which were of bacterial origin, and 8 distinct LDH proteins produced from bacterial species were identified in these samples (Fig. 1B). Bacterial LDH abundance was analysed as an indicator of LA levels given our previously findings that LDH relative abundance produced by clinical Lactobacillus isolates correlates with lactate concentrations in vitro [59]. High levels of host-derived LDH, which has only 39% sequence homology with LDH derived from gram-positive bacteria [60], were also detected but were excluded from the analysis. As expected, bacterial LDH protein was substantially less abundant in cervicovaginal samples from women with BV as compared to those with a microbiota dominated by Lactobacillus spp. (Fig. 1). Cervicovaginal secretions from women with a microbiota dominated by L. iners contained primarily L-LDH, whilst both L- and D-LDH were detected in women with L. crispatus-dominated microbiome (Fig. 1B), consistent with the known differential production of LA isomers by these species [34, 43]. L. iners was detected in a portion of women with Nugent BV at a relative abundance of up to 86% (Fig. 1A), and in these women, L. iners-derived L-LDH was also detected. Consistent with previous findings [2, 61], L. crispatus was rare in women with BV, and consequently, minimal D-LDH was detected.

Table 1 Demographic and clinical characteristics of cohort
Fig. 1
figure 1

Higher relative abundance of bacterial LDH in women with a Lactobacillus-dominated microbiome. The vaginal microbiome composition of 113 young women from Cape Town, South Africa, was determined using liquid chromatography-tandem mass spectrometry, and participants were classified as being bacterial vaginosis (BV) negative, positive, or intermediate (Int) using Nugent-BV criteria. Proteins were identified using MaxQuant, and a custom database was generated using de novo sequencing to filter the UniProt database. Taxonomy was assigned using UniProt, and relative abundance of each taxon was determined by aggregating the intensity-based absolute quantification (iBAQ) values of all proteins identified for each taxon. The relative abundance of the 20 most abundant bacterial species is indicated for each participant in (A). The relative abundance of L- and D-lactate dehydrogenase (LDH) protein derived from various bacterial sources was measured in vaginal secretions using metaproteomic analyses for each participant (B)

We found that total bacterial LDH relative abundance correlated inversely with cervicovaginal pH (Spearman’s Rho = −0.44; p < 0.001; Supplementary Fig. 1), which agrees with previous observations [46] and is consistent with the known effect of LA in maintaining low vaginal pH [36]. In contrast, host L-LDH showed a weak positive association with vaginal pH (Spearman’s Rho = 0.265, p = 0.007; data not shown). Taken together, these data indicate a higher abundance of bacterial-derived LDH in cervicovaginal samples from women with Lactobacillus-dominated microbiota where vaginal LA levels are known to be high and pH levels low.

High LDH relative abundance is associated with increased expression of epithelial barrier proteins in vivo

To examine the relationship between bacterial LDH abundance and host protein expression, we stratified women into high and low bacterial LDH categories based on the median relative LDH abundance of the cohort respectively. An analysis of host protein relative abundance identified 96 proteins which showed a significant differential abundance between high and low L-LDH groups (FDR adjusted p < 0.05; Supplementary Table 1), with 31 of these proteins (32%) having epithelial barrier function-related gene ontologies (Fig. 2). Although detection of L. crispatus, and consequently D-LDH protein, was lower in this cohort, 4 proteins were differentially abundant between women with high and low D-LDH (FDR adjusted p < 0.05; Supplementary Table 1). These included two proteins with gene ontology terms related to epithelial barrier function: keratin 8 and dermokine. Logistic regression analysis was performed to assess the relationship between L-LDH groups and expression of these epithelial barrier-associated proteins after adjusting for multiple potential confounders including Lactobacillus relative abundance, the presence of STIs, contraceptive use, and detection of prostate-specific antigen (PSA, an indicator of recent vaginal intercourse). We found that 11 of the 31 proteins remained significantly associated with L-LDH abundance, including tight junction protein 1 (TJP1/ZO-1; Fig. 2 and Supplementary Table 2). Similarly, 10 of 31 proteins remained significantly associated with L-LDH abundance following adjustment for BV status, STIs, contraceptive use and PSA. To account for protease production by nonoptimal bacteria that may contribute to epithelial barrier degradation and dysfunction, we aggregated the relative abundance of the 23 proteins with protease/peptidase activity produced by G. vaginalis, Megasphaera spp., Mobiluncus mulieris, Prevotella spp. and Sneathia amnii that were detected in our samples. Whilst total protease/peptidase relative abundance was inversely associated with epithelial barrier proteins, when we evaluated the association between L-LDH and epithelial barrier proteins, adjusting for protease/peptidase relative abundance, we found that 30/31 remained significantly associated with L-LDH. These analyses provide in vivo evidence of an association between cervicovaginal relative abundance of the key enzyme responsible for LA production and levels of epithelial junctional proteins that govern barrier integrity.

Fig. 2
figure 2

L-LDH relative abundance in the FRT is associated with expression of epithelial junction proteins. Unsupervised hierarchical clustering of intensity-based absolute quantification (iBAQ) values for epithelial barrier-related proteins that were found to be differentially abundant between women (n = 113) with levels of L-LDH protein above (high, purple) or below (low, yellow) the median value for the cohort. The moderated t-test (limma package, R) was used to identify proteins significantly associated with L-LDH relative abundance following false discovery rate adjustment for multiple comparisons. *Proteins that remained significantly associated with L-LDH after adjusting for Lactobacillus relative abundance, as well as the presence of sexually transmitted infections, contraceptive use, and detection of prostate-specific antigen

Lactic acid treatment of ectocervical epithelial cells enhances epithelial barrier integrity

Given the in vivo observations supporting an association between LA and cervicovaginal epithelial integrity, we assessed whether LA produced by vaginal Lactobacilli spp. was able to directly enhance epithelial barrier integrity in vitro. As discussed above, lactobacilli can produce both L and D isomers of LA dependent on species, with the biologically active protonated form of LA predominating at pH below 3.9 [38, 41, 42]. Ectocervical epithelial cells were therefore apically treated with a physiologically relevant concentration of 0.3% L- or D-LA at pH 3.9 and the change in epithelial barrier integrity assessed by TEER. Treatment of cells with either isoform of LA at pH 3.9 elicited a significant enhancement of barrier integrity as compared to untreated cells (Fig. 3), which was not observed when cells were treated with media acidified to the same pH with HCl. Ethylenediamine tetraacetic acid (EDTA) served as a positive control for barrier disruption in these experiments and elicited an expected decrease in TEER. Similar to our previous findings regarding the pH dependency of LA’s bioactivity [38, 41], no effect on TEER was observed when L- and D-LA containing media were neutralised to pH 7.0, indicating the epithelial barrier-enhancing effect of LA is mediated by the protonated form and not the lactate anion.

Fig. 3
figure 3

Protonated LA increases the barrier integrity of ectocervical epithelial cells. Ectocervical cells (Ect1) seeded into transwells were apically treated for 1 h with media alone (untreated, UT) or media containing 0.3% L- or D-lactic acid (LA) at either low (pH 3.9) or neutral (pH 7) pH, low pH media alone (HCl, pH 3.9), or 2 mM EDTA. Epithelial barrier integrity was assessed by transepithelial electrical resistance (TEER) measured at baseline (T0) and 24 h after treatment (T24) expressed as the TEER ratio (T24/T0). Graph shows median and interquartile range for each condition from 4 independent experiments. *p < 0.05 compared to untreated control, as determined by Mann-Whitney U-test

LA elicits altered expression of genes related to epithelial barrier integrity

To identify the mechanisms of LA-mediated barrier enhancement, we performed transcriptomic analysis of ectocervical epithelial cells treated with L- and D-LA (pH 3.9) using RNA-Seq. This analysis revealed 383 and 540 DEGs between untreated cells and those treated with either L-LA or D-LA, respectively (Fig. 4A and Supplementary Table 3). In contrast, treatment of ectocervical cells with media acidified to pH 3.9 with HCl had negligible effect on gene expression, with only one gene (HSPA6) found to be significantly differentially expressed compared to untreated cells (orange box in Fig. 4A). These data suggest that changes to gene transcription elicited by LA are unlikely to be due to low pH alone. No significant differences were observed in DEGs detected between L- and D-LA-treated cells using our stringent analysis criteria (upper boxes, Fig. 4A), suggesting that both isomers of LA mediate an analogous effect on the transcriptome of ectocervical epithelial cells, consistent with their similar effects on increasing epithelial barrier integrity (Fig. 3). Indeed, of the 383 DEGs elicited by L-LA treatment, 231 (60.3%) of the same genes were also significantly altered by D-LA treatment (Fig. 4A). This was further supported by the clustering of L- and D-LA-treated samples on a multidimensional scaling (MDS) plot (Supplementary Fig. 2). These data indicate treatment with L- and D-LA, but not low pH alone, elicits substantial gene changes to immortalized ectocervical cells, and that overlapping transcriptomic signatures are observed with both isomers of LA.

Fig. 4
figure 4

LA treatment of ectocervical epithelial cells increases the expression of genes relevant to barrier integrity. Ectocervical epithelial cells were apically treated with 0.3% L- or D-lactic acid (LA) or low pH media alone (HCl) for 1 h and transcriptomic changes assessed at 4-h post treatment by RNA-Seq. A Number of genes differentially expressed between the treatments is shown next to arrows indicating comparison groups as assessed by voom-limma analysis (Log2 fold change > 0.5, false discovery rate [FDR] < 0.05). There were 231 genes which were differentially expressed by both L-LA and D-LA as compared to untreated cells. B Gene ontologies enriched by LA treatment. For visual simplicity, ontologies enriched with a FDR < 0.01 are shown; all ontologies enriched with a FDR < 0.05 are shown in Supplementary Table 2. C Expression of genes related to epithelial tight junctions relative to untreated cells ranked by magnitude of up- or downregulation (left panel) or adjusted p-value (FDR). Genes highlighted in orange were those selected for further characterisation

Gene ontology analysis was performed to investigate pathways which were targeted by LA treatment. This revealed that LA treatment altered a broad range of cellular components, functions and processes (Fig. 4B and Supplementary Table 4) with 32/334 (9.6%) of pathways that were significantly enriched by L-LA treatment being associated with barrier function (Supplementary Table 4). A similar proportion of pathways enriched by D-LA treatment were also barrier function related (15/191; 7.9%). To assess the impact of LA treatment on epithelial tight junctions specifically, we separately evaluated expression of key genes known to encode critical tight junction proteins and found a substantial portion of these factors were significantly altered by LA treatment (Fig. 4C and Supplementary Table 5). The most significantly DEGs included those with established roles in facilitating tight junctions including CLDN1, CLDN4, OCLN, TJP2 and F11R (Fig. 4C).

To validate the differential expression of tight junction genes by LA treatment observed in the RNA-Seq analysis, we performed qRT-PCR analysis of a subset of genes in ectocervical cells treated with LA (pH 3.9) under the same conditions used to generate the RNA-Seq data. Given the similar impact of both L- and D- LA on transcriptional changes observed above, we limited the analysis to L-LA for these validation experiments. Consistent with the RNA-Seq data, qRT-PCR analysis demonstrated significantly increased expression of CLDN1, CLDN4, OCLN and F11R following LA treatment (Fig. 5, p < 0.05 for all), whilst expression of TJP2 was also significantly increased by L-LA in the qRT-PCR analysis. Interestingly, qRT-PCR analysis also showed increased expression of CLDN1 and CLDN4 due to low pH alone, but the magnitude of this effect was lower than that observed for L-LA (Fig. 5). We also analysed expression of the two most significant DEGs, CLDN1 and CLDN4, in primary cervicovaginal cells and observed a similar effect of L-LA on expression as in ectocervical cells (Supplementary Fig. 3). Taken together, these molecular analyses indicate that physiological concentrations of both L- and D-LA at pH 3.9 upregulate the expression of genes involved in cell barrier and tight junction formation which likely contributes to enhanced barrier function of ectocervical epithelial cells.

Fig. 5
figure 5

qRT-PCR analysis of expression of tight junction genes in ectocervical cells following L-LA treatment. Ectocervical epithelial cells were apically treated with 0.3% L-lactic acid (L-LA) or low pH media alone (HCl) for 1 h, and gene expression of tight junction factors claudin-1 (CLDN1), claudin-4 (CLDN4), occludin (OCLN), zona occludens-2/tight junction protein-2 (TJP2) and F11R was assessed at 4 h post treatment by RNA-Seq (A) or qRT-PCR (B). A Fold change in expression of tight junction genes in ectocervical epithelial cells treated with L-LA, D-LA or HCl-treated cells relative to untreated, as determined by RNA-Seq analysis. * and *** represent false discovery rate < 0.05 and < 0.001, respectively. B Relative expression (fold change) of tight junction genes compared to untreated cells (indicated by dotted line) as determined by qRT-PCR. Graph shows mean +/− SEM from n = 5–6 independent experiments. * and ** represent p < 0.05 and < 0.01, respectively, compared to untreated cells as determined by the Mann-Whitney U-test. NS, not significant

LA upregulates the expression of tight junction proteins claudin-1 and claudin-4

To confirm the upregulation of tight junction gene expression was associated with increased levels of tight junction proteins, ectocervical cells were treated with 0.3% L-LA pH 3.9 and levels of the most highly expressed factors claudin-1 and -4 assessed by Western blot. Consistent with the gene expression data, protein levels of both claudin-1 and -4 were significantly higher in L-LA-treated cells as compared to untreated cells (Fig. 6A and B, p = 0.01 and 0.03, respectively), and higher levels of claudin-4 were also observed in cells treated with low pH media alone (Fig. 6B). Protein levels of TJP2/ZO-2 were also significantly increased in L-LA-treated cells (Fig. 6C), which was consistent with the significantly increased expression of TJP2 found in the qRT-PCR analysis. Finally, to confirm that the protein levels observed in Western blot reflected relevant proteins localised at regions of cell-cell contact, we analysed ectocervical cells by fluorescence microscopy and observed focussed localisation of claudin-4 at intercellular regions in ectocervical cells treated with L-LA and a significantly increased intensity of staining for claudin-4 in L-LA-treated cells as compared to untreated samples (Fig. 6D). Thus, LA acts directly on FRT epithelial cells to increase the expression of junctional proteins at sites of cell-cell attachment, which likely contributes to strengthening of the epithelial barrier.

Fig. 6
figure 6

Increased expression of tight junction proteins in ectocervical cells by L-LA. Protein levels of barrier proteins claudin-1 (A), claudin-4 (B) and tight junction protein-2 (TJP2) (C) in ectocervical epithelial cells stimulated apically with 0.3% L-LA pH 3.9 (LA) or low pH media (HCl, pH 3.9) for 1 h followed by lysis of cells at 24 h and assessed by Western blot analysis. Target protein levels were standardised to β-actin in the same sample. D Expression of claudin-4 was visualised in ectocervical epithelial cells treated as above by immunofluorescence and confocal microscopy. Images from a single experiment representative of n = 4 replicates are shown (left), and mean fluorescence intensity (MFI) of claudin-4 signal intensity was assessed using ImageJ software and expressed relative to untreated samples. All graphs show mean +/− SEM (expressed as fold change compared to untreated cells) from 3 to 5 independent experiments. Significance was assessed by the Mann-Whitney U-test; *, **p < 0.05 and 0.01, respectively

Bacterial supernatants containing LA increase expression of tight junction genes

To evaluate whether LA produced by vaginal Lactobacillus spp. could similarly increase expression of tight junction proteins in the context of a complex milieu of bacterial components, VK2 epithelial cells were cultured with media containing 20% (v/v) filtered supernatant from L. crispatus, L. jensenii, L. iners or G. vaginalis cultures and gene expression analysed by RNA-Seq. Expression analysis of key tight junction genes described in Fig. 4C revealed a substantial portion of the 54 tight junction genes in this list which were also detected within the bacterial culture-treated cells were differentially expressed following treatment of vaginal epithelial cells with culture supernatants from L. crispatus and L. jensenii cultures (59% and 56%, respectively, Fig. 7A and Supplementary Table 6). These DEGs included CLDN4, OCLN, TJP2 and F11R genes shown to be most substantially altered following LA treatment (Fig. 7B). In contrast, less substantial changes to tight junction gene expression were observed following treatment with L. iners culture supernatants, with only 19% of genes significantly differentially expressed following treatment and no change to TJP2 and F11R gene expression (Fig. 7 and Supplementary Table 6). Minimal changes to tight junction gene expression were observed following treatment with G. vaginalis cultures. Whilst the growth of L. iners cultures (4.2 × 108 CFU/mL after 48 h culture) was slightly less than that of L. crispatus and L. jensenii (5.2 and 5.8 × 108 CFU/ml, respectively), adjusting RNASeq data for these small differences did not materially alter the results (analysis not shown). Media containing 20% L. crispatus and L. jensenii culture supernatants contained substantially higher levels of both total lactate and protonated LA as compared to that derived from L. iners or G. vaginalis cultures, even after adjusting for bacterial growth differences (Fig. 7C). These data indicate that supernatants from Lactobacillus spp. which produce high levels of LA upregulate the expression of tight junction genes in cervicovaginal epithelial cells in vitro, supporting a specific role for LA in enhancing epithelial barrier integrity.

Fig. 7
figure 7

Lactobacillus cultures producing high levels of LA upregulate tight junction gene expression. VK2 epithelial cells were treated with media containing 20% (v/v) filtered supernatant from cultures of L. crispatus, L. jensenii, L. iners or G. vaginalis bacteria. The impact of treatment on expression of key tight junction genes was analysed by RNA-Seq, and volcano plots of changes in expression following treatment are shown in (A), with red dots indicating significantly DEGs (FDR < 0.05). B Fold change in expression of claudin-4 (CLDN4), occludin (OCLN), zona occludens-2/tight junction protein-2 (TJP2) and F11R genes relative to cells treated with media alone (mean, from n = 3 replicates; *, ** and ***false discovery rate < 0.5, 0.01 and 0.001, respectively). C Concentration of total lactate (black bars) and protonated LA (grey bars) in treatment media, the latter calculated using the Henderson-Hasselbalch equation and the treatment media pH (range 6.6–7.1). Values were normalised to account for slight variations in growth rate of bacterial cultures

Discussion

This study reveals a novel mechanism by which cervicovaginal Lactobacillus spp. enhance FRT epithelial barrier integrity through the action of the metabolite LA that directly increases epithelial barrier integrity and expression of tight junction proteins in vitro. The physiological relevance of these findings was indicated in a cohort study which demonstrated an association between high bacterial LDH levels and increased relative abundance of epithelial barrier proteins in young women. Notably, this association was observed after adjusting for several confounders including BV status, nonoptimal bacterial protease relative abundance and Lactobacillus spp. abundance, indicating a direct association between L-LDH, which produces lactic acid, and barrier function in vivo. Together, these data suggest LA produced by vaginal lactobacilli may play an important role in enhancing the integrity of the cervicovaginal epithelial barrier in vivo to protect against infection with HIV and other STIs.

Experiments utilising human cervical explant tissue indicate HIV can penetrate both the single-layer columnar epithelial cells of the endocervix and upper FRT as well as the squamous, multilayered ectocervical epithelium of the ectocervix and vagina [20, 62], although studies in macaques indicate preferential virus transduction through the ectocervix and vagina [63]. Cervicovaginal mucus, particularly that produced by women with an L. crispatus-dominated microbiota, is able to effectively trap HIV virions [64], and this may represent an important barrier to HIV invasion in regions such as the columnar epithelium of the endocervix [20]. In the squamous epithelium of the ectocervix and vagina, HIV can easily penetrate the superficial epithelial layers of the stratum corneum where intercellular junctional molecules are not intact. However, penetration of HIV beyond the lower layers of epithelial cells expressing intact tight and adherens junction proteins is highly restricted [20, 65]. This suggests the integrity of mucosal epithelial junctions plays an important role in preventing HIV transmission. The data presented here demonstrate a direct LA-mediated enhancement of tight junction protein expression and increased epithelial barrier integrity in an ectocervical epithelial cell line, which may be relevant to HIV infection risk in vivo.

Enhanced barrier integrity of cervicovaginal epithelial cells cocultured with probiotic lactobacilli or conditioned media from Lactobacillus spp. cultures in vitro has been demonstrated previously [31, 66]. However, here we show that culture supernatants from Lactobacillus spp. producing high levels of LA increased tight junction gene expression in cervicovaginal epithelial cells to a much greater extent than L. iners, which produced substantially lower levels of LA. This suggests LA may be a major mediator of these barrier-enhancing effects. Conditioned media from bacterial cultures contain many soluble and insoluble components that may modulate epithelial barriers. Therefore, we used a reductionist approach to examine the specific activities of LA, added directly to cell culture media, on cervicovaginal cells. These studies demonstrate that it is the protonated form of LA, which predominates at pH < 3.86 present in the FRT of women with optimal Lactobacillus-dominated microbiota, which possesses epithelial barrier-promoting activity. These effects of protonated LA are uniquely relevant to the FRT as LA produced by intestinal lactobacilli in the gut would be present mostly as lactate due to the near neutral pH of the gut [67]. Whilst it is possible other bacterial-derived products may also contribute to the barrier-promoting properties of Lactobacillus spp. culture supernatants, our data suggest LA may be a major bioactive component of this effect.

Here, we observed a similar effect of L- and D-isomers of LA on enhancing epithelial barrier integrity and junctional molecule expression, consistent with our previous findings of comparable immunomodulatory and virucidal properties of both isomers [38, 41]. Indeed, our RNA-Seq analysis indicated that L- and D-LA elicited overlapping transcriptomic changes to ectocervical epithelial cells. L. crispatus, which produces predominantly D-LA, has been suggested to be the more favourable commensal Lactobacillus spp. due to its greater apparent stability and association with lower levels of genital inflammation and HIV acquisition [10]. D-LA also elicits more potent inhibition of Chlamydia trachomatis infection of FRT epithelial cells compared to L-LA in vitro [29]. However, L-LA has been shown to be more potent than D-LA in inactivating HIV at threshold concentrations in vitro [38]. Therefore, whilst our in vitro analysis did not identify differential effects of L- and D-LA on epithelial barrier integrity it remains possible LA isomers may have context and pathogen-specific differential activity in vivo. It is important to also consider that the total concentration of LA (both L- and D-LA) produced by lactobacilli can differ between species [34], which in turn influences pH and thus levels of the biologically relevant protonated form of LA [39]. Indeed, our data indicating L. crispatus and L. jensenii cultures produced substantially more LA than L. iners in vitro was consistent with their ability to upregulate tight junction gene expression to a greater extent.

The in vitro analyses presented in this study demonstrate an epithelial barrier-enhancing effect of LA, whilst our cohort study findings support the in vivo relevance of these observations in the complex FRT environment. Our data indicate that women with a greater abundance of bacterial LDH proteins (measured in ex vivo cervicovaginal samples) have increased expression of proteins important for maintaining epithelial barrier function and tight junction integrity, suggesting that these are physiologically relevant observations. Additionally, the finding that non-Lactobacillus taxa were also producing LDH, and that the relationship between LDH and epithelial barrier proteins was independent of BV status, Lactobacillus spp. relative abundance and protease/peptidase production by nonoptimal bacteria, suggests a direct relationship between LDH and barrier function. This suggests that although BV-associated bacteria may produce proteases and elicit host inflammatory responses that influence epithelial barrier integrity, the relationship between LDH and barrier function was independent of these factors. Our findings from a larger cohort (n = 113) of young South African women at increased risk of HIV support and extend previous observations [46] in Kenyan (n = 10) and North American (n = 31) women. Zevin et al. demonstrated enrichment of proteins associated with epidermal development and barrier function in women with Lactobacillus-dominated compared to G. vaginalis-dominated microbiota as well as a significant positive association between L-LDH abundance and INVO and SPR1A, which are known to act as scaffolding for epidermal layers [46]. Taken together, these in vitro and cohort data suggest a mechanistic link between bacterial-derived LA and cervicovaginal epithelial barrier integrity.

In this study, LA concentrations were inferred from in vivo levels of LDH protein in the metaproteome of women, which enabled the discrimination between relative abundance of host and bacterial LDH. Although the contribution of host versus bacterial LDH to LA measurements is difficult to differentiate, we have previously shown LDH relative abundance correlates directly with levels of total lactate production by lactobacilli in vitro [59]. Furthermore, high bacterial LDH abundance was associated with low vaginal pH. This is in contrast with host LDH, which showed a weak positive association with vaginal pH, consistent with previous findings [68] and may suggest higher levels of cell death (and hence host LDH release) at higher vaginal pH. Furthermore, vaginal bacteria, and not host cells, are considered the primary source of vaginal LA in vivo [43, 69]. Accurate measurement of cervicovaginal LA levels in clinical studies is challenging due to the fact that exposure of cervicovaginal fluid samples to the aerobic environment during collection and processing results in loss of CO2 and a subsequent increase in pH that may underestimate the in vivo concentration of the biologically relevant, protonated form of LA [36]. The low prevalence of L. crispatus and other D-LDH-producing bacteria in the cohort included in this study may have limited the identification of proteins associated with D-LDH levels. We performed proteomic analysis of cervicovaginal fluid, and although there is a 63% overlap between the FRT tissue and mucus proteome [70], not all cellular junctional molecules are detected using this approach. The majority of the in vitro mechanistic data presented here were performed in an immortalised ectocervical epithelial cell line, and whilst heightened expression of CLDN1 and CLDN4 genes was also observed in primary cervicovaginal epithelial cells following LA treatment (Supplementary Fig. 3), confirmation of other results in primary cells, or ideally primary multilayer 3D epivaginal tissue models, is warranted.

Our findings support the potential for LA to be used therapeutically to enhance the natural defences of the FRT mucosa and protect against pathogen infection. The combined abilities of LA to protect against pathogen-mediated inflammation, enhance the epithelial barrier and kill both BV-associated bacterial and HIV and other viruses [71] suggest it may be an ideal agent to promote FRT health. Although a plethora of over-the-counter LA-containing formulations are available for use in the FRT, the formulations and LA content of these products vary considerably, and data on their ability to alter sexual and reproductive health outcomes is limited [72]. However, a recent clinical trial of LA-containing gels in women with recurrent BV demonstrated BV resolution in 47% of participants [73]. Whilst this effect was inferior to that mediated by the standard of care treatment metronidazole alone [74] (70% resolution), LA-containing gels were preferred by women, and side effects were less common [73, 75]. The ability of the probiotic L. crispatus CTV-05 to modulate FRT health was recently evaluated in a cohort of women with BV where it was found to significantly reduce the rate of BV-recurrence [76]. Thus, the potential for LA to mediate beneficial in vivo effects is an area warranting future research.

Conclusions

By utilising a combination of in vitro mechanistic approaches and ex vivo analyses of FRT samples collected from a clinical cohort of young women, this study has identified a novel effect of the Lactobacillus metabolite LA on enhancing the integrity of the epithelial barrier through directly increasing levels of junctional molecules. These findings provide important mechanistic insights into the association of Lactobacillus-dominated microbiota with favourable sexual health outcomes and indicate potential avenues for therapeutic modulation of the FRT environment to enhance genital health.

Availability of data and materials

All data generated during the current study are included in this article (and its supplementary information files). Raw sequence reads will be made available upon request. Metaproteomics data is derived from a previously published cohort study (Alisoltani et al. 2020 Microbiome 8:165) with data available online http://fgtdb.org/. The sequences generated and analysed during this current study were uploaded to the NCBI Sequence Read Archive (SRA) data repository, with project numbers PRJNA823848 and PRJNA835889.

References

  1. Tuddenham S, Ravel J, Marrazzo JM. Protection and risk: Male and female genital microbiota and sexually transmitted infections. J Infect Dis. 2021;223(Supplement_3):S222–S35.

    PubMed  PubMed Central  Article  Google Scholar 

  2. 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(Suppl 1):4680–7.

    CAS  PubMed  Article  Google Scholar 

  3. Fettweis JM, Brooks JP, Serrano MG, Sheth NU, Girerd PH, Edwards DJ, et al. Differences in vaginal microbiome in African American women versus women of European ancestry. Microbiology (Reading). 2014;160(Pt 10):2272–82.

    CAS  Article  Google Scholar 

  4. Aldunate M, Srbinovski D, Hearps AC, Latham CF, Ramsland PA, Gugasyan R, et al. Antimicrobial and immune modulatory effects of lactic acid and short chain fatty acids produced by vaginal microbiota associated with eubiosis and bacterial vaginosis. Front Physiol. 2015;6:164.

    PubMed  PubMed Central  Article  Google Scholar 

  5. McKinnon LR, Achilles SL, Bradshaw CS, Burgener A, Crucitti T, Fredricks DN, et al. The evolving facets of bacterial vaginosis: implications for HIV transmission. AIDS Res Hum Retroviruses. 2019;35:219–28.

    PubMed  PubMed Central  Article  Google Scholar 

  6. Elovitz MA, Gajer P, Riis V, Brown AG, Humphrys MS, Holm JB, et al. Cervicovaginal microbiota and local immune response modulate the risk of spontaneous preterm delivery. Nat Commun. 2019;10:1305.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  7. Fettweis JM, Serrano MG, Brooks JP, Edwards DJ, Girerd PH, Parikh HI, et al. The vaginal microbiome and preterm birth. Nat Med. 2019;25:1012–21.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  8. Payne MS, Newnham JP, Doherty DA, Keelan JA, Furfaro LL, Pendal NL, et al. A specific bacterial DNA signature in the vagina of Australian women in midpregnancy predicts high risk of spontaneous preterm birth (the Predict1000 study). Am J Obstet Gynecol. 2021;224:635–6.

    CAS  PubMed  Article  Google Scholar 

  9. 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;8:1781–93.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. Gosmann C, Anahtar MN, Handley SA, Farcasanu M, Abu-Ali G, Bowman BA, et al. Lactobacillus-deficient cervicovaginal bacterial communities are associated with increased HIV acquisition in young South African women. Immunity. 2017;46:29–37.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. Srinivasan S, Richardson BA, Wallis J, Fiedler TL, Dezzutti CS, Chirenje ZM, et al. Vaginal microbiota and HIV acquisition risk among African women. In: Conference on retroviruses and opportunistic infections. Boston; 2018.

  12. Hillier SL, Nugent RP, Eschenbach DA, Krohn MA, Gibbs RS, Martin DH, et al. Association between bacterial vaginosis and preterm delivery of a low-birth-weight infant. The Vaginal Infections and Prematurity Study Group. N Engl J Med. 1995;333:1737–42.

    CAS  PubMed  Article  Google Scholar 

  13. Brotman RM, Klebanoff MA, Nansel TR, Yu KF, Andrews WW, Zhang J, et al. Bacterial vaginosis assessed by gram stain and diminished colonization resistance to incident gonococcal, chlamydial, and trichomonal genital infection. J Infect Dis. 2010;202:1907–15.

    PubMed  PubMed Central  Article  Google Scholar 

  14. Atashili J, Poole C, Ndumbe PM, Adimora AA, Smith JS. Bacterial vaginosis and HIV acquisition: a meta-analysis of published studies. AIDS. 2008;22:1493–501.

    PubMed  Article  Google Scholar 

  15. Low N, Chersich MF, Schmidlin K, Egger M, Francis SC, van de Wijgert JH, et al. Intravaginal practices, bacterial vaginosis, and HIV infection in women: individual participant data meta-analysis. PLoS Med. 2011;8:e1000416.

    PubMed  PubMed Central  Article  Google Scholar 

  16. Cohen CR, Lingappa JR, Baeten JM, Ngayo MO, Spiegel CA, Hong T, et al. Bacterial vaginosis associated with increased risk of female-to-male HIV-1 transmission: a prospective cohort analysis among African couples. PLoS Med. 2012;9:e1001251.

    PubMed  PubMed Central  Article  Google Scholar 

  17. Kaushic C. HIV-1 infection in the female reproductive tract: role of interactions between HIV-1 and genital epithelial cells. Am J Reprod Immunol. 2011;65:253–60.

    CAS  PubMed  Article  Google Scholar 

  18. Blaskewicz CD, Pudney J, Anderson DJ. Structure and function of intercellular junctions in human cervical and vaginal mucosal epithelia. Biol Reprod. 2011;85:97–104.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. Pudney J, Quayle AJ, Anderson DJ. Immunological microenvironments in the human vagina and cervix: mediators of cellular immunity are concentrated in the cervical transformation zone. Biol Reprod. 2005;73:1253–63.

    CAS  PubMed  Article  Google Scholar 

  20. Carias AM, McCoombe S, McRaven M, Anderson M, Galloway N, Vandergrift N, et al. Defining the interaction of HIV-1 with the mucosal barriers of the female reproductive tract. J Virol. 2013;87:11388–400.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. Shattock RJ, Moore JP. Inhibiting sexual transmission of HIV-1 infection. Nat Rev Microbiol. 2003;1:25–34.

    CAS  PubMed  Article  Google Scholar 

  22. Yasen A, Herrera R, Rosbe K, Lien K, Tugizov SM. Release of HIV-1 sequestered in the vesicles of oral and genital mucosal epithelial cells by epithelial-lymphocyte interaction. PLoS Pathog. 2017;13:e1006247.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  23. Kranjec C, Banks L. A systematic analysis of human papillomavirus (HPV) E6 PDZ substrates identifies MAGI-1 as a major target of HPV type 16 (HPV-16) and HPV-18 whose loss accompanies disruption of tight junctions. J Virol. 2011;85:1757–64.

    CAS  PubMed  Article  Google Scholar 

  24. Nazli A, Chan O, Dobson-Belaire WN, Ouellet M, Tremblay MJ, Gray-Owen SD, et al. Exposure to HIV-1 directly impairs mucosal epithelial barrier integrity allowing microbial translocation. PLoS Pathog. 2010;6:e1000852.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  25. Nazli A, Kafka JK, Ferreira VH, Anipindi V, Mueller K, Osborne BJ, et al. HIV-1 gp120 induces TLR2- and TLR4-mediated innate immune activation in human female genital epithelium. J Immunol. 2013;191:4246–58.

    CAS  PubMed  Article  Google Scholar 

  26. Cherne MD, Cole AL, Newberry L, Schmidt-Owens M, Deichen M, Cole AM. Matrix metalloproteinases expressed in response to bacterial vaginosis disrupt the endocervical epithelium, increasing transmigration of HIV. Infect Immun. 2020;88(4):e00041–20.

    PubMed  PubMed Central  Article  Google Scholar 

  27. Graver MA, Wade JJ. The role of acidification in the inhibition of Neisseria gonorrhoeae by vaginal lactobacilli during anaerobic growth. Ann Clin Microbiol Antimicrob. 2011;10:8.

    PubMed  PubMed Central  Article  Google Scholar 

  28. Nahui Palomino RA, Zicari S, Vanpouille C, Vitali B, Margolis L. Vaginal Lactobacillus inhibits HIV-1 replication in human tissues ex vivo. Front Microbiol. 2017;8:906.

    PubMed  PubMed Central  Article  Google Scholar 

  29. Edwards VL, Smith SB, McComb EJ, Tamarelle J, Ma B, Humphrys MS, et al. The cervicovaginal microbiota-host interaction modulates Chlamydia trachomatis infection. mBio. 2019;10(4):e01548–19.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. Takada K, Komine-Aizawa S, Kuramochi T, Ito S, Trinh QD, Pham NTK, et al. Lactobacillus crispatus accelerates re-epithelialization in vaginal epithelial cell line MS74. Am J Reprod Immunol. 2018;80:e13027.

    PubMed  Article  CAS  Google Scholar 

  31. Anton L, Sierra LJ, DeVine A, Barila G, Heiser L, Brown AG, et al. Common cervicovaginal microbial supernatants alter cervical epithelial function: mechanisms by which Lactobacillus crispatus contributes to cervical health. Front Microbiol. 2018;9:2181.

    PubMed  PubMed Central  Article  Google Scholar 

  32. Lennard K, Dabee S, Barnabas SL, Havyarimana E, Blakney A, Jaumdally SZ, et al. Microbial composition predicts genital tract inflammation and persistent bacterial vaginosis in South African adolescent females. Infect Immun. 2018;86(1):e00410–7.

    CAS  PubMed  Article  Google Scholar 

  33. Rose WA 2nd, McGowin CL, Spagnuolo RA, Eaves-Pyles TD, Popov VL, Pyles RB. Commensal bacteria modulate innate immune responses of vaginal epithelial cell multilayer cultures. PLoS One. 2012;7:e32728.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. Witkin SS, Mendes-Soares H, Linhares IM, Jayaram A, Ledger WJ, Forney LJ. Influence of vaginal bacteria and D- and L-lactic acid isomers on vaginal extracellular matrix metalloproteinase inducer: implications for protection against upper genital tract infections. mBio. 2013;4(4):e00460–13.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  35. Chetwin E, Manhanzva MT, Abrahams AG, Froissart R, Gamieldien H, Jaspan H, et al. Antimicrobial and inflammatory properties of South African clinical Lactobacillus isolates and vaginal probiotics. Sci Rep. 2019;9:1917.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  36. O'Hanlon DE, Moench TR, Cone RA. Vaginal pH and microbicidal lactic acid when lactobacilli dominate the microbiota. PLoS One. 2013;8:e80074.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. O'Hanlon DE, Moench TR, Cone RA. In vaginal fluid, bacteria associated with bacterial vaginosis can be suppressed with lactic acid but not hydrogen peroxide. BMC Infect Dis. 2011;11:200.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. Aldunate M, Tyssen D, Johnson A, Zakir T, Sonza S, Moench T, et al. Vaginal concentrations of lactic acid potently inactivate HIV. J Antimicrob Chemother. 2013;68:2015–25.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  39. Tyssen D, Wang YY, Hayward JA, Agius PA, DeLong K, Aldunate M, et al. Anti-HIV-1 activity of lactic acid in human cervicovaginal fluid. mSphere. 2018;3(4):e00055–18.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. Isaacs CE, Xu W. Theaflavin-3,3’-digallate and lactic acid combinations reduce herpes simplex virus infectivity. Antimicrob Agents Chemother. 2013;57:3806–14.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. Hearps AC, Tyssen D, Srbinovski D, Bayigga L, Diaz DJD, Aldunate M, et al. Vaginal lactic acid elicits an anti-inflammatory response from human cervicovaginal epithelial cells and inhibits production of pro-inflammatory mediators associated with HIV acquisition. Mucosal Immunol. 2017;10:1480–90.

    CAS  PubMed  Article  Google Scholar 

  42. Delgado-Diaz DJ, Tyssen D, Hayward JA, Gugasyan R, Hearps AC, Tachedjian G. Distinct immune responses elicited from cervicovaginal epithelial cells by lactic acid and short chain fatty acids associated with optimal and non-optimal vaginal microbiota. Front Cell Infect Microbiol. 2020;9:446.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  43. Boskey ER, Cone RA, Whaley KJ, Moench TR. Origins of vaginal acidity: high D/L lactate ratio is consistent with bacteria being the primary source. Hum Reprod. 2001;16:1809–13.

    CAS  PubMed  Article  Google Scholar 

  44. Ilhan ZE, Laniewski P, Thomas N, Roe DJ, Chase DM, Herbst-Kralovetz MM. Deciphering the complex interplay between microbiota, HPV, inflammation and cancer through cervicovaginal metabolic profiling. EBioMedicine. 2019;44:675–90.

    PubMed  PubMed Central  Article  Google Scholar 

  45. Alisoltani A, Manhanzva MT, Potgieter M, Balle C, Bell L, Ross E, et al. Microbial function and genital inflammation in young South African women at high risk of HIV infection. Microbiome. 2020;8:165.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. Zevin AS, Xie IY, Birse K, Arnold K, Romas L, Westmacott G, et al. Microbiome composition and function drives wound-healing impairment in the female genital tract. PLoS Pathog. 2016;12:e1005889.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  47. Barnabas SL, Dabee S, Passmore JS, Jaspan HB, Lewis DA, Jaumdally SZ, et al. Converging epidemics of sexually transmitted infections and bacterial vaginosis in southern African female adolescents at risk of HIV. Int J STD AIDS. 2018;29:531–9.

    PubMed  Article  Google Scholar 

  48. Potgieter MG, Nel A, Fortuin S, Garnett S, Wendoh JM, Tabb DL, et al. MetaNovo: a probabilistic approach to peptide and polymorphism discovery in complex metaproteomic datasets. bioRxiv. 2022. https://doi.org/10.1101/605550.

  49. Afiuni-Zadeh S, Boylan KLM, Jagtap PD, Griffin TJ, Rudney JD, Peterson ML, et al. Evaluating the potential of residual Pap test fluid as a resource for the metaproteomic analysis of the cervical-vaginal microbiome. Sci Rep. 2018;8:10868.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  50. 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.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  51. Smyth GK. limma: linear models for microarray data. In: Gentleman R, Carey VJ, Huber W, Irizarry RA, Dudoit S, editors. Bioinformatics and computational biology solutions using R and bioconductor. New York: Springer; 2005.

    Google Scholar 

  52. Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/. Accessed 24 June 2020.

  53. Law CW, Chen Y, Shi W, Smyth GK. voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 2014;15:R29.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  54. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  55. Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005;21:3448–9.

    CAS  PubMed  Article  Google Scholar 

  56. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  57. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  58. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods. 2001;25:402–8.

    CAS  PubMed  Article  Google Scholar 

  59. Manhanzva MT, Abrahams AG, Gamieldien H, Froissart R, Jaspan H, Jaumdally SZ, et al. Inflammatory and antimicrobial properties differ between vaginal Lactobacillus isolates from South African women with non-optimal versus optimal microbiota. Sci Rep. 2020;10:6196.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  60. Yeswanth S, Nanda Kumar Y, Venkateswara Prasad U, Swarupa V, Koteswara rao V, et al. Cloning and characterization of l-lactate dehydrogenase gene of Staphylococcus aureus. Anaerobe. 2013;24:43–8.

    CAS  PubMed  Article  Google Scholar 

  61. Gajer P, Brotman RM, Bai G, Sakamoto J, Schutte UM, Zhong X, et al. Temporal dynamics of the human vaginal microbiota. Sci Transl Med. 2012;4:132ra152.

    Article  Google Scholar 

  62. Carias AM, Allen SA, Fought AJ, Kotnik Halavaty K, Anderson MR, Jimenez ML, et al. Increases in endogenous or exogenous progestins promote virus-target cell interactions within the non-human primate female reproductive tract. PLoS Pathog. 2016;12:e1005885.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  63. Stieh DJ, Maric D, Kelley ZL, Anderson MR, Hattaway HZ, Beilfuss BA, et al. Vaginal challenge with an SIV-based dual reporter system reveals that infection can occur throughout the upper and lower female reproductive tract. PLoS Pathog. 2014;10:e1004440.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  64. Nunn KL, Wang YY, Harit D, Humphrys MS, Ma B, Cone R, et al. Enhanced trapping of HIV-1 by human cervicovaginal mucus is associated with Lactobacillus crispatus-dominant microbiota. mBio. 2015;6:e01084–15.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  65. Carias AM, Hope TJ. Barriers of mucosal entry of HIV/SIV. Curr Immunol Rev. 2019;15:4–13.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. Dizzell S, Nazli A, Reid G, Kaushic C. Protective effect of probiotic bacteria and estrogen in preventing HIV-1-mediated impairment of epithelial barrier integrity in female genital tract. Cells. 2019;8(10):1120.

    CAS  PubMed Central  Article  Google Scholar 

  67. Cummings JH, Pomare EW, Branch WJ, Naylor CP, Macfarlane GT. Short chain fatty acids in human large intestine, portal, hepatic and venous blood. Gut. 1987;28:1221–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  68. Niklasson O, Skude G, Mardh PA. Lactate dehydrogenase and its isoenzymes in vaginal fluid in vaginitis/vaginosis cases and in healthy controls. Int J STD AIDS. 2003;14:270–3.

    CAS  PubMed  Article  Google Scholar 

  69. Boskey ER, Telsch KM, Whaley KJ, Moench TR, Cone RA. Acid production by vaginal flora in vitro is consistent with the rate and extent of vaginal acidification. Infect Immun. 1999;67:5170–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  70. Berard AR, Perner M, Mutch S, Farr Zuend C, McQueen P, Burgener AD. Understanding mucosal and microbial functionality of the female reproductive tract by metaproteomics: implications for HIV transmission. Am J Reprod Immunol. 2018;80:e12977.

    PubMed  Article  Google Scholar 

  71. Tachedjian G, Aldunate M, Bradshaw CS, Cone RA. The role of lactic acid production by probiotic Lactobacillus species in vaginal health. Res Microbiol. 2017;168:782–92.

    CAS  PubMed  Article  Google Scholar 

  72. Plummer EL, Bradshaw CS, Doyle M, Fairley CK, Murray GL, Bateson D, et al. Lactic acid-containing products for bacterial vaginosis and their impact on the vaginal microbiota: a systematic review. PLoS One. 2021;16:e0246953.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  73. Ross J, Brittain C, Kai J, David M, Ozolins M, Jackson L, et al. O03.3 Intravaginal lactic acid gel versus oral metronidazole for treating women with recurrent bacterial vaginosis: the VITA randomised controlled trial. In: STI & HIV World Congress, July 14–17 2021: 2021: sexually transmitted infections; 2021.

    Google Scholar 

  74. Sexually transmitted infections treatment guidelines, 2021. https://www.cdc.gov/std/treatment-guidelines/bv.htm. Accessed 11 May 2022.

  75. Anstey Watkins J, Ross JDC, Thandi S, Brittain C, Kai J, Griffiths F. Acceptability of and treatment preferences for recurrent bacterial vaginosis-topical lactic acid gel or oral metronidazole antibiotic: qualitative findings from the VITA trial. PLoS One. 2019;14:e0224964.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  76. Cohen CR, Wierzbicki MR, French AL, Morris S, Newmann S, Reno H, et al. Randomized trial of Lactin-V to prevent recurrence of bacterial vaginosis. N Engl J Med. 2020;382:1906–15.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

We would like to acknowledge and thank the participants of the Women’s Initiative in Health Study and the Monash University Bioinformatics platform for assistance with RNAseq analysis. We thank Michael France for uploading sequences to the NCBI Sequence Read Archive. We gratefully acknowledge the contribution to this work of the Victorian Operational Infrastructure Support Program received by the Burnet Institute.

Funding

This work was supported by Project Grant (GNT1164982) awarded to GT, ACH, JAH, JR and ADB from the National Health and Medical Research Council of Australia and NHMRC Senior Research Fellowship GNT1117748 to GT. DJD was supported by a MBio Postgraduate Discovery Scholarship from the Department of Microbiology at Monash University at Melbourne, Australia, and BJ by an Australian Government Research Training Program Scholarship from Monash University. The metaproteomics laboratory work, data analysis and writing were supported by the Carnegie Corporation of New York, South African National Research Foundation (NRF), the Poliomyelitis Research Foundation and the South African Medical Research Council (PI: LM). The WISH cohort was supported by the European and Developing Countries Clinical Trials Partnership Strategic Primer grant (SP.2011.41304.038; PI J.S. Passmore) and the South African Department of Science and Technology (DST/CON 0260/2012; PI: JSP). Research reported in this publication was supported in part by the National Institute for Allergy and Infectious Diseases of the National Institutes of Health under awards number U01AI070921, UH2AI083264, and U19AI084044. 

Author information

Authors and Affiliations

Authors

Contributions

DJD, BJ, DT and SBS generated in vitro experimental data, whilst LB, ER, SD, SB, HG and L-GB contributed data from the WISH cohort study. DJD, BJ, JAH, AA, MP, AI, IA, JMB, NM, SD, SB, VLE, LM, ACH and GT contributed to data analysis and interpretation. L-GB, JR, ADB, JSP, LM, ACH and GT were involved in study design/conception and study/site supervision. ACH wrote the manuscript with assistance from DJD, BJ, JAH, LM and GT. The authors read and approved the final manuscript.

Authors’ information

Not applicable.

Corresponding author

Correspondence to Gilda Tachedjian.

Ethics declarations

Ethics approval and consent to participate

The Universities of Witwatersrand and Cape Town Research Ethics Committees approved the recruitment of women into the WISH study. Consent for participation was obtained as previously described [47].

Consent for publication

Not applicable.

Competing interests

GT and ACH are coinventors on a patent application examining the immunomodulatory effects of L-lactic acid on cervicovaginal epithelial cells. The other authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Association between cervicovaginal pH and relative abundance of bacterial L- and D-lactate dehydrogenase [LDH; Log2 intensity based absolute quantification (iBAQ) values]. Figure S2. Multidimensional scaling (MDS) plot of genes differentially expressed by ectocervical epithelial cells treated with 0.3% L-lactic acid (L-LA) pH 3.9 or D-lactic acid (D-LA) pH 3.9 (yellow and blue respectively), HCl pH 3.9 (green) or untreated cells (red). Figure S3. Primary cervicovaginal epithelial cells were apically treated with 0.3% L-lactic acid pH (L-LA, pH 3.9), or low pH media alone (HCl, pH 3.9) for 1 h and gene expression of tight junction factors claudin-1 (CLDN1) and claudin-4 (CLDN4) were assessed 4 h post-treatment by qRT-PCR.

Additional file 2: Supplementary Table 1.

Proteins with significantly different abundance in women with high as compared to low L-LDH and D-LDH abundance*. Supplementary Table 2. Barrier-related proteins with significantly different abundance in women with high as compared to low L-LDH abundance*. Supplementary Table 3. Genes differentially expressed by L-LA, D-LA and HCL treatment of Ect cells*. Supplementary Table 4. Gene ontology pathways significantly enriched by L-LA and D-LA treatment of Ect cells. Supplementary Table 5. Tight junction genes differentially expressed by L-LA, D-LA and HCL treatment of Ect cells*. Supplementary Table 6. Tight junction genes differentially expressed by treatment of VK2 cells with bacterial culture supernatants*.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Delgado-Diaz, D.J., Jesaveluk, B., Hayward, J.A. et al. Lactic acid from vaginal microbiota enhances cervicovaginal epithelial barrier integrity by promoting tight junction protein expression. Microbiome 10, 141 (2022). https://doi.org/10.1186/s40168-022-01337-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-022-01337-5

Keywords

  • Female reproductive tract
  • Lactic acid
  • Vaginal microbiome
  • Tight junctions
  • Transcriptomics
  • Epithelial cells
  • Metabolites
  • Lactobacilli
  • STIs
  • HIV