A pilot study to identify blood-based markers associated with response to treatment with Vedolizumab in patients with Inflammatory Bowel Disease.

The inflammatory bowel diseases (IBD) known as Crohn's disease (CD) and ulcerative colitis (UC) are chronic inflammatory diseases of the gastrointestinal tract believed to arise because of an imbalance between the epithelial, immune and microbial systems. It has been shown that biological differences (genetic, epigenetic, microbial, environmental, etc.) exist between patients with IBD, with multiple risk factors been associated with disease susceptibility and IBD-related phenotypes (e.g. disease location). It is also known that there is heterogeneity in terms of response to therapy in patients with IBD, including to biological therapies that target very specific biological pathways (e.g. TNF-alpha signaling, IL-23R signaling, immune cell trafficking, etc.). It is hypothesized that the better the match between the biology targeted by these advanced therapies and the predominant disease-associated pathways at play in each patient will favor a beneficial response. The aim of this pilot study was to identify potential biological differences associated with differential treatment response to the anti α4ß7 integrin therapy known as Vedolizumab. Our approach was to measure a broad range of analytes in the serum of patients prior to initiation of therapy and at the first clinical assessment visit, to identify potential markers of biological differences between patients at baseline and to see which biomarkers are most affected by treatment in responders. Our focus on early clinical response was to study the most proximal effects of therapy and to minimize confounders such as loss of response that occurs further distal to treatment initiation. Specifically, we performed targeted analyses of >150 proteins and metabolites, and untargeted analyses of >1100 lipid entities, in serum samples from 92 IBD patients (42 CD, 50 UC) immediately prior to initiation of therapy with vedolizumab (baseline samples) and at their first clinical assessment (14-week samples). We found lower levels of SDF-1a, but higher levels of PDGF-ßß, lactate, lysine, phenylalanine, branched chain amino acids, alanine, short/medium chain acylcarnitines, and triglycerides containing myristic acid in baseline serum samples of responders as compared to non-responders. We also observed an increase in serum levels of CXCL9 and citrate, as well as a decrease in IL-10, between baseline and week 14 samples. In addition, we observed that a group of metabolites and protein analytes was strongly associated with both treatment response and BMI status, although BMI status was not associated with treatment response.

Chemical and topological design of multicapped mRNA and capped circular RNA to augment translation.

Protein and vaccine therapies based on mRNA would benefit from an increase in translation capacity. Here, we report a method to augment translation named ligation-enabled mRNA-oligonucleotide assembly (LEGO). We systematically screen different chemotopological motifs and find that a branched mRNA cap effectively initiates translation on linear or circular mRNAs without internal ribosome entry sites. Two types of chemical modification, locked nucleic acid (LNA) N7-methylguanosine modifications on the cap and LNA + 5 × 2' O-methyl on the 5' untranslated region, enhance RNA-eukaryotic translation initiation factor (eIF4E-eIF4G) binding and RNA stability against decapping in vitro. Through multidimensional chemotopological engineering of dual-capped mRNA and capped circular RNA, we enhanced mRNA protein production by up to tenfold in vivo, resulting in 17-fold and 3.7-fold higher antibody production after prime and boost doses in a severe acute respiratory syndrome coronavirus 2 vaccine setting, respectively. The LEGO platform opens possibilities to design unnatural RNA structures and topologies beyond canonical linear and circular RNAs for both basic research and therapeutic applications.

Commensal consortia decolonize Enterobacteriaceae via ecological control.

Persistent colonization and outgrowth of potentially pathogenic organisms in the intestine can result from long-term antibiotic use or inflammatory conditions, and may perpetuate dysregulated immunity and tissue damage1,2. Gram-negative Enterobacteriaceae gut pathobionts are particularly recalcitrant to conventional antibiotic treatment3,4, although an emerging body of evidence suggests that manipulation of the commensal microbiota may be a practical alternative therapeutic strategy5-7. Here we isolated and down-selected commensal bacterial consortia from stool samples from healthy humans that could strongly and specifically suppress intestinal Enterobacteriaceae. One of the elaborated consortia, comprising 18 commensal strains, effectively controlled ecological niches by regulating gluconate availability, thereby re-establishing colonization resistance and alleviating Klebsiella- and Escherichia-driven intestinal inflammation in mice. Harnessing these activities in the form of live bacterial therapies may represent a promising solution to combat the growing threat of proinflammatory, antimicrobial-resistant Enterobacteriaceae infection.

