Gut microbiota: methodological aspects to describe taxonomy and functionality.

The human gut is populated by a rich and diverse microbiota that has been described as a human metabolic organ. Its composition has an impact on the health status of the host and could be the target or effector of dietary health effects. The advent of high-throughput sequencing technologies and their integration with advanced analysis methods enables the development of new approaches to characterize the gut microbiota composition and initiates the understanding of its functionality. These technological aspects are not necessarily the limiting factor to successfully identify biological correlations between gut microbiota, diet and health status, as there are other very important aspects. The aim of this article is to address the gut microbiota characterization methodologies, highlight some advantages and challenges, and give our opinion on how critical the sampling and the study design are.

Genome sequence of the lantibiotic bacteriocin producer Streptococcus salivarius strain K12.

Streptococcus salivarius is a prevalent commensal species of the oropharyngeal tract. S. salivarius strain K12 is an isolate from the saliva of a healthy child, used as an oral probiotic. Here, we report its genome sequence, i.e., the full sequence of the 190-kb megaplasmid pSsal-K12 and a high-quality draft 2.2-Gb chromosomal sequence.

Campylobacter jejuni group III phage CP81 contains many T4-like genes without belonging to the T4-type phage group: implications for the evolution of T4 phages.

CP81 is a virulent Campylobacter group III phage whose linear genome comprises 132,454 bp. At the nucleotide level, CP81 differs from other phages. However, a number of its structural and replication/recombination proteins revealed a relationship to the group II Campylobacter phages CP220/CPt10 and to T4-type phages. Unlike the T4-related phages, the CP81 genome does not contain conserved replication and virion modules. Instead, the respective genes are scattered throughout the phage genome. Moreover, most genes for metabolic enzymes of CP220/CPt10 are lacking in CP81. On the other hand, the CP81 genome contains nine similar genes for homing endonucleases which may be involved in the attrition of the conserved gene order for the virion core genes of T4-type phages. The phage apparently possesses an unusual modification of C or G bases. Efficient cleavage of its DNA was only achieved with restriction enzymes recognizing pure A/T sites. Uncommonly, phenol extraction leads to a significant loss of CP81 DNA from the aqueous layer, a property not yet described for other phages belonging to the T4 superfamily.

Genetic changes are introduced by repeated exposure of Salmonella spiked in low water activity and high fat matrix to heat.

WGS is used to define if isolates are "in" or "out" of an outbreak and/or microbial root cause investigation. No threshold of genetic differences is fixed and the conclusions on similarity between isolates are mainly based on the knowledge generated from previous outbreak investigations and reported mutation rates. Mutation rates in Salmonella when exposed to food processing conditions are lacking. Thus, in this study, the ability of heat and dry stress to cause genetic changes in two Salmonella serotypes frequently isolated from low moisture foods was investigated. S. enterica serovars S. Agona ATCC 51,957 and S. Mbandaka NCTC 7892 (ATCC 51,958) were repeatedly exposed to heat (90 °C for 5 min) in a low water activity and high fat matrix. No increased fitness of the strains was observed after 10 repeated heat treatments. However, genetic changes were introduced and the number of genetic differences increased with every heat treatment cycle. The genetic changes appeared randomly in the genome and were responsible for a population of diverse isolates with 0 to 28 allelic differences (0 to 38 SNPs) between them. This knowledge is key to interpret WGS results for source tracking investigations as part of a root cause analysis in a contamination event as isolates are exposed to stress conditions.

Whole Genome Sequencing Applied to Pathogen Source Tracking in Food Industry: Key Considerations for Robust Bioinformatics Data Analysis and Reliable Results Interpretation.

Whole genome sequencing (WGS) has arisen as a powerful tool to perform pathogen source tracking in the food industry thanks to several developments in recent years. However, the cost associated to this technology and the degree of expertise required to accurately process and understand the data has limited its adoption at a wider scale. Additionally, the time needed to obtain actionable information is often seen as an impairment for the application and use of the information generated via WGS. Ongoing work towards standardization of wet lab including sequencing protocols, following guidelines from the regulatory authorities and international standardization efforts make the technology more and more accessible. However, data analysis and results interpretation guidelines are still subject to initiatives coming from distinct groups and institutions. There are multiple bioinformatics software and pipelines developed to handle such information. Nevertheless, little consensus exists on a standard way to process the data and interpret the results. Here, we want to present the constraints we face in an industrial setting and the steps we consider necessary to obtain high quality data, reproducible results and a robust interpretation of the obtained information. All of this, in a time frame allowing for data-driven actions supporting factories and their needs.

