Fig. 2

Differentially expressed genes (DEGs) analysis and gene functional enrichment analysis between CeD and HC tissues. A Volcano map of the top 20 upregulated genes and top 20 downregulated genes. B Heatmap of a total of 137 differentially expressed genes with unsupervised clustering. C The bar diagram showed Gene Ontology enrichment analysis of DEGs based on Z-score. “BP” referred to “Biological Process,” “CC” referred to “Cellular Component” and “MF” referred to “Molecular Function.” D The bubble plot exhibited KEGG pathway analysis of DEGs with 10 most activated and 10 most inhibited pathways. E The trigram array plot demonstrated the Reactome pathway analysis of DEGs and the top 20 enrichment pathways