Functional Anabolic Network Analysis of Human-Associated Lactobacillus Strains

Members of the Lactobacillus genus are frequently utilized in the probiotic industry with many species conferring demonstrated health benefits; however, these effects are largely strain-dependent. We designed a method called PROTEAN (Probabilistic Reconstruction Of constituent Anabolic Networks) to computationally analyze the genomic annotations and predicted metabolic production capabilities of 144 strains across 16 species of Lactobacillus isolated from human intestinal, oral, and vaginal body sites. Using PROTEAN we conducted a genome-scale metabolic network comparison between strains, revealing that metabolic capabilities differ by isolation site. Notably, PROTEAN does not require a well-curated genome-scale metabolic network reconstruction to provide biological insights. We found that predicted metabolic capabilities of lactobacilli isolated from the vaginal microbiota cluster separately from intestinal and oral isolates, and we also uncovered an overlap in the predicted metabolic production capabilities of intestinal and oral isolates. Using machine learning, we determined the most informative metabolic products driving the difference between predicted metabolic capabilities of intestinal, oral, and vaginal isolates. Notably, intestinal and oral isolates were predicted to have a higher likelihood of producing D-alanine, D/L-serine, and L-proline, while the vaginal isolates were distinguished by a higher predicted likelihood of producing L-arginine, citrulline, and D/L-lactate. We found the distinguishing products to be consistent with published experimental literature. This study showcases a systematic technique, PROTEAN, for comparing the predicted functional metabolic output of microbes using genome-scale metabolic network analysis and computational modeling and provides unique insight into human-associated Lactobacillus biology.