MSPminer: abundance-based reconstitution of microbial pan-genomes from shotgun metagenomic data.

MOTIVATION: Analysis toolkits for shotgun metagenomic data achieve strain-level characterization of complex microbial communities by capturing intra-species gene content variation. Yet, these tools are hampered by the extent of reference genomes that are far from covering all microbial variability, as many species are still not sequenced or have only few strains available. Binning co-abundant genes obtained from de novo assembly is a powerful reference-free technique to discover and reconstitute gene repertoire of microbial species. While current methods accurately identify species core parts, they miss many accessory genes or split them into small gene groups that remain unassociated to core clusters. RESULTS: We introduce MSPminer, a computationally efficient software tool that reconstitutes Metagenomic Species Pan-genomes (MSPs) by binning co-abundant genes across metagenomic samples. MSPminer relies on a new robust measure of proportionality coupled with an empirical classifier to group and distinguish not only species core genes but accessory genes also. Applied to a large scale metagenomic dataset, MSPminer successfully delineates in a few hours the gene repertoires of 1661 microbial species with similar specificity and higher sensitivity than existing tools. The taxonomic annotation of MSPs reveals microorganisms hitherto unknown and brings coherence in the nomenclature of the species of the human gut microbiota. The provided MSPs can be readily used for taxonomic profiling and biomarkers discovery in human gut metagenomic samples. In addition, MSPminer can be applied on gene count tables from other ecosystems to perform similar analyses. AVAILABILITY AND IMPLEMENTATION: The binary is freely available for non-commercial users at www.enterome.com/downloads. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Quality control of microbiota metagenomics by k-mer analysis.

BACKGROUND: The biological and clinical consequences of the tight interactions between host and microbiota are rapidly being unraveled by next generation sequencing technologies and sophisticated bioinformatics, also referred to as microbiota metagenomics. The recent success of metagenomics has created a demand to rapidly apply the technology to large case-control cohort studies and to studies of microbiota from various habitats, including habitats relatively poor in microbes. It is therefore of foremost importance to enable a robust and rapid quality assessment of metagenomic data from samples that challenge present technological limits (sample numbers and size). Here we demonstrate that the distribution of overlapping k-mers of metagenome sequence data predicts sequence quality as defined by gene distribution and efficiency of sequence mapping to a reference gene catalogue. RESULTS: We used serial dilutions of gut microbiota metagenomic datasets to generate well-defined high to low quality metagenomes. We also analyzed a collection of 52 microbiota-derived metagenomes. We demonstrate that k-mer distributions of metagenomic sequence data identify sequence contaminations, such as sequences derived from "empty" ligation products. Of note, k-mer distributions were also able to predict the frequency of sequences mapping to a reference gene catalogue not only for the well-defined serial dilution datasets, but also for 52 human gut microbiota derived metagenomic datasets. CONCLUSIONS: We propose that k-mer analysis of raw metagenome sequence reads should be implemented as a first quality assessment prior to more extensive bioinformatics analysis, such as sequence filtering and gene mapping. With the rising demand for metagenomic analysis of microbiota it is crucial to provide tools for rapid and efficient decision making. This will eventually lead to a faster turn-around time, improved analytical quality including sample quality metrics and a significant cost reduction. Finally, improved quality assessment will have a major impact on the robustness of biological and clinical conclusions drawn from metagenomic studies.

Driving gut microbiota enterotypes through host genetics.

BACKGROUND: Population stratification based on interindividual variability in gut microbiota composition has revealed the existence of several ecotypes named enterotypes in humans and various animal species. Enterotypes are often associated with environmental factors including diet, but knowledge of the role of host genetics remains scarce. Moreover, enterotypes harbor functionalities likely associated with varying abilities and susceptibilities of their host. Previously, we showed that under controlled conditions, 60-day-old pig populations consistently split into two enterotypes with either Prevotella and Mitsuokella (PM enterotype) or Ruminococcus and Treponema (RT enterotype) as keystone taxa. Here, our aim was to rely on pig as a model to study the influence of host genetics to assemble enterotypes, and to provide clues on enterotype functional differences and their links with growth traits. RESULTS: We established two pig lines contrasted for abundances of the genera pairs specifying each enterotype at 60 days of age and assessed them for fecal microbiota composition and growth throughout three consecutive generations. Response to selection across three generations revealed, per line, an increase in the prevalence of the selected enterotype and in the average relative abundances of directly and indirectly selected bacterial genera. The PM enterotype was found less diverse than the RT enterotype but more efficient for piglet growth during the post-weaning period. Shotgun metagenomics revealed differentially abundant bacterial species between the two enterotypes. By using the KEGG Orthology database, we show that functions related to starch degradation and polysaccharide metabolism are enriched in the PM enterotype, whereas functions related to general nucleoside transport and peptide/nickel transport are enriched in the RT enterotype. Our results also suggest that the PM and RT enterotypes might differ in the metabolism of valine, leucin, and isoleucine, favoring their biosynthesis and degradation, respectively. CONCLUSION: We experimentally demonstrated that enterotypes are functional ecosystems that can be selected as a whole by exerting pressure on the host genetics. We also highlight that holobionts should be considered as units of selection in breeding programs. These results pave the way for a holistic use of host genetics, microbiota diversity, and enterotype functionalities to understand holobiont shaping and adaptation. Video Abstract.

