High-density SNP genotyping to define beta-globin locus haplotypes

Evolutionary history of sickle-cell mutation: implications for global genetic medicine

Resistance afforded by the sickle-cell trait against severe malaria has led to high frequencies of the sickle-cell mutation [HBB; c.20T>A, p.Glu6Val; OMIM: 141900 (HBB-βS)] in most parts of Africa. High-coverage sequencing and genotype data have now confirmed the single African origin of the sickle-cell gene variant [HBB; c.20T>A, p.Glu6Val; OMIM: 141900 (HBB-βS)]. Nevertheless, the classical HBB-like genes cluster haplotypes remain a rich source of HBB-βS evolutionary information. The overlapping distribution of HBB-βS and other disease-associated variants means that their evolutionary genetics must be investigated concurrently. In this review: (1) we explore the evolutionary history of HBB-βS and its implications in understanding human migration within and out of Africa: e.g. HBB haplotypes and recent migration paths of the Bantu expansion, occurrence of ~7% of the Senegal haplotype in Angola reflecting changes in population/SCD dynamics, and existence of all five classical HBB haplotype in Cameroon and Egypt suggesting a much longer presence of HBB-βS in these regions; (2) we discuss the time estimates of the emergence of HBB-βS in Africa and finally, (3) we discuss implications for genetic medicine in understanding complex epistatic interactions between HBB-βS and other gene variants selected under environmental pressure in Africa e.g. variants in HBB, HBA, G6PD, APOL1, APOE, OSBPL10 and RXRA.

Refinement of evolutionary medicine predictions based on clinical evidence for the manifestations of Mendelian diseases

Prediction methods have become an integral part of biomedical and biotechnological research. However, their clinical interpretations are largely based on biochemical or molecular data, but not clinical data. Here, we focus on improving the reliability and clinical applicability of prediction algorithms. We assembled and curated two large non-overlapping large databases of clinical phenotypes. These phenotypes were caused by missense variations in 44 and 63 genes associated with Mendelian diseases. We used these databases to establish and validate the model, allowing us to improve the predictions obtained from EVmutation, SNAP2 and PoPMuSiC 2.1. The predictions of clinical effects suffered from a lack of specificity, which appears to be the common constraint of all recently used prediction methods, although predictions mediated by these methods are associated with nearly absolute sensitivity. We introduced evidence-based tailoring of the default settings of the prediction methods; this tailoring substantially improved the prediction outcomes. Additionally, the comparisons of the clinically observed and theoretical variations led to the identification of large previously unreported pools of variations that were under negative selection during molecular evolution. The evolutionary variation analysis approach described here is the first to enable the highly specific identification of likely disease-causing missense variations that have not yet been associated with any clinical phenotype.

Genetic Modifiers of Fetal Haemoglobin in Sickle Cell Disease

Fetal haemoglobin (HbF) levels have a clinically beneficial effect on sickle cell disease (SCD). Patients with SCD demonstrate extreme variability in HbF levels (1-30%), a large part of which is likely genetically determined. The main genetic modifier loci for HbF persistence, HBS1L-MYB, BCL11A and the β-globin gene cluster in adults also act in SCD patients. Their effects are, however, modified significantly by a disease pathology that includes a drastically shortened erythrocyte lifespan with an enhanced survival of those red blood cells that carry HbF (F cells). We propose a model of how HbF modifier genes and disease pathology interact to shape HbF levels measured in patients. We review current knowledge on the action of these loci in SCD, their genetic architecture, and their putative functional components. At each locus, one strong candidate for a causative, functional DNA change has been proposed: Xmn1-HBG2 at the β-globin cluster, rs1427407 at BCL11A and the 3 bp deletion rs66650371 at HBS1L-MYB. These, however, explain only part of the impact of these loci and additional variants are yet to be identified. Further progress in understanding the genetic control of HbF levels requires that confounding factors inherent in SCD, such as ethnic complexity, the role of F cells and the influence of drugs, are suitably addressed. This will depend on international collaboration and on large, well-characterised patient cohorts with genome-wide single-nucleotide polymorphism or sequence data.

Genetic comparison of sickle cell anaemia cohorts from Brazil and the United States reveals high levels of divergence

