Rome III Criteria Capture Higher Irritable Bowel Syndrome SNP-Heritability and Highlight a Novel Genetic Link with Cardiovascular Traits.

BACKGROUND & AIMS: Irritable bowel syndrome (IBS) shows genetic predisposition, and large-scale genome-wide association studies (GWAS) are emerging, based on heterogeneous disease definitions. We investigated the genetic architecture of IBS defined according to gold standard Rome Criteria. METHODS: We conducted GWAS meta-analyses of Rome III IBS and its subtypes in 24,735 IBS cases and 77,149 asymptomatic control subjects from 2 independent European cohorts (UK Biobank and Lifelines). Single-nucleotide polymorphism (SNP)-based heritability (h2SNP) and genetic correlations (rg) with other traits were calculated. IBS risk loci were functionally annotated to identify candidate genes. Sensitivity and conditional analyses were conducted to assess impact of confounders. Polygenic risk scores were computed and tested in independent datasets. RESULTS: Rome III IBS showed significant SNP-heritability (up to 13%) and similar genetic architecture across subtypes, including those with manifestations at the opposite ends of the symptom spectrum (rg = 0.48 between IBS-D and IBS-C). Genetic correlations with other traits highlighted commonalities with family history of heart disease and hypertension, coronary artery disease, and angina pectoris (rg = 0.20-0.45), among others. Four independent GWAS signals (P < 5×10-8) were detected, including 2 novel loci for IBS (rs2035380) and IBS-mixed (rs2048419) that had been previously associated with hypertension and coronary artery disease. Functional annotation of GWAS risk loci revealed genes implicated in circadian rhythm (BMAL1), intestinal barrier (CLDN23), immunomodulation (MFHAS1), and the cyclic adenosine monophosphate pathway (ADCY2). Polygenic risk scores allowed the identification of individuals at increased risk of IBS (odds ratio, 1.34; P = 1.1×10-3). CONCLUSIONS: Rome III Criteria capture higher SNP-heritability than previously estimated for IBS. The identified link between IBS and cardiovascular traits may contribute to the delineation of alternative therapeutic strategies, warranting further investigation.

Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

Transmission and dynamics of mother-infant gut viruses during pregnancy and early life.

Early development of the gut ecosystem is crucial for lifelong health. While infant gut bacterial communities have been studied extensively, the infant gut virome remains under-explored. To study the development of the infant gut virome over time and the factors that shape it, we longitudinally assess the composition of gut viruses and their bacterial hosts in 30 women during and after pregnancy and in their 32 infants during their first year of life. Using shotgun metagenomic sequencing applied to dsDNA extracted from Virus-Like Particles (VLPs) and bacteria, we generate 205 VLP metaviromes and 322 total metagenomes. With this data, we show that while the maternal gut virome composition remains stable during late pregnancy and after birth, the infant gut virome is dynamic in the first year of life. Notably, infant gut viromes contain a higher abundance of active temperate phages compared to maternal gut viromes, which decreases over the first year of life. Moreover, we show that the feeding mode and place of delivery influence the gut virome composition of infants. Lastly, we provide evidence of co-transmission of viral and bacterial strains from mothers to infants, demonstrating that infants acquire some of their virome from their mother's gut.

Choice of DNA extraction method affects stool microbiome recovery and subsequent phenotypic association analyses.

The lack of standardization in the methods of DNA extraction from fecal samples represents the major source of experimental variation in the microbiome research field. In this study, we aimed to compare the metagenomic profiles and microbiome-phenotype associations obtained by applying two commercially available DNA extraction kits: the AllPrep DNA/RNA Mini Kit (APK) and the QIAamp Fast DNA Stool Mini Kit (FSK). Using metagenomic sequencing data available from 745 paired fecal samples from two independent population cohorts, Lifelines-DEEP (LLD, n = 292) and the 500 Functional Genomics project (500FG, n = 453), we confirmed significant differences in DNA yield and the recovered microbial communities between protocols, with the APK method resulting in a higher DNA concentration and microbial diversity. Further, we observed a massive difference in bacterial relative abundances at species-level between the APK and the FSK protocols, with > 75% of species differentially abundant between protocols in both cohorts. Specifically, comparison with a standard mock community revealed that the APK method provided higher accuracy in the recovery of microbial relative abundances, with the absence of a bead-beating step in the FSK protocol causing an underrepresentation of gram-positive bacteria. This heterogeneity in the recovered microbial composition led to remarkable differences in the association with anthropometric and lifestyle phenotypes. The results of this study further reinforce that the choice of DNA extraction method impacts the metagenomic profile of human gut microbiota and highlight the importance of harmonizing protocols in microbiome studies.

