Dans ce travail de recherche, nous sommes intéresses par la modélisation de l'hématopoïèse. Les cellules souches hématopoïétiques (CSH) sont des cellules indifférenciées de la moelle osseuse, possédant la capacité de se renouveler et de se différencier (pour la production des globules rouges, globules blancs et les plaquettes). Le processus de l'hématopoïèse souvent révèle des irrégularités qui causent les maladies hématologiques. En modélisant la leucémie myéloide chronique (LMC), une maladie hématologique fréquente, nous représentons l'hématopoïèse des cellules normales et cancéreuses par un système d'équations différentielles ordinaires (EDO). L'homéostasie des cellules normales et différente de l'homéostasie des cellules cancéreuses, et dépend de quelques lignées des cellules normales et cancéreuses. Nous analysons la dynamique globale du modèle pour obtenir les conditions de régénération de l'hématopoïèse ou bien la persistance de la LMC. Nous démontrons aussi que la coexistence des cellules normales et cancéreuses ne peut avoir lieu pour longtemps. Imatinib est un traitement de base de la LMC, avec un dosage variant de 400 à 1000 mg par jour. Certains patients présentent des réponses différentes à la thérapie, pouvant être hématologique, cytogénétique et moléculaire. La thérapie échoue dans deux cas: le patient demande un temps plus long pour réagir, alors il s'agit d'une réponse suboptimale; ou bien le patient résiste après une bonne réponse initiale. Pour déterminer le dosage optimal, nécessaire à la réduction des cellules cancéreuses, nous représentons les effets de la thérapie par un problème de contrôle optimal. Notre but est de minimiser le cout du traitement et le nombre des cellules cancéreuses. La réponse suboptimale, la résistance et le rétablissement sont alors obtenus suivant l'influence de l'imatinib sur les taux de division et de mortalité des cellules cancéreuses. Nous étudions par ailleurs l'hématopoïèse selon un modèle structuré en age, décrivant l'évolution des CSH normales et cancéreuses. Nous démontrons que le taux de division des CSH cancéreuses joue un rôle important dans la détermination du contrôle optimal. En contrôlant la croissance des cellules normales et cancéreuses avec compétition inter spécifique, nous démontrons que le dosage optimal dépend de l'homéostasie des CSH cancéreuses. / Modelling hematopoiesis represents a feature of our research. Hematopoietic stem cells (HSC) are undifferentiated cells, located in bone marrow, with unique abilities of self-renewal and differentiation (production of white cells, red blood cells and platelets).The process of hematopoiesis often exhibits abnormalities causing hematological diseases. In modelling Chronic Myeloid Leukemia (CML), a frequent hematological disease, we represent hematopoiesis of normal and leukemic cells by means of ordinary differential equations (ODE). Homeostasis of normal and leukemic cells are supposed to be different and depend on some lines of normal and leukemic HSC. We analyze the global dynamics of the model to obtain the conditions for regeneration of hematopoiesis and persistence of CML. We prove as well that normal and leukemic cells can not coexist for a long time. Imatinib is the main treatment of CML, with posology varying from 400 to 1000 mg per day. Some affected individuals respond to therapy with various levels being hematologic, cytogenetic and molecular. Therapy fails in two cases: the patient takes a long time to react, then suboptimal response occurs; or the patient resists after an initial response. Determining the optimal dosage required to reduce leukemic cells is another challenge. We approach therapy effects as an optimal control problem to minimize the cost of treatment and the level of leukemic cells. Suboptimal response, resistance and recovery forms are obtained through the influence of imatinib onto the division and mortality rates of leukemic cells. Hematopoiesis can be investigated according to age of cells. An age-structured system, describing the evolution of normal and leukemic HSC shows that the division rate of leukemic HSC plays a crucial role when determining the optimal control. When controlling the growth of cells under interspecific competition within normal and leukemic HSC, we prove that optimal dosage is related to homeostasis of leukemic HSC.

