The genetic dissection of fetal haemoglobin persistence in sickle cell disease in Nigeria.

The clinical severity of sickle cell disease (SCD) is strongly influenced by the level of fetal haemoglobin (HbF) persistent in each patient. Three major HbF loci (BCL11A, HBS1L-MYB, and Xmn1-HBG2) have been reported, but a considerable hidden heritability remains. We conducted a genome-wide association study for HbF levels in 1006 Nigerian patients with SCD (HbSS/HbSß0), followed by a replication and meta-analysis exercise in four independent SCD cohorts (3,582 patients). To dissect association signals at the major loci, we performed stepwise conditional and haplotype association analyses and included public functional annotation datasets. Association signals were detected for BCL11A (lead SNP rs6706648, ß = -0.39, P = 4.96 × 10-34) and HBS1L-MYB (lead SNP rs61028892, ß = 0.73, P = 1.18 × 10-9), whereas the variant allele for Xmn1-HBG2 was found to be very rare. In addition, we detected three putative new trait-associated regions. Genetically, dissecting the two major loci BCL11A and HBS1L-MYB, we defined trait-increasing haplotypes (P < 0.0001) containing so far unidentified causal variants. At BCL11A, in addition to a haplotype harbouring the putative functional variant rs1427407-'T', we identified a second haplotype, tagged by the rs7565301-'A' allele, where a yet-to-be-discovered causal DNA variant may reside. Similarly, at HBS1L-MYB, one HbF-increasing haplotype contains the likely functional small indel rs66650371, and a second tagged by rs61028892-'C' is likely to harbour a presently unknown functional allele. Together, variants at BCL11A and HBS1L-MYB SNPs explained 24.1% of the trait variance. Our findings provide a path for further investigation of the causes of variable fetal haemoglobin persistence in sickle cell disease.

Genome wide association study of silent cerebral infarction in sickle cell disease (HbSS and HbSC).

Not available.

F cell numbers are associated with an X-linked genetic polymorphism and correlate with haematological parameters in patients with sickle cell disease.

Patients with sickle cell disease (SCD) with high fetal haemoglobin (HbF) tend to have a lower incidence of complications and longer survival due to inhibition of deoxyhaemoglobin S (HbS) polymerisation by HbF. HbF-containing cells, namely F cells, are strongly influenced by genetic factors. We measured the percentage of F cells (Fcells%) in 222 patients with SCD to evaluate the association of (i) Fcells% with genetic HbF-modifier variants and (ii) Fcells% with haematological parameters. There was a different distribution of Fcells% in females compared to males. The association of the B-cell lymphoma/leukaemia 11A (BCL11A) locus with Fcells% (ß = 8·238; P < 0·001) and with HbF% (ß = 2·490; P < 0·001) was significant. All red cell parameters except for Hb and mean corpuscular Hb concentration levels in males and females were significantly different. The Fcells% was positively associated with mean cell Hb, mean cell volume and reticulocytes. To explain the significant gender difference in Fcells%, we tested for associations with single nucleotide polymorphisms on the X chromosomal region Xp22.2, where a genetic determinant of HbF had been previously hypothesised. We found in males a significant association with a SNP in FERM and PDZ domain-containing protein 4 (FRMPD4) and adjacent to male-specific lethal complex subunit 3 (MSL3). Thus, we have identified an X-linked locus that could account for a significant fraction of the Fcells% variation in patients with SCD.

A survey of genetic fetal-haemoglobin modifiers in Nigerian patients with sickle cell anaemia.

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.

g(HbF): a genetic model of fetal hemoglobin in sickle cell disease.

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 association of fetal-hemoglobin levels in individuals with sickle cell disease in Tanzania maps to conserved regulatory elements within the MYB core enhancer.

