Bifurcation Analysis of Cancer-Immunity Cycle

Despite advances in cancer immunotherapy, many patients still develop immune escape and disease progression. Understanding the dynamical mechanisms underlying the transition between tumor control and immune escape is therefore critical. Bifurcation theory offers valuable insights into how gradual changes in system parameters can lead to sudden qualitative shifts in dynamic behavior. In biological contexts such as cancer-immune interactions, these shifts can mark critical thresholds between tumor control by immunotherapy and uncontrollable progression. This study reviews bifurcation analysis for systems of difference equations, and demonstrates an application in immunotherapy of cancer. We apply the method of eigenvalues and Jacobian matrices to a previously developed model of the cancer-immunity cycle. This model examines how varying levels of immune suppression (related to PD-1-mediated immune evasion) affect system behavior. Along with a saddle-node bifurcation which is relatively more straightforward, the model reveals a stability bifurcation which depicts a transient stable equilibrium that appears (immune-limited state) and then disappears (immune-escape state). Since the potential to control cancer to a finite size depends on the existence of a stable equilibrium, understanding how to adjust immune reaction could be a game changer in immunotherapy. This framework therefore provides a mathematical basis for adaptive immunotherapy, in which treatment intensity is adjusted dynamically in response to evolving tumor-immune interactions. The stability bifurcation does not always coincide with the saddle-node bifurcation, instead it can occur at a lower immune suppression value (called an "early bifurcation"). Our case studies reveals that an early bifurcation arises only under a strong immune condition, whereas both bifurcations occur coincidentally under a weak immune condition. This discrepancy implies that relying only on the criterion with the saddle-node bifurcation may lead to an overestimate an immunotherapy effect under a strong immune environment. These results provide insight into the treatment strategy of anti-PD-1 based medicine. Along with the analysis of the immune suppression parameter, additional parameter analyses highlight that both the existence and stability of equilibria are highly sensitive near stability classification boundaries. Overall, this refined analysis strengthens the mathematical basis for modeling tumor-immune dynamics and may help optimize immuno-therapeutic strategies.