Multicenter evaluation of gut microbiome profiling by next-generation sequencing reveals major biases in partial-length metabarcoding approach.

Next-generation sequencing workflows, using either metabarcoding or metagenomic approaches, have massively contributed to expanding knowledge of the human gut microbiota, but methodological bias compromises reproducibility across studies. Where these biases have been quantified within several comparative analyses on their own, none have measured inter-laboratory reproducibility using similar DNA material. Here, we designed a multicenter study involving seven participating laboratories dedicated to partial- (P1 to P5), full-length (P6) metabarcoding, or metagenomic profiling (MGP) using DNA from a mock microbial community or extracted from 10 fecal samples collected at two time points from five donors. Fecal material was collected, and the DNA was extracted according to the IHMS protocols. The mock and isolated DNA were then provided to the participating laboratories for sequencing. Following sequencing analysis according to the laboratories' routine pipelines, relative taxonomic-count tables defined at the genus level were provided and analyzed. Large variations in alpha-diversity between laboratories, uncorrelated with sequencing depth, were detected among the profiles. Half of the genera identified by P1 were unique to this partner and two-thirds of the genera identified by MGP were not detected by P3. Analysis of beta-diversity revealed lower inter-individual variance than inter-laboratory variances. The taxonomic profiles of P5 and P6 were more similar to those of MGP than those obtained by P1, P2, P3, and P4. Reanalysis of the raw sequences obtained by partial-length metabarcoding profiling, using a single bioinformatic pipeline, harmonized the description of the bacterial profiles, which were more similar to each other, except for P3, and closer to the profiles obtained by MGP. This study highlights the major impact of the bioinformatics pipeline, and primarily the database used for taxonomic annotation. Laboratories need to benchmark and optimize their bioinformatic pipelines using standards to monitor their effectiveness in accurately detecting taxa present in gut microbiota.

The gut microbiota contributes to the pathogenesis of anorexia nervosa in humans and mice.

Anorexia nervosa (AN) is an eating disorder with a high mortality. About 95% of cases are women and it has a population prevalence of about 1%, but evidence-based treatment is lacking. The pathogenesis of AN probably involves genetics and various environmental factors, and an altered gut microbiota has been observed in individuals with AN using amplicon sequencing and relatively small cohorts. Here we investigated whether a disrupted gut microbiota contributes to AN pathogenesis. Shotgun metagenomics and metabolomics were performed on faecal and serum samples, respectively, from a cohort of 77 females with AN and 70 healthy females. Multiple bacterial taxa (for example, Clostridium species) were altered in AN and correlated with estimates of eating behaviour and mental health. The gut virome was also altered in AN including a reduction in viral-bacterial interactions. Bacterial functional modules associated with the degradation of neurotransmitters were enriched in AN and various structural variants in bacteria were linked to metabolic features of AN. Serum metabolomics revealed an increase in metabolites associated with reduced food intake (for example, indole-3-propionic acid). Causal inference analyses implied that serum bacterial metabolites are potentially mediating the impact of an altered gut microbiota on AN behaviour. Further, we performed faecal microbiota transplantation from AN cases to germ-free mice under energy-restricted feeding to mirror AN eating behaviour. We found that the reduced weight gain and induced hypothalamic and adipose tissue gene expression were related to aberrant energy metabolism and eating behaviour. Our 'omics' and mechanistic studies imply that a disruptive gut microbiome may contribute to AN pathogenesis.

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.

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.

Increasing the diversity of dietary fibers in a daily-consumed bread modifies gut microbiota and metabolic profile in subjects at cardiometabolic risk.

