Mathematical modeling of erythrocyte chimerism informs genetic intervention strategies for sickle cell disease

In vivo measurement of RBC survival in patients with sickle cell disease before or after hematopoietic stem cell transplantation

ABSTRACT

Stable, mixed-donor-recipient chimerism after allogeneic hematopoietic stem cell transplantation (HSCT) for patients with sickle cell disease (SCD) is sufficient for phenotypic disease reversal, and results from differences in donor/recipient-red blood cell (RBC) survival. Understanding variability and predictors of RBC survival among patients with SCD before and after HSCT is critical for gene therapy research which seeks to generate sufficient corrected hemoglobin to reduce polymerization thereby overcoming the red cell pathology of SCD. This study used biotin labeling of RBCs to determine the lifespan of RBCs in patients with SCD compared with patients who have successfully undergone curative HSCT, participants with sickle cell trait (HbAS), and healthy (HbAA) donors. Twenty participants were included in the analysis (SCD pre-HSCT: N = 6, SCD post-HSCT: N = 5, HbAS: N = 6, and HbAA: N = 3). The average RBC lifespan was significantly shorter for participants with SCD pre-HSCT (64.1 days; range, 35-91) compared with those with SCD post-HSCT (113.4 days; range, 105-119), HbAS (126.0 days; range, 119-147), and HbAA (123.7 days; range, 91-147) (P<.001). RBC lifespan correlated with various hematologic parameters and strongly correlated with the average final fraction of sickled RBCs after deoxygenation (P<.001). No adverse events were attributable to the use of biotin and related procedures. Biotin labeling of RBCs is a safe and feasible methodology to evaluate RBC survival in patients with SCD before and after HSCT. Understanding differences in RBC survival may ultimately guide gene therapy protocols to determine hemoglobin composition required to reverse the SCD phenotype as it relates directly to RBC survival. This trial was registered at www.clinicaltrials.gov as #NCT04476277.

A systems pharmacology model for gene therapy in sickle cell disease

We developed a mathematical model for autologous stem cell therapy to cure sickle cell disease (SCD). Experimental therapies using this approach seek to engraft stem cells containing a curative gene. These stem cells are expected to produce a lifelong supply of red blood cells (RBCs) containing an anti‐sickling hemoglobin. This complex, multistep treatment is expensive, and there is limited patient data available from early clinical trials. Our objective was to quantify the impact of treatment parameters, such as initial stem cell dose, efficiency of lentiviral transduction, and degree of bone marrow preconditioning on engraftment efficiency, peripheral RBC numbers, and anti‐sickling hemoglobin levels over time. We used ordinary differential equations to model RBC production from progenitor cells in the bone marrow, and hemoglobin assembly from its constituent globin monomers. The model recapitulates observed RBC and hemoglobin levels in healthy and SCD phenotypes. Treatment simulations predict dynamics of stem cell engraftment and RBC containing the therapeutic gene product. Post‐treatment dynamics show an early phase of reconstitution due to short lived stem cells, followed by a sustained RBC production from stable engraftment of long‐term stem cells. This biphasic behavior was previously reported in the literature. Sensitivity analysis of the model quantified relationships between treatment parameters and efficacy. The initial dose of transduced stem cells, and the intensity of myeloablative bone marrow preconditioning are predicted to most positively impact long‐term outcomes. The quantitative systems pharmacology approach used here demonstrates the value of model‐assisted therapeutic design for gene therapies in SCD.

Genetic therapies for the first molecular disease

Sickle cell disease (SCD) is a monogenic disorder characterized by recurrent episodes of severe bone pain, multi-organ failure, and early mortality. Although medical progress over the past several decades has improved clinical outcomes and offered cures for many affected individuals living in high-income countries, most SCD patients still experience substantial morbidity and premature death. Emerging technologies to manipulate somatic cell genomes and insights into the mechanisms of developmental globin gene regulation are generating potentially transformative approaches to cure SCD by autologous hematopoietic stem cell (HSC) transplantation. Key components of current approaches include ethical informed consent, isolation of patient HSCs, in vitro genetic modification of HSCs to correct the SCD mutation or circumvent its damaging effects, and reinfusion of the modified HSCs following myelotoxic bone marrow conditioning. Successful integration of these components into effective therapies requires interdisciplinary collaborations between laboratory researchers, clinical caregivers, and patients. Here we summarize current knowledge and research challenges for each key component, emphasizing that the best approaches have yet to be developed.