PICALO: principal interaction component analysis for the identification of discrete technical, cell-type, and environmental factors that mediate eQTLs.

Expression quantitative trait loci (eQTL) offer insights into the regulatory mechanisms of trait-associated variants, but their effects often rely on contexts that are unknown or unmeasured. We introduce PICALO, a method for hidden variable inference of eQTL contexts. PICALO identifies and disentangles technical from biological context in heterogeneous blood and brain bulk eQTL datasets. These contexts are biologically informative and reproducible, outperforming cell counts or expression-based principal components. Furthermore, we show that RNA quality and cell type proportions interact with thousands of eQTLs. Knowledge of hidden eQTL contexts may aid in the inference of functional mechanisms underlying disease variants.

Host genetic regulation of human gut microbial structural variation.

Although the impact of host genetics on gut microbial diversity and the abundance of specific taxa is well established1-6, little is known about how host genetics regulates the genetic diversity of gut microorganisms. Here we conducted a meta-analysis of associations between human genetic variation and gut microbial structural variation in 9,015 individuals from four Dutch cohorts. Strikingly, the presence rate of a structural variation segment in Faecalibacterium prausnitzii that harbours an N-acetylgalactosamine (GalNAc) utilization gene cluster is higher in individuals who secrete the type A oligosaccharide antigen terminating in GalNAc, a feature that is jointly determined by human ABO and FUT2 genotypes, and we could replicate this association in a Tanzanian cohort. In vitro experiments demonstrated that GalNAc can be used as the sole carbohydrate source for F. prausnitzii strains that carry the GalNAc-metabolizing pathway. Further in silico and in vitro studies demonstrated that other ABO-associated species can also utilize GalNAc, particularly Collinsella aerofaciens. The GalNAc utilization genes are also associated with the host's cardiometabolic health, particularly in individuals with mucosal A-antigen. Together, the findings of our study demonstrate that genetic associations across the human genome and bacterial metagenome can provide functional insights into the reciprocal host-microbiome relationship.

Health-related quality of life is linked to the gut microbiome in kidney transplant recipients.

Kidney transplant recipients (KTR) have impaired health-related quality of life (HRQoL) and suffer from intestinal dysbiosis. Increasing evidence shows that gut health and HRQoL are tightly related in the general population. Here, we investigate the association between the gut microbiome and HRQoL in KTR, using metagenomic sequencing data from fecal samples collected from 507 KTR. Multiple bacterial species are associated with lower HRQoL, many of which have previously been associated with adverse health conditions. Gut microbiome distance to the general population is highest among KTR with an impaired physical HRQoL (R = -0.20, P = 2.3 × 10-65) and mental HRQoL (R = -0.14, P = 1.3 × 10-3). Physical and mental HRQoL explain a significant part of variance in the gut microbiome (R2 = 0.58%, FDR = 5.43 × 10-4 and R2 = 0.37%, FDR = 1.38 × 10-3, respectively). Additionally, multiple metabolic and neuroactive pathways (gut brain modules) are associated with lower HRQoL. While the observational design of our study does not allow us to analyze causality, we provide a comprehensive overview of the associations between the gut microbiome and HRQoL while controlling for confounders.

The gut microbiome across the cardiovascular risk spectrum.

