Analysis of the Lactobacillus casei supragenome and its influence in species evolution and lifestyle adaptation

Phylogenomics and comparative genomics of Lactobacillus salivarius, a mammalian gut commensal

The genus Lactobacillus is a diverse group with a combined species count of over 200. They are the largest group within the lactic acid bacteria and one of the most important bacterial groups involved in food microbiology and human nutrition because of their fermentative and probiotic properties. Lactobacillus salivarius, a species commonly isolated from the gastrointestinal tract of humans and animals, has been described as having potential probiotic properties and results of previous studies have revealed considerable functional diversity existing on both the chromosomes and plasmids. Our study consists of comparative genomic analyses of the functional and phylogenomic diversity of 42 genomes of strains of L. salivarius using bioinformatic techniques. The main aim of the study was to describe intra-species diversity and to determine how this diversity is spread across the replicons. We found that multiple phylogenomic and non-phylogenomic methods used for reconstructing trees all converge on similar tree topologies, showing that different metrics largely agree on the evolutionary history of the species. The greatest genomic variation lies on the small plasmids, followed by the repA-type circular megaplasmid, with the chromosome varying least of all. Additionally, the presence of extra linear and circular megaplasmids is noted in several strains, while small plasmids are not always present. Glycosyl hydrolases, bacteriocins and proteases vary considerably on all replicons while two exopolysaccharide clusters and several clustered regularly interspaced short palindromic repeats-associated systems show a lot of variation on the chromosome. Overall, despite its reputation as a mammalian gastrointestinal tract specialist, the intra-specific variation of L. salivarius reveals potential strain-dependant effects on human health.

The Genomic Basis of Lactobacilli as Health-Promoting Organisms

Lactobacilli occupy a unique position in human culture and scientific history. Like brewer’s and baker’s yeast, lactobacilli have been associated with food production and preservation for thousands of years. Lactobacillus species are used in mixed microbial cultures, such as the classical Lactobacillus bulgaricus / Streptococcus thermophilus inoculum for yogurt fermentation, or combinations of diverse lactobacilli/yeasts in kefir grains. The association of lactobacilli consumption with greater longevity and improved health formed the basis for developing lactobacilli as probiotics, whose market has exploded worldwide in the past 10 years. The decade that followed the determination of the first genome sequence of a food-associated species, Lactobacillus plantarum , saw the application to lactobacilli of a full range of functional genomics methods to identify the genes and gene products that govern their distinctive phenotypes and health associations. In this review, we will briefly remind the reader of the range of beneficial effects attributed to lactobacilli, and then explain the phylogenomic basis for the distribution of these traits across the genus. Recognizing the strain specificity of probiotic effects, we review studies of intraspecific genomic variation and their contributions to identifying probiotic traits. Finally we offer a perspective on classification of lactobacilli into new genera in a scheme that will make attributing probiotic properties to clades, taxa, and species more logical and more robust.

Genetic, enzymatic and metabolite profiling of the Lactobacillus casei group reveals strain biodiversity and potential applications for flavour diversification

AIMS

The Lactobacillus casei group represents a widely explored group of lactic acid bacteria, characterized by a high level of biodiversity. In this study, the genetic and phenotypic diversity of a collection of more than 300 isolates of the Lact. casei group and their potential to produce volatile metabolites important for flavour development in dairy products, was examined.

METHODS AND RESULTS

Following confirmation of species by 16S rRNA PCR, the diversity of the isolates was determined by pulsed-field gel electrophoresis. The activities of enzymes involved in the proteolytic cascade were assessed and significant differences among the strains were observed. Ten strains were chosen based on the results of their enzymes activities and they were analysed for their ability to produce volatiles in media with increased concentrations of a representative aromatic, branched chain and sulphur amino acid. Volatiles were assessed using gas chromatography coupled with mass spectrometry. Strain-dependent differences in the range and type of volatiles produced were evident.

CONCLUSIONS

Strains of the Lact. casei group are characterized by genetic and metabolic diversity which supports variability in volatile production.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study provides a screening approach for the knowledge-based selection of strains potentially enabling flavour diversification in fermented dairy products.

CRISPR-Cas Technologies and Applications in Food Bacteria

Clustered regularly interspaced short palindromic repeats (CRISPRs) and CRISPR-associated (Cas) proteins form adaptive immune systems that occur in many bacteria and most archaea. In addition to protecting bacteria from phages and other invasive mobile genetic elements, CRISPR-Cas molecular machines can be repurposed as tool kits for applications relevant to the food industry. A primary concern of the food industry has long been the proper management of food-related bacteria, with a focus on both enhancing the outcomes of beneficial microorganisms such as starter cultures and probiotics and limiting the presence of detrimental organisms such as pathogens and spoilage microorganisms. This review introduces CRISPR-Cas as a novel set of technologies to manage food bacteria and offers insights into CRISPR-Cas biology. It primarily focuses on the applications of CRISPR-Cas systems and tools in starter cultures and probiotics, encompassing strain-typing, phage resistance, plasmid vaccination, genome editing, and antimicrobial activity.

Genotyping by PCR and High-Throughput Sequencing of Commercial Probiotic Products Reveals Composition Biases

Recent advances in microbiome research have brought renewed focus on beneficial bacteria, many of which are available in food and dietary supplements. Although probiotics have historically been defined as microorganisms that convey health benefits when ingested in sufficient viable amounts, this description now includes the stipulation “well defined strains,” encompassing definitive taxonomy for consumer consideration and regulatory oversight. Here, we evaluated 52 commercial dietary supplements covering a range of labeled species using plate counting and targeted genotyping. Strain identities were assessed using methods recently published by the United States Pharmacopeial Convention. We also determined the relative abundance of individual bacteria by high-throughput sequencing (HTS) of the 16S rRNA sequence using paired-end 2 × 250 bp Illumina MiSeq technology. Using these methods, we tested the hypothesis that products do contain the quantitative and qualitative list of labeled microbial species. We found that 17 samples (33%) were below label claim for CFU prior to their expiration dates. A multiplexed-PCR scheme showed that only 30/52 (58%) of the products contained a correctly labeled classification, with issues encompassing incorrect taxonomy, missing species, and un-labeled species. The HTS revealed that many blended products consisted predominantly of Lactobacillus acidophilus and Bifidobacterium animalis subsp. lactis. These results highlight the need for reliable methods to determine the correct taxonomy and quantify the relative amounts of mixed microbial populations in commercial probiotic products.

Functional analysis of a novel hydrogen peroxide resistance gene in Lactobacillus casei strain Shirota