Whole genome sequencing used in an industrial context reveals a Salmonella laboratory cross-contamination.

In 2013, during a routine laboratory analysis performed on food samples, one finished product from a European factory was tested positive for Salmonella Hadar. At the same period, one environmental isolate in the same laboratory was serotyped Salmonella Hadar. Prior to this event, the laboratory performed a proficiency testing involving a sample spiked with NCTC 9877 Salmonella Hadar. The concomitance of Salmonella Hadar detection led to the suspicion of a laboratory cross-contamination between the Salmonella Hadar isolate used in the laboratory proficiency testing and the Salmonella Hadar isolate found on the finished product by the same laboratory. Since the classical phenotypic serotyping method is able to attribute a serotype to Salmonella isolates with a common antigenic formula, but cannot differentiate strains of the same serotype within the subspecies, whole genome sequencing was used to test the laboratory cross-contamination hypothesis. Additionally, 12 Salmonella Hadar from public databases, available until the time of the event, were included in the whole genome sequencing analysis to better understand the genomic diversity of this serotype in Europe. The outcome of the analysis showed a maximum of ten single nucleotide polymorphisms (SNPs) between the isolates coming from the laboratory and the finished product, and thus confirmed the laboratory cross-contamination. These results combined with all additional investigations done at the factory, allowed to release finished product batches produced and thus circumvented unnecessary food waste and economic losses for the factory.

Salmonella Serotyping; Comparison of the Traditional Method to a Microarray-Based Method and an in silico Platform Using Whole Genome Sequencing Data.

Salmonella is one of the most common causes of food-borne diseases worldwide. While Salmonella molecular subtyping by Whole Genome Sequencing (WGS) is increasingly used for outbreak and source tracking investigations, serotyping remains as a first-line characterization of Salmonella isolates. The traditional phenotypic method for serotyping is logistically challenging, as it requires the use of more than 150 specific antisera and well trained personnel to interpret the results. Consequently, it is not a routine method for the majority of laboratories. Several rapid molecular methods targeting O and H loci or surrogate genomic markers have been developed as alternative solutions. With the expansion of WGS, in silico Salmonella serotype prediction using WGS data is available. Here, we compared a microarray method using molecular markers, the Check and Trace Salmonella assay (CTS) and a WGS-based serotype prediction tool that targets molecular determinants of serotype (SeqSero) to the traditional phenotypic method using 100 strains representing 45 common and uncommon serotypes. Compared to the traditional method, the CTS assay correctly serotyped 97% of the strains, four strains gave a double serotype prediction. Among the inconclusive data, one strain was not predicted and two strains were incorrectly identified. SeqSero was evaluated with two versions (SeqSero 1 and the alpha test version of SeqSero 2). The correct antigenic formula was predicted by SeqSero 1 for 96 and 95% of strains using raw reads and assembly, respectively. However, 34 and 33% of these predictions included multiple serotypes by raw reads and assembly. With raw reads, one strain was not identified and three strains were discordant with phenotypic serotyping result. With assembly, three strains were not predicted and two strains were incorrectly predicted. While still under development, SeqSero 2 maintained the accuracy of antigenic formula prediction at 98% and reduced multiple serotype prediction rate to 13%. One strain had no prediction and one strain was incorrectly predicted. Our study indicates that the CTS assay is a good alternative for routine laboratories as it is an easy to use method with a short turn-around-time. SeqSero is a reliable replacement for phenotypic serotyping if WGS is routinely implemented.

The role of prophage for genome diversification within a clonal lineage of Lactobacillus johnsonii: characterization of the defective prophage LJ771.

Two independent isolates of the gut commensal Lactobacillus johnsonii were sequenced. These isolates belonged to the same clonal lineage and differed mainly by a 40.8-kb prophage, LJ771, belonging to the Sfi11 phage lineage. LJ771 shares close DNA sequence identity with Lactobacillus gasseri prophages. LJ771 coexists as an integrated prophage and excised circular phage DNA, but phage DNA packaged into extracellular phage particles was not detected. Between the phage lysin gene and attR a likely mazE ("antitoxin")/pemK ("toxin") gene cassette was detected in LJ771 but not in the L. gasseri prophages. Expressed pemK could be cloned in Escherichia coli only together with the mazE gene. LJ771 was shown to be highly stable and could be cured only by coexpression of mazE from a plasmid. The prophage was integrated into the methionine sulfoxide reductase gene (msrA) and complemented the 5' end of this gene, creating a protein with a slightly altered N-terminal sequence. The two L. johnsonii strains had identical in vitro growth and in vivo gut persistence phenotypes. Also, in an isogenic background, the presence of the prophage resulted in no growth disadvantage.