B cell tolerance and autoimmunity: Lessons from repertoires.

Adaptive immune cell function is regulated by a highly diverse receptor recombined from variable germline-encoded segments that can recognize an almost unlimited array of epitopes. While this diversity enables the recognition of any pathogen, it also poses a risk of self-recognition, leading to autoimmunity. Many layers of regulation are present during both the generation and activation of B cells to prevent this phenomenon, although they are evidently imperfect. In recent years, our ability to analyze immune repertoires at scale has drastically increased, both through advances in sequencing and single-cell analyses. Here, we review the current knowledge on B cell repertoire analyses, focusing on their implication for autoimmunity. These studies demonstrate that a failure of tolerance occurs at multiple independent checkpoints in different autoimmune contexts, particularly during B cell maturation, plasmablast differentiation, and within germinal centers. These failures are marked by distinct repertoire features that may be used to identify disease- or patient-specific therapeutic approaches.

Emerging biochemical, microbial and immunological evidence in the search for why HLA-B∗27 confers risk for spondyloarthritis.

The strong association of the human leukocyte antigen B∗27 alleles (HLA-B∗27) with spondyloarthritis and related rheumatic conditions has long fascinated researchers, yet the precise mechanisms underlying its pathogenicity remain elusive. Here, we review how interplay between the microbiome, the immune system, and the enigmatic HLA-B∗27 could trigger spondyloarthritis, with a focus on whether HLA-B∗27 presents an arthritogenic peptide. We propose mechanisms by which the unique biochemical characteristics of the HLA-B∗27 protein structure, particularly its peptide binding groove, could dictate its propensity to induce pathological T cell responses. We further provide new insights into how TRBV9+ CD8+ T cells are implicated in the disease process, as well as how the immunometabolism of T cells modulates tissue-specific inflammatory responses in spondyloarthritis. Finally, we present testable models and suggest approaches to this problem in future studies given recent advances in computational biology, chemical biology, structural biology, and small-molecule therapeutics.

Small-molecule probe for IBD risk variant GPR65 I231L alters cytokine signaling networks through positive allosteric modulation.

The proton-sensing heterotrimeric guanine nucleotide-binding protein-coupled receptor GPR65 is expressed in immune cells and regulates tissue homeostasis in response to decreased extracellular pH, which occurs in the context of inflammation and tumorigenesis. Genome-wide association studies linked GPR65 to several autoimmune and inflammatory diseases such as multiple sclerosis and inflammatory bowel disease (IBD). The loss-of-function GPR65 I231L IBD risk variant alters cellular metabolism, impairs protective tissue functions, and increases proinflammatory cytokine production. Hypothesizing that a small molecule designed to potentiate GPR65 at subphysiological pH could decrease inflammatory responses, we found positive allosteric modulators of GPR65 that engage and activate both human and mouse orthologs of the receptor. We observed that the chemical probe BRD5075 alters cytokine and chemokine programs in dendritic cells, establishing that immune signaling can be modulated by targeting GPR65. Our investigation offers improved chemical probes to further interrogate the biology of human GPR65 and its clinically relevant genetic variants.

Correction of age-associated defects in dendritic cells enables CD4+ T cells to eradicate tumors.

Defective host defenses later in life are associated with changes in immune cell activities, suggesting that age-specific considerations are needed in immunotherapy approaches. In this study, we found that PD-1 and CTLA4-based cancer immunotherapies are unable to eradicate tumors in elderly mice. This defect in anti-tumor activity correlated with two known age-associated immune defects: diminished abundance of systemic naive CD8+ T cells and weak migratory activities of dendritic cells (DCs). We identified a vaccine adjuvant, referred to as a DC hyperactivator, which corrects DC migratory defects in the elderly. Vaccines containing tumor antigens and DC hyperactivators induced T helper type 1 (TH1) CD4+ T cells with cytolytic activity that drive anti-tumor immunity in elderly mice. When administered early in life, DC hyperactivators were the only adjuvant identified that elicited anti-tumor CD4+ T cells that persisted into old age. These results raise the possibility of correcting age-associated immune defects through DC manipulation.

