Distribution of Genes Related to Probiotic Effects Across Lacticaseibacillus rhamnosus Revealed by Population Structure.

Abstract

The Gram-positive Lacticaseibacillus rhamnosus has been broadly reported as capable of exerting beneficial health effects. Bacterial genomic diversity may promote niche specialization, thus creating subpatterns within populations. As L. rhamnosus advantageous effects have been widely reported at strain level and few is known regarding the distribution of beneficial genes among L. rhamnosus strains, we investigated all publicly available genomes of Lactobacillus and Lacticaseibacillus genera to study the pangenome and general population structure of L. rhamnosus. Core genome multilocus sequence typing detected eight L. rhamnosus phylogroups (PG1 to PG8). L. rhamnosus harbors an open pangenome; PG1, PG3, PG4, and PG5 exhibited highly conserved gene distribution patterns. Genes significantly associated to the PG1, which comprises L. rhamnosus GG, are mainly phage-related. The adhesion operon spaCBA-srtC1 was found in 44 (24.7%) genomes; however, considering only the PG1, the prevalence was of 65%. In PG2 the spaCBA-srtC1 prevalence was of 43%. Nevertheless, both human and milk-derived strains harbored this operon. Further, two main types of bacteriocin clusters were found (Bact1 and Bact2). Bact1 predictions indicate the presence of garQ, encoding the class II bacteriocin garvieacin Q, that is mainly present in the closely related PG8A and a PG2 subcluster. PG2 harbors two distinct subclusters, harboring either spaCBA-srtC1 or Bact1. Our findings provide novel insights on the distribution of biotechnological relevant genes across L. rhamnosus population, uncovering intra-species patterns that may bring forth the development of more efficient probiotic products.