Fig. 3

Computational modeling of the pyroptosis enrichment score (PES) and construction of immune infiltration co-expression network. A The venn diagram of pyroptosis related genes identified in three databases (GO database, Reactome database, and previous studies). B The boxplot showed the comparison of pyroptosis related genes between CeD and HC tissues in the derivation cohort. C The boxplot showed the differences in the PES between CeD and HC tissues in the derivation cohort. D The boxplot showed the differences in the PES between CeD and HC tissues in the validation cohort. E The boxplot showed the difference comparison of PES among three groups in dataset GSE125625. F The boxplot showed the difference comparison of PES among three groups from dataset GES191015. G The stacking histogram exhibited the infiltration proportion of different lymphoid cell types between CeD and HC in the derivation cohort. H Difference comparison of the infiltration score of various lymphoid cell types between CeD and HC in the derivation cohort. I The stacking histogram exhibited the infiltration proportion of different myeloid cell types between CeD and HC in the derivation cohort. J Difference comparison of the infiltrated score of various myeloid cell types between CeD and HC in the derivation cohort. K The correlation heatmap reflected co-expression patterns among immune cell infiltration score, PES, and response to IFN-γ pathway in the derivation cohort. L Establishment of the celiac disease rat model. M Representative immunofluorescence images showed the distribution of γδT cells in duodenal tissues of CeD rats and the statistical analysis of γδT cells infiltration degree quantified by the proportion of CD3⁺ γδTCR⁺ IELs to total IECs