Reduced variance of bone-marrow transit time of granulopoiesis—a possible pathomechanism of human cyclic neutropenia

Abstract. Human cyclic neutropenia (CN) is a haematological disorder characterized by oscillations in the numbers of neutrophilic granulocytes and other blood cells with a stable period of approximately 21 days. In most cases the neutrophils oscillate well below normal values such that these patients are chronically neutropenic. A comprehensive concept of the origin of CN is proposed. It assumes an abnormally small variance of the transit time of bone marrow cells (compared to normal human granulopoiesis) for the origin of the characteristic cycles. Furthermore, a reduced responsiveness of the immature granulopoietic bone marrow cells to the mitotic feedback stimuli is assumed to account for the subnormal neutrophil peaks. Together with feedback control provided in a simulation model of normal human granulopoiesis these two abnormalities can explain experimental and clinical cell kinetic data for bone marrow and blood in CN.     

Blast cell and granulocyte production in human leukemia: pathophysiological concepts based on computer simulation using discrete modeling techniques

The granulocyte production of two patients suffering from leukemia was studied extensively by means of the tritiated thymidine method of cellular kinetics. The data obtained (1-h labeling index, pattern of cell labeling, labeling intensity, as well as other conventional parameters of bone marrow and blood) were used to develop a computer model (GPSS-language) to fit the observations. From these models, it was concluded that patients with leukemia may have an abnormal granulopoiesis, characterized by a high degree of inefficiency (premature cell death, skipping of divisions with undisturbed maturation). However, the underlying mechanisms may be quite different. While it cannot be excluded that in acute myelocytic leukemia there is a stem and/or progenitor cell pool that is highly ineffective but still capable of feeding some cells into the granulocytic pathway, it is nevertheless possible, as shown in plasma cell leukemia, that the ineffective granulopoiesis may be the result of direct or indirect interaction between the "leukemic" and the "normal" cell clone.     

Response kinetics of radiation-induced micronucleated reticulocytes in human bone marrow culture

The frequency of micronucleated reticulocytes (MN-RETs) in the bone marrow or peripheral blood is a sensitive indicator of cytogenetic damage. While the kinetics of MN-RET induction in rodent models following irradiation has been investigated and reported, information about MN-RET induction of human bone marrow after radiation exposure is sparse. In this report, we describe a human long-term bone marrow culture (LTBMC), established in three-dimensional (3D) bioreactors, which sustains long-term erythropoiesis. Using this system, we measured the kinetics of human bone marrow red blood cell (RBC) and reticulocyte (RET) production, as well as the kinetics of human MN-RET induction following radiation exposure up to 6Gy. Human bone marrow established in the 3D bioreactor demonstrated an average percentage of RBCs among total viable cells peaking at 21% on day 21. The average percentage of RETs among total viable cells reached a maximum of 11% on day 14, and remained above 5% by day 28, suggesting that terminal erythroid differentiation was still active. Time- and dose-dependent induction of MN-RET by gamma radiation was observed in the human 3D LTBMC, with peak values occurring at approximately 3 days following 1Gy irradiation. A trend towards delayed peak to 3-5 days post-radiation was observed with radiation doses ≥2Gy. Our data reveal valuable information on the kinetics of radiation-induced MN-RET of human bone marrow cultured in the 3D bioreactor, a synthetic bioculture system, and suggest that this model may serve as a promising tool for studying MN-RET formation in human bone marrow, thereby providing opportunities to study bone marrow genotoxicity testing, mitigating agent effects, and other conditions that are not ordinarily feasible to experimental manipulation in vivo.     

Proliferation and degeneration of circulating CFU-GM in 15 adult normal subjects. Can any ontogenetic relationship be deduced from these data?