AIMS: Despite treatment advancements, cardiovascular disease remains a leading cause of death worldwide. Identifying new targets is crucial for enhancing preventive and therapeutic strategies. The gut microbiome has been associated with coronary artery disease (CAD), however our understanding of specific changes during CAD development remains limited. We aimed to investigate microbiome changes in participants without clinically manifest CAD with different cardiovascular risk levels and in patients with ST-elevation myocardial infarction (STEMI). METHODS: In this cross-sectional study, we characterized the gut microbiome using metagenomics of 411 faecal samples from individuals with low (n=130), intermediate (n=130) and high (n=125) cardiovascular risk based on the Framingham score, and STEMI patients (n=26). We analysed diversity, and differential abundance of species and functional pathways while accounting for confounders including medication and technical covariates. RESULTS: Collinsella stercoris, Flavonifractor plautii and Ruthenibacterium lactatiformans showed increased abundances with cardiovascular risk, while Streptococcus thermophilus was negatively associated. Differential abundance analysis revealed eight species and 49 predicted metabolic pathways that were differently abundant among the groups. In the gut microbiome of STEMI patients, there was a depletion of pathways linked to vitamin, lipid and amino-acid biosynthesis. CONCLUSION: We identified four microbial species showing a gradual trend in abundance from low-risk individuals to those with STEMI, and observed differential abundant species and pathways in STEMI patients compared to those without clinically manifest CAD. Further investigation is warranted to gain deeper understanding of their precise role in CAD progression and potential implications, with the ultimate goal of identifying novel therapeutic targets.

Despite previous studies demonstrating dysbiosis in STEMI patients, our understanding of the precise microbiome changes across the cardiovascular risk spectrum remains limited. This study addresses this knowledge gap by providing insights into the gut microbiome composition of individuals across varying cardiovascular risk levels and STEMI patients. By examining the gut microbiome of carefully selected participants from the general population with three different risk levels and a unique group of STEMI patients, we identified microbial species and pathways with differential abundance across the groups. Several of these species and pathways are associated with inflammation and lipid metabolism, which are key factors in CAD development. Collinsella stercoris, Flavonifractor plautii and Ruthenibacterium lactatiformans are increasingly abundant, while Streptococcus thermophilus is decreasingly abundant across the cardiovascular risk spectrum. The gut microbiome of STEMI patients showed eight differentially abundant species compared to groups at risk. Notably, four of these species, characterized by an elevated abundance in STEMI patients, have not been previously reported.

Analysis of microbial composition and sharing in low-biomass human milk samples: a comparison of DNA isolation and sequencing techniques.

Human milk microbiome studies are currently hindered by low milk bacterial/human cell ratios and often rely on 16S rRNA gene sequencing, which limits downstream analyses. Here, we aimed to find a method to study milk bacteria and assess bacterial sharing between maternal and infant microbiota. We tested four DNA isolation methods, two bacterial enrichment methods and three sequencing methods on mock communities, milk samples and negative controls. Of the four DNA isolation kits, the DNeasy PowerSoil Pro (PS) and MagMAX Total Nucleic Acid Isolation (MX) kits provided consistent 16S rRNA gene sequencing results with low contamination. Neither enrichment method substantially decreased the human metagenomic sequencing read-depth. Long-read 16S-ITS-23S rRNA gene sequencing biased the mock community composition but provided consistent results for milk samples, with little contamination. In contrast to 16S rRNA gene sequencing, 16S-ITS-23S rRNA gene sequencing of milk, infant oral, infant faecal and maternal faecal DNA from 14 mother-infant pairs provided sufficient resolution to detect significantly more frequent sharing of bacteria between related pairs compared to unrelated pairs. In conclusion, PS or MX kit-DNA isolation followed by 16S rRNA gene sequencing reliably characterises human milk microbiota, and 16S-ITS-23S rRNA gene sequencing enables studies of bacterial transmission in low-biomass samples.

Gut dysbiosis associates with cytokine production capacity in viral-suppressed people living with HIV.

Background: People living with human immunodeficiency virus (PLHIV) are exposed to chronic immune dysregulation, even when virus replication is suppressed by antiretroviral therapy (ART). Given the emerging role of the gut microbiome in immunity, we hypothesized that the gut microbiome may be related to the cytokine production capacity of PLHIV. Methods: To test this hypothesis, we collected metagenomic data from 143 ART-treated PLHIV and assessed the ex vivo production capacity of eight different cytokines [interleukin-1ß (IL-1ß), IL-6, IL-1Ra, IL-10, IL-17, IL-22, tumor necrosis factor, and interferon-γ] in response to different stimuli. We also characterized CD4+ T-cell counts, HIV reservoir, and other clinical parameters. Results: Compared with 190 age- and sex-matched controls and a second independent control cohort, PLHIV showed microbial dysbiosis that was correlated with viral reservoir levels (CD4+ T-cell-associated HIV-1 DNA), cytokine production capacity, and sexual behavior. Notably, we identified two genetically different P. copri strains that were enriched in either PLHIV or healthy controls. The control-related strain showed a stronger negative association with cytokine production capacity than the PLHIV-related strain, particularly for Pam3Cys-incuded IL-6 and IL-10 production. The control-related strain is also positively associated with CD4+ T-cell level. Conclusions: Our findings suggest that modulating the gut microbiome may be a strategy to modulate immune response in PLHIV.