Post-Transcriptional Genetic Silencing of BCL11A to Treat Sickle Cell Disease

BACKGROUND

Sickle cell disease is characterized by hemolytic anemia, pain, and progressive organ damage. A high level of erythrocyte fetal hemoglobin (HbF) comprising α- and γ-globins may ameliorate these manifestations by mitigating sickle hemoglobin polymerization and erythrocyte sickling. BCL11A is a repressor of γ-globin expression and HbF production in adult erythrocytes. Its down-regulation is a promising therapeutic strategy for induction of HbF.

METHODS

We enrolled patients with sickle cell disease in a single-center, open-label pilot study. The investigational therapy involved infusion of autologous CD34+ cells transduced with the BCH-BB694 lentiviral vector, which encodes a short hairpin RNA (shRNA) targeting BCL11A mRNA embedded in a microRNA (shmiR), allowing erythroid lineage-specific knockdown. Patients were assessed for primary end points of engraftment and safety and for hematologic and clinical responses to treatment.

RESULTS

As of October 2020, six patients had been followed for at least 6 months after receiving BCH-BB694 gene therapy; median follow-up was 18 months (range, 7 to 29). All patients had engraftment, and adverse events were consistent with effects of the preparative chemotherapy. All the patients who could be fully evaluated achieved robust and stable HbF induction (percentage HbF/(F+S) at most recent follow-up, 20.4 to 41.3%), with HbF broadly distributed in red cells (F-cells 58.9 to 93.6% of untransfused red cells) and HbF per F-cell of 9.0 to 18.6 pg per cell. Clinical manifestations of sickle cell disease were reduced or absent during the follow-up period.

CONCLUSIONS

This study validates BCL11A inhibition as an effective target for HbF induction and provides preliminary evidence that shmiR-based gene knockdown offers a favorable risk-benefit profile in sickle cell disease. (Funded by the National Institutes of Health; ClinicalTrials.gov number, NCT03282656).

Preclinical Evaluation of a Novel Lentiviral Vector Driving Lineage-Specific BCL11A Knockdown for Sickle Cell Gene Therapy

In this work we provide preclinical data to support initiation of a first-in-human trial for sickle cell disease (SCD) using an approach that relies on reversal of the developmental fetal-to-adult hemoglobin switch. Erythroid-specific knockdown of BCL11A via a lentiviral-encoded microRNA-adapted short hairpin RNA (shRNAmiR) leads to reactivation of the gamma-globin gene while simultaneously reducing expression of the pathogenic adult sickle β-globin. We generated a refined lentiviral vector (LVV) BCH-BB694 that was developed to overcome poor vector titers observed in the manufacturing scale-up of the original research-grade LVV. Healthy or sickle cell donor CD34+ cells transduced with Good Manufacturing Practices (GMP)-grade BCH-BB694 LVV achieved high vector copy numbers (VCNs) >5 and gene marking of >80%, resulting in a 3- to 5-fold induction of fetal hemoglobin (HbF) compared with mock-transduced cells without affecting growth, differentiation, and engraftment of gene-modified cells in vitro or in vivo. In vitro immortalization assays, which are designed to measure vector-mediated genotoxicity, showed no increased immortalization compared with mock-transduced cells. Together these data demonstrate that BCH-BB694 LVV is non-toxic and efficacious in preclinical studies, and can be generated at a clinically relevant scale in a GMP setting at high titer to support clinical testing for the treatment of SCD.

Unraveling Disease Pathophysiology with Mathematical Modeling

