Comparative Genomics Analysis of Habitat Adaptation by Lactobacillus kefiranofaciens

Lactobacillus kefiranofaciens is often found in fermented dairy products. Many strains of this species have probiotic properties, contributing to the regulation of immune metabolism and intestinal flora. This species was added to the list of lactic acid bacteria that can be added to food in China, in 2020. However, research on the genomics of this species is scarce. In this study we undertook whole genome sequencing analysis of 82 strains of L. kefiranofaciens from different habitats, of which 9 strains were downloaded from the NCBI RefSeq (National Center for Biotechnology Information RefSeq). The mean genome size of the 82 strains was 2.05 ± 0.25 Mbp, and the mean DNA G + C content was 37.47 ± 0.42%. The phylogenetic evolutionary tree for the core genes showed that all strains belonged to five clades with clear aggregation in relation to the isolation habitat; this indicated that the genetic evolution of L. kefiranofaciens was correlated to the isolation habitat. Analysis of the annotation results identified differences in the functional genes, carbohydrate active enzymes (CAZy) and bacteriocins amongst different isolated strains, which were related to the environment. Isolates from kefir grains had more enzymes for cellulose metabolism and a better ability to use vegetative substrates for fermentation, which could be used in feed production. Isolates from kefir grains also had fewer kinds of bacteriocin than isolates from sour milk and koumiss; helveticin J and lanthipeptide class I were not found in the isolates from kefir grains. The genomic characteristics and evolutionary process of L. kefiranofaciens were analyzed by comparative genomics and this paper explored the differences in the functional genes amongst the strains, aiming to provide a theoretical basis for the research and development of L. kefiranofaciens.