Impact of occupational pesticide exposure on the human gut microbiome.

The rising use of pesticides in modern agriculture has led to a shift in disease burden in which exposure to these chemicals plays an increasingly important role. The human gut microbiome, which is partially responsible for the biotransformation of xenobiotics, is also known to promote biotransformation of environmental pollutants. Understanding the effects of occupational pesticide exposure on the gut microbiome can thus provide valuable insights into the mechanisms underlying the impact of pesticide exposure on health. Here we investigate the impact of occupational pesticide exposure on human gut microbiome composition in 7198 participants from the Dutch Microbiome Project of the Lifelines Study. We used job-exposure matrices in combination with occupational codes to retrieve categorical and cumulative estimates of occupational exposures to general pesticides, herbicides, insecticides and fungicides. Approximately 4% of our cohort was occupationally exposed to at least one class of pesticides, with predominant exposure to multiple pesticide classes. Most participants reported long-term employment, suggesting a cumulative profile of exposure. We demonstrate that contact with insecticides, fungicides and a general "all pesticides" class was consistently associated with changes in the gut microbiome, showing significant associations with decreased alpha diversity and a differing beta diversity. We also report changes in the abundance of 39 different bacterial taxa upon exposure to the different pesticide classes included in this study. Together, the extent of statistically relevant associations between gut microbial changes and pesticide exposure in our findings highlights the impact of these compounds on the human gut microbiome.

Phage display sequencing reveals that genetic, environmental, and intrinsic factors influence variation of human antibody epitope repertoire.

Phage-displayed immunoprecipitation sequencing (PhIP-seq) has enabled high-throughput profiling of human antibody repertoires. However, a comprehensive overview of environmental and genetic determinants shaping human adaptive immunity is lacking. In this study, we investigated the effects of genetic, environmental, and intrinsic factors on the variation in human antibody repertoires. We characterized serological antibody repertoires against 344,000 peptides using PhIP-seq libraries from a wide range of microbial and environmental antigens in 1,443 participants from a population cohort. We detected individual-specificity, temporal consistency, and co-housing similarities in antibody repertoires. Genetic analyses showed the involvement of the HLA, IGHV, and FUT2 gene regions in antibody-bound peptide reactivity. Furthermore, we uncovered associations between phenotypic factors (including age, cell counts, sex, smoking behavior, and allergies, among others) and particular antibody-bound peptides. Our results indicate that human antibody epitope repertoires are shaped by both genetics and environmental exposures and highlight specific signatures of distinct phenotypes and genotypes.

Phage-display immunoprecipitation sequencing of the antibody epitope repertoire in inflammatory bowel disease reveals distinct antibody signatures.

Inflammatory bowel diseases (IBDs), e.g., Crohn's disease (CD) and ulcerative colitis (UC), are chronic immune-mediated inflammatory diseases. A comprehensive overview of an IBD-specific antibody epitope repertoire is, however, lacking. Using high-throughput phage-display immunoprecipitation sequencing (PhIP-Seq), we identified antibodies against 344,000 antimicrobial, immune, and food antigens in 497 individuals with IBD compared with 1,326 controls. IBD was characterized by 373 differentially abundant antibody responses (202 overrepresented and 171 underrepresented), with 17% shared by both IBDs, 55% unique to CD, and 28% unique to UC. Antibody reactivities against bacterial flagellins dominated in CD and were associated with ileal involvement, fibrostenotic disease, and anti-Saccharomyces cerevisiae antibody positivity, but not with fecal microbiome composition. Antibody epitope repertoires accurately discriminated CD from controls (area under the curve [AUC] = 0.89), and similar discrimination was achieved when using only ten antibodies (AUC = 0.87). Individuals with IBD thus show a distinct antibody repertoire against selected peptides, allowing clinical stratification and discovery of immunological targets.

The maternal gut microbiome during pregnancy and its role in maternal and infant health.