Lactic acid bacteria have a variety of mechanisms for tolerance to oxygen and reactive oxygen species, and these mechanisms differ among species. Lactobacillus casei strain Shirota grows well under aerobic conditions, indicating that the various systems involved in oxidative stress resistance function in this strain. To elucidate the mechanism of oxidative stress resistance in L. casei strain Shirota, we examined the transcriptome response to oxygen or hydrogen peroxide exposure. We then focused on an uncharacterized gene that was found to be up-regulated by both oxygen and hydrogen peroxide stress; we named the gene hprA1 (hydrogen peroxide resistance gene). This gene is widely distributed among lactobacilli. We investigated the involvement of this gene in oxidative stress resistance, as well as the mechanism of tolerance to hydrogen peroxide. Growth of L. casei MS105, an hprA1-disrupted mutant, was not affected by oxygen stress, whereas the survival rate of MS105 after hydrogen peroxide treatment was markedly reduced compared to that of the wild-type. However, the activity of MS105 in eliminating hydrogen peroxide was similar to that of the wild-type. We cloned hprA1 from L. caseiShirota and purified recombinant HprA1 protein from Escherichia coli. We demonstrated that the recombinant HprA1 protein bound to iron and prevented the formation of a hydroxyl radical in vitro. Thus, HprA1 protein probably contributes to hydrogen peroxide tolerance in L. casei strain Shirota by binding to iron in the cells and preventing the formation of a hydroxyl radical.

CRISPR technologies for bacterial systems: Current achievements and future directions

Throughout the decades of its history, the advances in bacteria-based bio-industries have coincided with great leaps in strain engineering technologies. Recently unveiled clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated proteins (Cas) systems are now revolutionizing biotechnology as well as biology. Diverse technologies have been derived from CRISPR/Cas systems in bacteria, yet the applications unfortunately have not been actively employed in bacteria as extensively as in eukaryotic organisms. A recent trend of engineering less explored strains in industrial microbiology-metabolic engineering, synthetic biology, and other related disciplines-is demanding facile yet robust tools, and various CRISPR technologies have potential to cater to the demands. Here, we briefly review the science in CRISPR/Cas systems and the milestone inventions that enabled numerous CRISPR technologies. Next, we describe CRISPR/Cas-derived technologies for bacterial strain development, including genome editing and gene expression regulation applications. Then, other CRISPR technologies possessing great potential for industrial applications are described, including typing and tracking of bacterial strains, virome identification, vaccination of bacteria, and advanced antimicrobial approaches. For each application, we note our suggestions for additional improvements as well. In the same context, replication of CRISPR/Cas-based chromosome imaging technologies developed originally in eukaryotic systems is introduced with its potential impact on studying bacterial chromosomal dynamics. Also, the current patent status of CRISPR technologies is reviewed. Finally, we provide some insights to the future of CRISPR technologies for bacterial systems by proposing complementary techniques to be developed for the use of CRISPR technologies in even wider range of applications.

Stress Physiology of Lactic Acid Bacteria

Lactic acid bacteria (LAB) are important starter, commensal, or pathogenic microorganisms. The stress physiology of LAB has been studied in depth for over 2 decades, fueled mostly by the technological implications of LAB robustness in the food industry. Survival of probiotic LAB in the host and the potential relatedness of LAB virulence to their stress resilience have intensified interest in the field. Thus, a wealth of information concerning stress responses exists today for strains as diverse as starter (e.g., Lactococcus lactis), probiotic (e.g., several Lactobacillus spp.), and pathogenic (e.g., Enterococcus and Streptococcus spp.) LAB. Here we present the state of the art for LAB stress behavior. We describe the multitude of stresses that LAB are confronted with, and we present the experimental context used to study the stress responses of LAB, focusing on adaptation, habituation, and cross-protection as well as on self-induced multistress resistance in stationary phase, biofilms, and dormancy. We also consider stress responses at the population and single-cell levels. Subsequently, we concentrate on the stress defense mechanisms that have been reported to date, grouping them according to their direct participation in preserving cell energy, defending macromolecules, and protecting the cell envelope. Stress-induced responses of probiotic LAB and commensal/pathogenic LAB are highlighted separately due to the complexity of the peculiar multistress conditions to which these bacteria are subjected in their hosts. Induction of prophages under environmental stresses is then discussed. Finally, we present systems-based strategies to characterize the "stressome" of LAB and to engineer new food-related and probiotic LAB with improved stress tolerance.

The Variable Regions of Lactobacillus rhamnosus Genomes Reveal the Dynamic Evolution of Metabolic and Host-Adaptation Repertoires

Lactobacillus rhamnosus is a diverse Gram-positive species with strains isolated from different ecological niches. Here, we report the genome sequence analysis of 40 diverse strains of L. rhamnosus and their genomic comparison, with a focus on the variable genome. Genomic comparison of 40 L. rhamnosus strains discriminated the conserved genes (core genome) and regions of plasticity involving frequent rearrangements and horizontal transfer (variome). The L. rhamnosus core genome encompasses 2,164 genes, out of 4,711 genes in total (the pan-genome). The accessory genome is dominated by genes encoding carbohydrate transport and metabolism, extracellular polysaccharides (EPS) biosynthesis, bacteriocin production, pili production, the cas system, and the associated clustered regularly interspaced short palindromic repeat (CRISPR) loci, and more than 100 transporter functions and mobile genetic elements like phages, plasmid genes, and transposons. A clade distribution based on amino acid differences between core (shared) proteins matched with the clade distribution obtained from the presence–absence of variable genes. The phylogenetic and variome tree overlap indicated that frequent events of gene acquisition and loss dominated the evolutionary segregation of the strains within this species, which is paralleled by evolutionary diversification of core gene functions. The CRISPR-Cas system could have contributed to this evolutionary segregation. Lactobacillus rhamnosus strains contain the genetic and metabolic machinery with strain-specific gene functions required to adapt to a large range of environments. A remarkable congruency of the evolutionary relatedness of the strains’ core and variome functions, possibly favoring interspecies genetic exchanges, underlines the importance of gene-acquisition and loss within the L. rhamnosus strain diversification.

Lactobacillus casei as a biocatalyst for biofuel production

Microbial fermentation of sugars from plant biomass to alcohols represents an alternative to petroleum-based fuels. The optimal biocatalyst for such fermentations needs to overcome hurdles such as high concentrations of alcohols and toxic compounds. Lactic acid bacteria, especially lactobacilli, have high innate alcohol tolerance and are remarkably adaptive to harsh environments. This study assessed the potential of five Lactobacillus casei strains as biocatalysts for alcohol production. L. casei 12A was selected based upon its innate alcohol tolerance, high transformation efficiency and ability to utilize plant-derived carbohydrates. A 12A derivative engineered to produce ethanol (L. casei E1) was compared to two other bacterial biocatalysts. Maximal growth rate, maximal optical density and ethanol production were determined under conditions similar to those present during alcohol production from lignocellulosic feedstocks. L. casei E1 exhibited higher innate alcohol tolerance, better growth in the presence of corn stover hydrolysate stressors, and resulted in higher ethanol yields.

The Impact of Lactobacillus casei on the Composition of the Cecal Microbiota and Innate Immune System Is Strain Specific

