Abstract
Salmonella negatively impacts the poultry industry and threatens animals' and humans' health. The gastrointestinal microbiota and its metabolites can modulate the host's physiology and immune system. Recent research demonstrated the role of commensal bacteria and short-chain fatty acids (SCFAs) in developing resistance to Salmonella infection and colonization. However, the complex interactions among chicken, Salmonella, host-microbiome, and microbial metabolites remain unelucidated. Therefore, this study aimed to explore these complex interactions by identifying the driver and hub genes highly correlated with factors that confer resistance to Salmonella. Differential gene expression (DEGs) and dynamic developmental genes (DDGs) analyses and weighted gene co-expression network analysis (WGCNA) were performed using transcriptome data from the cecum of Salmonella Enteritidis-infected chicken at 7 and 21 days after infection. Furthermore, we identified the driver and hub genes associated with important traits such as the heterophil/lymphocyte (H/L) ratio, body weight post-infection, bacterial load, propionate and valerate cecal contents, and Firmicutes, Bacteroidetes, and Proteobacteria cecal relative abundance. Among the multiple genes detected in this study, EXFABP, S100A9/12, CEMIP, FKBP5, MAVS, FAM168B, HESX1, EMC6, and others were found as potential candidate gene and transcript (co-) factors for resistance to Salmonella infection. In addition, we found that the PPAR and oxidative phosphorylation (OXPHOS) metabolic pathways were also involved in the host's immune response/defense against Salmonella colonization at the earlier and later stage post-infection, respectively. This study provides a valuable resource of transcriptome profiles from chicken cecum at the earlier and later stage post-infection and mechanistic understanding of the complex interactions among chicken, Salmonella, host-microbiome, and associated metabolites.