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.

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.

Dense sampling of bird diversity increases power of comparative genomics.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

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.

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.

Oral Microbiome in Relation to Periodontitis Severity and Systemic Inflammation.

Systemic inflammation induced by periodontitis is suggested to be the link between periodontitis and cardiovascular disease. The aim of this work was to explore the oral microbiome in periodontitis in relation to disease severity and systemic inflammation. The saliva and subgingival microbiome from periodontal pocket samples of patients with severe (n = 12) and mild periodontitis (n = 13) were analyzed using metagenomic shotgun sequencing. The taxa and pathways abundances were quantified. The diversity was assessed and the abundances to phenotype associations were performed using ANCOM and linear regression. A panel of inflammatory markers was measured in blood and was associated with taxa abundance. The microbial diversity and species richness did not differ between severe and mild periodontitis in either saliva or periodontal pockets. However, there were significant differences in the microbial composition between severe and mild periodontitis in the subgingival microbiome (i.e., pocket samples) and, in a lower grade, in saliva, and this is positively associated with systemic inflammatory markers. The "red complex" and "cluster B" abundances in periodontal pockets were strongly associated with inflammatory markers interleukin-6 and the white blood cell count. Our data suggest that systemic inflammation in severe periodontitis may be driven by the oral microbiome and may support the indirect (inflammatory) mechanism for the association between periodontitis and cardiovascular disease.

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.

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).

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.

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.

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.

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.

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.

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.

Western and Carnivorous Dietary Patterns are Associated with Greater Likelihood of IBD Development in a Large Prospective Population-based Cohort.

OBJECTIVE: Nutrition plays a role in the development of Crohn's disease [CD] and ulcerative colitis [UC]. However, prospective data on nutrition and disease onset are limited. Here, we analysed dietary patterns and scores in relation to inflammatory bowel disease [IBD] development in a prospective population-based cohort. METHODS: We analysed 125 445 participants of whom 224 individuals developed de novo UC and 97 CD over a maximum 14-year follow-up period. Participants answered health-related [also prospectively] and dietary questionnaires [FFQ] at baseline. Principal component analysis [PCA] was conducted deriving a-posteriori dietary patterns. Hypotheses-based a-priori dietary scores were also calculated, including the protein score, Healthy Eating Index, LifeLines Diet Score [LLDS], and alternative Mediterranean Diet Score. Logistic regression models were performed between dietary patterns, scores, and IBD development. RESULTS: PCA identified five dietary patterns. A pattern characterised by high intake of snacks, prepared meals, non-alcoholic beverages, and sauces along with low vegetables and fruit consumption was associated with higher likelihood of CD development (odds ratio [OR]: 1.16, 95% confidence interval [CI]: 1.03-1.30, p = 0.013). A pattern comprising red meat, poultry, and processed meat, was associated with increased likelihood of UC development [OR: 1.11, 95% CI: 1.01-1.20, p = 0.023]. A high diet quality score [LLDS] was associated with decreased risk of CD [OR: 0.95, 95% CI: 0.92-0.99, p = 0.009]. CONCLUSIONS: A Western dietary pattern was associated with a greater likelihood of CD development and a carnivorous pattern with UC development, whereas a relatively high diet quality [LLDS] was protective for CD development. Our study strengthens the importance of evaluating dietary patterns to aid prevention of IBD in the general population.

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.