Immune Gene Expression Covaries with Gut Microbiome Composition in Stickleback

Commensal microbial communities have immense effects on their vertebrate hosts, contributing to a number of physiological functions, as well as host fitness. In particular, host immunity is strongly linked to microbiota composition through poorly understood bi-directional links. Gene expression may be a potential mediator of these links between microbial communities and host function. However, few studies have investigated connections between microbiota composition and expression of host immune genes in complex systems. Here, we leverage a large study of laboratory-raised fish from the species Gasterosteus aculeatus (three-spined stickleback) to document correlations between gene expression and microbiome composition. First, we examined correlations between microbiome alpha diversity and gene expression. Our results demonstrate robust positive associations between microbial alpha diversity and expression of host immune genes. Next, we examined correlations between host gene expression and abundance of microbial taxa. We identified 15 microbial families that were highly correlated with host gene expression. These families were all tightly correlated with host expression of immune genes and processes, falling into one of three categories-those positively correlated, negatively correlated, and neutrally related to immune processes. Furthermore, we highlight several important immune processes that are commonly associated with the abundance of these taxa, including both macrophage and B cell functions. Further functional characterization of microbial taxa will help disentangle the mechanisms of the correlations described here. In sum, our study supports prevailing hypotheses of intimate links between host immunity and gut microbiome composition. IMPORTANCE Here, we document associations between host gene expression and gut microbiome composition in a nonmammalian vertebrate species. We highlight associations between expression of immune genes and both microbiome diversity and abundance of specific microbial taxa. These findings support other findings from model systems which have suggested that gut microbiome composition and host immunity are intimately linked. Furthermore, we demonstrate that these correlations are truly systemic; the gene expression detailed here was collected from an important fish immune organ (the head kidney) that is anatomically distant from the gut. This emphasizes the systemic impact of connections between gut microbiota and host immune function. Our work is a significant advancement in the understanding of immune-microbiome links in nonmodel, natural systems.