A tiny fraction of all species forms most of nature: Rarity as a sticky state.

Using data from a wide range of natural communities including the human microbiome, plants, fish, mushrooms, rodents, beetles, and trees, we show that universally just a few percent of the species account for most of the biomass. This is in line with the classical observation that the vast bulk of biodiversity is very rare. Attempts to find traits allowing the tiny fraction of abundant species to escape rarity have remained unsuccessful. Here, we argue that this might be explained by the fact that hyper-dominance can emerge through stochastic processes. We demonstrate that in neutrally competing groups of species, rarity tends to become a trap if environmental fluctuations result in gains and losses proportional to abundances. This counter-intuitive phenomenon arises because absolute change tends to zero for very small abundances, causing rarity to become a "sticky state", a pseudoattractor that can be revealed numerically in classical ball-in-cup landscapes. As a result, the vast majority of species spend most of their time in rarity leaving space for just a few others to dominate the neutral community. However, fates remain stochastic. Provided that there is some response diversity, roles occasionally shift as stochastic events or natural enemies bring an abundant species down allowing a rare species to rise to dominance. Microbial time series spanning thousands of generations support this prediction. Our results suggest that near-neutrality within niches may allow numerous rare species to persist in the wings of the dominant ones. Stand-ins may serve as insurance when former key species collapse.

Frailty is related to serum inflammageing markers: results from the VITAL study.

Frailty describes an age-associated state in individuals with an increased vulnerability and less resilience against adverse outcomes. To score frailty, studies have employed the questionnaires, such as the SF-36 and EQ-5D-3L, or the Frailty Index, a composite score based on deficit accumulation. Furthermore, ageing of the immune system is often accompanied by a state of low-grade inflammation (inflammageing). Here, we aimed to associate 29 circulating markers of inflammageing with frailty measures in a prospective cohort study to understand the mechanisms underlying ageing.Frailty measures and inflammageing markers were assessed in 317 participants aged 25-90. We determined four different measures of frailty: the Frailty Index based on 31 deficits, the EQ-5D-3L and two physical domains of the SF-36. Serum/plasma levels of inflammageing markers and CMV/EBV seropositivity were measured using different techniques: Quanterix, Luminex or ELISA.All four measures of frailty strongly correlated with age and BMI. Nineteen biomarkers correlated with age, some in a linear fashion (IL-6, YKL-40), some only in the oldest age brackets (CRP), and some increased at younger ages and then plateaued (CCL2, sIL-6R). After correcting for age, biomarkers, such as IL-6, CRP, IL-1RA, YKL-40 and elastase, were associated with frailty. When corrected for BMI, the number of associations reduced further.In conclusion, inflammageing markers, particularly markers reflecting innate immune activation, are related to frailty. These findings indicate that health decline and the accumulation of deficits with age is accompanied with a low-grade inflammation which can be detected by specific inflammatory markers.

Modulation of human intestinal microbiota in a clinical trial by consumption of a ß-D-glucan-enriched extract obtained from Lentinula edodes.

PURPOSE: The aim of this study was to evaluate the hypocholesterolemic, immune- and microbiota-modulatory effect of a mushroom extract in hypercholesterolemic subjects. METHODS: A randomized, controlled, double-blind, and parallel clinical trial was carried out with subjects from 18 to 65 years old (n = 52) with untreated mild hypercholesterolemia. Volunteers consumed a ß-D-glucan-enriched (BGE) mixture (10.4 g/day) obtained from shiitake mushrooms (Lentinula edodes) ensuring a 3.5 g/day of fungal ß-D-glucans or a placebo incorporated in three different commercial creams. RESULTS: This mixture showed hypocholesterolemic activities in vitro and in animal studies. After eight weeks intervention, no significant differences in lipid- or cholesterol-related parameters were found compared to placebo subjects as well as before and after the BGE mixture administration. No inflammatory or immunomodulatory responses were noticed and no changes in IL-1ß, IL-6, TNF-α or oxLDL were recorded. However, consumption of the BGE mixture was safe and managed to achieve the dietary fibre intake recommended as cardiovascular protective diet. Moreover, the BGE mixture modulated the colonic microbiota differently compared to placebo. Microbial community composition varied from before to after the intervention with several genera being positively or negatively correlated with some biomarkers related to cholesterol metabolism. CONCLUSION: These results suggested a relation between cholesterol metabolism, microbiota and BGE administration. Nevertheless, the precise significance of this differential modulation was not fully elucidated and requires further studies.