Some cardiometabolic risk factors such as dyslipidemia and insulin resistance are known to be associated with low gut microbiota richness. A link between gut microbiota richness and the diversity of consumed dietary fibers (DF) has also been reported. We introduced a larger diversity of consumed DF by using a daily consumed bread in subjects at cardiometabolic risk and assessed the impacts on the composition and functions of gut microbiota as well as on cardiometabolic profile. Thirty-nine subjects at cardiometabolic risk were included in a double-blind, randomized, cross-over, twice 8-week study, and consumed daily 150 g of standard bread or enriched with a 7-dietary fiber mixture (5.55 g and 16.05 g of fibers, respectively). Before and after intervention, stool samples were collected for gut microbiota analysis from species determination down to gene-level abundance using shotgun metagenomics, and cardiometabolic profile was assessed. Multi-fiber bread consumption significantly decreased Bacteroides vulgatus, whereas it increased Parabacteroides distasonis, Fusicatenibacter saccharivorans, an unclassified Acutalibacteraceae and an unclassified Eisenbergiella (q < 0.1). The fraction of gut microbiota carrying the gene coding for five families/subfamilies of glycoside hydrolases (CAZymes) were also increased and negatively correlated with peaks and total/incremental area under curve (tAUC/iAUC) of postprandial glycemia and insulinemia. Compared to control bread, multi-fiber bread decreased total cholesterol (-0.42 mM; q < 0.01), LDL cholesterol (-0.36 mM; q < 0.01), insulin (-2.77 mIU/l; q < 0.05), and HOMA (-0.78; q < 0.05). In conclusion, increasing the diversity of DF in a daily consumed product modifies gut microbiota composition and function and could be a relevant nutritional tool to improve cardiometabolic profile.

Long-term diosmectite use does not alter the gut microbiota in adults with chronic diarrhea.

BACKGROUND: Diosmectite, a natural colloidal clay, has been used worldwide for a number of approved indications, including the treatment of chronic functional diarrhea. Here, we used high-resolution whole metagenome shotgun sequencing to assess the impact of a 5 weeks administration of diosmectite (3 g/sachet, 3 sachets/day) on the fecal microbiota of 35 adults with functional chronic diarrhea. RESULTS: Gut microbiota was not impacted by diosmectite administration. In particular, richness remained stable and no microbial species displayed a significant evolution. Segregating patients either by diosmectite response (non responder, early responder, late responder) or by nationality (Great-Britain or Netherlands) yielded the same results. CONCLUSION: We concluded that no microbiota-related physiological alterations are expected upon long-term treatment with diosmectite. TRIAL REGISTRATION: Clinicaltrials.gov NCT03045926.

Microbiome and metabolome features of the cardiometabolic disease spectrum.

Previous microbiome and metabolome analyses exploring non-communicable diseases have paid scant attention to major confounders of study outcomes, such as common, pre-morbid and co-morbid conditions, or polypharmacy. Here, in the context of ischemic heart disease (IHD), we used a study design that recapitulates disease initiation, escalation and response to treatment over time, mirroring a longitudinal study that would otherwise be difficult to perform given the protracted nature of IHD pathogenesis. We recruited 1,241 middle-aged Europeans, including healthy individuals, individuals with dysmetabolic morbidities (obesity and type 2 diabetes) but lacking overt IHD diagnosis and individuals with IHD at three distinct clinical stages-acute coronary syndrome, chronic IHD and IHD with heart failure-and characterized their phenome, gut metagenome and serum and urine metabolome. We found that about 75% of microbiome and metabolome features that distinguish individuals with IHD from healthy individuals after adjustment for effects of medication and lifestyle are present in individuals exhibiting dysmetabolism, suggesting that major alterations of the gut microbiome and metabolome might begin long before clinical onset of IHD. We further categorized microbiome and metabolome signatures related to prodromal dysmetabolism, specific to IHD in general or to each of its three subtypes or related to escalation or de-escalation of IHD. Discriminant analysis based on specific IHD microbiome and metabolome features could better differentiate individuals with IHD from healthy individuals or metabolically matched individuals as compared to the conventional risk markers, pointing to a pathophysiological relevance of these features.

Intestinal Akkermansia muciniphila predicts clinical response to PD-1 blockade in patients with advanced non-small-cell lung cancer.

Aside from PD-L1 expression, biomarkers of response to immune checkpoint inhibitors (ICIs) in non-small-cell lung cancer (NSCLC) are needed. In a previous retrospective analysis, we documented that fecal Akkermansia muciniphila (Akk) was associated with clinical benefit of ICI in patients with NSCLC or kidney cancer. In the current study, we performed shotgun-metagenomics-based microbiome profiling in a large cohort of patients with advanced NSCLC (n = 338) treated with first- or second-line ICIs to prospectively validate the predictive value of fecal Akk. Baseline stool Akk was associated with increased objective response rates and overall survival in multivariate analyses, independent of PD-L1 expression, antibiotics, and performance status. Intestinal Akk was accompanied by a richer commensalism, including Eubacterium hallii and Bifidobacterium adolescentis, and a more inflamed tumor microenvironment in a subset of patients. However, antibiotic use (20% of cases) coincided with a relative dominance of Akk above 4.8% accompanied with the genus Clostridium, both associated with resistance to ICI. Our study shows significant differences in relative abundance of Akk that may represent potential biomarkers to refine patient stratification in future studies.

