BackgroundThe implantation of mechanical circulatory support devices in heart failure patients is associated with a systemic inflammatory response, potentially leading to death from multiple organ dysfunction syndrome. Previous studies point to the involvement of many mechanisms, but an integrative hypothesis does not yet exist. Using time-dependent whole-genome mRNA expression in circulating leukocytes, we constructed a systems-model to improve mechanistic understanding and prediction of adverse outcomes.MethodsWe sampled peripheral blood mononuclear cells from 22 consecutive patients undergoing mechanical circulatory support device (MCS) surgery, at 5 timepoints: day −1 preoperative, and postoperative days 1, 3, 5, and 8. Clinical phenotyping was performed using 12 clinical parameters, 2 organ dysfunction scoring systems, and survival outcomes. We constructed a strictly phenotype-driven time-dependent non-supervised systems-representation using weighted gene co-expression network analysis, and annotated eigengenes using gene ontology, pathway, and transcription factor binding site enrichment analyses. Genes and eigengenes were mapped to the clinical phenotype using a linear mixed-effect model, with Cox models also fit at each timepoint to survival outcomes.ResultsWe inferred a 19-module network, in which most module eigengenes correlated with at least one aspect of the clinical phenotype. We observed a response of advanced heart failure patients to surgery orchestrated into stages: first, activation of the innate immune response, followed by anti-inflammation, and finally reparative processes such as mitosis, coagulation, and apoptosis. Eigengenes related to red blood cell production and extracellular matrix degradation became predictors of survival late in the timecourse corresponding to multiorgan dysfunction and disseminated intravascular coagulation.ConclusionsOur model provides an integrative representation of leukocyte biology during the systemic inflammatory response following MCS device implantation. It demonstrates consistency with previous hypotheses, identifying a number of known mechanisms. At the same time, it suggests novel hypotheses about time-specific targets.Electronic supplementary materialThe online version of this article (doi:10.1186/s12920-017-0288-8) contains supplementary material, which is available to authorized users.