The probiotic function to impact human health is thought to be related to their ability to alter the composition of the gut microbiota and modulate the human innate immune system. The ability to function as a probiotic is believed to be strain specific. Strains of Lactobacillus casei are commonly utilized as probiotics that when consumed alter the composition of the gut microbiota and modulate the host immune response. L. casei strains are known to differ significantly in gene content. The objective of this study was to investigate seven different L. casei strains for their ability to alter the murine gut microbiota and modulate the murine immune system. C57BL/6 mice were fed L. casei strains at a dose of 10⁸ CFU/day/mouse for seven days and sacrificed 3.5h after the last administration. The cecal content and the ileum tissue were collected for microbiota analysis and immune profiling, respectively. While 5 of the L. casei strains altered the gut microbiota in a strain specific manner, two of the strains did not alter the overall cecal microbiota composition. The observed changes cluster into three groups containing between 1 and 2 strains. Two strains that did not affect the gut microbiota composition cluster together with the control in their impact on pattern recognition receptors (PRRs) expression, suggesting that the ability to alter the cecal microbiota correlates with the ability to alter PRR expression. They also cluster together in their impact on the expression of intestinal antimicrobial peptides (AMPs). This result suggests that a relationship exists between the capability of a L. casei strains to alter the composition of the gut microbiota, PRR regulation, and AMP regulation.

Draft Genome Sequence of the Respiration-Competent Strain Lactobacillus casei N87

Lactobacillus casei is used as a starter, adjunct, and/or probiotic culture in the production of fermented and functional foods. Here, we report the draft genome sequence of the respiration-competent strain L. casei N87, isolated from infant feces. This genome information may be useful for the study of respiratory metabolism in lactic acid bacteria.

Draft Genome Sequence of Lactobacillus casei DPC6800, an Isolate with the Potential to Diversify Flavor in Cheese

Lactobacillus casei is a nonstarter lactic acid bacterium commonly present in various types of cheeses. It is believed that strains of this species have a significant impact on the development of cheese flavor. The draft genome sequence of L. casei DPC6800, isolated from a semi-hard Dutch cheese, is reported.

CRISPR-Based Typing and Next-Generation Tracking Technologies

Bacteria occur ubiquitously in nature and are broadly relevant throughout the food supply chain, with diverse and variable tolerance levels depending on their origin, biological role, and impact on the quality and safety of the product as well as on the health of the consumer. With increasing knowledge of and accessibility to the microbial composition of our environments, food supply, and host-associated microbiota, our understanding of and appreciation for the ratio of beneficial to undesirable bacteria are rapidly evolving. Therefore, there is a need for tools and technologies that allow definite, accurate, and high-resolution identification and typing of various groups of bacteria that include beneficial microbes such as starter cultures and probiotics, innocuous commensals, and undesirable pathogens and spoilage organisms. During the transition from the current molecular biology-based PFGE (pulsed-field gel electrophoresis) gold standard to the increasingly accessible omics-level whole-genome sequencing (WGS) N-gen standard, high-resolution technologies such as CRISPR-based genotyping constitute practical and powerful alternatives that provide valuable insights into genome microevolution and evolutionary trajectories. Indeed, several studies have shown potential for CRISPR-based typing of industrial starter cultures, health-promoting probiotic strains, animal commensal species, and problematic pathogens. Emerging CRISPR-based typing methods open new avenues for high-resolution typing of a broad range of bacteria and constitute a practical means for rapid tracking of a diversity of food-associated microbes.

The Effect of Lactobacillus casei 32G on the Mouse Cecum Microbiota and Innate Immune Response Is Dose and Time Dependent

Lactobacilli have been associated with a variety of immunomodulatory effects and some of these effects have been related to changes in gastrointestinal microbiota. However, the relationship between probiotic dose, time since probiotic consumption, changes in the microbiota, and immune system requires further investigation. The objective of this study was to determine if the effect of Lactobacillus casei 32G on the murine gastrointestinal microbiota and immune function are dose and time dependent. Mice were fed L. casei 32G at doses of 10⁶, 10⁷, or 10⁸ CFU/day/mouse for seven days and were sacrificed 0.5h, 3.5h, 12h, or 24h after the last administration. The ileum tissue and the cecal content were collected for immune profiling by qPCR and microbiota analysis, respectively. The time required for L. casei 32G to reach the cecum was monitored by qPCR and the 32G bolus reaches the cecum 3.5h after the last administration. L. casei 32G altered the cecal microbiota with the predominance of Lachnospiraceae IS, and Oscillospira decreasing significantly (p < 0.05) in the mice receiving 10⁸ CFU/mouse 32G relative to the control mice, while a significant (p < 0.05) increase was observed in the prevalence of lactobacilli. The lactobacilli that increased were determined to be a commensal lactobacilli. Interestingly, no significant difference in the overall microbiota composition, regardless of 32G doses, was observed at the 12h time point. A likely explanation for this observation is the level of feed derived-nutrients resulting from the 12h light/dark cycle. 32G results in consistent increases in Clec2h expression and reductions in TLR-2, alpha-defensins, and lysozyme. Changes in expression of these components of the innate immune system are one possible explanation for the observed changes in the cecal microbiota. Additionally, 32G administration was observed to alter the expression of cytokines (IL-10rb and TNF-α) in a manner consistent with an anti-inflammatory response.

Human Gut-Commensalic Lactobacillus ruminis ATCC 25644 Displays Sortase-Assembled Surface Piliation: Phenotypic Characterization of Its Fimbrial Operon through In Silico Predictive Analysis and Recombinant Expression in Lactococcus lactis

Sortase-dependent surface pili (or fimbriae) in Gram-positive bacteria are well documented as a key virulence factor for certain harmful opportunistic pathogens. However, it is only recently known that these multi-subunit protein appendages are also belonging to the “friendly” commensals and now, with this new perspective, they have come to be categorized as a niche-adaptation factor as well. In this regard, it was shown earlier that sortase-assembled piliation is a native fixture of two human intestinal commensalics (i.e., Lactobacillus rhamnosus and Bifidobacterium bifidum), and correspondingly where the pili involved have a significant role in cellular adhesion and immunomodulation processes. We now reveal that intestinal indigenous (or autochthonous) Lactobacillus ruminis is another surface-piliated commensal lactobacillar species. Heeding to in silico expectations, the predicted loci for the LrpCBA-called pili are organized tandemly in the L. ruminis genome as a canonical fimbrial operon, which then encodes for three pilin-proteins and a single C-type sortase enzyme. Through electron microscopic means, we showed that these pilus formations are a surface assemblage of tip, basal, and backbone pilin subunits (respectively named LrpC, LrpB, and LrpA) in L. ruminis, and also when expressed recombinantly in Lactococcus lactis. As well, by using the recombinant-piliated lactococci, we could define certain ecologically relevant phenotypic traits, such as the ability to adhere to extracellular matrix proteins and gut epithelial cells, but also to effectuate an induced dampening on Toll-like receptor 2 signaling and interleukin-8 responsiveness in immune-related cells. Within the context of the intestinal microcosm, by wielding such niche-advantageous cell-surface properties the LrpCBA pilus would undoubtedly have a requisite functional role in the colonization dynamics of L. ruminis indigeneity. Our study provides only the second description of a native-piliated Lactobacillus species, but at the same time also involves the structural and functional characterization of a third type of lactobacillar pilus.

Diversity of CRISPR-Cas immune systems and molecular machines

Bacterial adaptive immunity hinges on CRISPR-Cas systems that provide DNA-encoded, RNA-mediated targeting of exogenous nucleic acids. A plethora of CRISPR molecular machines occur broadly in prokaryotic genomes, with a diversity of Cas nucleases that can be repurposed for various applications.

CRISPR-Based Technologies and the Future of Food Science

