1 Unité de Recherche Génomes et Milieux, Faculté des Sciences et Techniques, Université de Nouakchott, Nouveau Campus Universitaire, BP 5026, Nouakchott, Mauritanie
2 Department of Molecular Research, Maurilab Medical Analysis and Research Institute, 2434 Nouakchott, Mauritania
3 The African Centre of Excellence for Genomics of Infectious Diseases, Redeemer’s University, 232101 Ede, Nigeria
4 Faculté de Médecine, de Pharmacie et d’Odonto-stomatologie, Université de Nouakchott, Nouveau Campus Universitaire, BP 880 Nouakchott, Mauritanie
Abstract
Background: Sickle cell disease (SCD) is a major heritable genetic
disease in sub-Saharan Africa, including Mauritania. Fetal hemoglobin (HbF) can
affect the pathophysiology, moderate the clinical course, and offer prospects for
curative treatment of SCD. This study aimed to investigate the influence of
single nucleotide polymorphisms (SNPs) in the BCL11A gene on the levels
of HbF and hematological parameters in Mauritanian sickle cell (HbSS)
patients. Methods: Complete blood count was assessed in 565 patients
suspected to have SCD. Polymerase chain reaction (PCR)–restriction fragment
length polymorphism was performed to identify the HbSS, and sequencing
was used for genotyping three SNPs: rs4671393 (A
Keywords
- sickle cell disease
- fetal hemoglobin
- BCL11A gene
- single nucleotide polymorphisms
- Mauritania
Hemoglobinopathies are genetic disorders affecting the structure or production of the protein portion of hemoglobin [1]. About 5% of the global population carries an abnormal hemoglobin variant; 300,000 to 400,000 infants are born with severe hemoglobinopathies worldwide, and 80% of these cases occur in developing countries [2].
Among the most widespread hemoglobinopathies throughout the world, sickle cell
disease (SCD) is caused by an autosomal recessive inheritance of a single
nucleotide base substitution (HBB c.20T
Although it is caused by a single amino acid substitution, Glu6Val, SCD is characterized by a wide range of clinical manifestations [8, 9]. The cause of this clinical heterogeneity still needs to be fully understood. However, a relatively higher production of fetal hemoglobin (HbF) in adulthood and the consequent increased number of erythrocytes with HbF are likely to be the major factors that modulate the disease phenotype observed in SCD patients [9, 10]. Therefore, high HbF levels in adults are among the most common factors that can lead to a considerable decrease in disease severity by mitigating polymerization and precipitation of sickle cell hemoglobin in red blood cells [11].
HbF synthesis is generally reduced to very low levels, i.e., less than 0.6% of
the total hemoglobin in normal adults [12, 13]. However, the HbF level is
generally elevated in patients with SCD, ranging from 1% to 30% of the total
hemoglobin [11]. HbF expression is regulated by elements linked to the
A follow-up study involving six SCD patients recently validated the BCLA11 gene as an effective target for HbF induction therapy [16]. Several single nucleotide polymorphisms (SNPs) influencing the HbF level have been described at this QTL in patients with SCD or thalassemia worldwide [15, 17]. These QTLs represent 15–30% of the variation in HbF in sickle cell patients of African origin [11]. Subsequent studies on large populations of patients with thalassemia or SCD showed the distribution of several SNPs, which are significantly associated with high expression of HbF in subjects from various geographical origins [12, 13, 14].
SCD is a major heritable genetic disease that affects Africa. With a population
consisting of ‘white Moors’ of Arab descent and Black Africans, Mauritania is
also of concern. An earlier study conducted on hemoglobinopathies among 700 blood
donors in Nouakchott, the capital city of Mauritania, showed that 16.6% of
healthy blood donors were carriers (i.e., heterozygous) of abnormal hemoglobin
(21.4% of Mauritanians of Black African origin and 4.5% of white Moors) [18].
Sickle cell trait (HbA/S) was the most common hemoglobin disorder (8.7%),
followed by hemoglobin C (HbA/C) trait (3%) and
In the present study, three SNPs, rs4671393 (A
The study was conducted from June 2020 to December 2022 at the Maurilab Medical Analysis and Research Institute (MMARI) in Nouakchott. A total of 565 unrelated subjects suspected to be affected by hemoglobinopathy and referred to the MMARI by the Mauritanian SCD association were enrolled. The Mauritanian SCD association supports people with sickle cell disorder in improving their overall quality of life, including free diagnostic testing to detect the presence of abnormal hemoglobin in suspected patients with clinical symptoms suggestive of severe anemia or a family history of anemia. The basic sociodemographic information, including age and sex, was obtained from each participant and recorded in a questionnaire.
The associations between the presence of genetic polymorphisms, HbF levels, and hematologic parameters were analyzed in patients with HbSS aged at least 3 years and who had not received hydroxyurea treatment and red blood cell transfusions during the last three months before inclusion.
Venous blood samples were collected from each patient in K2-ethylenediaminetetraacetic acid (K2-EDTA) tubes (BD Vacutainer, Franklin Lakes, NJ, USA) and immediately analyzed. A complete blood count was obtained by a CELL-DYN Ruby automated Hematology analyzer (Abbott Laboratories, Diagnostics Division, Santa Clara, CA, USA). The blood parameters that were assessed included red blood cell count (RBC), hemoglobin concentration (Hb), sickle hemoglobin (HbS), mean cell (corpuscular) volume (MCV), mean cell hemoglobin concentration (MCHC), fetal hemoglobin (HbF), and hemoglobin A2 (HbA2). Capilarys™ 2 Flex piercing systems (Sebia electrophoresis, Lisses, France) were used to screen for hemoglobin fractions.
Only patients with homozygous sickle cell disease (HbSS) were included
in the genotyping analysis of this study. DNA samples were extracted from
peripheral blood using a DNeasy blood and tissue kit (Qiagen, Redwood, CA, USA),
according to the manufacturer’s instructions. Exon 1 of the
To genotype the SNPs in the BCL11A locus associated with HbF levels in
homozygous sickle cell patients, three regions containing the SNPs specific to
the BCL11A gene (rs11886868, rs4671393, and rs105252)
were amplified using the primer pairs shown in Table 1. PCR amplification was
carried out at a final volume of 25 µL containing 0.04 U Taq DNA
polymerase (Qiagen, Chatsworth, CA, USA), 1
| SNP | Localization | Primer code | Primer sequence 5 |
PCR fragment size |
| rs4671393 | Intron 2 | IVS-II-B_2F | CCAGGGAGCCCTCATTTTGT | 752 bp |
| IVS-II-B_2R | CCCATGGGGCTTCACTTTCT | |||
| rs11886868 | Intron 2 | IVS-II-B_1F | TGCGCTTCTCAGACCCAAAT | 830 bp |
| IVS-II-B_1R | GGCTCCCTGGACTCAAACTC | |||
| rs1052520 | 3 |
3 |
CACAGGCAGAGTCAAGTGCT | 621 bp |
| 3 |
GGGCTGTTTTGCCCAAAGTT |
SNP, single nucleotide polymorphism; F, forward primer; R, reverse primer; PCR, polymerase chain reaction.
The Pandas (https://pandas.pydata.org/), Scipy (https://scipy.org/), and statsmodels packages (https://www.statsmodels.org/stable/index.html) of Python software version 3 (https://www.python.org) were used to conduct statistical analysis. The relationship between HbF levels and hematological parameters was carried out using the Mann–Whitney U test, and the association between HbF levels and genotypes was established using both the Mann–Whitney U test and linear regression.
In the Mann–Whitney U test analysis, three different cut-off values were
investigated: 10%, 15%, and 20%. However, 15% was retained since the results
of this cut-off for HbF level were quite similar to those obtained from linear
regression, which employed the continuous version of HbF as a dependent variable.
Before building a linear regression model for genotypes as independent variables,
dummy variables were used to convert genotype values into binary variables. To
avoid perfect collinearity, we omitted one dummy variable per genotype (e.g.,
AA/AG in s4671393). Allelic and genotype frequencies were calculated and tested
for Hardy–Weinberg equilibrium (HWE) using Pearson’s chi-square test.
Statistical significance was set at p
Five hundred and sixty-five participants were enrolled in the study. The
participants consisted of 51.3% (290/565) females and 48.7% (275/565) males,
with a ratio of 0.95. Most participants (n = 382; 67.6%) had the normal
HbAA. Further, HbS variants were identified in 112 individuals:
heterozygous HbAS (10.9%, n = 62) and homozygous HbSS (8.8%,
n = 50). Seventy-one individuals (12.2%) were characterized as being affected by
other hemoglobinopathies, including 23 (4.1%) with minor
Based on the mean values, the analysis of the hematological parameters of the
HbSS patients showed normocytic anemia with a low mean (
| Parameters | HbF |
HbF |
p-value | ||
| Mean | Mean | ||||
| Age (year) | 13.85 | 10.51 | 19.6 | 14.11 | 0.26 |
| RBC ( |
2.75 | 0.96 | 3.11 | 0.67 | 0.19 |
| Hb (g/dL) | 7.39 | 2.91 | 8.69 | 1.32 | 0.02* |
| MCV (fL) | 87.1 | 14.1 | 85.92 | 10.35 | 0.98 |
| MCHC (g/dL) | 31.1 | 2 | 32.35 | 2.5 | 0.14 |
| HbF (%) | 8.51 | 3.82 | 21.06 | 3.15 | |
| HbS (%) | 88.12 | 3.75 | 76.3 | 2.98 | |
| HbA2 (%) | 3.3 | 0.59 | 2.64 | 0.51 | 0.0001* |
MCV, mean cellular (corpuscular) volume; Hb, hemoglobin; HbF, fetal hemoglobin; HbA2, hemoglobin A2; HbS, sickle hemoglobin; RBC, red blood cell; MCHC, mean cellular hemoglobin concentration; SD, standard deviation. *Statistically significant.
Blood samples from the 50 individuals with the HbSS genotype were
further genotyped for molecular analysis. The sequencing results of the
rs4671393 (A
| SNP | Genotype | Frequency | Allele | Frequency | p-value* | |
| rs4671393 | AA | 0.14 | A | 0.37 | 4.07 | 0.065 |
| AG | 0.47 | G | 0.63 | |||
| GG | 0.39 | |||||
| rs11886868 | CC | 0.17 | C | 0.39 | 11.35 | 0.001 |
| CT | 0.49 | T | 0.61 | |||
| TT | 0.34 |
All 50 samples from HbSS individuals displayed a wild-type genotype in
the rs1052520 allele. SNP, single nucleotide polymorphism; HWE,
Hardy–Weinberg equilibrium. *Significant at p
The distribution of genotypic and allelic frequencies did not deviate from the
HWE (p
Furthermore, relatively high frequencies of minor alleles A (0.37) and C (0.39) for rs4671393 and rs11886868, respectively, were observed among SCD patients.
Statistical analysis showed a significant correlation between rs4671393
and the level of HbF (Mann–Whitney test, p = 0.034) (Table 4). The
presence of the minor allele A in rs4671393 was associated with
an elevated level of HbF (mean 12.72
| Locus/chromosome location | SNP | Genotype | N | HbF (%) | Median | p-value | |
| Mean | |||||||
| BCL11A/Chr2 | rs4671393 A |
AA+AG | 27 | 12.72 | 6.26 | 12.1 | 0.034* |
| GG | 23 | 9.03 | 5.77 | 7.5 | |||
| rs11886868 C |
CC+CT | 26 | 11.83 | 6.42 | 10.25 | 0.25 | |
| TT | 21 | 10.01 | 6.41 | 8.9 | |||
Chr2, chromosome 2; SNP, single nucleotide polymorphism. *Statistically significant.
In SCD, a single
HbF is known as the principal SCD genetic modifier, and the presence of
This is the first study focusing on the genotyping of three SNPs (rs118868686, rs4671393, and rs1052520) in the BCL11A gene and their correlation with HbF levels in Mauritanian subjects with SCD. In a previous study that genotyped 95 SNPs, including 17 at the BCL11A locus, in 1032 African–Americans with sickle cell anemia, Galarneau et al. [31] found that rs4671393 was the genetic marker that was most strongly associated with HbF levels. In that study, they also found two other SNPs in intron 2 of BCL11A, rs7599488 and rs10189857, which were associated with HbF levels. In the present study, these latter SNPs were not included in our molecular analysis.
Our results showed the presence of three different genotypes for rs4671393 and rs11886868, with a slightly higher prevalence of minor alleles than those reported in Cameroonians [32]. This slight difference can be explained by the ethnic diversity of our study population, which includes a high proportion of non-Black African populations and the relatively small number of SCD patients included and analyzed. Regarding rs1052520, a single homozygous genotype (GG) was observed in all individuals with HbSS (n = 50), which suggests that this SNP may be fixed and may, therefore, not be suitable as a genetic marker in Mauritania.
In 2007, Menzel et al. [13] reported that intron 2 of the BCL11A gene has the strongest associations with HbF variation in a 14 kb sequence. In addition, genome-wide association studies identified the rs4671393 SNP as the determinant marker of the HbF level [32]. This SNP alone explained 10% of the variation in HbF levels in patients with sickle cell anemia from northern Brazil [33]. In addition, the minor A allele was also associated with significantly higher HbF levels [32]. It is worth noting that our results revealed a statistically significant correlation between rs4671393 and HbF level in individuals with HbSS. Indeed, the presence of the minor (A) allele was associated with the HbF level increase in our study population. These results are similar to those found in an Angolan pediatric population [34].
The main limitations of the present study are the inclusion of a small number of individuals with SCD and a limited number of SNPs analyzed in the BCL11A locus or reported to modulate HbF levels.
The results of the present study highlight the differential frequencies of three SNPs in the BCL11A locus in Mauritanian individuals with HbSS that were previously reported to be associated with HbF levels. The HbF level is crucial in determining the severity of SCD, and genetic modification of HbF expression may offer clinical benefits in diagnosis and disease monitoring. The findings of the present study have important implications for genetic counseling and clinical management of SCD patients.
However, further research is needed to explore the association between genetic markers and HbF levels to improve the characterization of sickle cell disease and potentially improve the management of this life-threatening condition.
The authors confirm that the data supporting the findings of this study are available within the article and its supplementary materials.
ATB collected data, performed blood and molecular analyses, interpreted the data, wrote the first draft of the manuscript. MT collected data and performed blood analysis. HS performed statistical analysis. SMG performed statistical analysis, interpreted data and supervised laboratory analyses. AM and AOMSB conceived, designed, and supervised the research, wrote, critically reviewed and edited the manuscript. All authors have participated sufficiently in the work to take public responsibility for appropriate portions of the content and agreed to be accountable for all aspects of the work in ensuring that questions related to its accuracy or integrity. All authors read and approved the final manuscript. All authors contributed to editorial changes in the manuscript.
The Institutional review board of the University of Nouakchott reviewed and
approved the study protocol (approval no. 00000062/19/UNA). The study was
conducted in accordance with the declaration of Helsinki and approved by the
ethics Committee of the University of Nouakchott. The objective of the study was
explained to all participants in their respective local mother dialect,
Hassaniya, Pular, Soninke or Wolof. Written informed consent was obtained from
each adult participant or, in the case of children
Not applicable.
This research received no external funding.
The authors declare no conflict of interest.
References
Publisher’s Note: IMR Press stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
