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.

Prehospital hydroxocobalamin for inhalation injury and cyanide toxicity in the United States - analysis of a database and survey of ems providers

Prehospital use of hydroxocobalamin as an antidote for cyanide toxicity, a serious complication of smoke inhalation, has yet to be universally adopted in the United States though its efficacy and safety have been demonstrated since 2006. The purpose of this study was to characterize practices of prehospital hydroxocobalamin administration via a survey of emergency medical services (EMS) and to report a case series from an EMS database to track use of hydroxocobalamin. The Fire Smoke Coalition Newsletter emailed a voluntary survey to EMS subscribers regarding hydroxocobalamin use. Survey responses were analyzed in addition to survival data from the Smoke Inhalation Treatment Database (SITD), a publically available, self-reported, online database for EMS regarding smoke inhalation patient outcomes. Analysis was compared to current published data from PubMed. The survey had a 14% response rate (284/2000). Only 38% reported prehospital utilization of a hydrogen cyanide antidote with 46% using hydroxocobalamin. 20% of responders reported a formal ALS protocol was in place for hydroxocobalamin use. For the SITD, 12 of 13 (92%) patients who received hydroxocobalamin for suspected inhalation survived. Other studies found a survival rate of 72% and 42% after administration of hydroxocobalamin for smoke inhalation. Prehospital administration of hydroxocobalamin for cyanide toxicity is uncommon in the United States, as evidenced by this analysis, despite well-documented safety and efficacy. Although a small sample, patients who received prehospital hydroxocobalamin had improved survival. This survival rate is significantly greater than those reported previously.

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.

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.

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β⁰) 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⁻⁶. 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.

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.

Heterogeneity of the egammadeltabeta-thalassaemias: characterization of three novel English deletions

We have characterized three novel epsilon gamma delta beta-thalassaemia deletions in three English families. Two of the deletions, 114 and 439 kb, removed the entire beta-globin gene complex, including a variable number of flanking olfactory receptor (HOR) genes. The 98-kb deletion extended 90-kb upstream of the epsilon gene to 8 kb upstream of the G gamma-gene, leaving the gamma,delta and beta-genes intact. The 439 kb deletion is the largest deletion reported so far to cause epsilon gamma delta beta-thalassaemia; heterozygotes for this deletion were variably affected by neonatal haemolytic anaemia. Two of the deletions were de novo. Breakpoints of all three deletions occurred within regions of L1 or Alu repeats and contained short regions of direct homology between the flanking sequences, a feature that is likely to have contributed to the illegitimate recombinations.