Comparative Genomics of Completely Sequenced Lactobacillus helveticus Genomes Provides Insights into Strain-Specific Genes and Resolves Metagenomics Data Down to the Strain Level

Microbial composition and viability of natural whey starters used in PDO Comté cheese-making

Natural whey starters (NWS) are cultures with undefined multiple-strains species commonly used to speed up the fermentation process of cheeses. The aim of this study was to explore the diversity and the viability of Comté cheese NWS microbiota. Culture-dependent methods, i.e. plate counting and genotypic characterization, and culture-independent methods, i.e. qPCR, viability-qPCR, fluorescence microscopy and DNA metabarcoding, were combined to analyze thirty-six NWS collected in six Comté cheese factories at two seasons. Our results highlighted that NWS were dominated by Streptococcus thermophilus (ST) and thermophilic lactobacilli. These species showed a diversity of strains based on Rep-PCR. The dominance of Lactobacillus helveticus (LH) over Lactobacillus delbrueckii (LD) varied depending on the factory and the season. This highlighted two types of NWS: the type-ST/LD (LD > LH) and the type-ST/LH (LD < LH). The microbial composition varied depending on cheese factory. One factory was distinguished by its level of culturable microbial groups (ST, enterococci and yeast) and its fungi diversity. The approaches used to estimate the viability showed that most NWS cells were viable. Further investigations are needed to understand the microbial diversity of these NWS.

Intraspecific microdiversity and ecological drivers of lactic acid bacteria in naturally fermented milk ecosystem

Traditional fermented milks are produced by inoculating technique, which selects well-adapted microorganisms that have been passed on through generations. Few reports have used naturally fermented milks as model ecosystems to investigate the mechanism of formation of intra-species microbial diversity. Here, we isolated and whole-genome-sequenced a total of 717 lactic acid bacterial isolates obtained from 12 independent naturally fermented milks collect from 12 regions across five countries. We further analyzed the within-sample intra-species phylogenies of 214 Lactobacillus helveticus isolates, 97 Lactococcus lactis subsp. lactis isolates, and 325 Lactobacillus delbrueckii subsp. bulgaricus isolates. We observed a high degree of intra-species genomic and functional gene diversity within-/between-sample(s). Single nucleotide polymorphism-based phylogenetic reconstruction revealed great within-sample intra-species heterogeneity, evolving from multiple lineages. Further phylogenetic reconstruction (presence-absence gene profile) revealed within-sample inter-clade functional diversity (based on carbohydrate-active enzyme- and peptidase-encoding genes) in all three investigated species/subspecies. By identifying and mapping clade-specific genes of intra-sample clades of the three species/subspecies to the respective fermented milk metagenome, we found extensive potential inter-/intra-species horizontal gene transfer events. Finally, the microbial composition of the samples is closely linked to the nucleotide diversity of the respective species/subspecies. Overall, our results contribute to the conservation of lactic acid bacteria resources, providing ecological insights into the microbial ecosystem of naturally fermented dairy products.

Assessment of beneficial effects and identification of host adaptation-associated genes of Ligilactobacillus salivarius isolated from badgers

BACKGROUND

Ligilactobacillus salivarius has been frequently isolated from the gut microbiota of humans and domesticated animals and has been studied as a candidate probiotic. Badger (Meles meles) is known as a "generalist" species that consumes complex foods and exhibits tolerance and resistance to certain pathogens, which can be partly attributed to the beneficial microbes such as L. salivarius in the gut microbiota. However, our understanding of the beneficial traits and genomic features of badger-originated L. salivarius remains elusive.

RESULTS

In this study, nine L. salivarius strains were isolated from wild badgers' feces, one of which exhibited good probiotic properties. Complete genomes of the nine L. salivarius strains were generated, and comparative genomic analysis was performed with the publicly available complete genomes of L. salivarius obtained from humans and domesticated animals. The strains originating from badgers harbored a larger genome, a higher number of protein-coding sequences, and functionally annotated genes than those originating from humans and chickens. The pan-genome phylogenetic tree demonstrated that the strains originating from badgers formed a separate clade, and totally 412 gene families (12.6% of the total gene families in the pan-genome) were identified as genes gained by the last common ancestor of the badger group. The badger group harbored significantly more gene families responsible for the degradation of complex carbohydrate substrates and production of polysaccharides than strains from other hosts; many of these were acquired by gene gain events.

CONCLUSIONS

A candidate probiotic and nine L. salivarius complete genomes were obtained from the badgers' gut microbiome, and several beneficial genes were identified to be specifically present in the badger-originated strains that were gained in the evolution. Our study provides novel insights into the adaptation of L. salivarius to the intestinal habitat of wild badgers and provides valuable strain and genome resources for the development of L. salivarius as a probiotic.

Bioinformatics and its role in the study of the evolution and probiotic potential of lactic acid bacteria

Due to their numerous well-established applications in the food industry, there have been many studies regarding the adaptation and evolution of lactic acid bacteria (LAB) in a wide variety of hosts and environments. Progress in sequencing technology and continual decreases in its costs have led to the availability of LAB genome sequence data. Bioinformatics has been central to the extraction of valuable information from these raw genome sequence data. This paper presents the roles of bioinformatics tools and databases in understanding the adaptation and evolution of LAB, as well as the bioinformatics methods used in the initial screening of LAB for probiotic potential. Moreover, the advantages, challenges, and limitations of employing bioinformatics for these purposes are discussed.

Comparative Genomics Analysis Provides New Insights into High Ethanol Tolerance of Lactiplantibacillus pentosus LTJ12, a Novel Strain Isolated from Chinese Baijiu