The centrosomal protein FGFR1OP controls myosin function in murine intestinal epithelial cells.

Recent advances in human genetics have shed light on the genetic factors contributing to inflammatory diseases, particularly Crohn's disease (CD), a prominent form of inflammatory bowel disease. Certain risk genes associated with CD directly influence cytokine biology and cell-specific communication networks. Current CD therapies primarily rely on anti-inflammatory drugs, which are inconsistently effective and lack strategies for promoting epithelial restoration and mucosal balance. To understand CD's underlying mechanisms, we investigated the link between CD and the FGFR1OP gene, which encodes a centrosome protein. FGFR1OP deletion in mouse intestinal epithelial cells disrupted crypt architecture, resulting in crypt loss, inflammation, and fatality. FGFR1OP insufficiency hindered epithelial resilience during colitis. FGFR1OP was crucial for preserving non-muscle myosin II activity, ensuring the integrity of the actomyosin cytoskeleton and crypt cell adhesion. This role of FGFR1OP suggests that its deficiency in genetically predisposed individuals may reduce epithelial renewal capacity, heightening susceptibility to inflammation and disease.

Discovery of disease-adapted bacterial lineages in inflammatory bowel diseases.

Gut bacteria are implicated in inflammatory bowel disease (IBD), but the strains driving these associations are unknown. Large-scale studies of microbiome evolution could reveal the imprint of disease on gut bacteria, thus pinpointing the strains and genes that may underlie inflammation. Here, we use stool metagenomes of thousands of IBD patients and healthy controls to reconstruct 140,000 strain genotypes, revealing hundreds of lineages enriched in IBD. We demonstrate that these strains are ancient, taxonomically diverse, and ubiquitous in humans. Moreover, disease-associated strains outcompete their healthy counterparts during inflammation, implying long-term adaptation to disease. Strain genetic differences map onto known axes of inflammation, including oxidative stress, nutrient biosynthesis, and immune evasion. Lastly, the loss of health-associated strains of Eggerthella lenta was predictive of fecal calprotectin, a biomarker of disease severity. Our work identifies reservoirs of strain diversity that may impact inflammatory disease and can be extended to other microbiome-associated diseases.

Efficient and accurate mixed model association tool for single-cell eQTL analysis.

Understanding the genetic basis of gene expression can help us understand the molecular underpinnings of human traits and disease. Expression quantitative trait locus (eQTL) mapping can help in studying this relationship but have been shown to be very cell-type specific, motivating the use of single-cell RNA sequencing and single-cell eQTLs to obtain a more granular view of genetic regulation. Current methods for single-cell eQTL mapping either rely on the "pseudobulk" approach and traditional pipelines for bulk transcriptomics or do not scale well to large datasets. Here, we propose SAIGE-QTL, a robust and scalable tool that can directly map eQTLs using single-cell profiles without needing aggregation at the pseudobulk level. Additionally, SAIGE-QTL allows for testing the effects of less frequent/rare genetic variation through set-based tests, which is traditionally excluded from eQTL mapping studies. We evaluate the performance of SAIGE-QTL on both real and simulated data and demonstrate the improved power for eQTL mapping over existing pipelines.

Oral and non-oral lichen planus show genetic heterogeneity and differential risk for autoimmune disease and oral cancer.