Impairment of gut microbial biotin metabolism and host biotin status in severe obesity: effect of biotin and prebiotic supplementation on improved metabolism.

OBJECTIVES: Gut microbiota is a key component in obesity and type 2 diabetes, yet mechanisms and metabolites central to this interaction remain unclear. We examined the human gut microbiome's functional composition in healthy metabolic state and the most severe states of obesity and type 2 diabetes within the MetaCardis cohort. We focused on the role of B vitamins and B7/B8 biotin for regulation of host metabolic state, as these vitamins influence both microbial function and host metabolism and inflammation. DESIGN: We performed metagenomic analyses in 1545 subjects from the MetaCardis cohorts and different murine experiments, including germ-free and antibiotic treated animals, faecal microbiota transfer, bariatric surgery and supplementation with biotin and prebiotics in mice. RESULTS: Severe obesity is associated with an absolute deficiency in bacterial biotin producers and transporters, whose abundances correlate with host metabolic and inflammatory phenotypes. We found suboptimal circulating biotin levels in severe obesity and altered expression of biotin-associated genes in human adipose tissue. In mice, the absence or depletion of gut microbiota by antibiotics confirmed the microbial contribution to host biotin levels. Bariatric surgery, which improves metabolism and inflammation, associates with increased bacterial biotin producers and improved host systemic biotin in humans and mice. Finally, supplementing high-fat diet-fed mice with fructo-oligosaccharides and biotin improves not only the microbiome diversity, but also the potential of bacterial production of biotin and B vitamins, while limiting weight gain and glycaemic deterioration. CONCLUSION: Strategies combining biotin and prebiotic supplementation could help prevent the deterioration of metabolic states in severe obesity. TRIAL REGISTRATION NUMBER: NCT02059538.

Combinatorial, additive and dose-dependent drug-microbiome associations.

During the transition from a healthy state to cardiometabolic disease, patients become heavily medicated, which leads to an increasingly aberrant gut microbiome and serum metabolome, and complicates biomarker discovery1-5. Here, through integrated multi-omics analyses of 2,173 European residents from the MetaCardis cohort, we show that the explanatory power of drugs for the variability in both host and gut microbiome features exceeds that of disease. We quantify inferred effects of single medications, their combinations as well as additive effects, and show that the latter shift the metabolome and microbiome towards a healthier state, exemplified in synergistic reduction in serum atherogenic lipoproteins by statins combined with aspirin, or enrichment of intestinal Roseburia by diuretic agents combined with beta-blockers. Several antibiotics exhibit a quantitative relationship between the number of courses prescribed and progression towards a microbiome state that is associated with the severity of cardiometabolic disease. We also report a relationship between cardiometabolic drug dosage, improvement in clinical markers and microbiome composition, supporting direct drug effects. Taken together, our computational framework and resulting resources enable the disentanglement of the effects of drugs and disease on host and microbiome features in multimedicated individuals. Furthermore, the robust signatures identified using our framework provide new hypotheses for drug-host-microbiome interactions in cardiometabolic disease.

Imidazole propionate is increased in diabetes and associated with dietary patterns and altered microbial ecology.

Microbiota-host-diet interactions contribute to the development of metabolic diseases. Imidazole propionate is a novel microbially produced metabolite from histidine, which impairs glucose metabolism. Here, we show that subjects with prediabetes and diabetes in the MetaCardis cohort from three European countries have elevated serum imidazole propionate levels. Furthermore, imidazole propionate levels were increased in subjects with low bacterial gene richness and Bacteroides 2 enterotype, which have previously been associated with obesity. The Bacteroides 2 enterotype was also associated with increased abundance of the genes involved in imidazole propionate biosynthesis from dietary histidine. Since patients and controls did not differ in their histidine dietary intake, the elevated levels of imidazole propionate in type 2 diabetes likely reflects altered microbial metabolism of histidine, rather than histidine intake per se. Thus the microbiota may contribute to type 2 diabetes by generating imidazole propionate that can modulate host inflammation and metabolism.

Long inverted repeats around the chromosome replication terminus in the model strain Bacillus thuringiensis serovar israelensis BGSC 4Q7.