A Validation Approach of an End-to-End Whole Genome Sequencing Workflow for Source Tracking of Listeria monocytogenes and Salmonella enterica.

Whole genome sequencing (WGS), using high throughput sequencing technology, reveals the complete sequence of the bacterial genome in a few days. WGS is increasingly being used for source tracking, pathogen surveillance and outbreak investigation due to its high discriminatory power. In the food industry, WGS used for source tracking is beneficial to support contamination investigations. Despite its increased use, no standards or guidelines are available today for the use of WGS in outbreak and/or trace-back investigations. Here we present a validation of our complete (end-to-end) WGS workflow for Listeria monocytogenes and Salmonella enterica including: subculture of isolates, DNA extraction, sequencing and bioinformatics analysis. This end-to-end WGS workflow was evaluated according to the following performance criteria: stability, repeatability, reproducibility, discriminatory power, and epidemiological concordance. The current study showed that few single nucleotide polymorphism (SNPs) were observed for L. monocytogenes and S. enterica when comparing genome sequences from five independent colonies from the first subculture and five independent colonies after the tenth subculture. Consequently, the stability of the WGS workflow for L. monocytogenes and S. enterica was demonstrated despite the few genomic variations that can occur during subculturing steps. Repeatability and reproducibility were also demonstrated. The WGS workflow was shown to have a high discriminatory power and has the ability to show genetic relatedness. Additionally, the WGS workflow was able to reproduce published outbreak investigation results, illustrating its capability of showing epidemiological concordance. The current study proposes a validation approach comprising all steps of a WGS workflow and demonstrates that the workflow can be applied to L. monocytogenes or S. enterica.

Genome analysis of phage JS98 defines a fourth major subgroup of T4-like phages in Escherichia coli.

Numerous T4-like Escherichia coli phages were isolated from human stool and environmental wastewater samples in Bangladesh and Switzerland. The sequences of the major head gene (g23) revealed that these coliphages could be placed into four subgroups, represented by the phages T4, RB69, RB49, and JS98. Thus, JS98 defines a new major subgroup of E. coli T4-like phages. We conducted an analysis of the 169-kb JS98 genome sequence. Overall, 198 of the 266 JS98 open reading frames (ORFs) shared amino acid sequence identity with the reference T4 phage, 41 shared identity with other T4-like phages, and 27 ORFs lacked any database matches. Genes on the plus strand encoded virion proteins, which showed moderate to high sequence identity with T4 proteins. The right genome half of JS98 showed a higher degree of sequence conservation with T4 and RB69, even for the nonstructural genes, than did the left genome half, containing exclusively nonstructural genes. Most of the JS98-specific genes were found in the left genome half. Two came as a hypervariability cluster, but most represented isolated genes, suggesting that they were acquired separately in multiple acquisition events. No evidence for DNA exchange between JS98 phage and the E. coli host genome or coliphages other than T4 was observed. No undesired genes which could compromise its medical use were detected in the JS98 genome sequence.

Gene expression of commensal Lactobacillus johnsonii strain NCC533 during in vitro growth and in the murine gut.

Work with pathogens like Vibrio cholerae has shown major differences between genes expressed in bacteria grown in vitro and in vivo. To explore this subject for commensals, we investigated the transcription of the Lactobacillus johnsonii NCC533 genome during in vitro and in vivo growth using the microarray technology. During broth growth, 537, 626, and 277 of the 1,756 tested genes were expressed during exponential phase, "adaptation" (early stationary phase), and stationary phase, respectively. One hundred one, 150, and 33 genes, respectively, were specifically transcribed in these three phases. To explore the in vivo transcription program, we fed L. johnsonii containing a resistance plasmid to antibiotic-treated mice. After a 2-day washout phase, we determined the viable-cell counts of lactobacilli that were in the lumina and associated with the mucosae of different gut segments. While the cell counts showed a rather uniform distribution along the gut, we observed marked differences with respect to the expression of the Lactobacillus genome. The largest number of transcribed genes was in the stomach (n = 786); the next-largest numbers occurred in the cecum (n = 391) and the jejunum (n = 296), while only 26 Lactobacillus genes were transcribed in the colon. In vitro and in vivo transcription programs overlapped only partially. One hundred ninety-one of the transcripts from the lactobacilli in the stomach were not detected during in vitro growth; 202 and 213 genes, respectively, were transcribed under all in vitro and in vivo conditions; but the core transcriptome for all growth conditions comprised only 103 genes. Forty-four percent of the NCC533 genes were not detectably transcribed under any of the investigated conditions. Nontranscribed genes were clustered on the genome and enriched in the variable-genome part. Our data revealed not only major differences between in vitro- and in vivo-expressed genes in a Lactobacillus gut commensal organism but also marked changes in the expression of genes along the digestive tract.