The on-going CRISPR craze is focused on the use of Cas9-based technologies for genome editing applications in eukaryotes, with high potential for translational medicine and next-generation gene therapy. Nevertheless, CRISPR-Cas systems actually provide adaptive immunity in bacteria, and have much promise for various applications in food bacteria that include high-resolution typing of pathogens, vaccination of starter cultures against phages, and the genesis of programmable and specific antibiotics that can selectively modulate bacterial population composition. Indeed, the molecular machinery from these DNA-encoded, RNA-mediated, DNA-targeting systems can be harnessed in native hosts, or repurposed in engineered systems for a plethora of applications that can be implemented in all organisms relevant to the food chain, including agricultural crops trait-enhancement, livestock breeding, and fermentation-based manufacturing, and for the genesis of next-generation food products with enhanced quality and health-promoting functionalities. CRISPR-based applications are now poised to revolutionize many fields within food science, from farm to fork. In this review, we describe CRISPR-Cas systems and highlight their potential for the development of enhanced foods.

The vaginal isolate Lactobacillus paracasei LPC-S01 (DSM 26760) is suitable for oral administration

Bacterial vaginosis is one of the most common urogenital diseases affecting women in reproductive age. The administration of probiotics as vaginal suppository has been proposed as a strategy to cure this condition and reduce its recurrence. Nonetheless, also oral consumption of probiotics, which is a more practical route of administration, proved to be an efficient strategy. In this perspective, we studied Lactobacillus paracasei LPC-S01 (DSM 26760), a human vaginal isolate included in commercial probiotic preparations for topical use, in order to assess if this bacterium can also perform as gastrointestinal probiotic. Comparative genomics revealed the presence of several accessory genes suggesting that LPC-S01 is a niche-generalist member of its species. According to a procedure conventionally used to predict the probiotic potential, we demonstrated that the probiotic properties of strain LPC-S01, with respect to those of the well-known probiotic references L. paracasei Shirota and DG, are equal for the bile tolerance and the reduction of NF-κB activation in Caco-2 cells, or superior for the tolerance to gastric juice and the adhesion to Caco-2 epithelial cells. We then demonstrated that LPC-S01 is susceptible to antibiotics indicated by EFSA and does not produce biogenic amines. Finally, a double-blind cross-over pilot intervention trial on healthy human volunteers showed that, after a 7-days oral consumption of capsules containing about 24 billion live cells, the fecal cell concentrations of strains LPC-S01 and DG (evaluated by qPCR) were not dissimilar. Specifically, both probiotics' cell concentrations were above the detection limit for an average of 5 days from the end of the treatment, corresponding to a mean number of evacuations of 7 ± 2. Taken together, these data demonstrate that the vaginal isolate L. paracasei LPC-S01 possesses safety and functional properties that may support its use as probiotic to be administered per os for potential intestinal as well as vaginal applications.

The bacterial pangenome as a new tool for analysing pathogenic bacteria

The bacterial pangenome was introduced in 2005 and, in recent years, has been the subject of many studies. Thanks to progress in next-generation sequencing methods, the pangenome can be divided into two parts, the core (common to the studied strains) and the accessory genome, offering a large panel of uses. In this review, we have presented the analysis methods, the pangenome composition and its application as a study of lifestyle. We have also shown that the pangenome may be used as a new tool for redefining the pathogenic species. We applied this to the Escherichia coli and Shigella species, which have been a subject of controversy regarding their taxonomic and pathogenic position.

•

Pangenome is a new way of studying pathogenic bacteria.

•

Pangenome can be used as a taxonomic tool.

•

This review describes pangenome in the world of pathogenic bacteria.

Correlation of Lactobacillus rhamnosus Genotypes and Carbohydrate Utilization Signatures Determined by Phenotype Profiling

Lactobacillus rhamnosus is a bacterial species commonly colonizing the gastrointestinal (GI) tract of humans and also frequently used in food products. While some strains have been studied extensively, physiological variability among isolates of the species found in healthy humans or their diet is largely unexplored. The aim of this study was to characterize the diversity of carbohydrate utilization capabilities of human isolates and food-derived strains of L. rhamnosus in relation to their niche of isolation and genotype. We investigated the genotypic and phenotypic diversity of 25 out of 65 L. rhamnosus strains from various niches, mainly human feces and fermented dairy products. Genetic fingerprinting of the strains by amplified fragment length polymorphism (AFLP) identified 11 distinct subgroups at 70% similarity and suggested niche enrichment within particular genetic clades. High-resolution carbohydrate utilization profiling (OmniLog) identified 14 carbon sources that could be used by all of the strains tested for growth, while the utilization of 58 carbon sources differed significantly between strains, enabling the stratification of L. rhamnosus strains into three metabolic clusters that partially correlate with the genotypic clades but appear uncorrelated with the strain's origin of isolation. Draft genome sequences of 8 strains were generated and employed in a gene-trait matching (GTM) analysis together with the publicly available genomes of L. rhamnosus GG (ATCC 53103) and HN001 for several carbohydrates that were distinct for the different metabolic clusters: l-rhamnose, cellobiose, l-sorbose, and α-methyl-d-glucoside. From the analysis, candidate genes were identified that correlate with l-sorbose and α-methyl-d-glucoside utilization, and the proposed function of these genes could be confirmed by heterologous expression in a strain lacking the genes. This study expands our insight into the phenotypic and genotypic diversity of the species L. rhamnosus and explores the relationships between specific carbohydrate utilization capacities and genotype and/or niche adaptation of this species.

An L-Fucose Operon in the Probiotic Lactobacillus rhamnosus GG Is Involved in Adaptation to Gastrointestinal Conditions

L-Fucose is a sugar present in human secretions as part of human milk oligosaccharides, mucins, and other glycoconjugates in the intestinal epithelium. The genome of the probiotic Lactobacillus rhamnosus GG (LGG) carries a gene cluster encoding a putative L-fucose permease (fucP), L-fucose catabolic pathway (fucI, fucK, fucU, and fucA), and a transcriptional regulator (fucR). The metabolism of L-fucose in LGG results in 1,2-propanediol production, and their fucI and fucP mutants displayed a severe and mild growth defect on L-fucose, respectively. Transcriptional analysis revealed that the fuc genes are induced by L-fucose and subject to a strong carbon catabolite repression effect. This induction was triggered by FucR, which acted as a transcriptional activator necessary for growth on L-fucose. LGG utilized fucosyl-α1,3-N-acetylglucosamine and contrarily to other lactobacilli, the presence of fuc genes allowed this strain to use the L-fucose moiety. In fucI and fucR mutants, but not in fucP mutant, L-fucose was not metabolized and it was excreted to the medium during growth on fucosyl-α1,3-N-acetylglucosamine. The fuc genes were induced by this fucosyl-disaccharide in the wild type and the fucP mutant but not in a fucI mutant, showing that FucP does not participate in the regulation of fuc genes and that L-fucose metabolism is needed for FucR activation. The l-fucose operon characterized here constitutes a new example of the many factors found in LGG that allow this strain to adapt to the gastrointestinal conditions.

Characterizing Diversity of Lactobacilli Associated with Severe Early Childhood Caries: A Study Protocol