Modeling has enabled fundamental advances in our understanding of the mechanisms of health and disease for centuries, since at least the time of William Harvey almost 500 years ago. Recent technological advances in molecular methods, computation, and imaging generate optimism that mathematical modeling will enable the biomedical research community to accelerate its efforts in unraveling the molecular, cellular, tissue-, and organ-level processes that maintain health, predispose to disease, and determine response to treatment. In this review, we discuss some of the roles of mathematical modeling in the study of human physiology and pathophysiology and some challenges and opportunities in general and in two specific areas: in vivo modeling of pulmonary function and in vitro modeling of blood cell populations.

Editing a γ-globin repressor binding site restores fetal hemoglobin synthesis and corrects the sickle cell disease phenotype

Sickle cell disease (SCD) is caused by a single amino acid change in the adult hemoglobin (Hb) β chain that causes Hb polymerization and red blood cell (RBC) sickling. The co-inheritance of mutations causing fetal γ-globin production in adult life hereditary persistence of fetal Hb (HPFH) reduces the clinical severity of SCD. HPFH mutations in the HBG γ-globin promoters disrupt binding sites for the repressors BCL11A and LRF. We used CRISPR-Cas9 to mimic HPFH mutations in the HBG promoters by generating insertions and deletions, leading to disruption of known and putative repressor binding sites. Editing of the LRF-binding site in patient-derived hematopoietic stem/progenitor cells (HSPCs) resulted in γ-globin derepression and correction of the sickling phenotype. Xenotransplantation of HSPCs treated with gRNAs targeting the LRF-binding site showed a high editing efficiency in repopulating HSPCs. This study identifies the LRF-binding site as a potent target for genome-editing treatment of SCD.

Genome editing of HBG1 and HBG2 to induce fetal hemoglobin

Induction of fetal hemoglobin (HbF) via clustered regularly interspaced short palindromic repeats/Cas9-mediated disruption of DNA regulatory elements that repress γ-globin gene (HBG1 and HBG2) expression is a promising therapeutic strategy for sickle cell disease (SCD) and β-thalassemia, although the optimal technical approaches and limiting toxicities are not yet fully defined. We disrupted an HBG1/HBG2 gene promoter motif that is bound by the transcriptional repressor BCL11A. Electroporation of Cas9 single guide RNA ribonucleoprotein complex into normal and SCD donor CD34+ hematopoietic stem and progenitor cells resulted in high frequencies of on-target mutations and the induction of HbF to potentially therapeutic levels in erythroid progeny generated in vitro and in vivo after transplantation of hematopoietic stem and progenitor cells into nonobese diabetic/severe combined immunodeficiency/Il2rγ-/-/KitW41/W41 immunodeficient mice. On-target editing did not impair CD34+ cell regeneration or differentiation into erythroid, T, B, or myeloid cell lineages at 16 to 17 weeks after xenotransplantation. No off-target mutations were detected by targeted sequencing of candidate sites identified by circularization for in vitro reporting of cleavage effects by sequencing (CIRCLE-seq), an in vitro genome-scale method for detecting Cas9 activity. Engineered Cas9 containing 3 nuclear localization sequences edited human hematopoietic stem and progenitor cells more efficiently and consistently than conventional Cas9 with 2 nuclear localization sequences. Our studies provide novel and essential preclinical evidence supporting the safety, feasibility, and efficacy of a mechanism-based approach to induce HbF for treating hemoglobinopathies.

Replicative cellular age distributions in compartmentalized tissues

The cellular age distribution of hierarchically organized tissues can reveal important insights into the dynamics of cell differentiation and self-renewal and associated cancer risks. Here, we examine the effect of progenitor compartments with varying differentiation and self-renewal capacities on the resulting observable distributions of replicative cellular ages. We find that strongly amplifying progenitor compartments, i.e. compartments with high self-renewal capacities, substantially broaden the age distributions which become skewed towards younger cells with a long tail of few old cells. For several of these strongly amplifying compartments, the age distribution becomes virtually independent of the influx from the stem cell compartment. By contrast, if tissues are organized into many downstream compartments with low self-renewal capacity, the shape of the replicative cell distribution in more differentiated compartments is dominated by stem cell dynamics with little added variation. In the limiting case of a strict binary differentiation tree without self-renewal, the shape of the output distribution becomes indistinguishable from that of the input distribution. Our results suggest that a comparison of cellular age distributions between healthy and cancerous tissues may inform about dynamical changes within the hierarchical tissue structure, i.e. an acquired increased self-renewal capacity in certain tumours. Furthermore, we compare our theoretical results to telomere length distributions in granulocyte populations of 10 healthy individuals across different ages, highlighting that our theoretical expectations agree with experimental observations.

