Multilocus sequence typing of Leuconostoc gelidum subsp. gasicomitatum, a Psychrotrophic lactic acid bacterium causing spoilage of packaged perishable foods

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.

Multilocus Sequence Typing of Leuconostoc mesenteroides Strains From the Qinghai-Tibet Plateau

Leuconostoc mesenteroides strains were a type of epiphytic bacterium widely used in fermented foods and products in the biochemical and pharmaceutical industries but data on its presence in foods from Qinghai-Tibet Plateau in China was scarce. In this study, molecular analysis based on multilocus sequence typing (MLST) with eight housekeeping genes (pyrG, groeL, rpoB, recA, uvrC, murC, carB, and pheS) was carried out on 45 L. mesenteroides strains isolated from different plants and dairy products from Qinghai-Tibet Plateau in China. The objective of this study was to perform genetic diversity analysis and explore the relationship between strains and isolate samples or separate regions. A total of 25 sequence types (STs) were identified with a diversity of up to 55.6%, which were grouped into one clonal complexes (CCs), 3 doublets and 17 singletons by eBURST. The results of minimum spanning tree and clustering analysis indicated these L. mesenteroides strains from the Qinghai-Tibet Plateau were relatively weakly related to the isolated region. However, there was a close relationship between the genotypes of L. mesenteroides strains and the type of the isolated sample, which was consistent with the results of API 50CH. The MLST scheme presented in this study provides a shareable and comparable sequence database and enhances our knowledge of the population diversity of L. mesenteroides strains which will be further used for the selection of industrial strains.

Microbiomes Associated With Foods From Plant and Animal Sources

Food microbiome composition impacts food safety and quality. The resident microbiota of many food products is influenced throughout the farm to fork continuum by farming practices, environmental factors, and food manufacturing and processing procedures. Currently, most food microbiology studies rely on culture-dependent methods to identify bacteria. However, advances in high-throughput DNA sequencing technologies have enabled the use of targeted 16S rRNA gene sequencing to profile complex microbial communities including non-culturable members. In this study we used 16S rRNA gene sequencing to assess the microbiome profiles of plant and animal derived foods collected at two points in the manufacturing process; post-harvest/pre-retail (cilantro) and retail (cilantro, masala spice mixes, cucumbers, mung bean sprouts, and smoked salmon). Our findings revealed microbiome profiles, unique to each food, that were influenced by the moisture content (dry spices, fresh produce), packaging methods, such as modified atmospheric packaging (mung bean sprouts and smoked salmon), and manufacturing stage (cilantro prior to retail and at retail). The masala spice mixes and cucumbers were comprised mainly of Proteobacteria, Firmicutes, and Actinobacteria. Cilantro microbiome profiles consisted mainly of Proteobacteria, followed by Bacteroidetes, and low levels of Firmicutes and Actinobacteria. The two brands of mung bean sprouts and the three smoked salmon samples differed from one another in their microbiome composition, each predominated by either by Firmicutes or Proteobacteria. These data demonstrate diverse and highly variable resident microbial communities across food products, which is informative in the context of food safety, and spoilage where indigenous bacteria could hamper pathogen detection, and limit shelf life.

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.

Fragment of the aspartyl-tRNA synthetase applicable as a shared classification and phylogenetic marker in particular representatives of the order Lactobacillales

The order Lactobacillales represents a morphologically, metabolically, and physiologically diverse group of bacteria. Lactic acid bacteria represent the core of this phylogenetic group. They are a part of epiphytic microflora, fermented dairy, meat, fruit and vegetable products, and the digestive tract of humans and animals. Despite the fact that these bacteria form a phenotypically and genotypically heterogeneous group, their phylogenetic relationship enables to propose a common genetic marker usable in classification, typing, and phylogeny. By creation of consensus sequence based on available genomic sequences of some representatives of order Lactobacillales, a specific primer-pair binding variable region of aspS gene (length of 615 nts) encoding the aspartyl-tRNA synthetase was designed. This gene has not yet been used in classification and phylogeny of the order Lactobacillales, although it meets the requirements of molecular markers (distribution and single copy in bacterial genomes, functional constancy and genetic stability, sequence variability among taxonomic units, irreplaceable role in proteosynthesis). Primers were applied on 54 type and wild Lactobacillales strains. Obtained sequences allowed to provide alignments for purpose of phylogenetic tree reconstructions that uncovered particular phylogenetic clusters of vagococci/enterococci, obligately homofermentative and heterofermentative lactobacilli. Although a relatively short fragment of the aspS gene (approximately 33% of the complete gene sequence) was evaluated, much higher sequence variability (61.8% of pairwise identity) among strains examined compared with 16S rRNA gene (90.7%, length of 1318 nt) provides a relatively simple and effective tool for classification and typing of selected representatives of the order Lactobacillales.

