Skip to main content

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 12 January 2022
Sec. Microbiome in Health and Disease
This article is part of the Research Topic The Female Reproductive Tract Microbiome - Gatekeeper for Sexual and Reproductive Health View all 20 articles

Vaginal Dysbiotic Microbiome in Women With No Symptoms of Genital Infections

  • 1Department of Molecular Immunology and Microbiology, Indian Council of Medical Research (ICMR)-National Institute for Research in Reproductive and Child Health, Mumbai, India
  • 2School of Applied Sciences & Technology (SAST-GTU), Gujarat Technological University, Ahmedabad, India
  • 3Department of Obstetrics and Gynecology, King Edward Memorial Hospital and Seth Gordhandas Sunderdas Medical College, Mumbai, India

The vaginal microbiome plays a critical role in determining the progression of female genital tract infections; however, little is known about the vaginal microbiota of Indian women. We aimed to investigate the vaginal microbial architecture of women with asymptomatic bacterial vaginosis (BV) (n=20) and normal microbiota (n=19). Microbial diversity was analyzed in vaginal swabs from regularly menstruating women (18-45yrs) by 16S rRNA V3-V4 amplicon (MiSeq Illumina) sequencing. Rarefaction analysis showed a higher number of species in normal flora compared to BV. Alpha diversity as measured by Pielou’s evenness revealed microbial diversity was significantly greater in BV samples than normal microbiota (p= 0.0165). Beta diversity comparison using UniFrac metrics indicated distinct microbial communities clustering between normal and BV flora. Firmicutes were the major phyla observed in vaginal specimens of normal microbiota whereas Actinobacteria, Fusobacteria, Bacteroidetes were significantly abundant in BV samples. Notably, the relative abundance of Lactobacillus was significantly high in normal microbiota. Conversely Gardnerella, Sneathia, Prevotella, Atopobium, Ureaplasma, Dialister significantly dominated dysbiotic microbiota. Relative frequency of Lactobacillus decreased significantly in BV (6%) as compared to normal microbiota (35.2%). L. fermentum, L. gasseri, L. iners, L. jensenii, L. mucosae, L. ruminis, L. salivarius, L. coleohominis was more exclusively present in normal microbiota. L. iners was detected from both the groups with a relative frequency of 50.4% and 17.2% in normal and BV microbiota respectively. Lefse analysis indicated Atopobium vaginae, Sneathia amnii, Mycoplasma hominis Prevotella disiens in the vaginal microbiota as a biomarker for dysbiosis and L. jensenii as a biomarker of a healthy microbiota. Firmicutes were negatively correlated to Tenericutes, Actinobacteria, Bacteroidetes, and Fusobacteria. Proteobacteria positively correlated to Tenericutes, and Bacteroidetes were shown to be positively correlated to Fusobacteria. Predicted functional analysis indicated differences in the functional profiles between BV and normal microbiota. Normal microbiota utilized pathways essential for phosphatidylglycerol biosynthesis I & II, peptidoglycan biosynthesis, geranylgeranyl diphosphate biosynthesis I, mevalonate pathway, CoA biosynthesis pathway I and pyrimidine nucleotide salvage; whereas BV bacteria had characteristic aromatic amino acid biosynthesis, pentose phosphate pathway, carbohydrate degradation. In conclusion, women with asymptomatic BV have vaginal microbiota significantly different than women with normal microbiota. Furthermore, the study provides insights into the vaginal microbial structure of Indian women that will enable us to explore the prospective candidates for restoring the vaginal microbiota.

Introduction

Reproductive Tract Infections (RTIs) caused by bacteria, viruses or protozoa have a profound impact on the reproductive and sexual health of men and women. The consequences of sexually transmitted infections/reproductive Tract Infections (STIs/RTIs) for reproductive health include pelvic inflammatory disease (PID), infertility (in women and men), ectopic pregnancy, and adverse pregnancy outcomes including miscarriage, stillbirth, preterm birth, and congenital infection (Mullick et al., 2005; Odogwu et al., 2021). STIs and RTIs are also known to increase the risk of HIV transmission (Taha and Gray, 2000; van de Wijgert et al., 2008; Ward and Rönn, 2010).

The healthy vaginal microbiota is predominantly dominated by lactobacilli which help in maintaining a healthy vaginal microenvironment and protecting the host from urogenital infections. Women with a healthy vaginal microbiome, are known to mostly harbor L. crispatus, L. gasseri, L. jensenii and L. iners which have been reported previously (Nunn and Forney, 2016; Pramanick et al., 2019). Based on the predominance of these species the vaginal microbiota has been classified as CST I, II, III and V, respectively (Ravel et al., 2011).