Extensive multilineage analysis in patients with mixed chimerism after allogeneic transplantation for sickle cell disease: insight into hematopoiesis and engraftment thresholds for gene therapy

Although studies of mixed chimerism following hematopoietic stem cell transplantation in patients with sickle cell disease (SCD) may provide insights into the engraftment needed to correct the disease and into immunological reconstitution, an extensive multilineage analysis is lacking. We analyzed chimerism simultaneously in peripheral erythroid and granulomonocytic precursors/progenitors, highly purified B and T lymphocytes, monocytes, granulocytes and red blood cells (RBC). Thirty-four patients with mixed chimerism and ≥12 months of follow-up were included. A selective advantage of donor RBC and their progenitors/precursors led to full chimerism in mature RBC (despite partial engraftment of other lineages), and resulted in the clinical control of the disease. Six patients with donor chimerism <50% had hemolysis (reticulocytosis) and higher HbS than their donor. Four of them had donor chimerism <30%, including a patient with AA donor (hemoglobin >10 g/dL) and three with AS donors (hemoglobin <10 g/dL). However, only one vaso-occlusive crisis occurred with 68.7% HbS. Except in the patients with the lowest chimerism, the donor engraftment was lower for T cells than for the other lineages. In a context of mixed chimerism after hematopoietic stem cell transplantation for SCD, myeloid (rather than T cell) engraftment was the key efficacy criterion. Results show that myeloid chimerism as low as 30% was sufficient to prevent a vaso-occlusive crisis in transplants from an AA donor but not constantly from an AS donor. However, the correction of hemolysis requires higher donor chimerism levels (i.e ≥50%) in both AA and AS recipients. In the future, this group of patients may need a different therapeutic approach.

Insights into Gene Therapy for Urea Cycle Defects by Mathematical Modeling

Metabolic liver diseases are attractive gene therapy targets that necessitate reconstitution of enzymatic activity in functionally complex biochemical pathways. The levels of enzyme activity required in individual hepatocytes and the proportion of the hepatic cell mass that must be gene corrected for therapeutic benefit vary in a disease-dependent manner that is difficult to predict. While empirical evaluation is inevitably required, useful insights can nevertheless be gained from knowledge of disease pathophysiology and theoretical approaches such as mathematical modeling. Urea cycle defects provide an excellent example. Building on a previously described one-compartment model of the urea cycle, we have constructed a two-compartment model that can simulate liver-targeted gene therapy interventions using the computational program Mathematica. The model predicts that therapeutically effective reconstitution of ureagenesis will correlate most strongly with the proportion of the hepatic cell mass transduced rather than the level of enzyme-encoding transgene expression achieved in individual hepatocytes. Importantly, these predictions are supported by experimental data in mice and human genotype/phenotype correlations. The most notable example of the latter is ornithine transcarbamylase deficiency (X-linked) where impairment of ureagenesis in male and female patients is closely simulated by the one- and two-compartment models, respectively. Collectively, these observations support the practical value of mathematical modeling in evaluation of the disease-specific gene transfer challenges posed by complex metabolic phenotypes.

New insights into RAS biology reinvigorate interest in mathematical modeling of RAS signaling

RAS is the most frequently mutated gene across human cancers, but developing inhibitors of mutant RAS has proven to be challenging. Given the difficulties of targeting RAS directly, drugs that impact the other components of pathways where mutant RAS operates may potentially be effective. However, the system-level features, including different localizations of RAS isoforms, competition between downstream effectors, and interlocking feedback and feed-forward loops, must be understood to fully grasp the opportunities and limitations of inhibiting specific targets. Mathematical modeling can help us discern the system-level impacts of these features in normal and cancer cells. New technologies enable the acquisition of experimental data that will facilitate development of realistic models of oncogenic RAS behavior. In light of the wealth of empirical data accumulated over decades of study and the advancement of experimental methods for gathering new data, modelers now have the opportunity to advance progress toward realization of targeted treatment for mutant RAS-driven cancers.