Effective identification of Lactobacillus casei group species: genome-based selection of the gene mutL as the target of a novel multiplex PCR assay

Lactobacillus casei,Lactobacillus paracasei and Lactobacillusrhamnosus form a closely related taxonomic group (the L. casei group) within the facultatively heterofermentative lactobacilli. Strains of these species have been used for a long time as probiotics in a wide range of products, and they represent the dominant species of nonstarter lactic acid bacteria in ripened cheeses, where they contribute to flavour development. The close genetic relationship among those species, as well as the similarity of biochemical properties of the strains, hinders the development of an adequate selective method to identify these bacteria. Despite this being a hot topic, as demonstrated by the large amount of literature about it, the results of different proposed identification methods are often ambiguous and unsatisfactory. The aim of this study was to develop a more robust species-specific identification assay for differentiating the species of the L. casei group. A taxonomy-driven comparative genomic analysis was carried out to select the potential target genes whose similarity could better reflect genome-wide diversity. The gene mutL appeared to be the most promising one and, therefore, a novel species-specific multiplex PCR assay was developed to rapidly and effectively distinguish L. casei, L. paracasei and L. rhamnosus strains. The analysis of a collection of 76 wild dairy isolates, previously identified as members of the L. casei group combining the results of multiple approaches, revealed that the novel designed primers, especially in combination with already existing ones, were able to improve the discrimination power at the species level and reveal previously undiscovered intraspecific biodiversity.

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

Leuconostoc gelidum subsp. gasicomitatum is a predominant lactic acid bacterium (LAB) in spoilage microbial communities of different kinds of modified-atmosphere packaged (MAP) food products. So far, only one genome sequence of a poultry-originating type strain of this bacterium (LMG 18811(T)) has been available. In the current study, we present the completely sequenced and functionally annotated genome of strain KG16-1 isolated from a vegetable-based product. In addition, six other vegetable-associated strains were sequenced to study possible "niche" specificity suggested by recent multilocus sequence typing. The genome of strain KG16-1 consisted of one circular chromosome and three plasmids, which together contained 2,035 CDSs. The chromosome carried at least three prophage regions and one of the plasmids encoded a galactan degradation cluster, which might provide a survival advantage in plant-related environments. The genome comparison with LMG 18811(T) and six other vegetable strains suggests no major differences between the meat- and vegetable-associated strains that would explain their "niche" specificity. Finally, the comparison with the genomes of other leuconostocs highlights the distribution of functionally interesting genes across the L. gelidum strains and the genus Leuconostoc.

Emergence of a New Population of Rathayibacter toxicus: An Ecologically Complex, Geographically Isolated Bacterium

Rathayibacter toxicus is a gram-positive bacterium that infects the floral parts of several Poaceae species in Australia. Bacterial ooze is often produced on the surface of infected plants and bacterial galls are produced in place of seed. R. toxicus is a regulated plant pathogen in the U.S. yet reliable detection and diagnostic tools are lacking. To better understand this geographically-isolated plant pathogen, genetic variation as a function of geographic location, host species, and date of isolation was determined for isolates collected over a forty-year period. Discriminant analyses of recently collected and archived isolates using Multi-Locus Sequence Typing (MLST) and Inter-Simple Sequence Repeats (ISSR) identified three populations of R. toxicus; RT-I and RT-II from South Australia and RT-III from Western Australia. Population RT-I, detected in 2013 and 2014 from the Yorke Peninsula in South Australia, is a newly emerged population of R. toxicus not previously reported. Commonly used housekeeping genes failed to discriminate among the R. toxicus isolates. However, strategically selected and genome-dispersed MLST genes representing an array of cellular functions from chromosome replication, antibiotic resistance and biosynthetic pathways to bacterial acquired immunity were discriminative. Genetic variation among isolates within the RT-I population was less than the within-population variation for the previously reported RT-II and RT-III populations. The lower relative genetic variation within the RT-I population and its absence from sampling over the past 40 years suggest its recent emergence. RT-I was the dominant population on the Yorke Peninsula during the 2013–2014 sampling period perhaps indicating a competitive advantage over the previously detected RT-II population. The potential for introduction of this bacterial plant pathogen into new geographic areas provide a rationale for understanding the ecological and evolutionary trajectories of R. toxicus.