Fig. 5

Identification of hub genes in neutrophil-associated DEGs. A Venn diagram representing the intersection genes taken between the differential gene set of MMD and neutrophil-associated gene set. A total of 30 intersection genes were identified. B Neutrophil-associated DEGs were extracted, and drawn as a heatmap. The heat map row names are the gene symbols intersected by the list of DEGs and neutrophil-associated genes. C Neutrophil-associated DEGs were subjected to GO functional enrichment analysis. The diagram shows the network between enriched GO terms, each dot being a pathway of enriched biological processes. The dot color indicates the terms of the cluster, and the most significant terms in the cluster are identified in bold font with color. D PPI network relationship graph of neutrophil-associated DEGs predicted by STRING. Each dot represents the protein molecule ultimately transcribed by each DEG and is distinguished by a different color. Connecting lines between proteins indicate the presence of their interactions, and the line thickness indicates the strength of interactions. The bottom row of dots shows protein molecules for which no network relationships were identified. E Identification of hub genes from the PPI network using MCODE. Cluster sub-networks were identified using Cytoscape’s MCODE plugin, and the top two cluster subnetworks are shown. Nodes and edges indicate the cluster members and their interconnections. F Heat map of the correlation between hub genes and immune cell ratios. The X-axis indicates the expression of hub genes; Y-axis shows the immune cell. Pie chart size represents the correlation coefficient, blue being positive and red negative correlation, color saturation indicating the strength of the correlation as shown on the bar to the right