Unravelling lactate-acetate and sugar conversion into butyrate by intestinal Anaerobutyricum and Anaerostipes species by comparative proteogenomics.

The d- and l-forms of lactate are important fermentation metabolites produced by intestinal bacteria but are found to negatively affect mucosal barrier function and human health. Both enantiomers of lactate can be converted with acetate into the presumed beneficial butyrate by a phylogenetically related group of anaerobes, including Anaerobutyricum and Anaerostipes spp. This is a low energy yielding process with a partially unknown pathway in Anaerobutyricum and Anaerostipes spp. and hence, we sought to address this via a comparative genomics, proteomics and physiology approach. We compared growth of Anaerobutyricum soehngenii on lactate with that on sucrose and sorbitol. Comparative proteomics revealed complete pathway of butyrate formation from sucrose, sorbitol and lactate. Notably, a gene cluster, lctABCDEF was abundantly expressed when grown on lactate. This gene cluster encodes a lactate dehydrogenase (lctD), electron transport proteins A and B (lctCB), nickel-dependent racemase (lctE), lactate permease (lctF) and short-chain acyl-CoA dehydrogenase (lctG). Investigation of available genomes of intestinal bacteria revealed this new gene cluster to be highly conserved in only Anaerobutyricum and Anaerostipes spp. Present study demonstrates that A. soehngenii and several related Anaerobutyricum and Anaerostipes spp. are highly adapted for a lifestyle involving lactate plus acetate utilization in the human intestinal tract.

Molecular ecology of the yet uncultured bacterial Ct85-cluster in the mammalian gut.

In our previous studies on irritable bowel syndrome (IBS) -associated microbiota by molecular methods, we demonstrated that a particular 16S rRNA gene amplicon was more abundant in the feces of healthy subjects or mixed type IBS (IBS-M) -sufferers than in the feces of individuals with diarrhea-type IBS (IBS-D). In the current study, we demonstrated that this, so called Ct85-amplicon, consists of a cluster of very heterogeneous 16S rRNA gene sequences, and defined six 16S rRNA gene types, a to f, within this cluster, each representing a novel species-, genus- or family level taxon. We then designed specific PCR primers for these sequence types, mapped the distribution of the Ct85-cluster sequences and that of the newly defined sequence types in several animal species and compared the sequence types present in the feces of healthy individuals and IBS sufferers using two IBS study cohorts, Finnish and Dutch. Various Ct85-cluster sequence types were detected in the fecal samples of several companion and production animal species with remarkably differing prevalences and abundances. The Ct85 sequence type composition of swine closely resembled that of humans. One of the five types (d) shared between humans and swine was not present in any other animals tested, while one sequence type (b) was found only in human samples. In both IBS study cohorts, one type (e) was more prevalent in healthy individuals than in the IBS-M group. By revealing various sequence types in the widespread Ct85-cluster and their distribution, the results improve our understanding of these uncultured bacteria, which is essential for future efforts to cultivate representatives of the Ct85-cluster and reveal their roles in IBS.

Genomic and physiological analyses of an indigenous strain, Enterococcus faecium 17OM39.