Emerging Genetic Therapy for Sickle Cell Disease

The genetic basis of sickle cell disease (SCD) was elucidated >60 years ago, yet current therapy does not rely on this knowledge. Recent advances raise prospects for improved, and perhaps curative, treatment. First, transcription factors, BCL11A and LRF/ZBTB7A, that mediate silencing of the β-like fetal (γ-) globin gene after birth have been identified and demonstrated to act at the γ-globin promoters, precisely at recognition sequences disrupted in rare individuals with hereditary persistence of fetal hemoglobin. Second, transformative advances in gene editing and progress in lentiviral gene therapy provide diverse opportunities for genetic strategies to cure SCD. Approaches include hematopoietic gene therapy by globin gene addition, gene editing to correct the SCD mutation, and genetic manipulations to enhance fetal hemoglobin production, a potent modifier of the clinical phenotype. Clinical trials may soon identify efficacious and safe genetic approaches to the ultimate goal of cure for SCD.

An Optimized Lentiviral Vector Efficiently Corrects the Human Sickle Cell Disease Phenotype

Autologous transplantation of hematopoietic stem cells transduced with a lentiviral vector (LV) expressing an anti-sickling HBB variant is a potential treatment for sickle cell disease (SCD). With a clinical trial as our ultimate goal, we generated LV constructs containing an anti-sickling HBB transgene (HBBAS3), a minimal HBB promoter, and different combinations of DNase I hypersensitive sites (HSs) from the locus control region (LCR). Hematopoietic stem progenitor cells (HSPCs) from SCD patients were transduced with LVs containing either HS2 and HS3 (β-AS3) or HS2, HS3, and HS4 (β-AS3 HS4). The inclusion of the HS4 element drastically reduced vector titer and infectivity in HSPCs, with negligible improvement of transgene expression. Conversely, the LV containing only HS2 and HS3 was able to efficiently transduce SCD bone marrow and Plerixafor-mobilized HSPCs, with anti-sickling HBB representing up to ∼60% of the total HBB-like chains. The expression of the anti-sickling HBB and the reduced incorporation of the β^S-chain in hemoglobin tetramers allowed up to 50% reduction in the frequency of RBC sickling under hypoxic conditions. Together, these results demonstrate the ability of a high-titer LV to express elevated levels of a potent anti-sickling HBB transgene ameliorating the SCD cell phenotype.

Genetic therapies for sickle cell disease

After decades with few novel therapeutic options for sickle cell disease (SCD), autologous hematopoietic stem cell (HSC) based genetic therapies including lentiviral gene therapy (GT), and genome editing (GE) now appear imminent. Lentiviral GT has advanced considerably in the past decade with promising clinical trial results in multiple disorders. For β-hemoglobinopathies, GT strategies of gene addition and fetal hemoglobin induction through BCL11A regulation are both being evaluated in open clinical trials. GE techniques offer the possibility of a nonviral curative approach, either through sickle hemoglobin mutation repair or fetal hemoglobin elevation. Although GE currently remains at the preclinical stage, multiple clinical trials will likely open soon. In addition to reviewing current strategies for GT and GE, this review highlights important next steps toward optimization of these therapies. All autologous cell-based genetic therapies rely on safely obtaining an adequate yield of autologous HSCs for genetic modification and transplantation. HSC collection is uniquely challenging in SCD. Peripheral mobilization with plerixafor has recently emerged as a promising approach. The acute and long-term toxicities associated with myeloablative conditioning are risks that may not be acceptable to a significant number of SCD patients, highlighting the need for novel conditioning regimens. Finally, increasing availability of autologous genetic therapies will require comprehensive and collaborative discussions regarding cost and access for SCD patients, at individual centers and worldwide.

Induction of fetal hemoglobin synthesis by CRISPR/Cas9-mediated editing of the human β-globin locus