Genetic analysis of admixed populations raises special concerns with regard to study design and data processing, particularly to avoid population stratification biases. The point mutation responsible for sickle cell anaemia codes for a variant hemoglobin, sickle hemoglobin or HbS, whose presence drives the pathophysiology of disease. Here we propose to explore ancestry and population structure in a genome-wide study with particular emphasis on chromosome 11 in two SCA admixed cohorts obtained from urban populations of Brazil (Pernambuco and São Paulo) and the United States (Pennsylvania). Ancestry inference showed different proportions of European, African and American backgrounds in the composition of our samples. Brazilians were more admixed, had a lower African background (43% vs. 78% on the genomic level and 44% vs. 76% on chromosome 11) and presented a signature of positive selection and Iberian introgression in the HbS region, driving a high differentiation of this locus between the two cohorts. The genetic structures of the SCA cohorts from Brazil and US differ considerably on the genome-wide, chromosome 11 and HbS mutation locus levels.

Haplotype Analysis of β-Thalassaemia Major and Carriers with Filipino β°-Deletion in Sabah, Malaysia

Objective

The Filipino β°-deletion has been reported as a unique mutation in East Malaysia with a severe phenotype due to the complete absence of β-globin chain synthesis. In this study, the haplotype patterns of the β-globin gene cluster were used to relate the human genetic variation to this specific β-thalassaemia mutation.

Methods

The 376 study subjects included 219 β-thalassaemia major (β-TM) patients with homozygous Filipino β°-deletion and 157 carriers with heterozygous Filipino β°-deletion from 10 government hospitals in different regions of Sabah. Genomic DNA was isolated from whole blood using silica membrane based DNA purification protocol. Polymerase chain reaction restriction fragment length polymorphism analysis (PCR-RFLP) was conducted on five markers within the β-globin gene cluster to construct the haplotype patterns.

Results

Four haplotypes (Haplotype I–IV) were identified with Haplotype I as the predominant haplotype with the highest frequency of 0.98, followed by Haplotype II, III and Haplotype IV with 0.02. Haplotype I was strongly linked with the Filipino β°-deletion among the indigenous population.

Conclusion

Haplotype I as the predominant haplotype suggests the patients with the Filipino β°-deletion in Sabah have a similar origin.

Haplotype Analysis of Three Common β-Thalassemia Mutations in Syrian Patients

β-Globin haplotypes were used to investigate the origin of three common β-globin mutations, IVS-I-110 (G>A); HBB: c.93-21G>A, IVS-I-1 (G>A); HBB: c.92 + 1G>A and codon 39 (C>T); HBB: c.118C > T in Syrian patients. Haplotype analysis was done for 49 unrelated patients with β-thalassemia (β-thal) and 20 unrelated healthy subjects by polymerase chain reaction (PCR)-based restriction fragment length polymorphism (RFLP) for the β-globin gene cluster of the following polymorphic restriction sites: HincII 5' to ε, HindIII 5' to ^Gγ, HindIII 5' to ^Aγ, HincII in ψβ, HincII 3' to ψβ, AvaII in β, and HinfI 3' to β. The IVS-I-110 mutation was associated with three haplotypes: I [+ - - - - + +] (79.4%), V [+ - - - - + -] (5.9%) and VII [+ - - - - - +] (14.7%), while, the two mutations IVS-I-1 and codon 39 were be linked to a single haplotype V (100.0%) and II [- + + - + + +] (100.0%), respectively. The normal chromosomes (β^A/β^A) were associated with four haplotypes, I (50.0%), II (7.5%), V (32.5%) and VII (10.0%). In the Syrian population, the IVS-I-110 mutation was associated with multi haplotypes, whereas the IVS-I-1 and codon 39 mutations have a single origin. More studies for the other mutations will be very useful for genetic epidemiological studies in Syria.

Presencia de haplotipos no africanos incrementa la diversidad genética en pacientes con anemia falciforme en Colombia

