Multi-level aggregation analysis of microbiome composition and host gene expression reveals associations with systemic and local immunity

Abstract

ABSTRACT The human gut microbiome plays a critical role in immune regulation, yet the molecular links between microbiome composition and host gene expression remain incompletely understood. We analyzed associations between host gene expression and microbiome composition in a cohort of 315 healthy individuals, integrating microarray-based gene expression data from three intestinal sites (ileum, transverse colon, and rectum) and six immune cell types with microbiome sequencing data. Using a hierarchical feature aggregation strategy combining principal component analysis, clustering, and covariate correction, we discovered significant associations primarily related to immunity. While microbial profiles were similar across the three intestinal sites, the transverse colon yielded the most “microbiome-host gene expression” associations. Among the immune cell types, CD8+ cells showed the highest number of associations. The first principal component of microbiome composition, reflecting a gradient from commensals (e.g., Ruminococcaceae and Christensenellaceae) to proinflammatory taxa ([Ruminococcus] gnavus and Lachnoclostridium), correlated with the expression of TNF-α-linked genes (HMOX1, CPI17, HSD3B2, and SLC5A1). Among individual genera, Catenibacterium abundance was associated with gene expression in both intestinal and immune cells, including negative associations with MRPS21 (related to mitochondrial function) in the transverse colon and with CD8+ gene programs related to T cell differentiation. These findings align with emerging evidence implicating mitochondrial dysfunction in intestinal inflammation. Our results identify multi-level associations between the gut microbiome and host gene expression, suggesting potential mechanisms by which microbiota shape local and systemic immunity and vice versa. The implicated genes and taxa represent candidates for experimental validation to improve understanding of host-microbiome homeostasis and its disruption in disease. IMPORTANCE The gut microbiome and immune system are engaged in a complex interplay throughout human life. While most associative studies focus on case-control comparisons—typically examining patients with conditions such as inflammatory bowel disease or metabolic diseases—less is known about the molecular links between the microbiome and immune system in healthy individuals. In this study of a large cohort of healthy individuals, we addressed this gap by applying multiscale modeling to tackle the high dimensionality of host-microbiome data. We identified multi-level associations between microbiome composition and host gene expression in both intestinal tissues and immune cells. These findings offer a valuable reference for understanding baseline host-microbiome communication and highlight molecular candidates—such as TNF-α-related genes and mitochondrial pathways—for future experimental validation. The gut microbiome and immune system are engaged in a complex interplay throughout human life. While most associative studies focus on case-control comparisons—typically examining patients with conditions such as inflammatory bowel disease or metabolic diseases—less is known about the molecular links between the microbiome and immune system in healthy individuals. In this study of a large cohort of healthy individuals, we addressed this gap by applying multiscale modeling to tackle the high dimensionality of host-microbiome data. We identified multi-level associations between microbiome composition and host gene expression in both intestinal tissues and immune cells. These findings offer a valuable reference for understanding baseline host-microbiome communication and highlight molecular candidates—such as TNF-α-related genes and mitochondrial pathways—for future experimental validation.