Bacillus thuringiensis serovar israelensis is the most widely used natural biopesticide against mosquito larvae worldwide. Its lineage has been actively studied and a plasmid-free strain, B. thuringiensis serovar israelensis BGSC 4Q7 (4Q7), has been produced. Previous sequencing of the genome of this strain has revealed the persistent presence of a 235 kb extrachromosomal element, pBtic235, which has been shown to be an inducible prophage, although three putative chromosomal prophages have been lost. Moreover, a 492 kb region, potentially including the standard replication terminus, has also been deleted in the 4Q7 strain, indicating an absence of essential genes in this area. We reanalysed the genome coverage distribution of reads for the previously sequenced variant strain, and sequenced two independently maintained samples of the 4Q7 strain. A 553 kb area, close to the 492 kb deletion, was found to be duplicated. This duplication presumably restored the equal sizes of the replichores, and a balanced functioning of replication termination. An analysis of genome assembly graphs revealed a transient association of the host chromosome with the pBtic235 element. This association may play a functional role in the replication of the bacterial chromosome, and the termination of this process in particular. The genome-restructuring events detected may modify the genetic status of cytotoxic or haemolytic toxins, potentially influencing strain virulence. Twelve of the single-nucleotide variants identified in 4Q7 were probably due to the procedure used for strain construction or were present in the precursor of this strain. No sequence variants were found in pBtic235, but the distribution of the corresponding 4Q7 reads indicates a significant difference from counterparts in natural B. thuringiensis serovar israelensis strains, suggesting a duplication or over-replication in 4Q7. Thus, the 4Q7 strain is not a pure plasmid-less offshoot, but a highly genetically modified derivative of its natural ancestor. In addition to potentially influencing virulence, genome-restructuring events can modify the replication termination machinery. These findings have potential implications for the conclusions of virulence studies on 4Q7 as a model, but they also raise interesting fundamental questions about the functioning of the Bacillus genome.

Integration of time-series meta-omics data reveals how microbial ecosystems respond to disturbance.

The development of reliable, mixed-culture biotechnological processes hinges on understanding how microbial ecosystems respond to disturbances. Here we reveal extensive phenotypic plasticity and niche complementarity in oleaginous microbial populations from a biological wastewater treatment plant. We perform meta-omics analyses (metagenomics, metatranscriptomics, metaproteomics and metabolomics) on in situ samples over 14 months at weekly intervals. Based on 1,364 de novo metagenome-assembled genomes, we uncover four distinct fundamental niche types. Throughout the time-series, we observe a major, transient shift in community structure, coinciding with substrate availability changes. Functional omics data reveals extensive variation in gene expression and substrate usage amongst community members. Ex situ bioreactor experiments confirm that responses occur within five hours of a pulse disturbance, demonstrating rapid adaptation by specific populations. Our results show that community resistance and resilience are a function of phenotypic plasticity and niche complementarity, and set the foundation for future ecological engineering efforts.

Correction to: Microbiome definition re-visited: old concepts and new challenges.

An amendment to this paper has been published and can be accessed via the original article.

Microbiome definition re-visited: old concepts and new challenges.

The field of microbiome research has evolved rapidly over the past few decades and has become a topic of great scientific and public interest. As a result of this rapid growth in interest covering different fields, we are lacking a clear commonly agreed definition of the term "microbiome." Moreover, a consensus on best practices in microbiome research is missing. Recently, a panel of international experts discussed the current gaps in the frame of the European-funded MicrobiomeSupport project. The meeting brought together about 40 leaders from diverse microbiome areas, while more than a hundred experts from all over the world took part in an online survey accompanying the workshop. This article excerpts the outcomes of the workshop and the corresponding online survey embedded in a short historical introduction and future outlook. We propose a definition of microbiome based on the compact, clear, and comprehensive description of the term provided by Whipps et al. in 1988, amended with a set of novel recommendations considering the latest technological developments and research findings. We clearly separate the terms microbiome and microbiota and provide a comprehensive discussion considering the composition of microbiota, the heterogeneity and dynamics of microbiomes in time and space, the stability and resilience of microbial networks, the definition of core microbiomes, and functionally relevant keystone species as well as co-evolutionary principles of microbe-host and inter-species interactions within the microbiome. These broad definitions together with the suggested unifying concepts will help to improve standardization of microbiome studies in the future, and could be the starting point for an integrated assessment of data resulting in a more rapid transfer of knowledge from basic science into practice. Furthermore, microbiome standards are important for solving new challenges associated with anthropogenic-driven changes in the field of planetary health, for which the understanding of microbiomes might play a key role. Video Abstract.