Lactic acid bacteria have received a significant amount of attention due to their probiotic characteristics. The species Lactiplantibacillus plantarum and Lactiplantibacillus pentosus are genotypically closely related, and their phenotypes are so similar that they are easily confused and mistaken. In the previous study, an ethanol-resistant strain, LTJ12, isolated from the fermented grains of soy sauce aroma type baijiu in North China, was originally identified as L. plantarum through a 16S rRNA sequence analysis. Here, the genome of strain LTJ12 was further sequenced using PacBio and Illumina sequencing technology to obtain a better understanding of the metabolic pathway underlying its resistance to ethanol stress. The results showed that the genome of strain LTJ12 was composed of one circular chromosome and three circular plasmids. The genome size is 3,512,307 bp with a GC content of 46.37%, and the number of predicted coding genes is 3248. Moreover, by comparing the coding genes with the GO (Gene Ontology), COG (Cluster of Orthologous Groups) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases, the functional annotation of the genome and an assessment of the metabolic pathways were performed, with the results showing that strain LTJ12 has multiple genes that may be related to alcohol metabolism and probiotic-related genes. Antibiotic resistance gene analysis showed that there were few potential safety hazards. Further, after conducting the comparative genomics analysis, it was found that strain LTJ12 is L. pentosus but not L. plantarum, but it has more functional genes than other L. pentosus strains that are mainly related to carbohydrate transport and metabolism, transcription, replication, recombination and repair, signal transduction mechanisms, defense mechanisms and cell wall/membrane/envelope biogenesis. These unique functional genes, such as gene 2754 (encodes alcohol dehydrogenase), gene 3093 (encodes gamma-D-glutamyl-meso-diaminopimelate peptidase) and some others may enhance the ethanol tolerance and alcohol metabolism of the strain. Taken together, L. pentosus LTJ12 might be a potentially safe probiotic with a high ethanol tolerance and alcohol metabolism. The findings of this study will also shed light on the accurate identification and rational application of the Lactiplantibacillus species.

Insights on the bacterial composition of Parmigiano Reggiano Natural Whey Starter by a culture-dependent and 16S rRNA metabarcoding portrait

Natural whey starters (NWS) are undefined bacterial communities produced daily from whey of the previous cheese-making round, by application of high temperature. As a result, in any dairy plant, NWS are continuously evolving, undefined mixtures of several strains and/or species of lactic acid bacteria, whose composition and performance strongly depend on the selective pressure acting during incubation. While NWS is critical to assure consistency to cheese-making process, little is known about the composition, functional features, and plant-to-plant fluctuations. Here, we integrated 16S rRNA metabarcoding and culture-dependent methods to profile bacterial communities of 10 NWS sampled in the production area of Parmigiano Reggiano cheese. 16S rRNA metabarcoding analysis revealed two main NWS community types, namely NWS type-H and NWS type-D. Lactobacillus helveticus was more abundant in NWS type-H, whilst Lactobacillus delbrueckii/St. thermophilus in NWS type-D, respectively. Based on the prediction of metagenome functions, NWS type-H samples were enriched in functional pathways related to galactose catabolism and purine metabolism, while NWS type-D in pathways related to aromatic and branched chain amino acid biosynthesis, which are flavor compound precursors. Culture-dependent approaches revealed low cultivability of individual colonies as axenic cultures and high genetic diversity in the pool of cultivable survivors. Co-culturing experiments showed that fermentative performance decreases by reducing the bacterial complexity of inoculum, suggesting that biotic interactions and cross-feeding relationships could take place in NWS communities, assuring phenotypic robustness. Even though our data cannot directly predict these ecological interactions, this study provides the basis for experiments targeted at understanding how selective regime affects composition, bacterial interaction, and fermentative performance in NWS.

The association of cervicovaginal Langerhans cells with clearance of human papillomavirus

Human papillomavirus (HPV) clearance is important in eliminating cervical cancer which contributes to high morbidity and mortality in women. Nevertheless, it remains largely unknown about key players in clearing pre-existing HPV infections. HPV antigens can be detected by the most important cervical antigen-presenting cells (Langerhans cells, LCs), of which the activities can be affected by cervicovaginal microbiota. In this review, we first introduce persistent HPV infections and then describe HPV-suppressed LCs activities, including but not limited to antigen uptake and presentation. Given specific transcriptional profiling of LCs in cervical epithelium, we also discuss the impact of cervicovaginal microbiota on LCs activation as well as the promise of exploring key microbial players in activating LCs and HPV-specific cellular immunity.

Functional characterization of prokaryotic dark matter: the road so far and what lies ahead

Eight-hundred thousand to one trillion prokaryotic species may inhabit our planet. Yet, fewer than two-hundred thousand prokaryotic species have been described. This uncharted fraction of microbial diversity, and its undisclosed coding potential, is known as the "microbial dark matter" (MDM). Next-generation sequencing has allowed to collect a massive amount of genome sequence data, leading to unprecedented advances in the field of genomics. Still, harnessing new functional information from the genomes of uncultured prokaryotes is often limited by standard classification methods. These methods often rely on sequence similarity searches against reference genomes from cultured species. This hinders the discovery of unique genetic elements that are missing from the cultivated realm. It also contributes to the accumulation of prokaryotic gene products of unknown function among public sequence data repositories, highlighting the need for new approaches for sequencing data analysis and classification. Increasing evidence indicates that these proteins of unknown function might be a treasure trove of biotechnological potential. Here, we outline the challenges, opportunities, and the potential hidden within the functional dark matter (FDM) of prokaryotes. We also discuss the pitfalls surrounding molecular and computational approaches currently used to probe these uncharted waters, and discuss future opportunities for research and applications.

Cervicovaginal microbiota significantly changed for HPV-positive women with high-grade squamous intraepithelial lesion

Lower female genital tract is colonized by a variety of microbes (cervicovaginal microbiota, CVM) which associate with the risk of genital infection. This study characterized CVM for 149 Chinese women with different status of human papillomavirus (HPV) infection and squamous intraepithelial lesion (SIL): no HPV infection (HPV-), HPV infection without significant SIL (HPV+NoSIL), HPV infection with low-grade SIL (HPV+LSIL) and HPV infection with high-grade SIL (HPV+HSIL). Analysis results showed CVM has dramatically changed in HPV+HSIL group when compared to HPV+LSIL group, but it exhibited no significant differences between HPV- and HPV+NoSIL groups as well as between HPV+NoSIL and HPV+LSIL groups. In consistence, random forest analysis found more notable differences in HPV+HSIL vs HPV+LSIL comparison than in other comparisons. In addition, depletion of Lactobacillus in CVM was more to be frequently identified in SIL-positive women as compared to SIL-negative individuals. Our findings suggested that significant CVM differences occurred when SIL developed to HSIL which was caused by persistent HPV infection.

