Fetal hemoglobin in sickle cell anemia: Molecular characterization of the unusually high fetal hemoglobin phenotype in African Americans

HbF inhibits deoxy-HbS polymerization. Patients with elevated HbF have fewer vasoocclusive complications and prolonged survival [1]. Three major HbF QTL are known. The C>T polymorphism (rs7482144) at nucleotide –158 upstream of HBG2 is associated with increased HbF in some SCA patients [2]. Polymorphisms in intron 2 of BCL11A represented by rs766432 was associated with HbF in healthy Northern Europeans [3], African Americans with SCA [4, 5], Chinese with β-thalassemia trait and Thai’s with HbE-β thalassemia [5]. BCL11A polymorphisms correlate highly with HbF levels in SCA, accounting for 7–12% of the HbF variance [6]. The HMIP polymorphisms are distributed in three LD blocks [7]. HMIP block 2 represented by rs9399137 is most significantly associated with HbF expression and might function as a distal regulatory element [8,9]. We studied a selected group of 20 African American SCA patients with exceptionally high HbF (mean 17.2%) which differed by more than 4 times the standard deviation of 30 other patients with low HbF (mean 5.0%; Table IA). All study subjects’ HBB underwent nucleotide sequencing to ascertain that they were HbS homozygotes. Multiplex ligation-dependent probe amplification (MLPA) was carried out to ensure that they did not harbor hereditary persistence of fetal hemoglobin (HPFH) 1, HPFH 2, Black (δβ)0- and Black (Aγδβ)0-thalassemia deletions [10]. Furthermore, their HBG2 and HBG1 promoters were also sequenced to be certain that they did not have promoter HPFH single nucleotide mutations [11,12]. They were unlikely to be on hydroxyurea based on their MCV being less than 100 fL. Table IA Hematologic results of high and low HbF study groups. In addition we conducted a subset analyses in 56 patients with unusually high HbF (mean 20.7%) which differed by more than 11 times the standard deviation of 489 patients with low HbF (mean 3.1%; table IB). These patients were selected from 1,086 subjects from the Cooperative Study of Sickle Cell Disease (CSSCD) who were previously investigated in a genome-wide association study (GWAS) of HbF [13]. Table IB Hematologic results of high and low HbF CSCCD groups. The MAF of rs7482144, also known as the Xmn I site, on chromosome 11p15 in the unusually high HbF group (10%) is not significantly different from that in the low HbF group (8%) as shown in Table II. The SNP rs5006884, a missense mutation (CTC>TTC or Leu172Phe) in OR51B6 (http://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi?rs=5006884) on chromosome 11p15 was reported to be associated with HbF in SCA in a GWAS [13]. In the present study, the MAF of rs5006884 in the unusually high HbF group (2%) is actually less than that in the low HbF group (10%), even though the difference is not statistically significant (Table II). It should be recognized that the small sample size in the current study does not afford sufficient power to discern possible MAF differences if present at these 2 SNPs between these 2 groups of patients. Table II HbF QTL minor allele frequencies in high and low HbF groups The SNP rs7482144 did not reach genome-wide significance in previous GWAS of HbF in the full CSCCD cohort [13]. In the present study with subsets from the CSCCD cohort, rs7482144 has a significantly higher MAF in the unusually high HbF group (30%) compared to the low HbF group (10%), P = 0.002. The MAF of rs5006884 in OR51B6 in the high HbF group (22%) is significantly higher than that in the low HbF group (10%), P = 0.00055. This SNP reached genome-wide significance in the GWAS of the full CSSCD cohort, however the mean HbF levels in homozygotes for this variant was only 10.6%. The present subset analysis shows that for SCA patients with unusually high HbF the frequency of this mutation is comparable to the general population (Table II). To examine the BCL11A QTL, three SNPs, rs766432, rs4671393 and rs11886868 were chosen for genotyping [6]. These 3 SNPs are in strong LD. Only data on rs766432 that is most highly correlated with HbF in SCA patients are presented. The MAF in the unusually high HbF group (45%) is significantly higher than that in the low HbF group (25%), P = 0.05 (Table II). The MAF in the subgroups from CSCCD were similar: 47% in the high and 20% in the low HbF groups, P = 6.399×10−10. The allele C of this SNP was associated with increased levels of HbF [13] and the mean HbF levels in subjects homozygous for the C allele was 8.26%. The present subset analysis in the CSSCD subgroups shows that in patients with unusually high HbF the frequency of this mutation is almost twice that of the full cohort. The QTL in HMIP is best represented by rs9399137 [9]. The MAF of rs9399137 among SCA patients of African descent without European admixture was reported to be 1–2% [14,15]. A GWAS on over 800 African American SCA patients did not detect genome-wide significance of association of this QTL with HbF [13]. The functional motif for this QTL is most likely a 3-bp (TAC) deletion which is in complete LD with the minor allele of rs9399137 [9]. The frequency of this 3-bp deletion is 23% in non-African HapMap populations, but only 5% in Africans [9]. In the current study, the MAF of rs9399137 in the African American SCA patients with unusually high HbF is 18%, significantly higher than that with low HbF (3%), P = 0.02 (Table II). Furthermore, the 3-bp deletion as reported by Farrell et al [9] was found for the first time in these African American SCA patients and it is in complete LD with the minor allele of rs9399137. Among the subset of CSSCD patients, the MAF of rs9399137 in the unusually high HbF group (9%) is also significantly higher than that in the low HbF group (3%), P = 0.006 (Table II). These results raise the possibility that some African American SCA patients with markedly elevated HbF might have inherited the minor allele of chromosome 6q23 QTL due to European genetic admixture. The minor T allele of rs7775698 tags either an ancestral T nucleotide found mostly in African populations, or a 3-bp deletion often found in European and Chinese populations [9]. In the present study, the MAF of rs7775698 in the high HbF group (40%) is higher than that in the low HbF group (21%). But the difference is not statistically significant (P = 0.07). Among the CSCCD patients, the MAF of rs7775698 in the high HbF group is 16% compared to 19% in the low HbF group (Table II). We found a 2-bp (CC) deletion plus an (A) insertion 19 bp downstream of rs9399137. This deletion/insertion was present in 33% of chromosome 6 in both high and low HbF groups. It is unlikely that it plays a significant functional role in modulating HbF expression. Its relatively high frequency in both groups makes it a probable haplotype marker in African American SCA patients. Two of the 3 known HbF major QTLs are present in the unusually high HbF patient study group at a frequency significantly higher than those with low HbF. Their cumulative effect was estimated by assigning to each minor allele at each locus a score of 1, and adding the scores to generate the Number of Minor Allele Present (NOMAP) [16]. There was an average of 2 minor alleles in the unusually high HbF group compared with 1 in the low HbF group (P = 0.001) as shown in Figure 1. Figure 1 To explore other possible genetic loci that may modulate HbF expression, we undertook nucleotide sequencing of a 14.1 kb DNA fragment between HBG1 and HBD in 15 high and 15 low HbF patients. This DNA fragment was chosen because BCL11A binds to this intergenic region and it also encompasses the Corfu deletion that in homozygotes is characterized by markedly elevated HbF [17, 18]. Thirty eight SNPs were found in both the high and low HbF groups (Supplementary Table I). SNP rs10128558 as described by Galarneau et al [19] was present in our study groups and found to be in complete LD with Xmn I polymorphism. Its MAF in the high HbF group (10%) is not different from the low HbF group (8%). Based on in silico analysis by TFSEARCH (threshold score 85), 22 of the 38 SNPs are associated with alterations in transcription factor binding including erythroid specific transcription factors. In addition, 4 SNPs have significantly higher major allele frequency in the high compared to the low HbF groups (P < 0.05). Three of these 4 SNPs result in alteration in transcription factor binding sites (Supplementary Table I). This study based on a small cohort of carefully selected African American SCA patients with unusually high HbF revealed that the MAF for rs766432 (BCL11A), rs9399137 and 3-bp deletion both within HMIP are much higher than that found in patients with low HbF. These findings also raise the possibility that some African American patients with markedly elevated HbF might have inherited the minor allele of chromosome 6q23 QTL due to European genetic admixture. Validation of these findings in a larger patient cohort and further functional investigations into these and other polymorphisms are warranted.