BACKGROUND: Common genetic variants residing near upstream regulatory elements for MYB, the gene encoding transcription factor cMYB, promote the persistence of fetal hemoglobin (HbF) into adulthood. While they have no consequences in healthy individuals, high HbF levels have major clinical benefits in patients with sickle cell disease (SCD) or ß thalassemia. Here, we present our detailed investigation of HBS1L-MYB intergenic polymorphism block 2 (HMIP-2), the central component of the complex quantitative-trait locus upstream of MYB, in 1,022 individuals with SCD in Tanzania. METHODS: We have looked at 1022 individuals with HbSS or HbS/ß(0) in Tanzania. In order to achieve a detailed analysis of HMIP-2, we performed targeted genotyping for a total of 10 SNPs and extracted additional 528 SNPs information from a genome wide scan involving the same population. Using MACH, we utilized the existing YRI data from 1000 genomes to impute 54 SNPs situated within HIMP-2. RESULTS: Seven HbF-increasing, low-frequency variants (ß > 0.3, p < 10(-5), f ≤ 0.05) were located in two partially-independent sub-loci, HMIP-2A and HMIP-2B. The spectrum of haplotypes carrying such alleles was diverse when compared to European and West African reference populations: we detected one such haplotype at sub-locus HMIP-2A, two at HMIP-2B, and a fourth including high-HbF alleles at both sub-loci ('Eurasian' haplotype clade). In the region of HMIP-2A a putative functional variant (a 3-bp indel) has been described previously, but no such candidate causative variant exists at HMIP-2B. Extending our dataset through imputation with 1000 Genomes, whole-genome-sequence data, we have mapped peak association at HMIP-2B to an 11-kb region around rs9494145 and rs9483788, flanked by two conserved regulatory elements for MYB. CONCLUSIONS: Studies in populations from the African continent provide distinct opportunities for mapping disease-modifying genetic loci, especially for conditions that are highly prevalent there, such as SCD. Population-genetic characteristics of our cohort, such as ethnic diversity and the predominance of shorter, African-type haplotypes, can add to the power of such studies.

Next generation sequencing identifies a novel rearrangement in the HBB cluster permitting to-the-base characterization.

Genetic testing for hemoglobinopathies is required for prenatal diagnosis, understanding complex cases where multiple pathogenic variants may be present or investigating cases of unexplained anemia. Characterization of disease causing variants that range from single base changes to large rearrangements may require several different labor-intensive methodologies. Multiplex ligation probe amplification analysis is the current method used to detect indels, but the technique does not characterize the breakpoints or detect balanced translocations. Here, we describe a next-generation sequencing (NGS) method that is able to identify and characterize a novel rearrangement of the HBB cluster responsible for ÎµÎ³Î´ß thalassemia in an English family. The structural variant involved a 59.0 kb inversion encompassing HBG2 exon 3, HBG1, HBD, HBB, and OR51V1, juxtaposed by a deletion of 122.6 kb including 82 bp of the inverted sequence, HBG2 exon 1 and 2, HBE, and the ß-locus control region. Identification of reads spanning the breakpoints provided to-the-base resolution of the rearrangement, subsequently confirmed by gap-PCR and Sanger sequence analysis. The same rearrangement, termed Inv-Del English V ÎµÎ³Î´ß thalassemia (HbVar 2935), was identified in two other unrelated English individuals with a similar hematological phenotype. Our NGS approach should be applicable as a diagnostic tool for other disorders.

Genome wide association study of fetal hemoglobin in sickle cell anemia in Tanzania.

BACKGROUND: Fetal hemoglobin (HbF) is an important modulator of sickle cell disease (SCD). HbF has previously been shown to be affected by variants at three loci on chromosomes 2, 6 and 11, but it is likely that additional loci remain to be discovered. METHODS AND FINDINGS: We conducted a genome-wide association study (GWAS) in 1,213 SCA (HbSS/HbSß0) patients in Tanzania. Genotyping was done with Illumina Omni2.5 array and imputation using 1000 Genomes Phase I release data. Association with HbF was analysed using a linear mixed model to control for complex population structure within our study. We successfully replicated known associations for HbF near BCL11A and the HBS1L-MYB intergenic polymorphisms (HMIP), including multiple independent effects near BCL11A, consistent with previous reports. We observed eight additional associations with P<10(-6). These associations could not be replicated in a SCA population in the UK. CONCLUSIONS: This is the largest GWAS study in SCA in Africa. We have confirmed known associations and identified new genetic associations with HbF that require further replication in SCA populations in Africa.

Global genetic architecture of an erythroid quantitative trait locus, HMIP-2.

HMIP-2 is a human quantitative trait locus affecting peripheral numbers, size and hemoglobin composition of red blood cells, with a marked effect on the persistence of the fetal form of hemoglobin, HbF, in adults. The locus consists of multiple common variants in an enhancer region for MYB (chr 6q23.3), which encodes the hematopoietic transcription factor cMYB. Studying a European population cohort and four African-descended groups of patients with sickle cell anemia, we found that all share a set of two spatially separate HbF-promoting alleles at HMIP-2, termed "A" and "B." These typically occurred together ("A-B") on European chromosomes, but existed on separate homologous chromosomes in Africans. Using haplotype signatures for "A" and "B," we interrogated public population datasets. Haplotypes carrying only "A" or "B" were typical for populations in Sub-Saharan Africa. The "A-B" combination was frequent in European, Asian, and Amerindian populations. Both alleles were infrequent in tropical regions, possibly undergoing negative selection by geographical factors, as has been reported for malaria with other hematological traits. We propose that the ascertainment of worldwide distribution patterns for common, HbF-promoting alleles can aid their further genetic characterization, including the investigation of gene-environment interaction during human migration and adaptation.

