Complete genome sequence of Leuconostoc gelidum subsp. gasicomitatum KG16-1, isolated from vacuum-packaged vegetable sausages

Pangenome and genomic taxonomy analyses of Leuconostoc gelidum and Leuconostoc gasicomitatum

BACKGROUND

Leuconostoc gelidum and Leuconostoc gasicomitatum have dual roles in foods. They may spoil cold-stored packaged foods but can also be beneficial in kimchi fermentation. The impact in food science as well as the limited number of publicly available genomes prompted us to create pangenomes and perform genomic taxonomy analyses starting from de novo sequencing of the genomes of 37 L. gelidum/L. gasicomitatum strains from our culture collection. Our aim was also to evaluate the recently proposed change in taxonomy as well as to study the genomes of strains with different lifestyles in foods.

METHODS

We selected as diverse a set of strains as possible in terms of sources, previous genotyping results and geographical distribution, and included also 10 publicly available genomes in our analyses. We studied genomic taxonomy using pairwise average nucleotide identity (ANI) and calculation of digital DNA-DNA hybridisation (dDDH) scores. Phylogeny analyses were done using the core gene set of 1141 single-copy genes and a set of housekeeping genes commonly used for lactic acid bacteria. In addition, the pangenome and core genome sizes as well as some properties, such as acquired antimicrobial resistance (AMR), important due to the growth in foods, were analysed.

RESULTS

Genome relatedness indices and phylogenetic analyses supported the recently suggested classification that restores the taxonomic position of L. gelidum subsp. gasicomitatum back to the species level as L. gasicomitatum. Genome properties, such as size and coding potential, revealed limited intraspecies variation and showed no attribution to the source of isolation. The distribution of the unique genes between species and subspecies was not associated with the previously documented lifestyle in foods. None of the strains carried any acquired AMR genes or genes associated with any known form of virulence.

CONCLUSION

Genome-wide examination of strains confirms that the proposition to restore the taxonomic position of L. gasicomitatum is justified. It further confirms that the distribution and lifestyle of L. gelidum and L. gasicomitatum in foods have not been driven by the evolution of functional and phylogenetic diversification detectable at the genome level.

Bacterial Communities and Prediction of Microbial Metabolic Pathway in Rice Wine Koji From Different Regions in China

Rice wine koji, a traditional homemade starter culture in China, is nutritious and delicious. The final quality of rice wine koji is closely related to the structure of its microbial community. However, the diversity of natural microorganisms in rice wine koji from different regions has not been evaluated. In this study, the microbial population of 92 naturally fermented rice koji samples collected from Hubei, Guangxi, and Sichuan was systematically analyzed by high-throughput sequencing. From all the rice wine koji samples, 22 phyla and 479 bacterial genera were identified. Weissella, Pediococcus, Lactobacillus, Enterobacter, Lactococcus, Pantoea, Bacillus, Staphylococcus, and Leuconostoc were the dominant genera in rice wine koji. The bacterial community structure of rice wine koji samples from different regions was significantly different (p < 0.05). The bacterial community composition of the samples from Hubei and Guangxi was similar, but significantly different from that of SC samples (p < 0.05). These differences may be caused by variations in geography, environment, or manufacturing. In addition, the results of microbial phenotype prediction by BugBase and bacterial functional potential prediction by PICRUSt showed that eight of the nine predicted phenotypic functions of rice wine koji samples from different regions were significantly different (p < 0.05) and that vigorous bacterial metabolism occurred in rice wine koji samples.

Transcriptomic time-series analysis of cold- and heat-shock response in psychrotrophic lactic acid bacteria

BACKGROUND

Psychrotrophic lactic acid bacteria (LAB) species are the dominant species in the microbiota of cold-stored modified-atmosphere-packaged food products and are the main cause of food spoilage. Despite the importance of psychrotrophic LAB, their response to cold or heat has not been studied. Here, we studied the transcriptome-level cold- and heat-shock response of spoilage lactic acid bacteria with time-series RNA-seq for Le. gelidum, Lc. piscium, and P. oligofermentans at 0 °C, 4 °C, 14 °C, 25 °C, and 28 °C.

RESULTS

We observed that the cold-shock protein A (cspA) gene was the main cold-shock protein gene in all three species. Our results indicated that DEAD-box RNA helicase genes (cshA, cshB) also play a critical role in cold-shock response in psychrotrophic LAB. In addition, several RNase genes were involved in cold-shock response in Lc. piscium and P. oligofermentans. Moreover, gene network inference analysis provided candidate genes involved in cold-shock response. Ribosomal proteins, tRNA modification, rRNA modification, and ABC and efflux MFS transporter genes clustered with cold-shock response genes in all three species, indicating that these genes could be part of the cold-shock response machinery. Heat-shock treatment caused upregulation of Clp protease and chaperone genes in all three species. We identified transcription binding site motifs for heat-shock response genes in Le. gelidum and Lc. piscium. Finally, we showed that food spoilage-related genes were upregulated at cold temperatures.

