A Generalist Lifestyle Allows Rare Gardnerella spp. to Persist at Low Levels in the Vaginal Microbiome

Gardnerella diversity and ecology in pregnancy and preterm birth

The vaginal microbiome has been linked to negative health outcomes including preterm birth. Specific taxa, including Gardnerella spp., have been identified as risk factors for these conditions. Historically, microbiome analysis methods have treated all Gardnerella spp. as one species, but the broad diversity of Gardnerella has become more apparent. We explore the diversity of Gardnerella clades and genomic species in the vaginal microbiome of pregnant women and their associations with microbiome composition and preterm birth. Relative abundance of Gardnerella clades and genomic species and other taxa was quantified in shotgun metagenomic sequencing data from three distinct cohorts of pregnant women. We also assessed the diversity and abundance of Gardnerella variants in 16S rRNA gene amplicon sequencing data from seven previously conducted studies in differing populations. Individual microbiomes often contained multiple Gardnerella variants, and the number of clades was associated with increased microbial load, or the ratio of non-human reads to human reads. Taxon co-occurrence patterns were largely consistent across Gardnerella clades and among cohorts. Some variants previously described as rare were prevalent in other cohorts, highlighting the importance of surveying a diverse set of populations to fully capture the diversity of Gardnerella. The diversity of Gardnerella both across populations and within individual vaginal microbiomes has long been unappreciated, as has been the intra-species diversity of many other members of the vaginal microbiome. The broad genomic diversity of Gardnerella has led to its reclassification as multiple species; here we demonstrate the diversity of Gardnerella found within and between vaginal microbiomes.IMPORTANCEThe present study shows that single microbiomes can contain all currently known species of Gardnerella and that multiple similar species can exist within the same environment. Furthermore, surveys of demographically distinct populations suggest that some species appear more commonly in certain populations. Further studies in broad and diverse populations will be necessary to fully understand the ecological roles of each Gardnerella sp., how they can co-exist, and their distinct impacts on microbial communities, preterm birth, and other health outcomes.

Vaginal microbiomes show ethnic evolutionary dynamics and positive selection of Lactobacillus adhesins driven by a long-term niche-specific process

The vaginal microbiome's composition varies among ethnicities. However, the evolutionary landscape of the vaginal microbiome in the multi-ethnic context remains understudied. We perform a systematic evolutionary analysis of 351 vaginal microbiome samples from 35 multi-ethnic pregnant women, in addition to two validation cohorts, totaling 462 samples from 90 women. Microbiome alpha diversity and community state dynamics show strong ethnic signatures. Lactobacillaceae have a higher ratio of non-synonymous to synonymous polymorphism and lower nucleotide diversity than non-Lactobacillaceae in all ethnicities, with a large repertoire of positively selected genes, including the mucin-binding and cell wall anchor genes. These evolutionary dynamics are driven by the long-term evolutionary process unique to the human vaginal niche. Finally, we propose an evolutionary model reflecting the environmental niches of microbes. Our study reveals the extensive ethnic signatures in vaginal microbial ecology and evolution, highlighting the importance of studying the host-microbiome ecosystem from an evolutionary perspective.

Transport and Utilization of Glycogen Breakdown Products by Gardnerella spp. from the Human Vaginal Microbiome