The vaginal microbiome has been studied in pregnancy (Romero et al., 2014: Mehta et al., 2020), cervical cancer (Champer et al., 2018; Chambers et al., 2021; and infertility (Zhang et al., 2021). Alteration in the vaginal ecology during which the commensal lactobacilli are reduced and displaced by the anaerobic species of the genera Gardnerella, Sneathia, Prevotella, Atopobium. can lead to bacterial vaginosis. The vaginal dysbiosis not only affects the quality of life in women but could lead to poor reproductive sequelae in IVF patients (Haahr et al., 2016; Skafte-Holm et al., 2021), adverse pregnancy outcomes and increase risk of cancer (Yang et al., 2020). This change in the microbial ecology has been attributed to various factors such as sex hormones (Brotman et al., 2014), hormonal contraception (Achilles and Hillier, 2013), sexual practices, personal hygiene and diet (Lewis et al., 2017). Significant differences in the vaginal microbiome composition related to races/ethnicity have been reported (Peipert et al., 2008; Fettweis et al., 2014). However, in 15- 84% of apparently healthy women BV could be asymptomatic (Koumans et al., 2007; Pramanick et al., 2019). While many studies have characterized the vaginal microbiota of symptomatic cases of BV, asymptomatic BV remained uncharacterized.

India has a varied population and diversified ethnic composition. The limited vaginal microbiota studies on Indian women have focused only on lactobacilli (Garg et al., 2009; Madhivanan et al., 2015; Das Purkayastha et al., 2019). These studies have often reported contrasting observations. Furthermore, the microbiome of non-pregnant women of reproductive age in India remains unexplored. In this work, we have characterized the vaginal microbiota of a healthy and asymptomatic BV microbiota.

Methods

Study Participants and Recruitment

Vaginal swabs from 39 married, non-pregnant regularly menstruating women of the reproductive age group (18-45yrs) were collected for the study. The participants were recruited from the Gynecology and Obstetrics Out-Patients Department of King Edward Memorial Hospital and Seth Gordhandas Sunderdas Medical College, Mumbai. The study had the approval of the human ethics review board at the institute (Protocol Number 215/2012) and the collaborating institute (Protocol No EC/GOV-5/2012). Informed consent was obtained from all the participants. Inclusion criteria were women who were clinically healthy and in the age group of 18 - 45 years, had maintained sexual abstinence of at least 5 days, not using any vaginal hygienic products or lubricants, willing to get screened for STIs and other infections of the lower genital tract (vagina and cervix), have been diagnosed as normal with no infections of the lower genital tract based on the laboratory evidence, had not taken any antibiotics in the last six weeks, not menstruating and preferably 5 days following LMP, did not have any dysfunctional uterine bleeding or genital malignancy or genital prolapse. Pregnant women, women using hormonal contraceptives or OC’s, women with intrauterine devices, chronic drug intake, or with chronic medical illness were excluded from the study.

Nugent Scoring and Amsel Criteria

One swab was used to assess vaginal health by Amsel’s criteria (Amsel et al., 1983) and Nugent scoring (Nugent et al., 1991). The second swab was used for DNA extraction for sequencing. Nugent scoring was based on the microscopic examination of different vaginal bacterial morphotypes and their enumeration on gram stained vaginal smears. Normal microbiota (score 0–3), intermediate (score 4–6), or BV microbiota (score 7–10) were the scores assigned by Nugent scoring. Amsel criteria were based on clinical symptoms and signs that include vaginal pH>4.5, the presence of clue cells (vaginal epithelial cells laden with coccobacilli) using wet mount microscopy, homogeneous white vaginal discharge and fishy odor of discharge (10% KOH amine test). A patient who satisfied three of these four criteria was diagnosed as positive for BV.

Sample Processing and DNA Extraction

Genomic DNA was extracted according to the instruction on MO BIO Powersoil DNA extraction kit (Qiagen, Catalog No. 12888). The extracted DNA was quantified and purity was determined on the NanoDrop reader. The samples were stored at -200C until further use.

Amplicon Sequencing

The amplicons libraries were prepared using the Nextera XT Index kit (Illumina Inc.) as per the 16S metagenomic sequencing library preparation protocol. Primers for the amplification of the V3-V4 region of 16S rDNA were GCCTACGGGNGGCWGCAG and reverse ACTACHVGGGTATCTAATCC. The amplicon libraries were purified by AMPureXP beads and quantified using a Qubit fluorometer. The libraries were sequenced on MiSeq using 2x300 bp chemistry.

Data Analysis and Statistics

Sequence read quality was assessed using Fastqc. Based on the observed quality, the reads were filtered out with the following parameters: all forward/reverse reads were trimmed by 10 bases from the left; from the right end, the reads were trimmed at length 280 for forward and length 220 for reverse. Reads were further pre-processed with DADA2 for denoising and eliminating chimera sequences and duplicates (Supplementary Table 1). The predicted ASVs were normalized by the minimum number of feature sequences in a sample from each of the study groups, respectively. Read pre-processing and taxonomic classification were performed in QIIME2 framework (Bolyen et al., 2019) using the pre-trained classifier for V3-V4 region of 16S rRNA genes of the SILVA database. Beta diversity was calculated with unweighted and weighted UniFrac metrics to compute the distribution of the samples of studied groups, viz., BV and Normal in Principal Coordinates Analysis (PCoA). Alpha diversity was estimated using the number of observed amplicon sequence variants (ASVs) using Pielou’s evenness index. The microbial taxonomy was studied for statistical differences (p < 0.05) between the studied groups and the core microbiota was compared using Microbiome Analyst (Dhariwal et al., 2017). Functional attributes corresponding to the observed microbiota were assessed using the q2-PICRUST plugin (Douglas et al., 2020). The statistical differences between the normal and BV were studied by performing Kruskal–Wallis test in the STAMP v2.1 software (Parks et al., 2014).

Linear discriminant analysis (LDA) Effect Size (LEfSe) method (Segata et al., 2011) was employed to predict biomarker species and functional attributes (LDA score cutoff value of 4) with significant differences in abundance between the two groups. Linear discriminant analysis (LDA) Effect Size (LEfSe) method was employed to predict biomarker species and functional attributes (LDA score cutoff value of 4) with significant differences in abundance between the two groups.

Data Availability

The sequences were submitted in the NCBI database under the Bio project Accession No.: PRJNA674451. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA674451.

Results

For this study, we recruited 39 regularly menstruating women of the reproductive age group of 18-45yrs, who were not on antibiotics. Nugent scoring of the vaginal swabs revealed, 19 women had normal microbiota and 20 had BV. The recruited participant’s characteristics are summarized in Table 1. There was no significant difference in age and the menstrual phases of participants recruited for each group (Table 1). Due to its cosmopolitan nature the population of Mumbai is characterized by people from different parts of the country. About 15.38% (6) participants were nonresidents of Mumbai whereas 10.26% (4) have migrated in the last one year from sample collection.

TABLE 1
www.frontiersin.org

Table 1 Characteristics of participants recruited for the study.

The vaginal microbial profile of these women was characterized using Illumina MiSeq sequencing of the V3-V4 region of the 16s rRNA gene with an average read length of 301bp. A total of 7041,312 reads were generated from 39 swabs with an average of 352,006 sequences per sample. On average 92.38% of sequences had a Phred score of >20.

Comparison of Vaginal Microbiota Between Healthy and Asymptomatic BV Samples

The rarefaction analysis revealed a higher number of observed species in Normal microbiota compared to BV samples. However, the comparison difference revealed an FDR Adjusted P-value = 0.19 showing a less significant difference between Normal and BV samples. The rarefaction curve generated by the OTUs obtained from both groups indicated sufficient sequencing depth (Supplementary Figure 1). The alpha diversity measured by the Pielou evennessindex revealed women with normal microbiota had less diverse communities as compared to women with asymptomatic BV. As observed from Figure 1A, a significant (p=0.0165) difference exists between the community evenness of BV and Normal samples. Beta diversity measures the amount of dissimilarity between the samples. Beta diversity comparison using UniFrac metrics, used to measure differences in microbial abundances between the groups, demonstrated distinct microbial communities between normal and BV samples (p = 0.0127). The first two components of the PCoA explained 78.32% of the variations among the samples in axes 1 (49.95) and 2 (28.41%) respectively. The cluster corresponding to the normal microbiota showed variation along axis 1, from the BV population (Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1 (A) Alpha-diversity represented by the Pielou evenness index in the BV and normal samples indicating significant difference (p-value = 0.0165) between the groups based on the Kruskal-Wallis test. The boxes represent the distributions of the alpha diversity index and show the median for each condition. Whiskers extend to the furthest data point. (B) PCoA generated using weighted Unifrac distance indicates distinct clustering of samples from each group. Each point corresponds to an individual sample. For each experimental group, an ellipse around the centroid is depicted. The first two components of the variance are represented by Plotting BV(blue) vs Normal (Red) samples with significant separation between the two groups. The first two components represent 78.32 percent of variance, individually depicted in parentheses next to Axis1 and Axis2.

Firmicutes were the most abundant phyla in normal microbiota whereas Bacteroidetes, Actinobacteria, Fusobacteria were significantly abundant in the asymptomatic BV group (Figure 2A). Lactobacillus was the predominant genus in normal microbiota whereas Gardnerella, Sneathia, Prevotella, Atopobium, significantly dominated the BV microbiota (Figure 2B). The heatmap in Figure 3 depicts the abundance of the major genus present in each sample. The mean proportion of Lactobacillus and Gardnerella in normal microbiota was 87.2% and 3.1% as compared to 24% and 28.2% in BV samples respectively. Sneathia and Atopobium were absent in normal microbiota but constituted 12.1% and 4.9% in asymptomatic BV samples.

FIGURE 2
www.frontiersin.org

Figure 2 Relative abundance of bacteria at the (A) phylum and (B) genus level in the normal and BV groups as computed by Welch’s t-test using STAMP software. The middle shows the difference between proportions of abundance in the 95% confidence interval, and the value at the right is the P-value. P < 0.05 represents the significant difference between the two groups.

FIGURE 3
www.frontiersin.org

Figure 3 Genus level composition of the vaginal microbiota in normal and BV samples. The heatmap depicts the relative abundances of the most abundant genera. Each column represents one sample. Orange are samples with normal microbiota and Blue are BV samples.

Lactobacillus was detected in all the samples of normal microbiota whereas it was completely missing from two samples with asymptomatic BV. L. iners (84.21%), (65%) and unidentified Lactobacillus spps. (84.21%), (80%) were the most frequently detected lactobacilli in both normal and BV samples respectively. Additionally, most of the women with normal microbiota harbored L. jensenii (36.84%) and L. gasseri (31.57%). Women with asymptomatic BV were frequently detected with L. gasseri (15%), L. salivarius (15%) and L. fermentum (15%). Other lactobacilli exclusively present in normal microbiota were L. coleohominis (10.52%) and L. ruminis (5.26%).

Classification of Microbiota Based on Bacterial Markers

Taxonomic biomarker identification in each group was carried out using LefSe analysis. Discriminate analysis using LefSe showed Lactobacillus jensenii, Comamonas, Weissella were enriched in normal samples while Atopobium vaginae, Sneathia amnii, Mycoplasma hominis were enriched in BV samples (Supplementary Figure 2).

Moreover, Cladogram represents differentially abundant genus and species (Figure 4). Lactobacillus jensenii was significantly abundant in normal microbiota as compared to asymptomatic BV. On the other hand, species such as Atopobium vaginae, Coriobacteriales bacterium DNF00809, Sneathia amnii were differentially elevated in asymptomatic BV samples.

FIGURE 4
www.frontiersin.org

Figure 4 Cladogram representation of differentially abundant bacterial families detected using LEfSe. The cladogram diagram shows the microbial species with significant differences in the two groups. Red and green, indicate BV and normal microbiota groups respectively, with the species classification at the level of phylum, class, order, family, and genus shown from the outside to the inside. The red and green nodes in the phylogenetic tree represent microbial species that play an important role in the BV and normal microbiota groups, respectively. Yellow nodes represent species with no significant difference. Significantly abundant bacterial groups observed in the study are shown in the list on the right hand side.

Correlation of Microbial Members of the Vaginal Ecology

Firmicutes were negatively correlated to Tenericutes, Actinobacteria, Bacteroidetes and Fusobacteria. Proteobacteria positively correlated to Tenericutes, and Bacteroidetes were shown to be positively correlated to Fusobacteria (Figure 5A). Lactobacillus highly negatively correlated to Gardnerella and other BV-related bacteria, but positively correlated to Pseudomonas. Sneathia positively correlated with Enterococcus, Mycoplasma and other pathogens such as Dialister, Prevotella, Atopobium, Streptococcus (Figure 5B) (Supplementary Table 2).

FIGURE 5
www.frontiersin.org

Figure 5 Bacterial correlation (A) at the phylum level. Positive correlations are represented as red and negative correlations as blue (B) Network of associations between the bacteria at the genus level. Red are samples with normal microbiota and green are BV samples. The nodes represent genera of bacteria; the edges represent the correlation coefficients between genera. Node size indicates the number of subjects in which the genus was seen.

Functional Profile of Vaginal Microbiota

The normal microbiota was enriched in pathways involved in phosphatidylglycerol biosynthesis I & II, peptidoglycan biosynthesis, geranylgeranyl diphosphate biosynthesis I, mevalonate pathway, CoA biosynthesis pathway I and pyrimidine nucleotide salvage while the microbiota in BV women was enriched with aromatic amino acid biosynthesis, pentose phosphate pathway, carbohydrate degradation (Figure 6 and Table 2).

FIGURE 6
www.frontiersin.org

Figure 6 Predicted functional profiling: Differential pathway abundance predicted by PICRUSt between normal and BV group. Pathways (BioCyc IDs) with an LDA score of 4 and above in normal (Green) and BV (Red) samples are represented.

TABLE 2
www.frontiersin.org

Table 2 The enriched pathways in vaginal microbiota of Indian women.

Discussion

Lactobacillus which is the keystone bacteria of vaginal microbiota was significantly present in the normal sample as compared to asymptomatic BV microbiota. In contrast to our observation, a study on Estonia women of reproductive age has reported lactobacilli dominance in healthy as well as asymptomatic BV women (Drell et al., 2013). Lactobacillus, Comamonas, Weissella were identified as the biomarkers for a healthy microbiota. Lactobacillus, Leuconostoc, Weissella, and Streptococcus have been commonly identified in vaginal samples of other populations (Jin et al., 2007). Along with Lactobacillus Weissella spp. have been isolated in fermented food items and feces (Ruiz-Moyano et al., 2011; Gomathi et al., 2014; Zhang et al., 2014; Kang et al., 2020; Pabari et al., 2020) and demonstrated to have probiotic and anti-inflammatory properties (Fatmawati et al, 2020). Weissella is another lactic acid bacterium known to produce H2O2 in the genital tract of women (Jin et al., 2007) and shown to be a potential probiotic for vaginal health (Lee, 2005). Atopobium, Gardnerella, Sneathia were the biomarker genus for asymptomatic BV. The taxonomic biomarker identified in BV samples was similar to previously reported in Chinese women (Ling et al., 2010). All the enriched genera in asymptomatic BV samples have been reported as vaginal pathogens in various studies and contribute to the sequelae of spontaneous abortions, preterm birth, infertility and cancer (Agarwal et al., 2020). The presence of dysbiotic microbiota in vagina of asymptomatic women addresses the need for the characterization of the microbiome in those women who may have spontaneous abortions or infertility problems or those who go in for in vitro fertilization and embryo transfer (IVF-ET) procedures (Koedooder et al., 2019; Kong et al., 2020; Wang et al., 2021). Our previous study on cultivable microbiota also reported a distinct microbial diversity of asymptomatic BV from the normal samples (Pramanick et al., 2019). The translational aspects of the vaginal microbiome and metabolome data could be exploited and based on the identification of specific biomarkers from BV and normal samples, new diagnostic point of care tests and treatment modalities could be developed.

The predicted key functional pathways in the normal microbiota were peptidoglycan biosynthesis, phosphatidylglycerol biosynthesis I and II, peptidoglycan biosynthesis IV and pyrimidine nucleotide salvage pathways. Increased cell wall organization and peptidoglycan biosynthesis in Lactobacillus dominated microbiomes have been associated with reduced FGT inflammation (Alisoltani et al., 2020) and in modulating the immune system (Song et al., 2018). Peptidoglycan biosynthesis IV and pyrimidine nucleotide salvage pathways have been further described during the proliferative phase of menstrual cycle and associated with increased bacterial proliferation (Chen et al., 2017). Our studies show that normal microbiota dominant in lactobacillus species showed mevalonate pathways. These strains are demonstrated to alleviate hyperlipidemia by modulating AMPK and downregulating cholesterol biosynthesis via the mevalonate pathway and Bloch pathway (Lew et al., 2020) and also increase the resilience of tissue cells to cholesterol-dependent cytolysins (Griffin et al., 2017). Geranylgeranyl diphosphate biosynthesis I (via mevalonate) forms geranylgeranyl diphosphate which is a crucial compound involved in the biosynthesis of a variety of terpenes and terpenoids, including central compounds such as ubiquinones and menaquinones. In addition, GGPP is also used in posttranslational modifications of proteins (geranylgeranylation) which is important for membrane adhesion as well as the function of some proteins (Jiang et al., 2017). Aromatic amino acid biosynthesis, carbohydrate degradation, pentose phosphate pathway that may cause preterm birth were enriched in asymptomatic BV women. Women who gave preterm birth had a vaginal microbiome enriched with the nonoxidative branch of the pentose phosphate pathway (Odogwu et al., 2021). Thus, functionally, many of the enriched BV pathways were also reflective of vaginal microbiota associated with preterm birth. Our findings of functionally enhanced pentose phosphate pathway in asymptomatic BV women have implications of the possibility of preterm deliveries in these women if not treated. Though the need to treat or not treat BV has been controversial our study supports that women with asymptomatic BV need to be treated to prevent adverse pregnancy outcomes.

From this study, we were able to establish the vaginal microbiome of non-pregnant Indian women for the first time. Our study shows that women with no genital symptoms of bacterial vaginosis may have a dysbiotic microbiome. The prevalence of asymptomatic BV in otherwise healthy women and functional pathways that may cause preterm birth highlights the need for further research on its pathogenesis. The presence of these bacteria in women with asymptomatic BV is worrisome. Hence it is imperative to maintain homeostasis to prevent any future episodes of infections.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA674451.

Ethics Statement

The studies involving human participants were reviewed and approved by human ethics review board at the ICMR-National Institute for Research in Reproductive and Child Health (Protocol Number 215/2012) and the collaborating institute King Edward Hospital and Seth Gordhandas Sunderdas Medical CollegeGor(Protocol No EC/GOV-5/2012). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

RP and CA designed the study. NM and HW recruited the participants and carried out sample collection. RP carried out sample collection and its processing. NN performed the bioinformatic analysis. RP and CA interpreted the data. RP and CA wrote the paper. CA was involved in the acquisition of funding and review of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Indian Council of Medical Research (Grant number 5/7/726/2012-RHN).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We express our gratitude to the Director (ICMR-NIRRCH) for constant guidance and encouragement. The laboratory is funded by grants from Indian Council of Medical Research (ICMR), Govt. of India. RP is thankful to ICMR JRF on the project. The manuscript bears the NIRRCH ID.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2021.760459/full#supplementary-material

Supplementary Figure 1 | Rarefaction curves showing the number of observed species for samples belonging to BV and Normal condition at each rarefaction depth starting from 100 seqs/sample to 80000 seqs/sample

Supplementary Figure 2 | Differentially abundant bacterial taxa between normal and BV microbiota by LEfSe analyses Significant bacterial genera were determined by the Kruskal-Wallis test (P<0.05) with LDA score greater than 2. Normal and BV microbiota are represented by green and red bars respectively.

Supplementary Table 1 | Sequence statistics of the samples after deionizing by DADA2

Supplementary Table 2 | Pearson correlation coefficients and P-values between bacterial taxa.

References

Achilles, S. L., Hillier, S. L. (2013). The Complexity of Contraceptives: Understanding Their Impact on Genital Immune Cells and Vaginal Microbiota. AIDS 27 Suppl 1 (01), S5–s15. doi: 10.1097/QAD.0000000000000058

PubMed Abstract | CrossRef Full Text | Google Scholar

Agarwal, K., Robinson, L. S., Aggarwal, S., Foster, L. R., Hernandez-Leyva, A., Lin, H., et al. (2020). Glycan Cross-Feeding Supports Mutualism Between Fusobacterium and the Vaginal Microbiota. PloS Biol. 18 (8), e3000788. doi: 10.1371/journal.pbio.3000788

PubMed Abstract | CrossRef Full Text | Google Scholar

Alisoltani, A., Manhanzva, M. T., Potgieter, M., Balle, C., Bell, L., Ross, E., et al. (2020). Microbial Function and Genital Inflammation in Young South African Women at High Risk of HIV Infection. Microbiome 8 (1), 165. doi: 10.1186/s40168-020-00932-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Amsel, R., Totten, P. A., Spiegel, C. A., Chen, K. C., Eschenbach, D., Holmes, K. K. (1983). Nonspecific Vaginitis. Diagnostic Criteria and Microbial and Epidemiologic Associations. Am. J. Med. 74 (1), 14–22. doi: 10.1016/0002-9343(83)91112-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, Interactive, Scalable and Extensible Microbiome Data Science Using QIIME 2. Nat. Biotechnol. 37 (8), 852–857. doi: 10.1038/s41587-019-0209-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Brotman, R. M., Ravel, J., Bavoil, P. M., Gravitt, P. E., Ghanem, K. G. (2014). Microbiome, Sex Hormones, and Immune Responses in the Reproductive Tract: Challenges for Vaccine Development Against Sexually Transmitted Infections. Vaccine 32 (14), 1543–1552. doi: 10.1016/j.vaccine.2013.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Chambers, L. M., Bussies, P., Vargas, R., Esakov, E., Tewari, S., Reizes, O., et al. (2021). The Microbiome and Gynecologic Cancer: Current Evidence and Future Opportunities. Curr. Oncol. Rep. 23 (8), 92. doi: 10.1007/s11912-021-01079-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Champer, M., Wong, A. M., Champer, J., Brito, I. L., Messer, P. W., Hou, J. Y., et al. (2018). The Role of the Vaginal Microbiome in Gynaecological Cancer. BJOG: Int. J. Obstetrics Gynaecol. 125 (3), 309–315. doi: 10.1111/1471-0528.14631

CrossRef Full Text | Google Scholar

Chen, C., Song, X., Wei, W., Zhong, H., Dai, J., Lan, Z., et al. (2017). The Microbiota Continuum Along the Female Reproductive Tract and its Relation to Uterine-Related Diseases. Nat. Commun. 8(1), 875. doi: 10.1038/s41467-017-00901-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Das Purkayastha, S., Bhattacharya, M. K., Prasad, H. K., Upadhyaya, H., Lala, S. D., Pal, K., et al. (2019). Contrasting Diversity of Vaginal Lactobacilli Among the Females of Northeast India. BMC Microbiol. 19 (1), 198. doi: 10.1186/s12866-019-1568-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Dhariwal, A., Chong, J., Habib, S., King, I. L., Agellon, L. B., Xia, J. (2017). MicrobiomeAnalyst: A Web-Based Tool for Comprehensive Statistical, Visual and Meta-Analysis of Microbiome Data. Nucleic Acids Res. 45 (W1), W180–W188. doi: 10.1093/nar/gkx295

PubMed Abstract | CrossRef Full Text | Google Scholar

Douglas, G. M., Maffei, V. J., Zaneveld, J. R., Yurgel, S. N., Brown, J. R., Taylor, C., et al. (2020). PICRUSt2 for Prediction of Metagenome Functions. Nat. Biotechnol. 38(6), 685–688. doi: 10.1038/s41587-020-0548-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Drell, T., Lillsaar, T., Tummeleht, L., Simm, J., Aaspõllu, A., Väin, E., et al. (2013). Characterization of the Vaginal Micro- and Mycobiome in Asymptomatic Reproductive-Age Estonian Women. PloS One 8 (1), e54379. doi: 10.1371/journal.pone.0054379

PubMed Abstract | CrossRef Full Text | Google Scholar

Fatmawati, N., Gotoh, K., Mayura, I., Nocianitri, K. A., Suwardana, G., Komalasari, N., et al. (2020). Enhancement of Intestinal Epithelial Barrier Function by Weissella Confusa F213 and Lactobacillus Rhamnosus FBB81 Probiotic Candidates in an In Vitro Model of Hydrogen Peroxide-Induced Inflammatory Bowel Disease. BMC Res. Notes 13 (1), 489. doi: 10.1186/s13104-020-05338-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Fettweis, J. M., Brooks, J. P., Serrano, M. G., Sheth, N. U., Girerd, P. H., Edwards, D. J., et al. (2014). Differences in Vaginal Microbiome in African American Women Versus Women of European Ancestry. Microbiol. (Reading) 160 (Pt 10), 2272–2282. doi: 10.1099/mic.0.081034-0

CrossRef Full Text | Google Scholar

Garg, K. B., Ganguli, I., Das, R., Talwar, G. P. (2009). Spectrum of Lactobacillus Species Present in Healthy Vagina of Indian Women. Ind. J. Med. Res. 129, 652–657.

Google Scholar

Gomathi, S., Sasikumar, P., Anbazhagan, K., Sasikumar, S., Kavitha, M., Selvi, M. S., et al. (2014). Screening of Indigenous Oxalate Degrading Lactic Acid Bacteria From Human Faeces and South Indian Fermented Foods: Assessment of Probiotic Potential. Sci. World J. 2014, 1–11. doi: 10.1155/2014/648059

CrossRef Full Text | Google Scholar

Griffin, S., Preta, G., Sheldon, I. M. (2017). Inhibiting Mevalonate Pathway Enzymes Increases Stromal Cell Resilience to a Cholesterol-Dependent Cytolysin. Sci. Rep. 7, 17050. doi: 10.1038/s41598-017-17138-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Haahr, T., Jensen, J. S., Thomsen, L., Duus, L., Rygaard, K., Humaidan, P. (2016). Abnormal Vaginal Microbiota may be Associated With Poor Reproductive Outcomes: A Prospective Study in IVF Patients. Hum. Reprod. 31 (4), 795–803. doi: 10.1093/humrep/dew026

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, C., Diao, F., Sang, Y.-J., Xu, N., Zhu, R.-L., Wang, X.-X., et al. (2017). GGPP-Mediated Protein Geranylgeranylation in Oocyte Is Essential for the Establishment of Oocyte-Granulosa Cell Communication and Primary-Secondary Follicle Transition in Mouse Ovary. PloS Genet. 13 (1), e1006535. doi: 10.1371/journal.pgen.1006535

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, L., Tao, L., Pavlova, S. I., So, J.-S., Kiwanuka, N., Namukwaya, Z., et al. (2007). Species Diversity and Relative Abundance of Vaginal Lactic Acid Bacteria From Women in Uganda and Korea. J. Appl. Microbiol. 102, 1107–1115. doi: 10.1111/j.1365-2672.2006.03147.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, M. S., Lee, D. S., Lee, S. A., Kim, M. S., Nam, S. H. (2020). Effects of Probiotic Bacterium Weissella Cibaria CMU on Periodontal Health and Microbiota: A Randomised, Double-Blind, Placebo-Controlled Trial. BMC Oral. Health 20 (1), 243. doi: 10.1186/s12903-020-01231-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Koedooder, R., Singer, M., Schoenmakers, S., Savelkoul, P., Morré, S. A., de Jonge, J. D., et al. (2019). The Vaginal Microbiome as a Predictor for Outcome of In Vitro Fertilization With or Without Intracytoplasmic Sperm Injection: A Prospective Study. Hum. Reprod. (Oxford England) 34 (6), 1042–1054. doi: 10.1093/humrep/dez065

CrossRef Full Text | Google Scholar

Kong, Y., Liu, Z., Shang, Q., Gao, Y., Li, X., Zheng, C., et al. (2020). The Disordered Vaginal Microbiota Is a Potential Indicator for a Higher Failure of In Vitro Fertilization. Front. Med. 7, 217. doi: 10.3389/fmed.2020.00217

CrossRef Full Text | Google Scholar

Koumans, E. H., Sternberg, M., Bruce, C., McQuillan, G., Kendrick, J., Sutton, M., et al. (2007). The Prevalence of Bacterial Vaginosis in the United State-2004; Associations With Symptoms, Sexual Behaviors, and Reproductive Health. Sexually Transmitted Dis. 34 (11), 864–869. doi: 10.1097/OLQ.0b013e318074e565

CrossRef Full Text | Google Scholar

Lee, Y. (2005). Characterization of Weissella Kimchii PL9023 as a Potential Probiotic for Women. FEMS Microbiol. Lett. 250 (1), 157–162. doi: 10.1016/j.femsle.2005.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Lew, L. C., Hor, Y. Y., Jaafar, M. H., Lau, A. S., Lee, B. K., Chuah, L. O., et al. (2020). Lactobacillus Strains Alleviated Hyperlipidemia and Liver Steatosis in Aging Rats via Activation of AMPK. Int. J. Mol. Sci. 21 (16):5872. doi: 10.3390/ijms21165872

CrossRef Full Text | Google Scholar

Lewis, F., Bernstein, K. T., Aral, S. O. (2017). Vaginal Microbiome and Its Relationship to Behavior, Sexual Health, and Sexually Transmitted Diseases. Obstetrics Gynecol. 129 (4), 643–654. doi: 10.1097/AOG.0000000000001932

CrossRef Full Text | Google Scholar

Ling, Z., Kong, J., Liu, F., Zhu, H., Chen, X., Wang, Y., et al. (2010). Molecular Analysis of the Diversity of Vaginal Microbiota Associated With Bacterial Vaginosis. BMC Genomics 11, 488. doi: 10.1186/1471-2164-11-488

PubMed Abstract | CrossRef Full Text | Google Scholar

Madhivanan, P., Alleyn, H. N., Raphael, E., Krupp, K., Ravi, K., Nebhrajani, R., et al. (2015). Identification of Culturable Vaginal Lactobacillus Species Among Reproductive Age Women in Mysore, India. J. Med. Microbiol. 64 (6), 636–641. doi: 10.1099/jmm.0.000070

PubMed Abstract | CrossRef Full Text | Google Scholar

Mehta, O., Ghosh, T. S., Kothidar, A., Gowtham, M., Mitra, R., Kshetrapal, R., et al. (2020). Vaginal Microbiome of Pregnant Indian Women: Insights Into the Genome of Dominant Lactobacillus Species. Microb. Ecol. 80(2), 487–499. doi: 10.1007/s00248-020-01501-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mullick, S., Watson-Jones, D., Beksinska, M., Mabey, D. (2005). Sexually Transmitted Infections in Pregnancy: Prevalence, Impact on Pregnancy Outcomes, and Approach to Treatment in Developing Countries. Sex Transm. Infect. 81, 294–302. doi: 10.1136/sti.2002.004077

PubMed Abstract | CrossRef Full Text | Google Scholar

Nugent, R. P., Krohn, M. A., Hillier, S. L. (1991). Reliability of Diagnosing Bacterial Vaginosis is Improved by Standardized Method of Gram Stain Interpretation. J. Clin. Microbiol. 29, 297–301. doi: 10.1128/jcm.29.2.297-301.1991

PubMed Abstract | CrossRef Full Text | Google Scholar

Nunn, K. L., Forney, L. J. (2016). Unraveling the Dynamics of the Human Vaginal Microbiome. Yale J. Biol. Med. 89 (3), 331–337.

PubMed Abstract | Google Scholar

Odogwu, N. M., Chen, J., Onebunne, C. A., Jeraldo, P., Yang, L., Johnson, S., et al. (2021). Predominance of Atopobium Vaginae at Midtrimester: A Potential Indicator of Preterm Birth Risk in a Nigerian Cohort. mSphere 6 (1), e01261–e01220. doi: 10.1128/mSphere.01261-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Pabari, K., Pithva, S., Kothari, C., Purama, R. K., Kondepudi, K. K., Vyas, B., et al. (2020). Evaluation of Probiotic Properties and Prebiotic Utilization Potential of Weissella Paramesenteroides Isolated From Fruits. Probiotics Antimicrobial Proteins 12 (3), 1126–1138. doi: 10.1007/s12602-019-09630-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Parks, D. H., Tyson, G. W., Hugenholtz, P., Beiko, R. G. (2014). STAMP: Statistical Analysis of Taxonomic and Functional Profiles. Bioinf. (Oxford England) 30 (21), 3123–3124. doi: 10.1093/bioinformatics/btu494

CrossRef Full Text | Google Scholar

Peipert, J. F., Lapane, K. L., Allsworth, J. E., Redding, C. A., Blume, J. D., Stein, M. D. (2008). Bacterial Vaginosis, Race, and Sexually Transmitted Infections: Does Race Modify the Association? Sex Transm. Dis. 35 (4), 363–367. doi: 10.1097/OLQ.0b013e31815e4179

PubMed Abstract | CrossRef Full Text | Google Scholar

Pramanick, R., Mayadeo, N., Warke, H., Begum, S., Aich, P., Aranha, C. (2019). Vaginal Microbiota of Asymptomatic Bacterial Vaginosis and Vulvovaginal Candidiasis: Are They Different From Normal Microbiota? Microbial Pathogenesis 134, 103599. doi: 10.1016/j.micpath.2019.103599

PubMed Abstract | CrossRef Full Text | Google Scholar

Ravel, J., Gajer, P., Abdo, Z., Schneider, G. M., Koenig, S. S., McCulle, S. L., et al. (2011). Vaginal Microbiome of Reproductive-Age Women. Proc. Natl. Acad. Sci. U. S. A. 108(Suppl 1), 4680–4687. doi: 10.1073/pnas.1002611107

PubMed Abstract | CrossRef Full Text | Google Scholar

Romero, R., Hassan, S. S., Gajer, P., Tarca, A. L., Fadrosh, D. W., Nikita, L., et al. (2014). The Composition and Stability of the Vaginal Microbiota of Normal Pregnant Women is Different From That of non-Pregnant Women. Microbiome 2 (1), 4. doi: 10.1186/2049-2618-2-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruiz-Moyano, S., Martín, A., Benito, M. J., Hernández, A., Casquete, R., de Guia Córdoba, M. (2011). Application of Lactobacillus Fermentum HL57 and Pediococcus Acidilactici SP979 as Potential Probiotics in the Manufacture of Traditional Iberian Dry-Fermented Sausages. Food Microbiol. 28 (5), 839–847. doi: 10.1016/j.fm.2011.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic Biomarker Discovery and Explanation. Genome Biol. 12 (6), R60. doi: 10.1186/gb-2011-12-6-r60

PubMed Abstract | CrossRef Full Text | Google Scholar

Skafte-Holm, A., Humaidan, P., Bernabeu, A., Lledo, B., Jensen, J. S., Haahr, T. (2021). The Association Between Vaginal Dysbiosis and Reproductive Outcomes in Sub-Fertile Women Undergoing IVF-Treatment: A Systematic PRISMA Review and Meta-Analysis. Pathogens 10 (3), 295. doi: 10.3390/pathogens10030295

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, J., Lang, F., Zhao, N., Guo, Y., Zhang, H. (2018). Vaginal Lactobacilli Induce Differentiation of Monocytic Precursors Toward Langerhans-Like Cells: In Vitro Evidence. Front. Immunol. 9, 2437. doi: 10.3389/fimmu.2018.02437

PubMed Abstract | CrossRef Full Text | Google Scholar

Taha, T. E., Gray, R. (2000). Genital Tract Infections and Perinatal Transmission of HIV. Ann. N. Y. Acad. Sci. 918, 84–98. doi: 10.1111/j.1749-6632.2000.tb05477.x

PubMed Abstract | CrossRef Full Text | Google Scholar

van de Wijgert, J. H., Morrison, C. S., Cornelisse, P. G., Munjoma, M., Moncada, J., Awio, P., et al. (2008). Bacterial Vaginosis and Vaginal Yeast, But Not Vaginal Cleansing, Increase HIV-1 Acquisition in African Women. J. Acquir. Immune Defic. Syndr. 48 (2), 203–210. doi: 10.1097/QAI.0b013e3181743936

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., Zhou, G., Wu, L., Huang, X., Li, Y., Luo, B., et al. (2021). The Microbial Composition of Lower Genital Tract May Affect the Outcome of In Vitro Fertilization-Embryo Transfer. Front. Microbiol. 12, 729744. doi: 10.3389/fmicb.2021.729744

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, H., Rönn, M. (2010). Contribution of Sexually Transmitted Infections to the Sexual Transmission of HIV. Curr. Opin. HIV AIDS 5 (4), 305–310. doi: 10.1097/COH.0b013e32833a8844

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Q., Wang, Y., Wei, X., Zhu, J., Wang, X., Xie, X., et al. (2020). The Alterations of Vaginal Microbiome in HPV16 Infection as Identified by Shotgun Metagenomic Sequencing Front. Cell. Infect. Microbiol. 23, 1–14. doi: 10.3389/fcimb.2020.00286

CrossRef Full Text | Google Scholar

Zhang, F., Dai, J., Chen, T. (2021). Role of Lactobacillus in Female Infertility Via Modulating Sperm Agglutination and Immobilization. Front. Cell Infect. Microbiol. 10, 620529. doi: 10.3389/fcimb.2020.620529

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Peng, X., Zhang, N., Liu, L., Wang, Y., Ou, S. (2014). Cytotoxicity Comparison of Quercetin and its Metabolites From In Vitro Fermentation of Several Gut Bacteria. Food Funct. 5, 2152–2156. doi: 10.1039/C4FO00418C

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: microbiome & dysbiosis, vaginal microbiota (VMB), asymptomatic BV, Lactobacillus, bacterial vaginosis (BV)

Citation: Pramanick R, Nathani N, Warke H, Mayadeo N and Aranha C (2022) Vaginal Dysbiotic Microbiome in Women With No Symptoms of Genital Infections. Front. Cell. Infect. Microbiol. 11:760459. doi: 10.3389/fcimb.2021.760459

Received: 18 August 2021; Accepted: 14 December 2021;
Published: 12 January 2022.

Edited by:

Lindi Masson, Burnet Institute, Australia

Reviewed by:

Frank Schuren, Netherlands Organisation for Applied Scientific Research, Netherlands
Michele Tanturli, Università degli Studi di Firenze, Italy

Copyright © 2022 Pramanick, Nathani, Warke, Mayadeo and Aranha. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Clara Aranha, YXJhbmhhY0BuaXJyaC5yZXMuaW4=

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.