We studied, in 15 normal adults, the "in vitro" proliferation and differentiation of circulating CFU-GM, in order to assess their implication in the processes that regulate the dynamic equilibrium of granulopoiesis, analogous to bone marrow CFU-GM, and to deduce by their growth behaviour the ontogenetic relationship between CFU-GM subpopulations in the circulating and bone marrow compartments respectively. We found that "in vitro" proliferation of circulating CFU-GM predominates over their degeneration. We believe circulating CFU-GM and bone marrow CFU-GM are not implicated in granulopoiesis regulation in the same manner, and that circulating CFU-GM are more immature than bone marrow CFU-GM when taking proliferation and GM-CSF response into account. One cannot ignore the hypothesis that cells that "in vivo" are quiescent, are recruited "in vitro" with GM-CSF. Finally we would like to draw attention to the parallelism between CFU-GM classification in types 1 and 2 using monoclonal antibodies to track surface antigens, and our classification obtained by using a new mathematical model that takes the "birth" and "dead" of cellular aggregates into account.     

GROWING BONES ON THE COMPUTER—SOME PITFALLS OF A COMPUTER SIMULATION OF THE EFFECTS OF RADIATION ON BONE GROWTH

The linear arrangement of the proliferating cells in the growth cartilage of bone provides an apparently simple system for studying the relationship of cell kinetics to a gross endpoint, i.e. bone growth. The system has been reduced to a simple computer model which follows the effects of radiation on the growth of fifty or one hundred individual cartilage columns. When the computed data were compared with experimental curves some factors were discovered that had been omitted from the original model. A revised model proved of limited value in sorting out the relative importance of competing kinetic parameters but indicated some directions for future experimental studies.     

No effects of levamisole on cytotoxic drug-induced changes of human granulopoiesis.

The effect of Levamisole on the human granulopoiesis was studied in patients randomized to receive, in addition to adjuvant chemotherapy for primary breast cancer, either no other treatment or additional unspecific immune therapy with Levamisole. The reaction of granulopoiesis to the cytostatic drugs, as characterized by changes of peripheral blood polymorphonuclear neutrophils (PMN), functional bone marrow granulocyte reserve, serial bone marrow cytology, and granulopoietic stem cells (CFU-C) in marrow and blood, was not affected by administration of Levamisole. The data support the concept that Levamisole has no direct effect on human bone marrow granulopoiesis, but that an allergic mechanism is involved in the pathogenesis of Levamisole-induced agranulocytosis. The expectation that Levamisole exerts a beneficial effect by stimulation of the granulopoiesis, as previously suggested for BCG and Corynebacterium parvum, could not be substantiated in our studies.     

Neutrophil marrow release. A system analysis.

Neutrophil marrow egress is governed by several processes. The most important are cell maturation, functional behavior of marrow sinusoids and humoral or neuro-vascular factors. Neutrophil release cannot be observed directly but is reflected in the size, cellular composition and kinetics of the nonproliferating pool of granulocytopoiesis in bone marrow and of blood neutrophil pool. These experimentally determined parameters were used as the basis of a mathematical model study. The model describes two catenated compartments, the nonproliferating pool of granulocytopoiesis in marrow and the total blood granulocyte pool. Cell transit from one pool to the other was assumed to be age-dependent. It was expressed by a positive sloping sigmoidal function that defines the egress potential fo the cells that increases with cell maturation. During maturation granulocytopoietic cells develop intense motility which determines the morphology of the cells on smears. Relationship between cell motility and its morphology was defined by functions determining the age-dependent probabilities of cell fixation as metamyelocytes, band- and segmented forms, respectively. The parameters of this model could be so adjusted that all experimental data were matched within experimental errors. Thus, qualitative and quantitative information on neutrophil marrow egress was obtained for normal and pathological states of granulocytopoiesis.     

The complex effect of granulocyte colony-stimulating factor on human granulopoiesis analyzed by a new physiologically-based mathematical model

