High-quality metagenome assembly from long accurate reads with metaMDBG.

We introduce metaMDBG, a metagenomics assembler for PacBio HiFi reads. MetaMDBG combines a de Bruijn graph assembly in a minimizer space with an iterative assembly over sequences of minimizers to address variations in genome coverage depth and an abundance-based filtering strategy to simplify strain complexity. For complex communities, we obtained up to twice as many high-quality circularized prokaryotic metagenome-assembled genomes as existing methods and had better recovery of viruses and plasmids.

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.

Efficient High-Quality Metagenome Assembly from Long Accurate Reads using Minimizer-space de Bruijn Graphs.

We introduce a novel metagenomics assembler for high-accuracy long reads. Our approach, implemented as metaMDBG, combines highly efficient de Bruijn graph assembly in minimizer space, with both a multi-k' approach for dealing with variations in genome coverage depth and an abundance-based filtering strategy for simplifying strain complexity. The resulting algorithm is more efficient than the state-of-the-art but with better assembly results. metaMDBG was 1.5 to 12 times faster than competing assemblers and requires between one-tenth and one-thirtieth of the memory across a range of data sets. We obtained up to twice as many high-quality circularised prokaryotic metagenome assembled genomes (MAGs) on the most complex communities, and a better recovery of viruses and plasmids. metaMDBG performs particularly well for abundant organisms whilst being robust to the presence of strain diversity. The result is that for the first time it is possible to efficiently reconstruct the majority of complex communities by abundance as near-complete MAGs.

Population-level impacts of antibiotic usage on the human gut microbiome.

The widespread usage of antimicrobials has driven the evolution of resistance in pathogenic microbes, both increased prevalence of antimicrobial resistance genes (ARGs) and their spread across species by horizontal gene transfer (HGT). However, the impact on the wider community of commensal microbes associated with the human body, the microbiome, is less well understood. Small-scale studies have determined the transient impacts of antibiotic consumption but we conduct an extensive survey of ARGs in 8972 metagenomes to determine the population-level impacts. Focusing on 3096 gut microbiomes from healthy individuals not taking antibiotics we demonstrate highly significant correlations between both the total ARG abundance and diversity and per capita antibiotic usage rates across ten countries spanning three continents. Samples from China were notable outliers. We use a collection of 154,723 human-associated metagenome assembled genomes (MAGs) to link these ARGs to taxa and detect HGT. This reveals that the correlations in ARG abundance are driven by multi-species mobile ARGs shared between pathogens and commensals, within a highly connected central component of the network of MAGs and ARGs. We also observe that individual human gut ARG profiles cluster into two types or resistotypes. The less frequent resistotype has higher overall ARG abundance, is associated with certain classes of resistance, and is linked to species-specific genes in the Proteobacteria on the periphery of the ARG network.

Crop management shapes the diversity and activity of DNA and RNA viruses in the rhizosphere.

BACKGROUND: The rhizosphere is a hotspot for microbial activity and contributes to ecosystem services including plant health and biogeochemical cycling. The activity of microbial viruses, and their influence on plant-microbe interactions in the rhizosphere, remains undetermined. Given the impact of viruses on the ecology and evolution of their host communities, determining how soil viruses influence microbiome dynamics is crucial to build a holistic understanding of rhizosphere functions. RESULTS: Here, we aimed to investigate the influence of crop management on the composition and activity of bulk soil, rhizosphere soil, and root viral communities. We combined viromics, metagenomics, and metatranscriptomics on soil samples collected from a 3-year crop rotation field trial of oilseed rape (Brassica napus L.). By recovering 1059 dsDNA viral populations and 16,541 ssRNA bacteriophage populations, we expanded the number of underexplored Leviviricetes genomes by > 5 times. Through detection of viral activity in metatranscriptomes, we uncovered evidence of "Kill-the-Winner" dynamics, implicating soil bacteriophages in driving bacterial community succession. Moreover, we found the activity of viruses increased with proximity to crop roots, and identified that soil viruses may influence plant-microbe interactions through the reprogramming of bacterial host metabolism. We have provided the first evidence of crop rotation-driven impacts on soil microbial communities extending to viruses. To this aim, we present the novel principal of "viral priming," which describes how the consecutive growth of the same crop species primes viral activity in the rhizosphere through local adaptation. CONCLUSIONS: Overall, we reveal unprecedented spatial and temporal diversity in viral community composition and activity across root, rhizosphere soil, and bulk soil compartments. Our work demonstrates that the roles of soil viruses need greater consideration to exploit the rhizosphere microbiome for food security, food safety, and environmental sustainability. Video Abstract.

Stimulation of Distinct Rhizosphere Bacteria Drives Phosphorus and Nitrogen Mineralization in Oilseed Rape under Field Conditions.