HBS1L-MYB intergenic variants modulate fetal hemoglobin via long-range MYB enhancers.

Genetic studies have identified common variants within the intergenic region (HBS1L-MYB) between GTP-binding elongation factor HBS1L and myeloblastosis oncogene MYB on chromosome 6q that are associated with elevated fetal hemoglobin (HbF) levels and alterations of other clinically important human erythroid traits. It is unclear how these noncoding sequence variants affect multiple erythrocyte characteristics. Here, we determined that several HBS1L-MYB intergenic variants affect regulatory elements that are occupied by key erythroid transcription factors within this region. These elements interact with MYB, a critical regulator of erythroid development and HbF levels. We found that several HBS1L-MYB intergenic variants reduce transcription factor binding, affecting long-range interactions with MYB and MYB expression levels. These data provide a functional explanation for the genetic association of HBS1L-MYB intergenic polymorphisms with human erythroid traits and HbF levels. Our results further designate MYB as a target for therapeutic induction of HbF to ameliorate sickle cell and ß-thalassemia disease severity.

Genetic determinants of haemolysis in sickle cell anaemia.

Haemolytic anaemia is variable among patients with sickle cell anaemia and can be estimated by reticulocyte count, lactate dehydrogenase, aspartate aminotransferase and bilirubin levels. Using principal component analysis of these measurements we computed a haemolytic score that we used as a subphenotype in a genome-wide association study. We identified in one cohort and replicated in two additional cohorts the association of a single nucleotide polymorphism in NPRL3 (rs7203560; chr16p13·3) (P = 6·04 × 10(-07) ). This association was validated by targeted genotyping in a fourth independent cohort. The HBA1/HBA2 regulatory elements, hypersensitive sites (HS)-33, HS-40 and HS-48 are located in introns of NPRL3. Rs7203560 was in perfect linkage disequilibrium (LD) with rs9926112 (r(2)  = 1) and in strong LD with rs7197554 (r(2)  = 0·75) and rs13336641 (r(2)  = 0·77); the latter is located between HS-33 and HS-40 sites and next to a CTCF binding site. The minor allele for rs7203560 was associated with the -∝(3·7) thalassaemia gene deletion. When adjusting for HbF and ∝ thalassaemia, the association of NPRL3 with the haemolytic score was significant (P = 0·00375) and remained significant when examining only cases without gene deletion∝ thalassaemia (P = 0·02463). Perhaps by independently down-regulating expression of the HBA1/HBA2 genes, variants of the HBA1/HBA2 gene regulatory loci, tagged by rs7203560, reduce haemolysis in sickle cell anaemia.

HbA2 levels in normal adults are influenced by two distinct genetic mechanisms.

Using a genome-wide association study, we found that common inter-individual differences in haemoglobin A(2) (HbA(2) , α(2) δ(2) ) levels are largely governed by genetic factors (42% of variability). The influence of age (1%) and sex (4%) was small. HbA(2) levels were influenced by two loci: the HBS1L-MYB locus on chromosome 6q, which has been shown to have pleiotropic effects on other haematological traits; and a second locus surrounding HBB, the gene encoding ß-globin. Our results suggest that HbA(2) levels in adults are influenced by two different biological processes: one via kinetics of erythropoiesis, and the other, via competition between HBB and HBD activity.

A novel 506kb deletion causing ÎµÎ³Î´ß thalassemia.

We describe a novel deletion causing ÎµÎ³Î´ß thalassemia in a Pakistani family. The Pakistani deletion is 506kb in length, and the second largest ÎµÎ³Î´ß thalassemia deletion reported to date. It removes the entire ß globin gene (HBB) cluster, extending from 431kb upstream to 75kb downstream of the ε globin gene (HBE). The breakpoint junction occurred within a 160bp palindrome embedded in LINE/LTR repeats, and contained a short (9bp) region of direct homology which may have contributed to the recombination event. Characterization of the deletion breakpoints has been particularly challenging due to the complexity of DNA deletion, insertion and inversion, involving a multitude of methodologies, mirroring the changing DNA analysis technologies.