The human gut microbiome plays a crucial role in human health and efforts need to be done for cultivation and characterisation of bacteria with potential health benefits. Here, we isolated a bacterium from a healthy Indian adult faeces and investigated its potential as probiotic. The cultured bacterial strain 17OM39 was identified as Enterococcus faecium by 16S rRNA gene sequencing. The strain 17OM39 exhibited tolerance to acidic pH, showed antimicrobial activity and displayed strong cell surface traits such as hydrophobicity and autoaggregation capacity. The strain was able to tolerate bile salts and showed bile salt hydrolytic (BSH) activity, exopolysaccharide production and adherence to human HT-29 cell line. Importantly, partial haemolytic activity was detected and the strain was susceptible to the human serum. Genomics investigation of strain 17OM39 revealed the presence of diverse genes encoding for proteolytic enzymes, stress response systems and the ability to produce essential amino acids, vitamins and antimicrobial compound Bacteriocin-A. No virulence factors and plasmids were found in this genome of the strain 17OM39. Collectively, these physiological and genomic features of 17OM39 confirm the potential of this strain as a candidate probiotic.

Reclassification of Eubacterium hallii as Anaerobutyricum hallii gen. nov., comb. nov., and description of Anaerobutyricum soehngenii sp. nov., a butyrate and propionate-producing bacterium from infant faeces.

A bacterial strain designated L2-7T, phylogenetically related to Eubacterium hallii DSM 3353T, was previously isolated from infant faeces. The complete genome of strain L2-7T contains eight copies of the 16S rRNA gene with only 98.0-98.5 % similarity to the 16S rRNA gene of the previously described type strain E. hallii. The next closest validly described species is Anaerostipes hadrus DSM 3319T (90.7 % 16S rRNA gene similarity). A polyphasic taxonomic approach showed strain L2-7T to be a novel species, related to type strain E. hallii DSM 3353T. The experimentally observed DNA-DNA hybridization value between strain L2-7T and E. hallii DSM 3353T was 26.25 %, close to that calculated from the genomes (34.3 %). The G+C content of the chromosomal DNA of strain L2-7T was 38.6 mol%. The major fatty acids were C16 : 0, C16 : 1cis9 and a component with summed feature 10 (C18 : 1c11/t9/t6c). Strain L2-7T had higher amounts of C16 : 0 (30.6 %) compared to E. hallii DSM 3353T (19.5 %) and its membrane contained phosphatidylglycerol and phosphatidylethanolamine, which were not detected in E. hallii DSM 3353T. Furthermore, 16S rRNA gene phylogenetic analysis advocates that E. hallii DSM 3353T is misclassified, and its reclassification as a member of the family Lachnospiraceae is necessary. Using a polyphasic approach, we propose that E. hallii (=DSM 3353T=ATCC 27751T) be reclassified as the type strain of a novel genus Anaerobutyricum sp. nov., comb. nov. and we propose that strain L2-7T should be classified as a novel species, Anaerobutyricum soehngenii sp. nov. The type strain is L2-7T (=DSM 17630T=KCTC 15707T).

An Inducible Operon Is Involved in Inulin Utilization in Lactobacillus plantarum Strains, as Revealed by Comparative Proteogenomics and Metabolic Profiling.