El objetivo de este estudio fue identificar la frecuencia de haplotipos dentro del cluster de Beta globina presente en pacientes con anemia falciforme en Colombia, establecer la presencia de haplotipos no africanos en esta población, así como verificar variaciones en el patrón de desequilibrio de ligamiento dentro del cluster de Beta globina. Se analizaron 83 individuos con anemia falciforme, los haplotipos se formaron utilizando cinco sitios de restricción dentro del cluster de Beta globina, se estableció la frecuencia de haplotipos, se calculó el grado de desequilibrio de ligamiento entre los sitios de restricción, así como la similitud genética de esta población con otra de afectados en América. Los haplotipos más frecuentes en la población fueron Benín ( 35,1 %) y Bantú (26, 5 %), ambos africanos. Sin embargo, haplotipos presentes en poblaciones indígenas americanas y europeas alcanzaron frecuencias entre el 2 – 10 %, así como haplotipos que no han sido reportados en otras poblaciones. Los sitios de restricción presentaron bajo o nulo desequilibrio de ligamiento entre ellos. Al compararse con otras poblaciones, la población colombiana presentó mayor similitud con la población de Venezuela en donde Benin y Bantú son también predominantes. Nuestros resultados muestran que el mestizaje ha facilitado el paso de la mutación para la anemia falciforme a un contexto genético no africano (amerindio y europeo). Además, el mestizaje también ha alterado el patrón de desequilibrio de ligamiento dentro del cluster de Beta globina generando modificaciones que pueden tener influencia en estudios de asociación dentro de esta población de afectados. 

Original Research: A case-control genome-wide association study identifies genetic modifiers of fetal hemoglobin in sickle cell disease

Sickle cell disease (SCD) is a group of inherited blood disorders that have in common a mutation in the sixth codon of the β-globin (HBB) gene on chromosome 11. However, people with the same genetic mutation display a wide range of clinical phenotypes. Fetal hemoglobin (HbF) expression is an important genetic modifier of SCD complications leading to milder symptoms and improved long-term survival. Therefore, we performed a genome-wide association study (GWAS) using a case-control experimental design in 244 African Americans with SCD to discover genetic factors associated with HbF expression. The case group consisted of subjects with HbF≥8.6% (133 samples) and control group subjects with HbF≤£3.1% (111 samples). Our GWAS results replicated SNPs previously identified in an erythroid-specific enhancer region located in the second intron of the BCL11A gene associated with HbF expression. In addition, we identified SNPs in the SPARC, GJC1, EFTUD2 and JAZF1 genes as novel candidates associated with HbF levels. To gain insights into mechanisms of globin gene regulation in the HBB locus, linkage disequilibrium (LD) and haplotype analyses were conducted. We observed strong LD in the low HbF group in contrast to a loss of LD and greater number of haplotypes in the high HbF group. A search of known HBB locus regulatory elements identified SNPs 5' of δ-globin located in an HbF silencing region. In particular, SNP rs4910736 created a binding site for a known transcription repressor GFi1 which is a candidate protein for further investigation. Another HbF-associated SNP, rs2855122 in the cAMP response element upstream of Gγ-globin, was analyzed for functional relevance. Studies performed with siRNA-mediated CREB binding protein (CBP) knockdown in primary erythroid cells demonstrated γ-globin activation and HbF induction, supporting a repressor role for CBP. This study identifies possible molecular determinants of HbF production.

Beta-globin gene haplotypes among cameroonians and review of the global distribution: is there a case for a single sickle mutation origin in Africa?

Studies of hemoglobin S haplotypes in African subpopulations have potential implications for patient care and our understanding of genetic factors that have shaped the prevalence of sickle cell disease (SCD). We evaluated HBB gene cluster haplotypes in SCD patients from Cameroon, and reviewed the literature for a global distribution. We reviewed medical records to obtain pertinent socio-demographic and clinical features for 610 Cameroonian SCD patients, including hemoglobin electrophoresis and full blood counts. RFLP-PCR was used to determine the HBB gene haplotype on 1082 chromosomes. A systematic review of the current literature was undertaken to catalogue HBB haplotype frequencies in SCD populations around the world. Benin (74%; n = 799) and Cameroon (19%; n = 207) were the most prevalent haplotypes observed among Cameroonian patients. There was no significant association between HBB haplotypes and clinical life events, anthropometric measures, hematological parameters, or fetal hemoglobin (HbF) levels. The literature review of the global haplotype distributions was consistent with known historical migrations of the people of Africa. Previously reported data from Sudan showed a distinctly unusual pattern; all four classical haplotypes were reported, with an exceptionally high proportion of the Senegal, Cameroon, and atypical haplotypes. We did not observe any significant associations between HBB haplotype and SCD disease course in this cohort. Taken together, the data from Cameroon and from the wider literature suggest that a careful reassessment of African HBB haplotypes may shed further light on the evolutionary dynamics of the sickle allele, which could suggest a single origin of the sickle mutation.

