Helen Rooks

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 β 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.

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.

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.

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.

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.

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 (r2 = 1) and in strong LD with rs7197554 (r2 = 0·75) and rs13336641 (r2 = 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·7thalassaemia 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 A2 (HbA2, α2δ2) levels are largely governed by genetic factors (42% of variability). The influence of age (1%) and sex (4%) was small. HbA2 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 HbA2 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.

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 ɛγδβ-thalassaemias: characterization of three novel English deletions

We have characterized three novel ɛγδβ‐thalassaemia deletions in three English families. Two of the deletions, 114 and 439 kb, removed the entire β‐globin gene complex, including a variable number of flanking olfactory receptor (HOR) genes. The 98‐kb deletion extended 90‐kb upstream of the ɛ gene to 8 kb upstream of the Gγ‐gene, leaving the γ,δ and β‐genes intact. The 439 kb deletion is the largest deletion reported so far to cause ɛγδβ‐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.