The host immune response against infection requires the coordinated action of many diverse cell subsets that dynamically adapt to a pathogen threat. Due to the complexity of such a response, most immunological studies have focused on a few genes, proteins, or cell types. With the development of "omic"-technologies and computational analysis methods, attempts to analyze and understand complex system dynamics are now feasible. However, the decomposition of transcriptomic data sets generated from complete organs remains a major challenge. Here, we combined Weighted Gene Coexpression Network Analysis (WGCNA) and Digital Cell Quantifier (DCQ) to analyze time-resolved mouse splenic transcriptomes in acute and chronic Lymphocytic Choriomeningitis Virus (LCMV) infections. This enabled us to generate hypotheses about complex immune functioning after a virus-induced perturbation. This strategy was validated by successfully predicting several known immune phenomena, such as effector cytotoxic T lymphocyte (CTL) expansion and exhaustion. Furthermore, we predicted and subsequently verified experimentally macrophage-CD8 T cell cooperativity and the participation of virus-specific CD8+ T cells with an early effector transcriptome profile in the host adaptation to chronic infection. Thus, the linking of gene expression changes with immune cell kinetics provides novel insights into the complex immune processes within infected tissues.

This work is supported by a grant from the Spanish Ministry of Economy, Industry and Competitiveness and FEDER grant no. SAF2016-75505-R (AEI/MINEICO/FEDER, UE) and the María de Maeztu Programme for Units of Excellence in R&D (MDM-2014-0370). GB and AM are also supported by the Russian Science Foundation (grant 18-11-00171). AE-C and SH are also supported by Instituto de Salud Carlos III (ISCIII) grant from the Spanish Ministry of Economy, Industry and Competitiveness and FEDER grant no. PT17/0009/0019. IG-V is supported by European Research Council (637885)