Multiple Gardnerella species frequently cooccur in vaginal microbiomes, and several factors, including competition for nutrients such as glycogen could determine their population structure. Although Gardnerella spp. can hydrolyze glycogen to produce glucose, maltose, maltotriose, and maltotetraose, how these sugars are transported and utilized for growth is unknown. We determined the distribution of genes encoding transporter proteins associated with the uptake of glucose, maltose, and malto-oligosaccharides and maltodextrins among Gardnerella species. A total of five different ABC transporters were identified in Gardnerella spp. of which MusEFGK2I and MalXFGK were conserved across all 15 Gardnerella isolates. RafEFGK and TMSP (trehalose, maltose, sucrose, and palatinose) operons were specific to G. vaginalis while the MalEFG transporter was identified in G. leopoldii only. Although no glucose specific sugar-symporters were identified, putative "glucose/galactose porters" and components of a phosphotransferase system were identified. In laboratory experiments, all Gardnerella isolates grew more in the presence of glucose, maltose, maltotriose, and maltotetraose compared to unsupplemented media. In addition, most isolates (10/15) showed significantly more growth on maltotetraose compared to glucose (Kruskal Wallis, P < 0.05) suggesting their preference for longer chain malto-oligosaccharides. Our findings show that although putative MusEFGK2I and MalXFGK transporters are found in all Gardnerella spp., some species-specific transporters are also present. Observed distribution of genes encoding transporter systems was consistent with laboratory observations that Gardnerella spp. grow better on longer chain malto-oligosaccharides. IMPORTANCE Increased abundance of Gardnerella spp. is a diagnostic characteristic of bacterial vaginosis, an imbalance in the human vaginal microbiome associated with troubling symptoms and negative reproductive health outcomes, including increased transmission of sexually transmitted infections and preterm birth. Competition for nutrients is likely an important factor in causing dramatic shifts in the vaginal microbial community. Gardnerella produces enzymes to digest glycogen, an important nutrient source for vaginal bacteria, but little is known about the mechanisms in Gardnerella for uptake of the products of this digestion, or whether Gardnerella use some or all of the products. Our results indicate that Gardnerella may have evolved to preferentially use a subset of the glycogen breakdown products, which would help them reduce direct competition with some other bacteria in the vagina.

Metabolic Network Models of the Gardnerella Pangenome Identify Key Interactions with the Vaginal Environment

Gardnerella is the primary pathogenic bacterial genus present in the polymicrobial condition known as bacterial vaginosis (BV). Despite BV's high prevalence and associated chronic and acute women's health impacts, the Gardnerella pangenome is largely uncharacterized at both the genetic and functional metabolic levels. Here, we used genome-scale metabolic models to characterize in silico the Gardnerella pangenome metabolic content. We also assessed the metabolic functional capacity in a BV-positive cervicovaginal fluid context. The metabolic capacity varied widely across the pangenome, with 38.15% of all reactions being core to the genus, compared to 49.60% of reactions identified as being unique to a smaller subset of species. We identified 57 essential genes across the pangenome via in silico gene essentiality screens within two simulated vaginal metabolic environments. Four genes, gpsA, fas, suhB, and psd, were identified as core essential genes critical for the metabolic function of all analyzed bacterial species of the Gardnerella genus. Further understanding these core essential metabolic functions could inform novel therapeutic strategies to treat BV. Machine learning applied to simulated metabolic network flux distributions showed limited clustering based on the sample isolation source, which further supports the presence of extensive core metabolic functionality across this genus. These data represent the first metabolic modeling of the Gardnerella pangenome and illustrate strain-specific interactions with the vaginal metabolic environment across the pangenome. IMPORTANCE Bacterial vaginosis (BV) is the most common vaginal infection among reproductive-age women. Despite its prevalence and associated chronic and acute women's health impacts, the diverse bacteria involved in BV infection remain poorly characterized. Gardnerella is the genus of bacteria most commonly and most abundantly represented during BV. In this paper, we use metabolic models, which are a computational representation of the possible functional metabolism of an organism, to investigate metabolic conservation, gene essentiality, and pathway utilization across 110 Gardnerella strains. These models allow us to investigate in silico how strains may differ with respect to their metabolic interactions with the vaginal-host environment.

Review of the evidence of the effects of human papillomavirus infection and Gardnerella vaginalis, and their co-infection on infertility