Statin therapy is associated with lower prevalence of gut microbiota dysbiosis.

Microbiome community typing analyses have recently identified the Bacteroides2 (Bact2) enterotype, an intestinal microbiota configuration that is associated with systemic inflammation and has a high prevalence in loose stools in humans1,2. Bact2 is characterized by a high proportion of Bacteroides, a low proportion of Faecalibacterium and low microbial cell densities1,2, and its prevalence varies from 13% in a general population cohort to as high as 78% in patients with inflammatory bowel disease2. Reported changes in stool consistency3 and inflammation status4 during the progression towards obesity and metabolic comorbidities led us to propose that these developments might similarly correlate with an increased prevalence of the potentially dysbiotic Bact2 enterotype. Here, by exploring obesity-associated microbiota alterations in the quantitative faecal metagenomes of the cross-sectional MetaCardis Body Mass Index Spectrum cohort (n = 888), we identify statin therapy as a key covariate of microbiome diversification. By focusing on a subcohort of participants that are not medicated with statins, we find that the prevalence of Bact2 correlates with body mass index, increasing from 3.90% in lean or overweight participants to 17.73% in obese participants. Systemic inflammation levels in Bact2-enterotyped individuals are higher than predicted on the basis of their obesity status, indicative of Bact2 as a dysbiotic microbiome constellation. We also observe that obesity-associated microbiota dysbiosis is negatively associated with statin treatment, resulting in a lower Bact2 prevalence of 5.88% in statin-medicated obese participants. This finding is validated in both the accompanying MetaCardis cardiovascular disease dataset (n = 282) and the independent Flemish Gut Flora Project population cohort (n = 2,345). The potential benefits of statins in this context will require further evaluation in a prospective clinical trial to ascertain whether the effect is reproducible in a randomized population and before considering their application as microbiota-modulating therapeutics.

Gut Bacteria Composition Drives Primary Resistance to Cancer Immunotherapy in Renal Cell Carcinoma Patients.

BACKGROUND: The development of immune checkpoint blockade (ICB) has revolutionized the clinical outcome of renal cell carcinoma (RCC). Nevertheless, improvement of durability and prediction of responses remain unmet medical needs. While it has been recognized that antibiotics (ATBs) decrease the clinical activity of ICB across various malignancies, little is known about the direct impact of distinct intestinal nonpathogenic bacteria (commensals) on therapeutic outcomes of ICB in RCC. OBJECTIVE: To evaluate the predictive value of stool bacteria composition for ICB efficacy in a cohort of advanced RCC patients. DESIGN, SETTING, AND PARTICIPANTS: We prospectively collected fecal samples from 69 advanced RCC patients treated with nivolumab and enrolled in the GETUG-AFU 26 NIVOREN microbiota translational substudy phase 2 trial (NCT03013335) at Gustave Roussy. We recorded patient characteristics including ATB use, prior systemic therapies, and response criteria. We analyzed 2994 samples of feces from healthy volunteers (HVs). In parallel, preclinical studies performed in RCC-bearing mice that received fecal transplant (FMT) from RCC patients resistant to ICB (NR-FMT) allowed us to draw a cause-effect relationship between gut bacteria composition and clinical outcomes for ICB. The influence of tyrosine kinase inhibitors (TKIs) taken before starting nivolumab on the microbiota composition has also been assessed. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Metagenomic data (MG) from whole genome sequencing (WGS) were analyzed by multivariate and pairwise comparisons/fold ratio to identify bacterial fingerprints related to ATB or prior TKI exposure and patients' therapeutic response (overall response and progression-free survival), and compared with the data from cancer-free donors. RESULTS AND LIMITATIONS: Recent ATB use (n = 11; 16%) reduced objective response rates (from 28% to 9%, p < 0.03) and markedly affected the composition of the microbiota, facilitating the dominance of distinct species such as Clostridium hathewayi, which were also preferentially over-represented in stools from RCC patients compared with HVs. Importantly, TKIs taken prior to nivolumab had implications in shifting the microbiota composition. To establish a cause-effect relationship between gut bacteria composition and ICB efficacy, NR-FMT mice were successfully compensated with either FMT from responding RCC patients or beneficial commensals identified by WGS-MG (Akkermansia muciniphila and Bacteroides salyersiae). CONCLUSIONS: The composition of the microbiota is influenced by TKIs and ATBs, and impacts the success of immunotherapy. Future studies will help sharpen the role of these specific bacteria and their potential as new biomarkers. PATIENT SUMMARY: We used quantitative shotgun DNA sequencing of fecal microbes as well as preclinical models of fecal or bacterial transfer to study the association between stool composition and (pre)clinical outcome to immune checkpoint blockade. Novel insights into the pathophysiological relevance of intestinal dysbiosis in the prognosis of kidney cancer may lead to innovative therapeutic solutions, such as supplementation with probiotics to prevent primary resistance to therapy.