Advances in DNA sequencing technologies have drastically changed our perception of the structure and complexity of the plant microbiome. By comparison, our ability to accurately identify the metabolically active fraction of soil microbiota and its specific functional role in augmenting plant health is relatively limited. Important ecological interactions being performed by microbes can be investigated by analyzing the extracellular protein fraction. Here, we combined a unique protein extraction method and an iterative bioinformatics pipeline to capture and identify extracellular proteins (metaexoproteomics) synthesized in the rhizosphere of Brassica spp. We first validated our method in the laboratory by successfully identifying proteins related to a host plant (Brassica rapa) and its bacterial inoculant, Pseudomonas putida BIRD-1. This identified numerous rhizosphere specific proteins linked to the acquisition of plant-derived nutrients in P. putida. Next, we analyzed natural field-soil microbial communities associated with Brassica napus L. (oilseed rape). By combining metagenomics with metaexoproteomics, 1,885 plant, insect, and microbial proteins were identified across bulk and rhizosphere samples. Metaexoproteomics identified a significant shift in the metabolically active fraction of the soil microbiota responding to the presence of B. napus roots that was not apparent in the composition of the total microbial community (metagenome). This included stimulation of rhizosphere-specialized bacteria, such as Gammaproteobacteria, Betaproteobacteria, and Flavobacteriia, and the upregulation of plant beneficial functions related to phosphorus and nitrogen mineralization. Our metaproteomic assessment of the "active" plant microbiome at the field-scale demonstrates the importance of moving beyond metagenomics to determine ecologically important plant-microbe interactions underpinning plant health. IMPORTANCE Plant-microbe interactions are critical to ecosystem function and crop production. While significant advances have been made toward understanding the structure of the plant microbiome, learning about its full functional role is still in its infancy. This is primarily due to an incomplete ability to determine in situ plant-microbe interactions actively operating under field conditions. Proteins are the functional entities of the cell. Therefore, their identification and relative quantification within a microbial community provide the best proxy for which microbes are the most metabolically active and which are driving important plant-microbe interactions. Here, we provide the first metaexoproteomics assessment of the plant microbiome using field-grown oilseed rape as the model crop species, identifying key taxa responsible for specific ecological interactions. Gaining a mechanistic understanding of the plant microbiome is central to developing engineered plant microbiomes to improve sustainable agricultural approaches and reduce our reliance on nonrenewable resources.

Impact of sulfamethoxazole on a riverine microbiome.

The continued emergence of bacterial pathogens presenting antimicrobial resistance is widely recognised as a global health threat and recent attention focused on potential environmental reservoirs of antibiotic resistance genes (ARGs). Freshwater environments such as rivers represent a potential hotspot for ARGs and antibiotic resistant bacteria as they are receiving systems for effluent discharges from wastewater treatment plants (WWTPs). Effluent also contains low levels of different antimicrobials including antibiotics and biocides. Sulfonamides are antibacterial chemicals widely used in clinical, veterinary and agricultural settings and are frequently detected in sewage sludge and manure in addition to riverine ecosystems. The impact of such exposure on ARG prevalence and diversity is unknown, so the aim of this study was to investigate the release of a sub-lethal concentration of the sulfonamide compound sulfamethoxazole (SMX) on the river bacterial microbiome using a flume system. This system was a semi-natural in vitro flume using river water (30 L) and sediment (6 kg) with circulation to mimic river flow. A combination of 'omics' approaches were conducted to study the impact of SMX exposure on the microbiomes within the flumes. Metagenomic analysis showed that the addition of low concentrations of SMX (<4 µg L-1) had a limited effect on the bacterial resistome in the water fraction only, with no impact observed in the sediment. Metaproteomics did not show differences in ARGs expression with SMX exposure in water. Overall, the river bacterial community was resilient to short term exposure to sub-lethal concentrations of SMX which mimics the exposure such communities experience downstream of WWTPs throughout the year.

STRONG: metagenomics strain resolution on assembly graphs.

We introduce STrain Resolution ON assembly Graphs (STRONG), which identifies strains de novo, from multiple metagenome samples. STRONG performs coassembly, and binning into metagenome assembled genomes (MAGs), and stores the coassembly graph prior to variant simplification. This enables the subgraphs and their unitig per-sample coverages, for individual single-copy core genes (SCGs) in each MAG, to be extracted. A Bayesian algorithm, BayesPaths, determines the number of strains present, their haplotypes or sequences on the SCGs, and abundances. STRONG is validated using synthetic communities and for a real anaerobic digestor time series generates haplotypes that match those observed from long Nanopore reads.

Gene duplication drives genome expansion in a major lineage of Thaumarchaeota.

Ammonia-oxidising archaea of the phylum Thaumarchaeota are important organisms in the nitrogen cycle, but the mechanisms driving their radiation into diverse ecosystems remain underexplored. Here, existing thaumarchaeotal genomes are complemented with 12 genomes belonging to the previously under-sampled Nitrososphaerales to investigate the impact of lateral gene transfer (LGT), gene duplication and loss across thaumarchaeotal evolution. We reveal a major role for gene duplication in driving genome expansion subsequent to early LGT. In particular, two large LGT events are identified into Nitrososphaerales and the fate of these gene families is highly lineage-specific, being lost in some descendant lineages, but undergoing extensive duplication in others, suggesting niche-specific roles. Notably, some genes involved in carbohydrate transport or coenzyme metabolism were duplicated, likely facilitating niche specialisation in soils and sediments. Overall, our results suggest that LGT followed by gene duplication drives Nitrososphaerales evolution, highlighting a previously under-appreciated mechanism of genome expansion in archaea.

