Heme-stress activated NRF2 skews fate trajectories of bone marrow cells from dendritic cells towards red pulp-like macrophages in hemolytic anemia

Heme is an erythrocyte-derived toxin that drives disease progression in hemolytic anemias, such as sickle cell disease. During hemolysis, specialized bone marrow-derived macrophages with a high heme-metabolism capacity orchestrate disease adaptation by removing damaged erythrocytes and heme-protein complexes from the blood and supporting iron recycling for erythropoiesis. Since chronic heme-stress is noxious for macrophages, erythrophagocytes in the spleen are continuously replenished from bone marrow-derived progenitors. Here, we hypothesized that adaptation to heme stress progressively shifts differentiation trajectories of bone marrow progenitors to expand the capacity of heme-handling monocyte-derived macrophages at the expense of the homeostatic generation of dendritic cells, which emerge from shared myeloid precursors. This heme-induced redirection of differentiation trajectories may contribute to hemolysis-induced secondary immunodeficiency. We performed single-cell RNA-sequencing with directional RNA velocity analysis of GM-CSF-supplemented mouse bone marrow cultures to assess myeloid differentiation under heme stress. We found that heme-activated NRF2 signaling shifted the differentiation of bone marrow cells towards antioxidant, iron-recycling macrophages, suppressing the generation of dendritic cells in heme-exposed bone marrow cultures. Heme eliminated the capacity of GM-CSF-supplemented bone marrow cultures to activate antigen-specific CD4 T cells. The generation of functionally competent dendritic cells was restored by NRF2 loss. The heme-induced phenotype of macrophage expansion with concurrent dendritic cell depletion was reproduced in hemolytic mice with sickle cell disease and spherocytosis and associated with reduced dendritic cell functions in the spleen. Our data provide a novel mechanistic underpinning of hemolytic stress as a driver of hyposplenism-related secondary immunodeficiency.


INTRODUCTION
The hemolytic destruction of red blood cells (RBCs) generates large quantities of cell-free hemoglobin and heme, which accumulate in the blood and tissues. These RBC-derived toxins drive morbidity and mortality in hemolytic anemias via a number of pathophysiological pathways that cause nitric oxide depletion, lipid peroxidation, and heme engagement of cellular signaling receptors [1][2][3][4]. While primary heme toxicities are localized in the cardiovascular system and kidneys, more recent research has identified heme as a context-dependent stimulator and suppressor of immunity and inflammation [3][4][5]. Immunosuppressive functions of heme may contribute towards a poorly understood secondary immunodeficiency, called hyposplenism, which is a critical contributor to mortality in patients with genetic hemolytic anemia [6][7][8][9].
Hemolytic anemias such as sickle cell disease (SCD), thalassemia, and spherocytosis are the most frequent monogenic diseases worldwide and have exerted strong evolutionary pressure on physiological adaptation enhancing disease tolerance to anemia and heme-stress [10,11]. A key role of this adaptation is provided by mononuclear phagocytes in the liver and spleen. These macrophages remove damaged RBCs and toxic heme-protein complexes from the blood, which limits hemolysis-related tissue toxicity and ultimately sustains iron metabolism and erythropoiesis [12][13][14][15][16][17]. However, chronic heme stress, such as that in hemolytic anemias, is noxious for erythrophagocytic macrophages, and their continuous replacement by bone marrow (BM)-derived monocytes is thus needed to maintain homeostasis [18][19][20][21][22][23]. This replacement is driven by a regulatory network that includes the heme-activated transcription factor Spi-C [22] and PPARγ [24]. The antioxidant and iron-recycling functions are further supported by ATF1 [25] and the basic leucine zipper protein transcription factor NFE2L2 (NRF2), which is known as the master regulator of the cellular response to oxidative stress [26].
Like monocyte-derived macrophages, dendritic cells (DCs) emerge from myeloid-rooted BM progenitors. DCs are the professional antigen-presenting cells and are essential for the initiation of effector T cell responses [27]. Despite their key role at the interface of innate and adaptive immunity, there is a lack of research defining the influence of hemolysis and heme stress on DCs and how a potential deregulation of DC homeostasis could contribute to hemolysis-associated immunodeficiency. Similar to heme-handling macrophages in the hemolytic spleen, DCs are short-lived cells that continuously emerge from the BM and differentiate into mature antigen-presenting cells in tissues [27][28][29]. Monocyte-derived macrophages and DCs are progenies of common myeloid BM-precursors, which implies a potential competition of macrophage and DC generation [30][31][32]. Based on this prior knowledge, we reasoned that adaptation to chronic heme stress may progressively shift the differentiation trajectories of myeloid BM progenitors to expand the capacity of hemehandling macrophages at the expense of homeostatic DC production. Although such a hematopoietic deviation may support adaptive physiology in hemolytic stress, the resulting deficiency in DCs may contribute to hyposplenism. Therefore, to explore how heme-stress could perturbate macrophage and DC homeostasis, we performed time-resolved single-cell RNA-sequencing (scRNA-seq) experiments with mouse BM cultures. We found that heme redirected the differentiation trajectory of BM cells toward antioxidant and iron-recycling macrophages rather than DCs via an NRF2-dependent signaling pathway. Heme stress in mice with SCD replicated this phenotype in vivo and showed the expansion of heme-handling macrophages and a concomitant reduction in DCs with an impaired capacity to prime antigen-specific T cells. Our data reveal a new regulatory function of heme in the immune system and points to a novel mechanistic understanding of secondary immunodeficiency in hemolytic anemia.