Naturally occurring, large deletions in the β-globin locus result in hereditary persistence of fetal hemoglobin, a condition that mitigates the clinical severity of sickle cell disease (SCD) and β-thalassemia. We designed a clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated protein 9 (Cas9) (CRISPR/Cas9) strategy to disrupt a 13.6-kb genomic region encompassing the δ- and β-globin genes and a putative γ-δ intergenic fetal hemoglobin (HbF) silencer. Disruption of just the putative HbF silencer results in a mild increase in γ-globin expression, whereas deletion or inversion of a 13.6-kb region causes a robust reactivation of HbF synthesis in adult erythroblasts that is associated with epigenetic modifications and changes in chromatin contacts within the β-globin locus. In primary SCD patient-derived hematopoietic stem/progenitor cells, targeting the 13.6-kb region results in a high proportion of γ-globin expression in erythroblasts, increased HbF synthesis, and amelioration of the sickling cell phenotype. Overall, this study provides clues for a potential CRISPR/Cas9 genome editing approach to the therapy of β-hemoglobinopathies.

At least 20% donor myeloid chimerism is necessary to reverse the sickle phenotype after allogeneic HSCT

Novel curative therapies using genetic transfer of normal globin-producing genes into autologous hematopoietic stem cells (HSCs) are in clinical trials for patients with sickle cell disease (SCD). The percentage of transferred globin necessary to cure SCD is currently not known. In the setting of allogeneic nonmyeloablative HSC transplants (HSCTs), stable mixed chimerism is sufficient to reverse the disease. We regularly monitored 67 patients after HSCT. After initially robust engraftment, 3 of these patients experienced declining donor myeloid chimerism (DMC) levels with eventual return of disease. From this we discovered that 20% DMC is necessary to reverse the sickle phenotype. We subsequently developed a mathematical model to test the hypothesis that the percentage of DMC necessary is determined solely by differences between donor and recipient red blood cell (RBC) survival times. In our model, the required 20% DMC can be entirely explained by the large differences between donor and recipient RBC survival times. Our model predicts that the requisite DMC and therefore necessary level of transferred globin is lowest in patients with the highest reticulocyte counts and concomitantly shortened RBC lifespans.

Curative approaches for sickle cell disease: A review of allogeneic and autologous strategies

Despite sickle cell disease (SCD) first being reported >100years ago and molecularly characterized >50years ago, patients continue to experience severe morbidity and early mortality. Although there have been substantial clinical advances with immunizations, penicillin prophylaxis, hydroxyurea treatment, and transfusion therapy, the only cure that can be offered is hematopoietic stem cell transplantation (HSCT). In this work, we summarize the various allogeneic curative approaches reported to date and discuss open and upcoming clinical research protocols. Then we consider gene therapy and gene editing strategies that may enable cure based on autologous HSCs.

Extinction rates in tumour public goods games

Cancer evolution and progression are shaped by cellular interactions and Darwinian selection. Evolutionary game theory incorporates both of these principles, and has been proposed as a framework to understand tumour cell population dynamics. A cornerstone of evolutionary dynamics is the replicator equation, which describes changes in the relative abundance of different cell types, and is able to predict evolutionary equilibria. Typically, the replicator equation focuses on differences in relative fitness. We here show that this framework might not be sufficient under all circumstances, as it neglects important aspects of population growth. Standard replicator dynamics might miss critical differences in the time it takes to reach an equilibrium, as this time also depends on cellular turnover in growing but bounded populations. As the system reaches a stable manifold, the time to reach equilibrium depends on cellular death and birth rates. These rates shape the time scales, in particular, in coevolutionary dynamics of growth factor producers and free-riders. Replicator dynamics might be an appropriate framework only when birth and death rates are of similar magnitude. Otherwise, population growth effects cannot be neglected when predicting the time to reach an equilibrium, and cell-type-specific rates have to be accounted for explicitly.

Gene Therapy for β-Hemoglobinopathies