(ENGLISH) Hematopoiesis is dependent on the actions of hematopoietic stem cells (HSCs). This process is tightly controlled through a complex array of extrinsic and intrinsic factors. Even though the hematopoiesis seems to be well conserved across the disparate vertebrate animals, erythroid and thrombocytic differentiation have changed during the evolution of mammals. Specifically, adult mammalian red blood cells have the unique feature of being enucleated, and mammalian thrombocytes are not individual cells, but fragments of megakaryocytes, instead. It is likely that these enhancements provided a survival advantage to early mammalian species; however, they also bring up the question of evolutionary origin of these cells that studied using zebrafish (Danio rerio) model. First, it was necessary to generate a toolbox of a recombinant cytokines and optimized culture media that allowed us to manipulate zebrafish hematopoietic cells ex vivo in liquid and clonal cultures. Interestingly, teleost species underwent an extra duplication event during their evolution and as a result, two copies (paralogs) of some of the genes are present in zebrafish. This was also the case for majority of the cytokines from our toolbox and here, we provide functional characterization of these paralogs. Strikingly, our results...

Indiana University-Purdue University Indianapolis (IUPUI) / Juvenile Myelomonocytic Leukemia (JMML) is rare, fatal myeloproliferative disease (MPD) affecting young children, and is characterized by expansion of monocyte lineage cells and hypersensitivity to Granulocyte Macrophage-Colony Stimulating Factor (GM-CSF) stimulation. JMML is frequently associated with gain-of-function mutations in the PTPN11 gene, which encodes the protein tyrosine phosphatase, Shp2. Activating Shp2 mutations are known to promote hyperactivation of the Ras-Erk signaling pathway, but Akt is also observed to have enhanced phosphorylation, suggesting a potential role for Phosphatidylinositol-3-Kinase (PI3K)-Akt signaling in mutant Shp2-induced GM-CSF hypersensitivity and leukemogenesis. Having demonstrated that Class IA PI3K is hyperactivated in the presence of mutant Shp2 and contributes to GM-CSF hypersensitivity, I hypothesized the hematopoietic-specific Class IA PI3K catalytic subunit p110δ is a crucial mediator of mutant Shp2-induced PI3K hyperactivation and GM-CSF hypersensitivity in vitro and MPD development in vivo. I crossed gain-of-function mutant Shp2 D61Y inducible knockin mice, which develop fatal MPD, with mice expressing kinase-dead mutant p110δ D910A to evaluate p110δ’s role in mutant Shp2-induced GM-CSF hypersensitivity in vitro and MPD development in vivo. As a comparison, I also crossed Shp2 D61Y inducible knockin mice with mice bearing inducible knockout of the ubiquitously expressed Class IA PI3K catalytic subunit, p110α. I found that genetic interruption of p110δ, but not p110α, significantly reduced GM-CSF-stimulated hyperactivation of both the Ras-Erk and PI3K-Akt signaling pathways, and as a consequence, resulted in reduced GM-CSF-stimulated hyper-proliferation in vitro. Furthermore, I found that mice bearing genetic disruption of p110δ, but not p110α, in the presence of gain-of-function mutant Shp2 D61Y, had on average, smaller spleen sizes, suggesting that loss of p110δ activity reduced MPD severity in vivo. I also investigated the effects of three PI3K inhibitors with high specificity for p110δ, IC87114, GDC-0941, and GS-9820 (formerly known as CAL-120), on mutant Shp2-induced GM-CSF hypersensitivity. These inhibitors with high specificity for p110δ significantly reduced GM-CSF-stimulated hyperactivation of PI3K-Akt and Ras-Erk signaling and reduced GM-CSF-stimulated hyperproliferation in cells expressing gain-of-function Shp2 mutants. Collectively, these findings show that p110δ-dependent PI3K hyperactivation contributes to mutant Shp2-induced GM-CSF hypersensitivity and MPD development, and that p110δ represents a potential novel therapeutic target for JMML.

