Fetal hemoglobin (HbF) is a strong modifier of sickle cell disease (SCD) severity and is associated with 3 common genetic loci. Quantifying the genetic effects of the 3 loci would specifically address the benefits of HbF increases in patients. Here, we have applied statistical methods using the most representative variants: rs1427407 and rs6545816 in BCL11A, rs66650371 (3-bp deletion) and rs9376090 in HMIP-2A, rs9494142 and rs9494145 in HMIP-2B, and rs7482144 (Xmn1-HBG2 in the β-globin locus) to create g(HbF), a genetic quantitative variable for HbF in SCD. Only patients aged ≥5 years with complete genotype and HbF data were studied. Five hundred eighty-one patients with hemoglobin SS (HbSS) or HbSβ0 thalassemia formed the discovery cohort. Multiple linear regression modeling rationalized the 7 variants down to 4 markers (rs6545816, rs1427407, rs66650371, and rs7482144) each independently contributing HbF-boosting alleles, together accounting for 21.8% of HbF variability (r2) in the HbSS or HbSβ0 patients. The model was replicated with consistent r2 in 2 different cohorts: 27.5% in HbSC patients (N = 186) and 23% in 994 Tanzanian HbSS patients. g(HbF), our 4-variant model, provides a robust approach to account for the genetic component of HbF in SCD and is of potential utility in sickle genetic and clinical studies.

Genetic variants at three quantitative trait loci (QTL) for fetal haemoglobin (HbF), BCL11A, HBS1L-MYB and the β-globin gene cluster, have attracted interest as potential targets of therapeutic strategies for HbF reactivation in sickle cell anaemia (SCA). We carried out the first systematic evaluation of critical single nucleotide polymorphisms at these disease modifier loci in Nigerian patients with SCA. Common variants for BCL11A and HBS1L-MYB were strongly associated with HbF levels. At both loci, secondary association signals were detected, illustrating the mapping resolution attainable in this population. For BCL11A, the two independent sites of association were represented by rs1427407 (primary site, p = 7.0 x 10(-10)) and rs6545816 (secondary site, conditioned on rs1427407: p = 0.02) and for HBS1L-MYB by rs9402686 (HMIP-2B, p = 1.23 x 10(-4)) and rs66650371 (HMIP-2A, p = 0.002). Haplotype analysis revealed similarities in the genetic architecture of BCL11A and HBS1L-MYB in Nigerian patients. Variants at both loci also alleviated anaemia. The variant allele for the γ globin gene promoter polymorphism XmnI-HBG2 was too infrequent in our patients to be evaluated in this relatively small study. Studying the large and diverse SCA patient populations in African countries such as Nigeria will be key for a clearer understanding of how these loci work and for the discovery of new disease modifier genes.