Lichen planus (LP) is a T-cell-mediated inflammatory disease affecting squamous epithelia in many parts of the body, most often the skin and oral mucosa. Cutaneous LP is usually transient and oral LP (OLP) is most often chronic, so we performed a large-scale genetic and epidemiological study of LP to address whether the oral and non-oral subgroups have shared or distinct underlying pathologies and their overlap with autoimmune disease. Using lifelong records covering diagnoses, procedures, and clinic identity from 473,580 individuals in the FinnGen study, genome-wide association analyses were conducted on carefully constructed subcategories of OLP (n = 3,323) and non-oral LP (n = 4,356) and on the combined group. We identified 15 genome-wide significant associations in FinnGen and an additional 12 when meta-analyzed with UKBB (27 independent associations at 25 distinct genomic locations), most of which are shared between oral and non-oral LP. Many associations coincide with known autoimmune disease loci, consistent with the epidemiologic enrichment of LP with hypothyroidism and other autoimmune diseases. Notably, a third of the FinnGen associations demonstrate significant differences between OLP and non-OLP. We also observed a 13.6-fold risk for tongue cancer and an elevated risk for other oral cancers in OLP, in agreement with earlier reports that connect LP with higher cancer incidence. In addition to a large-scale dissection of LP genetics and comorbidities, our study demonstrates the use of comprehensive, multidimensional health registry data to address outstanding clinical questions and reveal underlying biological mechanisms in common but understudied diseases.

Blended Genome Exome (BGE) as a Cost Efficient Alternative to Deep Whole Genomes or Arrays.

Genomic scientists have long been promised cheaper DNA sequencing, but deep whole genomes are still costly, especially when considered for large cohorts in population-level studies. More affordable options include microarrays + imputation, whole exome sequencing (WES), or low-pass whole genome sequencing (WGS) + imputation. WES + array + imputation has recently been shown to yield 99% of association signals detected by WGS. However, a method free from ascertainment biases of arrays or the need for merging different data types that still benefits from deeper exome coverage to enhance novel coding variant detection does not exist. We developed a new, combined, "Blended Genome Exome" (BGE) in which a whole genome library is generated, an aliquot of that genome is amplified by PCR, the exome regions are selected and enriched, and the genome and exome libraries are combined back into a single tube for sequencing (33% exome, 67% genome). This creates a single CRAM with a low-coverage whole genome (2-3x) combined with a higher coverage exome (30-40x). This BGE can be used for imputing common variants throughout the genome as well as for calling rare coding variants. We tested this new method and observed >99% r 2 concordance between imputed BGE data and existing 30x WGS data for exome and genome variants. BGE can serve as a useful and cost-efficient alternative sequencing product for genomic researchers, requiring ten-fold less sequencing compared to 30x WGS without the need for complicated harmonization of array and sequencing data.

An esophagus cell atlas reveals dynamic rewiring during active eosinophilic esophagitis and remission.

Coordinated cell interactions within the esophagus maintain homeostasis, and disruption can lead to eosinophilic esophagitis (EoE), a chronic inflammatory disease with poorly understood pathogenesis. We profile 421,312 individual cells from the esophageal mucosa of 7 healthy and 15 EoE participants, revealing 60 cell subsets and functional alterations in cell states, compositions, and interactions that highlight previously unclear features of EoE. Active disease displays enrichment of ALOX15+ macrophages, PRDM16+ dendritic cells expressing the EoE risk gene ATP10A, and cycling mast cells, with concomitant reduction of TH17 cells. Ligand-receptor expression uncovers eosinophil recruitment programs, increased fibroblast interactions in disease, and IL-9+IL-4+IL-13+ TH2 and endothelial cells as potential mast cell interactors. Resolution of inflammation-associated signatures includes mast and CD4+ TRM cell contraction and cell type-specific downregulation of eosinophil chemoattractant, growth, and survival factors. These cellular alterations in EoE and remission advance our understanding of eosinophilic inflammation and opportunities for therapeutic intervention.

Gut microbiome and metabolome profiling in Framingham heart study reveals cholesterol-metabolizing bacteria.

Accumulating evidence suggests that cardiovascular disease (CVD) is associated with an altered gut microbiome. Our understanding of the underlying mechanisms has been hindered by lack of matched multi-omic data with diagnostic biomarkers. To comprehensively profile gut microbiome contributions to CVD, we generated stool metagenomics and metabolomics from 1,429 Framingham Heart Study participants. We identified blood lipids and cardiovascular health measurements associated with microbiome and metabolome composition. Integrated analysis revealed microbial pathways implicated in CVD, including flavonoid, γ-butyrobetaine, and cholesterol metabolism. Species from the Oscillibacter genus were associated with decreased fecal and plasma cholesterol levels. Using functional prediction and in vitro characterization of multiple representative human gut Oscillibacter isolates, we uncovered conserved cholesterol-metabolizing capabilities, including glycosylation and dehydrogenation. These findings suggest that cholesterol metabolism is a broad property of phylogenetically diverse Oscillibacter spp., with potential benefits for lipid homeostasis and cardiovascular health.

