Cytokine and molecular networks in sepsis cases: a network biology approach

Background: Sepsis is a life-threatening condition of organ dysfunction caused by a dysregulated host immune response to infection. We performed network analysis of cytokine molecules and compared network structures between a systematic inflammatory response syndrome (SIRS) or normal control (NC) group and a sepsis group. Results: We recruited SIRS (n = 33) and sepsis (n = 89) patients from electronic medical records (EMR) according to whether data on PCT, CRP, interleukin (IL)-1β, IL-2, IL-4, IL-5, IL-6, IL-9, IL-10, IL-12p70, IL-13, IL-17, IL-22, TNF-α, and IFN-γ levels were available. From the public GEO dataset, GSE66099, GSE9960, GSE95233, GSE57065 were downloaded. Genes corresponding to 15 molecules were extracted from an expression array. A correlation matrix was formed for the 15 molecules and statistically significant molecular pairs were used as pairs for network analysis of coexpression. The number of molecular or gene expression pairs significantly correlated among the SIRS or control and sepsis groups are as follows for datasets: EMR, 15 and 15; GEO66099-1, 13 and 15; GEO9960, 13 and 11; GSE95233, 13 and 8; GSE66099-2, 15 and 14; GSE57065, 14 and 13, respectively. Network analysis revealed that network diameter, number of nodes and shortest path were equal to or lower in the sepsis group. Conclusions: The coexpression network in sepsis patients was relatively small sized and had lower shortest paths compared with the SIRS group or healthy control group. Cytokines with one degree (k = 1) are increased in sepsis group compared with SIRS or healthy control group. IL-9 and IL-2 were not included in network of sepsis group indicating that these cytokines showed no correlation with other cytokines. These data might imply that cytokines tend to be dysregulated in the sepsis group compared to that of SIRS or normal control groups.