Lactobacilli have been consistently associated with dental caries for decades; however, knowledge of this group of bacteria in the etiology of the disease is limited to quantitative elucidation. Nowadays, explicit identification of oral Lactobacillus species is possible, despite their taxonomic complexity. Here we describe a combined approach involving both cultivation and genetic methods to ascertain and characterize the diversity and abundance of the Lactobacillus population in the oral cavities of children with severe early childhood caries (S-ECC). Eighty 3- to 6-year-old children (40 S-ECC and 40 caries free) who were seeking dental care at the Pediatric Dental Clinic of Bellevue Hospital in New York City were invited to participate in this study. Clinical data on socio-demographic information and oral health behavior were obtained from the primary caregiver. The data included a detailed dental examination, children's medical history, and a questionnaire survey. Combined non-stimulated saliva and supra-gingival plaque samples were collected from each child and cultivated on selective media for quantitative measures of lactobacilli levels. The procedure for Lactobacillus species screening will include the random selection of 50 colonies per plate, extraction of DNA from each colony, and genotyping by arbitrarily primed polymerase chain reaction (AP-PCR). Each unique Lactobacillus AP-PCR genotype will be selected for taxonomic assessment by 16S rRNA gene sequencing analysis. Lactobacillus species will be identified by comparing the 16S rRNA sequences with the Ribosomal Database and the Human Oral Microbiome Database. Meanwhile, the same set of clinical samples will be independently subjected to genomic DNA isolation, 16S rRNA amplification with Lactobacillus genus-specific primers, sequencing, and taxonomic identification, both at genus and species levels with a customized pipeline. The distribution and phylogenetic differences of these Lactobacillus species will be compared between children with or without S-ECC. One of the main objectives of this study is to establish a study protocol for the identification and characterization of lactobacilli in the oral cavity. Future caries risk assessments can include lactobacilli counts (quantitative) and the presence/absence of specific cariogenic genetic signatures of a Lactobacillus species (qualitative) associated with S-ECC.

Ten years of pan-genome analyses

Next generation sequencing technologies have engendered a genome sequence data deluge in public databases. Genome analyses have transitioned from single or few genomes to hundreds to thousands of genomes. Pan-genome analyses provide a framework for estimating the genomic diversity of the dataset at hand and predicting the number of additional whole genomes sequences that would be necessary to fully characterize that diversity. We review recent implementations of the pan-genome approach, its impact and limits, and we propose possible extensions, including analyses at the whole genome multiple sequence alignment level.

Comparative genomics of Lactobacillus crispatus suggests novel mechanisms for the competitive exclusion of Gardnerella vaginalis

BACKGROUND

Lactobacillus crispatus is a ubiquitous micro-organism encountered in a wide range of host-associated habitats. It can be recovered from the gastrointestinal tract of animals and it is a common constituent of the vaginal microbiota of humans. Moreover, L. crispatus can contribute to the urogenital health of the host through competitive exclusion and the production of antimicrobial agents. In order to investigate the genetic diversity of this important urogenital species, we performed a comparative genomic analysis of L. crispatus.

RESULTS

Utilizing the completed genome sequence of a strain ST1 and the draft genome sequences of nine other L. crispatus isolates, we defined the scale and scope of the pan- and core genomic potential of L. crispatus. Our comparative analysis identified 1,224 and 2,705 ortholog groups present in all or only some of the ten strains, respectively. Based on mathematical modeling, sequencing of additional L. crispatus isolates would result in the identification of new genes and functions, whereas the conserved core of the ten strains was a good representation of the final L. crispatus core genome, estimated to level at about 1,116 ortholog groups. Importantly, the current core was observed to encode bacterial components potentially promoting urogenital health. Using antibody fragments specific for one of the conserved L. crispatus adhesins, we demonstrated that the L. crispatus core proteins have a potential to reduce the ability of Gardnerella vaginalis to adhere to epithelial cells. These findings thereby suggest that L. crispatus core proteins could protect the vagina from G. vaginalis and bacterial vaginosis.

CONCLUSIONS

Our pan-genome analysis provides insights into the intraspecific genome variability and the collective molecular mechanisms of the species L. crispatus. Using this approach, we described the differences and similarities between the genomes and identified features likely to be important for urogenital health. Notably, the conserved genetic backbone of L. crispatus accounted for close to 60% of the ortholog groups of an average L. crispatus strain and included factors for the competitive exclusion of G. vaginalis, providing an explanation on how this urogenital species could improve vaginal health.

High efficiency electrotransformation of Lactobacillus casei

We investigated whether protocols allowing high efficiency electrotransformation of other lactic acid bacteria were applicable to five strains of Lactobacillus casei (12A, 32G, A2-362, ATCC 334 and BL23). Addition of 1% glycine or 0.9 M NaCl during cell growth, limitation of the growth of the cell cultures to OD600 0.6-0.8, pre-electroporation treatment of cells with water or with a lithium acetate (100 mM)/dithiothreitol (10 mM) solution and optimization of electroporation conditions all improved transformation efficiencies. However, the five strains varied in their responses to these treatments. Transformation efficiencies of 10(6) colony forming units μg(-1) pTRKH2 DNA and higher were obtained with three strains which is sufficient for construction of chromosomal gene knock-outs and gene replacements.

Phenotypical analysis of the Lactobacillus rhamnosus GG fimbrial spaFED operon: surface expression and functional characterization of recombinant SpaFED pili in Lactococcus lactis

A noticeable genomic feature of many piliated Gram-positive bacterial species is the presence of more than one pilus-encoding operon. Paradigmatically, the gut-adapted Lactobacillus rhamnosus GG strain contains two different fimbrial operons in its genome. However, whereas one of these operons (called spaCBA) is encoding for the functionally mucus-/collagen-binding SpaCBA pilus, for the other operon (called spaFED) any native expression of the SpaFED-called pili is still the subject of some uncertainty. Irrespective of such considerations, we decided it would be of relevance or interest to decipher the gross structure of this pilus type, and as well assess its functional capabilities for cellular adhesion and immunostimulation. For this, and by following the approach we had used previously to explicate the immuno-properties of SpaCBA pili, we constructed nisin-inducible expression clones producing either wild-type or SpaF pilin-deleted surface-assembled L. rhamnosus GG SpaFED pili on Lactococcus lactis cells. Using these piliated lactococcal constructs, we found that the pilin-polymerized architecture of a recombinant-produced SpaFED pilus coincides with sequence-based functional predictions of the related pilins, and in fact is prototypical of those other sortase-dependent pilus-like structures thus far characterized for piliated Gram-positive bacteria. Moreover, we confirmed that among the different pilin subunits encompassing spaFED operon-encoded pili, the SpaF pilin is a main adhesion determinant, and when present in the assembled structure can mediate pilus binding to mucus, certain extracellular matrix proteins, and different gut epithelial cell lines. However, somewhat unexpectedly, when recombinant SpaFED pili are surface-attached, we found that they could not potentiate the existing lactococcal cell-induced immune responses so elicited from intestinal- and immune-related cells, but rather instead, they could dampen them. Accordingly, we have now provided the first phenotypical description of a spaFED pilus operon, and with that furthered the functional understanding of surface piliation for a particular gut-commensalic genre of piliated Gram-positive bacteria.