Four functional profiles for fibre and mucin metabolism in the human gut microbiome.

BACKGROUND: With the emergence of metagenomic data, multiple links between the gut microbiome and the host health have been shown. Deciphering these complex interactions require evolved analysis methods focusing on the microbial ecosystem functions. Despite the fact that host or diet-derived fibres are the most abundant nutrients available in the gut, the presence of distinct functional traits regarding fibre and mucin hydrolysis, fermentation and hydrogenotrophic processes has never been investigated. RESULTS: After manually selecting 91 KEGG orthologies and 33 glycoside hydrolases further aggregated in 101 functional descriptors representative of fibre and mucin degradation pathways in the gut microbiome, we used nonnegative matrix factorization to mine metagenomic datasets. Four distinct metabolic profiles were further identified on a training set of 1153 samples, thoroughly validated on a large database of 2571 unseen samples from 5 external metagenomic cohorts and confirmed with metatranscriptomic data. Profiles 1 and 2 are the main contributors to the fibre-degradation-related metagenome: they present contrasted involvement in fibre degradation and sugar metabolism and are differentially linked to dysbiosis, metabolic disease and inflammation. Profile 1 takes over Profile 2 in healthy samples, and unbalance of these profiles characterize dysbiotic samples. Furthermore, high fibre diet favours a healthy balance between profiles 1 and profile 2. Profile 3 takes over profile 2 during Crohn's disease, inducing functional reorientations towards unusual metabolism such as fucose and H2S degradation or propionate, acetone and butanediol production. Profile 4 gathers under-represented functions, like methanogenesis. Two taxonomic makes up of the profiles were investigated, using either the covariation of 203 prevalent genomes or metagenomic species, both providing consistent results in line with their functional characteristics. This taxonomic characterization showed that profiles 1 and 2 were respectively mainly composed of bacteria from the phyla Bacteroidetes and Firmicutes while profile 3 is representative of Proteobacteria and profile 4 of methanogens. CONCLUSIONS: Integrating anaerobic microbiology knowledge with statistical learning can narrow down the metagenomic analysis to investigate functional profiles. Applying this approach to fibre degradation in the gut ended with 4 distinct functional profiles that can be easily monitored as markers of diet, dysbiosis, inflammation and disease. Video Abstract.

In Vitro Modelling of Oral Microbial Invasion in the Human Colon.

Recent advances in the human microbiome characterization have revealed significant oral microbial detection in stools of dysbiotic patients. However, little is known about the potential interactions of these invasive oral microorganisms with commensal intestinal microbiota and the host. In this proof-of-concept study, we proposed a new model of oral-to-gut invasion by the combined use of an in vitro model simulating both the physicochemical and microbial (lumen- and mucus-associated microbes) parameters of the human colon (M-ARCOL), a salivary enrichment protocol, and whole-metagenome shotgun sequencing. Oral invasion of the intestinal microbiota was simulated by injection of enriched saliva in the in vitro colon model inoculated with a fecal sample from the same healthy adult donor. The mucosal compartment of M-ARCOL was able to retain the highest species richness levels over time, while species richness levels decreased in the luminal compartment. This study also showed that oral microorganisms preferably colonized the mucosal microenvironment, suggesting potential oral-to-intestinal mucosal competitions. This new model of oral-to-gut invasion can provide useful mechanistic insights into the role of oral microbiome in various disease processes. IMPORTANCE Here, we propose a new model of oral-to-gut invasion by the combined use of an in vitro model simulating both the physicochemical and microbial (lumen- and mucus-associated microbes) parameters of the human colon (M-ARCOL), a salivary enrichment protocol, and whole-metagenome shotgun sequencing. Our study revealed the importance of integrating the mucus compartment, which retained higher microbial richness during fermentation, showed the preference of oral microbial invaders for the mucosal resources, and indicated potential oral-to-intestinal mucosal competitions. It also underlined promising opportunities to further understand mechanisms of oral invasion into the human gut microbiome, define microbe-microbe and mucus-microbe interactions in a compartmentalized fashion, and help to better characterize the potential of oral microbial invasion and their persistence in the gut.

