Markov field network model of multi-modal data predicts effects of immune system perturbations on intravenous BCG vaccination in macaques.

Analysis of multi-modal datasets can identify multi-scale interactions underlying biological systems but can be beset by spurious connections due to indirect impacts propagating through an unmapped biological network. For example, studies in macaques have shown that Bacillus Calmette-Guerin (BCG) vaccination by an intravenous route protects against tuberculosis, correlating with changes across various immune data modes. To eliminate spurious correlations and identify critical immune interactions in a public multi-modal dataset (systems serology, cytokines, and cytometry) of vaccinated macaques, we applied Markov fields (MFs), a data-driven approach that explains vaccine efficacy and immune correlations via multivariate network paths, without requiring large numbers of samples (i.e., macaques) relative to multivariate features. We find that integrating multiple data modes with MFs helps remove spurious connections. Finally, we used the MF to predict outcomes of perturbations at various immune nodes, including an experimentally validated B cell depletion that induced network-wide shifts without reducing vaccine protection.