Extensive genomic variability of knops blood group polymorphisms is associated with sickle cell disease in Africa

Sickle cell disease (SCD) is a multisystem disorder characterized by chronic hemolytic anemia, vaso-occlusive crises, and marked variability in disease severity. Patients require transfusions to manage disease complications, with complements, directed by complement regulatory genes (CR1) and its polymorphisms, implicated in the development of alloantibodies. We hypothesize that CR1 polymorphisms affect complement regulation and function, leading to adverse outcome in SCD. To this end, we determined the genomic diversity of complement regulatory genes by examining single nucleotide polymorphisms associated with Knops blood group antigens. Genomic DNA samples from 130 SCD cases and 356 control Africans, 331 SCD cases and 497 control African Americans, and 254 Caucasians were obtained and analyzed, utilizing a PCR-RFLP (polymerase chain reaction-restriction fragment length polymorphism) assay. Analyzing for ethnic diversity, we found significant differences in the genotypic and allelic frequencies of Sl1/Sl2 (rs17047661) and McCa/b (rs17047660) polymorphisms between Africans, African Americans, and Caucasians (P < 0.05). The homozygote mutant variants had significantly higher frequencies in Africans and African Americans but were insignificant in Caucasians (80.2% and 59.6% vs 5.9% for Sl1/2; and 36% and 24% vs 1.8% for McCa/b). With SCD, we did not detect any difference among cases and controls either in Africa or in the United States. However, we found significant difference in genotypic (P < 0.0001) and allelic frequencies (P < 0.0001) of Sl1/Sl2 (rs17047661) and McCa/b (rs17047660) polymorphisms between SCD groups from Africa and the United States. There was no difference in haplotype frequencies of these polymorphisms among or between groups. The higher frequency of CR1 homozygote mutant variants in Africa but not United States indicates a potential pathogenic role, possibly associated with complicated disease pathophysiology in the former and potentially protective in the latter. The difference between sickle cell groups suggests potential genetic drift or founder effect imposed on the disease in the United States, but not in Africa, and a possible confirmation of the ancestral susceptibility hypothesis. The lower haplotype frequencies among sickle cell and control populations in the United States may be due to the admixture and the dilution of African genetic ancestry in the African American population.

Frequency and origin of haplotypes associated with the beta-globin gene cluster in individuals with trait and sickle cell anemia in the Atlantic and Pacific coastal regions of Colombia

Sickle cell anemia is a genetic disease with high prevalence in people of African descent. There are five typical haplotypes associated with this disease and the haplotypes associated with the beta-globin gene cluster have been used to establish the origin of African-descendant people in America. In this work, we determined the frequency and the origin of haplotypes associated with hemoglobin S in a sample of individuals with sickle cell anemia (HbSS) and sickle cell hemoglobin trait (HbAS) in coastal regions of Colombia. Blood samples from 71 HbAS and 79 HbSS individuals were obtained. Haplotypes were determined based on the presence of variable restriction sites within the β-globin gene cluster. On the Pacific coast of Colombia the most frequent haplotype was Benin, while on the Atlantic coast Bantu was marginally higher than Benin. Eight atypical haplotypes were observed on both coasts, being more diverse in the Atlantic than in the Pacific region. These results suggest a differential settlement of the coasts, dependent on where slaves were brought from, either from the Gulf of Guinea or from Angola, where the haplotype distributions are similar. Atypical haplotypes probably originated from point mutations that lost or gained a restriction site and/or by recombination events.

Acute chest syndrome is associated with single nucleotide polymorphism-defined beta globin cluster haplotype in children with sickle cell anaemia