Lactococcus lactis metabolism and gene expression during growth on plant tissues

Lactic acid bacteria have been isolated from living, harvested, and fermented plant materials; however, the adaptations these bacteria possess for growth on plant tissues are largely unknown. In this study, we investigated plant habitat-specific traits of Lactococcus lactis during growth in an Arabidopsis thaliana leaf tissue lysate (ATL). L. lactis KF147, a strain originally isolated from plants, exhibited a higher growth rate and reached 7.9-fold-greater cell densities during growth in ATL than the dairy-associated strain L. lactis IL1403. Transcriptome profiling (RNA-seq) of KF147 identified 853 induced and 264 repressed genes during growth in ATL compared to that in GM17 laboratory culture medium. Genes induced in ATL included those involved in the arginine deiminase pathway and a total of 140 carbohydrate transport and metabolism genes, many of which are involved in xylose, arabinose, cellobiose, and hemicellulose metabolism. The induction of those genes corresponded with L. lactis KF147 nutrient consumption and production of metabolic end products in ATL as measured by gas chromatography-time of flight mass spectrometry (GC-TOF/MS) untargeted metabolomic profiling. To assess the importance of specific plant-inducible genes for L. lactis growth in ATL, xylose metabolism was targeted for gene knockout mutagenesis. Wild-type L. lactis strain KF147 but not an xylA deletion mutant was able to grow using xylose as the sole carbon source. However, both strains grew to similarly high levels in ATL, indicating redundancy in L. lactis carbohydrate metabolism on plant tissues. These findings show that certain strains of L. lactis are well adapted for growth on plants and possess specific traits relevant for plant-based food, fuel, and feed fermentations.

Genome-scale reconstruction of metabolic networks of Lactobacillus casei ATCC 334 and 12A

Lactobacillus casei strains are widely used in industry and the utility of this organism in these industrial applications is strain dependent. Hence, tools capable of predicting strain specific phenotypes would have utility in the selection of strains for specific industrial processes. Genome-scale metabolic models can be utilized to better understand genotype-phenotype relationships and to compare different organisms. To assist in the selection and development of strains with enhanced industrial utility, genome-scale models for L. casei ATCC 334, a well characterized strain, and strain 12A, a corn silage isolate, were constructed. Draft models were generated from RAST genome annotations using the Model SEED database and refined by evaluating ATP generating cycles, mass-and-charge-balances of reactions, and growth phenotypes. After the validation process was finished, we compared the metabolic networks of these two strains to identify metabolic, genetic and ortholog differences that may lead to different phenotypic behaviors. We conclude that the metabolic capabilities of the two networks are highly similar. The L. casei ATCC 334 model accounts for 1,040 reactions, 959 metabolites and 548 genes, while the L. casei 12A model accounts for 1,076 reactions, 979 metabolites and 640 genes. The developed L. casei ATCC 334 and 12A metabolic models will enable better understanding of the physiology of these organisms and be valuable tools in the development and selection of strains with enhanced utility in a variety of industrial applications.

Estimation of prokaryotic supergenome size and composition from gene frequency distributions

BACKGROUND

Because prokaryotic genomes experience a rapid flux of genes, selection may act at a higher level than an individual genome. We explore a quantitative model of the distributed genome whereby groups of genomes evolve by acquiring genes from a fixed reservoir which we denote as supergenome. Previous attempts to understand the nature of the supergenome treated genomes as random, independent collections of genes and assumed that the supergenome consists of a small number of homogeneous sub-reservoirs. Here we explore the consequences of relaxing both assumptions.

RESULTS

We surveyed several methods for estimating the size and composition of the supergenome. The methods assumed that genomes were either random, independent samples of the supergenome or that they evolved from a common ancestor along a known tree via stochastic sampling from the reservoir. The reservoir was assumed to be either a collection of homogeneous sub-reservoirs or alternatively composed of genes with Gamma distributed gain probabilities. Empirical gene frequencies were used to either compute the likelihood of the data directly or first to reconstruct the history of gene gains and then compute the likelihood of the reconstructed numbers of gains.

CONCLUSIONS

Supergenome size estimates using the empirical gene frequencies directly are not robust with respect to the choice of the model. By contrast, using the gene frequencies and the phylogenetic tree to reconstruct multiple gene gains produces reliable estimates of the supergenome size and indicates that a homogeneous supergenome is more consistent with the data than a supergenome with Gamma distributed gain probabilities.

Functional genomics of lactic acid bacteria: from food to health

Genome analysis using next generation sequencing technologies has revolutionized the characterization of lactic acid bacteria and complete genomes of all major groups are now available. Comparative genomics has provided new insights into the natural and laboratory evolution of lactic acid bacteria and their environmental interactions. Moreover, functional genomics approaches have been used to understand the response of lactic acid bacteria to their environment. The results have been instrumental in understanding the adaptation of lactic acid bacteria in artisanal and industrial food fermentations as well as their interactions with the human host. Collectively, this has led to a detailed analysis of genes involved in colonization, persistence, interaction and signaling towards to the human host and its health. Finally, massive parallel genome re-sequencing has provided new opportunities in applied genomics, specifically in the characterization of novel non-GMO strains that have potential to be used in the food industry. Here, we provide an overview of the state of the art of these functional genomics approaches and their impact in understanding, applying and designing lactic acid bacteria for food and health.

Motif depletion in bacteriophages infecting hosts with CRISPR systems

BACKGROUND

CRISPR is a microbial immune system likely to be involved in host-parasite coevolution. It functions using target sequences encoded by the bacterial genome, which interfere with invading nucleic acids using a homology-dependent system. The system also requires protospacer associated motifs (PAMs), short motifs close to the target sequence that are required for interference in CRISPR types I and II. Here, we investigate whether PAMs are depleted in phage genomes due to selection pressure to escape recognition.

RESULTS

To this end, we analyzed two data sets. Phages infecting all bacterial hosts were analyzed first, followed by a detailed analysis of phages infecting the genus Streptococcus, where PAMs are best understood. We use two different measures of motif underrepresentation that control for codon bias and the frequency of submotifs. We compare phages infecting species with a particular CRISPR type to those infecting species without that type. Since only known PAMs were investigated, the analysis is restricted to CRISPR types I-C and I-E and in Streptococcus to types I-C and II. We found evidence for PAM depletion in Streptococcus phages infecting hosts with CRISPR type I-C, in Vibrio phages infecting hosts with CRISPR type I-E and in Streptococcus thermopilus phages infecting hosts with type II-A, known as CRISPR3.

CONCLUSIONS

The observed motif depletion in phages with hosts having CRISPR can be attributed to selection rather than to mutational bias, as mutational bias should affect the phages of all hosts. This observation implies that the CRISPR system has been efficient in the groups discussed here.

A comparative pan-genome perspective of niche-adaptable cell-surface protein phenotypes in Lactobacillus rhamnosus