β-Thalassemia and sickle cell disease (SCD) are the world's two most widely disseminated hereditary hemoglobinopathies. β-Thalassemia originated in the Mediterranean, Middle Eastern, and Asian regions, and SCD originated in central Africa. However, subsequent population migration means that these two diseases are now global and thus constitute a growing health problem in many countries. Despite remarkable improvements in medical care for patients with β-hemoglobinopathies, there is still only one definitive treatment option: allogeneic hematopoietic stem cell (HSC) transplantation. The development of gene therapy for β-hemoglobinopathies has been justified by (1) the limited availability of human leukocyte antigen (HLA)-identical donors, (2) the narrow window of application of HSC transplantation to the youngest patients, and (3) recent advances in HSC-based gene therapy. The huge ongoing efforts in translational medicine and the high number of related publications show that gene therapy has the potential to become the treatment of choice for patients who lack either an HLA genoidentical sibling or an alternative, medically acceptable donor. In this dynamic scientific context, we first summarize the main steps toward clinical translation of this therapeutic approach and then discuss novel lentiviral- and genome editing-based treatment strategies for β-hemoglobinopathies.

Long-Term Engraftment and Fetal Globin Induction upon BCL11A Gene Editing in Bone-Marrow-Derived CD34+ Hematopoietic Stem and Progenitor Cells

To develop an effective and sustainable cell therapy for sickle cell disease (SCD), we investigated the feasibility of targeted disruption of the BCL11A gene, either within exon 2 or at the GATAA motif in the intronic erythroid-specific enhancer, using zinc finger nucleases in human bone marrow (BM) CD34⁺ hematopoietic stem and progenitor cells (HSPCs). Both targeting strategies upregulated fetal globin expression in erythroid cells to levels predicted to inhibit hemoglobin S polymerization. However, complete inactivation of BCL11A resulting from bi-allelic frameshift mutations in BCL11A exon 2 adversely affected erythroid enucleation. In contrast, bi-allelic disruption of the GATAA motif in the erythroid enhancer of BCL11A did not negatively impact enucleation. Furthermore, BCL11A exon 2-edited BM-CD34⁺ cells demonstrated a significantly reduced engraftment potential in immunodeficient mice. Such an adverse effect on HSPC function was not observed upon BCL11A erythroid-enhancer GATAA motif editing, because enhancer-edited CD34⁺ cells achieved robust long-term engraftment and gave rise to erythroid cells with elevated levels of fetal globin expression when chimeric BM was cultured ex vivo. Altogether, our results support further clinical development of the BCL11A erythroid-specific enhancer editing in BM-CD34⁺ HSPCs as an autologous stem cell therapy in SCD patients.

Lineage-specific BCL11A knockdown circumvents toxicities and reverses sickle phenotype

Reducing expression of the fetal hemoglobin (HbF) repressor BCL11A leads to a simultaneous increase in γ-globin expression and reduction in β-globin expression. Thus, there is interest in targeting BCL11A as a treatment for β-hemoglobinopathies, including sickle cell disease (SCD) and β-thalassemia. Here, we found that using optimized shRNAs embedded within an miRNA (shRNAmiR) architecture to achieve ubiquitous knockdown of BCL11A profoundly impaired long-term engraftment of both human and mouse hematopoietic stem cells (HSCs) despite a reduction in nonspecific cellular toxicities. BCL11A knockdown was associated with a substantial increase in S/G2-phase human HSCs after engraftment into immunodeficient (NSG) mice, a phenotype that is associated with HSC exhaustion. Lineage-specific, shRNAmiR-mediated suppression of BCL11A in erythroid cells led to stable long-term engraftment of gene-modified cells. Transduced primary normal or SCD human HSCs expressing the lineage-specific BCL11A shRNAmiR gave rise to erythroid cells with up to 90% reduction of BCL11A protein. These erythrocytes demonstrated 60%-70% γ-chain expression (vs. < 10% for negative control) and a corresponding increase in HbF. Transplantation of gene-modified murine HSCs from Berkeley sickle cell mice led to a substantial improvement of sickle-associated hemolytic anemia and reticulocytosis, key pathophysiological biomarkers of SCD. These data form the basis for a clinical trial application for treating sickle cell disease.