Most genome editing analyses to date are based on quantifying small insertions and deletions. Here, we show that CRISPR-Cas9 genome editing can induce large gene modifications, such as deletions, insertions, and complex local rearrangements in different primary cells and cell lines. We analyzed large deletion events in hematopoietic stem and progenitor cells (HSPCs) using different methods, including clonal genotyping, droplet digital polymerase chain reaction, single-molecule real-time sequencing with unique molecular identifier, and long-amplicon sequencing assay. Our results show that large deletions of up to several thousand bases occur with high frequencies at the Cas9 on-target cut sites on the HBB (11.7 to 35.4%), HBG (14.3%), and BCL11A (13.2%) genes in HSPCs and the PD-1 (15.2%) gene in T cells. Our findings have important implications to advancing genome editing technologies for treating human diseases, because unintended large gene modifications may persist, thus altering the biological functions and reducing the available therapeutic alleles.

Red blood cell (RBC) deformability is a valuable hemorheological biomarker that can be used to assess the clinical status and response to therapy of individuals with sickle cell disease (SCD). RBC deformability has been measured by ektacytometry for decades, which uses shear or osmolar stress. However, ektacytometry is a population based measurement that does not detect small-fractions of abnormal RBCs. A single cell-based, functional RBC deformability assay would complement ektacytometry and provide additional information. Here, we tested the relative merits of the OcclusionChip, which measures RBC deformability by microcapillary occlusion, and ektacytometry. We tested samples containing glutaraldehyde-stiffened RBCs for up to 1% volume fraction; ektacytometry detected no significant change in Elongation Index (EI), while the OcclusionChip showed significant differences in Occlusion Index (OI). OcclusionChip detected a significant increase in OI in RBCs from an individual with sickle cell trait (SCT) and from a subject with SCD who received allogeneic hematopoietic stem cell transplant (HSCT), as the sample was taken from normoxic (pO2:159 mmHg) to physiologic hypoxic (pO2:45 mmHg) conditions. Oxygen gradient ektacytometry detected no difference in EI for SCT or HSCT. These results suggest that the single cell-based OcclusionChip enables detection of sickle hemoglobin (HbS)-related RBC abnormalities in SCT and SCD, particularly when the HbS level is low. We conclude that the OcclusionChip is complementary to the population based ektacytometry assays, and providing additional sensitivity and capacity to detect modest abnormalities in red cell function or small populations of abnormal red cells.

Human genetic variation has enabled the identification of several key regulators of fetal-to-adult hemoglobin switching, including BCL11A, resulting in therapeutic advances. However, despite the progress made, limited further insights have been obtained to provide a fuller accounting of how genetic variation contributes to the global mechanisms of fetal hemoglobin (HbF) gene regulation. Here, we have conducted a multi-ancestry genome-wide association study of 28,279 individuals from several cohorts spanning 5 continents to define the architecture of human genetic variation impacting HbF. We have identified a total of 178 conditionally independent genome-wide significant or suggestive variants across 14 genomic windows. Importantly, these new data enable us to better define the mechanisms by which HbF switching occurs in vivo. We conduct targeted perturbations to define BACH2 as a new genetically-nominated regulator of hemoglobin switching. We define putative causal variants and underlying mechanisms at the well-studied BCL11A and HBS1L-MYB loci, illuminating the complex variant-driven regulation present at these loci. We additionally show how rare large-effect deletions in the HBB locus can interact with polygenic variation to influence HbF levels. Our study paves the way for the next generation of therapies to more effectively induce HbF in sickle cell disease and β-thalassemia.

We read with great interest the recently published report of Nardo‐Marino et al. 1 The authors report results of a study designed to: (1) investigate the relationship between oxygen gradient ektacytometry‐derived biomarkers and two established markers of disease severity in children with sickle cell anaemia (SCA), (2) examine the impact of hydroxycarbamide (hydroxyurea [HU]) on red blood cell (RBC) deformability measured with oxygen gradient ektacytometry and (3) investigate the ability of oxygen gradient ektacytometry‐derived parameters to predict pain as a clinical outcome in SCA. The authors found that all oxygen gradient ektacytometry‐derived parameters correlated with known modifiers of disease severity, and that specifically the point of sickling (PoS) correlated with transcranial Doppler velocities. On the other hand, they did not find that oxygen gradient ektacytometry parameters independently predicted pain events or response to HU treatment. Therefore, the authors concluded that oxygen gradient ektacytometry could not predict pain events in SCA and questioned its use in clinical practice.