Metagenomics and artificial intelligence in the context of human health

Human microbiome is ubiquitous, dynamic, and site-specific consortia of microbial communities. The pathogenic nature of microorganisms within human tissues has led to an increase in microbial studies. Characterization of genera, like Streptococcus, Cutibacterium, Staphylococcus, Bifidobacterium, Lactococcus and Lactobacillus through culture-dependent and culture-independent techniques has been reported. However, due to the unique environment within human tissues, it is difficult to culture these microorganisms making their molecular studies strenuous. MGs offer a gateway to explore and characterize hidden microbial communities through a culture-independent mode by direct DNA isolation. By function and sequence-based MGs, Scientists can explore the mechanistic details of numerous microbes and their interaction with the niche. Since the data generated from MGs studies is highly complex and multi-dimensional, it requires accurate analytical tools to evaluate and interpret the data. Artificial intelligence (AI) provides the luxury to automatically learn the data dimensionality and ease its complexity that makes the disease diagnosis and disease response easy, accurate and timely. This review provides insight into the human microbiota and its exploration and expansion through MG studies. The review elucidates the significance of MGs in studying the changing microbiota during disease conditions besides highlighting the role of AI in computational analysis of MG data.

Functional strain redundancy and persistent phage infection in Swiss hard cheese starter cultures

Undefined starter cultures are poorly characterized bacterial communities from environmental origin used in cheese making. They are phenotypically stable and have evolved through domestication by repeated propagation in closed and highly controlled environments over centuries. This makes them interesting for understanding eco-evolutionary dynamics governing microbial communities. While cheese starter cultures are known to be dominated by a few bacterial species, little is known about the composition, functional relevance, and temporal dynamics of strain-level diversity. Here, we applied shotgun metagenomics to an important Swiss cheese starter culture and analyzed historical and experimental samples reflecting 82 years of starter culture propagation. We found that the bacterial community is highly stable and dominated by only a few coexisting strains of Streptococcus thermophilus and Lactobacillus delbrueckii subsp. lactis. Genome sequencing, metabolomics analysis, and co-culturing experiments of 43 isolates show that these strains are functionally redundant, but differ tremendously in their phage resistance potential. Moreover, we identified two highly abundant Streptococcus phages that seem to stably coexist in the community without any negative impact on bacterial growth or strain persistence, and despite the presence of a large and diverse repertoire of matching CRISPR spacers. Our findings show that functionally equivalent strains can coexist in domesticated microbial communities and highlight an important role of bacteria-phage interactions that are different from kill-the-winner dynamics.

Functional strain redundancy and persistent phage infection in Swiss hard cheese starter cultures

Undefined starter cultures are poorly characterized bacterial communities from environmental origin used in cheese making. They are phenotypically stable and have evolved through domestication by repeated propagation in closed and highly controlled environments over centuries. This makes them interesting for understanding eco-evolutionary dynamics governing microbial communities. While cheese starter cultures are known to be dominated by a few bacterial species, little is known about the composition, functional relevance, and temporal dynamics of strain-level diversity. Here, we applied shotgun metagenomics to an important Swiss cheese starter culture and analyzed historical and experimental samples reflecting 82 years of starter culture propagation. We found that the bacterial community is highly stable and dominated by only a few coexisting strains of Streptococcus thermophilus and Lactobacillus delbrueckii subsp. lactis. Genome sequencing, metabolomics analysis, and co-culturing experiments of 43 isolates show that these strains are functionally redundant, but differ tremendously in their phage resistance potential. Moreover, we identified two highly abundant Streptococcus phages that seem to stably coexist in the community without any negative impact on bacterial growth or strain persistence, and despite the presence of a large and diverse repertoire of matching CRISPR spacers. Our findings show that functionally equivalent strains can coexist in domesticated microbial communities and highlight an important role of bacteria-phage interactions that are different from kill-the-winner dynamics.

Reconstitution and Transmission of Gut Microbiomes and Their Genes between Generations

Microbiomes are transmitted between generations by a variety of different vertical and/or horizontal modes, including vegetative reproduction (vertical), via female germ cells (vertical), coprophagy and regurgitation (vertical and horizontal), physical contact starting at birth (vertical and horizontal), breast-feeding (vertical), and via the environment (horizontal). Analyses of vertical transmission can result in false negatives (failure to detect rare microbes) and false positives (strain variants). In humans, offspring receive most of their initial gut microbiota vertically from mothers during birth, via breast-feeding and close contact. Horizontal transmission is common in marine organisms and involves selectivity in determining which environmental microbes can colonize the organism's microbiome. The following arguments are put forth concerning accurate microbial transmission: First, the transmission may be of functions, not necessarily of species; second, horizontal transmission may be as accurate as vertical transmission; third, detection techniques may fail to detect rare microbes; lastly, microbiomes develop and reach maturity with their hosts. In spite of the great variation in means of transmission discussed in this paper, microbiomes and their functions are transferred from one generation of holobionts to the next with fidelity. This provides a strong basis for each holobiont to be considered a unique biological entity and a level of selection in evolution, largely maintaining the uniqueness of the entity and conserving the species from one generation to the next.

Antimicrobial Susceptibility of Lactobacillus delbrueckii subsp. lactis from Milk Products and Other Habitats

As components of many cheese starter cultures, strains of Lactobacillus delbrueckii subsp. lactis (LDL) must be tested for their antimicrobial susceptibility to avoid the potential horizontal transfer of antibiotic resistance (ABR) determinants in the human body or in the environment. To this end, a phenotypic test, as well as a screening for antibiotic resistance genes (ARGs) in genome sequences, is commonly performed. Historically, microbiological cutoffs (MCs), which are used to classify strains as either 'sensitive' or 'resistant' based on the minimal inhibitory concentrations (MICs) of a range of clinically-relevant antibiotics, have been defined for the whole group of the obligate homofermentative lactobacilli, which includes LDL among many other species. This often leads to inaccuracies in the appreciation of the ABR status of tested LDL strains and to false positive results. To define more accurate MCs for LDL, we analyzed the MIC profiles of strains originating from various habitats by using the broth microdilution method. These strains' genomes were sequenced and used to complement our analysis involving a search for ARGs, as well as to assess the phylogenetic proximity between strains. Of LDL strains, 52.1% displayed MICs that were higher than the defined MCs for kanamycin, 9.9% for chloramphenicol, and 5.6% for tetracycline, but no ARG was conclusively detected. On the other hand, all strains displayed MICs below the defined MCs for ampicillin, gentamycin, erythromycin, and clindamycin. Considering our results, we propose the adaptation of the MCs for six of the tested clinically-relevant antibiotics to improve the accuracy of phenotypic antibiotic testing.