Both Disease Activity and HLA-B27 Status Are Associated With Gut Microbiome Dysbiosis in Spondyloarthritis Patients.

OBJECTIVE: Gut microbiome dysbiosis has previously been reported in spondyloarthritis (SpA) patients and could be critically involved in the pathogenesis of this disorder. The objectives of this study were to further characterize the microbiota structure in SpA patients and to investigate the relationship between dysbiosis and disease activity in light of the putative influence of the genetic background. METHODS: Shotgun sequencing was performed on fecal DNA isolated from stool samples from 2 groups of adult volunteers: SpA patients (n = 102) and healthy controls (n = 63). A subset of the healthy controls comprised the age-matched siblings of patients whose HLA-B27 status was known. Changes in gut microbiota composition were assessed based on species diversity, enterotypes, and taxonomic and functional differences. RESULTS: Dysbiosis was confirmed in SpA patients as compared to healthy controls. The restriction of microbiota diversity was detected in patients with the most active disease, and the abundance of several bacterial species was correlated with Bath Ankylosing Spondylitis Disease Activity Index score. Among healthy controls, significant differences in microbiota composition were also detected between the HLA-B27-positive and the HLA-B27-negative siblings of SpA patients. We highlighted a decreased abundance of several species of bacteria in SpA patients, especially those bacteria belonging to the Clostridiales order. Among the few species of bacteria showing increased abundance, Ruminococcus gnavus was one of the top differentiating species. CONCLUSION: These findings reveal that genetic background and level of disease activity are likely to influence the composition of the gut microbiota of patients with SpA. It may be appropriate for further research on chronic arthritis to focus on these key parameters.

Alteration of Gut Microbiome in Patients With Schizophrenia Indicates Links Between Bacterial Tyrosine Biosynthesis and Cognitive Dysfunction.

Background: Schizophrenia (SCZ) is a heterogeneous neuropsychiatric disorder for which current treatment has insufficient efficacy and severe adverse effects. The modifiable gut microbiome might be a potential target for intervention to improve neurobiological functions through the gut-microbiome-brain axis. Methods: In this case-control study, gut microbiota of 132 patients with SCZ and increased waist circumference were compared with gut microbiota of two age- and sex-matched control groups, composed of 132 healthy individuals and 132 individuals with metabolic syndrome. Shotgun sequencing was used to characterize fecal samples at the taxonomic and functional levels. Cognition of the patients with SCZ was evaluated using the Brief Assessment of Cognition instrument. Results: SCZ gut microbiota differed significantly from those of healthy control subjects and individuals with metabolic syndrome in terms of richness and global composition. SCZ gut microbiota were notably enriched in Flavonifractor plautii, Collinsella aerofaciens, Bilophila wadsworthia, and Sellimonas intestinalis, while depleted in Faecalibacterium prausnitzii, Ruminococcus lactaris, Ruminococcus bicirculans, and Veillonella rogosae. Functional potential of the gut microbiota accounted for 11% of cognition variability. In particular, the bacterial functional module for synthesizing tyrosine, a precursor for dopamine, was in SCZ cases positively associated with cognitive score (ρ = 0.34, q ≤ .1). Conclusions: Overall, this study shows that the gut microbiome of patients with SCZ differs greatly from that of healthy control subjects or individuals with metabolic syndrome. Cognitive function of patients with SCZ is associated with the potential for gut bacterial biosynthesis of tyrosine, a precursor for dopamine, suggesting that gut microbiota might be an intervention target for alleviation of cognitive dysfunction in SCZ.

Recovery of Metagenome-Assembled Genomes from a Human Fecal Sample with Pacific Biosciences High-Fidelity Sequencing.

Here, we report the recovery of 89 metagenome-assembled genomes (MAGs) derived from a human fecal sample subjected to Pacific Biosciences (PacBio) high-fidelity (HiFi) sequencing. A total of 9 MAGs consisted of complete circular contigs, and 45 MAGs were high-quality draft genomes according to the minimum information about a metagenome-assembled genome (MIMAG) standards.

Draft Genome Sequence of Isolate POC01, a Novel Anaerobic Member of the Oscillospiraceae Family, Isolated from Human Feces.