The draft genomes of Lactobacillus plantarum strains isolated from Asian fermented foods, infant feces, and shrimp intestines were sequenced and compared to those of well-studied strains. Among 28 strains of L. plantarum, variations in the genomic features involved in ecological adaptation were elucidated. The genome sizes ranged from approximately 3.1 to 3.5 Mb, of which about 2,932 to 3,345 protein-coding sequences (CDS) were predicted. The food-derived isolates contained a higher number of carbohydrate metabolism-associated genes than those from infant feces. This observation correlated to their phenotypic carbohydrate metabolic profile, indicating their ability to metabolize the largest range of sugars. Surprisingly, two strains (P14 and P76) isolated from fermented fish utilized inulin. ß-Fructosidase, the inulin-degrading enzyme, was detected in the supernatants and cell wall extracts of both strains. No activity was observed in the cytoplasmic fraction, indicating that this key enzyme was either membrane-bound or extracellularly secreted. From genomic mining analysis, a predicted inulin operon of fosRABCDXE, which encodes ß-fructosidase and many fructose transporting proteins, was found within the genomes of strains P14 and P76. Moreover, pts1BCA genes, encoding sucrose-specific IIBCA components involved in sucrose transport, were also identified. The proteomic analysis revealed the mechanism and functional characteristic of the fosRABCDXE operon involved in the inulin utilization of L. plantarum The expression levels of the fos operon and pst genes were upregulated at mid-log phase. FosE and the LPXTG-motif cell wall anchored ß-fructosidase were induced to a high abundance when inulin was present as a carbon source. IMPORTANCE: Inulin is a long-chain carbohydrate that may act as a prebiotic, which provides many health benefits to the host by selectively stimulating the growth and activity of beneficial bacteria in the colon. While certain lactobacilli can catabolize inulin, this has not yet been described for Lactobacillus plantarum, and an associated putative inulin operon has not been reported in this species. By using comparative and functional genomics, we showed that two L. plantarum strains utilized inulin and identified functional inulin operons in their genomes. The proteogenomic data revealed that inulin degradation and uptake routes, which related to the fosRABCDXE operon and pstBCA genes, were widely expressed among L. plantarum strains. The present work provides a novel understanding of gene regulation and mechanisms of inulin utilization in probiotic L. plantarum generating opportunities for synbiotic product development.

Draft Genome Sequences of Yersinia pestis Strains from the 1994 Plague Epidemic of Surat and 2002 Shimla Outbreak in India.

We report the first draft genome sequences of the strains of plague-causing bacteria, Yersinia pestis, from India. These include two strains from the Surat epidemic (1994), one strain from the Shimla outbreak (2002) and one strain from the plague surveillance activity in the Deccan plateau region (1998). Genome size for all four strains is ~4.49 million bp with 139-147 contigs. Average sequencing depth for all four genomes was 21x.

Pre-vaccination immunotypes reveal weak and robust antibody responders to influenza vaccination.

Effective vaccine-induced immune responses are particularly essential in older adults who face an increased risk of immunosenescence. However, the complexity and variability of the human immune system make predicting vaccine responsiveness challenging. To address this knowledge gap, our study aimed to characterize immune profiles that are predictive of vaccine responsiveness using "immunotypes" as an innovative approach. We analyzed an extensive set of innate and adaptive immune cell subsets in the whole blood of 307 individuals (aged 25-92) pre- and post-influenza vaccination which we associated with day 28 hemagglutination inhibition (HI) antibody titers. Building on our previous work that stratified individuals into nine immunotypes based on immune cell subsets, we identified two pre-vaccination immunotypes associated with weak and one showing robust day 28 antibody response. Notably, the weak responders demonstrated HLA-DR+ T-cell signatures, while the robust responders displayed a high naïve-to-memory T-cell ratio and percentage of nonclassical monocytes. These specific signatures deepen our understanding of the relationship between the baseline of the immune system and its functional potential. This approach could enhance our ability to identify individuals at risk of immunosenescence. Our findings highlight the potential of pre-vaccination immunotypes as an innovative tool for informing personalized vaccination strategies and improving health outcomes, particularly for aging populations.

Profiling the fecal microbiome and its modulators across the lifespan in the Netherlands.

Defining what constitutes a healthy microbiome throughout our lives remains an ongoing challenge. Understanding to what extent host and environmental factors can influence it has been the primary motivation for large population studies worldwide. Here, we describe the fecal microbiome of 3,746 individuals (0-87 years of age) in a nationwide study in the Netherlands, in association with extensive questionnaires. We validate previous findings, such as infant-adult trajectories, and explore the collective impact of our variables, which explain over 40% of the variation in microbiome composition. We identify associations with less explored factors, particularly those ethnic related, which show the largest impact on the adult microbiome composition, diversity, metabolic profiles, and CAZy (carbohydrate-active enzyme) repertoires. Understanding the sources of microbiome variability is crucial, given its potential as a modifiable target with therapeutic possibilities. With this work, we aim to serve as a foundational element for the design of health interventions and fundamental research.