Massive Survey on Bacterial-Bacteriophages Biodiversity and Quality of Natural Whey Starter Cultures in Trentingrana Cheese Production

This study focused on the microbial and bacteriophages identification and characterization in cheese-production facilities that use natural whey starter (NWS) cultures for Trentingrana production. Bacterial and phage screening was carried out on cooked not acidified whey and NWS samples isolated from six dairy factories, for 4 consecutive days in four different months. By means of a combined approach, using plate counts, bacterial isolation, and metataxonomic analysis Lactobacillus helveticus was found occurring as the dominant species in NWS cultures and Levilactobacillus brevis as codominant in the cheese factories where the temperature of NWS production was mainly lower than 40°C, suggesting that the variability in the parameters of the NWS culture preparation could differently modulate the bacterial species in NWS cultures. Using turbidity test approach on 303 bacterial isolates from the NWS cultures, 120 distinct phages were identified. L. helveticus phage contamination of NWS cultures was revealed in most of the analyzed samples, but despite the great recovery of bacteriophage contamination cases, the microbial quality of NWS cultures was high. Our results support the presence of natural bacteriophage resistance mechanisms in L. helveticus. The use of NWS cultures probably creates an ideal environment for the proliferation of different L. helveticus strains balanced with their phages without a clear dominance. It is evident, from this study, that the presence of a high biodiversity of NWS bacterial strains is relevant to avoid phages dominance in NWS cultures and consequently to keep a good acidification ability.

Genome and Pangenome Analysis of Lactobacillus hilgardii FLUB-A New Strain Isolated from Mead

The production of mead holds great value for the Polish liquor industry, which is why the bacterium that spoils mead has become an object of concern and scientific interest. This article describes, for the first time, Lactobacillus hilgardii FLUB newly isolated from mead, as a mead spoilage bacteria. Whole genome sequencing of L. hilgardii FLUB revealed a 3 Mbp chromosome and five plasmids, which is the largest reported genome of this species. An extensive phylogenetic analysis and digital DNA-DNA hybridization confirmed the membership of the strain in the L. hilgardii species. The genome of L. hilgardii FLUB encodes 3043 genes, 2871 of which are protein coding sequences, 79 code for RNA, and 93 are pseudogenes. L. hilgardii FLUB possesses three clustered regularly interspaced short palindromic repeats (CRISPR), eight genomic islands (44,155 bp to 6345 bp), and three (two intact and one incomplete) prophage regions. For the first time, the characteristics of the genome of this species were described and a pangenomic analysis was performed. The concept of the pangenome was used not only to establish the genetic repertoire of this species, but primarily to highlight the unique characteristics of L. hilgardii FLUB. The core of the genome of L. hilgardii is centered around genes related to the storage and processing of genetic information, as well as to carbohydrate and amino acid metabolism. Strains with such a genetic constitution can effectively adapt to environmental changes. L. hilgardii FLUB is distinguished by an extensive cluster of metabolic genes, arsenic detoxification genes, and unique surface layer proteins. Variants of MRS broth with ethanol (10-20%), glucose (2-25%), and fructose (2-24%) were prepared to test the strain's growth preferences using Bioscreen C and the PYTHON script. L. hilgardii FLUB was found to be more resistant than a reference strain to high concentrations of alcohol (18%) and sugars (25%). It exhibited greater preference for fructose than glucose, which suggests it has a fructophilic nature. Comparative genomic analysis supported by experimental research imitating the conditions of alcoholic beverages confirmed the niche specialization of L. hilgardii FLUB to the mead environment.

Niche-specific adaptation of Lactobacillus helveticus strains isolated from malt whisky and dairy fermentations

Lactobacillus helveticus is a well characterized lactobacillus for dairy fermentations that is also found in malt whisky fermentations. The two environments contain considerable differences related to microbial growth, including the presence of different growth inhibitors and nutrients. The present study characterized L. helveticus strains originating from dairy fermentations (called milk strains hereafter) and malt whisky fermentations (called whisky strains hereafter) by in vitro phenotypic tests and comparative genomics. The whisky strains can tolerate ethanol more than the milk strains, whereas the milk strains can tolerate lysozyme and lactoferrin more than the whisky strains. Several plant-origin carbohydrates, including cellobiose, maltose, sucrose, fructooligosaccharide and salicin, were generally metabolized only by the whisky strains, whereas milk-derived carbohydrates, i.e. lactose and galactose, were metabolized only by the milk strains. Milk fermentation properties also distinguished the two groups. The general genomic characteristics, including genomic size, number of coding sequences and average nucleotide identity values, differentiated the two groups. The observed differences in carbohydrate metabolic properties between the two groups correlated with the presence of intact specific enzymes in glycoside hydrolase (GH) families GH1, GH4, GH13, GH32 and GH65. Several GHs in the milk strains were inactive due to the presence of stop codon(s) in genes encoding the GHs, and the inactivation patterns of the genes encoding specific enzymes assigned to GH1 in the milk strains suggested a possible diversification manner of L. helveticus strains. The present study has demonstrated how L. helveticus strains have adapted to their habitats.

An integrated model system to gain mechanistic insights into biofilm-associated antimicrobial resistance in Pseudomonas aeruginosa MPAO1