Mediterranean diet intervention in overweight and obese subjects lowers plasma cholesterol and causes changes in the gut microbiome and metabolome independently of energy intake.

OBJECTIVES: This study aimed to explore the effects of an isocaloric Mediterranean diet (MD) intervention on metabolic health, gut microbiome and systemic metabolome in subjects with lifestyle risk factors for metabolic disease. DESIGN: Eighty-two healthy overweight and obese subjects with a habitually low intake of fruit and vegetables and a sedentary lifestyle participated in a parallel 8-week randomised controlled trial. Forty-three participants consumed an MD tailored to their habitual energy intakes (MedD), and 39 maintained their regular diets (ConD). Dietary adherence, metabolic parameters, gut microbiome and systemic metabolome were monitored over the study period. RESULTS: Increased MD adherence in the MedD group successfully reprogrammed subjects' intake of fibre and animal proteins. Compliance was confirmed by lowered levels of carnitine in plasma and urine. Significant reductions in plasma cholesterol (primary outcome) and faecal bile acids occurred in the MedD compared with the ConD group. Shotgun metagenomics showed gut microbiome changes that reflected individual MD adherence and increase in gene richness in participants who reduced systemic inflammation over the intervention. The MD intervention led to increased levels of the fibre-degrading Faecalibacterium prausnitzii and of genes for microbial carbohydrate degradation linked to butyrate metabolism. The dietary changes in the MedD group led to increased urinary urolithins, faecal bile acid degradation and insulin sensitivity that co-varied with specific microbial taxa. CONCLUSION: Switching subjects to an MD while maintaining their energy intake reduced their blood cholesterol and caused multiple changes in their microbiome and metabolome that are relevant in future strategies for the improvement of metabolic health.

Identification of rumen microbial biomarkers linked to methane emission in Holstein dairy cows.

Mitigation of greenhouse gas emissions is relevant for reducing the environmental impact of ruminant production. In this study, the rumen microbiome from Holstein cows was characterized through a combination of 16S rRNA gene and shotgun metagenomic sequencing. Methane production (CH4 ) and dry matter intake (DMI) were individually measured over 4-6 weeks to calculate the CH4 yield (CH4 y = CH4 /DMI) per cow. We implemented a combination of clustering, multivariate and mixed model analyses to identify a set of operational taxonomic unit (OTU) jointly associated with CH4 y and the structure of ruminal microbial communities. Three ruminotype clusters (R1, R2 and R3) were identified, and R2 was associated with higher CH4 y. The taxonomic composition on R2 had lower abundance of Succinivibrionaceae and Methanosphaera, and higher abundance of Ruminococcaceae, Christensenellaceae and Lachnospiraceae. Metagenomic data confirmed the lower abundance of Succinivibrionaceae and Methanosphaera in R2 and identified genera (Fibrobacter and unclassified Bacteroidales) not highlighted by metataxonomic analysis. In addition, the functional metagenomic analysis revealed that samples classified in cluster R2 were overrepresented by genes coding for KEGG modules associated with methanogenesis, including a significant relative abundance of the methyl-coenzyme M reductase enzyme. Based on the cluster assignment, we applied a sparse partial least-squares discriminant analysis at the taxonomic and functional levels. In addition, we implemented a sPLS regression model using the phenotypic variation of CH4 y. By combining these two approaches, we identified 86 discriminant bacterial OTUs, notably including families linked to CH4 emission such as Succinivibrionaceae, Ruminococcaceae, Christensenellaceae, Lachnospiraceae and Rikenellaceae. These selected OTUs explained 24% of the CH4 y phenotypic variance, whereas the host genome contribution was ~14%. In summary, we identified rumen microbial biomarkers associated with the methane production of dairy cows; these biomarkers could be used for targeted methane-reduction selection programmes in the dairy cattle industry provided they are heritable.