Indiana University-Purdue University Indianapolis (IUPUI) / Skeletal manifestations including short stature, osteoporosis, kyphoscoliosis, and tibial dysplasia cumulatively affect approximately 70% of patients with neurofibromatosis type 1 (NF1). Tibial pseudarthrosis, the chronic non-union of a spontaneous fracture, is a debilitating skeletal malady affecting young children with NF1. These non-healing fractures respond poorly to treatment and often require amputation of the affected limb due to limited understanding of the causative mechanisms. To better understand the cellular and molecular pathogenesis of these osseous defects, we have established a new mouse model which recapitulates a spectrum of skeletal pathologies frequently observed in patients with NF1. Nf1flox/-;Col2.3Cre mice, harboring Nf1 nullizygous osteoblasts on a Nf1+/- background, exhibit multiple osseous defects which are closely reminiscent of those found in NF1 patients, including runting (short stature), bone mass deficits, spinal deformities, and tibial fracture non-union. Through adoptive bone marrow transfer studies, we have demonstrated that the Nf1 haploinsufficient hematopoietic system pivotally mediates the pathogenesis of bone loss and fracture non-union in Nf1flox/-;Col2.3Cre mice. By genetic ablation of a single Nf1 allele in early myeloid development, under the control of LysMCre, we have further delineated that Nf1 haploinsufficient myeloid progenitors and osteoclasts are the culprit lineages mediating accelerated bone loss. Interestingly, conditional Nf1 haploinsufficiency in mature osteoclasts, induced by CtskCre, was insufficient to trigger enhanced lytic activity. These data provide direct genetic evidence for Nf1’s temporal significance as a gatekeeper of the osteoclast progenitor pool in primitive myelopoiesis. On the molecular level, we found that transforming growth factor-beta1 (TGF-β1), a primary mediator in the spatiotemporal coupling of bone remodeling, is pathologically overexpressed by five- to six- fold in both NF1 patients and in mice. Nf1 deficient osteoblasts, the principal source of TGF-β1 in the bone matrix, overexpress TGF-β1 in a gene dosage dependent fashion. Moreover, p21Ras dependent hyperactivation of the Smad pathway accentuates responses to pathological TGF-β1 signals in Nf1 deficient bone cells. As a proof of concept, we demonstrate that pharmacologic TβRI kinase inhibition can rescue bone mass defects and prevent tibial fracture non-union in Nf1flox/-;Col2.3Cre mice, suggesting that targeting TGF-β1 signaling in myeloid lineages may provide therapeutic benefit for treating NF1 skeletal defects.

Stem cells differentiate through an organized hierarchy of intermediate cell types to terminally differentiated cell types. This process is largely guided by master transcriptional regulators, but it also depends on the expression of many other types of genes. The discrete cell types in the differentiation hierarchy are often identified based on the expression or non-expression of certain marker genes. Historically, these have often been various cell-surface proteins, which are fairly easy to assay biochemically but are not necessarily causative of the cell type, in the sense of being master transcriptional regulators. This raises important questions about how gene expression across the whole genome controls or reflects cell state, and in particular, differentiation hierarchies. Traditional approaches to understanding gene expression patterns across multiple conditions, such as principal components analysis or K-means clustering, can group cell types based on gene expression, but they do so without knowledge of the differentiation hierarchy. Hierarchical clustering and maximization of parsimony can organize the cell types into a tree, but in general this tree is different from the differentiation hierarchy. Using hematopoietic differentiation as an example, we demonstrate how many genes other than marker genes are able to discriminate between different branches of the differentiation tree by proposing two models for detecting genes that are up-regulated or down-regulated in distinct lineages. We then propose a novel approach to solving the following problem: Given the differentiation hierarchy and gene expression data at each node, construct a weighted Euclidean distance metric such that the minimum spanning tree with respect to that metric is precisely the given differentiation hierarchy. We provide a set of linear constraints that are provably sufficient for the desired construction and a linear programming framework to identify sparse sets of weights, effectively identifying genes that are most relevant for discriminating different parts of the tree. We apply our method to microarray gene expression data describing 38 cell types in the hematopoiesis hierarchy, constructing a sparse weighted Euclidean metric that uses just 175 genes. These 175 genes are different than the marker genes that were used to identify the 38 cell types, hence offering a novel alternative way of discriminating different branches of the tree. A DAVID functional annotation analysis shows that the 175 genes reflect major processes and pathways active in different parts of the tree. However, we find that there are many alternative sets of weights that satisfy the linear constraints. Thus, in the style of random-forest training, we also construct metrics based on random subsets of the genes and compare them to the metric of 175 genes. Our results show that the 175 genes frequently appear in the random metrics, implicating their significance from an empirical point of view as well. Finally, we show how our linear programming method is able to identify columns that were selected to build minimum spanning trees on the nodes of random variable-size matrices.