Lactobacillus rhamnosus is a ubiquitously adaptable Gram-positive bacterium and as a typical commensal can be recovered from various microbe-accessible bodily orifices and cavities. Then again, other isolates are food-borne, with some of these having been long associated with naturally fermented cheeses and yogurts. Additionally, because of perceived health benefits to humans and animals, numerous L. rhamnosus strains have been selected for use as so-called probiotics and are often taken in the form of dietary supplements and functional foods. At the genome level, it is anticipated that certain genetic variances will have provided the niche-related phenotypes that augment the flexible adaptiveness of this species, thus enabling its strains to grow and survive in their respective host environments. For this present study, we considered it functionally informative to examine and catalogue the genotype-phenotype variation existing at the cell surface between different L. rhamnosus strains, with the presumption that this might be relatable to habitat preferences and ecological adaptability. Here, we conducted a pan-genomic study involving 13 genomes from L. rhamnosus isolates with various origins. In using a benchmark strain (gut-adapted L. rhamnosus GG) for our pan-genome comparison, we had focused our efforts on a detailed examination and description of gene products for certain functionally relevant surface-exposed proteins, each of which in effect might also play a part in niche adaptability among the other strains. Perhaps most significantly of the surface protein loci we had analyzed, it would appear that the spaCBA operon (known to encode SpaCBA-called pili having a mucoadhesive phenotype) is a genomic rarity and an uncommon occurrence in L. rhamnosus. However, for any of the so-piliated L. rhamnosus strains, they will likely possess an increased niche-specific fitness, which functionally might presumably be manifested by a protracted transient colonization of the gut mucosa or some similar microhabitat.

Assessment of aerobic and respiratory growth in the Lactobacillus casei group

One hundred eighty four strains belonging to the species Lactobacillus casei, L. paracasei and L. rhamnosus were screened for their ability to grow under aerobic conditions, in media containing heme and menaquinone and/or compounds generating reactive oxygen species (ROS), in order to identify respiratory and oxygen-tolerant phenotypes. Most strains were able to cope with aerobic conditions and for many strains aerobic growth and heme or heme/menaquinone supplementation increased biomass production compared to anaerobic cultivation. Only four L. casei strains showed a catalase-like activity under anaerobic, aerobic and respiratory conditions and were able to survive in presence of H2O2 (1 mM). Almost all L. casei and L. paracasei strains tolerated menadione (0.2 mM) and most tolerated pyrogallol (50 mM), while L. rhamnosus was usually resistant only to the latter compound. This is the first study in which an extensive screening of oxygen and oxidative stress tolerance of members of the L. casei group has been carried out. Results allowed the selection of strains showing the typical traits of aerobic and respiratory metabolism (increased pH and biomass under aerobic or respiratory conditions) and unique oxidative stress response properties. Aerobic growth and respiration may confer technological and physiological advantages in the L. casei group and oxygen-tolerant phenotypes could be exploited in several food industry applications.

Production of aroma compounds in lactic fermentations

This review describes recent scientific research on the production of aroma compounds by lactic acid bacteria (LAB) in fermented food products. We discuss the various precursor molecules for the formation of aroma compounds in connection with the metabolic pathways involved. The roles of nonmetabolic properties such as cell lysis are also described in relation to aroma formation. Finally, we provide an overview of the literature on methods to steer and control aroma formation by LAB in mixed culture fermentations. We demonstrate that the technological progress made recently in high-throughput analysis methods has been driving the development of new approaches to understand, control, and steer aroma formation in (dairy) fermentation processes. This currently entails proposing new rules for designing stable, high-performance mixed cultures constituting a selection of strains, which in concert and on the basis of their individual predicted gene contents deliver the required functionalities.

Comparative functional genomics of Lactobacillus spp. reveals possible mechanisms for specialization of vaginal lactobacilli to their environment

Lactobacilli are found in a wide variety of habitats. Four species, Lactobacillus crispatus, L. gasseri, L. iners, and L. jensenii, are common and abundant in the human vagina and absent from other habitats. These may be adapted to the vagina and possess characteristics enabling them to thrive in that environment. Furthermore, stable codominance of multiple Lactobacillus species in a single community is infrequently observed. Thus, it is possible that individual vaginal Lactobacillus species possess unique characteristics that confer to them host-specific competitive advantages. We performed comparative functional genomic analyses of representatives of 25 species of Lactobacillus, searching for habitat-specific traits in the genomes of the vaginal lactobacilli. We found that the genomes of the vaginal species were significantly smaller and had significantly lower GC content than those of the nonvaginal species. No protein families were found to be specific to the vaginal species analyzed, but some were either over- or underrepresented relative to nonvaginal species. We also found that within the vaginal species, each genome coded for species-specific protein families. Our results suggest that even though the vaginal species show no general signatures of adaptation to the vaginal environment, each species has specific and perhaps unique ways of interacting with its environment, be it the host or other microbes in the community. These findings will serve as a foundation for further exploring the role of lactobacilli in the ecological dynamics of vaginal microbial communities and their ultimate impact on host health.

Comparative genome analysis of Lactobacillus rhamnosus clinical isolates from initial stages of dental pulp infection: identification of a new exopolysaccharide cluster

The human oral microbiome has a major role in oral diseases including dental caries. Our studies on progression of caries infection through dentin and more recently, the invasion of vital dental pulp, detected Lactobacillus rhamnosus in the initial stages of infection of vital pulp tissue. In this study employing current high-throughput next generation sequencing technology we sought to obtain insight into genomic traits of tissue invasive L. rhamnosus, to recognise biomarkers that could provide an understanding of pathogenic potential of lactobacilli, generally regarded as safe. Roche GS FLX+ technology was used to generate whole genome sequences of two clinical isolates of L. rhamnosus infecting vital pulp. Detailed genome-wide comparison of the genetic profiles of tissue invasive L. rhamnosus with probiotic L. rhamnosus was performed to test the hypothesis that specific strains of L. rhamnosus possessing a unique gene complement are selected for the capacity to invade vital pulp tissue. Analysis identified 264 and 258 genes respectively, from dental pulp-invasive L. rhamnosus strains LRHMDP2 and LRHMDP3 isolated from two different subjects that were not present in the reference probiotic L. rhamnosus strain ATCC 53103 (GG). Distinct genome signatures identified included the presence of a modified exopolysaccharide cluster, a characteristic confirmed in a further six clinical isolates. Additional features of LRHMDP2 and LRHMDP3 were altered transcriptional regulators from RpoN, NtrC, MutR, ArsR and zinc-binding Cro/CI families, as well as changes in the two-component sensor kinase response regulator and ABC transporters for ferric iron. Both clinical isolates of L. rhamnosus contained a single SpaFED cluster, as in L. rhamnosus Lc705, instead of the two Spa clusters (SpaCBA and SpaFED) identified in L. rhamnosus ATCC 53103 (GG). Genomic distance analysis and SNP divergence confirmed a close relationship of the clinical isolates but segregation from the reference probiotic L. rhamnosus strain ATCC 53103 (GG).

Lactobacillus buchneri genotyping on the basis of clustered regularly interspaced short palindromic repeat (CRISPR) locus diversity

