Cross-Talk Between Intestinal Microbiota and Host Gene Expression in Gilthead Sea Bream (Sparus aurata) Juveniles: Insights in Fish Feeds for Increased Circularity and Resource Utilization

© 2021 Naya-Català, do Vale Pereira, Piazzon, Fernandes, Calduch-Giner, Sitjà-Bobadilla, Conceição and Pérez-Sánchez.New types of fish feed based on processed animal proteins (PAPs), insect meal, yeast, and microbial biomasses have been used with success in gilthead sea bream. However, some drawback effects on feed conversion and inflammatory systemic markers were reported in different degrees with PAP- and non-PAP-based feed formulations. Here, we focused on the effects of control and two experimental diets on gut mucosal-adherent microbiota, and how it correlated with host transcriptomics at the local (intestine) and systemic (liver and head kidney) levels. The use of tissue-specific PCR-arrays of 93 genes in total rendered 13, 12, and 9 differentially expressed (DE) genes in the intestine, liver, and head kidney, respectively. Illumina sequencing of gut microbiota yielded a mean of 125,350 reads per sample, assigned to 1,281 operational taxonomic unit (OTUs). Bacterial richness and alpha diversity were lower in fish fed with the PAP diet, and discriminant analysis displayed 135 OTUs driving the separation between groups with 43 taxa correlating with 27 DE genes. The highest expression of intestinal pcna and alpi was achieved in PAP fish with intermediate values in non-PAP, being the pro-inflammatory action of alpi associated with the presence of Psychrobacter piscatorii. The intestinal muc13 gene was down-regulated in non-PAP fish, with this gene being negatively correlated with anaerobic (Chloroflexi and Anoxybacillus) and metal-reducing (Pelosinus and Psychrosinus) bacteria. Other inflammatory markers (igm, il8, tnfα) were up-regulated in PAP fish, positively correlating the intestinal igm gene with the inflammasome activator Escherichia/Shigella, whereas the systemic expression of il8 and tnfα was negatively correlated with the Bacilli class in PAP fish and positively correlated with Paracoccus yeei in non-PAP fish. Overall changes in the expression pattern of il10, galectins (lgals1, lgals8), and toll-like receptors (tlr2, tlr5, tlr9) reinforced the anti-inflammatory profile of fish fed with the non-PAP diet, with these gene markers being associated with a wide range of OTUs. A gut microbiota-liver axis was also established, linking the microbial generation of short chain fatty acids with the fueling of scd1- and elovl6-mediated lipogenesis. In summary, by correlating the microbiome with host gene expression, we offer new insights in the evaluation of fish diets promoting gut and metabolism homeostasis, and ultimately, the health of farmed fish.This study was supported by the EU H2020 GAIN Project (Aquaculture intensification in Europe, contract 773330). This study reflects only the views of the authors, and the European Union cannot be held responsible for any use that may be made of the information contained herein. Additional funding was obtained by the EU H2020 Research Innovation Program under the TNA Program (project AE150004) at the IATS-CSIC Research Infrastructure within the AQUAEXCEL2020 Project (652831), and the Spanish MICINN Project (Bream-AquaINTECH, RTI2018–094128-B-I00). MP was funded by a Ramón y Cajal Postdoctoral Research Fellowship (RYC2018-024049-I/AEI/10.13039/501100011033, co-funded by the European Social Fund (ESF) and ACOND/2020 Generalitat Valenciana)