There is growing knowledge that the maternal gut microbiome undergoes substantial changes during pregnancy. However, despite the recognition that the maternal gut microbiome influences maternal and infant health, we still have a limited understanding of the clinical and environmental factors that can impact the maternal gut microbiome during pregnancy and the consequences of these changes. Here, we review the current body of knowledge about factors shaping the maternal gut microbiome during pregnancy and its role in the development of pregnancy complications and infant health.

Faecal metabolome and its determinants in inflammatory bowel disease.

OBJECTIVE: Inflammatory bowel disease (IBD) is a multifactorial immune-mediated inflammatory disease of the intestine, comprising Crohn's disease and ulcerative colitis. By characterising metabolites in faeces, combined with faecal metagenomics, host genetics and clinical characteristics, we aimed to unravel metabolic alterations in IBD. DESIGN: We measured 1684 different faecal metabolites and 8 short-chain and branched-chain fatty acids in stool samples of 424 patients with IBD and 255 non-IBD controls. Regression analyses were used to compare concentrations of metabolites between cases and controls and determine the relationship between metabolites and each participant's lifestyle, clinical characteristics and gut microbiota composition. Moreover, genome-wide association analysis was conducted on faecal metabolite levels. RESULTS: We identified over 300 molecules that were differentially abundant in the faeces of patients with IBD. The ratio between a sphingolipid and L-urobilin could discriminate between IBD and non-IBD samples (AUC=0.85). We found changes in the bile acid pool in patients with dysbiotic microbial communities and a strong association between faecal metabolome and gut microbiota. For example, the abundance of Ruminococcus gnavus was positively associated with tryptamine levels. In addition, we found 158 associations between metabolites and dietary patterns, and polymorphisms near NAT2 strongly associated with coffee metabolism. CONCLUSION: In this large-scale analysis, we identified alterations in the metabolome of patients with IBD that are independent of commonly overlooked confounders such as diet and surgical history. Considering the influence of the microbiome on faecal metabolites, our results pave the way for future interventions targeting intestinal inflammation.

gutSMASH predicts specialized primary metabolic pathways from the human gut microbiota.

The gut microbiota produce hundreds of small molecules, many of which modulate host physiology. Although efforts have been made to identify biosynthetic genes for secondary metabolites, the chemical output of the gut microbiome consists predominantly of primary metabolites. Here we introduce the gutSMASH algorithm for identification of primary metabolic gene clusters, and we used it to systematically profile gut microbiome metabolism, identifying 19,890 gene clusters in 4,240 high-quality microbial genomes. We found marked differences in pathway distribution among phyla, reflecting distinct strategies for energy capture. These data explain taxonomic differences in short-chain fatty acid production and suggest a characteristic metabolic niche for each taxon. Analysis of 1,135 individuals from a Dutch population-based cohort shows that the level of microbiome-derived metabolites in plasma and feces is almost completely uncorrelated with the metagenomic abundance of corresponding metabolic genes, indicating a crucial role for pathway-specific gene regulation and metabolite flux. This work is a starting point for understanding differences in how bacterial taxa contribute to the chemistry of the microbiome.

Microbial vitamin production mediates dietary effects on diabetic risk.

Adequate levels of essential vitamins are important for the prevention of diabetes. While the main efforts to address this are currently focused on the intake of vitamin supplements, improving and maintaining intrinsic vitamin production capacity, which is determined by gut microbes, has received insufficient attention. In this study, we systematically investigated the relationship between gut microbial vitamin production and factors related to diabetes and cardiometabolic health in a deeply phenotyped cohort, Lifelines-DEEP (N = 1,135). We found that blood glucose-related factors, lipids, circulating inflammation, and fecal short-chain fatty acids are associated with gut microbial vitamin production. Use of laxatives and metformin are associated with increased levels of vitamin B1/B6 biosynthesis pathways. We further reveal a mediatory role for microbial vitamin B1/B2 production on the influence of fruit intake on diabetes risk. This study provides preliminary evidence for microbiome-targeted vitamin metabolism interventions to promote health.

Systematic analysis of relationships between plasma branched-chain amino acid concentrations and cardiometabolic parameters: an association and Mendelian randomization study.