CONCLUSIONS

The results of this study provide new insights on the cold- and heat-shock response of psychrotrophic LAB. In addition, candidate genes involved in cold- and heat-shock response predicted using gene network inference analysis could be used as targets for future studies.

Comparative Genomics of Streptococcus thermophilus Support Important Traits Concerning the Evolution, Biology and Technological Properties of the Species

Streptococcus thermophilus is a major starter for the dairy industry with great economic importance. In this study we analyzed 23 fully sequenced genomes of S. thermophilus to highlight novel aspects of the evolution, biology and technological properties of this species. Pan/core genome analysis revealed that the species has an important number of conserved genes and that the pan genome is probably going to be closed soon. According to whole genome phylogeny and average nucleotide identity (ANI) analysis, most S. thermophilus strains were grouped in two major clusters (i.e., clusters A and B). More specifically, cluster A includes strains with chromosomes above 1.83 Mbp, while cluster B includes chromosomes below this threshold. This observation suggests that strains belonging to the two clusters may be differentiated by gene gain or gene loss events. Furthermore, certain strains of cluster A could be further subdivided in subgroups, i.e., subgroup I (ASCC 1275, DGCC 7710, KLDS SM, MN-BM-A02, and ND07), II (MN-BM-A01 and MN-ZLW-002), III (LMD-9 and SMQ-301), and IV (APC151 and ND03). In cluster B certain strains formed one distinct subgroup, i.e., subgroup I (CNRZ1066, CS8, EPS, and S9). Clusters and subgroups observed for S. thermophilus indicate the existence of lineages within the species, an observation which was further supported to a variable degree by the distribution and/or the architecture of several genomic traits. These would include exopolysaccharide (EPS) gene clusters, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs)-CRISPR associated (Cas) systems, as well as restriction-modification (R-M) systems and genomic islands (GIs). Of note, the histidine biosynthetic cluster was found present in all cluster A strains (plus strain NCTC12958^T) but was absent from all strains in cluster B. Other loci related to lactose/galactose catabolism and urea metabolism, aminopeptidases, the majority of amino acid and peptide transporters, as well as amino acid biosynthetic pathways were found to be conserved in all strains suggesting their central role for the species. Our study highlights the necessity of sequencing and analyzing more S. thermophilus complete genomes to further elucidate important aspects of strain diversity within this starter culture that may be related to its application in the dairy industry.

Deciphering intra-species bacterial diversity of meat and seafood spoilage microbiota using gyrB amplicon sequencing: A comparative analysis with 16S rDNA V3-V4 amplicon sequencing

Meat and seafood spoilage ecosystems harbor extensive bacterial genomic diversity that is mainly found within a small number of species but within a large number of strains with different spoilage metabolic potential. To decipher the intraspecies diversity of such microbiota, traditional metagenetic analysis using the 16S rRNA gene is inadequate. We therefore assessed the potential benefit of an alternative genetic marker, gyrB, which encodes the subunit B of DNA gyrase, a type II DNA topoisomerase. A comparison between 16S rDNA-based (V3-V4) amplicon sequencing and gyrB-based amplicon sequencing was carried out in five types of meat and seafood products, with five mock communities serving as quality controls. Our results revealed that bacterial richness in these mock communities and food samples was estimated with higher accuracy using gyrB than using16S rDNA. However, for Firmicutes species, 35% of putative gyrB reads were actually identified as sequences of a gyrB paralog, parE, which encodes subunit B of topoisomerase IV; we therefore constructed a reference database of published sequences of both gyrB and pare for use in all subsequent analyses. Despite this co-amplification, the deviation between relative sequencing quantification and absolute qPCR quantification was comparable to that observed for 16S rDNA for all the tested species. This confirms that gyrB can be used successfully alongside 16S rDNA to determine the species composition (richness and evenness) of food microbiota. The major benefit of gyrB sequencing is its potential for improving taxonomic assignment and for further investigating OTU richness at the subspecies level, thus allowing more accurate discrimination of samples. Indeed, 80% of the reads of the 16S rDNA dataset were represented by thirteen 16S rDNA-based OTUs that could not be assigned at the species-level. Instead, these same clades corresponded to 44 gyrB-based OTUs, which differentiated various lineages down to the subspecies level. The increased ability of gyrB-based analyses to track and trace phylogenetically different groups of strains will generate improved resolution and more reliable results for studies of the strains implicated in food processes.