A prevalent sexually transmitted infection, the human papillomavirus (HPV) is typically obtained just after the first sexual activity. The majority of HPV infections are asymptomatic and temporary. Cervical, anal, penile, vaginal, vulvar, and oropharyngeal cancers can occur due to recurrent infections with high-risk (hr)-HPV strains, generally decades later. Infections with HPV are significantly associated with reproductive function abnormalities. Per recent research, HPV infections may result in male infertility by reducing sperm motility. The hr-HPV infection was a risk factor for miscarriage, and the indiscriminate HPV genotype increased the probability of premature labor unexpectedly. Women's endometrial trophoblastic cell implantation is decreased by HPV. Gardnerella vaginalis (GV), an anaerobic bacterium that is a component of the natural vaginal flora, can be associated with bacterial vaginosis (BV) when it starts to overgrow and emerge as the dominant species. Reduced Lactobacillus species abundance and GV are linked to female infertility. Data from in vitro studies suggests that sialidase produced by GV may facilitate the entry and growth of papilloma and other sexually transmitted viruses. Also, based on some studies conducted in the past, it can be said that GV and BV are associated with the development of uterine cancer. However, there is still not enough information about the exact mechanism of GV and HPV in causing infertility, which requires more research.

Gardnerella subgroup dominant microbiomes are associated with divergent cervicovaginal immune responses in a longitudinal cohort of Kenyan women

Most cervicovaginal microbiome-immunology studies to date have relied on 16S rDNA microbial profiling which does not resolve the molecular subgroups of Gardnerella, believed to be central to the pathogenesis of bacterial vaginosis (BV) and subsequent risk of HIV acquisition. Here we used the cpn60 universal target which in addition to other microbial taxa, resolves four Gardnerella subgroups, for cervicovaginal microbial profiling in a longitudinal cohort of Kenyan women to examine associations with cellular and soluble markers of inflammation and HIV susceptibility. Participants (N = 41) were sampled, contributing 362 samples for microbiome analysis. All non-Lactobacillus dominant microbial communities were associated with high pro-inflammatory cytokine levels. Divergent associations were observed among different Gardnerella subgroup dominated communities with respect to the chemokine IP-10. Specifically, Gardnerella subgroup A dominant and polymicrobial communities were associated with reduced concentrations of IP-10 in adjusted linear mixed models (p<0.0001), compared to microbial communities dominated by Lactobacillus (non-iners) species. However, these associations did not translate to significant differences in the proportion or absolute number of CCR5, HLA-DR and CD38 expressed on cervical CD4⁺ T- cells. These findings suggest that some associations between Gardnerella subgroup dominant microbiomes and mucosal immunity differ and are relevant for the study of BV-pathogenesis and understanding the mechanisms of BV-associated HIV risk.

Vaginal Lactobacillus iners abundance is associated with outcome in antibiotic treatment of bacterial vaginosis and capable of inhibiting Gardnerella

Bacterial vaginosis is characterized as a polymicrobial dysbiosis with the loss of Lactobacillus spp. and growth of multiple anerobic bacteria, including Gardnerella, Prevotella and Atopobium ranked as the top three most abundant. A total of nine Gardnerella genomospecies have been identified, yet the association between their distribution or any exact Lactobacillus species with BV occurrence or prognosis remains controversial. A total of 308 patients and 62 healthy women who sought annual examinations were recruited, with 130 BV patients and 41 healthy women who met our inclusion criteria finally included. Vaginal samples were used for microscopic examination, 16S rRNA sequencing, bacterial culture and isolation. Isolates of Gardnerella vaginalis, Fannyhessae vaginae (used to be called Atopobium vaginae) and Lactobacillus iners were used for competition tests. We found that the relative abundances of Gardnerella, Prevotella and Atopobium were elevated in BV patients compared to healthy people (p<0.0001), yet no significant differences were found among patients with different clinical outcomes (p>0.05). Seven out of nine Gardnerella genomospecies were present in both BV patients and healthy women, and the relative abundances of all detected genomospecies were higher in BV patients (p<0.05). Cured patients possessed higher GS03 than intermediate and failed patients (p=0.005, 0.0337). L. iners was significantly higher in cured patients than in the other two groups (p=0.0021, p<0.0001), and its ability to inhibit the growth of G. vaginalis and F. vaginae was validated. In summary, seven Gardnerella genomospecies were detected in Chinese BV patients, but no association of its distribution and BV occurrence or prognosis was found. The relative abundance of L. iners was higher in cured patients, and its antimicrobial activity against G. vaginalis and F. vaginae was validated through in vitro inhibition experiment. L. iners could become a predictive indicator of clinical outcomes of BV patients, and its antimicrobial function might be beneficial to BV patients.