Genetics of fetal hemoglobin in Tanzanian and British patients with sickle cell anemia.

Fetal hemoglobin (HbF, α(2)γ(2)) is a major contributor to the remarkable phenotypic heterogeneity of sickle cell anemia (SCA). Genetic variation at 3 principal loci (HBB cluster on chromosome 11p, HBS1L-MYB region on chromosome 6q, and BCL11A on chromosome 2p) have been shown to influence HbF levels and disease severity in ß-thalassemia and SCA. Previous studies in SCA, however, have been restricted to populations from the African diaspora, which include multiple genealogies. We have investigated the influence of these 3 loci on HbF levels in sickle cell patients from Tanzania and in a small group of African British sickle patients. All 3 loci have a significant impact on the trait in both patient groups. The results suggest the presence of HBS1L-MYB variants affecting HbF in patients who are not tracked well by European-derived markers, such as rs9399137. Additional loci may be identified through independent genome-wide association studies in African populations.

The HBS1L-MYB intergenic interval associated with elevated HbF levels shows characteristics of a distal regulatory region in erythroid cells.

HBS1L-MYB intergenic polymorphism (HMIP) on chromosome 6q23 is associated with elevated fetal hemoglobin levels and has pleiotropic effects on several hematologic parameters. To investigate potential regulatory activity in the region, we have measured sensitivity of the sequences to DNase I cleavage that identified 3 tissue-specific DNase I hypersensitive sites in the core intergenic interval. Chromatin immunoprecipitation with microarray (ChIP-chip) analysis showed strong histone acetylation in a defined interval of 65 kb corresponding to the core HBS1L-MYB intergenic region in primary human erythroid cells but not in non-MYB-expressing HeLa cells. ChIP-chip analysis also identified several potential cis-regulatory elements as strong GATA-1 signals that coincided with the DNase I hypersensitive sites present in MYB-expressing erythroid cells. We suggest that HMIP contains regulatory sequences that could be important in hematopoiesis by controlling MYB expression. This study provides the functional link between genetic association of HMIP with control of fetal hemoglobin and other hematologic parameters. We also present a large-scale analysis of histone acetylation as well as RNA polymerase II and GATA-1 interactions on chromosome 6q, and alpha and beta globin gene loci. The data suggest that GATA-1 regulates numerous genes of various functions on chromosome 6q.

A QTL influencing F cell production maps to a gene encoding a zinc-finger protein on chromosome 2p15.

F cells measure the presence of fetal hemoglobin, a heritable quantitative trait in adults that accounts for substantial phenotypic diversity of sickle cell disease and beta thalassemia. We applied a genome-wide association mapping strategy to individuals with contrasting extreme trait values and mapped a new F cell quantitative trait locus to BCL11A, which encodes a zinc-finger protein, on chromosome 2p15. The 2p15 BCL11A quantitative trait locus accounts for 15.1% of the trait variance.

Two candidate genes for low platelet count identified in an Asian Indian kindred by genome-wide linkage analysis: glycoprotein IX and thrombopoietin.

A genome-wide linkage analysis of platelet count was carried out in a large Asian Indian kindred. Linkage analysis showed one marker (D3S1309) on chromosome 3q with a lod score of 3.26 and another (D3S1282) approximately 30 cM centromeric, with a lod score of 2.52. Multipoint analysis of chromosome 3q identified two peaks with maximum multipoint lod scores of 3.52 and 4.11 under markers D3S1309 and D3S1282, respectively. Two strong candidate genes for platelet variation were identified in the linked region; thrombopoietin (THPO) and glycoprotein IX (GPIX). Resequencing of four individuals revealed five single-nucleotide polymorphisms (SNPs) in THPO and one mutation in the transmembrane region of GPIX. Analysis of variance showed that the GPIX mutation and one THPO SNP accounted for 6 and 4% of the variation in platelet count, respectively. The THPO SNP lies in the 3' untranslated region of the gene and has not been previously reported. The G to A transition at nucleotide 653 resulted in an Ala 156 (GCC) to Thr (ACC) replacement in the GPIX protein. The GPIX mutation was recently identified in a Chinese patient with Bernard-Soulier syndrome (BSS), a rare recessive bleeding disorder characterized by thrombocytopenia and giant platelets. One copy of the GPIX mutation was found in 300 European individuals with platelet counts within the normal range. The results suggest that two QTLs on chromosome 3q influence platelet count variation in the Asian Indian kindred, with the GPIX transmembrane mutation and the 3' UTR SNP in THPO being strong candidates.