Pseudomonas aeruginosa MPAO1 is the parental strain of the widely utilized transposon mutant collection for this important clinical pathogen. Here, we validate a model system to identify genes involved in biofilm growth and biofilm-associated antibiotic resistance. Our model employs a genomics-driven workflow to assemble the complete MPAO1 genome, identify unique and conserved genes by comparative genomics with the PAO1 reference strain and genes missed within existing assemblies by proteogenomics. Among over 200 unique MPAO1 genes, we identified six general essential genes that were overlooked when mapping public Tn-seq data sets against PAO1, including an antitoxin. Genomic data were integrated with phenotypic data from an experimental workflow using a user-friendly, soft lithography-based microfluidic flow chamber for biofilm growth and a screen with the Tn-mutant library in microtiter plates. The screen identified hitherto unknown genes involved in biofilm growth and antibiotic resistance. Experiments conducted with the flow chamber across three laboratories delivered reproducible data on P. aeruginosa biofilms and validated the function of both known genes and genes identified in the Tn-mutant screens. Differential protein abundance data from planktonic cells versus biofilm confirmed the upregulation of candidates known to affect biofilm formation, of structural and secreted proteins of type VI secretion systems, and provided proteogenomic evidence for some missed MPAO1 genes. This integrated, broadly applicable model promises to improve the mechanistic understanding of biofilm formation, antimicrobial tolerance, and resistance evolution in biofilms.

Acetate kinase and peptidases are associated with the proteolytic activity of Lactobacillus helveticus isolated from fermented food

Lactobacillus (L.) helveticus is widely used in food industry due to its high proteolytic activity. However, such activity varies greatly between isolates, and the determining factors regulating the strength of proteolytic activity in L. helveticus are unclear. This study sequenced the genomes of 60 fermented food-originated L. helveticus and systemically examined the proteolytic activity-determining factors. Our analyses found that the strength of proteolytic activity in L. helveticus was independent of the isolation source, geographic location, phylogenetic closeness between isolates, and distribution of cell envelope proteinases (CEPs). Genome-wide association study (GWAS) identified two genes, the acetate kinase (ackA) and a hypothetical protein, and 15 single nucleotide polymorphisms (SNPs) that were associated with the strength of the proteolytic activity. Further investigating the functions of these gene components revealed that ackA and two cysteine peptidases coding genes (pepC and srtA) rather than the highly heterogeneous and intraspecific CEPs were linked to the level of proteolytic activity. Moreover, the sequence type (ST) defined by SNP analysis revealed a total of ten STs, and significantly weaker proteolytic activity was observed among isolates of ST2. This study provides practical information for future selection of L. helveticus of strong proteolytic activity.

Core and Accessory Genome Comparison of Australian and International Strains of O157 Shiga Toxin-Producing Escherichia coli

Shiga toxin-producing Escherichia coli (STEC) is a foodborne pathogen, and serotype O157:H7 is typically associated with severe disease. Australian STEC epidemiology differs from many other countries, as severe outbreaks and HUS cases appear to be more often associated with non-O157 serogroups. It is not known why Australian strains of O157 STEC might differ in virulence to international strains. Here we investigate the reduced virulence of Australian strains. Multiple genetic analyses were performed, including SNP-typing, to compare the core genomes of the Australian to the international isolates, and accessory genome analysis to determine any significant differences in gene presence/absence that could be associated with their phenotypic differences in virulence. The most distinct difference between the isolates was the absence of the stx2a gene in all Australian isolates, with few other notable differences observed in the core and accessory genomes of the O157 STEC isolates analyzed in this study. The presence of stx1a in most Australian isolates was another notable observation. Acquisition of stx2a seems to coincide with the emergence of highly pathogenic STEC. Due to the lack of other notable genotypic differences observed between Australian and international isolates characterized as highly pathogenic, this may be further evidence that the absence of stx2a in Australian O157 STEC could be a significant characteristic defining its mild virulence. Further work investigating the driving force(s) behind Stx prophage loss and acquisition is needed to determine if this potential exists in Australian O157 isolates.

Linking Comparative Genomics of Nine Potato-Associated Pseudomonas Isolates With Their Differing Biocontrol Potential Against Late Blight

For plants, the advantages of associating with beneficial bacteria include plant growth promotion, reduction of abiotic and biotic stresses and enhanced protection against various pests and diseases. Beneficial bacteria rightly equipped for successful plant colonization and showing antagonistic activity toward plant pathogens seem to be actively recruited by plants. To gain more insights into the genetic determinants responsible for plant colonization and antagonistic activities, we first sequenced and de novo assembled the complete genomes of nine Pseudomonas strains that had exhibited varying antagonistic potential against the notorious oomycete Phytophthora infestans, placed them into the phylogenomic context of known Pseudomonas biocontrol strains and carried out a comparative genomic analysis to define core, accessory (i.e., genes found in two or more, but not all strains) and unique genes. Next, we assessed the colonizing abilities of these strains and used bioassays to characterize their inhibitory effects against different stages of P. infestans' lifecycle. The phenotype data were then correlated with genotype information, assessing over three hundred genes encoding known factors for plant colonization and antimicrobial activity as well as secondary metabolite biosynthesis clusters predicted by antiSMASH. All strains harbored genes required for successful plant colonization but also distinct arsenals of antimicrobial compounds. We identified genes coding for phenazine, hydrogen cyanide, 2-hexyl, 5-propyl resorcinol and pyrrolnitrin synthesis, as well as various siderophores, pyocins and type VI secretion systems. Additionally, the comparative genomic analysis revealed about a hundred accessory genes putatively involved in anti-Phytophthora activity, including a type II secretion system (T2SS), several peptidases and a toxin. Transcriptomic studies and mutagenesis are needed to further investigate the putative involvement of the novel candidate genes and to identify the various mechanisms involved in the inhibition of P. infestans by different Pseudomonas strains.

A Proteogenomic Resource Enabling Integrated Analysis of Listeria Genotype-Proteotype-Phenotype Relationships