We report the isolation, culture, and genome sequencing of isolate POC01, a strictly anaerobic bacterium isolated from a healthy donor, representing a previously uncultured member of the Oscillospiraceae family.

Benchmarking second and third-generation sequencing platforms for microbial metagenomics.

Shotgun metagenomic sequencing is a common approach for studying the taxonomic diversity and metabolic potential of complex microbial communities. Current methods primarily use second generation short read sequencing, yet advances in third generation long read technologies provide opportunities to overcome some of the limitations of short read sequencing. Here, we compared seven platforms, encompassing second generation sequencers (Illumina HiSeq 300, MGI DNBSEQ-G400 and DNBSEQ-T7, ThermoFisher Ion GeneStudio S5 and Ion Proton P1) and third generation sequencers (Oxford Nanopore Technologies MinION R9 and Pacific Biosciences Sequel II). We constructed three uneven synthetic microbial communities composed of up to 87 genomic microbial strains DNAs per mock, spanning 29 bacterial and archaeal phyla, and representing the most complex and diverse synthetic communities used for sequencing technology comparisons. Our results demonstrate that third generation sequencing have advantages over second generation platforms in analyzing complex microbial communities, but require careful sequencing library preparation for optimal quantitative metagenomic analysis. Our sequencing data also provides a valuable resource for testing and benchmarking bioinformatics software for metagenomics.

Changes in Gut Microbiota of Patients with Atopic Dermatitis During Balneotherapy.

INTRODUCTION: Balneotherapy (BT) is the treatment of disease through the use of thermal spring water (TSW). It has been used for centuries and remains a popular form of treatment for dermatologic diseases such as atopic dermatitis (AD). Recent findings highlighted the role of the gut microbiota in AD and the possible crosstalk between gut and skin microbiomes in this pathology. Nevertheless, changes in the composition of the gut microbiota after balneotherapy remain to be elucidated. METHODS: A total of 96 patients, with moderate to severe AD according to the SCORAD, were enrolled. Stool samples were collected prior and post a 3-week balneotherapy at the thermal care center of La Roche-Posay (France). Composition of the gut microbiota was assessed by shotgun metagenomic sequencing. RESULTS: Species associated with high gut microbiota richness tended to correlate negatively with disease severity (SCORAD) and positively with SCORAD reduction, while species associated with low richness displayed the opposite pattern. Relative abundance of 23 species was significantly altered during BT, these changes being significantly associated with SCORAD reduction during BT, suggesting that gut microbiota composition and AD progression were associated through the treatment. Microbial functions related to gut-brain axis such as GABA and tryptophan metabolism were also altered by the treatment. Long-standing AD patients exhibited a better gut microbial profile than recently diagnosed patients, with higher MSP richness and species associated with SCORAD reduction. CONCLUSION: In patients with AD, clinical disease parameters such as SCORAD or disease duration are intricately linked to the gut microbiota composition. SCORAD reduction occurring during BT was also associated with gut microbiota. The gut-brain-skin axis via neurotransmitter such as GABA should be further studied in diseases such as AD.

Diet Supplementation with NUTRIOSE, a Resistant Dextrin, Increases the Abundance of Parabacteroides distasonis in the Human Gut.

SCOPE: An imbalance of the gut microbiota ("dysbiosis") is associated with numerous chronic diseases, and its modulation is a promising novel therapeutic approach. Dietary supplementation with soluble fiber is one of several proposed modulation strategies. This study aims at confirming the impact of the resistant dextrin NUTRIOSE (RD), a soluble fiber with demonstrated beneficial health effects, on the gut microbiota of healthy individuals. METHODS AND RESULTS: Fifty healthy women are enrolled and supplemented daily with either RD (n = 24) or a control product (n = 26) during 6 weeks. Characterization of the fecal metagenome with shotgun sequencing reveals that RD intake dramatically increases the abundance of the commensal bacterium Parabacteroides distasonis. Furthermore, presence in metagenomes of accessory genes from P. distasonis, coding for susCD (a starch-binding membrane protein complex) is associated with a greater increase of the species. This suggests that response to RD might be strain-dependent. CONCLUSION: Supplementation with RD can be used to specifically increase P. distasonis in gut microbiota of healthy women. The magnitude of the response may be associated with fiber-metabolizing capabilities of strains carried by subjects. Further research will seek to confirm that P. distasonis directly modulates the clinical effects observed in other studies.

The effects of Aronia berry (poly)phenol supplementation on arterial function and the gut microbiome in middle aged men and women: Results from a randomized controlled trial.