BACKGROUND: Branched-chain amino acids (BCAAs; valine, leucine, and isoleucine) are essential amino acids that are associated with an increased risk of cardiometabolic diseases (CMD). However, there are still only limited insights into potential direct associations between BCAAs and a wide range of CMD parameters, especially those remaining after correcting for covariates and underlying causal relationships. METHODS: To shed light on these relationships, we systematically characterized the associations between plasma BCAA concentrations and a large panel of 537 CMD parameters (including atherosclerosis-related parameters, fat distribution, plasma cytokine concentrations and cell counts, circulating concentrations of cardiovascular-related proteins and plasma metabolites) in 1400 individuals from the Dutch population cohort LifeLines DEEP and 294 overweight individuals from the 300OB cohort. After correcting for age, sex, and BMI, we assessed associations between individual BCAAs and CMD parameters. We further assessed the underlying causality using Mendelian randomization. RESULTS: A total of 838 significant associations were detected for 409 CMD parameters. BCAAs showed both common and specific associations, with the most specific associations being detected for isoleucine. Further, we found that obesity status substantially affected the strength and direction of associations for valine, which cannot be corrected for using BMI as a covariate. Subsequent univariable Mendelian randomization (UVMR), after removing BMI-associated SNPs, identified seven significant causal relationships from four CMD traits to BCAA levels, mostly for diabetes-related parameters. However, no causal effects of BCAAs on CMD parameters were supported. CONCLUSIONS: Our cross-sectional association study reports a large number of associations between BCAAs and CMD parameters. Our results highlight some specific associations for isoleucine, as well as obesity-specific effects for valine. MR-based causality analysis suggests that altered BCAA levels can be a consequence of diabetes and alteration in lipid metabolism. We found no MR evidence to support a causal role for BCAAs in CMD. These findings provide evidence to (re)evaluate the clinical importance of individual BCAAs in CMD diagnosis, prevention, and treatment.

Influence of the microbiome, diet and genetics on inter-individual variation in the human plasma metabolome.

The levels of the thousands of metabolites in the human plasma metabolome are strongly influenced by an individual's genetics and the composition of their diet and gut microbiome. Here, by assessing 1,183 plasma metabolites in 1,368 extensively phenotyped individuals from the Lifelines DEEP and Genome of the Netherlands cohorts, we quantified the proportion of inter-individual variation in the plasma metabolome explained by different factors, characterizing 610, 85 and 38 metabolites as dominantly associated with diet, the gut microbiome and genetics, respectively. Moreover, a diet quality score derived from metabolite levels was significantly associated with diet quality, as assessed by a detailed food frequency questionnaire. Through Mendelian randomization and mediation analyses, we revealed putative causal relationships between diet, the gut microbiome and metabolites. For example, Mendelian randomization analyses support a potential causal effect of Eubacterium rectale in decreasing plasma levels of hydrogen sulfite-a toxin that affects cardiovascular function. Lastly, based on analysis of the plasma metabolome of 311 individuals at two time points separated by 4 years, we observed a positive correlation between the stability of metabolite levels and the amount of variance in the levels of that metabolite that could be explained in our analysis. Altogether, characterization of factors that explain inter-individual variation in the plasma metabolome can help design approaches for modulating diet or the gut microbiome to shape a healthy metabolome.

Diversity and Ecology of Caudoviricetes Phages with Genome Terminal Repeats in Fecal Metagenomes from Four Dutch Cohorts.

The human gut harbors numerous viruses infecting the human host, microbes, and other inhabitants of the gastrointestinal tract. Most of these viruses remain undiscovered, and their influence on human health is unknown. Here, we characterize viral genomes in gut metagenomic data from 1950 individuals from four population and patient cohorts. We focus on a subset of viruses that is highly abundant in the gut, remains largely uncharacterized, and allows confident complete genome identification­phages that belong to the class Caudoviricetes and possess genome terminal repeats. We detect 1899 species-level units belonging to this subset, 19% of which do not have complete representative genomes in major public gut virome databases. These units display diverse genomic features, are predicted to infect a wide range of microbial hosts, and on average account for <1% of metagenomic reads. Analysis of longitudinal data from 338 individuals shows that the composition of this fraction of the virome remained relatively stable over a period of 4 years. We also demonstrate that 54 species-level units are highly prevalent (detected in >5% of individuals in a cohort). Finally, we find 34 associations between highly prevalent phages and human phenotypes, 24 of which can be explained by the relative abundance of potential hosts.