An Integrative multi-lineage model of variation in leukopoiesis and acute myelogenous leukemia
© The Author(s) 2017
Received: 8 February 2017
Accepted: 11 August 2017
Published: 25 August 2017
Acute myelogenous leukemia (AML) progresses uniquely in each patient. However, patients are typically treated with the same types of chemotherapy, despite biological differences that lead to differential responses to treatment.
Here we present a multi-lineage multi-compartment model of the hematopoietic system that captures patient-to-patient variation in both the concentration and rates of change of hematopoietic cell populations. By constraining the model against clinical hematopoietic cell recovery data derived from patients who have received induction chemotherapy, we identified trends for parameters that must be met by the model; for example, the mitosis rates and the probability of self-renewal of progenitor cells are inversely related. Within the data-consistent models, we found 22,796 parameter sets that meet chemotherapy response criteria. Simulations of these parameter sets display diverse dynamics in the cell populations. To identify large trends in these model outputs, we clustered the simulated cell population dynamics using k-means clustering and identified thirteen ‘representative patient’ dynamics. In each of these patient clusters, we simulated AML and found that clusters with the greatest mitotic capacity experience clinical cancer outcomes more likely to lead to shorter survival times. Conversely, other parameters, including lower death rates or mobilization rates, did not correlate with survival times.
Using the multi-lineage model of hematopoiesis, we have identified several key features that determine leukocyte homeostasis, including self-renewal probabilities and mitosis rates, but not mobilization rates. Other influential parameters that regulate AML model behavior are responses to cytokines/growth factors produced in peripheral blood that target the probability of self-renewal of neutrophil progenitors. Finally, our model predicts that the mitosis rate of cancer is the most predictive parameter for survival time, followed closely by parameters that affect the self-renewal of cancer stem cells; most current therapies target mitosis rate, but based on our results, we propose that additional therapeutic targeting of self-renewal of cancer stem cells will lead to even higher survival rates.
Acute myelogenous leukemia (AML) is a cancer of the white blood cells stemming from the myeloid lineage that produces cells including neutrophils and monocytes . On average, four in 100,000 individuals will develop AML, with a median age of 67 years at diagnosis. AML has the highest mortality rate of all of the different types of leukemia . The French-American-British (FAB) co-operative group identified eight different subcategories of AML, M0-M7, based on the specific cell population from which AML arises . More recent classification systems now exist that include other criteria besides white cell morphology [4, 5]. Regardless, the majority of patients with AML receive essentially identical induction chemotherapy with drugs that target DNA replication . Unfortunately, for the patients over 55 years of age diagnosed with AML, the 5-year survival rate is 10% and the relapse rate is 80% . Given the variability of AML progression among patients, understanding determinants of disease progression will lead to therapeutic advances that result in lower rates of relapse and higher rates of survival.
Computational modeling is an effective tool to personalize therapies for improved patient outcomes. Model-based personalized therapy has improved several medical interventions including: glucose control in type I diabetes [8, 9]; automatic control of anesthesia dosage [10–12]; and pacemaker control of heart rate variability . Presently, AML treatment is not benefiting from model-based personalized approaches, in part due to a lack of multi-lineage AML computational models [14, 15]. Key features of a computational model of hematopoiesis includes: normal formation of mature cells from hematopoietic stem cells (HSCs) [16–24]; abnormal development leading to leukemia [18, 25–27]; treatment using chemotherapy or transplantation [17, 25, 26]; movement of cells between the bone marrow/tissue and the peripheral blood where AML is frequently measured [20, 28–35]; and the feedback signals from cytokines that participate in the regulation of hematopoiesis [20, 25, 26, 28]. Previous models that contain various elements of these key features exist (summarized in Additional file 1: Table S1 in Supplement 1). Here, we focus on the development of a multi-lineage AML model that captures patient variability and use the model to identify patient sub-types and identify potential novel therapeutic targets.
In this work, we develop a semi-mechanistic multi-lineage multi-compartment computational model of hematopoiesis and AML. The model describes the interactions of various lineages of hematopoietic stem cells, including neutrophils, lymphocytes, and monocytes. Patient data is used to constrain the model, creating a large number of data-consistent solutions that differ in their responses to chemotherapy. We apply the constraints from patient data to characterize normal hematopoiesis versus leukemic differentiation and growth in unique parameter sets representing 22,796 simulated patient dynamics. To identify trends in dynamical variability, we clustered the simulated dynamics into thirteen ‘representative patient’ groups. We find that our mathematical model can be used to demonstrate the typical variation in the dynamics of AML progression and patient clinical outcomes. Using sensitivity analysis and correlation analysis, we identified important parameters whose variation most affects the development of leukemia.
The leukopoiesis process consists of a set of discrete stages of cellular differentiation. There are three different stages of cell development representing the maturation of the stem cells into mature leukocytes (Fig. 1, solid black circles). The first stage contains the hematopoietic stem cells (HSCs) that either asymmetrically self-renew or differentiate. We assume this HSC compartment consists of many stages of cells that have the capability to differentiate into various lineages, as previous work has shown that the dynamics of mature cells are appropriately represented in a two-compartment system [36–38]. The stem cells proportionally differentiate into distinct progenitors, including common lymphoid progenitors and granulocyte/myeloid progenitors that then clonally amplify their branching lineages . All three lineages are necessary to demonstrate lineage dominance and cell-cell interactions amongst the three most-commonly measured white blood cell types (see Additional file 1: Supplement 3, Table S5 for further justification of these three lineages). Other multi-lineage models do not consider these three lineages [21–24]. Granulocyte/myeloid progenitors produce mature neutrophils and monocytes, which have relatively short life-spans and are replaced daily. Lymphoid progenitors branch into immature B-lymphocytes and T-lymphocytes, yielding more complex and less predictable dynamics than those of myeloid cells. Thus, we focused on studying myeloid cell production, destruction, and deviations associated with myeloid leukemias instead of the less predictable lymphoid leukemias. The process of proliferation and differentiation in our model ensures a physiological mechanism that repopulates the bone marrow after chemotherapy and ensures that the formation of a single mutated cell in the HSC compartment can result in AML.
Mature cells are maintained in three different compartments: the bone marrow or tissue (yellow compartment in Fig. 1), the circulating peripheral blood (red compartment in Fig. 1), and the non-circulating peripheral blood, or marginal pool (blue compartment in Fig. 1). Cells migrate from both the bone marrow and marginal pool to the peripheral blood and are recruited to the bone marrow or marginal pools from the peripheral blood [40–42]. Since the concentration of normal white blood cells in the bone marrow and peripheral blood may vary over several orders of magnitude , it is essential to compute the dynamics of the proliferation of cells within the bone marrow and the mobilization of cells. Other multi-lineage models do not demonstrate movement of cells amongst compartments [21–24], and mobilization is typically limited to neutrophils or monocytes [20, 27–31, 33–35]. Additionally, neutrophils and monocytes also remain in non-circulating peripheral blood pools, or marginal pools, to allow for massive demargination upon a trigger event, such as chemotherapy [40, 43].
Here, X 0 is the state of interest; p 0 is the proliferation term, typically described as asymmetric self-renewal (circular arrows in Fig. 1; based on ; described further in Additional file 1: Supplement 3); s −1 is the differentiation (specialization) term from a stem cell or progenitor state, X −1, to the state of interest, X 0; s 0 is the differentiation term from the state of interest, X 0, to a more mature state (straight solid arrows in Fig. 1); m −1 is the mobilization or margination term from the corresponding bone marrow or peripheral blood state, X b m/p b , to the state of interest, X 0 (forward dashed arrows in Fig. 1); r 0 is the recruitment term of X 0 to another compartment (backward dashed arrows in Fig. 1); and d 0 is the loss, or death rate (x’s in Fig. 1). Further details of this model are described in Additional file 1: Supplement 3. A few key physiological details, including feedback and cell-cell interactions in the multi-lineage model, are described below.
Here, k is the Michaelis-Menten constant; C is the concentration of the cell state that produces cytokines for feedback; r is the rate the feedback modifies; and X is the cell state targeted by cytokines. To model the formation of AML from monocytes, we simulated the multiple-hit hypothesis and assumed all feedback signals did not effect a single normal monocyte progenitor, as expected from literature .
Inhibitory feedback prevents cells from growing uncontrollably in any state. Thus, our model has several inhibition mechanisms to regulate the concentration of cells. In particular, inhibition feedback modifies every rate except cell death in our model. All movement is inhibited by the concentration of cells in the compartment to which cells are moving to prevent overcrowding using Michaelis-Menten type kinetics (Eq. 2). Proliferation and differentiation are hindered by inhibiting the associated self-renewal probability of stem cells . For progenitor cells, this inhibition occurs from cytokines produced by mature cells of the same lineage in the peripheral blood . However, for stem cell proliferation and differentiation, this inhibition occurs from a scaled combination of the concentration of stem cells to ensure that sufficient stem cells are in the hematopoietic system and the concentration of cells in the bone marrow (yellow compartment in Fig. 1) do not exceed the capacity of the marrow. All negative feedback signals are pictorially described in Additional file 1: Figure S1. The derivation of both the form for asymmetric self-renewal and the capacity of the bone marrow are in Additional file 1: Supplement 3.
In contrast, several physiological processes in hematopoiesis are triggered by positive feedback. In response to an inflammatory event, mature white blood cells proliferate to accommodate this response. We include two states in our model that demonstrate the effect of debris clearance due to chemotherapy. Monocytes are recruited into the bone marrow to become activated macrophages to help clear excessive cellular debris , which we have modeled as the ‘Apoptosis’ state. Macrophages assist in recruiting other white blood cells, such as lymphocytes, neutrophils, and monocytes, during inflammation or other diseases , which we have also incorporated into the model. Positive feedback is used in our model to demonstrate several processes, including inducing proliferation of macrophages by apoptotic debris; promoting the recruitment of neutrophils, lymphocytes, and other monocytes into the tissue in response to high levels of macrophages; and increasing the clearing rate of apoptosis due to a high level of macrophages. All positive feedback signals are pictorially described in Additional file 1: Figure S2. Additional file 1: Table S4 in Supplement 2 describes each of the feedback processes in detail.
Capabilities of multi-lineage hematopoiesis model
Healthy dynamical acceptability criteria
1. Peripheral blood recovers > 80 cells/ μL 
2. Stem cells recover > 1 cell/ μL 
4. Marginal pool is within one order of magnitude of peripheral blood 
5. Cell counts reduce with chemotherapy to < 20% of original value 
6. Recovery overshoot < 12 times value five days after overshoot [hematopoietic stem cell transplant patient data from Dr. Robert Nelson, shown in Additional file 1: Supplement 5].
7. Amplitude of the second peak of oscillations deviates < 18% from the amplitude of the first peak [cutoff calculated, explained in Additional file 5: Supplement 5].
Physiological capabilities of model
Peripheral blood recovers > 80 cells/ μ L: This ensures that the peripheral blood sufficiently recovers after chemotherapy treatment, and the final peripheral blood cell count must be > 80 cells/ μL .
Stem cells recovery > 1 cell/ μ L: The final stem cell concentration must be between 1 and 100 cells/ μL. If the stem cell concentration is within the acceptable bounds, the stem cells should also show recovery after chemotherapy occurring .
Final dynamic values are within acceptable ranges: The final uni-lineage peripheral blood cell counts after chemotherapy must be within the upper limits of expected normal peripheral blood cell counts [40, 59]. Neutropenia (low number of neutrophils in the peripheral blood) is very common after undergoing chemotherapy , but neutrophilia mostly occurs if another underlying disease exists to cause the high number of neutrophils. We do not want to model neutrophilia during recovery from chemotherapy because that is not a typical response in patients with AML. The values each cell state must be within are depicted in Additional file 1: Supplement 4, Table S6.
Marginal pool is within one order of magnitude of peripheral blood: The final concentration of the cells in the marginal pool must be within one order of magnitude of the final concentration of the peripheral blood because the size of the marginal pool and peripheral blood compartments are approximately the same [43, 61, 62]. Though a marginal pool can exist for lymphocytes, it is much smaller than that of neutrophils and monocytes . Thus, we do not incorporate a lymphocyte marginal pool in our model.
Cell counts reduce with chemotherapy to <20 % of original value: The concentration of cells remaining at the end of chemotherapy should be less than 20% of the concentration of cells before chemotherapy begins .
Recovery overshoot < 12 times value five days after overshoot: Patients experience lymphopenia, or low white blood cell concentrations, during chemotherapy. After the chemotherapy regimen is completed, patients’ cell counts almost always increase. In some instances, patient cell counts increase and peak above their final homeostatic cell concentration. Here, we define this overshoot level as the maximum concentration of cells after chemotherapy administration is completed. If the overshoot level is more than twelve times the concentration of cells five days after the overshoot (or the concentration of cells at the end of the simulation, whichever comes first), the solution is rejected. This overshoot was the maximum observed in the HSCT patient data from the Simon Cancer Center (Additional file 1: Figure S6 in Supplement 5).
Amplitude of the second peak of oscillations deviates < 18% from the amplitude of the first peak: Several simulations produced oscillations. Though some small oscillations are reasonable and can occur physiologically, large oscillations are unlikely because they will not maintain normal homeostasis in the body and most patients who are treated for AML do not report large oscillations in cell counts. We manually classified 100 simulations for sufficiently dampened oscillations. We calculated that a simulation was sufficiently dampened if the cell concentration of the second peak value post-chemotherapy was less than an 18% decrease from the first peak post-chemotherapy. Thus, we implemented this cutoff for all of our simulations. An example of a simulation that does not meet this criteria and a simulation that does meet this criteria are shown in Additional file 1: Figure S7 in Supplement 5.
For each uni-lineage parameter search, two criteria were salient in determining which models met the physiological features required based on responses to chemotherapy. The remainder of the criteria are mostly robust in constraining the parameter space. One of the salient criteria was the lack of oscillations for each uni-lineage model. This implies that the need to maintain a steady-state concentration of cells for each of the lineages independently is one of the key constraints in normal hematopoiesis and is consistent with parameters that lead to balanced feedback. The second criteria that largely determined the acceptable uni-lineage models were the final peripheral blood acceptable dynamics of the neutrophil uni-lineage model, stem cell recovery for the lymphocyte uni-lineage model, and sufficient response to chemotherapy for the monocyte uni-lineage model (indicated with asterisks in Fig. 3). This reflects the necessity of maintaining high concentrations of neutrophils in homeostasis and the ability of cells to repopulate cell populations for appropriate recovery. If we relaxed the constraints imposed on the model set by the features of normal physiological dynamics, the acceptable parameter sets would be larger as well.
Parameter constraints with separatrix method
The separatrix method distinguishes between parameter regions of acceptable solutions and regions of unacceptable solutions, and in log-space, this manifests as a clear separatrix. We first constrain to the acceptable range of each parameter (dashed black box in Fig. 4). For each unique pair of uni-lineage parameters, we then divide the parameter range rectangle into 10×10 equally sized bins. If the acceptable parameter set is uniformly distributed in parameter space, every bin has a cutoff of (number of acceptable simulations)/100 acceptable parameter sets. Thus, each of the 100 bins is either acceptable or unacceptable depending on whether more acceptable parameter sets existed in the bin than the cutoff value. We then further constrain this acceptable parameter density by removing corners of the range rectangle that did not maintain a cutoff density of acceptable parameter sets (dashed green lines in Fig. 4). Thus, we found a set of inequalities for each unique pair of parameter sets in the uni-lineage parameterization in log space. We used this separatrix constraint on the parameter search space for multi-lineage parameterization.
The separatrices of each pair-wise parameter relationship confirm parameter relationships from previous works, and we identified additional required relationships in the acceptable solutions. As Getto et al.  and Stiehl et al.  show, the mitosis rate of stem cells and the probability of self-renewal are inversely related. However, we found that this also extends to all progenitor cells (Fig. 4 a-c). Specifically, a linear relationship emerges between the log of the self-renewal probability and the log of the mitosis rate. Solutions that have larger mitosis rate or self-renewal probabilities produce oscillations (red x’s in Fig. 4), which we constrained against in our dynamical acceptability criteria (Criteria 7 in Table 1). We also found a similar inverse linear relationship between the log of the self-renewal probability of neutrophil progenitors and the log of the homeostatic constant of neutrophils in the peripheral blood; solutions that have lower values of either of these two parameters becomes unacceptable (Fig. 4 d). Overall, this finding means that a specific relationship is maintained between the probability of self-renewal and the mitosis rate for all cells that are capable of self-renewing to ensure that the cell population remains steady. Additionally, due to the high concentration of neutrophils, a constrained feedback mechanism exists between the self-renewal probability of neutrophil progenitors and the cytokines produced by neutrophils in the peripheral blood to maintain appropriate homeostatic concentrations of those cells. Overall, we confirmed the dependence on the mitosis rate and self-renewal probability of stem cells found by other groups (Fig. 4 e), but we also found additional dependencies on these same parameters of progenitor cells, as well as a specific dependence of the cytokines produced by neutrophils in the peripheral blood and the self-renewal probability of neutrophil progenitor cells.
Multi-lineage acceptable solutions
Sensitivity analysis reveals importance of self-renewal probability and cytokines
We carried out a sensitivity analysis on the multi-lineage solutions to determine the most important parameters whose variation most greatly affected output. Using a partial rank correlation coefficient (PRCC) method [66–68], we determined the parameters that were most impactful in determining concentrations of stem cells, neutrophils in the peripheral blood, lymphocytes in the peripheral blood, and monocytes in the peripheral blood. We found that, in general, the probability of self-renewal and the Michaelis-Menten constant that modifies both the rate of self-renewal probability and the movement of mature cells into the peripheral blood were the most important parameters for their respective cell states. This indicates that the probability of a cell to self-renew is very sensitive to changes in the system and can cause very drastic outcomes in the final cell concentration if this rate is affected, which corroborates with findings in Marciniak-Czochra et al. . Additionally, for stem cells, the feedback cytokines from the concentration of cells in the bone marrow and the concentration of stem cells were very important parameters. These constants also modify the probability of self-renewal in the stem cells. Finally, many parameters that are associated with neutrophil homeostasis (self-renewal probability of neutrophil progenitors, mitosis rate, Michaelis-Menten constant, and mobilization rate) all were very important in the stem cell concentration. This is probably due to the large concentration of neutrophils affecting the homeostasis values of stem cells. The PRCC results are in Additional file 1: Supplement 6.
Representative patient clusters
AML is a derivative of common granulocyte macrophage progenitor cells [3, 69]. Thus, leukemia is modeled parallel to the manner in which monocytes are modeled. We assume that the cancer stem cells self-renew and proliferate at the same rate as monocytic progenitor cells. These then differentiate into mature cancer within the bone marrow, which can mobilize into the peripheral blood. The only difference in the model between leukemic stem cells and monocytes is that the leukemic cells do not maintain homeostasis, so all feedback that regulates self-renewal [18, 70], differentiation, and movement were removed and the cancer cells cannot be recruited back to the tissue. Since we have modeled feedback from various signals, the mutated cancer cell that grows indefinitely is the product of several mutations that cause all feedback to be inhibited, as currently hypothesized in literature .
For a simple validation to ensure that cancer is not growing too slowly in our model, we can compare overall survival times of the simulations without treatments to trials in which patients received low-dose treatments. The majority of clinical patients who are unable to receive traditional chemotherapy and received hydroxyurea as a treatment post-diagnosis with AML died within one year of diagnosis even with favorable stratification . We find the survival time in all of our simulations to be less than four months in after diagnosis of AML >20% blasts (Fig. 8 b). Thus, cancer growth fits within an expected bound of patient survival.
Cancer parameter bounds
a M c2
m r M c2
0.5m r M2
m b Mc
0.0001m b M
m b M
When a patient is diagnosed with AML, physicians assess the white blood cell counts of the patient to guide the course of action for treating the leukemia. Many patients are successfully treated with a stem cell or bone marrow transplant from a matched donor. However, for those patients who are unable to receive a transplant, chemotherapy is administered periodically to lower the cancer load on the patient. Generally, a high white blood cell count makes the patient a candidate for chemotherapy. However, it is likely that factors other than the absolute cell count are important for determining a patient’s prognosis. Using a multi-lineage model of the formation of white blood cells, we find that (1) certain physiological rate relationships are necessary to prevent unstable cell population growth and (2) the rate of growth of the cancer is an important prognostic factor in determining the survival time of patients.
We determined that there are specific physiological rates and cell concentrations that are crucial in maintaining homeostasis. As expected, progenitor cells are very important in homeostatic mechanisms. A sensitivity analysis demonstrated that the probability of self-renewal of all progenitor cells and parameters that modify this probability are the most important parameters in hematopoietic dynamics. Other groups have found similar results  in the context cancer [75, 76] and specifically AML . Additionally, we found an inverse relationship between the self-renewal probability of stem cells and the mitosis rate of stem cells (Fig. 4 e), corroborating findings of other groups [38, 52]. This relationship between mitosis rate and self-renewal probability also exists for progenitor cells (Fig. 4 a-c), which was not found previously. In general, solutions that have high mitosis rates or self-renewal probabilities lead to oscillations in the cell concentrations. This oscillation is highly undesirable and can occur in diseases such as cyclic neutropenia. Previous work has shown that “re-entry" into the stem cell compartment was found to be one of the factors that control oscillation in cyclic neutropenia . Our work suggests that the analogous self-renewal probability and mitosis rate of the same cell, whether it is the stem cell or another progenitor cell, are both crucial to control to prevent diseases such as cyclic neutropenia. Chemotherapy that targets only one of these rates may not be sufficient in controlling physiological oscillations. Thus, we recommend that the feedback mechanisms that govern the self-renewal probabilities are explored as potential pharmaceutical targets in AML treatment.
We found that peripheral blood concentration levels of neutrophils and lymphocytes are important factors in maintaining homeostasis in our model (Fig. 6), specifically in the multi-lineage context. This was evidenced in three ways. First, when we identified the capabilities of our multi-lineage model, we found that lymphocyte dynamics constrained the range of acceptable solutions for all three lineages (Fig. 5). This is likely due to the high selection rate of lymphocytes in the thymus (dL3, which is modeled as the tissue/bone marrow compartment)  and low death rate of lymphocytes in the peripheral blood (dL3pb) ; this can lead to fast dynamic changes if not constrained appropriately. Second, the inverse relationship between the self-renewal probability of neutrophil progenitors and the homeostatic term for neutrophil progenitors is a novel relationship in hematopoietic modeling that tightly regulates healthy neutrophil behavior (Fig. 4 d). This means that in order to maintain healthy levels of neutrophils at homeostasis, either the fraction of neutrophil progenitor cells that self-renew has to be low or the feedback mechanism that maintains neutrophil homeostasis has to have a low threshold for turning self-renewal off. Thus, self-renewal is very tightly regulated to ensure that cells do not grow indefinitely. To reduce the cancer load for potential treatment, physicians could target the feedback mechanisms associated with self renewal probability, in particular focusing on cytokine signaling. Third, the large concentration of neutrophils cause the parameters associated with neutrophils to be important parameters in determining the concentration of stem cells, as determined by global sensitivity analysis. We confirmed the importance of neutrophils in relation to lymphocytes in the multi-lineage model by testing the robustness of the multi-lineage model in Additional file 1: Supplement 6. We also found that the neutrophils in our multi-lineage model might be over-constrained by comparing to clinical overshoot data in comparison to the uni-lineage models. Specifically, the neutrophil overshoot of the multi-lineage model only reflected the neutrophil overshoot in 85% of clinical data (Additional file 1: Figure S6 in Supplement 5). Thus, our model encompasses most of the dynamics of the patient population. Furthermore, this reflects the importance of the sensing mechanism of each of these cells to maintain appropriate homeostatic levels. When we modeled cancer, removing these homeostatic terms from our model allowed the cancer to grow indefinitely, as expected. Thus, treatment that targets the feedback receptors on cancer cells can drastically help reduce the cancer load.
Individual patients are likely to develop unique phenotypes of AML. The dynamics and relative ratios of the absolute concentrations of neutrophils, lymphocytes, and monocytes vary widely across the patient population. Furthermore, phenotypic variability of leukemia may depend on the specific source cell within the population of common granulocyte progenitors that first becomes cancerous. Here, we develop a mathematical model that is able to describe patient variability in AML. We find that though cancer will form when homeostatic mechanisms are altered, additional mutations that increase the mitosis rate of cancer will reduce the survival of the patient without intervention (Fig. 8 e). Current chemotherapy regimens inhibit DNA replication, which corresponds to inhibiting the mitosis rate of the cancer . However, we also find that if the mechanism that determines the probability of self-renewal of cancer cells is mutated, then this can be an additional target for pharmaceutical treatments (Fig. 8 f), similar to what has been found in other work [75, 77].
The multi-lineage model developed in this work can be modified to explore mechanisms governing hematopoiesis and leukopoiesis. For example, the model could be adapted to discriminate amongst chemotherapy regimens on simulated patients to identify characteristics of patients that would benefit from certain regimens. Then, this could lead to identification of treatment schedules optimized for individual patients. The patient clusters could be used to stratify real patients by using the dynamics of patients’ responses to chemotherapy to match to a specific simulation and its outcome. Additionally, since we found that the probability of self-renewal is a potential secondary target for chemotherapy, various specific feedback mechanisms could be incorporated to identify the exact cytokine signal that would be the best target for therapy. This could be further used to identify the feedback signals that are most sensitive to cause AML growth. Furthermore, the linear relationship in log-log space between mitosis rates and self-renewal probabilities that was found in our model could be tested experimentally to ensure these relationships exist, though mediating cytokines properly is very difficult. Though we did not find specific relationships amongst the different lineages, the signaling mechanisms that control bone marrow size could also be explored experimentally and compared with the model to regulate overproduction of a particular cell lineage with respect to other lineages. Specifically, GM-CSF is often used to regulate the stem cell production of neutrophils and macrophages . This could be manipulated to identify effects in lymphocyte counts and the production of AML. Overall, the multi-lineage model we present here can be extended to characterize many aspects of hematopoiesis.
The multi-lineage model of hematopoiesis and leukopoiesis developed in this work can be readily adapted and expanded to incorporate many other immunological effects. Various groups have already modeled bone marrow transplantation and its potential complications [31, 33, 34]. Using our model, bone marrow or stem cell transplantation and transplant rejection can be integrated to predict graft rejection, and different lymphocyte sub-types can be added to the model, such as natural killer cells, to model an alternative outcome of transplantation: graft versus host disease. Lymphocytic leukemias can be explored by altering the homeostatic mechanisms of lymphocytes; analogously, myelodysplastic syndromes and minimal residual disease could be further characterized by determining parameter changes that lead to appropriate model behavior. Additionally, a wide array of other immunological diseases or the wide array of patient responses (including a large T-cell repertoire) can be adapted into this multi-lineage model to characterize the complexity of normal function and other immunological diseases. More specifically, future work could characterize the mechanisms that cause spontaneous remission of AML in the presence of bacterial infection .
The multi-lineage mathematical model we have created of hematopoiesis and leukopoiesis can help identify how individuals differ in their white blood cell and leukemia production. This model is useful for several reasons. We have determined crucial parameters and parameter relationships that can be used as potential drug targets for both AML and other potential immunological disorders. For many chemotherapeutic drugs, DNA replication is targeted , which aligns with targeting the mitosis rate. However, a combination therapy that also addresses the probability of a cancer cell to self-renew could be potentially helpful for patients whose cancer becomes resistant to the initial therapy. There is no current way to experimentally determine how these different lineages interact and limit each other. Thus, this multi-lineage model is a very powerful tool that can aid in understanding how blood forms normally and misforms into AML in individual patients. In addition, the model captures multiple dynamics that represent specific patient subgroups. By clustering the wide population of individual differences, we find that one advantage of this multi-lineage model is that it can be readily extended to investigate personalization of treatment schedules of individual patients to prolong overall survival.
Support for this research comes from the National Institutes of Health under Grant R01HD073156 and the Indiana CTSI TL1 Award, funded in part by Grant UL1TR001108.
Availability of data and materials
The clinical data that support the findings of this study are available from Dr. Robert P. Nelson, Jr., but restrictions apply to the availability of these data, which were used under IRB approval for the current study, and so are not publicly available. The remaining data are however available from the authors upon reasonable request.
JMS and SR created the multi-lineage model. JMS completed the model analysis and drafted the written work. RPN contributed to the conception, interpretation, and revision of the work. AER, DMU, and TKU contributed to the study design, interpretation, and revision of the written work. All authors read and approved the final manuscript.
Ethics approval and consent to participate
De-identified clinical data was used under Purdue University’s IRB Protocol# 1011009928 from the Indiana University Simon Cancer Center from patients who received hematopoietic stem cell transplants during a study approved by the IRB of Indiana University-Purdue University Indianapolis.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Koeffler HP, Golde DW. Acute myelogenous leukemia: a human cell line responsive to colony-stimulating activity. Science. 1978; 200:1153–4.View ArticlePubMedGoogle Scholar
- Epidemiology S, Program ERS. Research Data (1973-2012), National cancer institute, DCCPS, Surveillance research program, Surveillance systems branch. 2015. www.seer.cancer.gov. Accessed 10 Jan 2017.
- Bennett JM, Catovsky D, Daniel MT, Flandrin G, Galton DA, Gralnick HR, Sultan C. Proposals for the classification of the acute leukaemias. french-american-british (fab) co-operative group. Br J Hematol. 1976; 33:451–8.View ArticleGoogle Scholar
- Grimwade D, Walker H, Harrison G, Oliver F, Chatters S, Harrison CJ, Wheatley K, Burnett AK, Goldstone AH. The predictive value of hierarchical cytogenetic classification in older adults with acute myeloid leukemia (aml): analysis of 1065 patients entered into the united kingdom medical research council aml11 trial. Blood. 2001; 98:1312–9.View ArticlePubMedGoogle Scholar
- Vardiman JW, Thiele J, Arber DA, Brunning RD, Borowitz MJ, Porwit A, Harris NL, Beau MML, Hellström-Lindberg E, Tefferi A, Bloomfield CD. The 2008 revision of the world health organization (who) classification of myeloid neoplasms and acute leukemia: rationale and important changes. Blood. 2009; 114:937–51.View ArticlePubMedGoogle Scholar
- Kimby E, Nygren P, Glimelius B. A systematic overview of chemotherapy effects in acute myeloid leukaemia. Acta Oncol. 2001; 40:231–52.View ArticlePubMedGoogle Scholar
- Baron F, Storb R. Hematopoietic cell transplantation after reduced-intensity conditioning for older adults with acute myeloid leukemia in complete remission. Curr Opinions Hematol. 2007; 14:145–51.View ArticleGoogle Scholar
- Parker RS, Doyle FJ, Peppas NA. A model-based algorithm for blood glucose control in type i diabetic patients. IEEE Trans Biomed Eng. 1999; 46:148–57.View ArticlePubMedGoogle Scholar
- Wong XW, Singh-Levett I, Hollingsworth LJ, Shaw GM, Hann CE, Lotz T, Lin J, Wong OS, Chase JG. A novel, model-based insulin and nutrition delivery controller for glycemic regulation in critically ill patients. Diabetes Technol Ther. 2006; 8:174–90.View ArticlePubMedGoogle Scholar
- Dua P, Dua V, Pistikopoulos EN. Model based parametric control in anesthesia. Eur Symp Comput Aided Process Eng. 2005; 20:1015–20.Google Scholar
- Mendez JA, Torres S, Reboso JA, Jagannivas HR, Hepsiba D. Model-based controller for anesthesia automation. IEEE Int Conf Autom Sci Eng. 2009;379–84.Google Scholar
- Jagannivas N, Hepsiba D. Control of anaesthesia concentration using model based controller. Int J Innov Res Technol. 2014; 1:237–44.Google Scholar
- Bogdan P, Jain S, Marculescu R. Pacemaker control of heart rate variability: A cyber physical system perspective. ACM Trans Embed Comput Syst (TECS); 12(1s):1–22.Google Scholar
- Chakrabarty A, Pearce SM, R P Nelson J, Rundell AE. Treating acute myeloid leukemia via hsc transplantation: A preliminary study of multi-objective personalization strategies. Am Control Conf (ACC). 2013:3790–5.Google Scholar
- Noble SL, Sherer E, Hannemann RE, Ramkrishna D, Vik T, Rundell AE. Using adaptive model predictive control to customize maintenance therapy chemotherapeutic dosing for childhood acute lymphoblastic leukemia. J Theor Biol. 2010; 264:990–1002.View ArticlePubMedGoogle Scholar
- Peng CA, Koller MR, Palsson BO. Unilineage model of hematopoiesis predicts self-renewal of stem and progenitor cells based on ex vivo growth data. Biotechnol Bioeng. 1996; 52:24–33.View ArticlePubMedGoogle Scholar
- Scholz M, Engel C, Loeffler M. Modelling human granulopoiesis under poly-chemotherapy with g-csf support. J Math Biol. 2005; 50:397–439.View ArticlePubMedGoogle Scholar
- Michor F, Hughes TP, Iwasa Y, Branford S, Shah NP, Sawyers CL, Nowak MA. Dynamics of chronic myeloid leukaemia. Nature. 2005; 435:1267–70.View ArticlePubMedGoogle Scholar
- Marciniak-Czochra A, Stiehl T, Ho AD, Jager W, Wagner W. Modeling of asymmetric cell division in hematopoietic stem cells–regulation of self-renewal is essential for efficient repopulation. Stem Cells Dev. 2009; 18:377–85.View ArticlePubMedGoogle Scholar
- Ho T, Clermont G, Parker RS. A model of neutrophil dynamics in response to inflammatory and cancer chemotherapy challenges. Comput Chem Eng. 2013; 51:187–96.View ArticleGoogle Scholar
- Colijn C, Mackey MC. A mathematical model of hematopoiesis–i. periodic chronic myelogenous leukemia. J Theor Biol. 2005; 237:117–32.View ArticlePubMedGoogle Scholar
- Colijn C, Mackey MC. A mathematical model of hematopoiesis–ii. cyclical neutropenia. J Theor Biol. 2005; 237:133–46.View ArticlePubMedGoogle Scholar
- Manesso E, Teles J, Bryder D, Peterson C. Dynamical modelling of haematopoiesis: an integrated view over the system in homeostasis and under perturbation. J R Soc Interface. 2012;10(80).Google Scholar
- Székely T, Burrage K, Mangel M, Bonsall MB. Stochastic dynamics of interacting haematopoietic stem cell niche lineages. PLOS Comput Biol. 2014; 10(9):e1003794.View ArticlePubMedPubMed CentralGoogle Scholar
- Moore H, Li NK. A mathematical model for chronic myelogenous leukemia (cml) and t cell interaction. J Theor Biol. 2004; 227:513–23.View ArticlePubMedGoogle Scholar
- DeConde R, Kim PS, Levy D, Lee PP. Post-transplantation dynamics of the immune response to chronic myelogenous leukemia. J Theor Biol. 2005; 236:39–59.View ArticlePubMedGoogle Scholar
- Stiehl T, Baran N, Ho AD, Marciniak-Czochra A. Clonal selection and therapy resistance in acute leukemias: mathematical modelling explains different proliferation patterns at diagnosis and relapse. J R Soc Interface. 2014; 11.Google Scholar
- Takumi K, Garssen J, de Jonge R, de Jong W, Havelaar A. Release kinetics and cell trafficking in relation to bacterial growth explain the time course of blood neutrophils and monocytes during primary salmonella infection. Int Immunol. 2005; 17:85–93.View ArticlePubMedGoogle Scholar
- Lahoz-Beneytez J, Elemans M, Zhang Y, Ahmed R, Salam A, Block M, Niederalt C, Asquith B, Macallan D. Human neutrophil kinetics: modeling of stable isotope labeling data supports short blood neutrophil half-lives. Blood. 2016; 127:3431–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Craig M, Humphries AR, Mackey MC. A mathematical model of granulopoiesis incorporating the negative feedback dynamics and kinetics of g-csf/neutrophil binding and internalization. Bull Math Biol. 2016; 78(12):2304–57.View ArticlePubMedGoogle Scholar
- Stiehl T, Ho A, Marciniak-Czochra A. The impact of cd34+ cell dose on engraftment after scts: personalized estimates based on mathematical modeling. Bone Marrow Transplant. 2014; 49(1):30–7.View ArticlePubMedGoogle Scholar
- Stiehl T, Ho AD, Marciniak-Czochra A. Assessing hematopoietic (stem-) cell behavior during regenerative pressure. In: A Systems Biology Approach to Blood. New York: Springer: 2014. p. 347–67.Google Scholar
- Ostby I, Rusten LS, Kvalheim G, Grottum P. A mathematical model for reconstitution of granulopoiesis after high dose chemotherapy with autologous stem cell transplantation. J Math Biol. 2003; 47:101–36.View ArticlePubMedGoogle Scholar
- Østby I, Kvalheim G, Rusten LS, Grottum P. Mathematical modeling of granulocyte reconstitution after high-dose chemotherapy with stem cell support: effect of post-transplant g-csf treatment. J Theor Biol. 2004; 231:69–83.View ArticlePubMedGoogle Scholar
- Engel C, Scholz M, Loeffler M. A computational model of human granulopoiesis to simulate the hematotoxic effects of multicycle polychemotherapy. Blood. 2004; 104:2323–31.View ArticlePubMedGoogle Scholar
- Stiehl T, Marciniak-Czochra A. Characterization of stem cells using mathematical models of multistage cell lineages. ath Comput Model. 2011; 53(7):1505–17.View ArticleGoogle Scholar
- Nakata Y, Getto P, Marciniak-Czochra A, Alarcón T. Stability analysis of multi-compartment models for cell production sytems. J Biol Dyn. 2012; 6:2–18.View ArticlePubMedGoogle Scholar
- Getto P, Marciniak-Czochra A, Nakata Y, dM. Vivanco M. Global dynamics of two-compartment models for cell production systems with regulatory mechanisms. Math Biosci. 2013; 245:258–68.View ArticlePubMedGoogle Scholar
- Metcalf D. Clonal analysis of proliferation and differentiation of paired daughter cells: action of granulocyte-macrophage colony-stimulating factor on granulocyte-macrophage precursors. Proc Natl Acad Sci. 1980; 77(9):5327–330.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith CW. Production, distribution, and fate of neutrophils In: Kaushansky K, Lichtman MA, Prchal JT, Levi MM, Press OW, Burns LJ, Caligiuri M, editors. Williams Hematology. New York: McGraw-Hill: 2015. Chap. 61.Google Scholar
- Douglas SD, Douglas AG. Production, distribution, and activation of monocytes and macrophages In: Kaushansky K, Lichtman MA, Prchal JT, Levi MM, Press OW, Burns LJ, Caligiuri M, editors. Williams Hematology. New York: McGraw-Hill: 2015. Chap. 68.Google Scholar
- Seet CS, Crooks GM. Lymphopoiesis In: Kaushansky K, Lichtman MA, Prchal JT, Levi MM, Press OW, Burns LJ, Caligiuri M, editors. Williams Hematology. New York: McGraw-Hill: 2015. Chap. 74.Google Scholar
- Klonz A, Wonigeit K, Pabst R, Westermann J. The marginal blood pool of the rat contains not only granulocytes, but also lymphocytes, nk-cells and monocytes: a second intravascular compartment, its cellular composition, adhesion molecule expression and interaction with the peripheral blood pool. Scandanavian J Immunol. 1996; 44:461–9.View ArticleGoogle Scholar
- Zhang CC, Lodish HF. Cytokines regulating hematopoietic stem cell function. Curr Opin Hematol. 2008; 15:307–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Gonda TJ, D’Andrea RJ. Activating mutations in cytokine receptors: Implications for receptor function and role in disease. Blood. 1997; 89(2):355–69.PubMedGoogle Scholar
- Mizuki M, Schw able J, Steura C, Choudhary C, Agrawal S, Sargin B, Steffen B, Matsumura I, Kanakura Y, B ohmer FD, M uller-Tidow C, Berdel WE, Serve H. Suppression of myeloid transcription factors and induction of stat response genes by aml-specific flt3 mutations. Blood. 2003; 101:3164–73.View ArticlePubMedGoogle Scholar
- Kohl TM, Ellwart SSJW, Hiddemann W, Spiekermann K. Kit exon 8 mutations associated with core-binding factor (cbf)–acute myeloid leukemia (aml) cause hyperactivation of the receptor in response to stem cell factor. Blood. 2005; 105:3319–21.View ArticlePubMedGoogle Scholar
- Levine RL, Huntly MLBJP, Loh ML, Beran M, Stoffregen E, Berger R, Clark JJ, Willis SG, Nguyen KT, Flores NJ, Estey E, Gattermann N, Armstrong S, Look AT, Griffin JD, Bernard OA, Heinrich MC, Gilliland DG, Druker B, Deininger MWN. The jak2v617f activating mutation occurs in chronic myelomonocytic leukemia and acute myeloid leukemia, but not in acute lymphoblastic leukemia or chronic lymphocytic leukemia. Blood. 2005; 106:3377–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Knudson AG. Mutation and cancer: statistical study of retinoblastoma. Proc Natl Acad Sci. 1971; 68(4):820–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Rubinow SI, Lebowitz JL. A mathematical model of the acute myeloblastic leukemic state in man. Biophys J. 1976; 16:897–910.View ArticlePubMedPubMed CentralGoogle Scholar
- van den Akker E, Satchwell TJ, Pellegrin S, Daniels G, Toye AM. The majority of the in vitro erythroid expansion potential resides in cd34– cells, outweighing the contribution of cd34+ cells and significantly increasing the erythroblast yield from peripheral blood samples. Haematologica. 2010; 95.Google Scholar
- Stiehl T, Baran N, Ho AD, Marciniak-Czochra A. Cell division patterns in acute myeloid leukemia stem-like cells determineclinical course:a model to predict patient survival. Cancer Res. 2015; 75:940–9.View ArticlePubMedGoogle Scholar
- Shochat E, Rom-Kedar V, Segel LA. G-csf control of neutrophils dynamics in the blood. Bull Math Biol. 2007; 69:2299–338.View ArticlePubMedGoogle Scholar
- Holländer GA, Widmer B, Burakoff SJ. loss of normal thymic repertoire selection and persistence of autoreactive t cells in graft vs host disease. J Immunol. 1994; 152:1609–17.PubMedGoogle Scholar
- Kim PS, Lee PP, Levy D. Mini-transplants for chronic myelogenous leukemia: A modeling perspective. Biol Control Theory: Curr Challenges. 2007; 357:3–20.Google Scholar
- van Furth R, Dulk MMCD-D, Mattie H. Quantitative study on the production and kinetics of mononuclear phagocytes during an acute inflammatory reaction. J Exp Med. 1973; 138(6):1314–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Whitelaw DM, Batho HF. The disribution of monocytes in the rat. Cell Tissue Kinet. 1972; 5:215–25.PubMedGoogle Scholar
- Schmaier AH, Lazarus HM, (eds).Concise Guide to Hematology. West Sussex, UK: Wiley-Blackwell; 2012.Google Scholar
- Le T, Bhushan V, (eds).First Aid for the USMLE Step 1 2012. New York: McGraw-Hill Education; 2012.Google Scholar
- Morstyn G, Souza L, Keech J, Sheridan W, Campbell L, Alton N, Green M, Metcalf D, Fox R. Effect of granulocyte colony stimulating factor on neutropenia induced by cytotoxic chemotherapy. Lancet. 1988; 331(8587):667–72.View ArticleGoogle Scholar
- Harlan JM. Leukocyte-endothelial interactions. Blood. 1985; 65:513–25.PubMedGoogle Scholar
- Jedrzejczak WW. Mobilization of the marginal pool of neutrophils with epinephrine. results in healthy persons, patients with neutropenias, patients with neutrophilias, and patients with changes in neutrophil count induced by cancer chemotherapy. Haematologica. 1979; 64:586–96.PubMedGoogle Scholar
- Wäsch R, Reisser S, Hahn J, Bertz H, Engelhardt M, Kunzmann R, Veelken H, Holler E, Finke J. Rapid achievement of complete donor chimerism and low regimen-related toxicity after reduced conditioning with fludarabine, carmustine, melphalan and allogeneic transplantation. Bone Marrow Transplant. 2000; 26:243–50.View ArticlePubMedGoogle Scholar
- Jain AK, Murty MN, Flynn PJ. Data clustering: a review. ACM Comput Surv (CSUR). 1999; 13:264–323.View ArticleGoogle Scholar
- Amari S, Wu S. Improving support vector machine classifiers by modifying kernal functions. Neural Netw. 1999; 12:783–9.View ArticlePubMedGoogle Scholar
- Zheng Y, Rundell A. Comparative study of parameter sensitivity analyses of the tcr-activated erk-ampk signaling pathway. IEEE Proc Syst Biol. 2006; 153:201–11.View ArticleGoogle Scholar
- Kinzer-Ursem TL, Linderman JJ. Both ligand- and cell-specific parameters control ligand agonism in a kinetic model of g protein–coupled receptor signaling. PLOS Comput Biol. 2007; 3(1):e6.View ArticlePubMedPubMed CentralGoogle Scholar
- Marino S, Hogue IB, Ray CJ, Kirschner DE. A methodology for performing global uncertainty and sensitivity analysis in systems biology. Theor Biol. 2008; 254:178–96.View ArticleGoogle Scholar
- Krivtsov AV, Twomey D, Feng Z, Stubbs MC, Wang Y, Faber J, Levine JE, Wang J, Hahn WC, Gilliland DG, et al. Transformation from committed progenitor to leukaemia stem cell initiated by mll–af9. Nature. 2006; 442(7104):818–22.View ArticlePubMedGoogle Scholar
- Haeno H, Levine RL, Gilliland DG, Michor F. A progenitor cell origin of myeloid malignancies. Proc Natl Acad Sci U S A. 2009; 106:16616–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Zeijlemaker W, Gratama JW, Schuurhuis G. Tumor heterogeneity makes aml a “moving target” for detection of residual disease. Cytom Part B: Clin Cytom. 2014; 86(1):3–14.View ArticleGoogle Scholar
- Döhner H, Estey EH, Amadori S, Appelbaum FR, Büchner T, Burnett AK, Dombret H, Fenaux P, Grimwade D, Larson RA, Lo-Coco F, Naoe T, Niederwieser D, Ossenkoppele GJ, Sanz MA, Sierra J, Tallman MS, Löwenberg B, Bloomfield CD. Diagnosis and management of acute myeloid leukemia in adults: recommendations from an international expert panel, on behalf of the european leukemianet. Blood. 2010; 115:453–74.View ArticlePubMedGoogle Scholar
- Schiffer CA, Gurbuxani S, Larson RA, Rosmarin AG. Clinical manifestations, pathologic features, and diagnosis of acute myeloid leukemia. Waltham. http://www.uptodate.com/contents/clinical-manifestations-pathologic-features-anddiagnosis-of-acute-myeloid-leukemia. Accessed 10 Jan 2017.
- Burnett AK, Milligan D, Prentice AG, Goldstone AH, McMullin MF, Hills RK, Wheatley K. A comparison of low-dose cytarabine and hydroxyurea with or without all-trans retinoic acid for acute myeloid leukemia and high-risk myelodysplastic syndrome in patients not considered fit for intensive treatment. Cancer. 2007; 109(6):1114–24.View ArticlePubMedGoogle Scholar
- Ashkenazi R, Gentry SN, Jackson TL. Pathways to tumorigenesis—modeling mutation acquisition in stem cells and their progeny. Neoplasia. 2008; 10(11):1170–111826.View ArticlePubMedPubMed CentralGoogle Scholar
- Gentry S, Ashkenazi R, Jackson T. A maturity-structured mathematical model of mutation, acquisition in the absence of homeostatic regulation. Math Model Nat Phenom. 2009; 4(3):156–82.View ArticleGoogle Scholar
- Stiehl T, Marciniak-Czochra A. Mathematical modeling of leukemogenesis and cancer stem cell dynamics. Math Model Nat Phenom. 2012; 7:166–202.View ArticleGoogle Scholar
- Cytarabine: Drug Information. Waltham. http://www.uptodate.com/contents/cytarabine-patient-drug-information. Accessed 10 Jan 2017.
- Metcalf D. Hematopoietic cytokines. Blood. 2008; 111(2):485–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Maywald O, Buchheidt D, Bergmann J, Schoch C, Ludwig WD, Reiter A, Hastka J, Lengfelder E, Hehlmann R. Spontaneous remission in adult acute myeloid leukemia in association with systemic bacterial infection—case report and review of the literature. Ann Hematol. 2004; 83(3):189–94.View ArticlePubMedGoogle Scholar