Chemical tools to expand the ligandable proteome: diversity-oriented synthesis-based photoreactive stereoprobes.

Chemical proteomics enables the global assessment of small molecule-protein interactions in native biological systems and has emerged as a versatile approach for ligand discovery. The range of small molecules explored by chemical proteomics has, however, been limited. Here, we describe a diversity-oriented synthesis (DOS)-inspired library of stereochemically-defined compounds bearing diazirine and alkyne units for UV light-induced covalent modification and click chemistry enrichment of interacting proteins, respectively. We find that these 'photo-stereoprobes' interact in a stereoselective manner with hundreds of proteins from various structural and functional classes in human cells and demonstrate that these interactions can form the basis for high-throughput screening-compatible nanoBRET assays. Integrated phenotypic analysis and chemical proteomics identified photo-stereoprobes that modulate autophagy by engaging the mitochondrial serine protease CLPP. Our findings show the utility of photo-stereoprobes for expanding the ligandable proteome, furnishing target engagement assays, and discovering and characterizing bioactive small molecules by cell-based screening.

Integrated annotation prioritizes metabolites with bioactivity in inflammatory bowel disease.

Microbial biochemistry is central to the pathophysiology of inflammatory bowel diseases (IBD). Improved knowledge of microbial metabolites and their immunomodulatory roles is thus necessary for diagnosis and management. Here, we systematically analyzed the chemical, ecological, and epidemiological properties of ~82k metabolic features in 546 Integrative Human Microbiome Project (iHMP/HMP2) metabolomes, using a newly developed methodology for bioactive compound prioritization from microbial communities. This suggested >1000 metabolic features as potentially bioactive in IBD and associated ~43% of prevalent, unannotated features with at least one well-characterized metabolite, thereby providing initial information for further characterization of a significant portion of the fecal metabolome. Prioritized features included known IBD-linked chemical families such as bile acids and short-chain fatty acids, and less-explored bilirubin, polyamine, and vitamin derivatives, and other microbial products. One of these, nicotinamide riboside, reduced colitis scores in DSS-treated mice. The method, MACARRoN, is generalizable with the potential to improve microbial community characterization and provide therapeutic candidates.

Comprehensive Association Analyses of Extraintestinal Manifestations in Inflammatory Bowel Disease.