Listeria monocytogenes is an opportunistic foodborne pathogen responsible for listeriosis, a potentially fatal foodborne disease. Many different Listeria strains and serotypes exist, but a proteogenomic resource that bridges the gap in our molecular understanding of the relationships between the Listeria genotypes and phenotypes via proteotypes is still missing. Here, we devised a next-generation proteogenomics strategy that enables the community to rapidly proteotype Listeria strains and relate this information back to the genotype. Based on sequencing and de novo assembly of the two most commonly used Listeria model strains, EGD-e and ScottA, we established two comprehensive Listeria proteogenomic databases. A genome comparison established core- and strain-specific genes potentially responsible for virulence differences. Next, we established a DIA/SWATH-based proteotyping strategy, including a new and robust sample preparation workflow, that enables the reproducible, sensitive, and relative quantitative measurement of Listeria proteotypes. This reusable and publicly available DIA/SWATH library covers 70% of open reading frames of Listeria and represents the most extensive spectral library for Listeria proteotype analysis to date. We used these two new resources to investigate the Listeria proteotype in states mimicking the upper gastrointestinal passage. Exposure of Listeria to bile salts at 37 °C, which simulates conditions encountered in the duodenum, showed significant proteotype perturbations including an increase of FlaA, the structural protein of flagella. Given that Listeria is known to lose its flagella above 30 °C, this was an unexpected finding. The formation of flagella, which might have implications on infectivity, was validated by parallel reaction monitoring and light and scanning electron microscopy. flaA transcript levels did not change significantly upon exposure to bile salts at 37 °C, suggesting regulation at the post-transcriptional level. Together, these analyses provide a comprehensive proteogenomic resource and toolbox for the Listeria community enabling the analysis of Listeria genotype-proteotype-phenotype relationships.

Oligotrophic bacterium Hymenobacter latericoloratus CGMCC 16346 degrades the neonicotinoid imidacloprid in surface water

The intensive and extensive application of imidacloprid in agriculture has resulted in water pollution and risks to aquatic invertebrates. However, pure bacteria remediation of imidacloprid in surface water environments has not been studied. Here, we isolated an imidacloprid-degrading bacterium from a water environment, examined its imidacloprid degradation in pure culture and surface water, sequenced its genome, and compared its Clusters of Orthologous Groups (COG) protein categorization with that for another imidacloprid-degrading bacterium. The isolate was an obligate oligotrophic bacterium, Hymenobacter latericoloratus CGMCC 16346, which degraded imidacloprid via hydroxylation by co-metabolism in pure culture. Resting cells degraded 64.4% of 100 mg/L imidacloprid in 6 days in the presence of co-substrate maltose, and growing culture degraded 40.8% of imidacloprid in 10 days. H. latericoloratus CGMCC 16346 degraded imidacloprid in surface water without co-substrate supplementation and retained imidacloprid-degrading activity after 30 days. The half-life of imidacloprid in surface water was decreased from 173.3 days in the control to 57.8 days by CGMCC 16346 inoculation. Genome sequencing and COG analysis indicated that carbohydrate metabolism and transport, cell wall/membrane biogenesis, and defense mechanisms are enriched in H. latericoloratus CGMCC 16346 compared with the copiotrophic imidacloprid-degrading Pseudoxanthomonas indica CGMCC 6648, indicating that H. latericoloratus CGMCC 16346 is adapted to live in oligotrophic water environments and biofilms. H. latericoloratus CGMCC 16346 is a promising bioremediation agent for elimination of imidacloprid contamination from surface water.

Draft Genome Sequences of Six Lactobacillus gasseri and Three Lactobacillus paragasseri Strains Isolated from the Female Bladder

Lactobacilli are dominant members of the healthy female bladder microbiota. Here, we report the complete genome sequences of six Lactobacillus gasseri and three Lactobacillus paragasseri strains isolated from catheterized urine samples. These L. paragasseri genomes are the first publicly available sequences of the species from the bladder.

Lactobacilli are dominant members of the healthy female bladder microbiota. Here, we report the complete genome sequences of six Lactobacillus gasseri and three Lactobacillus paragasseri strains isolated from catheterized urine samples. These L. paragasseri genomes are the first publicly available sequences of the species from the bladder.

Draft Genome Sequences of 11 Lactobacillus jensenii Strains Isolated from the Female Bladder

Lactobacillus jensenii, a protective bacterium in the vaginal microbiota, is also a member of the female urinary tract community. Here, we report 11 genome sequences of L. jensenii strains isolated from catheterized urine from women. This effort greatly increases our knowledge of the genetic diversity of this species within the bladder.

Comparative Genomics of Marine Sponge-Derived Streptomyces spp. Isolates SM17 and SM18 With Their Closest Terrestrial Relatives Provides Novel Insights Into Environmental Niche Adaptations and Secondary Metabolite Biosynthesis Potential

The emergence of antibiotic resistant microorganisms has led to an increased need for the discovery and development of novel antimicrobial compounds. Frequent rediscovery of the same natural products (NPs) continues to decrease the likelihood of the discovery of new compounds from soil bacteria. Thus, efforts have shifted toward investigating microorganisms and their secondary metabolite biosynthesis potential, from diverse niche environments, such as those isolated from marine sponges. Here we investigated at the genomic level two Streptomyces spp. strains, namely SM17 and SM18, isolated from the marine sponge Haliclona simulans, with previously reported antimicrobial activity against clinically relevant pathogens; using single molecule real-time (SMRT) sequencing. We performed a series of comparative genomic analyses on SM17 and SM18 with their closest terrestrial relatives, namely S. albus J1074 and S. pratensis ATCC 33331 respectively; in an effort to provide further insights into potential environmental niche adaptations (ENAs) of marine sponge-associated Streptomyces, and on how these adaptations might be linked to their secondary metabolite biosynthesis potential. Prediction of secondary metabolite biosynthetic gene clusters (smBGCs) indicated that, even though the marine isolates are closely related to their terrestrial counterparts at a genomic level; they potentially produce different compounds. SM17 and SM18 displayed a better ability to grow in high salinity medium when compared to their terrestrial counterparts, and further analysis of their genomes indicated that they possess a pool of 29 potential ENA genes that are absent in S. albus J1074 and S. pratensis ATCC 33331. This ENA gene pool included functional categories of genes that are likely to be related to niche adaptations and which could be grouped based on potential biological functions such as osmotic stress, defense; transcriptional regulation; symbiotic interactions; antimicrobial compound production and resistance; ABC transporters; together with horizontal gene transfer and defense-related features.

Genomic Comparison of Lactobacillus helveticus Strains Highlights Probiotic Potential