Clustered regularly interspaced short palindromic repeats (CRISPR) in combination with associated sequences (cas) constitute the CRISPR-Cas immune system, which uptakes DNA from invasive genetic elements as novel "spacers" that provide a genetic record of immunization events. We investigated the potential of CRISPR-based genotyping of Lactobacillus buchneri, a species relevant for commercial silage, bioethanol, and vegetable fermentations. Upon investigating the occurrence and diversity of CRISPR-Cas systems in Lactobacillus buchneri genomes, we observed a ubiquitous occurrence of CRISPR arrays containing a 36-nucleotide (nt) type II-A CRISPR locus adjacent to four cas genes, including the universal cas1 and cas2 genes and the type II signature gene cas9. Comparative analysis of CRISPR spacer content in 26 L. buchneri pickle fermentation isolates associated with spoilage revealed 10 unique locus genotypes that contained between 9 and 29 variable spacers. We observed a set of conserved spacers at the ancestral end, reflecting a common origin, as well as leader-end polymorphisms, reflecting recent divergence. Some of these spacers showed perfect identity with phage sequences, and many spacers showed homology to Lactobacillus plasmid sequences. Following a comparative analysis of sequences immediately flanking protospacers that matched CRISPR spacers, we identified a novel putative protospacer-adjacent motif (PAM), 5'-AAAA-3'. Overall, these findings suggest that type II-A CRISPR-Cas systems are valuable for genotyping of L. buchneri.

Genomic adaptation of the Lactobacillus casei group

Lactobacillus casei, L. paracasei, and L. rhamnosus form a closely related taxonomic group (Lactobacillus casei group) within the facultatively heterofermentative lactobacilli. Here, we report the complete genome sequences of L. paracasei JCM 8130 and L. casei ATCC 393, and the draft genome sequence of L. paracasei COM0101, all of which were isolated from daily products. Furthermore, we re-annotated the genome of L. rhamnosus ATCC 53103 (also known as L. rhamnosus GG), which we have previously reported. We confirmed that ATCC 393 is distinct from other strains previously described as L. paracasei. The core genome of 10 completely sequenced strains of the L. casei group comprised 1,682 protein-coding genes. Although extensive genome-wide synteny was found among the L. casei group, the genomes of ATCC 53103, JCM 8130, and ATCC 393 contained genomic islands compared with L. paracasei ATCC 334. Several genomic islands, including carbohydrate utilization gene clusters, were found at the same loci in the chromosomes of the L. casei group. The spaCBA pilus gene cluster, which was first identified in GG, was also found in other strains of the L. casei group, but several L. paracasei strains including COM0101 contained truncated spaC gene. ATCC 53103 encoded a higher number of proteins involved in carbohydrate utilization compared with intestinal lactobacilli, and extracellular adhesion proteins, several of which are absent in other strains of the L. casei group. In addition to previously fully sequenced L. rhamnosus and L. paracasei strains, the complete genome sequences of L. casei will provide valuable insights into the evolution of the L. casei group.

Comparative genomic and functional analysis of 100 Lactobacillus rhamnosus strains and their comparison with strain GG

Lactobacillus rhamnosus is a lactic acid bacterium that is found in a large variety of ecological habitats, including artisanal and industrial dairy products, the oral cavity, intestinal tract or vagina. To gain insights into the genetic complexity and ecological versatility of the species L. rhamnosus, we examined the genomes and phenotypes of 100 L. rhamnosus strains isolated from diverse sources. The genomes of 100 L. rhamnosus strains were mapped onto the L. rhamnosus GG reference genome. These strains were phenotypically characterized for a wide range of metabolic, antagonistic, signalling and functional properties. Phylogenomic analysis showed multiple groupings of the species that could partly be associated with their ecological niches. We identified 17 highly variable regions that encode functions related to lifestyle, i.e. carbohydrate transport and metabolism, production of mucus-binding pili, bile salt resistance, prophages and CRISPR adaptive immunity. Integration of the phenotypic and genomic data revealed that some L. rhamnosus strains possibly resided in multiple niches, illustrating the dynamics of bacterial habitats. The present study showed two distinctive geno-phenotypes in the L. rhamnosus species. The geno-phenotype A suggests an adaptation to stable nutrient-rich niches, i.e. milk-derivative products, reflected by the alteration or loss of biological functions associated with antimicrobial activity spectrum, stress resistance, adaptability and fitness to a distinctive range of habitats. In contrast, the geno-phenotype B displays adequate traits to a variable environment, such as the intestinal tract, in terms of nutrient resources, bacterial population density and host effects.

Some bacterial species are specialists and adapted to a single niche, while others are generalists and able to grow in various environmental conditions. Lactobacillus rhamnosus is a generalist and its members can often be found in different human cavities but also in various artisanal and industrial dairy products. To gain insights into the genetic complexity and ecological versatility of this species, we collected 100 L. rhamnosus strains from different niches. Genomic and functional analysis of these revealed a dichotomy within the species that reflected its adaptation to particular niches. The variable regions identified in the L. rhamnosus genome encode lifestyle traits that allowed us to demonstrate that some L. rhamnosus isolates possibly resided in multiple habitats. Our work brings valuable data on the ecological dynamics and adaptability of the species and provides a basis for a model explaining the ecology of L. rhamnosus in an anthropocentric perspective. Finally, we observed that a set of pheno-genomic markers, i.e. CRISPR oligotyping or carbohydrate metabolism, would be sufficient and among the best ways to differentiate the L. rhamnosus strains, providing a general approach to select the highest diversity in these and other bacterial species.

Comparative genome analysis of Lactobacillus casei strains isolated from Actimel and Yakult products reveals marked similarities and points to a common origin

The members of the Lactobacillus genus are widely used in the food and feed industry and show a remarkable ecological adaptability. Several Lactobacillus strains have been marketed as probiotics as they possess health-promoting properties for the host. In the present study, we used two complementary next-generation sequencing technologies to deduce the genome sequences of two Lactobacillus casei strains LcA and LcY, which were isolated from the products Actimel and Yakult, commercialized as probiotics. The LcA and LcY draft genomes have, respectively, an estimated size of 3067 and 3082 Mb and a G+C content of 46.3%. Both strains are close to identical to each other and differ by no more than minor chromosomal re-arrangements, substitutions, insertions and deletions, as evident from the verified presence of one insertion-deletion (InDel) and only 29 single-nucleotide polymorphisms (SNPs). In terms of coding capacity, LcA and LcY are predicted to encode a comparable exoproteome, indicating that LcA and LcY are likely to establish similar interactions with human intestinal cells. Moreover, both L. casei LcA and LcY harboured a 59.6 kb plasmid that shared high similarities with plasmids found in other L. casei strains, such as W56 and BD-II. Further analysis revealed that the L. casei plasmids constitute a good evolution marker within the L. casei species. The plasmids of the LcA and LcY strains are almost identical, as testified by the presence of only three verified SNPs, and share a 3.5 kb region encoding a remnant of a lactose PTS system that is absent from the plasmids of W56 and BD-II but conserved in another smaller L. casei plasmid (pLC2W). Our observations imply that the results obtained in animal and human experiments performed with the Actimel and Yakult strains can be compared with each other as these strains share a very recent common ancestor.

Funding Information The present work was supported by the Center of Excellence in Microbial Food Safety Research (Academy of Finland, Grant 141140), Grant ERC 250172 – Microbes Inside from the European Research Council and Grants 137389 and 141130 from the Academy of Finland. F.P.D. was funded by a postdoctoral research fellowship (Academy of Finland, Grant 252123).