The paper is devoted to mathematical modelling of erythropoiesis, production of red blood cells in the bone marrow. We discuss intra-cellular regulatory networks which determine self-renewal and differentiation of erythroid progenitors. In the case of excessive self-renewal, immature cells can fill the bone marrow resulting in the development of leukemia. We introduce a parameter characterizing the strength of mutation. Depending on its value, leukemia will or will not develop. The simplest...