Integrated omics analysis reveals the immunologic characteristics of cystic Peyer’s patches in the cecum of Bactrian camels

Bactrian camels have specific mucosa-associated lymphoid tissue (MALT) throughout the large intestine, with species-unique cystic Peyer’s patches (PPS) as the main type of tissue. However, detailed information about the molecular characteristics of PPS remains unclear. This study applied a transcriptomic analysis, untargeted metabolomics, and 16S rDNA sequencing to compare the significant differences between PPS and the adjacent normal intestine tissues (NPPS) during the healthy stage of three young Bactrian camels. The results showed that samples from PPS could be easily differentiated from the NPPS samples based on gene expression profile, metabolites, and microbial composition, separately indicated using dimension reduction methods. A total of 7,568 up-regulated and 1,266 down-regulated differentially expressed genes (DEGs) were detected, and an enrichment analysis found 994 DEGs that participated in immune-related functions, and a co-occurance network analysis identified nine hub genes (BTK, P2RX7, Pax5, DSG1, PTPN2, DOCK11, TBX21, IL10, and HLA-DOB) during multiple immunologic processes. Further, PPS and NPPS both had a similar pattern of most compounds among all profiles of metabolites, and only 113 differentially expressed metabolites (DEMs) were identified, with 101 of these being down-regulated. Deoxycholic acid (DCA; VIP = 37.96, log2FC = −2.97, P = 0), cholic acid (CA; VIP = 13.10, log2FC = −2.10, P = 0.01), and lithocholic acid (LCA; VIP = 12.94, log2FC = −1.63, P = 0.01) were the highest contributors to the significant dissimilarities between groups. PPS had significantly lower species richness (Chao1), while Firmicutes (35.92% ± 19.39%), Bacteroidetes (31.73% ± 6.24%), and Proteobacteria (13.96% ± 16.21%) were the main phyla across all samples. The LEfSe analysis showed that Lysinibacillus, Rikenellaceae_RC9_gut_group, Candidatus_Stoquefichus, Mailhella, Alistipes, and Ruminococcaceae_UCG_005 were biomarkers of the NPPS group, while Escherichia_Shigella, Synergistes, Pyramidobacter, Odoribacter, Methanobrevibacter, Cloacibacillus, Fusobacterium, and Parabacteroides were significantly higher in the PPS group. In the Procrustes analysis, the transcriptome changes between groups showed no significant correlations with metabolites or microbial communities, whereas the alteration of metabolites significantly correlated with the alteration of the microbial community. In the co-occurrence network, seven DEMs (M403T65-neg, M329T119-neg, M309T38-neg, M277T42-2-neg, M473T27-neg, M747T38-1-pos, and M482t187-pos) and 14 genera (e.g., Akkermansia, Candidatus-Stoquefichus, Caproiciproducens, and Erysipelatoclostridium) clustered much more tightly, suggesting dense interactions. The results of this study provide new insights into the understanding of the immune microenvironment of the cystic PPS in the cecum of Bactrian camels.