Research of cervical microbiota alterations with human papillomavirus infection status and women age in Sanmenxia area of China

Background

Human papillomavirus (HPV) infection is the leading cause of cervical cancer. More and more studies discovered that cervical microbiota (CM) composition correlated with HPV infection and the development of cervical cancer. However, more studies need to be implemented to clarify the complex interaction between microbiota and the mechanism of disease development, especially in a specific area of China.

Materials and methods

In this study, 16S rDNA sequencing was applied on 276 Thin-prep Cytologic Test (TCT) samples of patients from the Sanmenxia area. Systematical analysis of the microbiota structure, diversity, group, and functional differences between different HPV infection groups and age groups, and co-occurrence relationships of the microbiota was carried out.

Results

The major microbiota compositions of all patients include Lactobacillus iners, Escherichia coli, Enterococcus faecalis, and Atopobium vaginae at species level, and Staphylococcus, Lactobacillus, Gardnerella, Bosea, Streptococcus, and Sneathia in genus level. Microbiota diversity was found significantly different between HPV-positive (Chao1 index: 98.8869, p < 0.01), unique-268 infected (infections with one of the HPV genotype 52, 56, or 58, 107.3885, p < 0.01), multi-268 infected (infections with two or more of HPV genotype 52, 56, and 58, 97.5337, p = 0.1012), other1 (94.9619, p < 0.05) groups and HPV-negative group (83.5299). Women older than 60 years old have higher microbiota diversity (108.8851, p < 0.01, n = 255) than younger women (87.0171, n = 21). The abundance of Gardnerella and Atopobium vaginae was significantly higher in the HPV-positive group than in the HPV-negative group, while Burkholderiaceae and Mycoplasma were more abundant in the unique-268 group compared to the negative group. Gamma-proteobacteria and Pseudomonas were found more abundant in older than 60 patients than younger groups. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups (COG) analysis revealed the effects on metabolism by microbiota that the metabolism of cells, proteins, and genetic information-related pathways significantly differed between HPV-negative and positive groups. In contrast, lipid metabolism, signal transduction, and cell cycle metabolism pathway significantly differed between multi-268 and negative groups.

Conclusion

The HPV infection status and age of women were related to CM’s diversity and function pathways. The complex CM co-occurrent relationships and their mechanism in disease development need to be further investigated.

Population Density Affects the Outcome of Competition in Co-cultures of Gardnerella Species Isolated from the Human Vaginal Microbiome

Negative frequency-dependent selection is one possible mechanism for maintenance of rare species in communities, but the selective advantage of rare species may be checked at lower overall population densities where resources are abundant. Gardnerella spp. belonging to cpn60 subgroup D, are detected at low levels in vaginal microbiomes and are nutritional generalists relative to other more abundant Gardnerella spp., making them good candidates for negative frequency-dependent selection. The vaginal microbiome is a dynamic environment, and the resulting changes in density of the microbiota may explain why subgroup D never gains dominance. To test this, we co-cultured subgroup D isolates with isolates from the more common and abundant subgroup C. Deep amplicon sequencing of rpoB was used to determine proportional abundance of each isolate at 0 h and 72 h in 152 co-cultures and to calculate change in proportion. D isolates had a positive change in proportional abundance in most co-cultures regardless of initial proportion. Initial density affected the change in proportion of subgroup D isolates either positively or negatively depending on the particular isolates combined, suggesting that growth rate, population density and other intrinsic features of the isolates influenced the outcome. Our results demonstrate that population density is an important factor influencing the outcome of competition between Gardnerella spp. isolated from the human vaginal microbiome.

