 Research article
 Open Access
 Published:
Cancer immune control dynamics: a clinical data driven model of systemic immunity in patients with metastatic melanoma
BMC Bioinformatics volume 22, Article number: 197 (2021)
Abstract
Background
Recent clinical advances in cancer immunotherapeutics underscore the need for improved understanding of the complex relationship between cancer and the multiple, multifunctional, interdependent, cellular and humoral mediators/regulators of the human immune system. This interdisciplinary effort exploits engineering analysis methods utilized to investigate anomalous physical system behaviors to explore immune system behaviors. Cancer Immune Control Dynamics (CICD), a systems analysis approach, attempts to identify differences between systemic immune homeostasis of 27 healthy volunteers versus 14 patients with metastatic malignant melanoma based on daily serial measurements of conventional peripheral blood biomarkers (15 cell subsets, 35 cytokines). The modeling strategy applies engineering control theory to analyze an individual’s immune system based on the biomarkers’ dynamic nonlinear oscillatory behaviors. The reverse engineering analysis uses a Singular Value Decomposition (SVD) algorithm to solve the inverse problem and identify a solution profile of the active biomarker relationships. Herein, 28,605 biologically possible biomarker interactions are modeled by a set of matrix equations creating a system interaction model. CICD quantifies the model with a participant’s biomarker data then computationally solves it to measure each relationship’s activity allowing a visualization of the individual’s current state of immunity.
Results
CICD results provide initial evidence that this modelbased analysis is consistent with identified roles of biomarkers in systemic immunity of cancer patients versus that of healthy volunteers. The mathematical computations alone identified a plausible network of immune cells, including T cells, natural killer (NK) cells, monocytes, and dendritic cells (DC) with cytokines MCP1 [CXCL2], IP10 [CXCL10], and IL8 that play a role in sustaining the state of immunity in advanced cancer.
Conclusions
With CICD modeling capabilities, the complexity of the immune system is mathematically quantified through thousands of possible interactions between multiple biomarkers. Therefore, the overall state of an individual’s immune system regardless of clinical status, is modeled as reflected in their blood samples. It is anticipated that CICDbased capabilities will provide tools to specifically address cancer and treatment modulated (immune checkpoint inhibitors) parameters of human immunity, revealing clinically relevant biological interactions.
Background
In recent years, new insights into the state of the systemic immunity in cancer patients suggest a constellation of multiple abnormalities in the immune system that have the potential to directly impact not only the clinical response to immune checkpoint inhibitors (ICI) therapy but also the natural history of the malignant disease [1]. It is becoming increasingly clear that the panoply of multiple aberrancies in the immune system of cancer patients likely represent a manifestation of a complex set of biological processes that require an interrogative approach capable of complex systems analyses that take into account multiple variably interdependent parameters (biomarkers). The scientific community is increasingly recognizing that to analyze such complex systems, a challenging interdisciplinary approach to create meaningful biological computational tools is urgently needed [2,3,4,5,6,7]. A team of medical oncologists/immunologists and systems engineers have come together to overcome these challenges by applying wellestablished mathematical algorithms and engineering knowledge of physical nonlinear oscillations to biological oscillations of cells and cytokines in the immune system [8,9,10,11,12]. The product of this multiyear collaboration is the innovative biological mathematical tool, Cancer Immune Control Dynamics (CICD), a clinical data driven model of systemic immunity. CICD strives to bridge biology and engineering by generating an adaptable biological modelbased analysis program that incorporates multiple biomarkers. This endeavor attempts to understand and organize the individual roles of multiple measured, interacting, variably interdependent mediators (cellular and humoral) of systemic immune homeostasis (both up and down regulatory) through measurable biomarker concentrations in peripheral blood. The overall state of an individual’s immune system regardless of clinical status, is modeled as reflected in their blood samples.
Mathematical modeling concepts used in various engineering systems are the foundation of the CICD model presented herein. Physical system models simultaneously embed a multitude of complex and often illdefined representational features including system disturbances, environment, feedback relationships, modularity, adaptability, robustness, redundancy, etc. These same illdefined features are observed in the adaptive immune system that generate a dynamic response via a complex network of immune biomarker interactions that are heterogeneous, highly redundant and maintain homeostasis. Biological and engineering technologies appear very different but at the system level have characterizations that are mathematically equivalent and can potentially be modeled using the same principles [13,14,15]. To mathematically model the dynamics of the immune system, we assume that the oscillations of biomarkers’ concentrations will expose their interrelations, regardless of their disparate phenotypes and disparate functions (up vs. down regulatory). By mathematically determining which biomarkers are fluctuating together (and which are not), thereby possibly mediating each other via activation or suppression, the underlying biological complexity of the immune system is uncovered.
Mathematical modeling of biological systems, specifically the immune system has greatly increased in recent years [7]. CICD modeling utilizes the common ordinary differential equation (ODE) methods [16,17,18]. This approach both enables calibration against various data and is also computationally efficient making it the most widespread and flexible of models. However, its complexity grows with the addition of more equations and consequently an exponentially growing number of tobedefined parameter values. This size problem places limits on ODE model growth potential for predictive analysis applications. To date these models are effective for wellknown situations but do not fully include the ability to investigate the enormous number of possible immune interactions that can take place in the human system concurrently. CICD circumvents these limitations by applying the same ODE to each biomarker to represent all possible interactions and parameters that affect its rate of change. Patient immune parameter measurements (concentrations of peripheral blood biomarkers) populate these equations producing a mathematical representation of an individual’s immune system.
Once a system is modeled mathematically, reverse engineering strategies [19, 20] can decompose a physical system’s complex interactions to reveal otherwise hidden features, structure, and control principles. By exploiting these reverse engineering methodologies for immunology [8], the connectivity and underlying dynamics characterizations for biomarker interaction networks can potentially be exposed. CICD methodology, constructed on these principles, is able to quantify underlying biomarker relationships in the complex immune system network through an expandable, descriptive method of mathematical modeling based solely on an individual’s test period data record.
CICD results provide initial evidence that modelbased systems analysis leads to patient specific clinical insights pertaining to the biology of disease (antitumor immunity) in humans. Herein we present results of our initial effort to develop an interrogative approach that reverse engineers in silico a model of the human immune system in patients with metastatic melanoma (and healthy volunteers) utilizing peripheral blood derived measurements (timeseries) of immune function (biomarkers). This effort has adopted a datadriven approach that exploits engineering analysis methods utilized to investigate anomalous system behaviors. First verification of the computational process for this novel approach is presented in detail to validate the calculations used to identify and quantify biomarker pairings prevalent in clinical diagnoses. Next, confirmation of CICD modeling begins with the portrayal of the biological community’s accepted truths regarding homeostasis, redundancy, heterogeneity, and homogeneity.
CICD exposes the underlying causes consisting of the biomarker relationships that together affect the observed populations of the biomarkers in the blood, thereby providing an insight rich snapshot of an individual’s state of immune homeostasis. Its main objective is discovery, to help lead to improvements in the efficacy of existing immune therapies (patient selection; drug combinations) and insight into new therapeutic targets that may significantly reduce the time to discover new therapeutics capable of meaningful clinical impacts.
Methods
Engineering approach
CICD’s modeling methodology belongs to a specialized subdiscipline known as inverse problems, important mathematical problems in science and mathematics because they provide information about parameters that are not directly observed. In essence, an inverse problem in science consists of the process of calculating from a set of observations the causal factors that produced them. CICD’s fundamental assumption (Fig. 1a) is that a multidimensional cause and effect dynamic biomarker interaction network occurs between biomarker populations in the immune system. CICD uses the observed measurable output populations to mathematically calculate the activity of the underlying input interactions that are potential signals utilized to control the immune system. CICD takes advantage of the oscillatory nature of the changes in biomarker concentrations over time in order to resolve interparametric relationships that result in the maintenance of systemic immune homeostasis in cancer versus healthy.
CICD utilizes principles of control theory to examine the internal dynamics of the human immune system. As in mechanical systems, a sensor provides continuous feedback information to a controller to induce a response to maintain the stability (homeostasis) of the system. Prior to distinguishing biological embedded sensors and controllers, a model first needs to be developed through system identification which builds a mathematical model of the dynamical system from measured data. With a viable mathematical model, additional engineering principles and algorithms are applied to reveal the underlying causes or signals that provide information to the system [9]. CICD’s innovation is that it provides a quantitative math model representation for the dynamic interaction network of the immune system based on serially collected peripheral blood measurements. Figure 1b provides the foundational architecture of CICD’s analysis process. The system identification process builds the math model by (1) utilizing the serially collected peripheral blood samples (biomarker data) as measured data for patient biological status, (2) representing the biomarker population dynamics as mathematical relationships between measured data expressed as nonlinear ordinary differential Lotka–Volterra predator–prey equations, and (3) creating a biomarker interaction network, through a matrix generalization of these equations using a truncated Kolmogorov–Gabor polynomial to include all possible biomarker relationships as specified in the Knowledge Model. These components generate predator–prey equations for each biomarker studied that are combined into a system characterization matrix to provide a dynamic behavior model. This matrix is the foundation for the reverse engineering analysis. This math model is reverse engineered using methods of linear matrix computation [11], specifically Singular Value Decomposition (SVD) to approximate all biologically possible nonlinear cause and effect coupling mechanisms. The composite resultant data exposes the most active biomarker relationships or causes that associate with observed effects as measured in the blood data. This unique modal analysis solution process creates a database of information that can be mined for characterizing patient and multipatient cause and effect network interaction dynamics.
The development of CICD has evolved through collaboration between medical and engineering disciplines using extensive control theory experience in the computational modeling of illconditioned systems. The associated vision for the analysis framework for how CICD’s modelbased information should be formulated, collected and presented has been influenced by multiple contributing sources [10, 13, 16, 21,22,23,24,25].
CICD patient biological status
The data used consists of 35 cytokines (plasma concentrations): EGF, EOTAXIN, FGF2, FLT3L, FRACTALKINE, GCSF, GMCSF, IFNa2, IFNg, IL10, IL12p40, IL12p70, IL13, IL15, IL17A, IL1a, IL1b, IL1ra, IL2, IL3, IL4, IL5, IL6, IL7, IL8, IL9, IP10, MCP1, MCP3, MIP1a, MIP1b, TGFa, TNFa, TNFb, VEGF and 15 cell phenotypes (relative concentrations among peripheral blood mononuclear cells): CD11c+, CD11c/CD14+, CD11c/CD86+, CD11c/HLADR+, CD123/HLADR+, CD14/CD197+, CD16/CD56+, CD3+, CD3/CD4+, CD3/CD62L+, CD3/CD69+, CD3/CD8+, CD4/CD294+, CD4/TIM3+, and CD56+. Blood samples are obtained at approximately the same time of day Monday thru Friday over a 2week period in 14 patients with stage IV metastatic melanoma, not on active therapy as well as in 27 healthy volunteers. Cells and cytokines are collectively referred to as biomarkers. Computational analysis for system oscillatory dynamics builds upon the implicit assumption that data sampling is a continuous dynamic process for which characterization models are derived from available data records. Herein linear interpolation defines the biomarker rate of change, slope, as constant between blood draws providing an initial, simple, numerical approximate of the entire test period response. CICD analysis is valid at any instant between blood draws. Six equally timed sampling points are chosen between each draw to approximate dynamic biomarker actions executed between blood draw instances. The methods used to obtain the clinical data is found in Additional file 1.
CICD model of biomarker population dynamics
A mathematical expression that can model biomarkers’ fluctuations over time in relation to each other’s growth and decay is needed to determine the prominent relationships. The Lotka–Volterra predator–prey equations are firstorder nonlinear differential equations used to describe the dynamics of biological systems in which two species interact. To “translate” these mathematical equations into immune biology, each cell and cytokine is viewed as an individual species. The predator–prey approach merely provides a simple mathematical framework to organize the competing influences of the numerous immuneactivating and immunesuppressing sets of signals involving both cellular and cytokine biomarkers. CICD generalizes these established constitutive relationship equations to include multiple biomarkers to model their interacting populations. The state variables of immunology are biomarker population size observed in blood (measured concentrations). This modeling approach provides the ability to capture concurrently oscillatory, exponential growth and decay behaviors between interacting biomarker populations as they are continuously being used and replaced (concentrations going down and up, respectively). Each biomarker’s population change is expressed via an equation, which CICD then uses to build a large, generalized, biomarker interaction matrix model.
CICD model of biomarker interaction network
The fundamental CICD derived relationship consists of three components “Target”, “Source”, and “Modulator”, modeled as the population of the Target biomarker is stimulated or suppressed by the population of the Source biomarker that is either modulated by the population of the Modulator biomarker (bilinear) or is not modulated (linear) (Fig. 2). This relational concept provides a means to create a profile of the dynamic immune system by defining thousands of sets of potential immune signals (up and down regulatory). A major design goal for CICD is to identify and quantify only the biologically possible underlying networks of connectivity in the human immune system. To ensure the broad applicability of CICD’s analysis methodology, a systemsengineering approach is adopted where theoretically all combinations can be included while allowing the user an ability to reduce or expand the characterization model according to the needs of the study. Through quantification of the observed behaviors of these possible causes, CICD creates a descriptive mathematical model of what has occurred in contrast to predicative models that hypothesize what will occur.
CICD uses a template, referred to herein as the Knowledge Model, to generate a comprehensive list of biomarker relationships that may cause a change in a biomarker population. In order to maintain biological integrity, a template is incorporated to identify specific biologically possible combinations for the analysis according to the biomarker’s type (Fig. 2). The Knowledge Model template states which group, cells or cytokines, can populate each component. With these criteria, currently 7 types of relationships, the 35 cytokines and 15 cells studied generate 28,605 unique biomarker relationships. This allows the application of system identification concepts to generalize the predator–prey equations to a mathematical model that approximates specific biologically possible relationships according to the study’s needs and is consistent with wellestablished biological knowledge. The generalization provides the framework to include all possible causes, i.e. biomarker relationships that affect the biomarker population rate of change.
Once all relationships to be included are determined, the Kolmogorov–Gabor polynomial (Eq. 1) models all actions that drive cellcytokine interaction network dynamics. It provides a framework to express each possible population change cause as an atomic element in the equation, as a linear (nonmodulated source) or bilinear (modulated source) product of state variables, each with an associated unknown parameter. Based on CICD’s multidimensional cause and effect assumption, the rate of change of a Target biomarker is expressed as being equal to the algebraic sum of each atomic element for every CICD derived relationship as dictated by the Knowledge Model.
The generic format of the Kolmogorov–Gabor polynomial is:
Herein, each variable of state \(x_{1} , \ldots , x_{n}\) represents a biomarker population, (Sources and Modulators), the a’s are the unknowns of the atomic elements of population change, and \(Y\left( {x_{1} , \ldots , x_{n} } \right)\) is the resultant of all possible relationships as dictated by the Knowledge Model that impact one system biomarker (Target). The unknown (a’s) incorporate all known and unknown underlying mediating parameters. CICD truncates the polynomial (Eq. 1) at the bilinear term to maintain SVD solution accuracy consistent with the square root of machine precision as well as limit problem size. The polynomial models all biomarkers relationships by providing a framework to express the set of Target equations (herein, 50) in linear matrix form (Eq. 2). This matrix provides the characterization of the immune regulatory system interactions at any time instant during the test period.
CICD model of reverse engineering
The system dynamics analysis relationship is expressed in linear matrix form (Eq. 2) as:
Herein \(\left\{ {\left. {\frac{dx}{{dt}}} \right\}} \right.\) is the short (50 × 1) vector of biomarker state variables, \(\left\{ a \right\}\) is the long (28,605 × 1) vector of all linear and bilinear relationship unknowns and \(\left[ K \right]\) is the long thin (50 × 28,605) rectangular matrix of Kolmogorov polynomial state variable linear and bilinear scalar products. For all atomic element constitutive relationships implicit in equation (Eq. 1), the CICD characterization matrix \(\left[ K \right]\) is completely quantified by biomarker values. The participant’s biomarker data populates each of the Knowledge Model’s 28,605 relationships for every analysis time point, currently 48 for 10 days of serial blood collections. This step of data preparation prior to the SVD calculation generates 48 uniquely valued matrixes each containing 28,605 non–zero data points for every participant. This nontypical predator–prey equation format creates a linear matrix algebraic relationship between the time instant biomarker change vector \(\left\{ {\left. {\frac{dx}{{dt}}} \right\}} \right.\) and all Kolmogorov Polynomial unknowns. The net result is a SVDfriendly reformulation of the nonlinear predator–prey ODE dynamics problem into an identical linear ODE problem that has a welldefined matrix \(\left[ K \right]\) of time varying coefficients and a computable vector of unknowns, i.e. \(\left\{ a \right\}\) for every time instant for which \(\left\{ {\left. {\frac{dx}{{dt}}} \right\}} \right.\) is defined. By making the plausible linear interpolation hypothesis for all biomarker response signals inserted between measured data points, a welldefined constant value for \(\left\{ {\left. {\frac{dx}{{dt}}} \right\}} \right.\) is obtained. The second equality sign introduces the SVD (Eq. 3) relationship that equates matrix \(\left[ K \right]\) to the matrix triple product of Left Singular Vectors (LSV’s) in \(\left[ U \right]\), singular values in \(\left[ {\Sigma } \right]\) and Right Singular Vectors (RSV’s) in \(\left[ V \right]^{T}\). These orthogonal vectors decompose the time instant set of fully coupled biomarker system equations into an equally sized set of fully decoupled generalized coordinate system equations. Singular values are scale factors associated with all orthonormalized LSV’s and RSV’s.
From the mathematical viewpoint, CICD fundamental Eqs. (2) and (3) directly equates the blood sample input data to the CICD results output. Input data totally defines contents of \(\left\{ {\left. {\frac{dx}{{dt}}} \right\}} \right.\) and \(\left[ K \right]\). Computed results totally define the contents of \(\left\{ a \right.\}\), \(\left[ U \right], \left[ \Sigma \right]\) and \(\left[ V \right]^{T}\). CICD analysis ensures that the resultant data is accurate to machine precision by having a quality check incorporated into the program. The resultant data, the computed values for \(\left\{ a \right.\}\), \(\left[ U \right]\), \(\left[ \Sigma \right]\) and \(\left[ V \right]^{T}\) is used to back compute the input data to ensure both the original input data and recalculated input data are within expected machine precision bounds. This validation ensures the integrity and numerical stability of the computed results.
An SVD algorithm variant specialized to work with the very long thin (50 × 28,605) rectangular matrices is utilized by CICD to circumvent wellknown numerical computation problems associated with several nearly equal singular values [26]. Conceptually, the SVD iterative solution process creates a sequential product of similarity transformations that is stopped before an illconditioned computing step is executed. The net result is a minimum norm, least effort, solution that has machine solution accuracy and a LSV matrix \(\left[ U \right]\) that is not unit diagonal. This specialized SVD decomposition of the system characterization matrix \(\left[ K \right]\) is key to CICD’s ability to reverse engineer a patient’s test period record to quantify all active causes that result in a biomarker population change. CICD analysis uncovers and quantifies the active biomarker relationships specific to an individual’s immune network through the interaction strength values for all LSV’s and RSV’s elements.
CICD model of biomarker cause and effect interaction network
Once SVD calculates the resultant data structures, they must be interpreted to understand the biological interactions relative to the measured input data. Figure 3 provides a flow chart on how the large relational data structure is processed to interpret and extract biologically relevant insight.
CICD analysis decomposes a \(\left[ K \right]\) matrix by SVD at every analysis time point during the test period so that each LSV column of \(\left[ U \right]\) has an associated RSV row of \(\left[ V \right]^{T}\) that contains all possible relationships as defined by the Knowledge Model (Fig. 3a). This results in two additional dimensions to the current threedimensional CICD derived relationship. Therefore, at each time point, the SVD calculation generates a CICD relationship with 5 dimensional components and an associated value: (1) test period time instant, (2) LSV column, (3) Target biomarker, (4) Source biomarker, and (5) Modulator biomarker.
The SVD variant computes 50 LSVs corresponding directly to the number of biomarkers analyzed. To connect the math of Eq. (3) to biology it is necessary to label each specific LSV column of \(\left[ U \right]\) and RSV row of \(\left[ V \right]^{T}\). This effort is nontrivial. At every analysis time instant, a nametagging algorithm associates a unique biomarker name to each LSV according to its dominant biomarker component. This associates each LSV natural mode of behavior with a specific biomarker. Ambiguous cases exist per LSV matrix \(\left[ U \right]\) but are relatively few in which additional steps are utilized. This labeling method ensures that every LSV is associated with a unique biomarker nametag according its behavioral activity.
CICD methodology includes a small set of response threshold rules to reflect the biology response expectation that LSV solutions will have a clearly identified dominator component and a modest number of supporter components, while removing very small valued interactions. For the resultant data analysis, all LSV components with values above the threshold are considered active. Therefore, for each 5dimensional relationship if its value is greater than the threshold limit it is assumed active and counted as a single Occurrence Count (Fig. 3b).
Each CICD calculated unknown value is assumed to represent all parameters acting on its specific relationship. Herein we removed those relationships with an unknown value of zero, per the original equation in which each term is the product of the unknown and the state variables representing the total causal effect of the specific biomarker relationship. If this effect is equal to zero, then it is assumed not active and is not included in the current investigation.
A key data processing step of the result data is a series of (Σ) summations to expose the highly active biomarker networks. Figure 3c illustrates the multiple views of the data that are generated by combining CICD relationships with like components. Each level provides insight into an individual’s immune state via several resultant test period relationships that can be explored for possible clinical insight.
CICD via SVD, computes 50 Singular Values with an associated set of Left and Right Singular Vectors (LSV, RSV). Each LSV has 50 components, one for each biomarker. Each RSV has 28,605 components, one for each possible biomarker relationship defined by the Knowledge Model in the current study. For modal analysis, this generates 50 × 28,605 unique 4dimensional relationships, each with a unique Knowledge Model relationship and an associated LSV across the time period creating CICD’s fundamental data structure, the Quad relationship (Fig. 3d). Therefore, the Quad relationship contains a unique Modulator, Source, Target, and LSV combination with the number of Occurrence Counts over the test period.
The CICD analysis establishes a new unique measurement, Occurrence Count that quantifies the influence of these biomarkers relationships or signals in an individual’s immune system. The total Occurrence Count is the number of times each specific Quad is active during the test period. The maximum Occurrence Count for each Quad corresponds directly to the number of blood draws and number of time instances analyzed, herein 48 for the 10 days of serial peripheral blood sample collections. To aid current patient comparisons as well as potential future studies (reducing the number of days for serial peripheral blood collections), each Quad’s Occurrence Count is divided by the maximum count as dictated by the number of time points then multiplied by 100. This normalization provides a consistent maximum value independent of study length as well as a representation of the percentage of activity for the Quad across the study period. Therefore, active Quads have a range of Occurrence Counts from 0 to 100. The CICD resultant data analyzed herein is the normalized Occurrence Count value of each unique Quad and is the basis for all result figures presented. CICD analysis presented herein originates with the Quads that are combined by LSV and Targets to obtain the SourceModulator Doublets in which high variation was observed across individuals. The principal results illustrate the variability of the Doublet Occurrence Counts across patient and clinical groups.
Results
Patient peripheral blood derived serial biomarker data
A sample of the peripheral blood data is presented in Fig. 4. 50 biomarkers are analyzed for both cancer patients and healthy volunteers across the 2week test period (consecutive weekdays only). Biomarker values vary both across time and between individuals. The highly dynamic and fluctuating nature of the immune system is demonstrated through a sample of the CICD biomarker input data for an individual cancer patient and a healthy donor. These graphs illustrate the overall general similarities and the deviations found between the study groups. It is these variations that CICD seek to quantify via underlying unobservable interactions that manifest in the observed and measurable concentrations. Each participant’s graphs of cell and cytokines concentrations as shown in Fig. 4 are found in Additional file 2.
CICD data analysis
A verification of methodology is presented through the findings that only 133,000 Quads out of a possible 1.4 million are nonzero across all participants and only several thousand for an individual. The CICD final results consist of a list of active Quad relationships and their normalized Occurrence Count for each patient.
The total Occurrence Count for each Quad, Triplet, Doublet and Singlet is explored and compared between patients and patient groups to reveal active biomarker interactions. Table 1 provides a small sample of the thousands of active Quads for a cancer patient and a healthy donor, furthermore it illustrates the summation process as depicted in Fig. 3c for Triplets (Modulator, Source, and Target) and Doublets (Modulator and Source). Due to the large dataset only one Target, IL8 for the ModulatorSource Doublet, IP10/CD3 + is provided herein, but it has been observed throughout datasets that a consistent pattern of activity is repeated for Targets/LSVs in an active Doublet. Three levels of Target/LSV activity have been continuously observed throughout all Quad datasets (Fig. 5). First, if the Doublet is active it will have at a minimum an Occurrence Count for the Quad with the same cell in the Source, Target, and LSV (Quad: IP10/CD3+/CD3+/CD3 +). The Source cell stimulates/suppresses itself. The next level of observed activity is in the Dominator LSV, the Quad with the same cytokine in the Target and LSV (Quad: IP10/CD3+/IL8/IL8). The redundancy of the system is illustrated by the final level of activity, the Supporter LSVs, the Quads with a different LSV cytokine from the Target (Quad: IP10/CD3+/IL8/nonIL8 Cytokines). These additional Occurrence Counts are viewed as supporting activity for the Triplet, increasing the effects on the Target cytokine. It has been observed that if a specific Doublet is found active in a specific patient, the Occurrence Counts for the cell Target, Dominator LSV and Supporter LSV activity presents in a similar pattern. The supporting activity varies across the LSVs but is only observed if both the cell Target and the Dominator LSV Quad are active for the Triplet providing the justification for the summation process. This cumulative pattern of Occurrence Counts manifests in the total ModulatorSource Doublet Occurrence Count, the higher the sum value, the higher the level of supporting LSVs across all Targets, and therefore the higher its overall effects on all biomarkers Targets. Much of the CICD result data reflects observations seen in the biomarker data that cannot be currently quantified. For the two study participants provided in the paper, the IL8 input data (blood biomarker concentrations) in the cancer patient is higher and appears to be oscillating in the same pattern as IP10 compared to the presented healthy volunteer (Fig. 4). For the same two study participants, both similarities and differences for their Occurrence Counts are observed for these biomarkers in the specific sample of Quads provided in Table 1. Here IP10 is modulating CD3 + which is targeting CD3 + and IL8 across several biomarkers’ LSV. The cancer patient in this example has a much higher level of activity for this IP10/CD3+/IL8 Triplet and IP10/CD3 + Doublet generated from various supporting biomarker LSVs as compared to this specific healthy volunteer. Therefore, according to CICD analysis the relationship IP10/CD3 + is affecting multiple cell and cytokine biomarkers more intensely in this cancer patient as compared to the healthy volunteer. Also, IL8 is acting as a nonmodulator source on several cellular targets but only for the cancer patient. These specific Knowledge Model relationships are now quantified by CICD result data. The Triplets and Doublets compared here are only a few of the thousands of observable comparisons found between the study participants and clinical groups (cancer vs. healthy) with CICD, illustrating that variations observed in the biomarker blood data can have many underlying causes that can be mathematically uncovered. This summation process is repeated for all Quads and additional software has been developed specifically tailored to explore CICD data both within and across patients.
CICD modeling of biologic interactions
CICD begins to reveal the subtle differences that are masked by the commonalities in all donors (cancer and healthy). A main goal of CICD is to uncover the specific networks of active interactions in cancer patients that contribute to their clinical status, relative to healthy donors. The presented flow diagrams visualize the Modulator/Source Doublets and their total Occurrence Counts for both the Doublet and Singlet summations (Fig. 6). Each flow diagram begins with a Quad dataset that is reduced to Doublets (edges) and Singlets (nodes). As explained previously with Table 1, the patterns of Occurrence Counts repeat for all Targets and LSVs across a specific Doublet. Therefore, it can be inferred that active Doublets affect all biomarker Targets to varying degrees, the larger the total Doublet Occurrence Count value the greater the effects on all its biomarkers Targets. The Targets and LSVs of each Doublet provides additional insight and is an area of continued analysis however, is beyond the scope of this paper.
To provide an initial comprehensive view of the CICD resultant data, two cohorts, Healthy or Cancer were averaged across each unique active Quad (133,000 total) generating an average Cancer and an average Healthy CICD Quad dataset. These two representative datasets are compared to highlight the information obtained via CICD modeling (Fig. 6). Both complete Quad datasets Fig. 6a, b and three Quad comparison datasets Fig. 6c, d, e are portrayed (note, dominant signals of MCP1 [CCL2], IP10 [CXCL10], IL8 and CD11c +). For each cancer patient and healthy volunteer, the same complete flow diagrams are generated and found in Additional file 2.
The prominent Doublet activity levels (edges) are emphasized in Fig. 7. The biggest differences found between cancer patients and healthy volunteers were modulated by MCP1, IP10, and IL8, seen amongst the top 10 modulatorsource doublets enhanced in cancer depicted in Fig. 7a. These chemokines were found to modulate several types of immune cells, including T cells, natural killer (NK) cells, monocytes, and dendritic cells (DC) (Fig. 7b).
The CICD model’s Quad dataset and resultant flow diagrams portray many of the biology’s accepted fundamental truths. The immune regulatory system is highly redundant and tends toward a relatively stable equilibrium among these multiple interdependent elements as seen through the multiple CICD defined relationships driving specific biomarkers activity. This homeostasis is maintained by numerous underlying physiological processes, many that can be possibly defined through CICD. The homogenous and heterogeneous characteristics of the system is modeled through the multiple identical relationships, yet different values found in the CICD datasets between both cancer and healthy cohorts as well as among individuals. This is first observed when examining the full Quad and Triplet data, the same relationship yet difference values are found across patients and cohorts that is then visualized by the ModulatorSource Doublets flow diagrams. The calculations of the subtle differences underneath the seemly similar raw data found with CICD provides the individualized modeling needed to unlock clinical pathways to tailor treatment to a patient’s specific immune status.
Discussion
The main premise of this work is that singular biomarker discriminators (single measurement at a single time) have not been successful in describing the interconnected nature of different regulators of immune function in peripheral blood of humans with cancer. One of the main challenges in clinical cancer immunology is that, despite multiple efforts, it has been difficult to define a single measurement/predictor of immune cell function that encompasses the state of immunity in a patient with cancer. This has been particularly challenging in efforts of developing predictive biomarkers for cancer immunotherapy, where the only biomarker of outcomes for antiPD1 therapy (tumor associated PDL1 expression) continues to be highly debated and remains of notable but limited clinical value. Part of the reason for this challenge is the acknowledged complexity of the numerous interacting components of human immune system regulation, not all of which play the exact same role in different patients. Thus, we undertook an effort to collectively examine the multitude of commonly used biomarkers of the state of systemic immunity, analyze their interactive properties based on their temporal variability and mathematically define (CICD) the overall state of systemic immunity, shared by cohorts of individuals, but, attained via similar but not identical pathways. Herein we present an integrated analysis of a multitude of preselected “biomarkers” of immunity that describe the difference between the state of systemic immunity in cancer patients, versus that of healthy volunteers.
The data reveals a dominant signal by two chemokines MCP1 (CCL2) and IP10 (CXCL10) that are known to affect multiple cell types, as suggested by our data, with a somewhat more prominent impact on CD11c + cells in peripheral blood. These data are consistent with our appreciation of the role of chemokines in advanced cancer [27] as well as the putative role of CD11c + immune cells in maintenance of chronic inflammation [28]. When looking at the biggest differences between cancer patients versus healthy controls, we noticed that MCP1, IP10, and IL8 were amongst the top modulators in cancer, acting on several types of immune cells, shown in Fig. 7. These results suggest that the main mechanisms of immunologic homeostasis of cancer are driven by a subset of chemokines for which there is published data demonstrating their role in cancer immunity. MCP1 (CCL2) and IL8 (CXCL8) are chemokines known to allow tumor progression by promoting tumor angiogenesis [29, 30]. However, MCP1 is also known to attract monocytes, NK cells, and T cells [31, 32]. Furthermore, IL8 may also modulate the immune system in a favorable manner, as it has been shown to act as a chemoattractant for T cells [33]. In contrast to MCP1 and IL8, IP10 (CXCL10) is able to inhibit angiogenesis [34, 35]. IP10 is secreted by cells in the context if inflammation and promotes a Th1 immune response involving T cells, NK cells, and DCs [36,37,38,39]. Thus, the modulatorsource interactions that CICD associated from peripheral blood samples aligns with observations in the literature.
Hence, in a purely mathematical way, the current analyses have identified a plausible network of cells/cytokines (Modulator/Source) that play a role in sustaining the state of immunity in advanced cancer. This provides a starting point to a more focused analyses regarding the role of these dominant immune mediators in patients with cancer, and offers an analytical tool that may be used to test additional biomarkers, relevant to MCP1 (CCL2) and IP10 (CXCL10) and in more detail describe the network of interactions in the absence or presence of therapeutic interventions, and resultant clinical outcomes, possibly identifying novel therapeutic targets and/or prognostic/predictive biomarkers.
At the current juncture, the CICD model is only asking what biomarkers (measured in peripheral blood only) are working together regardless of previous knowledge of their function. CICD observes the interaction behaviors of the immune system, without considering the many possible underlying mechanisms that combine to produce an active interaction. Identification of the subgroups (sub systems) of interacting biomarkers that discriminate cancer patients from healthy volunteers, elucidate critical differences in the states of immune homeostasis that may be relevant to the natural history of the malignancy as some of these elements of homeostasis (PD1 + immune cells) are targets of modern cancer immunotherapy.
A fundamental feature of CICD is the Knowledge Model. Unlike common approaches that focus on a limited number of biomarker interactions, CICD’s goal was to include all possible immune biomarker interactions to begin to unravel the immune regulatory complexity. The plasticity of the Knowledge Model allows future studies to expand and/or modify its interactions and biomarkers to assist in revealing biological interactions thru mathematics. The Knowledge Model is the unique feature of CICD that allows a flexible analysis of biomarkers’ interaction network without any experimental constraints based on previously investigated interactions. The innovative Knowledge Model lists all biologically possible causes or relationships between the biomarkers studied, whether or not they have previously been studied. Future studies look to modify the Knowledge Model by either narrowing or expanding the number and/or type of relationships to further explore the immune interaction network are currently under way. The flexibility of the Knowledge Model is fundamental in uncovering the complexity of the system.
The second essential requirement for CICD is the sequential collection of peripheral blood samples. This provides the necessary information to specifically assess an individual without any theoretical estimation of parameters. As seen reflected in the CICD results data, the raw data also contains similarities and differences among participants. Through CICD, the specific underlying causes that manifest in these observed biomarker oscillations can be quantified. Future studies will focus on determining the number blood draw data points required for optimal analysis reliability.
It is well established that biological systems are complex and redundant systems that will effectively maintain functionality if a failure occurs within the system. The human immune system illustrates this concept through numerous immunological mediators (biomarkers) working side by side to protect the body from disease. This presents a most interesting math modeling challenge to engineers that expect laws of behavior to be known exactly and solutions to be predictable, unique and verifiable. Rather than hypothesize ODEs to solve, a Knowledge Model is used to specify equations for all possible biomarker interaction response, and a set of characterization parameters defined as a vector of unknowns. The net result is an underdetermined system with far fewer equations than unknowns [40]. Such systems have an infinite number of equally valid solutions, solvable via SVD methods and the pseudoinverse. It is wellknown that for this class of SVD problems, singular values are unique while singular vectors are not. Recent work [41,42,43] is enhancing SVD capability to compute one solution that is least effort and which satisfies additional modeling objectives and application needs. Through the SVD algorithm, one of an infinite number of equally valid and least effort solutions is computed, inherently providing a model of the redundant, regulatory/system homeostasis mechanisms of the human immune system. The redundancy and robustness utilized by the immune system to maintain functionality despite changes in the system due to disease can now be mathematically modeled via the power of SVD.
A comprehensive literature review of current biological math modeling techniques was completed to assess CICD’s approach [2,3,4,5,6,7, 13,14,15,16,17,18, 44,45,46,47,48,49,50,51,52,53,54,55,56,57]. A sidebyside comparison of several aspects of current techniques vs. the alternative CICD methodology is presented in Table 2 to elaborate the advantages of CICD. Most current models tend to speculate on a limited subset of potential interactions and attempt to recover model parameters to match a set of observed data, a predicative forward modeling strategy. CICD exploits the uncommon reverse modeling strategy to determine underlying causalities that is expressed in the clinical data to understand the complexity of the immune system. Overall CICD does not assume a postulated mathematical model; instead it measures outcomes in the sample data which provide critical insights into the generally unknown internal working of the immune system response. Mathematical modeling converts assumptions into conclusions with certainty always relative to choice of assumptions [44]. To have coincidence in conclusions it has been necessary for current models to keep to simple and limited number of assumptions. These lowresolution models provide limited descriptions of the immune system and the need to shift to dynamic comprehensive modeling is apparent [45].
Current math models utilized multiple mathematical approaches such as various differential equations and agentbased modeling to model a limited number of specific biological assumptions by measurable or estimated parameters. CICD revises this current structure by allowing thousands of assumptions, both known and unknown relationships to be quantified using one generalized ODE equation with only clinical data and without parameterization. Current models’ results are highly reliant on multiple parameter values. A major advantage of CICD is that its algorithms and thereby its results are not dependent on the accuracy of estimated or theorized parameters. CICD allows all underlying parameters to be represented by the unknown value. By isolating the values for all parameters, the computational burden of the model is dramatically reduced, and the system of equations simplified, thereby enabling CICD analysis to focus on the information obtained via SVD. This vector of calculated parameter values is another avenue for investigation in future studies. An additional unique and superior characteristic of CICD modeling is the ability to change how many and which specific biomarkers will be analyzed without an assumed knowledge of a biomarker’s mechanism in the immune system. The flexibility of CICD is boundless and enables the ability to tailor the biomarkers and relationships to the focus of the study. CICD results are based on a patient’s blood draw measurements alone. This individualized immune math model can potentially unlock tailored treatment plans for a person’s specific clinical status and the ability to monitor their immune system over the course of treatment as well as beyond. CICD’s application of the SVD algorithm in this unique way is the foundation of its approach. The SVD algorithm has enabled considerable advances in the modern world, such as facial recognition [58]. Now CICD uses the power of SVD to continue to advance understanding of the biological world.
Conclusions
CICD equations expose underlying causes, previously imperceptible nonlinear coupling biomarker interactions that together affect the observed populations of the biomarkers in the blood, thereby providing an insight rich snapshot of an individual’s state of immune homeostasis. Its main objective is discovery, to help lead to improvements in the efficacy of existing immune therapies (patient selection; drug combinations) and insight into new therapeutic targets that may significantly reduce the time to discover new therapeutics capable of meaningful clinical impacts.
This new perspective of biological modeling from the current predictive approach to the CICD descriptive approach unleashes vast potentials not only in understanding of the complexity of the immune system but continues the pathway to individual treatment plans. Possible applications of CICD are many, not only in other cancers but in other diseases as well, utilizing this model of complex, redundant, illdefined, biological systems. This collaboration between biology and engineering disciplines exemplifies how “the purpose of computation is insight, not numbers” is pivotal in the quest for a cure [59].
Availability of data and materials
The data that support the findings of this study are available from Math for Medicine, Inc. but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of Math for Medicine, Inc.
Abbreviations
 CICD:

Cancer immune control dynamics
 SVD:

Singular value decomposition
 ICI:

Immune checkpoint inhibitors
 ODE:

Ordinary differential equation
 LSV:

Left singular vectors
 RSV:

Right singular vectors
References
McKinney EF, Smith KG. T cell exhaustion and immunemediated diseasethe potential for therapeutic exhaustion. Curr Opin Immunol. 2016;43:74–80.
Bianca C, Pennisi M. Immune system modelling by topdown and bottomup approaches. Int Math Forum. 2012;7:109–28.
Callard RE, Yates AJ. Immunology and mathematics: crossing the divide. Immunology. 2005;115:21–33.
Eftimie R, Gillard JJ, Cantrell DA. Mathematical models for immunology: current state of the art and future research directions. Bull Math Biol. 2016;78:2091–134.
Germain RN, MeierSchellersheim M, NitaLazar A, Fraser ID. Systems biology in immunology: a computational modeling perspective. Annu Rev Immunol. 2011;29:527–85.
Motta S, Pappalardo F. Mathematical modeling of biological systems. Brief Bioinform. 2013;14:411–22.
Shou W, Bergstrom CT, Chakraborty AK, Skinner FK. Theory, models and biology. Elife. 2015;4:e07158. https://doi.org/10.7554/eLife.
Juang JN. Applied system identification. Englewood Cliffs: Prentice Hall; 1994. p. 394.
Junkins JL, Kim Y. Introduction to dynamics and control of flexible structures: American Institute of Aeronautics and Astronautics. https://books.google.com/books?id=WbS6K_LU7DoC; 1993.
Sowa JF. Knowledge representation: logical, philosophical, and computational foundations (ed. Przemieniecki, J.S.). Przemieniecki JS, editor: Pacific Grove: Brooks/Cole 2000.
Golub GH, Van Loan CF. Matrix computations. London: Johns Hopkins University Press; 1983.
Zadeh L, Desoer C. Linear system theory: the state space approach. McGrawHill, New York. https://books.google.com/books?id=CEVdyQEACAAJ; 1963.
Csete ME, Doyle JC. Reverse engineering of biological complexity. Science. 2002;295:1664–9.
Ludewig B, Stein JV, Sharpe J, CervantesBarragan L, Thiel V, Bocharov G. A global “imaging’’’ view on systems approaches in immunology.” Eur J Immunol. 2012;42:3116–25.
Villaverde AF, Banga JR. Reverse engineering and identification in systems biology: strategies, perspectives and challenges. J R Soc Interface. 2014;11:20130505.
Eftimie R, Bramson JL, Earn DJ. Interactions between the immune system and cancer: a brief review of nonspatial mathematical models. Bull Math Biol. 2011;73:2–32.
Kim PS, Levy D, Lee PP. Modeling and simulation of the immune system as a selfregulating network. Methods Enzymol. 2009;467:79–109.
Vodovotz Y, Xia A, Read EL, BassaganyaRiera J, Hafler DA, Sontag E, et al. Solving immunology? Trends Immunol. 2017;38:116–27.
Strang G, Aarikka K. Introduction to applied mathematics. Wellesley: WellesleyCambridge Press; 1986. p. 16.
Tarantola A. Inverse problem theory and methods for model parameter estimation. London: Society for Industrial and Applied Mathematics; 2005. p. 89.
Josephson J, Josephson S. Abductive inference: computation, philosophy, technology. Cambridge: Cambridge University Press; 1996.
Barker S. Aircraft as a system of systems: a business process perspective. Warrendale: SAE International; 2019. https://doi.org/10.4271/T139.
Thakar J, Poss M, Albert R, Long GH, Zhang R. Dynamic models of immune responses: what is the ideal level of detail? Theor Biol Med Model. 2010;7:35.
Brusic V, Petrovsky N. Immunoinformatics and its relevance to understanding human immune disease. Expert Rev Clin Immunol. 2005;1:145–57.
GómezPérez A, FernándezLópez M, Corcho O. Ontological engineering. Advanced information and knowledge processing. Heidelberg: Springer; 2003.
Chan TF. An improved algorithm for computing the singular valure decomposition. ACM Trans Math Softw. 1982;8:72–83.
Nagarsheth N, Wicha MS, Zou W. Chemokines in the cancer microenvironment and their relevance in cancer immunotherapy. Nat Rev Immunol. 2017;17:559–72.
Winslow GM, Papillion AM, Kenderes KJ, Levack RC. CD11c+ Tbet+ memory B cells: immune maintenance during chronic infection and inflammation? Cell Immunol. 2017;321:8–17.
Chow MT, Luster AD. Chemokines in cancer. Cancer Immunol Res. 2014;2:1125–31.
Payne AS, Cornelius LA. The role of chemokines in melanoma tumor growth and metastasis. J Invest Dermatol. 2002;118:915–22.
Conti I, Rollins BJ. CCL2 (monocyte chemoattractant protein1) and cancer. Semin Cancer Biol. 2004;14:149–54.
Fuentes ME, Durham SK, Swerdel MR, Lewin AC, Barton DS, Megill JR, et al. Controlled recruitment of monocytes and macrophages to specific organs through transgenic expression of monocyte chemoattractant protein1. J Immunol. 1995;155:5769–76.
Taub DD, Anver M, Oppenheim JJ, Longo DL, Murphy WJ. T lymphocyte recruitment by interleukin8 (IL8). IL8induced degranulation of neutrophils releases potent chemoattractants for human T lymphocytes both in vitro and in vivo. J Clin Invest. 1996;97:1931–41.
Luster AD, Greenberg SM, Leder P. The IP10 chemokine binds to a specific cell surface heparan sulfate site shared with platelet factor 4 and inhibits endothelial cell proliferation. J Exp Med. 1995;182:219–31.
Campanella GS, Colvin RA, Luster AD. CXCL10 can inhibit endothelial cell proliferation independently of CXCR3. PLoS ONE. 2010;5:e12700.
Groom JR, Richmond J, Murooka TT, Sorensen EW, Sung JH, Bankert K, et al. CXCR3 chemokine receptorligand interactions in the lymph node optimize CD4+ T helper 1 cell differentiation. Immunity. 2012;37:1091–103.
Wendel M, Galani IE, SuriPayer E, Cerwenka A. Natural killer cell accumulation in tumors is dependent on IFNgamma and CXCR3 ligands. Cancer Res. 2008;68:8437–45.
Liu M, Guo S, Stiles JK. The emerging role of CXCL10 in cancer (review). Oncol Lett. 2011;2:583–9.
Bonecchi R, Bianchi G, Bordignon PP, D’Ambrosio D, Lang R, Borsatti A, et al. Differential expression of chemokine receptors and chemotactic responsiveness of type 1 T helper cells (Th1s) and Th2s. J Exp Med. 1998;187:129–34.
Boyd S. Introduction to linear dynamical systems. Stanford: Stanford Engineering Everywhere; 2007.
Stojnic M. Linear underdetermined systems with sparse solutions: redirecting a challenge? arXiv preprint. arXiv:1304.0004; 2013.
Qian J, Andrew AL, Chu D, Tan RCE. Methods for solving underdetermined systems. Numer Linear Algebra Appl. 2017;25:e2127. https://doi.org/10.1002/nla.2127.
Pena JM, Sauer T. SVD update methods for large matrices and applications. Linear Algebra Appl. 2019;561:41–62.
Gunawardena J. Models in biology: “accurate descriptions of our pathetic thinking.” BMC Biol. 2014;12:29.
Bocharov G, Volpert V, Ludewig B, Meyerhans A. Editorial: mathematical modeling of the immune system in homeostasis. Infection Disease Front Immunol. 2019;10:2944.
Azeloglu EU, Iyengar R. Good practices for building dynamical models in systems biology. Sci Signal. 2015;8:fs8.
Kabiri Chimeh M, Heywood P, Pennisi M, Pappalardo F, Richmond P. Parallelisation strategies for agent based simulation of immune systems. BMC Bioinform. 2019;20:579.
Castiglione F, Ghersi D, Celada F. Computer modeling of clonal dominance: memoryantinaive and its curbing by attrition. Front Immunol. 2019;10:1513.
Figueredo GP, Siebers PO, Aickelin U. Investigating mathematical models of immunointeractions with earlystage cancer under an agentbased modelling perspective. BMC Bioinform. 2013;14(Suppl 6):S6.
Grossman Z. Immunological paradigms, mechanisms, and models: conceptual understanding is a prerequisite to effective modeling. Front Immunol. 2019;10:2522.
MeierSchellersheim M, Varma R, Angermann BR. Mechanistic models of cellular signaling, cytokine crosstalk, and cell–cell communication in immunology. Front Immunol. 2019;10:2268.
Peskov K, Azarov I, Chu L, Voronova V, Kosinsky Y, Helmlinger G. Quantitative mechanistic modeling in support of pharmacological therapeutics development in immunooncology. Front Immunol. 2019;10:924.
Belfiore M, Pennisi M, Arico G, Ronsisvalle S, Pappalardo F. In silico modeling of the immune system: cellular and molecular scale approaches. Biomed Res Int Int. 2014;2014:371809.
Bianca C, Chiacchio F, Pappalardo F, Pennisi M. Mathematical modeling of the immune system recognition to mammary carcinoma antigen. BMC Bioinform. 2012;13(Suppl 17):S21.
Cappuccio A, Tieri P, Castiglione F. Multiscale modelling in immunology: a review. Brief Bioinform. 2016;17:408–18.
Depillis L, Gallegos A, Radunskaya A. A model of dendritic cell therapy for melanoma. Front Oncol. 2013;3:56.
Norton KA, Gong C, Jamalian S, Popel AS. Multiscale agentbased and hybrid modeling of the tumor immune microenvironment. Processes (Basel). 2019. https://doi.org/10.3390/pr7010037.
Zeng G. Facial recognition with singular value decomposition. In: Elleithy K, editor. Advances and innovations in systems, computing sciences and software engineering. Dordrecht: Springer; 2007.
Hamming R. Numerical methods for scientists and engineers. New York: Courier Corporation; 2012.
Acknowledgements
This manuscript is dedicated to the memory of our dear Harry Frisch, the heart and soul behind this effort, whose passion to reveal the dynamic immune system through mathematics has made this work possible. We thank Vera J. Suman, Ph.D. Mayo Clinic and Wolfgang Kliemann, Ph.D. Iowa State University for their valuable assistance. Figures 1, 2, 3, 5 and 7 were created with Biorender.com.
Funding
Mayo Foundation for Medical Education and Research. Math for Medicine. M4M did not provide any monetary funding to the project but did make its personnel and software available to the project. M4M personnel did not participate in collection, analysis, and interpretation of the processed data. M4M did contribute ideas, reports and diagrams to the project.
Author information
Authors and Affiliations
Contributions
HF, JT, SM conceptualization. HF, AS, VM, WN data curation. HF, AS, VM, AL, SM formal analysis. AS and SM funding acquisition. HF, AS, VM, WN, SM investigation. HF, AS, VM, J.T, SM methodology. AS and SM project administration. WN and SM resources. HF, AS, VM software. SM supervision. HF, AS, VM, AL, SM validation. AS and VM visualization. HF, AS, VM, SM writingoriginal draft preparation. HF, AS, VM, JT, LB, SM writingreview and editing. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All research has been performed in accordance with relevant guidelines and regulations. This study has been approved by the Institutional Review Board of Mayo Clinic. All participants provided written informed consent for study participation.
Consent for publication
Not applicable.
Competing interests
SNM has intellectual property that is involved in this work. Math for Medicine, Inc. (M4M) is the owner of the proprietary processing and analysis software used to produce the reports and data presented in the paper (licensed from the Mayo Clinic). There was no funding provided by M4M for the execution of these studies. Coinvestigators from the Mayo Clinic were supported by institutional funds for the purpose of the study. All other authors declare no potential conflicts of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1. Flow Cytometry and Cytokine Measurement
: Detailed methods used to obtain biomarker measurements.
Additional file 2. Individual Peripheral Blood Derived Serial Biomarker Data and Flow Diagrams
: A graph of peripheral blood biomarker measurements and flow diagram is provided for each cancer patient and each healthy volunteer. All Quads for the individual are included in each diagram. Graphs are in the same format as Fig. 4. Flow diagrams are in the same format as Fig. 6.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Frisch, H.P., Sprau, A., McElroy, V.F. et al. Cancer immune control dynamics: a clinical data driven model of systemic immunity in patients with metastatic melanoma. BMC Bioinformatics 22, 197 (2021). https://doi.org/10.1186/s12859021040257
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12859021040257
Keywords
 Math modeling
 Oncology
 Melanoma
 Systemic immunity
 Peripheral blood biomarkers