Genome-wide association analysis reveals specialization to hosts and niches in multiple species of the <i>Lactobacillaceae</i>

The Lactobacillaceae inhabit diverse environments, but the extent of their habitat adaptation remains unclear and the colonization factors unknown. First, we applied multiple machine learning models to determine if we can distinguish strains of the same species isolated from two different habitats based on their gene content. Surprisingly, we show that no species is differentially adapted to the oral cavity versus the human gut, or food versus the human gut, while only Lactobacillus crispatus showed specialization to the human urogenital system versus human gut. We then asked which species of Lactobacillaceae are habitat-specialized and how they could be identified. Using multiple lifestyle predictors incorporated in logistic regression models, we found that Limosilactobacillus reuteri, Ligilactobacillus ruminis, Ligilactobacillus salivarius, L. crispatus, and Limosilactobacillus mucosae displayed the highest degrees of host specialization. Applying our microbial genome-wide association study tool, aurora, to these species identified genes encoding adhesins and bacteriocins as the strongest and most common adaptation factors. This work establishes a generalizable framework for identifying novel species-habitat pairs with strong evidence of specialization and for uncovering the genomic features underlying within-species host and habitat adaptation.