Divergent erythroid megakaryocyte fates in Blvrb-deficient mice establish non-overlapping cytoprotective functions during stress hematopoiesis

Cytoprotective mechanisms of heme oxygenases function by derivatizing heme to generate carbon monoxide, ferrous iron, and isomeric biliverdins, followed by rapid NAD(P)H-dependent biliverdin reduction to the antioxidant bilirubin using two non-overlapping biliverdin reductases that display biliverdin isomer-restricted redox activity. Although cytoprotective functions of heme oxygenases are widely recognized, concomitant effects of downstream biliverdin reductases remain incomplete. A computational model predicated on murine hematopoietic single-cell transcriptomic data identified Blvrb as a biological driver linked to the tumor necrosis factor stress pathway as a predominant source of variation defining hematopoietic cell heterogeneity. In vivo studies using Blvrb-deficient mice established the dispensable role of Blvrb in steady-state hematopoiesis, although model validation using aged Blvrb-deficient mice established an important cytoprotective function in stress hematopoiesis with dichotomous megakaryocyte-biased hematopoietic recovery. Defective stress erythropoiesis was evident in Blvrb^−/− spleens and in bone marrow erythroid development, occurring in conjunction with defective lipid peroxidation as a marker of oxidant mishandling. Cell autonomous effects on megakaryocyte lineage bias were documented using multipotential progenitor assays. These data provide the first physiological function of murine Blvrb in a non-redundant pathway of stress cytoprotection. Divergent effects on erythroid/megakaryocyte lineage speciation impute a novel redox-regulated mechanism for lineage partitioning.