Enterococcus faecalis is widely studied as a common gut commensal and a nosocomial pathogen. In fact, Enterococcus faecalis is ubiquitous in nature, and it has been isolated from various niches, including the gastrointestinal tract, faeces, blood, urine, water, and fermented foods (such as dairy products). In order to elucidate the role of habitat in shaping the genome of Enterococcus faecalis, we performed a comparative genomic analysis of 78 strains of various origins. Although no correlation was found between the strain isolation habitat and the phylogeny of Enterococcus faecalis from our whole genome-based phylogenetic analysis, our results revealed some environment-associated features in the analysed Enterococcus faecalis genomes. Significant differences were found in the genome size and the number of predicted open reading frames (ORFs) between strains originated from different environments. In general, strains from water sources had the smallest genome size and the least number of predicted ORFs. We also identified 293 environment-specific genes, some of which might link to the adaptive strategies for survival in particular environments. In addition, the number of antibiotic resistance genes was significantly different between strains isolated from dairy products, water, and blood. Strains isolated from blood had the largest number of antibiotic resistance genes. These findings improve our understanding of the role of habitat in shaping the genomes of Enterococcus faecalis.

As an important nosocomial pathogen, Enterococcus faecium has received increasing attention in recent years. However, a large number of studies have focused on the hospital-associated isolates and ignored isolates originated from the natural environments. In this study, comparative genomic analysis was conducted on 161 isolates originated from human, animal, and naturally fermented dairy products. The results showed that the environment played an important role in shaping the genomes of Enterococcus faecium. The isolates from human had the largest average genome size, while the isolates from dairy products had the smallest average genome size and fewest antibiotic resistance genes. A phylogenetic tree was reconstructed based on the genomes of these isolates, which revealed new insights into the phylogenetic relationships among the dairy isolates and those from hospitals, communities, and animals. Furthermore, 202 environment-specific genes were identified, including 136 dairy-specific, 31 human blood-specific, and 35 human gastrointestinal-specific genes. Interestingly, five dairy-specific genes (namely lacF, lacA/B, lacD, lacG, and lacC) that constituted an integrated lactose metabolism pathway existed in almost all dairy isolates. The pathway conservation demonstrated an active role of the environment in shaping the genomes of Enterococcus faecium. This study shows that the Enterococcus faecium species has great genomic plasticity and high versatility to occupy broad ecological roles, dwelling as non-harmful dairy and animal gut commensals as well as significant nosocomial pathogens that disseminate antibiotic resistance genes.

Lactobacillus delbrueckii ssp. bulgaricus is one of the most commonly used starter cultures for yogurt production. However, how its genetic background affects acid production capacity is unclear. This study investigated the industrial potential of L. delbrueckii ssp. bulgaricus using population genomics and GWAS analysis. To meet our goal, population genetics and functional genomics analyses were performed on 188 newly sequenced L. delbrueckii ssp. bulgaricus strains isolated from naturally fermented dairy products together with 19 genome sequences retrieved from the NCBI database. Four distinct clusters were identified, and they were correlated with the geographical sites where the samples were collected. The GWAS analysis about acidification fermentation results with sucrose-fortified reconstituted skim milk revealed a significant association between l-lactate dehydrogenase (lldD; Ldb2036) and the bacterial acid production rate. Our study has broadened the understanding of the population structure and genetic diversity of L. delbrueckii ssp. bulgaricus by identifying potential targets for further research, development, and use of lactic acid bacteria in the dairy industry.

Background: The two subspecies of Bifidobacterium catenulatum, B. catenulatum subsp. kashiwanohense and B. catenulatum subsp. catenulatum, are usually from the infant and adult gut, respectively. However, the genomic analysis of their functional difference and genetic divergence has been rare. Here, 16 B. catenulatum strains, including 2 newly sequenced strains, were analysed through comparative genomics. Results: A phylogenetic tree based on 785 core genes indicated that the two subspecies of B. catenulatum were significantly separated. The comparison of genomic characteristics revealed that the two subspecies had significantly different genomic sizes (p < 0.05) but similar GC contents. The functional comparison revealed the most significant difference in genes of carbohydrate utilisation. Carbohydrate-active enzymes (CAZyme) present two clustering patterns in B. catenulatum. The B. catenulatum subsp. kashiwanohense specially including the glycoside hydrolases 95 (GH95) and carbohydrate-binding modules 51 (CBM51) families involved in the metabolism of human milk oligosaccharides (HMO) common in infants, also, the corresponding fucosylated HMO gene clusters were detected. Meanwhile, B. catenulatum subsp. catenulatum rich in GH3 may metabolise more plant-derived glycan in the adult intestine. Conclusions: These findings provide genomic evidence of carbohydrate utilisation bias, which may be a key cause of the genetic divergence of two B. catenulatum subspecies.

In the ecosystem of spontaneously fermented cow milk, the characteristics and relationship of bacterial communities and nonvolatile components at different scales of geographical distances (provincial, county, and village levels) are unclear. Here, 25 sampling sites from Xin Jiang and Tibet, 2 provinces of China, were selected based on the distribution of spontaneously fermented cow milk and used for metagenomic and metabolomic analysis. At the provincial geographical distance, the same predominant species, Lactobacillus delbrueckii ssp. bulgaricus and Streptococcus thermophilus, were detected in Xin Jiang and Tibet. Further, the richness of the bacterial composition of samples from Tibet was higher than those from Xin Jiang; specifically, at the species level, 28 species were identified in Tibet samples but only 7 species in Xin Jiang samples. At the provincial geographical level, we detected significant differences in bacterial structure, shown in principal coordinate analysis plots, and significant differences (Simpson index) in bacterial diversity were also detected. However, at the county and village levels, no significant differences were detected in bacterial communities and diversity, but a difference in bacterial compositions was detectable. This indicates that bacterial communities and diversity of spontaneously fermented milk dissimilarity significantly increased with geographic distance. For the nonvolatile component profiles, the partial least squares discriminant analysis plot (R²Y > 0.5 and Q² > 0.5 for the goodness-of-fit and predictive ability parameter, respectively) showed that samples from different geographical distances (provincial, county, and village) were all separated, which indicated that all the discriminations in nonvolatile components profiles were from different geographical distances. Investigating relationships between lactic acid bacteria and discriminatory nonvolatile components at the county level showed that 9 species were positively correlated with 16 discriminatory nonvolatile components, all species with low abundance rather than the predominant species L. delbrueckii ssp. bulgaricus and Strep. thermophilus, which indicates the importance of the selection of autochthonous nonpredominant bacteria.