Food Spoilage-Associated Leuconostoc, Lactococcus, and Lactobacillus Species Display Different Survival Strategies in Response to Competition

Psychrotrophic lactic acid bacteria (LAB) are the prevailing spoilage organisms in packaged cold-stored meat products. Species composition and metabolic activities of such LAB spoilage communities are determined by the nature of the meat product, storage conditions, and interspecies interactions. Our knowledge of system level responses of LAB during such interactions is very limited. To expand it, we studied interactions between three common psychrotrophic spoilage LAB (Leuconostoc gelidum, Lactococcus piscium, and Lactobacillus oligofermentans) by comparing their time course transcriptome profiles obtained during their growth in individual, pairwise, and triple cultures. The study revealed how these LAB employed different strategies to cope with the consequences of interspecies competition. The fastest-growing bacterium, Le. gelidum, attempted to enhance its nutrient-scavenging and growth capabilities in the presence of other LAB through upregulation of carbohydrate catabolic pathways, pyruvate fermentation enzymes, and ribosomal proteins, whereas the slower-growing Lc. piscium and Lb. oligofermentans downregulated these functions. These findings may explain the competitive success and predominance of Le. gelidum in a variety of spoiled foods. Peculiarly, interspecies interactions induced overexpression of prophage genes and restriction modification systems (mechanisms of DNA exchange and protection against it) in Lc. piscium and Lb. oligofermentans but not in Le. gelidum Cocultivation induced also overexpression of the numerous putative adhesins in Lb. oligofermentans These adhesins might contribute to the survival of this slowly growing bacterium in actively growing meat spoilage communities.IMPORTANCE Despite the apparent relevance of LAB for biotechnology and human health, interactions between members of LAB communities are not well known. Knowledge of such interactions is crucial for understanding how these communities function and, consequently, whether there is any possibility to develop new strategies to interfere with their growth and to postpone spoilage of packaged and refrigerated foods. With the help of controlled experiments, detailed regulation events can be observed. This study gives an insight into the system level interactions and the different competition-induced survival strategies related to enhanced uptake and catabolism of carbon sources, overexpression of adhesins and putative bacteriocins, and the induction of exchange of genetic material. Even though this experiment dealt with only three LAB strains in vitro, these findings agreed well with the relative abundance patterns typically reported for these species in natural food microbial communities.

Complete genome sequence of Leuconostoc suionicum DSM 20241T provides insights into its functional and metabolic features

The genome of Leuconostoc suionicum DSM 20241^T (=ATCC 9135^T = LMG 8159^T = NCIMB 6992^T) was completely sequenced and its fermentative metabolic pathways were reconstructed to investigate the fermentative properties and metabolites of strain DSM 20241^T during fermentation. The genome of L. suionicum DSM 20241^T consists of a circular chromosome (2026.8 Kb) and a circular plasmid (21.9 Kb) with 37.58% G + C content, encoding 997 proteins, 12 rRNAs, and 72 tRNAs. Analysis of the metabolic pathways of L. suionicum DSM 20241^T revealed that strain DSM 20241^T performs heterolactic acid fermentation and can metabolize diverse organic compounds including glucose, fructose, galactose, cellobiose, mannose, sucrose, trehalose, arbutin, salcin, xylose, arabinose and ribose.

Metagenomic approach reveals microbial diversity and predictive microbial metabolic pathways in Yucha, a traditional Li fermented food

Yucha is a typical traditional fermented food of the Li population in the Hainan province of China, and it is made up of cooked rice and fresh fish. In the present study, metagenomic approach and culture-dependent technology were applied to describe the diversity of microbiota and identify beneficial microbes in the Yucha. At the genus level, Lactobacillus was the most abundant genus (43.82% of the total reads), followed by Lactococcus, Enterococcus, Vibrio, Weissella, Pediococcus, Enterobacter, Salinivibrio, Acinetobacter, Macrococcus, Kluyvera and Clostridium; this result was confirmed by q-PCR. PCoA based on Weighted UniFrac distances showed an apparent clustering pattern for Yucha samples from different locations, and Lactobacillus sakei, Lactobacillus saniviri and Staphylococcus sciuri represented OTUs according to the major identified markers. At the microbial functional level, it was observed that there was an enrichment of metabolic functional features, including amino acid and carbohydrate metabolism, which implied that the microbial metabolism in the Yucha samples tended to be vigorous. Accordingly, we further investigated the correlation between the predominant microbes and metabolic functional features. Thirteen species of Lactobacillus (147 strains) were isolated, and Lactobacillus plantarum (60 isolates) and Lactobacillus pentosus (34 isolates) were isolated from every sample.