In sickle cell disease (SCD), higher whole blood viscosity is a risk factor for vaso-occlusive crisis, avascular necrosis, and proliferative retinopathy. Blood viscosity is strongly impacted by hemoglobin (Hb) levels and red blood cell (RBC) deformability. Voxelotor is a hemoglobin S (HbS) polymerization inhibitor with anti-sickling properties that increases the Hb affinity for oxygen, thereby reducing HbS polymerization. In clinical trials, voxelotor increased Hb by an average of 1g/dl, creating concern that this rise in Hb could increase viscosity, particularly when the drug was cleared. To investigate this potential rebound hyperviscosity effect, we treated SCD mice with GBT1118, a voxelotor analog, and stopped the treatment to determine the effect on blood viscosity and RBC deformability under a range of oxygen concentrations. GBT1118 treatment increased Hb, improved RBC deformability by increasing the elongation index under normoxic (EImax) and hypoxic conditions (EImin), and decreased the point of sickling (PoS) without increasing blood viscosity. The anti-sickling effects and improvement of RBC deformability balanced the effect of increased Hb such that there was no increase in blood viscosity. Forty-eight hours after ceasing GBT1118, Hb declined from the rise induced by treatment, viscosity did not increase, and EImin remained elevated compared to control animals. Hb and PoS were not different from control animals, suggesting a return to native oxygen affinity and clearance of the drug. RBC deformability did not return to baseline, suggesting some residual rheological improvement. These data suggest that concerns regarding viscosity rise above pre-treatment levels upon sudden cessation of voxelotor are not warranted.

Individuals with monogenic disorders can experience variable phenotypes that are influenced by genetic variation. To investigate this in sickle cell disease (SCD), we performed whole-genome sequencing (WGS) of 722 individuals with hemoglobin HbSS or HbSb0-thalassemia from Baylor College of Medicine and from the St. Jude Children's Research Hospital Sickle Cell Clinical Research and Intervention Program (SCCRIP) longitudinal cohort study. We developed pipelines to identify genetic variants that modulate sickle hemoglobin polymerization in red blood cells and combined these with painassociated variants to build a polygenic score (PGS) for acute vaso-occlusive pain (VOP). Overall, we interrogated the a-thalassemia deletion -α3.7 and 133 candidate singlenucleotide polymorphisms (SNPs) across 66 genes for associations with VOP in 327 SCCRIP participants followed longitudinally over 6 years. Twenty-one SNPs in 9 loci were associated with VOP, including 3 (BCL11A, MYB, and the b-like globin gene cluster) that regulate erythrocyte fetal hemoglobin (HbF) levels and 6 (COMT, TBC1D1, KCNJ6, FAAH, NR3C1, and IL1A) that were associated previously with various pain syndromes. An unweighted PGS integrating all 21 SNPs was associated with the VOP event rate (estimate, 0.35; standard error, 0.04; P = 5.9 × 10-14) and VOP event occurrence (estimate, 0.42; standard error, 0.06; P = 4.1 × 10-13). These associations were stronger than those of any single locus. Our findings provide insights into the genetic modulation of VOP in children with SCD. More generally, we demonstrate the utility of WGS for investigating genetic contributions to the variable expression of SCD-associated morbidities.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Sickle cell disease (SCD) is associated with severe morbidity and early mortality. Two large population studies found an increased risk for leukemia in individuals with SCD. Notably, while the relative risk of leukemia development is high, the absolute risk is low in individuals with SCD who do not receive cell-based therapies. However, the risk of leukemia in SCD is high after graft rejection and with gene therapy. Clonal hematopoiesis (CH) is a well-recognized premalignant condition in the general population and in patients after high-dose myelotoxic therapies. Recent studies suggest that CH may be more common in SCD than in the general population, outside the cell-based therapy setting. Here, we review risk factors for CH and progression to leukemia in SCD. We surmise why patients with SCD are at an increased risk for CH and why leukemia incidence is unexpectedly high after graft rejection and gene therapy for SCD. Currently, we are unable to reliably assess genetic risk factors for leukemia development after curative therapies for SCD. Given our current knowledge, we recommend counseling patients about leukemia risk and discussing the importance of an individualized benefit/risk assessment that incorporates leukemia risk in patients undergoing curative therapies for SCD.