Cytokine modules defined in combined cytokine profiles of Il10^+/+ and Il10^−/− mice are associated with B-cell outcomes. (A) Heat map of cytokine modules displaying pairwise reliability scores of adjusted cytokines over 1000 random samplings of subjects. The optimal number of clusters was determined by the Tibshirani gap statistic method. (B) Table of associations between PM5 module and immune cell outcomes with differential abundance between Il10^+/+ and Il10^−/− bone marrow (BM) or peripheral blood (PB). Log-cytokine regression coefficient, P value, FWER, and false-discovery rate (FDR) were calculated as described.²⁹