BACKGROUND AND AIMS: Berry (poly)phenol consumption has been associated with cardioprotective benefits, however little is known on the role the gut microbiome may play on such health benefits. Our objective was to investigate the effects of aronia berry (poly)phenol consumption on cardiometabolic health and gut microbiome richness and composition in prehypertensive middle-aged men and women. METHODS: A total of 102 prehypertensive participants were included in a parallel 12-week randomized double-blind placebo-controlled trial. Volunteers were randomly allocated to daily consume an encapsulated (poly)phenol-rich aronia berry extract (Aronia, n = 51) or a matched maltodextrin placebo (Control, n = 51). Blood pressure (BP) and arterial function (office and 24 h), endothelial function (measured as flow-mediated dilation), serum biochemistry (including blood lipids), plasma and urine (poly)phenol metabolites as well as gut microbiome composition through shotgun metagenomic sequencing were monitored over the study period. Relationships between vascular outcomes, (poly)phenol metabolites and gut microbiome were investigated using an integrated multi-levels approach. RESULTS: A significant improvement in arterial indices measured as augmentation index (AIx) and pulse wave velocity (PWV) was found in the Aronia compared to Control group (awake Δ PWV = -0.24 m/s; 95% CI: -0.79, -0.01 m/s, P < 0.05; 24 h peripheral Δ AIx = -6.8; -11.2, -2.3, %, P = 0.003; 24 h central Δ AIx = -3.3; -5.5, -1.0, %, P = 0.006). No changes in BP, endothelial function or blood lipids were found following the intervention. Consumption of aronia (poly)phenols led to a significant increase in gut microbiome gene richness and in the abundance of butyrate-producing species such as Lawsonibacter asaccharolyticus and Intestinimonas butyriciproducens species, compared to Control group. Results from an approach including metabolomic, metagenomic and clinical outcomes highlighted associations between aronia-derived phenolic metabolites, arterial stiffness, and gut microbiome. CONCLUSIONS: Aronia berry (poly)phenol consumption improved arterial function in prehypertensive middle-aged individuals, possibly via modulation of gut microbiome richness and composition based on the associations observed between these parameters. CLINICAL TRIAL REGISTRY: The National Institutes of Health (NIH)-randomized trial records held on the NIH ClinicalTrials.gov website (NCT03434574). Aronia Berry Consumption on Blood Pressure.

StrainFLAIR: strain-level profiling of metagenomic samples using variation graphs.

Current studies are shifting from the use of single linear references to representation of multiple genomes organised in pangenome graphs or variation graphs. Meanwhile, in metagenomic samples, resolving strain-level abundances is a major step in microbiome studies, as associations between strain variants and phenotype are of great interest for diagnostic and therapeutic purposes. We developed StrainFLAIR with the aim of showing the feasibility of using variation graphs for indexing highly similar genomic sequences up to the strain level, and for characterizing a set of unknown sequenced genomes by querying this graph. On simulated data composed of mixtures of strains from the same bacterial species Escherichia coli, results show that StrainFLAIR was able to distinguish and estimate the abundances of close strains, as well as to highlight the presence of a new strain close to a referenced one and to estimate its abundance. On a real dataset composed of a mix of several bacterial species and several strains for the same species, results show that in a more complex configuration StrainFLAIR correctly estimates the abundance of each strain. Hence, results demonstrated how graph representation of multiple close genomes can be used as a reference to characterize a sample at the strain level.

High engraftment capacity of frozen ready-to-use human fecal microbiota transplants assessed in germ-free mice.

The number of indications for fecal microbiota transplantation is expected to rise, thus increasing the needs for production of readily available frozen or freeze-dried transplants. Using shotgun metagenomics, we investigated the capacity of two novel human fecal microbiota transplants prepared in maltodextrin-trehalose solutions (abbreviated MD and TR for maltodextrin:trehalose, 3:1, w/w, and trehalose:maltodextrin 3:1, w/w, respectively), to colonize a germ-free born mouse model. Gavage with frozen-thawed MD or TR suspensions gave the taxonomic profiles of mouse feces that best resembled those obtained with the fresh inoculum (Spearman correlations based on relative abundances of metagenomic species around 0.80 and 0.75 for MD and TR respectively), while engraftment capacity of defrosted NaCl transplants most diverged (Spearman correlations around 0.63). Engraftment of members of the family Lachnospiraceae and Ruminoccocaceae was the most challenging in all groups of mice, being improved with MD and TR transplants compared to NaCl, but still lower than with the fresh preparation. Improvement of engraftment of this important group in maintaining health represents a challenge that could benefit from further research on fecal microbiota transplant manufacturing.