RESULTS
Heme exposure suppressed the phenotype of functional DCs in mouse BM cultures To investigate the impact of heme-exposure on myeloid cell differentiation trajectories in vitro, we cultured lineage-depleted (Lin − ) BM cells from C57BL/6 mice supplemented with recombinant M-CSF or GM-CSF, which are known to induce BM cells into macrophages and DCs, respectively. Within 7 days, these M-CSF-and GM-CSF BM cultures generated populations of CD11b + /CD115 + / MHC-class II low macrophages and CD11c + /CD115 − /MHC-class II high DC-like cells, respectively (Supplementary Fig. 1A) [30,33]. We characterized the potential heme-handling function of these cells by measuring the uptake of SnMP, a fluorescent heme-analog, using a spectral flow analyzer. In contrast to the DC-like cells, we detected a strong SnMP-typical fluorescence signal in the macrophages, suggesting functional heterogeneity ( Supplementary Fig. 1B).
To fully characterize the cell-type heterogeneity and changes induced by heme under the two culture conditions, we exposed M-CSF BM cultures and GM-CSF BM cultures to heme-albumin (heme) or albumin (vehicle control) and then performed a multiplexed scRNA-seq experiment (Fig. 1A). Cells under each condition were labeled with DNA-barcoded antibodies, pooled and processed for sequencing. Figure 1A shows a uniform manifold approximation and projection (UMAP) of the 4852 cells merged across the four treatment conditions (growth factor × treatment). The 15 clusters were assigned to cell types using an unbiased approach with an algorithm that performs a per-cluster search for populations with similar gene expression in a database of >15 million annotated single-cell transcriptomes [34] (Supplementary Fig. 1C). The M-CSF BM cultures generated only macrophages, with minor database hits consistent with monocytes and microglia. In contrast, the GM-CSF BM cultures were more heterogeneous and consisted of macrophages, DCs and neutrophils (Fig. 1B). Split UMAP plots for the four conditions showed clear segregation of the GM-CSF and M-CSF BM cultures and also of the control and heme-exposed cells, which implies a strong effect of heme-exposure on the overall gene expression profile (Fig. 1C).
The key strength of this experiment was the precise definition of how the presence of heme changed the cell types that were generated. In the M-CSF BM cultures all cells were assigned a macrophage phenotype regardless of whether heme was present in the culture or not. In contrast and unexpectedly, in the GM-CSF BM cultures, we found that the prevailing DC population represented by clusters 12 and 13 completely disappeared in the presence of heme, whereas a small population exhibiting neutrophil-specific gene expression emerged in only the hemetreated cells (cluster 14) (Fig. 1C, Supplementary Fig. 1C). The expression of canonical marker genes for DCs [35,36], macrophages, and neutrophils is shown in Supplementary Fig. 1D.
To infer potential functional implications from the gene expression phenotype induced by heme, we performed gene-set enrichment analysis (GSEA) of differentially expressed genes identified from a pseudo-bulk analysis of heme-treated versus control samples ( Fig. 1D and Supplementary Fig. 2A). In both GM-CSF and M-CSF BM cultures, heme strongly induced the hallmark gene set for the reactive oxygen species (ROS) pathway including the antioxidant genes Gclm, Prdx1, Gsr, Txnrd1, and Cat (Fig. 1E). We also found that heme enhanced expression of the ironrecycling machinery, defined by the heme transporter Slc48a1, the heme-degradation enzyme heme oxygenase (Hmox1), and genes related to iron storage (Ftl1, Fth1) and iron export (Slc40a1) (Fig. 1E). In contrast, we found consistent suppression (i.e., negative enrichment) of genes classified as IFN alpha and IFN gamma responses in the heme-treated compared to the control cells ( Fig. 1D and Supplementary Fig. 2A). Collectively, these GSEA results suggest that heme enhances antioxidant and iron-recycling functions, while specific immune functions may be suppressed.
The defining immune function of DCs is their capacity for antigen-specific priming of naive T cells [27,37]. Therefore, to functionally validate our scRNA-seq study, we analyzed the activation and proliferation of T cells isolated from OVA-specific recombinant T cell receptor transgenic mice (OT-2 and OT-1 mice) in cocultures with antigen-pulsed BM cells by evaluating carboxyfluorescein succinimidyl ester (CFSE) dilution and cytokine production. We performed these experiments under the assumption that the observed T cell response would directly reflect the output of mature and functionally competent DCs in the different cell cultures. Importantly, before the coculture with T cells was initiated, the BM cultures were washed to avoid direct interference of heme with T cell activation and proliferation. Consistent with our phenotype-based definition derived from the scRNA-seq data (Fig. 1C), coculture of OT-2 CD4 + T cells with OVA 323-339 peptidepulsed GM-CSF BM cultures, but not M-CSF BM cultures resulted in the appearance of multiple CD69 + T cell populations with progressively reduced (i.e., diluted) CFSE fluorescence, indicating T cell activation and division (Fig. 1F, G). This activation was accompanied by the secretion of IL-2, which was dose-dependent in respect to the OVA 323-339 peptide concentration used to pulse BM cultures (Fig. 1H). Heme exposure completely suppressed the CD4 + T cell-priming capacity of GM-CSF BM cultures (Fig. F, G, H), which was consistent with the heme-suppressed DC phenotype (Fig. 1D).
The ability of heme-treated GM-CSF BM culture to cross-present ovalbumin-derived antigen to CD8 + T was also impaired as indicated by the reduced proliferation of OT-1 cells in coculture experiments ( Supplementary Fig. 1E).
To verify that this effect of heme is maintained in vivo, we adoptively transferred CFSE-labeled OT-2 CD4 + T cells and OVA 323-339 peptide-pulsed GM-CSF BM cells that were cultured with heme or vehicle into CD45.1 Rag1 −/− mice and analyzed the proliferation of the T cells in the spleen after 3 days by flow cytometry. In this experiment, DC retained a strong antigen-presentation capacity that was profoundly suppressed by heme Fig. 1 Heme exposure suppressed CD11c + /MHC class II high /CD115 (Csf1R) − DC phenotype and function in GM-CSF mouse BM cultures. A Left: Outline of the multiplexed scRNA-seq experiment. M-CSF-and GM-CSF BM cultures were stimulated with heme (300 µM) or vehicle on day 3. On day 7, the cells were harvested, tagged with oligo-barcoded antibodies, and pooled into a single scRNA-seq sample for sequencing. Demultiplexing and reassignment of treatment groups were performed after quality assessment, filtering, and cell-type attribution. Right: UMAP plot showing cells colored by DNA-barcoded antibodies (hashtag labeling). B UMAP plot showing cells colored by cell type. C Factorial UMAP plots of the dataset separated by culture supplement (M-CSF versus GM-CSF) and treatment condition (vehicle versus heme exposure). The data points are colored by cell type (orange = macrophages, purple = neutrophils, green = DCs). D GSEA of differentially expressed genes (heme exposure versus vehicle) in GM-CSF or M-CSF BM cultures. The heatmaps represent the running enrichment score per gene-set category (red for positive enrichment, blue for negative enrichment). Each row represents one gene-set category. E UMAP plots showing the expression levels of the top five enriched genes defining the ROS pathway and genes associated with iron recycling in the whole data set. F M-CSF-and GM-CSF BM cells were cultured with heme (300 µM) or vehicle. On day 7, the cells were harvested, pulsed with OVA 323-339 (1 µg/ ml), washed to remove any remaining heme, and then incubated with CFSE-labeled naive CD4 + T cells isolated from OT-2 mice in a 1/5 BM cells to T cells ratio. The proliferation of CD4 + T cells was assessed by evaluating the degree of CFSE dilution measured by flow cytometry after 3 days in coculture. The cells were also stained for the T cell activation marker CD69. Representative data for cells gated based on positive CD4 expression are shown. G Cumulative T cell proliferation data from three independent coculture experiments. H GM-CSF BM cells were cultured with heme (300 µM) or vehicle. On day 7, the cells were harvested, pulsed with OVA 323-339 (1 or 10 µg/ml), washed to remove any remaining heme, and then incubated with CFSE-labeled naive CD4 + T cells. The proliferation of CD4 + T cells was assessed after 3 days in coculture by flow cytometry. Cumulative T cell proliferation data from three independent coculture experiments are shown. The data in G, H are presented as the means ± SDs. Each dot represents one independent experiment. t-test (G); one-way ANOVA with Tukey's multiple comparison test (H); n. s. not significant, *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001, ****P ≤ 0.0001. Fig. 1F, G). In summary, these experiments indicated that heme-stress suppressed the capacity of GM-CSF BM cultures to induce activation and proliferation of antigenspecific T cells consistent with the disappearance of DCs.