Genetic diversity at the human β-globin locus has been implicated as a modifier of sickle cell anaemia (SCA) severity. However, haplotypes defined by restriction fragment length polymorphism sites across the β-globin locus have not been consistently associated with clinical phenotypes. To define the genetic structure at the β-globin locus more thoroughly, we performed high-density single nucleotide polymorphism (SNP) mapping in 820 children who were homozygous for the sickle cell mutation (HbSS). Genotyping results revealed very high linkage disequilibrium across a large region spanning the locus control region and the HBB (β-globin gene) cluster. We identified three predominant haplotypes accounting for 96% of the β(S) -carrying chromosomes in this population that could be distinguished using a minimal set of common SNPs. Consistent with previous studies, fetal haemoglobin level was significantly associated with β(S) -haplotypes. After controlling for covariates, an association was detected between haplotype and rate of hospitalization for acute chest syndrome (ACS) (incidence rate ratio 0·51, 95% confidence interval 0·29-0·89) but not incidence rate of vaso-occlusive pain or presence of silent cerebral infarct (SCI). Our results suggest that these SNP-defined β(S) -haplotypes may be associated with ACS, but not pain or SCI in a study population of children with SCA.

Stroke in sickle cell anemia patients: a need for multidisciplinary approaches

Sickle cell anemia (SCA) is an autosomal recessive disorder, with Mendelian inheritance pattern, caused by a missense mutation in the β-polypeptide chain of the hemoglobin B. SCA preferentially affects populations in countries where malaria was/is present (e.g. Africa, USA, Brazil). Thereby, in USA, the incidence of SCA is relatively high, around 1/500, and the prevalence is about 1/1000. In Brazil, SCA represents a major public health problem with an incidence ranging from 1/2000 to 1/600 depending on the regions. Homozygotic patients present more severe medical conditions and reduced life expectancy than heterozygous individuals who generally are asymptomatic. Eventually, this life-threatening disease displays a complex etiology owing to heterogeneous phenotypes and clinical outcomes, subsequently affecting the management of the patients. One of the most critical complications associated with SCA is stroke, a leading neurologic cause of death and disability. About 24% of SCA patients have a stroke by the age of 45 and 11% by the age of 20. From the general population, twin and familial aggregation studies as well as genome-wide association studies (GWAS), mostly in pediatric populations with ischemic stroke, showed that the risk of stroke has a substantial genetic component. Nevertheless, to fully characterize genomic contributors of stroke and permit reliable personalized medicine, multidisciplinary studies incorporating knowledge from clinical medicine, epidemiology, genetics, and molecular biology, are required. In this manuscript, stroke in SCA patients is extensively reviewed with emphasis to the US and Brazilian populations. Recent advances in genomics analysis of stroke in SCA patients are highlighted.

Evolutionary constraints in the β-globin cluster: the signature of purifying selection at the δ-globin (HBD) locus and its role in developmental gene regulation

Human hemoglobins, the oxygen carriers in the blood, are composed by two α-like and two β-like globin monomers. The β-globin gene cluster located at 11p15.5 comprises one pseudogene and five genes whose expression undergoes two critical switches: the embryonic-to-fetal and fetal-to-adult transition. HBD encodes the δ-globin chain of the minor adult hemoglobin (HbA2), which is assumed to be physiologically irrelevant. Paradoxically, reduced diversity levels have been reported for this gene. In this study, we sought a detailed portrait of the genetic variation within the β-globin cluster in a large human population panel from different geographic backgrounds. We resequenced the coding and noncoding regions of the two adult β-globin genes (HBD and HBB) in European and African populations, and analyzed the data from the β-globin cluster (HBE, HBG2, HBG1, HBBP1, HBD, and HBB) in 1,092 individuals representing 14 populations sequenced as part of the 1000 Genomes Project. Additionally, we assessed the diversity levels in nonhuman primates using chimpanzee sequence data provided by the PanMap Project. Comprehensive analyses, based on classic neutrality tests, empirical and haplotype-based studies, revealed that HBD and its neighbor pseudogene HBBP1 have mainly evolved under purifying selection, suggesting that their roles are essential and nonredundant. Moreover, in the light of recent studies on the chromatin conformation of the β-globin cluster, we present evidence sustaining that the strong functional constraints underlying the decreased contemporary diversity at these two regions were not driven by protein function but instead are likely due to a regulatory role in ontogenic switches of gene expression.