Global transcriptome analysis of the heat shock response of Bifidobacterium longum.

Bifidobacteria are natural inhabitants of the human gastrointestinal tract and have been widely used as functional foods in different products. During industrial processing, bacterial cells undergo several stresses that can limit large-scale production and stability of the final product. To better understand the stress-response mechanisms of bifidobacteria, microarrays were used to obtain a global transcriptome profile of Bifidobacterium longum NCC2705 exposed to a heat shock treatment at 50 degrees C for 3, 7 and 12 min. Gene expression data highlighted a profound modification of gene expression, with 46% of the genes being altered. This analysis revealed a slow-down of Bi. longum general metabolic activity during stress with a simultaneous activation of the classical heat shock stimulon. Moreover, the expression of several genes with unknown function was highly induced under stress conditions. Three of these were conserved in other bacteria species where they were also previously shown to be induced by high temperature, suggesting their widespread role in the heat stress response. Finally, the implication of the trans-translation machinery in the response of Bi. longum cells to heat shock was suggested by the induction of the gene encoding the tmRNA-associated small protein B (SmpB) with concomitant high constitutive expression of the tmRNA gene.

Fully Closed Genome Sequences of Five Type Strains of the Genus Cronobacter and One Cronobacter sakazakii Strain.

Cronobacteris associated with infant infections and the consumption of reconstituted infant formula. Here we sequenced and closed six genomes ofC. condimenti(T),C. muytjensii(T),C. universalis(T),C. malonaticus(T),C. dublinensis(T), andC. sakazakiithat can be used as reference genomes in single nucleotide polymorphism (SNP)-based next-generation sequencing (NGS) analysis for source tracking investigations.

Genome Sequence of Lactobacillus fermentum Strain NCC2970 (CNCM I-5068).

Lactobacillus fermentum NCC2970 (CNCM I-5068) is a lactic acid bacterium originating from the Nestle Culture Collection. Here, we disclose its full 1.9-Gb genome sequence comprising one chromosome with no plasmid.

Safety analysis of a Russian phage cocktail: from metagenomic analysis to oral application in healthy human subjects.

Phage therapy has a long tradition in Eastern Europe, where preparations are comprised of complex phage cocktails whose compositions have not been described. We investigated the composition of a phage cocktail from the Russian pharmaceutical company Microgen targeting Escherichia coli/Proteus infections. Electron microscopy identified six phage types, with numerically T7-like phages dominating over T4-like phages. A metagenomic approach using taxonomical classification, reference mapping and de novo assembly identified 18 distinct phage types, including 7 genera of Podoviridae, 2 established and 2 proposed genera of Myoviridae, and 2 genera of Siphoviridae. De novo assembly yielded 7 contigs greater than 30 kb, including a 147-kb Myovirus genome and a 42-kb genome of a potentially new phage. Bioinformatic analysis did not reveal undesired genes and a small human volunteer trial did not associate adverse effects with oral phage exposure.

Oral T4-like phage cocktail application to healthy adult volunteers from Bangladesh.

The genomic diversity of 99 T4-like coliphages was investigated by sequencing an equimolar mixture with Illumina technology and screening them against different databases for horizontal gene transfer and undesired genes. A 9-phage cocktail was given to 15 healthy adults from Bangladesh at a dose of 3×10(9) and 3×10(7) plaque-forming units and placebo respectively. Phages were detected in 64% of the stool samples when subjects were treated with higher titer phage, compared to 30% and 28% with lower-titer phage and placebo, respectively. No Escherichia coli was present in initial stool samples, and no amplification of phage was observed. One percent of the administered oral phage was recovered from the feces. No adverse events were observed by self-report, clinical examination, or from laboratory tests for liver, kidney, and hematology function. No impact of oral phage was seen on the fecal microbiota composition with respect to bacterial 16S rRNA from stool.