DESMAN: a new tool for de novo extraction of strains from metagenomes.

We introduce DESMAN for De novo Extraction of Strains from Metagenomes. Large multi-sample metagenomes are being generated but strain variation results in fragmentary co-assemblies. Current algorithms can bin contigs into metagenome-assembled genomes but are unable to resolve strain-level variation. DESMAN identifies variants in core genes and uses co-occurrence across samples to link variants into haplotypes and abundance profiles. These are then searched for against non-core genes to determine the accessory genome of each strain. We validated DESMAN on a complex 50-species 210-genome 96-sample synthetic mock data set and then applied it to the Tara Oceans microbiome.

Inferring Aggregated Functional Traits from Metagenomic Data Using Constrained Non-negative Matrix Factorization: Application to Fiber Degradation in the Human Gut Microbiota.

Whole Genome Shotgun (WGS) metagenomics is increasingly used to study the structure and functions of complex microbial ecosystems, both from the taxonomic and functional point of view. Gene inventories of otherwise uncultured microbial communities make the direct functional profiling of microbial communities possible. The concept of community aggregated trait has been adapted from environmental and plant functional ecology to the framework of microbial ecology. Community aggregated traits are quantified from WGS data by computing the abundance of relevant marker genes. They can be used to study key processes at the ecosystem level and correlate environmental factors and ecosystem functions. In this paper we propose a novel model based approach to infer combinations of aggregated traits characterizing specific ecosystemic metabolic processes. We formulate a model of these Combined Aggregated Functional Traits (CAFTs) accounting for a hierarchical structure of genes, which are associated on microbial genomes, further linked at the ecosystem level by complex co-occurrences or interactions. The model is completed with constraints specifically designed to exploit available genomic information, in order to favor biologically relevant CAFTs. The CAFTs structure, as well as their intensity in the ecosystem, is obtained by solving a constrained Non-negative Matrix Factorization (NMF) problem. We developed a multicriteria selection procedure for the number of CAFTs. We illustrated our method on the modelling of ecosystemic functional traits of fiber degradation by the human gut microbiota. We used 1408 samples of gene abundances from several high-throughput sequencing projects and found that four CAFTs only were needed to represent the fiber degradation potential. This data reduction highlighted biologically consistent functional patterns while providing a high quality preservation of the original data. Our method is generic and can be applied to other metabolic processes in the gut or in other ecosystems.

Identification of metabolic signatures linked to anti-inflammatory effects of Faecalibacterium prausnitzii.

UNLABELLED: Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified on the basis of human clinical data. The mechanisms underlying its beneficial effects are still unknown. Gnotobiotic mice harboring F. prausnitzii (A2-165) and Escherichia coli (K-12 JM105) were subjected to 2,4,6-trinitrobenzenesulfonic acid (TNBS)-induced acute colitis. The inflammatory colitis scores and a gas chromatography-time of flight (GC/TOF) mass spectrometry-based metabolomic profile were monitored in blood, ileum, cecum, colon, and feces in gnotobiotic mice. The potential anti-inflammatory metabolites were tested in vitro. We obtained stable E. coli and F. prausnitzii-diassociated mice in which E. coli primed the gastrointestinal tract (GIT), allowing a durable and stable establishment of F. prausnitzii. The disease activity index, histological scores, myeloperoxidase (MPO) activity, and serum cytokine levels were significantly lower in the presence of F. prausnitzii after TNBS challenge. The protective effect of F. prausnitzii against colitis was correlated to its implantation level and was linked to overrepresented metabolites along the GIT and in serum. Among 983 metabolites in GIT samples and serum, 279 were assigned to known chemical reactions. Some of them, belonging to the ammonia (α-ketoglutarate), osmoprotective (raffinose), and phenolic (including anti-inflammatory shikimic and salicylic acids) pathways, were associated with a protective effect of F. prausnitzii, and the functional link was established in vitro for salicylic acid. We show for the first time that F. prausnitzii is a highly active commensal bacterium involved in reduction of colitis through in vivo modulation of metabolites along the GIT and in the peripheral blood. IMPORTANCE: Inflammatory bowel diseases (IBD) are characterized by low proportions of F. prausnitzii in the gut microbiome. This commensal bacterium exhibits anti-inflammatory effects through still unknown mechanisms. Stable monoassociated rodents are actually not a reproducible model to decipher F. prausnitzii protective effects. We propose a new gnotobiotic rodent model providing mechanistic clues. In this model, F. prausnitzii exhibits protective effects against an acute colitis and a protective metabolic profile is linked to its presence along the digestive tract. We identified a molecule, salicylic acid, directly involved in the protective effect of F. prausnitzii. Targeting its metabolic pathways could be an attractive therapeutic strategy in IBD.