Neutropenia, frequently a side effect of chemo- and radiotherapy, increases susceptibility to microbial infections and is a life-threatening condition. For realistically predicting drug treatment effects on granulopoiesis, we have constructed a new mathematical model of granulopoiesis in the bone marrow and in the peripheral blood, featuring cell cycle phase transition and detailed granulocyte-colony stimulating factor (G-CSF) pharmacokinetics (PK) and pharmacodynamics (PD), including intracellular second messenger. Using this model, in conjunction with clinical results, we evaluated the system parameters, implemented those in the model and successfully retrieved the results of several independent clinical experiments under a wide range of G-CSF regimens. Our results show that the introduction of G-CSF-controlled intracellular second messenger is indispensable for precise retrieval of the clinical results, and suggest that the half-life of this messenger varies between a single and multiple G-CSF administration schedules. In addition, our model provided reliable steady-state, as well as dynamic, estimations of human granulopoiesis parameters. These included an estimation of apoptosis index in the post-mitotic compartment, which corroborates previous results. At present the model is used for suggesting improved drug regimens.     

Granulopoiesis in long-term culture by marrow from mice with busulfan-induced chronic latent aplasia

Mice given high-dose busulfan therapy develop a chronic latent marrow aplasia characterized by normal peripheral blood neutrophil numbers, hematocrits and marrow cellularity but reduced numbers of pluripotent hemopoietic stem cells (CFU-s) and granulocyte-monocyte progenitor cells (CFU-gm). To study the pathogenesis of this lesion, bone marrow was propagated in long-term marrow cultures (LTMC). Small amounts of normal marrow readily established and sustained long-term granulopoiesis in vitro. In contrast, inocula of marrow from busulfan-treated animals containing three to five times as many stem and progenitor cells failed to establish long-term granulopoiesis in vitro. These results suggest that high-dose busulfan therapy produces a qualitative defect in either the hemopoietic stem cells, the stromal-forming elements, or both, rendering them incapable of establishing long-term granulopoiesis in vitro. Furthermore, mixing experiments employing normal and busulfan-damaged marrow demonstrate that this qualitative defect is not due to the emergence of a suppressor cell population. LTMC can show types of marrow damage not detectable by other techniques currently available and represent a powerful tool for studying latent bone marrow failure.     

Modeling kinetics of a large‐scale fed‐batch CHO cell culture by Markov chain Monte Carlo method

Markov chain Monte Carlo (MCMC) method was applied to model kinetics of a fed‐batch Chinese hamster ovary cell culture process in 5,000‐L bioreactors. The kinetic model consists of six differential equations, which describe dynamics of viable cell density and concentrations of glucose, glutamine, ammonia, lactate, and the antibody fusion protein B1 (B1). The kinetic model has 18 parameters, six of which were calculated from the cell culture data, whereas the other 12 were estimated from a training data set that comprised of seven cell culture runs using a MCMC method. The model was confirmed in two validation data sets that represented a perturbation of the cell culture condition. The agreement between the predicted and measured values of both validation data sets may indicate high reliability of the model estimates. The kinetic model uniquely incorporated the ammonia removal and the exponential function of B1 protein concentration. The model indicated that ammonia and lactate play critical roles in cell growth and that low concentrations of glucose (0.17 mM) and glutamine (0.09 mM) in the cell culture medium may help reduce ammonia and lactate production. The model demonstrated that 83% of the glucose consumed was used for cell maintenance during the late phase of the cell cultures, whereas the maintenance coefficient for glutamine was negligible. Finally, the kinetic model suggests that it is critical for B1 production to sustain a high number of viable cells. The MCMC methodology may be a useful tool for modeling kinetics of a fed‐batch mammalian cell culture process. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

MicroRNA Profiling in Human Neutrophils during Bone Marrow Granulopoiesis and In Vivo Exudation

The purpose of this study was to describe the microRNA (miRNA) expression profiles of neutrophils and their precursors from the initiation of granulopoiesis in the bone marrow to extravasation and accumulation in skin windows. We analyzed three different cell populations from human bone marrow, polymorphonuclear neutrophil (PMNs) from peripheral blood, and extravasated PMNs from skin windows using the Affymetrix 2.0 platform. Our data reveal 135 miRNAs differentially regulated during bone marrow granulopoiesis. The majority is differentially regulated between the myeloblast/promyelocyte (MB/PM) and myelocyte/metamyelocyte (MC/MM) stages of development. These 135 miRNAs were divided into six clusters according to the pattern of their expression. Several miRNAs demonstrate a pronounced increase or reduction at the transition between MB/PM and MC/MM, which is associated with cell cycle arrest and the initiation of terminal differentiation. Seven miRNAs are differentially up-regulated between peripheral blood PMNs and extravasated PMNs and only one of these (miR-132) is also differentially regulated during granulopoiesis. The study indicates that several different miRNAs participate in the regulation of normal granulopoiesis and that miRNAs might also regulate activities of extravasated neutrophils. The data present the miRNA profiles during the development and activation of the neutrophil granulocyte in healthy humans and thus serves as a reference for further research of normal and malignant granulocytic development.     