(Supplementary
Collectively, based on these experimental results, we hypothesized that heme exposure redirects BM cell differentiation trajectories in favor of antioxidant and iron-recycling macrophages at the expense of functionally mature DCs. Furthermore, we found that cell-type heterogeneity in the GM-CSF BM cultures provides a unique opportunity to explore heme-deviated myeloid cell differentiation in a single culture dish, which prompted the use of only GM-CSF BM cultures in the remainder of the study. Heme skewed the cell differentiation fates of mouse BM cultures from DCs towards macrophages To capture the cell differentiation trajectories of GM-CSF BM cultures and the effect of heme in more detail, we extended our studies with another scRNA-seq experiment using a time-resolved scRNA-seq experiment, as outlined in Fig. 2A. Briefly, we cultured Lin − BM cells with GM-CSF, added heme or vehicle on day 3 of culture, and cultured the cells for an additional 24, 72, or 96 h. Using CellRank [38], we combined transcriptional similarity data with directional RNA velocity information to delineate three terminal differentiation states, or macrostates (Fig. 2B). CellSearch similarity matching attributed these macrostates to DCs, macrophages, and neutrophils ( Fig. 2C), which is consistent with the expression of canonical marker genes für these cell types shown in Fig. 2D. Additionally, we identified three transitional cell states consistent with progenitors and cells at an intermediate differentiation state (Fig. 2E). Figure 2F shows the UMAP plots obtained in the experiment with data from cells under each treatment condition and exposure time viewed separately. Over time, the transcriptome profiles progressed along the RNA velocity vectors calculated for the whole dataset (Fig. 2E), with the hemetreated and control cultures showing very different outcomes. Cells in the control culture initially progressed on day 3 along two trajectories towards neutrophils and DCs. At the end of the study, on day 7, the neutrophil arm had disappeared, leaving a prevailing population of DCs (brown) alongside an intermediate state population with marker gene expression consistent with immature DCs and macrophages/monocytes. In the heme-exposed BM culture, the propagation and collapse of the neutrophils were comparable to those in the control culture, but the final differentiation fate of the macrophage/DC population was strongly skewed away from DCs towards Spic + red pulp-like macrophages, with high expression of key heme and iron metabolism genes such as the heme transporter Slc48a1, Hmox1, ferritin genes (Fth1, Flt1), and the iron exporter Slc40a1 ( Supplementary Fig. 2B). This qualitative interpretation was supported by a quantitative analysis of the discrete cell state proportions (Fig. 2G) and by the day-wise cumulative absorption probabilities (Fig. 2H), which represented the probabilistic fate that any cell within the experiment would differentiate into one of the three terminal states representing neutrophils, macrophages, or DCs. Here, the highest absorption probability on day 7 suggested a preferential DC cell fate for cells in the control condition and a preferential macrophage cell fate for cells in the heme-treated condition. Individual cell absorption probabilities for all cells across timepoints and treatments are shown in Fig. 2I [38,39]. A GSEA of differentially expressed genes in heme versus control cells is presented for each time point in Supplementary Fig. 3.
Collectively, these results reveal that heme exposure skewed the probabilistic differentiation fates of GM-CSF BM cultures from DCs and towards macrophages with a red pulp macrophage-like gene expression profile supporting heme metabolism and iron recycling.
Suppression of DCs in heme-exposed BM cultures proceeds via the transcription factor NRF2 Next, we aimed to identify the signaling pathway through which heme suppresses the generation of mature DCs. We analysed the genome-wide chromatin accessibility in single cells by ATAC-seq (scATAC-seq) of GM-CSF supplemented BM cells that were exposed on day 3 of the culture to heme or vehicle for 6 h. Figure 3A shows UMAPs of 23,691 analyzed nuclei, which could be classified as neutrophil or macrophage/DC precursors based on the pseudo gene expression of cell-type specific genes. We subsequently performed a motif enrichment analysis on the differentially accessible regions in heme versus vehicle-treated BM cells across the whole dataset. We found the most enriched motif corresponded to the Nfe2l2 motif (MA0150.2), with a cell specific Nfe2l2 enrichment in the heme-exposed macrophage/DC precursors.
NRF2 activation proceeds via nuclear translocation of the transcription factor after dissociation from its cytoplasmic captor protein KEAP1. This model was consistent with a strong nuclear accumulation of NRF2 antigen in HEK293 cells that were treated with heme for 2, 4, and 6 h, suggesting that rapid NRF2 activation and nuclear translocation is a highly conserved response to hemestress ( Supplementary Fig. 4A).
To validate the function of NRF2, we performed a two-factorial multiplexed scRNA-seq experiment using GM-CSF BM cultures from Nrf2 −/− mice and wild-type littermates that were treated with heme or with the selective NRF2 activator ML-334 ( Supplementary Fig. 4B). As an initial screen, we compared differential gene expression data of the compound-treated (heme or ML-334) versus vehicle-treated wild-type cells after functional consolidation of the data by GSEA. Supplementary Fig. 4C demonstrates the high correlation of the ranked GSEA results for heme and ML-334, respectively. Under both treatment conditions, genes involved in heme metabolism and the ROS pathway were strongly upregulated, while IFN response genes were suppressed to the greatest degree, suggesting that selective NRF2 activation replicates the effect of heme on GM-CSF BM cultures. Figure 3B shows UMAP plots of the scRNA-seq experiment for the cell densities under each experimental condition separated by genotype and treatment. Expression heatmaps of canonical marker genes of mature DCs (Ccr7, Cd83), monocytesmacrophages (Csf1r), and neutrophils (Csf3r) and that of representative genes known to be regulated by NRF2 (Gclm, Nqo1, Scl7a11, and Hmox1) are highlighted in Fig. 3C, and these results are accompanied by a cumulative NRF2 target gene expression score, which compares the mean expression of the Nrf2 target genes from the TRRUST gene-set library to a randomly selected reference set. The data indicated that both heme and ML-334 activated Nrf2-related gene expression in wild-type cells and almost completely blocked the emergence of Ccr7 + Cd83 + Csf1r − DCs. Nrf2 deficiency rescued the generation of DCs in hemetreated and ML-334-treated cells, which suggests that a linked heme-NRF2 signaling pathway suppressed DC generation in the heme-exposed BM cultures. In addition, we found that the hemeinduced suppression of DC generation was more pronounced in Hmox1 knockout GM-CSF BM cell cultures ( Supplementary Fig. 4D). Hmox1 knockout cells accumulate intracellular heme, resulting in stronger induction of NRF2 target genes compared to wild-type BM cells ( Supplementary Fig. 4E). These studies confirmed that the DC-suppressive signaling was directly mediated by intracellular heme and not by the primary heme-catabolites.
Collectively, the scRNA-seq results indicate that heme-activated NRF2 suppresses the emergence of mature DCs from BM progenitors.
Cell-intrinsic heme-NRF2 signaling suppressed DC and functional antigen-presentation capacity in mouse BM cultures In the next step, we validated the function of the heme-NRF2 signaling pathway as a suppressor of DC generation by flow cytometry and T cell proliferation assays. For these studies we treated GM-CSF BM cultures with different concentrations of heme Fig. 2 Heme redirects the cell differentiation pathways in GM-CS mouse BM cultures from DCs towards macrophages. A Outline of the multiplexed scRNA-seq experiment. GM-CSF BM-cultures were initiated on day 0, 1, 3, and 4. Each culture was treated with heme (300 µM) or vehicle control (black arrows) on day 3 and cells were harvested on day 7 providing heme exposure periods of 0, 24, 72, and 96 h. After, the cells were harvested, tagged with DNA-barcoded antibodies, pooled, and processed for sequencing. After quality assessment, filtering, and cell-type attribution, the experiment was demultiplexed to reassign individual cells to their sample of origin. A UMAP plot of the scRNA-seq results highlighting the cells that were automatically selected by the CellRank algorithm as most representative of one of three, yet undefined, terminal differentiation states (1 brown, 2 blue, 3 violet). The terminal state-typical cell populations selected in B were used in a similarity search against >15 million annotated cells from published studies. Similarity matching was performed using CellSearch (Bioturing). For each macrostate (1, 2, or 3), the heatmap shows the number of high-quality matches for DCs, macrophages, and neutrophils that were then used for cell-type attribution. C Expression intensity projections of the scRNA-seq data showing canonical marker genes for cDCs, neutrophils, and macrophages. D UMAP plots of the scRNA-seq data separated by treatment status (control and heme). The superimposed arrows indicate RNA velocity vectors that were calculated for each treatment condition across time points. The colors represent the three terminal states and three additional transition states. The macrostate colors will be used for the remainder of the figure. E UMAP plots of the scRNA-seq data separated by heme treatment status and time point illustrate the treatment-dependent differentiation pathways across the four time points. F Stack bar charts illustrating discrete cell-type attributions as a proportion of the whole number of cells per sample, for each time and treatment condition. Numbers are percentages. G Cumulative absorption probabilities summarizing the probability of each cell reaching each terminal state defined as DCs, macrophages, and neutrophils given for each time point for both experimental conditions. H UMAP plots with the cells colored according to their probability of reaching the terminal state defined as DCs, macrophages, or neutrophils. Data of all treatment conditions (heme × time) are merged. and selective NRF2 activators for 4 days. The results demonstrated that the absence of the Nrf2 gene blocked the heme-induced loss of CD11c + /MHC class II high /CD115 (Csf1R) − DCs (Fig. 4A, B). Additionally, the two NRF2 activators ML-334 and RA-839 did not suppress differentiation into the DCs in Nrf2 −/− cells (Fig. 4C, D), confirming the target-specificity of the observed effect. DC generation was also reduced in Keap1 GM-CSF BM-cultures, which have constitutively active NRF2, independent of the presence of a pharmacological activator ( Supplementary Fig. 4F). We also performed heme treatment experiments with mixed cultures of Nrf2 −/− and wild-type BM cells in the same cell culture dish, exploiting the congenic CD45.1/CD45.2 system, which can be used with flow cytometry to discriminate cell origin [40] (Fig. 4E). The results showed that in a mixed culture heme selectively  The IL-2 concentration was measured using a Bio-Plex assay. The data in B and D, and H are presented as the means ± SDs. Each dot represents one independent experiment. The data in F are representative of two independent experiments, and the data in G are representative of six independent experiments. The data presented in B, D, and H were analyzed by one-way ANOVA with Tukey's multiple comparison test. n.s. not significant, *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001, ****P ≤ 0.0001.
Taken together, these data support the existence of a cellintrinsic heme-NRF2 signaling pathway that prevents the generation of functionally competent DCs.
Phenotypical and functional DC depletion in the spleen of heme-exposed mice with genetic hemolytic anemia or Hmox1 deficiency To investigate the potential involvement of heme-mediated effects on DCs in a pathophysiological context, we explored mouse models of diseases related to heme exposure. These investigations focused on the spleen, which has a dual function in mice during hemolysis, acting as a clearance organ for damaged RBCs and cell-free heme and as an extramedullary hematopoietic organ [41,42]. We thus expected that the spleen would mirror a heme-exposed BM-like microenvironment and replicate our cell culture model in vivo.
Intracellular heme levels are defined by the influx of extracellular heme into the cell (e.g., via erythrophagocytosis or the uptake of heme-protein complexes) and by intracellular heme degradation through heme oxygenases [43]. The first model investigated in this study represents a condition of cell-intrinsic heme excess and is induced by conditional Hmox1 knockout. In these model mice, we deleted Hmox1 by tamoxifen treatment, which led to the progressive depletion of CD11c + /MHC class II high DCs in the spleen. In contrast to the DCs, the pool of neutrophils was markedly enhanced 9 weeks after tamoxifen treatment (Fig. 5A), indicating that the observed depletion of DCs was not caused by nonspecific toxicity but rather by more specific hematopoietic deregulation.
The second model used, which represents cell-extrinsic heme excess, is a mouse model of SCD. Berkeley SCD mice exhibited pronounced hemolysis with lower hematocrit and more reticulocytes than their control littermates (Fig. 5B). Intravascular hemolysis in SCD causes considerable free heme exposure, which leads to the accumulation of heme-adducted albumin in the plasma (Fig. 5B) [44]. In the spleen of SCD mice spectral flow cytometry revealed expanded populations of F4/80 + macrophages and CD11b + /Ly6G + neutrophils, whereas the CD11c + / MHC class II high DC population was reduced (Fig. 5C). These data suggest that the SCD mice exhibited deregulated homeostasis of the main phagocyte populations in the spleen through hemeinduced expansion of macrophages and neutrophils and depletion of DCs.
To more specifically characterize these changes in the SCD mouse spleen, we performed scRNA-seq experiments with spleen cell suspensions from wild-type and SCD mice that were (1) enriched for macrophages using anti-F4/80 antibody-coated magnetic beads or (2) enriched for DCs using a DC isolation kit based on the depletion of unwanted cells (mainly B cells and erythropoietic cells). Details of the F4/80 + macrophage scRNA-seq analysis of 21,170 cells are provided in Supplementary Figs. 5 and 6. As expected from the results of the flow cytometry experiment (Fig. 5C), red pulp and erythroblastic island macrophages were expanded in the spleen of the SCD mouse (Supplementary Figs. 5 and 6). Consistent with an adaptive macrophage response, we found strongly increased expression of genes related to heme and iron metabolism and enhanced expression of the erythrophagocytic marker MARCO [26] in the hemolytic SCD mouse spleen (Fig. 5D).
In the scRNA-seq analysis of enriched DCs, we used CellRank [38] to assign 9419 wild-type and 10,660 SCD splenic cells to one of 10 macrostates (Supplementary Fig. 7A). Based on a set of marker genes, every macrostate was assigned to one of six hematopoietic cell types ( Fig. 5E and Supplementary Fig. 7B/C). Figure 5F shows UMAP plots separated by genotype with arrows showing the RNA velocity vectors, which represent the potential differentiation trajectories. A significantly lower proportion of conventional DCs (cDCs) was found in the spleen of the SCD mouse than in the wild-type mouse (5.9% vs. 27.3%). Consistent with neutrophilia in the peripheral blood, we found a strongly leftshifted population of neutrophils in the SCD mouse spleen. This is visualized by the scattered distribution of the cells along the RNA velocity vector with more cells positioned close to the vetor origins indicating less mature cells (Fig. 5F). Because categorial cell-type attribution based on marker genes as described above can be ambiguous and prone to bias, we analyzed splenic cells using a second method based on the use of RNA velocity data to assign each cell n probabilities of reaching a predefined number of n-terminal differentiation states (i.e., macrostates) [38,39]. We set the parameters to six terminal differentiation states, which reflected neutrophils, monocytes, plasmacytoid DCs (pDCs), cDCs, B cells, and red pulp macrophages. Supplementary Fig. 7D shows the calculated probability for every cell in the data set to reach any of the six macrostates. Supplementary Fig. 7E shows the cumulative absorbance probabilities split by genotype confirming that the SCD mouse spleen had a contracted cDC compartment.
To further define the subtype of DCs depleted from the SCD mice, we selected all DCs from the scRNA-seq dataset described in Fig. 5 and regrouped them through another round of unsupervised Leiden clustering, which showed five distinct clusters (Fig. 6A). We matched the gene expression profiles of each cluster with subtype-discriminating DC markers that were defined using the population comparison browser from ImmGen [45], which allowed us to attribute clusters 2, 0, and 1 to cDC1s, cDC2s, and pDCs respectively (Fig. 6B). Additionally, we identified cluster 3 as a population of DCs highly similar to our in vitro-generated DCs in GM-CSF BM cultures, which showed notably strong expression of Ccr7 and Cd83 (Fig. 6C). Figure 6C shows gene expression intensity UMAP projections for canonical DC genes that were not used for our Immgen based classification [27,37,[46][47][48][49][50], demonstrating cDC1s with high expression of Xcr1, cDC2s expressing Itgam and Sirpa, and pDCs expressing Siglech. The data indicated that the observed depletion of DCs in the SCD mouse model could be attributed to a deficiency in cDC2s and the strong depletion of Ccr7 + Cd83 + mature DCs. We validated the scRNA-seq data across multiple mice by qRT-PCR analysis of enriched DCs from eight control and eight SCD animals, which revealed profound suppression of the signals for Cd74, H2-Ab1, and the DC maturation markers Cd83, Ccr7, and Cd40 in the SCD mouse spleens compared with the spleens of control littermates. The different genotype resulted in a clear segregation of wild-type and SCD mice in the hierarchical clustering analysis shown in Fig. 6D. Collectively, these phenotype data suggest that SCD mice have severe deficiency of cDC2 and Ccr7 + Cd83 + mature DCs in the spleen.
To determine whether the quantitative DC deficiency in SCD mice observed by flow cytometry and scRNA-seq reflects an impairment of DC function, we administered an agonistic anti-CD40 antibody, which provides a very potent activation signal to DCs [51], and measured the cytokine response in the blood. Consistent with the decreased expression of marker genes for mature CD40 + DCs, we detected significantly attenuated plasma levels of IL12p70 and CXCL10 in SCD mice compared with littermates (Fig. 6E), indicating lower DC activity.
We then aimed to determine whether the depletion of cDCs in the spleens of hemolytic mice coincided with a defective CD4 T cell priming. Because Berkley SCD mice have a mixed genetic background, the validity of functional immunological assays performed with this strain might be limited. Therefore, we used a mouse strain with hereditary spherocytosis that exhibits severe intravascular hemolysis resulting from a mutation in the spectrin alpha (Spta) gene [52]. As shown in Fig. 6F and Supplementary  Fig. 8A, Spta sph/sph mice also demonstrated a reduced number of CD11c + /MHC class II high DCs, and significantly attenuated anti-CD40 antibody-induced plasma levels of IL12p70 and CXCL10 compared to littermate control animals ( Supplementary Fig. 8B), suggesting that DC depletion in the spleen is a generic sequela of hemolysis.
Furthermore, we measured the antigen-presentation capacity of OVA 323-339 -pulsed splenic DCs from Spta sph/sph and wild-type littermates in a coculture assay with CFSE-labeled naive CD4 + The cells were gated from live CD45 + cells. Right: Flow cytometry density plots of the spleen cell suspensions from the mice analyzed at 9 weeks after tamoxifen induction. The suspensions were stained for CD45 and Ly6G. B Hematocrit, reticulocytes, and plasma concentration of heme-adducted albumin in SCD mice and wild-type littermates. C Left: Flow cytometry contour plots of spleen cell suspensions from SCD and wild-type mice stained for neutrophils (CD11b + Ly6G + ), macrophages (F4/80 + ), and DCs (CD11c + MHC class II high ). The cells were run on a spectral flow cytometer and gated as live CD45 + CD3 _ CD19 _ cells. Right: Percentages of neutrophils, macrophages, and DCs in three SCD and control spleens. D Expression heatmap of heme and iron metabolism genes macrophages enriched populations isolated from the spleen of SCD and wild-type mice using anti-F4/80 antibody-coated magnetic dynabeads. The dot size indicates the fraction of cells expressing the genes, and the color indicates the normalized Z-score. E UMAP plot of pooled scRNA-seq data from two sequencing experiments with spleen cell suspensions from SCD and wild-type mice. Before sequencing, the spleen cells were enriched for DCs by negative selection using a DC isolation kit. The colors indicate the marker gene-derived cell-type attribution delineated in Supplementary Fig. 5. F UMAP plot of the scRNA-seq dataset split by genotype and colored by cell type. RNA velocity vectors are superimposed onto the UMAP plots with arrows indicating the velocity direction and magnitude (arrow length). G The stack bar charts illustrate each cell type as a proportion of the whole number of cells per sample. The data in B, C are presented as the means ± SDs. Each dot represents one mouse. t-test (B, C); n.s. not significant, *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001, ****P ≤ 0.0001.
T cells from OT-2 mice. The proliferation of CFSE-labeled lymphocytes was markedly impaired after coculture with DCs enriched from the spleens of Spta sph/sph mice compared with coculture with DCs from the spleens of wild-type mice ( Fig. 6G and Supplementary Fig. 8C). We subsequently evaluated the ability of splenic DCs to initiate a T cell response in vivo. To this end, we transferred CFSE-labeled CD4 + T cells from CD45.1 × OT-2 congenic mice into Spta sph/sph or wild-type littermates and immunized the recipients with OVA protein mixed with an agonistic anti-CD40 antibody as the adjuvant (Fig. 6H). After 3 days, we found that the proliferation of CD45.1 + CD4 + T cells in the spleens of Spta sph/sph mice was significantly decreased compared to that of wild-type animals (Fig. 6I, K).
Overall, the results demonstrate that genetic hemolytic anemia alters the homeostasis of DCs in the spleen, leading to the depletion of cDC2s and mature Ccr7 + Cd83 + DCs. This altered phenotype is accompanied by reductions in cytokine release and antigen-specific CD4 + T cell expansion in the spleen, providing a novel mechanistic underpinning of hyposplenism-related secondary immunodeficiency in hemolytic anemia.

DISCUSSION
Macrophages are pivotal defenders against the consequences of hemolysis [15,20]. This protection relies on the high capacity of specialized macrophages in the spleen and liver to remove damaged RBCs and cell-free heme from the blood and on functionally linked iron-recycling pathways that promote adaptation to anemia [15,16,53,54]. In this study, we explored the hypothesis that hemolysis-enhanced heme stress promotes the generation of homeostatic macrophages at the expense of DCs, which emerge from the same BM progenitors in a potentially competitive way. We found that the activation of NRF2 by heme skewed the differentiation trajectories of BM cultures towards red pulp macrophage-like cells that strongly expressed genes with an antioxidant function and roles in heme and iron metabolism. This physiologically meaningful adaptation-designed to increase disease tolerance during hemolytic stress-deterred the generation of mature DCs and impaired the capacity to prime antigenspecific CD4 T cells. Translational models of genetic hemolytic anemia reproduced the outcome of heme-deregulated myelopoiesis. In models of SCD and spherocytosis, we found that the expansion of heme-and iron-metabolizing macrophages in the spleen was accompanied by a profound deficiency in phenotypically and functionally mature DCs and impaired T cell priming capacity in the spleen.
We obtained the key finding of our study that heme skews myeloid differentiation trajectories in growth-factor supplemented mouse BM cultures. The supplementation of BM cultures with GM-CSF constitutes the first in vitro technique established for the production of DC-like cells, and these cells serve as a key model for basic science and therapeutic innovation in immunology [33]. More recently, the phenotypically and functionally very diverse cellular output of GM-CSF supported BM cultures, consisting of neutrophils, macrophages, and DCs, has called into question the utility of this system for producing mature DCs for research and therapeutic purposes [30,[55][56][57]. However, for our studies, the pronounced cell-type heterogeneity of GM-CSF BM cultures provided a unique opportunity for investigating how heme exposure impacts differentiation trajectories across a broad spectrum of potential cell fates within a single-cell culture dish.
To avoid the potential ambiguities of flow cytometry-based celltype delineation, we decided to utilize a focused scRNA-seq strategy to validate our hypothesis that heme-exposure redirects the differentiation of myeloid precursors towards antioxidant and iron-recycling macrophages at the expense of DC generation. We leveraged the scRNA-seq approach with unbiased cell-type attribution based on RNA velocity analysis, which accounts for the continuity of differentiation states across divergent differentiation trajectories [38,39]. In our GM-CSF BM cell culture model, this high-dimensional analysis confirmed three trajectories for the generation of DCs, macrophages, and neutrophils. Heme exposure fundamentally redirected the differentiation pathway of GM-CSF BM cultures by inhibiting the production of functionally competent DCs and favoring macrophages with a Spic + Hmox1 + Slc40a1 + red pulp macrophage-like phenotype. Interestingly, we found that the heme-suppressible DC population generated in vitro shared prominent phenotypic characteristics, such as high Ccr7 and Cd83 expression, with the DC population that was distinctly depleted in the SCD mouse spleen. Although the gene expression profile of the in vitro DCs did not precisely match the full gene expression profile defined for any of the canonical DC subtypes found in vivo, it was previously noted that DCs produced with this culture method share principal characteristics with Ccr7 high migratory DCs and cDC2s, notably their strong capacity for the antigen-specific priming of naive CD4 + T cells [30,57]. In a broader perspective, our study not only defines a new mechanism by which heme-exposure distorts DCs and immune homeostasis but it also provides a unique high-dimensional analysis of a classical cell culture model, which was and still is very popular amongst immunologists. As such, our data may provide a resource supporting novel interpretations of past and future studies.
We provide substantial evidence showing that the suppression of functionally competent DCs in our BM culture model was mediated by NRF2. NRF2 is the master regulator of the cellular response to oxidative stress. The finding that NRF2 activation by heme or synthetic agonists exerts such a strong effect on cell fate decisions was, therefore, unexpected. Based on prior literature, two principal mechanisms could explain this observation. By supporting the expression of genes involved in heme and iron metabolism, heme-activated NRF2 in myeloid precursors could push probabilistic cell fate decisions toward red pulp macrophage-like phagocytes [26]. In addition, NRF2 suppresses several inflammatory signaling pathways that cooperatively promote the generation of mature DCs, such as the NFkB [58], STAT-1 [59], and STING-related interferon signaling pathways [60]. Through this suppressive activity, heme-activated NRF2 could block the formation of mature DCs at any stage between early progenitor commitment and final maturation.
The intracellular translocation of heme by erythrophagocytosis or the receptor-mediated uptake of heme-protein complexes can suppress proinflammatory transcriptional networks in macrophages [26,59], disrupt actin cytoskeletal dynamics to cause dysfunctional chemotaxis and phagocytosis [61], inhibit Fig. 6 Functional DC depletion in the SCD spleen is caused by a loss of cDC2s and Ccr7 + Cd83 + mature DCs. A Left: All DCs from the scRNA-seq dataset in Fig. 5E/F were extracted and clustered using an unsupervised Leiden algorithm (numbers in parenthesis are cluster identities). For illustration of the differences, the dataset was separated by genotype (wild-type versus SCD). Right: The stack bar charts represent the number of cells per stratified for genotype. Numbers in parenthesis are cell counts. B Cell-type attribution. The dot plot shows the normalized expression intensity (color) and the fraction of positive cells (size) for the top five discriminating ImmGen marker genes for cDC1s, cDC2s, and pDCs in each cluster. The dot size indicates the fraction of cells, the code color the mean expression. The dendrogram shows the hierarchical clustering based on the PCA components between the five clusters. C UMAP intensity projections illustrate the expression of Xcr1, Siglech, Itgam, Sirpa, Ccr7, and Cd83 across the whole dataset of wild-type and SCD DCs. D Hierarchical clustering analysis of the normalized mRNA expression intensities for Cd74, H2-Ab1, Ccr7, Cd40, and Cd83 in isolated splenic DC suspensions. Isolated DC samples were collected from the spleen of eight SCD mice and eight wild-type mice using a DC isolation kit and analyzed by qRT-PCR. Each column indicates one mouse. The almost perfect clustering within genotypes and the long dendrogram distance suggest highly different expressions across the set of marker genes in wild-type versus SCD mice. E Plasma concentrations of CXCL10 and IL12p70 in SCD mice and wild-type littermates treated with saline (control) or an agonistic anti-CD40 antibody measured with a Bio-Plex assay. F Percentage of live CD45 + cells splenic CD11c + MHC class II high DCs in SCD mice, Spta sph/sph mice and their littermates measured by flow cytometry in DC enriched splenic cell suspension. G Enriched DC suspensions from the spleen of Spta sph/sph mice or wild-type littermates were pulsed with OVA 323-339 (10 µg/ml) and cocultured with CFSE-labeled naive CD4 + T cells isolated from OT-2 mice. The proliferation of CD4 + T cells was assessed with flow cytometry by evaluating the degree of CFSE dilution after 4 days in coculture in a 1/3 enriched DCs to T cells ratio. The cells were also stained for the T cell activation marker CD25. Representative data for cells gated based on positive CD4 expression are shown. H Schematic representation of the adoptive T cell transfer experiments. CD4 + T cells were isolated from the spleen of CD45.1 × OT-2 mice, labeled with CFSE and injected i.v. into Spta sph/sph mice or wild-type littermates. Subsequently, the mice were intravenously immunized with ovalbumin protein (75 µg) mixed with an agonistic anti-CD40 antibody. Three days later, the spleens were harvested, CD4 + T cells were negatively enriched and proliferation of CD4 + T cells was assessed by evaluating the degree of CFSE dilution. I Top row: representative density plots of the CD45.1 and CD45. 2 expression in the enriched CD4 fraction isolated from the spleen of immunized mice. Cells were gated based on positive CD4 expression. Bottom row: CFSE dilution in cells gated based on positive CD45.1 and CD4 expression. K Cumulative T cell proliferation data from three independent experiments. The data in F, E, and K are presented as the means ± SDs, and each dot represents one mouse. The data in G are representative of three independent experiments. t-test (F, K); one-way ANOVA with Tukey's multiple comparison test (E); n.s. not significant, *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001, ****P ≤ 0.0001.
proteasome function [62], and enhance immunosuppressive heme metabolites produced by Hmox1 [63]. In light of these prior studies and our own observations, it is likely that heme is a multifactorial driver of hyposplenism-related secondary immunodeficiency and, collectively, these results may stimulate exploring heme-directed therapeutics to reconstitute immune defense in patients with hemolytic anemia.
In summary, our results reveal a role of free heme as a disruptor of immune homeostasis at the interface of innate and adaptive immunity and delineate a novel mechanism underlying hyposplenism.
Mice of both sexes that were 12-16 weeks old were used for all experiments. All experimental protocols were reviewed and approved by the Veterinary Office of the Canton of Zurich (ZH161 2020). All animals were maintained at the animal facility of the University of Zurich (LASC) and treated in accordance with guidelines provided by the Swiss Federal Veterinary Office. For comparisons of genotypes, we always selected matched wild-type and knockout mice from the same litter. For CD40 treatment experiments matched animals were allocated to control and treatment groups randomly. Sample analysis was performed by a person blinded to the treatment group. Mice have a numerical identifier provided by the animal facility.

BM cell cultures
BM cells were isolated by flushing the femurs and tibias of 8-to 10-weekold mice, followed by straining the BM through a 70-μm filter. Lin − progenitor cells were isolated using a Lin + cell depletion kit (Lin: CD5, CD11b, CD45R (B220), Gr-1 (Ly-6G/C), 7-4, and Ter-119; Miltenyi Biotec) and 4 × 10 5 cells pro well were plated in tissue culture-treated six-well plates (Nunc Multidishes with UpCell Surface, Thermo Fisher) in 4 ml of complete RPMI-1640 medium (10% fetal calf serum (FCS), 1% L-glutamine) supplemented with 1% penicillin/streptomycin (P/S) and 20 ng/ml recombinant mouse GM-CSF or 100 ng/ml recombinant mouse M-CSF (PeproTech). For the M-CSF-supplemented cultures, half of the medium was removed on day 3 and new medium supplemented with M-CSF (100 ng/ml) was added. For the GM-CSF-supplemented cultures, fresh medium containing GM-CSF (2×, 40 ng/ml) was added on day 2. On day 3, half of the medium was removed and new medium supplemented with GM-CSF (20 ng/ml) was added. Hmox1 and Keap1 knockout were obtained by treating the cells with 4-hydroxytamoxifen (1 μg/ml) (Sigma). Cultures were treated on day 3 as indicated in the results section with heme-albumin, albumin (as vehicle/control), 10 μM ML-314 (Tocris), or 15 μM RA-839 (Tocris). The BM cells were harvested for analysis on day 7 (unless mentioned otherwise) from the thermosensitive cell culture plates after cooling to room temperature for 30 min.

Porphyrin preparation for cell culture
Hemin (heme-chloride) and SnMP (tin mesoporphyrin) were obtained from Frontier Scientific (Newark). The porphyrins were dissolved in 10 ml of NaOH (100 mM) at 37°C, and 10 ml of 20% human serum albumin (CSL Behring AG) was then added. After 1 h of incubation at 37°C, the pH of the solution was adjusted to pH 7.4 using ortho-phosphoric acid, and the final volume was adjusted to 25 ml with a saline solution. The porphyrinalbumin solutions were filter-sterilized (0.22 µm) and used immediately [67].

Tamoxifen treatment of mice
Tamoxifen was prepared by dissolving in corn oil for a final concentration of 80 mg/ml. Hmox1 tm1.1Hes and littermates were given 60 μl (4.8 mg) tamoxifen solution once a day for 5 days by oral gavage.

Preparation of spleen cell suspensions and cell enrichment
Spleens were harvested and mechanically disrupted in PBS and passed through a 70-μm cell strainer. The cell suspensions were then centrifuged, incubated in RBC lysis buffer (BioLegend) for 2 min at 37°C and centrifuged once more to obtain spleen cell populations devoid of mature erythrocytes. DCs were negatively enriched from spleen single-cell suspensions using a Dynabeads DC enrichment kit (Invitrogen) according to the manufacturer's instructions. Macrophages were positively enriched from spleen single-cell suspensions using anti-F4/80 antibody (rat, Bioscience, Clone T45-2342)-coated Dynabeads (anti-rat IgG, Invitrogen) according to the manufacturer's instructions.
CD4 + T and CD8 + T cell isolation and CFSE labeling T cells were negatively enriched from spleen single-cell suspensions using a MagniSort CD4 or CD8 enrichment kit (Invitrogen) according to the manufacturer's instructions. Isolated CD4 + T or CD8 + T cells were labeled with CFSE (Thermo Fisher) at 37°C for 20 min, washed with PBS and counted before use. The final purity confirmed by flow cytometry to be >95%.

OT-2 assay
A total of 10 4 M-CSF-or GM-CSF BM cells or DCs enriched from spleen cell suspensions were plated in 96-well round-bottom plates, pulsed or not with 1 or 10 μg/ml OVA peptide 323-339 (Sigma) for 30 min at 37°C, and washed three times with PBS. Subsequently, 5 × 10 4 CFSE-labeled naive CD4 + T cells isolated from spleens of OT-2 mice were added to the BM cells or enriched splenic DCs in complete RPMI-1640 medium and cocultured at 37°C. CFSE dilution was assessed by flow cytometry and IL-2 concentration was measured in the supernatant by Bio-Plex assay after 3 or 4 days as specified.

OT-1 assay
A total of 10 4 GM-CSF BM cells were plated in 96-well round-bottom plates, pulsed or not with 150 μg/ml ovalbumin (InvivoGen) for 30 min at 37°C, and washed three times with PBS. Subsequently, 5 × 10 4 CFSE-labeled naive CD8 + T cells isolated from spleens of OT-1 mice were added to the BM cells in complete RPMI-1640 medium and cocultured at 37°C. CFSE dilution was assessed by flow cytometry.
Transfer of GM-CSF phagocytes and OT-2-specific CD4 + T cells into CD45.1 Rag1 −/− mice GM-CSF BM cells that were exposed to heme-albumin or albumin (vehicle/ control) as described in the BM-derived cell cultures section were F. Vallelian et al. incubated in complete RPMI-1640 medium with 10 μg/ml OVA 323-339 for 30 min at 37°C. The cells were washed three times in PBS and injected intravenously via the tail vein into Rag1 −/− mice expressing the CD45.1 antigen 2 h after the intravenous transfer of 1.5 × 10 6 CFSE-labeled CD4 + T cells isolated from OT-2 mouse spleens. Three days later, the CFSE dilution of splenic CD45.2 + CD4 + T cells was evaluated by flow cytometry.
Transfer of OT-2 CD4 + T cells into Spta sph/sph and Spta wt/wt mice First, 1.5 × 10 6 CFSE-labeled CD4 + T cells isolated from OT-2 × CD45.1 mouse spleens were injected intravenously via the tail vein into Spta sph/sph mice and wild-type littermates (F1 C57BL/6J × WB). Two hours later, the mice were challenged intravenously with 75 μg of OVA protein (InvivoGen) mixed with 75 μg of agonistic anti-CD40 antibodies. The mice were sacrificed 3 days later, CD4 + T cells were enriched from the spleen by negative selection and the CFSE dilution of CD45.1 + CD4 + T cells was assessed by flow cytometry.

Bio-Plex cytokine assays
The concentrations of IL-2, IL12p70, and CXCL10 in coculture supernatants or plasma were determined using Bio-Plex Cytokine Assays (Bio-Rad). The assays were performed with a Bio-Plex 200 system (Bio-Rad), and the results were analyzed using Bio-Plex Data Pro software (Bio-Rad).

Flow cytometry
Cells were preincubated with a LIVE/DEAD Fixable Near-IR cell stain kit (Invitrogen) and with Mouse BD Fc Block™ (≤1 μg/million cells in 100 μl, BD Biosciences) at 4°C for 10 min.

Western blot analysis and nuclei isolation procedure
Total cellular and nuclear proteins were extracted from cultured human HEK293 cells grown in 10 cm cell culture dishes under standard growth conditions. For total cellular protein analysis, cells were lysed in in RIPA buffer containing 50 mM Tris-HCl pH 7.6, 150 mM NaCl, 1% Nonidet P-40 substitute, 0.5% sodium deoxycholate, 0.1% SDS and 0.2 U/μl RNase inhibitor. Nuclear proteins were analysed after a nuclei isolation procedure according to a modified 10× Genomics protocol. In brief, cells were lysed on ice for 5 min in chilled (4°C) lysis buffer containing 10 mM Tris-HCl, 10 mM NaCl, 3 mM MgCl 2 , and 0.1% Nonidet P-40 substitute in nuclease-free water. Lysis efficiency was assessed by microscopy and proved to be over 80%. Nuclei were then washed twice with chilled (4°C) nuclei wash buffer containing 1% BSA and 0.2 U/μl RNase inhibitor in phosphate-buffered saline, and finally resuspended in RIPA buffer. After three freeze-thaw cycles, cellular debris was removed by centrifugation at 16,000 g for 15 min. Protein concentrations of all samples were determined using a Pierce BCA protein assay (Thermo Fisher Scientific, Rockford, IL, USA). A reducing SDS-PAGE was loaded with 3 μg total protein/well for subsequent Western blot analysis using a Criterion TGX Stain-Free Precast Gel (Bio-Rad Laboratories, Hercules, CA USA). Fluorescent gel imaging for total protein analysis was performed according to the manufacturer's protocol with a ChemiDoc MP imager (Bio-Rad Laboratories). Proteins were transferred to a polyvinylidene difluoride membrane (Bio-Rad Laboratories), which was blocked by incubating for 1 h at room temperature in 10% goat serum (Thermo Fisher Scientific) and 1% BSA in PBS on an orbital shaker. For immunodetection, a polyclonal primary antibody directed against NRF2 (rabbit anti-human, Thermo Fisher Scientific; 1:2000) and an appropriate HRP-conjugated secondary antibody (goat anti-rabbit, Bio-Rad Laboratories, 1:3000) were used. The blot was developed with Clarity Max Western ECL Substrate (Bio-Rad Laboratories) and analyzed on the ChemiDoc MP imager with Image Lab Analysis software, Version 5.2.1 (Bio-Rad Laboratories). Normalized band intensities were calculated by division of the background-subtracted NRF2 band volumes by the corresponding background-subtracted total protein lane volumes. For all experiments, the sample volume was adjusted to a target capture of 10,000 cells and loaded on the 10× Genomics chromium next-GEM chip G to generate gel-beads-in-emulsion (GEMs). The GEM solution was placed in a Applied Biosystems Veriti 96 well thermal cycler for reverse transcription as described by the 10× Genomics instruction guide (53:00 min at 53°C followed by 5:00 min at 85°C). The resulting barcoded cDNA was then cleaned using Dynabeads MyOne Silane and amplified for 11 cycles (as recommended by the 10× Genomics user guide for a target cell recovery of >6000 cells). After amplification, for multiplexed experiments, cDNA generated from polyadenylated mRNA for the 3' gene expression library was separated from DNA from the Cell Surface Protein Feature Barcode for the Cell Surface Protein library with Dynabeads MyOne Silane and SPRIselect reagents based on size. The quality and concentration of both cDNA and DNA were assessed using High-Sensitivity D5000 ScreenTape (Agilent). All samples presented product sizes with a narrow distribution centered around 2000 pb and yielded between 50 and 800 ng cDNA (manually selecting products between 100-250 and 5000-6000 bp). cDNA and DNA were then subjected to enzymatic fragmentation, end repair and A-tailing. Adaptors were ligated to the fragmented cDNA and DNA, and the sample index was added during sample index PCR (set for 12 cycles, as recommended by the 10× Genomics user guide to correlate with a cDNA/DNA input of 12-150 ng). Library quality and concentration were assessed using High-Sensitivity D5000 ScreenTape (Agilent). All libraries showed an average fragment size of around 400 pb. For multiplexed runs, 3ʹ Gene Expression and Cell Surface Protein libraries were pooled at a ratio of 4:1 and sequenced using the Illumina NovaSeq 6000 system with a sequencing depth of 50,000 and 12,500 reads per cell, respectively, following the recommendations of 10× Genomics (paired-end reads, single indexing, read 1 = 28 cycles, i7 = 8 cycles, i5 = 0 cycles and read 2 = 91 cycles). For non multiplexed runs, 3ʹ Gene Expression libraries for each sample were pooled at an equimolar amount and sequenced using the Illumina NovaSeq 6000 system with a sequencing depth of 50,000 reads per cell, following the same recommendations of 10× Genomics.

Sequencing-based workflows and data analysis
scRNA-seq data analysis. The Cell Ranger Single-Cell Software Suite (version 4.0.0) was used for cDNA oligopeptide alignment, barcode assignment and UMI counting from fastq data from the Illumina sequencing. For each sample, the cell-containing droplets were filtered from the empty droplets, and this step was followed by the generation of an expression matrix using Cell Ranger Count (version 4.0.0). Demultiplexing of the cells within each sample was performed with the filtered matrix produced by Cell Ranger in R (version 4.4) using Seurat (version 3.2.3) and the HTODemux function (positive quantile set at 0.99). The resulting gene expression matrices were further analyzed with Python (version 3.8.8) using the Scanpy (version 1.6.0) [68], scVelo (version 0.2.2) [39] and CellRank (version 1.1.0) [38] libraries. Cells with a total number of expressed genes >5000 or <500, a proportion of mitochondrial genes <15%, or a proportion of ribosomal genes <30% and genes expressed in <20 cells were excluded from downstream analyses. After filtering, normalization was performed using the scran normalization algorithm implemented in the scran R package (version 1.18.3), which was followed by log normalization. For dimensional reduction, principal component analysis (PCA) of each sample was performed using the tl.pca function with the default settings from the Scanpy library. For cell culture samples, a shared nearest neighbor graph was built using the pp.neighbors function based on the first 15 PCAs or first 30 PCAs, respectively. Using the Leiden algorithm [69], cells were clustered and visualized in 2D using UMAP. For expression intensity projections, the cells were colored according to their log-normalized gene expression. Differentially expressed genes (p-value <0.05 and a log2-fold change >0.5) were determined using the tl. rank_genes_groups function from the Scanpy library with a pairwise Wilcoxon rank-sum test. For RNA velocity analysis, an expression matrix with spliced and unspliced transcript reads was calculated using the velocyto pipeline [70] by applying the Python script velocyto.py (version 0.17.17) to the Cell Ranger output folder. Using the scVelo library, RNA velocities were calculated by stochastically modeling the transcriptional dynamics of splicing kinetics, and the velocities were then projected as streamlines onto the low-dimensional UMAP embedding. With CellRank, cells were assigned to a predefined number of macrostates combining RNA velocity and transcriptional similarity information.
Pathway and transcription factor enrichment analysis of scRNA-seq data. GSEA was performed using the GSEApy library [71]. As the input, we used a gene list consisting of differentially expressed genes between the different conditions. Using gene lists ranked according to a combined score calculated with the tl.rank_genes_groups function from Scanpy, we performed pre ranked GSEA using the hallmark gene sets from the Molecular Signature Database (MSigDB) [72]. For visualization, the terms were ordered according to their normalized enrichment score, and the running enrichment score for each term was plotted along the ranked genes as a heatmap.
Nrf2 target gene expression scores were calculated for each cell using the average expression of the Nrf2 gene set from the TRRUST gene-set library subtracted with the average expression of a randomly sampled reference set of similar length.
scATAC-seq data acquisition. One million cells per condition were treated for nuclei isolation following 10× genomics recommendation. Briefly, the cells were incubated in lysis buffer (Tris-HCL 10 mM, NaCl 10 mM, MgCl2 3 mM, BSA 1%, Tween-20 0.1%, Nonidet P40 substitute 0.1% and Digitonin 0.01%) for 3 min and washed in washing buffer (Tris-HCL 10 mM, NaCl 10 mM, MgCl2 3 mM, BSA 1%, Tween-20 0.1%) before resuspension in PBS supplemented with BSA. Nuclei quality was assessed using trypan blue staining (brightfield microscopy) and hoechst staining (fluorescence microscopy). Nuclei were then processed according to the chromium next GEM single-cell ATAC reagents kit v1.1 (10× genomics). Sample volume was adjusted for a target capture of 10,000 nuclei and incubated with Tn5 transposase for 60 min at 37°C. Transposase-treated nuclei were then loaded on the chromium H chip to generate GEMs. The GEM solution was placed in a Applied Biosystems Veriti 96 well thermal cycler for 12 cycles to generate 10× barcoded single strand DNA. DNA solutions were cleaned from leftover reagents using silane magnetic beads and solid phase reversible immobilization (SPRI) beads were used to eliminate unused barcodes. Sample index PCR was performed with 9 cycles as recommended by 10× genomics for an input of 6000-10,000 nuclei and the products were cleaned using SPRI beads. Library quality and concentration were assessed using High-Sensitivity D5000 ScreenTape (Agilent). All libraries showed fragment size between 150 and 700 pb. Libraries were sequenced using the Illumina NovaSeq 6000 system with a sequencing depth of 25,000 reads per nuclei, following the recommendations of 10× Genomics (paired-end reads, dual indexing, read 1 = 50 cycles, i7 = 8 cycles, i5 = 16 cycles and read 2 = 50 cycles).
scATAC-seq data analysis. Demultiplexing, read filtering, barcode counting and alignment on the Grcm38p6 (mm10) reference genome was performed using CellRangerATAC (version2.0.0) and the two samples aggregated using cellrangerATAC aggr command with correction for sequencing depth. Data was then analyzed using R (version 4.1.2) and the packages Seurat (version 4.0.5) [73], Signac (version 1.4.0) [74], JASPAR2020 (version 0.99.10) [75] and chromVAR (version 1.16.0) [76]. Peak calling from cellranger output was used for quality control metrics. Quality control metrics were obtained from the output of Cellranger (percentage read in peaks and ratio of fragments in blacklist regions) and calculated using the NucleosomeSignal and TSSenrichement functions. Nuclei with >3000 and <100,000 peak fragments, >40% reads in peaks, <0.025% peaks in blacklist region, <4 nucleosome signal score and >2 TSS enrichment score were kept. Peak calling was performed using Macs2 [77] and the output was used for further analysis. Normalization, highly variable feature finding and dimension reduction were carried out using the RunTFIDF, FindTopfeatures and RunSVD functions, respectively. RunUMAP, Findneighbors were then run on the latent semantic indexing dimensions 2-30 to avoid accounting technical difference due to sequencing depth. Differential accessible regions were calculated using the FindMarkers function using a logistic regression framework and a likelihood ratio test. Genes closest to differentially accessible peaks were determined using the ClosestFeature function. Motif analysis on differentially accessible peaks was performed based on the JASPAR2020 [75] database using the FindMotifs function. Per-cell motif activity score was computed using the ChromVAR function. A pseudo-gene activity matrix was calculated using the Geneactivity function which sums the number of fragments overlapping the promoter and gene body region for each cell (limited to the region 2 kbp upstream and downstream from gene coordinates). Cell-type attribution was based on this gene activity matrix.
RT-qPCR analysis. Total RNA was isolated from spleen-cell suspensions using the RNeasy Mini Kit (Qiagen) according to the manufacturer's instructions. Reverse transcription was performed with TaqMan reverse transcription reagents (Life Technologies). Real-time PCR was performed using Fast SYBR™ Green Master Mix (Applied Biosystems) to determine the expression levels of target genes using the primers listed in Table 1 shown below. Relative mRNA levels for experimental samples were calculated with 7500 Fast System Sequence Detection Software version 1.4 (Applied Biosystems) after normalization to the Hprt levels.

Free heme measurement
Mouse plasma samples were analyzed by SEC-high-performance liquid chromatography using an Agilent 1260 II HPLC attached to a quaternary pump and a photodiode array detector (DAD) (Agilent). Plasma samples and hemoglobin standards were separated on a Diol-300 (3 µm, 300 × 8.0 mm) column (YMC Co., Ltd.) with PBS (pH 7.4, Bichsel) as the mobile phase at a flow rate of 1 mL/min. For all samples, two wavelengths were recorded (λ = 280 nm and λ = 414 nm). From these data, we quantified peak areas of hemoglobin-haptoglobin complexes, free hemoglobin (dimers and tetramers), and other heme-adducted proteins [78].

Blood analysis
Hematocrit levels of blood samples of SCD and wild-type mice were measured by the Veterinary Laboratory of the University of Zurich.

Statistical analysis
Data plotting and statistical analysis were performed with Prism 9 (GraphPad), JMP 15 (SAS). For intergroup comparisons, we used a t-test or ANOVA with Tukey's posttest as indicated in the figure legends. All data points are displayed in the graphs as mean ± standard deviation (n.s. = not significant, *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001, ****P ≤ 0.0001). Normal distribution with similar variance between groups was assumed. No samples were excluded except samples were physically lost during the experiment (e.g., failure of cardiac puncture and blood sampling in a mouse).