Lactobacillus helveticus belongs to the large group of lactic acid bacteria (LAB), which are the major players in the fermentation of a wide range of foods. LAB are also present in the human gut, which has often been exploited as a reservoir of potential novel probiotic strains, but several parameters need to be assessed before establishing their safety and potential use for human consumption. In the present study, six L. helveticus strains isolated from natural whey cultures were analyzed for their phenotype and genotype in exopolysaccharide (EPS) production, low pH and bile salt tolerance, bile salt hydrolase (BSH) activity, and antibiotic resistance profile. In addition, a comparative genomic investigation was performed between the six newly sequenced strains and the 51 publicly available genomes of L. helveticus to define the pangenome structure. The results indicate that the newly sequenced strain UC1267 and the deposited strain DSM 20075 can be considered good candidates for gut-adapted strains due to their ability to survive in the presence of 0.2% glycocholic acid (GCA) and 1% taurocholic and taurodeoxycholic acid (TDCA). Moreover, these strains had the highest bile salt deconjugation activity among the tested L. helveticus strains. Considering the safety profile, none of these strains presented antibiotic resistance phenotypically and/or at the genome level. The pangenome analysis revealed genes specific to the new isolates, such as enzymes related to folate biosynthesis in strains UC1266 and UC1267 and an integrated phage in strain UC1035. Finally, the presence of maltose-degrading enzymes and multiple copies of 6-phospho-β-glucosidase genes in our strains indicates the capability to metabolize sugars other than lactose, which is related solely to dairy niches.

Pushing the limits of de novo genome assembly for complex prokaryotic genomes harboring very long, near identical repeats

Generating a complete, de novo genome assembly for prokaryotes is often considered a solved problem. However, we here show that Pseudomonas koreensis P19E3 harbors multiple, near identical repeat pairs up to 70 kilobase pairs in length, which contained several genes that may confer fitness advantages to the strain. Its complex genome, which also included a variable shufflon region, could not be de novo assembled with long reads produced by Pacific Biosciences’ technology, but required very long reads from Oxford Nanopore Technologies. Importantly, a repeat analysis, whose results we release for over 9600 prokaryotes, indicated that very complex bacterial genomes represent a general phenomenon beyond Pseudomonas. Roughly 10% of 9331 complete bacterial and a handful of 293 complete archaeal genomes represented this ‘dark matter’ for de novo genome assembly of prokaryotes. Several of these ‘dark matter’ genome assemblies contained repeats far beyond the resolution of the sequencing technology employed and likely contain errors, other genomes were closed employing labor-intense steps like cosmid libraries, primer walking or optical mapping. Using very long sequencing reads in combination with assembly algorithms capable of resolving long, near identical repeats will bring most prokaryotic genomes within reach of fast and complete de novo genome assembly.

Long-read based de novo assembly of low-complexity metagenome samples results in finished genomes and reveals insights into strain diversity and an active phage system

BACKGROUND

Complete and contiguous genome assemblies greatly improve the quality of subsequent systems-wide functional profiling studies and the ability to gain novel biological insights. While a de novo genome assembly of an isolated bacterial strain is in most cases straightforward, more informative data about co-existing bacteria as well as synergistic and antagonistic effects can be obtained from a direct analysis of microbial communities. However, the complexity of metagenomic samples represents a major challenge. While third generation sequencing technologies have been suggested to enable finished metagenome-assembled genomes, to our knowledge, the complete genome assembly of all dominant strains in a microbiome sample has not been demonstrated. Natural whey starter cultures (NWCs) are used in cheese production and represent low-complexity microbiomes. Previous studies of Swiss Gruyère and selected Italian hard cheeses, mostly based on amplicon metagenomics, concurred that three species generally pre-dominate: Streptococcus thermophilus, Lactobacillus helveticus and Lactobacillus delbrueckii.

RESULTS

Two NWCs from Swiss Gruyère producers were subjected to whole metagenome shotgun sequencing using the Pacific Biosciences Sequel and Illumina MiSeq platforms. In addition, longer Oxford Nanopore Technologies MinION reads had to be generated for one to resolve repeat regions. Thereby, we achieved the complete assembly of all dominant bacterial genomes from these low-complexity NWCs, which was corroborated by a 16S rRNA amplicon survey. Moreover, two distinct L. helveticus strains were successfully co-assembled from the same sample. Besides bacterial chromosomes, we could also assemble several bacterial plasmids and phages and a corresponding prophage. Biologically relevant insights were uncovered by linking the plasmids and phages to their respective host genomes using DNA methylation motifs on the plasmids and by matching prokaryotic CRISPR spacers with the corresponding protospacers on the phages. These results could only be achieved by employing long-read sequencing data able to span intragenomic as well as intergenomic repeats.

CONCLUSIONS

Here, we demonstrate the feasibility of complete de novo genome assembly of all dominant strains from low-complexity NWCs based on whole metagenomics shotgun sequencing data. This allowed to gain novel biological insights and is a fundamental basis for subsequent systems-wide omics analyses, functional profiling and phenotype to genotype analysis of specific microbial communities.

Improved genome of Agrobacterium radiobacter type strain provides new taxonomic insight into Agrobacterium genomospecies 4

The reported Agrobacterium radiobacter DSM 30174^T genome is highly fragmented, hindering robust comparative genomics and genome-based taxonomic analysis. We re-sequenced the Agrobacterium radiobacter type strain, generating a dramatically improved genome with high contiguity. In addition, we sequenced the genome of Agrobacterium tumefaciens B6^T, enabling for the first time, a proper comparative genomics of these contentious Agrobacterium species. We provide concrete evidence that the previously reported Agrobacterium radiobacter type strain genome (Accession Number: ASXY01) is contaminated which explains its abnormally large genome size and fragmented assembly. We propose that Agrobacterium tumefaciens be reclassified as Agrobacterium radiobacter subsp. tumefaciens and that Agrobacterium radiobacter retains it species status with the proposed name of Agrobacterium radiobacter subsp. radiobacter. This proposal is based, first on the high pairwise genome-scale average nucleotide identity supporting the amalgamation of both Agrobacterium radiobacter and Agrobacterium tumefaciens into a single species. Second, maximum likelihood tree construction based on the concatenated alignment of shared genes (core genes) among related strains indicates that Agrobacterium radiobacter NCPPB3001 is sufficiently divergent from Agrobacterium tumefaciens to propose two independent sub-clades. Third, Agrobacterium tumefaciens demonstrates the genomic potential to synthesize the L configuration of fucose in its lipid polysaccharide, fostering its ability to colonize plant cells more effectively than Agrobacterium radiobacter.

