Epigenome-wide association study on the plasma metabolome suggests self-regulation of the glycine and serine pathway through DNA methylation.

BACKGROUND: The plasma metabolome reflects the physiological state of various biological processes and can serve as a proxy for disease risk. Plasma metabolite variation, influenced by genetic and epigenetic mechanisms, can also affect the cellular microenvironment and blood cell epigenetics. The interplay between the plasma metabolome and the blood cell epigenome remains elusive. In this study, we performed an epigenome-wide association study (EWAS) of 1183 plasma metabolites in 693 participants from the LifeLines-DEEP cohort and investigated the causal relationships in DNA methylation-metabolite associations using bidirectional Mendelian randomization and mediation analysis. RESULTS: After rigorously adjusting for potential confounders, including genetics, we identified five robust associations between two plasma metabolites (L-serine and glycine) and three CpG sites located in two independent genomic regions (cg14476101 and cg16246545 in PHGDH and cg02711608 in SLC1A5) at a false discovery rate of less than 0.05. Further analysis revealed a complex bidirectional relationship between plasma glycine/serine levels and DNA methylation. Moreover, we observed a strong mediating role of DNA methylation in the effect of glycine/serine on the expression of their metabolism/transport genes, with the proportion of the mediated effect ranging from 11.8 to 54.3%. This result was also replicated in an independent population-based cohort, the Rotterdam Study. To validate our findings, we conducted in vitro cell studies which confirmed the mediating role of DNA methylation in the regulation of PHGDH gene expression. CONCLUSIONS: Our findings reveal a potential feedback mechanism in which glycine and serine regulate gene expression through DNA methylation.

Intestinal Blastocystis is linked to healthier diets and more favorable cardiometabolic outcomes in 56,989 individuals from 32 countries.

Diet impacts human health, influencing body adiposity and the risk of developing cardiometabolic diseases. The gut microbiome is a key player in the diet-health axis, but while its bacterial fraction is widely studied, the role of micro-eukaryotes, including Blastocystis, is underexplored. We performed a global-scale analysis on 56,989 metagenomes and showed that human Blastocystis exhibits distinct prevalence patterns linked to geography, lifestyle, and dietary habits. Blastocystis presence defined a specific bacterial signature and was positively associated with more favorable cardiometabolic profiles and negatively with obesity (p < 1e-16) and disorders linked to altered gut ecology (p < 1e-8). In a diet intervention study involving 1,124 individuals, improvements in dietary quality were linked to weight loss and increases in Blastocystis prevalence (p = 0.003) and abundance (p < 1e-7). Our findings suggest a potentially beneficial role for Blastocystis, which may help explain personalized host responses to diet and downstream disease etiopathogenesis.

Reproducible gut microbial signatures in bipolar and schizophrenia spectrum disorders: A metagenome-wide study.

BACKGROUND: Numerous studies report gut microbiome variations in bipolar disorder (BD) and schizophrenia spectrum disorders (SSD) compared to healthy individuals, though, there is limited consensus on which specific bacteria are associated with these disorders. METHODS: In this study, we performed a comprehensive metagenomic shotgun sequencing analysis in 103 Dutch patients with BD/SSD and 128 healthy controls matched for age, sex, body mass index and income, while accounting for diet quality, transit time and technical confounders. To assess the replicability of the findings, we used two validation cohorts (total n = 203), including participants from a distinct population with a different metagenomic isolation protocol. RESULTS: The gut microbiome of the patients had a significantly different ß-diversity, but not α-diversity nor neuroactive potential compared to healthy controls. Initially, twenty-six bacterial taxa were identified as differentially abundant in patients. Among these, the previously reported genera Lachnoclostridium and Eggerthella were replicated in the validation cohorts. Employing the CoDaCoRe learning algorithm, we identified two bacterial balances specific to BD/SSD, which demonstrated an area under the receiver operating characteristic curve (AUC) of 0.77 in the test dataset. These balances were replicated in the validation cohorts and showed a positive association with the severity of psychiatric symptoms and antipsychotic use. Last, we showed a positive association between the relative abundance of Klebsiella and Klebsiella pneumoniae with antipsychotic use and between the Anaeromassilibacillus and lithium use. CONCLUSIONS: Our findings suggest that microbial balances could be a reproducible method for identifying BD/SSD-specific microbial signatures, with potential diagnostic and prognostic applications. Notably, Lachnoclostridium and Eggerthella emerge as frequently occurring bacteria in BD/SSD. Last, our study reaffirms the previously established link between Klebsiella and antipsychotic medication use and identifies a novel association between Anaeromassilibacillus and lithium use.

Antibody signatures against viruses and microbiome reflect past and chronic exposures and associate with aging and inflammation.

Encounters with pathogens and other molecules can imprint long-lasting effects on our immune system, influencing future physiological outcomes. Given the wide range of microbes to which humans are exposed, their collective impact on health is not fully understood. To explore relations between exposures and biological aging and inflammation, we profiled an antibody-binding repertoire against 2,815 microbial, viral, and environmental peptides in a population cohort of 1,443 participants. Utilizing antibody-binding as a proxy for past exposures, we investigated their impact on biological aging, cell composition, and inflammation. Immune response against cytomegalovirus (CMV), rhinovirus, and gut bacteria relates with telomere length. Single-cell expression measurements identified an effect of CMV infection on the transcriptional landscape of subpopulations of CD8 and CD4 T-cells. This examination of the relationship between microbial exposures and biological aging and inflammation highlights a role for chronic infections (CMV and Epstein-Barr virus) and common pathogens (rhinoviruses and adenovirus C).