Multiple vaccine comparison in the same adults reveals vaccine-specific and age-related humoral response patterns: an open phase IV trial.

Vaccine responsiveness is often reduced in older adults. Yet, our lack of understanding of low vaccine responsiveness hampers the development of effective vaccination strategies to reduce the impact of infectious diseases in the ageing population. Young-adult (25-49 y), middle-aged (50-64 y) and older-adult ( ≥ 65 y) participants of the VITAL clinical trials (n = 315, age-range: 28-98 y), were vaccinated with an annual (2019-2020) quadrivalent influenza (QIV) booster vaccine, followed by a primary 13-valent pneumococcal-conjugate (PCV13) vaccine (summer/autumn 2020) and a primary series of two SARS-CoV-2 mRNA-1273 vaccines (spring 2021). This unique setup allowed investigation of humoral responsiveness towards multiple vaccines within the same individuals over the adult age-range. Booster QIV vaccination induced comparable H3N2 hemagglutination inhibition (HI) titers in all age groups, whereas primary PCV13 and mRNA-1273 vaccination induced lower antibody concentrations in older as compared to younger adults (primary endpoint). The persistence of humoral responses, towards the 6 months timepoint, was shorter in older adults for all vaccines (secondary endpoint). Interestingly, highly variable vaccine responder profiles overarching multiple vaccines were observed. Yet, approximately 10% of participants, mainly comprising of older male adults, were classified as low responders to multiple vaccines. This study aids the identification of risk groups for low vaccine responsiveness and hence supports targeted vaccination strategies. Trial number: NL69701.041.19, EudraCT: 2019-000836-24.

A tool to assess the mock community samples in 16S rRNA gene-based microbiota profiling studies.

Inclusion and investigation of technical controls in microbiome sequencing studies is important for understanding technical biases and errors. Here, we present chkMocks, a general R-based tool that allows researchers to compare the composition of mock communities that are processed along with samples to their theoretical composition. A visual comparison between experimental and theoretical community composition and their correlation is provided for researchers to assess the quality of their sample processing workflows.

Reducing bias in microbiome research: Comparing methods from sample collection to sequencing.

Background: Microbiota profiles are strongly influenced by many technical aspects that impact the ability of researchers to compare results. To investigate and identify potential biases introduced by technical variations, we compared several approaches throughout the entire workflow of a microbiome study, from sample collection to sequencing, using commercially available mock communities (from bacterial strains as well as from DNA) and multiple human fecal samples, including a large set of positive controls created as a random mix of several participant samples. Methods: Human fecal material was sampled, and aliquots were used to test two commercially available stabilization solutions (OMNIgene·GUT and Zymo Research) in comparison to samples frozen immediately upon collection. In addition, the methodology for DNA extraction, input of DNA, or the number of PCR cycles were analyzed. Furthermore, to investigate the potential batch effects in DNA extraction, sequencing, and barcoding, we included 139 positive controls. Results: Samples preserved in both the stabilization buffers limited the overgrowth of Enterobacteriaceae when compared to unpreserved samples stored at room temperature (RT). These stabilized samples stored at RT were different from immediately frozen samples, where the relative abundance of Bacteroidota was higher and Actinobacteriota and Firmicutes were lower. As reported previously, the method used for cell disruption was a major contributor to variation in microbiota composition. In addition, a high number of cycles during PCR lead to an increase in contaminants detected in the negative controls. The DNA extraction had a significant impact on the microbial composition, also observed with the use of different Illumina barcodes during library preparation and sequencing, while no batch effect was observed in replicate runs. Conclusion: Our study reaffirms the importance of the mechanical cell disruption method and immediate frozen storage as critical aspects in fecal microbiota studies. A comparison of storage conditions revealed that the bias was limited in RT samples preserved in stabilization systems, and these may be a suitable compromise when logistics are challenging due to the size or location of a study. Moreover, to reduce the effect of contaminants in fecal microbiota profiling studies, we suggest the use of ~125 pg input DNA and 25 PCR cycles as optimal parameters during library preparation.

