Fig. 3
IgA-production-related pathways are upregulated in the duodenum of WLMB. A Sample correlation heatmap based on the Pearson’s correlation coefficients of all gene expression between every two samples. B PCA plot of gene expression between AA and WLMB. C Volcano plot of genes, Q value < 0.05 and |log2(fold change)| > 0.5850 as the threshold parameters. D Heatmap of DEGs in KEGG pathways related to IgA production, based on fragments per kilobase per million (FPKM). E Enrichment plots of DEGs from the gene set enrichment analysis (GSEA). NES, normalized enrichment score