Phenotypic characterization of Gardnerella vaginalis subgroups suggests differences in their virulence potential

The well-known genotypic and phenotypic diversity of G. vaginalis resulted in its classification into at least four subgroups (clades) with diverse genomic properties. To evaluate the virulence potential of G. vaginalis subgroups, we analyzed the virulence-related phenotypic characteristics of 14 isolates of Glade 1, 12 isolates of Glade 2, 8 isolates of Glade 4 assessing their in vitro ability to grow as a biofilm, produce the toxin vaginolysin, and express sialidase activity. Significant differences in VLY production were found (p = 0.023), but further analysis of Glade pairs did not confirm this finding. The amount of biofim did not differ significantly among the clades. Analysis of sialidase activity indicated statistically significant differences among the clades (p < 0.001). Production of active recombinant G. vaginalis sialidase demonstrated the link between the sld gene and enzymatic activity, which may be differentially regulated at the transcriptional level. Statistical classification analysis (random forests algorithm) showed that G. vaginalis clades could be best defined by the profiles of two phenotypic characteristics: sialidase activity and vaginolysin production. The results of principal component analysis and hierarchical clustering suggested that all isolates can be subgrouped into three clusters, the structures of which are determined based on phenotypic characteristics of the isolates. Clade 4 was the most homogenous group, as all isolates were found in the same cluster, which is characterized by low production of all studied virulence factors. Clade 2 isolates were mainly distributed between two clusters, whereas Glade 1 isolates were found in all three clusters that were characterized by a distinct profile of phenotypic characteristics. Our findings suggest that G. vaginalis subgroups with different virulence potential might play distinct roles in vaginal microbiota.