Wavelet clustering analysis as a tool for characterizing community structure in the human microbiome.

Human microbiome research is helped by the characterization of microbial networks, as these may reveal key microbes that can be targeted for beneficial health effects. Prevailing methods of microbial network characterization are based on measures of association, often applied to limited sampling points in time. Here, we demonstrate the potential of wavelet clustering, a technique that clusters time series based on similarities in their spectral characteristics. We illustrate this technique with synthetic time series and apply wavelet clustering to densely sampled human gut microbiome time series. We compare our results with hierarchical clustering based on temporal correlations in abundance, within and across individuals, and show that the cluster trees obtained by using either method are significantly different in terms of elements clustered together, branching structure and total branch length. By capitalizing on the dynamic nature of the human microbiome, wavelet clustering reveals community structures that remain obscured in correlation-based methods.

Unveiling and harnessing the human gut microbiome in the rising burden of non-communicable diseases during urbanization.

The world is witnessing a global increase in the urban population, particularly in developing Asian and African countries. Concomitantly, the global burden of non-communicable diseases (NCDs) is rising, markedly associated with the changing landscape of lifestyle and environment during urbanization. Accumulating studies have revealed the role of the gut microbiome in regulating the immune and metabolic homeostasis of the host, which potentially bridges external factors to the host (patho-)physiology. In this review, we discuss the rising incidences of NCDs during urbanization and their links to the compositional and functional dysbiosis of the gut microbiome. In particular, we elucidate the effects of urbanization-associated factors (hygiene/pollution, urbanized diet, lifestyles, the use of antibiotics, and early life exposure) on the gut microbiome underlying the pathogenesis of NCDs. We also discuss the potential and feasibility of microbiome-inspired and microbiome-targeted approaches as novel avenues to counteract NCDs, including fecal microbiota transplantation, diet modulation, probiotics, postbiotics, synbiotics, celobiotics, and precision antibiotics.

Changes in human gut flora with age: an Indian familial study.

BACKGROUND: The gut micro flora plays vital role in health status of the host. The majority of microbes residing in the gut have a profound influence on human physiology and nutrition. Different human ethnic groups vary in genetic makeup as well as the environmental conditions they live in. The gut flora changes with genetic makeup and environmental factors and hence it is necessary to understand the composition of gut flora of different ethnic groups. Indian population is different in physiology from western population (YY paradox) and thus the gut flora in Indian population is likely to differ from the extensively studied gut flora in western population. In this study we have investigated the gut flora of two Indian families, each with three individuals belonging to successive generations and living under the same roof. RESULTS: Denaturation gradient gel electrophoresis analysis showed age-dependant variation in gut microflora amongst the individuals within a family. Different bacterial genera were dominant in the individual of varying age in clone library analysis. Obligate anaerobes isolated from individuals within a family showed age related differences in isolation pattern, with 27% (6 out of 22) of the isolates being potential novel species based on 16S rRNA gene sequence. In qPCR a consistent decrease in Firmicutes number and increase in Bacteroidetes number with increasing age was observed in our subjects, this pattern of change in Firmicutes / Bacteroidetes ratio with age is different than previously reported in European population. CONCLUSION: There is change in gut flora with age amongst the individuals within a family. The isolation of high percent of novel bacterial species and the pattern of change in Firmicutes /Bacteroidetes ratio with age suggests that the composition of gut flora in Indian individuals may be different than the western population. Thus, further extensive study is needed to define the gut flora in Indian population.

Dynamic metabolic interactions and trophic roles of human gut microbes identified using a minimal microbiome exhibiting ecological properties.