T4 phages against Escherichia coli diarrhea: potential and problems.

A combination of in vitro and in vivo experiments with comparative phage genomics was used for the rational design of a phage cocktail against E. coli diarrhea. Orally applied T4 coliphages representing three different subgroups (T4-, RB49- and JS98-like phages) had no negative impact on the murine gut microbiota. T4 phages were found with high titers in the cecum and colon and lower titers in the small intestine, but were not detected in the blood, liver or spleen. No adverse effects were observed after one-month exposure to phage nor were serum anti-T4 antibodies detected. T4 phages belonging to the same subgroup showed closely related genomes that differed by 12 (phage JS10 vs. JS98 reference) to 17 (phage JSE vs. RB49 reference) insertion/deletions mostly representing single small ORFs. Bioinformatic analysis did not reveal undesired genes in the T4 genomes. Sequence variability was seen over the tail fibre genes, but the variability did not correlate with phage host range. The investigated T4 phages were not only species- but also strain-specific, necessitating the use of phage cocktails consisting of 10 and 16 T4 phage isolates to cover half to two thirds of E. coli strains representing the five main pathotypes isolated from diarrhea patients.

Similarity and differences in the Lactobacillus acidophilus group identified by polyphasic analysis and comparative genomics.

A set of lactobacilli were investigated by polyphasic analysis. Multilocus sequence analysis, DNA typing, microarray analysis, and in silico whole-genome alignments provided a remarkably consistent pattern of similarity within the Lactobacillus acidophilus complex. On microarray analysis, 17 and 5% of the genes from Lactobacillus johnsonii strain NCC533 represented variable and strain-specific genes, respectively, when tested against four independent isolates of L. johnsonii. When projected on the NCC533 genome map, about 10 large clusters of variable genes were identified, and they were enriched around the terminus of replication. A quarter of the variable genes and two-thirds of the strain-specific genes were associated with mobile DNA. Signatures for horizontal gene transfer and modular evolution were found in prophages and in DNA from the exopolysaccharide biosynthesis cluster. On microarray hybridizations, Lactobacillus gasseri strains showed a shift to significantly lower fluorescence intensities than the L. johnsonii test strains, and only genes encoding very conserved cellular functions from L. acidophilus hybridized to the L. johnsonii array. In-silico comparative genomics showed extensive protein sequence similarity and genome synteny of L. johnsonii with L. gasseri, L. acidophilus, and Lactobacillus delbrueckii; moderate synteny with Lactobacillus casei; and scattered X-type sharing of protein sequence identity with the other sequenced lactobacilli. The observation of a stepwise decrease in similarity between the members of the L. acidophilus group suggests a strong element of vertical evolution in a natural phylogenetic group. Modern whole-genome-based techniques are thus a useful adjunct to the clarification of taxonomical relationships in problematic bacterial groups.

The genome sequence of the probiotic intestinal bacterium Lactobacillus johnsonii NCC 533.

Lactobacillus johnsonii NCC 533 is a member of the acidophilus group of intestinal lactobacilli that has been extensively studied for their "probiotic" activities that include, pathogen inhibition, epithelial cell attachment, and immunomodulation. To gain insight into its physiology and identify genes potentially involved in interactions with the host, we sequenced and analyzed the 1.99-Mb genome of L. johnsonii NCC 533. Strikingly, the organism completely lacked genes encoding biosynthetic pathways for amino acids, purine nucleotides, and most cofactors. In apparent compensation, a remarkable number of uncommon and often duplicated amino acid permeases, peptidases, and phosphotransferase-type transporters were discovered, suggesting a strong dependency of NCC 533 on the host or other intestinal microbes to provide simple monomeric nutrients. Genome analysis also predicted an abundance (>12) of large and unusual cell-surface proteins, including fimbrial subunits, which may be involved in adhesion to glycoproteins or other components of mucin, a characteristic expected to affect persistence in the gastrointestinal tract (GIT). Three bile salt hydrolases and two bile acid transporters, proteins apparently critical for GIT survival, were also detected. In silico genome comparisons with the >95% complete genome sequence of the closely related Lactobacillus gasseri revealed extensive synteny punctuated by clear-cut insertions or deletions of single genes or operons. Many of these regions of difference appear to encode metabolic or structural components that could affect the organisms competitiveness or interactions with the GIT ecosystem.