BACKGROUND & AIMS: Patients with inflammatory bowel disease (IBD) frequently develop extraintestinal manifestations (EIMs) that contribute substantially to morbidity. We assembled the largest multicohort data set to date to investigate the clinical, serologic, and genetic factors associated with EIM complications in IBD. METHODS: Data were available in 12,083 unrelated European ancestry IBD cases with presence or absence of EIMs (eg, ankylosing spondylitis [ankylosing spondylitis and sacroiliitis], primary sclerosing cholangitis [PSC], peripheral arthritis, and skin and ocular manifestations) across 4 cohorts (Cedars-Sinai Medical Center, National Institute for Diabetes and Digestive and Kidney Diseases IBD Genetics Consortium, Sinai Helmsley Alliance for Research Excellence Consortium, and Risk Stratification and Identification of Immunogenetic and Microbial Markers of Rapid Disease Progression in Children with Crohn's Disease cohort). Clinical and serologic parameters were analyzed by means of univariable and multivariable regression analyses using a mixed-effects model. Within-case logistic regression was performed to assess genetic associations. RESULTS: Most EIMs occurred more commonly in female subjects (overall EIM: P = 9.0E-05, odds ratio [OR], 1.2; 95% CI, 1.1-1.4), with CD (especially colonic disease location; P = 9.8E-09, OR, 1.7; 95% CI, 1.4-2.0), and in subjects who required surgery (both CD and UC; P = 3.6E-19, OR, 1.7; 95% CI, 1.5-1.9). Smoking increased risk of EIMs except for PSC, where there was a "protective" effect. Multiple serologic associations were observed, including with PSC (anti-nuclear cytoplasmic antibody; IgG and IgA, anti-Saccharomyces cerevisiae antibodies; and anti-flagellin) and any EIM (anti-nuclear cytoplasmic antibody; IgG and IgA, anti-Saccharomyces cerevisiae antibodies; and anti-Pseudomonas fluorescens-associated sequence). We identified genome-wide significant associations within major histocompatibility complex (ankylosing spondylitis and sacroiliitis, P = 1.4E-15; OR, 2.5; 95% CI, 2.0-3.1; PSC, P = 2.7E-10; OR, 2.8; 95% CI, 2.0-3.8; ocular, P = 2E-08, OR, 3.6; 95% CI, 2.3-5.6; and overall EIM, P = 8.4E-09; OR, 2.2; 95% CI, 1.7-2.9) and CPEB4 (skin, P = 2.7E-08; OR, 1.5; 95% CI, 1.3-1.8). Genetic associations implicated tumor necrosis factor, JAK-STAT, and IL6 as potential targets for EIMs. Contrary to previous reports, only 2% of our subjects had multiple EIMs and most co-occurrences were negatively correlated. CONCLUSIONS: We have identified demographic, clinical, and genetic associations with EIMs that revealed underlying mechanisms and implicated novel and existing drug targets-important steps toward a more personalized approach to IBD management.

Absence of a pancreatic microbiome in intraductal papillary mucinous neoplasm.

OBJECTIVE: This study aims to validate the existence of a microbiome within intraductal papillary mucinous neoplasm (IPMN) that can be differentiated from the taxonomically diverse DNA background of next-generation sequencing procedures. DESIGN: We generated 16S rRNA amplicon sequencing data to analyse 338 cyst fluid samples from 190 patients and 19 negative controls, the latter collected directly from sterile syringes in the operating room. A subset of samples (n=20) and blanks (n=5) were spiked with known concentrations of bacterial cells alien to the human microbiome to infer absolute abundances of microbial traces. All cyst fluid samples were obtained intraoperatively and included IPMNs with various degrees of dysplasia as well as other cystic neoplasms. Follow-up culturing experiments were conducted to assess bacterial growth for microbiologically significant signals. RESULTS: Microbiome signatures of cyst fluid samples were inseparable from those of negative controls, with no difference in taxonomic diversity, and microbial community composition. In a patient subgroup that had recently undergone invasive procedures, a bacterial signal was evident. This outlier signal was not characterised by higher taxonomic diversity but by an increased dominance index of a gut-associated microbe, leading to lower taxonomic evenness compared with the background signal. CONCLUSION: The 'microbiome' of IPMNs and other pancreatic cystic neoplasms does not deviate from the background signature of negative controls, supporting the concept of a sterile environment. Outlier signals may appear in a small fraction of patients following recent invasive endoscopic procedures. No associations between microbial patterns and clinical or cyst parameters were apparent.

Translational genetics identifies a phosphorylation switch in CARD9 required for innate inflammatory responses.

Population genetics continues to identify genetic variants associated with diseases of the immune system and offers a unique opportunity to discover mechanisms of immune regulation. Multiple genetic variants linked to severe fungal infections and autoimmunity are associated with caspase recruitment domain-containing protein 9 (CARD9). We leverage the CARD9 R101C missense variant to uncover a biochemical mechanism of CARD9 activation essential for antifungal responses. We demonstrate that R101C disrupts a critical signaling switch whereby phosphorylation of S104 releases CARD9 from an autoinhibited state to promote inflammatory responses in myeloid cells. Furthermore, we show that CARD9 R101C exerts dynamic effects on the skin cellular contexture during fungal infection, corrupting inflammatory signaling and cell-cell communication circuits. Card9 R101C mice fail to control dermatophyte infection in the skin, resulting in high fungal burden, yet show minimal signs of inflammation. Together, we demonstrate how translational genetics reveals molecular and cellular mechanisms of innate immune regulation.

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.