Characterization of an α-Glucosidase Enzyme Conserved in Gardnerella spp. Isolated from the Human Vaginal Microbiome

Gardnerella spp. in the vaginal microbiome are associated with bacterial vaginosis, in which a lactobacillus-dominated community is replaced with mixed bacteria, including Gardnerella species. Co-occurrence of multiple Gardnerella species in the vaginal environment is common, but different species are dominant in different women. Competition for nutrients, including glycogen, could play an important role in determining the microbial community structure. Digestion of glycogen into products that can be taken up and further processed by bacteria requires the combined activities of several enzymes collectively known as amylases, which belong to glycoside hydrolase family 13 (GH13) within the CAZy classification system. GH13 is a large and diverse family of proteins, making prediction of their activities challenging. SACCHARIS annotation of the GH13 family in Gardnerella resulted in identification of protein domains belonging to eight subfamilies. Phylogenetic analysis of predicted amylase sequences from 26 genomes demonstrated that a putative α-glucosidase-encoding sequence, CG400_06090, was conserved in all Gardnerella spp. The predicted α-glucosidase enzyme was expressed, purified, and functionally characterized. The enzyme was active on a variety of maltooligosaccharides with maximum activity at pH 7. K_m, k_cat, and k_cat/K_m values for the substrate 4-nitrophenyl α-d-glucopyranoside were 8.3 μM, 0.96 min^-1, and 0.11 μM^-1 min^-1, respectively. Glucose was released from maltose, maltotriose, maltotetraose, and maltopentaose, but no products were detected when the enzyme was incubated with glycogen. Our findings show that Gardnerella spp. produce an α-glucosidase enzyme that may contribute to the multistep process of glycogen metabolism by releasing glucose from maltooligosaccharides. IMPORTANCE Increased abundance of Gardnerella spp. is a diagnostic characteristic of bacterial vaginosis, an imbalance in the human vaginal microbiome associated with troubling symptoms, and negative reproductive health outcomes, including increased transmission of sexually transmitted infections and preterm birth. Competition for nutrients is likely an important factor in causing dramatic shifts in the vaginal microbial community but little is known about the contribution of bacterial enzymes to the metabolism of glycogen, a major carbon source available to vaginal bacteria. The significance of our research is characterizing the activity of an enzyme conserved in Gardnerella species that likely contributes to the ability of these bacteria to utilize glycogen.

Slipped-Strand Mispairing in the Gene Encoding Sialidase NanH3 in Gardnerella spp

Cell wall proteins with sialidase activity are involved in carbohydrate assimilation, adhesion to mucosal surfaces, and biofilm formation. Gardnerella spp. inhabit the human vaginal microbiome and encode up to three sialidase enzymes, two of which are suspected to be cell wall associated. Here, we demonstrate that the gene encoding extracellular sialidase NanH3 is found almost exclusively in Gardnerella piotii and the closely related species Gardnerella genome sp. 3, and its presence correlates with a sialidase-positive phenotype in a collection of 112 Gardnerella isolates. The nanH3 gene sequence includes a homopolymeric repeat of cytosines that varies in length within cell populations, indicating that this gene is subject to slipped-strand mispairing, a mechanism of phase variation in bacteria. Variation in the length of the homopolymer sequence results in production of either the full-length sialidase protein or truncated peptides lacking the sialidase domain due to introduction of reading-frame shifts and premature stop codons. Phase variation in NanH3 may be involved in immune evasion or modulation of adhesion to host epithelial cells and formation of biofilms characteristic of the vaginal dysbiosis known as bacterial vaginosis.