Microbe-microbe interactions in the human gut are influenced by host-derived glycans and diet. The high complexity of the gut microbiome poses a major challenge for unraveling the metabolic interactions and trophic roles of key microbes. Synthetic minimal microbiomes provide a pragmatic approach to investigate their ecology including metabolic interactions. Here, we rationally designed a synthetic microbiome termed Mucin and Diet based Minimal Microbiome (MDb-MM) by taking into account known physiological features of 16 key bacteria. We combined 16S rRNA gene-based composition analysis, metabolite measurements and metatranscriptomics to investigate community dynamics, stability, inter-species metabolic interactions and their trophic roles. The 16 species co-existed in the in vitro gut ecosystems containing a mixture of complex substrates representing dietary fibers and mucin. The triplicate MDb-MM's followed the Taylor's power law and exhibited strikingly similar ecological and metabolic patterns. The MDb-MM exhibited resistance and resilience to temporal perturbations as evidenced by the abundance and metabolic end products. Microbe-specific temporal dynamics in transcriptional niche overlap and trophic interaction network explained the observed co-existence in a competitive minimal microbiome. Overall, the present study provides crucial insights into the co-existence, metabolic niches and trophic roles of key intestinal microbes in a highly dynamic and competitive in vitro ecosystem.

Inter-species Metabolic Interactions in an In-vitro Minimal Human Gut Microbiome of Core Bacteria.

Knowledge of the functional roles and interspecies interactions are crucial for improving our understanding of the human intestinal microbiome in health and disease. However, the complexity of the human intestinal microbiome and technical challenges in investigating it pose major challenges. In this proof-of-concept study, we rationally designed, assembled and experimentally tested a synthetic Diet-based Minimal Microbiome (Db-MM) consisting of ten core intestinal bacterial species that together are capable of efficiently converting dietary fibres into short chain fatty acids (SCFAs). Despite their genomic potential for metabolic competition, all ten bacteria coexisted during growth on a mixture of dietary fibres, including pectin, inulin, xylan, cellobiose and starch. By integrated analyses of metabolite production, community composition and metatranscriptomics-based gene expression data, we identified interspecies metabolic interactions leading to production of key SCFAs such as butyrate and propionate. While public goods, such as sugars liberated from colonic fibres, are harvested by non-degraders, some species thrive by cross-feeding on energetically challenging substrates, including the butyrogenic conversion of acetate and lactate. Using a reductionist approach in an in-vitro system combined with functional measurements, our study provides key insights into the complex interspecies metabolic interactions between core intestinal bacterial species.

An in vitro fermentation model to study the impact of bacteriophages targeting Shiga toxin-encoding Escherichia coli on the colonic microbiota.

Lytic bacteriophages are considered safe for human consumption as biocontrol agents against foodborne pathogens, in particular in ready-to-eat foodstuffs. Phages could, however, evolve to infect different hosts when passing through the gastrointestinal tract (GIT). This underlines the importance of understanding the impact of phages towards colonic microbiota, particularly towards bacterial families usually found in the colon such as the Enterobacteriaceae. Here we propose in vitro batch fermentation as model for initial safety screening of lytic phages targeting Shiga toxin-producing Escherichia coli (STEC). As inoculum we used faecal material of three healthy donors. To assess phage safety, we monitored fermentation parameters, including short chain fatty acid production and gas production/intake by colonic microbiota. We performed shotgun metagenomic analysis to evaluate the outcome of phage interference with colonic microbiota composition and functional potential. During the 24 h incubation, concentrations of phage and its host were also evaluated. We found the phage used in this study, named E. coli phage vB_EcoS_Ace (Ace), to be safe towards human colonic microbiota, independently of the donors' faecal content used. This suggests that individuality of donor faecal microbiota did not interfere with phage effect on the fermentations. However, the model revealed that the attenuated STEC strain used as phage host perturbed the faecal microbiota as based on metagenomic analysis, with potential differences in metabolic output. We conclude that the in vitro batch fermentation model used in this study is a reliable safety screening for lytic phages intended to be used as biocontrol agents.