Phylogenomic Analysis of Lactobacillus curvatus Reveals Two Lineages Distinguished by Genes for Fermenting Plant-Derived Carbohydrates

Lactobacillus curvatus is a lactic acid bacterium encountered in many different types of fermented food (meat, seafood, vegetables, and cereals). Although this species plays an important role in the preservation of these foods, few attempts have been made to assess its genomic diversity. This study uses comparative analyses of 13 published genomes (complete or draft) to better understand the evolutionary processes acting on the genome of this species. Phylogenomic analysis, based on a coalescent model of evolution, revealed that the 6,742 sites of single nucleotide polymorphism within the L. curvatus core genome delineate two major groups, with lineage 1 represented by the newly sequenced strain FLEC03, and lineage 2 represented by the type-strain DSM20019. The two lineages could also be distinguished by the content of their accessory genome, which sheds light on a long-term evolutionary process of lineage-dependent genetic acquisition and the possibility of population structure. Interestingly, one clade from lineage 2 shared more accessory genes with strains of lineage 1 than with other strains of lineage 2, indicating recent convergence in carbohydrate catabolism. Both lineages had a wide repertoire of accessory genes involved in the fermentation of plant-derived carbohydrates that are released from polymers of α/β-glucans, α/β-fructans, and N-acetylglucosan. Other gene clusters were distributed among strains according to the type of food from which the strains were isolated. These results give new insight into the ecological niches in which L. curvatus may naturally thrive (such as silage or compost heaps) in addition to fermented food.

Symposium review: Genomic investigations of flavor formation by dairy microbiota

Flavor is one of the most important attributes of any fermented dairy product. Dairy consumers are known to be willing to experiment with different flavors; thus, many companies producing fermented dairy products have looked at culture manipulation as a tool for flavor diversification. The development of flavor is a complex process, originating from a combination of microbiological, biochemical, and technological aspects. A key driver of flavor is the enzymatic activities of the deliberately inoculated starter cultures, in addition to the environmental or "nonstarter" microbiota. The contribution of microbial metabolism to flavor development in fermented dairy products has been exploited for thousands of years, but the availability of the whole genome sequences of the bacteria and yeasts involved in the fermentation process and the possibilities now offered by next-generation sequencing and downstream "omics" technologies is stimulating a more knowledge-based approach to the selection of desirable cultures for flavor development. By linking genomic traits to phenotypic outputs, it is now possible to mine the metabolic diversity of starter cultures, analyze the metabolic routes to flavor compound formation, identify those strains with flavor-forming potential, and select them for possible commercial application. This approach also allows for the identification of species and strains not previously considered as potential flavor-formers, the blending of strains with complementary metabolic pathways, and the potential improvement of key technological characteristics in existing strains, strains that are at the core of the dairy industry. An in-depth knowledge of the metabolic pathways of individual strains and their interactions in mixed culture fermentations can allow starter blends to be custom-made to suit industry needs. Applying this knowledge to starter culture research programs is enabling research and development scientists to develop superior starters, expand flavor profiles, and potentially develop new products for future market expansion.

Linezolid resistance genes and genetic elements enhancing their dissemination in enterococci and streptococci

Linezolid is considered a last resort drug in treatment of severe infections caused by Gram-positive pathogens, resistant to other antibiotics, such as vancomycin-resistant enterococci (VRE), methicillin-resistant staphylococci and multidrug resistant pneumococci. Although the vast majority of Gram-positive pathogenic bacteria remain susceptible to linezolid, resistant isolates of enterococci, staphylococci and streptococci have been reported worldwide. In these bacteria, apart from mutations, affecting mostly the 23S rRNA genes, acquisition of such genes as cfr, cfr(B), optrA and poxtA, often associated with mobile genetic elements (MGE), plays an important role for resistance. The purpose of this paper is to provide an overview on diversity and epidemiology of MGE carrying linezolid-resistance genes among clinically-relevant Gram-positive pathogens such as enterococci and streptococci.

Conserved Genome Organization and Core Transcriptome of the Lactobacillus acidophilus Complex

The Lactobacillus genus encompasses a genetically and functionally diverse group of species, and contains many strains widely formulated in the human food supply chain as probiotics and starter cultures. Within this genetically expansive group, there are several distinct clades that have high levels of homology, one of which is the Lactobacillus acidophilus group. Of the uniting features, small genomes, low GC content, adaptation to dairy environments, and fastidious growth requirements, are some of the most defining characteristics of this group. To better understand what truly links and defines this clade, we sought to characterize the genomic organization and content of the genomes of several members of this group. Through core genome analysis we explored the synteny and intrinsic genetic underpinnings of the L. acidophilus clade, and observed key features related to the evolution and adaptation of these organisms. While genetic content is able to provide a large map of the potential of each organism, it does not always reflect their functionality. Through transcriptomic data we inferred the core transcriptome of the L. acidophilus complex to better define the true metabolic capabilities that unite this clade. Using this approach we have identified seven small ORFs that are both highly conserved and transcribed in diverse members of this clade and could be potential novel small peptide or untranslated RNA regulators. Overall, our results reveal the core features of the L. acidophilus complex and open new avenues for the enhancement and formulation and of next generation probiotics and starter cultures.

Population Dynamics of Lactobacillus helveticus in Swiss Gruyère-Type Cheese Manufactured With Natural Whey Cultures

Lactobacillus helveticus, a ubiquitous bacterial species in natural whey cultures (NWCs) used for Swiss Gruyère cheese production, is considered to have crucial functions for cheese ripening such as enhancing proteolysis. We tracked the diversity and abundance of L. helveticus strains during 6 months of ripening in eight Swiss Gruyère-type cheeses using a culture-independent typing method. The study showed that the L. helveticus population present in NWCs persisted in cheese and demonstrated a stable multi-strain coexistence during cheese ripening. With regard to proteolysis, one of the eight L. helveticus populations exhibited less protein degradation during ripening.