Cytoskeletal influences on nuclear shape in granulocytic HL-60 cells

Background  

During granulopoiesis in the bone marrow, the nucleus differentiates from ovoid to lobulated shape. Addition of retinoic acid (RA) to leukemic HL-60 cells induces development of lobulated nuclei, furnishing a convenient model system for nuclear differentiation during granulopoiesis. Previous studies from our laboratory have implicated nuclear envelope composition as playing important roles in nuclear shape changes. Specifically noted were: 1) a paucity of lamins A/C and B1 in the undifferentiated and RA treated cell forms; 2) an elevation of lamin B receptor (LBR) during induced granulopoiesis.     

In vivo regulation of granulocyte protluction in acute inflammation

In order to study in vivo the enhanced granulopoiesis that occurs during acute inflammation, 1-3 sterile metallic copper rods were inserted subcutaneously into mice either at the same place (one abscess) or at different sites (multiple abscesses). Diffusion chambers filled with bone marrow cells were implanted intraperitoneally for 3 days. When a single abscess was created, the granulocytic content of the diffusion chamber increased similarly whatever the number of inserted copper rods. However, there was a direct relationship between the number of abscesses and the number of granulocytic cells harvested from the diffusion chambers. In order to investigate the role of T-lymphocytes in the production of diffusible stimulating factors that act on diffusion chamber granulopoiesis, cyclosporin A (CyA) was given to the mice with implanted copper rods. CyA abrogated the induced enhancement of CFU-S, CFU-GM and mature granulocyte numbers inside the diffusion chamber. The stimulatory effect of inflammation on diffusion chamber granulopoiesis was not observed in T-lymphocyte-deficient nude mice. These data suggest that in vivo stimulation of granulopoiesis is related to the level of inflammation, and that this effect requires the functional integrity of T-lymphocytes.     

A histomorphometric analysis of trephine biopsies of bone marrow from 65 patients with chronic myeloid leukemia. Classification of patients into subgroups with different survival patterns