Unraveling interindividual variation of trimethylamine N-oxide and its precursors at the population level.

Trimethylamine N-oxide (TMAO) is a circulating microbiome-derived metabolite implicated in the development of atherosclerosis and cardiovascular disease (CVD). We investigated whether plasma levels of TMAO, its precursors (betaine, carnitine, deoxycarnitine, choline), and TMAO-to-precursor ratios are associated with clinical outcomes, including CVD and mortality. This was followed by an in-depth analysis of their genetic, gut microbial, and dietary determinants. The analyses were conducted in five Dutch prospective cohort studies including 7834 individuals. To further investigate association results, Mendelian Randomization (MR) was also explored. We found only plasma choline levels (hazard ratio [HR] 1.17, [95% CI 1.07; 1.28]) and not TMAO to be associated with CVD risk. Our association analyses uncovered 10 genome-wide significant loci, including novel genomic regions for betaine (6p21.1, 6q25.3), choline (2q34, 5q31.1), and deoxycarnitine (10q21.2, 11p14.2) comprising several metabolic gene associations, for example, CPS1 or PEMT. Furthermore, our analyses uncovered 68 gut microbiota associations, mainly related to TMAO-to-precursors ratios and the Ruminococcaceae family, and 16 associations of food groups and metabolites including fish-TMAO, meat-carnitine, and plant-based food-betaine associations. No significant association was identified by the MR approach. Our analyses provide novel insights into the TMAO pathway, its determinants, and pathophysiological impact on the general population.

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.

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.

Selective enrichment of founding reproductive microbiomes allows extensive vertical transmission in a fungus-farming termite.

Mutualistic coevolution can be mediated by vertical transmission of symbionts between host generations. Termites host complex gut bacterial communities with evolutionary histories indicative of mixed-mode transmission. Here, we document that vertical transmission of gut bacterial strains is congruent across parent to offspring colonies in four pedigrees of the fungus-farming termite Macrotermes natalensis. We show that 44% of the offspring colony microbiome, including more than 80 bacterial genera and pedigree-specific strains, are consistently inherited. We go on to demonstrate that this is achieved because colony-founding reproductives are selectively enriched with a set of non-random, environmentally sensitive and termite-specific gut microbes from their colonies of origin. These symbionts transfer to offspring colony workers with high fidelity, after which priority effects appear to influence the composition of the establishing microbiome. Termite reproductives thus secure transmission of complex communities of specific, co-evolved microbes that are critical to their offspring colonies. Extensive yet imperfect inheritance implies that the maturing colony benefits from acquiring environmental microbes to complement combinations of termite, fungus and vertically transmitted microbes; a mode of transmission that is emerging as a prevailing strategy for hosts to assemble complex adaptive microbiomes.

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.

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.

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.

Beyond personal space: Unveiling the transmission pattern of the human gut and oral microbiome.

Valles-Colomer et al. evaluated the rate of microbial strain sharing to infer person-to-person microbial transmission. The presence of shared strains does not necessarily indicate direct transmission between individuals. Future research should investigate the potential for noncommunicable diseases to be transmitted through microbial strain transmission. Identification of harmful strains with high transmissibility and transmission routes can enable the development of targeted interventions.

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.

Genetic, parental and lifestyle factors influence telomere length.

The average length of telomere repeats (TL) declines with age and is considered to be a marker of biological ageing. Here, we measured TL in six blood cell types from 1046 individuals using the clinically validated Flow-FISH method. We identified remarkable cell-type-specific variations in TL. Host genetics, environmental, parental and intrinsic factors such as sex, parental age, and smoking are associated to variations in TL. By analysing the genome-wide methylation patterns, we identified that the association of maternal, but not paternal, age to TL is mediated by epigenetics. Single-cell RNA-sequencing data for 62 participants revealed differential gene expression in T-cells. Genes negatively associated with TL were enriched for pathways related to translation and nonsense-mediated decay. Altogether, this study addresses cell-type-specific differences in telomere biology and its relation to cell-type-specific gene expression and highlights how perinatal factors play a role in determining TL, on top of genetics and lifestyle.

Effect of host genetics on the gut microbiome in 7,738 participants of the Dutch Microbiome Project.

Host genetics are known to influence the gut microbiome, yet their role remains poorly understood. To robustly characterize these effects, we performed a genome-wide association study of 207 taxa and 205 pathways representing microbial composition and function in 7,738 participants of the Dutch Microbiome Project. Two robust, study-wide significant (P < 1.89 × 10-10) signals near the LCT and ABO genes were found to be associated with multiple microbial taxa and pathways and were replicated in two independent cohorts. The LCT locus associations seemed modulated by lactose intake, whereas those at ABO could be explained by participant secretor status determined by their FUT2 genotype. Twenty-two other loci showed suggestive evidence (P < 5 × 10-8) of association with microbial taxa and pathways. At a more lenient threshold, the number of loci we identified strongly correlated with trait heritability, suggesting that much larger sample sizes are needed to elucidate the remaining effects of host genetics on the gut microbiome.