Disease progression in autoimmune conditions such as systemic lupus erythematosus (SLE) is highly heterogeneous, and the cellular and molecular mechanisms driving disease-onset dynamics remain incompletely understood. Here, based on single-cell transcriptomics data on lupus-prone NZB/W F1 mice, we derived a mathematical cell-cell interaction model recapitulating early dynamics of innate immune cells in lupus nephritis. We identified a diverse landscape of tissue-associated ILC and vessel-associated NK cell populations. We conceived a scalable mathematical framework for analysis of immune-cell interaction dynamics. A specific model formulation considers ILC as amplifiers of inflammatory processes in the presence of autoantibodies in lupus-prone individuals. Systematic model analyses highlight the impact of positive feedback loops and spontaneous inflammatory events or environmental stimuli, and the timing-dependent effectiveness of depletion therapies. Additionally, our model links the critical role of ILC populations to hallmarks of SLE such as highly heterogeneous disease dynamics. Overall, our findings lay the groundwork towards a mathematical model of immune-tissue cellular crosstalk, enabling quantification of disease severity and prediction of responses to biologic treatments in autoimmune diseases.
Kreider et al. (Thu,) studied this question.