A histomorphometric (planimetric) study was performed on trephine biopsies of the bone marrow taken at presentation from 65 patients (31 males and 34 females, with a median age of 48 years) with chronic myeloid leukemia (CML). Specimens from 20 patients (9 males and 11 females, with a median age of 53 years) without any hematologic disorders served as controls. Of the various histologic variables tested, only the counts of neutrophilic granulocytes per 1 sq mm, the ratio of granulocytopoiesis to megakaryopoiesis and the density of reticulin (argyrophilic) fibers revealed a significant correlation with the prognosis. The CML patients were separated into two groups with different survival patterns. Group I (34 patients with a median survival of 24 months) mostly contained cases with the so-called "megakaryocytic subtype" of CML, which is accompanied by variable degrees of fibrosis; group II (31 patients with a median survival of 36 months) mainly contained cases with the "granulocytic subtype," which is not accompanied by myelofibrosis. Among the morphometric parameters, a positive correlation existed between the megakaryocyte count and the reticulin fiber density, which underlines the important role of that cell lineage in fibrillogenesis. There were multiple interrelationships between the histomorphometric variables and the laboratory data. Consequently, multivariate regression methods (using Cox's proportional hazards model) were applied to assess the relative predictive value of the patient characteristics for survival. The derived prognostic model divided the patients into two risk groups, with median survivals of 14 and 41 months, respectively. In order of their entry into the regression model, these variables were percentage of neutrophils in the differential blood count, amount of granulopoiesis, liver size, percentage of peripheral myeloblasts and density of reticulin fibers in the bone marrow. In comparing the two patient groups, based on bone marrow histomorphometric parameters, this model revealed that two of those factors (amount of granulopoiesis and density of reticulin fibers) had a significant correlation with the prognosis.     

AN ANALYSIS OF BONE MARROW ERYTHROPOIESIS IN THE MOUSE

Following the model of the erythropoietic system developed in the rat by Tarbutt and Blackett, the authors have carried out a kinetic analysis of bone marrow erythropoiesis in the mouse.
Using ⁵⁹Fe labelling techniques the size of the recognizable precursor cells and of the functional cell compartments have been estimated, while the flow-rate from the unrecognized precursor cells to the recognizable cells and from the latter compartment to the circulating erythrocytes have been evaluated by ⁵⁵Fe autoradiography.
Differences in the kinetic parameters of the erythropoietic mouse bone marrow compared with the rat bone marrow are reported, whose interpretation has required a more detailed analysis of the original model.     

k-Cone analysis: determining all candidate values for kinetic parameters on a network scale

The absence of comprehensive measured kinetic values and the observed inconsistency in the available in vitro kinetic data has hindered the formulation of network-scale kinetic models of biochemical reaction networks. To meet this challenge we present an approach to construct a convex space, termed the k-cone, which contains all the allowable numerical values of the kinetic constants in large-scale biochemical networks. The definition of the k-cone relies on the incorporation of in vivo concentration data and a simplified approach to represent enzyme kinetics within an established constraint-based modeling approach. The k-cone approach was implemented to define the allowable combination of numerical values for a full kinetic model of human red blood cell metabolism and to study its correlated kinetic parameters. The k-cone approach can be used to determine consistency between in vitro measured kinetic values and in vivo concentration and flux measurements when used in a network-scale kinetic model. k-Cone analysis was successful in determining whether in vitro measured kinetic values used in the reconstruction of a kinetic-based model of Saccharomyces cerevisiae central metabolism could reproduce in vivo measurements. Further, the k-cone can be used to determine which numerical values of in vitro measured parameters are required to be changed in a kinetic model if in vivo measured values are not reproduced. k-Cone analysis could identify what minimum number of in vitro determined kinetic parameters needed to be adjusted in the S. cerevisiae model to be consistent with the in vivo data. Applying the k-cone analysis a priori to kinetic model development may reduce the time and effort involved in model building and parameter adjustment. With the recent developments in high-throughput profiling of metabolite concentrations at a whole-cell scale and advances in metabolomics technologies, the k-cone approach presented here may hold the promise for kinetic characterization of metabolic networks as well as other biological functions at a whole-cell level.     

Stimulating factors and cell recruitment in murine bone marrow stem cells and EMT6 tumours.

The role of a stimulating factor in cell recruitment and the kinetics of its secretion were investigated by in vivo and in vitro techniques. The association of these two methods made it possible to demonstrate that a non-cycling population liberates a factor which in turn stimulates quiescent bone marrow stem cells into DNA synthesis. Moreover, it seems that undamaged cells are capable of secreting this factor. A stimulating factor responsible for cell recruitment was also demonstrated in an experimental EMT6 tumour and the kinetics of its secretion reported.     

Analysis of Chinese hamster ovary cell metabolism through a combined computational and experimental approach

Optimization of cell culture processes can benefit from the systematic analysis of experimental data and their organization in mathematical models, which can be used to decipher the effect of individual process variables on multiple outputs of interest. Towards this goal, a kinetic model of cytosolic glucose metabolism coupled with a population-level model of Chinese hamster ovary cells was used to analyse metabolic behavior under batch and fed-batch cell culture conditions. The model was parameterized using experimental data for cell growth dynamics, extracellular and intracellular metabolite profiles. The results highlight significant differences between the two culture conditions in terms of metabolic efficiency and motivate the exploration of lactate as a secondary carbon source. Finally, the application of global sensitivity analysis to the model parameters highlights the need for additional experimental information on cell cycle distribution to complement metabolomic analyses with a view to parameterize kinetic models.