Skip to main content
PeerJ logoLink to PeerJ
. 2022 Jul 5;10:e13487. doi: 10.7717/peerj.13487

G6PD and HBB polymorphisms in the Senegalese population: prevalence, correlation with clinical malaria

Fatou Thiam 1,✉,#, Gora Diop 2,3,#, Cedric Coulonges 4, Céline Derbois 5, Babacar Mbengue 6, Alassane Thiam 3, Cheikh Momar Nguer 1, Jean Francois Zagury 4, Jean-Francois Deleuze 5, Alioune Dieye 6
Editor: Katherine Mitsouras
PMCID: PMC9266585  PMID: 35811813

Abstract

Background

Host genetic factors contribute to the variability of malaria phenotypes and can allow a better understanding of mechanisms involved in susceptibility and/or resistance to Plasmodium falciparum infection outcomes. Several genetic polymorphisms were reported to be prevalent among populations living in tropical malaria-endemic regions and induce protection against malaria. The present study aims to investigate the prevalence of HBB (chr11) and G6PD (chrX) deficiencies polymorphisms among Senegalese populations and their associations with the risk for severe Plasmodium falciparum malaria occurrence.

Methods

We performed a retrospective study with 437 samples, 323 patients recruited in hospitals located in three different endemic areas where malaria episodes were confirmed and 114 free malaria controls. The patients enrolled were classified into two groups: severe malaria (SM) (153 patients) and uncomplicated malaria (UM) (170 patients). PCR and DNA sequencing assessed host genetic polymorphisms in HBB and G6PD. Using a multivariate regression and additive model, estimates of the impact of human HBB and G6PD polymorphisms on malaria incidence were performed.

Results

Six frequent SNPs with minor allele frequencies (MAF) > 3% were detected in the HBB gene (rs7946748, rs7480526, rs10768683, rs35209591, HbS (rs334) and rs713040) and two in the G6PD gene (rs762515 and rs1050828 (G6PD-202 G > A). Analysis of selected HbS polymorphism showed significant association with protective effect against severe malaria with a significant p-value = 0.033 (OR 0.38, 95% CI [0.16–0.91]) for SM vs. UM comparison. Surprisingly, our study did not identify the protective effect of variant HbC polymorphism against severe malaria. Finally, we found some of the polymorphisms, like HbS (rs334), are associated with age and biological parameters like eosinophils, basophils, lymphocytes etc.

Conclusion

Our data report HBB and G6PD polymorphisms in the Senegalese population and their correlation with severe/mild malaria and outcome. The G6PD and HBB deficiencies are widespread in West Africa endemic malaria regions such as The Gambia, Mali, and Burkina Faso. The study shows the critical role of genetic factors in malaria outcomes. Indeed, genetic markers could be good tools for malaria endemicity prognosis.

Keywords: HBB, G6PD, Polymorphisms, Severe malaria, Senegal

Background

Malaria is caused by Plasmodium species infection and affects hundreds of millions of people per year. In 2020, an estimated 241 million malaria cases occurred in 85 malaria-endemic countries, increasing from 227 million in 2019, with most of this increase coming from countries in the African Region (95% of cases). Malaria disease remains a significant cause of death and is still the 4th leading cause of death from infectious diseases worldwide (World Health Organization, 2021; Papaioannou, Utzinger & Vounatsou, 2019; Moxon et al., 2020). The symptoms caused by Plasmodium falciparum infection are highly variable, ranging from asymptomatic (AM) and uncomplicated (UM) to severe clinical forms (SM) (Conway, 2007; Miller et al., 2002; Niang et al., 2016). Reasons for these different progressions from asymptomatic, uncomplicated to severe forms or fatal outcomes are not fully understood (Milner, 2018). Understanding these mechanisms is necessary for global malaria eradication (Mackinnon et al., 2005; Weatherall & Clegg, 2002; Griffiths et al., 2005).

Host genetic factors contribute to the variability of malaria phenotypes (Moxon et al., 2020) and may provide insight into the mechanisms involved in susceptibility and/or resistance to Plasmodium species infection outcome (Miller et al., 2002). Several genetic polymorphisms are prevalent among populations living in tropical malaria-endemic regions and induce protection against malaria (Miller et al., 2002; Niang et al., 2016). Glucose-6-phosphate dehydrogenase (G6PD) and β-chain of haemoglobin (HBB) deficiencies are more prevalent in malaria-endemic countries and are two of the essential loci conferring resistance to severe malaria in humans (Mackinnon et al., 2005; Grignard et al., 2019).

Mutations in the HBB gene, located on Chromosome 11 (Onda et al., 2005), are responsible for several serious hemoglobinopathies, such as sickle cell anaemia and β-thalassemia (Lopera-Mesa et al., 2015; Modiano et al., 2001; Taylor, Parobek & Fairhurst, 2012; Williams, 2005a; Williams, 2015; Williams, 2005b). To date, there are more than 900 variants in HBB (Aldakeel et al., 2020), and in malaria-endemic areas, several HBB variants are known to have a protective effect against P. falciparum malaria (Grant et al., 2015). HbAS variant has been associated with 50% and 80% reduced risks of developing uncomplicated and severe malaria (Williams, 2005b), HbAC and HbC variants have been associated with 30% and 93% reduced risks of developing mild and severe malaria, respectively (Modiano et al., 2001; Taylor, Parobek & Fairhurst, 2012; Agarwal et al., 2000; Travassos et al., 2015). HbC is prevalent only in West Africa  (Tetard et al., 2017; Piel et al., 2013a), and allele frequencies above 15% have been described in West African populations (Piel et al., 2013b). The sickle hemoglobin  HbS is encountered at frequencies of up to 18% across sub-Saharan Africa, the Middle East and South Asia. HbSS homozygous subjects develop sickle-cell anaemia disease associated with low life expectancy. Generally, subjects succumb in early life; unlike the HbCC homozygosity, the HbAS and HbAC heterozygous traits are clinically benign (Piel et al., 2013a).

G6PD is encoded by a 16.2 kb gene found on the X chromosome  (Gomez-Manzo et al., 2016). Approximately 160 genetic variants causing clinical deficiency of G6PD have been characterised (Mason, Bautista & Gilsanz, 2007). G6PD deficiency affects over 400 million people living in tropical and subtropical countries, 15–30% of whom are founded in sub-Saharan Africa  (Lwanira et al., 2017; Carter et al., 2011; Clarke et al., 2017). A high diversity of G6PD variants has been reported, including the most common form G6PD B (wild type enzyme), characterised by the 202G-376A allele. In sub-Saharan Africa, up to 40% of the population carries the G6PD A+ form (non-deficient type) characterised by 202G-376G allele with 85% enzyme activity (Battistuzzi et al., 1977). With the presence of an added mutation at the 202 nucleotides (G > A), the double mutant (202A-376G) is named G6PD A- with only 12% enzyme activity compared to the B allele (Beutler, 1989). The 202A- allele has been associated with protection against severe malaria in African populations in Mali, the Gambia, and Uganda (Guindo et al., 2007; Ruwende et al., 1995; Walakira et al., 2017). However, this association was less significant than sickle haemoglobin and α-thalassemia (Verra, Mangano & Modiano, 2009).

The high prevalence of malaria protective polymorphisms is clearly associated with endemicity, but the prevalence varies among populations living in endemic areas. Local malaria risk can explain some of these differences, as suggested by the decreased prevalence of protective polymorphisms with increasing altitude in endemic countries (Hedrick, 2011). The prevalence of some protective polymorphisms has also been shown to vary between ethnic groups (Hedrick, 2011).

In Senegal, the polymorphisms of HBB and G6PD have already been investigated. Main studies focus on a small number of participants and/or ethnic groups. Since 1960, some data have been obtained concerning the prevalence of HbS and G6PD polymorphisms in Bedik and Niokholonko populations in the Southeastern part (Mauran-Sendrail et al., 1975). The study of G6PD from the Niokholonko ethnic group reported the frequencies of two variants of G6PD: 0.18 for G6PD A+ and 0.07 for G6PD A- (Vergnes, Gherardi & Bouloux, 1975). It also appears that in the Sereer ethnic group, from the Niakhar area located 150 km southeast of Dakar, Senegal, the allele G6PD-376G/202A is not the common variant but rather the 376G/968C variant (De Araujo et al., 2006). A study performed in 2005 in the Dakar region focused on 300 individuals described the prevalence and morbidity of G6PD Deficiency in the Sickle Cell Disease homozygote group (Diop et al., 2005). In neighbouring countries such as the Gambia, it was reported that the minor allele frequency of G6PD-202A in healthy participants was 0.03, considerably lower than those reported in eastern African countries such as Kenya (0.18) and Malawi (0.19) (Clarke et al., 2017).

However, until now, there has been a lack of data about the frequencies of the general Senegalese population, from which full studies have not yet been performed. It is essential to perform additional screening of HBB and G6PD loci to identify potential variants associated with susceptibility to severe malaria and corresponding hematologic parameters.

Senegal is a West African country where malaria risk varies from one region to another. Our present study has two aims: (i) first, to investigate the prevalence of key malaria-protective polymorphisms in G6PD and HBB genes in the Senegalese population, (ii) secondly, we will check the correlation between main polymorphisms in HBB and G6PD genes and severe malaria outcome. In participants from three ecological zones, we characterised polymorphisms on HBB and G6PD genes where the malaria endemicity was different. More data on G6PD and HBB polymorphisms and their incidence on malaria outcome and corresponding haematological parameters are essential for different malaria eradication strategies in endemic West African countries.

Material and Methods

Senegalese population, ethnical distribution and study sites

The cohort included black Senegalese-born individuals recruited between 2003 and 2015, whose parents and grandparents were born in Senegal, a malaria-endemic country in the Sahelian zone of West Africa.

The Senegalese population (SP) recorded in 2013 is 13.508.715 inhabitants, including near equality between men 49.9% and women 50.1% (Agence Nationale de la Statistique et de la Démographie (ANSD), 2014). SP is unevenly distributed in space and is concentrated in the West and Center of the Country, while the East and the North are sparsely populated. Then the highest density of population is found in cities such as Dakar, the Capital, Thies, Touba, Saint-Louis-Louga, and Kaolack.

Senegalese ethnic groups are located almost everywhere in Senegal, but the ‘Wolof Ethnic group’ is concentrated in the West and West-Center; the ‘Serere populations’ are focused on Peanut-Basin, specifically in Kaolack and Fatick regions. The Peulh ethnic groups are scattered throughout Senegal areas but are more established in Ferlo along the river and the Kolda region. The ‘Sarakholes’ live in Bakel. Other ethnic groups, such as Diolas, Baïnouks, Balantes and Madingues Ethnic groups, are located in Casamance. Bedik and Niokholonko populations are located in the country’s Southeastern part, such as Tambacounda and Kedougou regions.

Clinical Malaria participants were enrolled from five “site-regions” in Senegalese regional hospitals in St Louis and Louga in the North, Dakar, Tambacounda and Kolda (Fig. 1). These targeted urban cities polarise different ethnic groups, particularly in Dakar, the capital of Senegal, even though there is a difference in the ethnic distribution between North to South.

Figure 1. The study sites.

Figure 1

Study region map. Black points mark the towns surveyed in Senegal and show the study sites. The yellow color areas keep the different ecological areas of endemicity, in north the preelimination area: St Louis-Louga regions, the capital Dakar-Diamniadio center and the south endemic region Tambacounda-Kolda.

For convenience, these five regions were divided into three ecological areas according to the climatic gradient and malaria endemicity, translated from North to South. The northern part is Sahelian climates with unstable or seasonal malaria, including St Louis and Louga regions. These sites belong to the ecological area where malaria is hypoendemic and where the annual incidence in 2019 was less than five cases per 1,000 inhabitants. Some sites are under 1/1000 with very low anopheles vector density. This region is located in the malaria prelimination area context (Thiam et al., 2021; Thiam et al., 2020). The ‘centre’ part is Dakar, an urban centre with low malaria prevalence rates and better health coverage due to better healthcare access. The southern region includes Kolda, and Tambacounda and corresponds to the Sudano-Guinean zone, with annual rainfall exceeding 800 mm (Sahondra Harisoa et al., 2001). Until now, malaria transmission has been hyperendemic in southern regions, and the high transmission season in Senegal occurs mainly between July and October (Jambou et al., 2001; Diouf et al., 2017).

Sampling, clinical and controls participants

Clinical malaria participants corresponded to participants with Plasmodium-positive QBC (Quantitative Buffy Coat). The QBC is a qualitative test, conceptualised as far back as 1974, and is used to identify the malarial parasites in the peripheral blood faster than the conventional thick smear method. The technique involves staining the centrifuged and compressed red cell layer with acridine orange, which is later examined under a fluorescence microscope (Benito et al., 1994; Prashanth, 2012).

According to the criteria defined by Saissy, Rouvin & Koulmann (2003) the malaria participants were assigned to two groups: uncomplicated malaria (UM) and severe malaria (SM). ‘Severe malaria’ phenotype is characterised by the cerebral form, an organ failure and/or metabolic dysfunctions secondary to P. falciparum infection. High parasitemia, Severe anaemia and cerebral form are the three most ‘severe forms’ currently explored. But other complications such as hypoglycaemia, thrombocytopenia, renal insufficiency, hepatic or even pulmonary oedema may appear alone or in combination. For uncomplicated forms, we considered patients who had fever with P. falciparum parasitaemia of <25,000 parasites/µL of blood with no evidence of impaired consciousness or seizures at the time of enrollment and no other past background of mental illness, meningitis or accidental head injury were included.

To ensure homogeneity of data, inclusion criteria were: (1) Only black Senegalese individuals with P. falciparum infection confirmed in diagnostic; (2) persons born in Senegal; whose parents and “grandparents” were born in Senegal. Exclusion criteria were: (1) other racial or ethnic groups living in Senegal; (2) participants with clinical signs of severity or any state that may interfere with the study, such as a recent pregnancy and childbirth. The control group (CTR) corresponded to unaffected participants, free of Plasmodium infection. These participants were recruited at Pasteur institute at the biomedical laboratory according to the same criteria used for clinical participants inclusion, and live in the same environmental and endemic area”.

A blood sample (10 mL) was taken from venous blood in heparin-containing vacutainer tubes for each patient. Hemogram and Several blood cell indices, such as red blood cells (RBCs) count, Hematocrit, MCV (mean corpuscular volume) and MCHC (mean corpuscular volume concentration), were taken in the different groups and evaluated according to references normal values suggested elsewhere (Lopera-Mesa et al., 2015; Modiano et al., 2001). In SM and UM participants, parasitaemia h been estimated by first counting the number of parasites per 200 white blood cells in a thick blood film and then calculating the parasite count/µL from the total white blood cell count µL (Ngom et al., 2007).

DNA extraction, PCR and sequencing

As previously described by Diop et al. (2018), genomic DNA was extracted from peripheral blood using a standard Qiagen Kit following the manufacturer’s recommendations (QIAmp kit Cat. No 51306).

Oligonucleotide primers were designed to amplify interest regions in HBB and G6PD by PCR. At all, PCR reactions were performed using a Gotaq®Green Master Mix (Promega, Germany) in a total volume of 25 µl containing 25 ng of genomic DNA (5 ng/µl) and 2.5 µl of each primer (10 µM). PCR conditions were initial denaturation at 95 °C for 5 min, 35 cycles at 95 °C for 30 s, 62 °C for 45 s, and 72 °C for 1 min, with a final extension at 72 °C for 10 min.

PCR amplification for the Hemoglobin HBB gene was performed using a reference sequence located at chromosome 11 (GRCh38:11:5226434-5227234). The sequence region was downloaded from the BLAST tools, and the primers positions of HbS—rs334 ([A/T]), HbC—rs33930165 ([C/T]) polymorphisms were designed. The reference sequence downloaded is orientated forward concerning the reference genome. The primers sequences designed for HBB amplicons were reversed concerning the reference, such as forward (5′-ACTCCTAAGCCAGTGCCAGA-3′) and reverse (5′-CGATCCTGAGACTTCCACAC-3′) primers. An amplicon of 801 bases was obtained, targeting mutations of deficiencies, including HbS (rs334) and HbC (rs33930165) and surrounding polymorphisms.

For the G6PD gene, the location of the amplicon was (GRCh38: X:154535674-154536374). PCR amplification was performed with forward (5′-GTCTTCTGGGTCAGGG AT-3′) and reverse (5′-GGAGAAAGCTCTCTCTCC-3′) primers, amplifying 701 bases long and targeting the G6PD-202 and their surrounding SNPs. Then the designed primers are reverse with respect to the human reference sequence (GRCh38: X:154535674-154536374) (http://www.ensembl.org/index.html). Then the amplicon left out SNP previously identified in the Senegalese population, such as G6PD +376, G6PD +376, G6PD +542, G6PD +680 and G6PD +968.

The PCR products were examined for specificity via 2% agarose gel electrophoresis, and the amplicons were purified using BioGel P100 gels (Bio-Rad). Sequencing reactions (2 µl of PCR product) were performed using the Dye Terminator v3.1 method in an ABI PRISMs 3730 DNA Analyzer (Applied Biosystems, Foster City, CA, USA). Sequencing conditions were: 96 °C for 5 min, 25 cycles of 96 °C for 10 s; 60 °C for 4 min and 15 °C forever, and PCR products were purified with Sephadex G50 superfine columns (GE Healthcare). Alignment and SNP discovery were performed using NG_009015.2 and NG_059281. as reference sequences for G6PD and HBB, respectively. Analysis was performed with Genalys version 2.0 software (Takahashi et al., 2003).

Statistical analysis

Statistical association tests Epi Info software (version 7.0.8.1) were performed to evaluate the association between malaria clinical status and G6PD and HBB polymorphisms. At first Allelic frequencies and Hardy-Weinberg equilibrium were calculated, as described (Rodriguez, Gaunt & Day, 2009). The differences in allelic frequencies between the three groups (SM, UM, CTR) were determined using the logistic regression analysis method, as reported previously (Barrett et al., 2005; Tregouet & Garelle, 2007).

Associations between SNPs and different parameters such as age, biological and hematological parameters were performed using the Mann–Whitney test, and then associations with p values < 0.05 were considered statistically significant.

Different approaches were explored for association analysis with G6PD loci polymorphism to assess hemizygosity. First, a random inactivation was coded accordingly to genotype as [0, 2] for males and [0, 1, 2] for females, accounting for the random inactivation and assuming an activation of one of the two X chromosomes. Second, an ‘SNP random inactivation escape’ was taken into account, and a test was performed in combined samples, given that a proportion of SNPs can escape from random inactivation used at first. Then gender was treated as a covariate, and the association tests were adjusted in the discovery sample. Third, the association test for males and females was performed separately to investigate sex-specific effects underlying X chromosome loci. An association test accounting for gender × SNP interaction was also performed. Meta-analyses were then conducted using the R package (version 1.9–8) for combined discovery analyses. For this model approach separating males and females would cause a lost statistical power according to the small size of our sampling. Finally, the heterogeneity among these datasets from the different models was evaluated by tests such as I2 and Q test > 50% or P < 0.05 were considered heterogeneous, and the random-effect model was applied; otherwise, the fixed-effect model would be used. Linkage disequilibrium (LD) was computed for each pair of polymorphisms within the HBB and G6PD genes using Haploview software (Barrett et al., 2005). The Lewontin’s D′ coefficient correlates to the level of recombination, which helps find a haploblock. The LD was calculated as a summary and according to all ethnic group populations. Haplotype estimates were obtained using Shapeit4 (Delaneau et al., 2019). Nominal p-values were corrected under the cofounders association effects HbS and HbC polymorphisms.

Ethics

The study’s objectives have been explained clearly using local dialect before performing inclusion of patients in hospitals centres. The protocol has been reviewed according to the rules issued by the National Committee for Ethics for Health Research (CNERS) of Senegal and according to the procedures established by the Cheikh Anta Diop University of Dakar (UCAD) for the ethical approval of any research involving human participants. Written informed consent was obtained from adult participants and parents or legal representatives of children. In addition, based on the information provided, UCAD’s Committee on Research and Ethics (CER) considers that the research proposed respects the appropriate ethical standard and, as a result, approves its execution under “Protocole 0344/2018/CER-UCAD”. All patients enrolled in the cohort/or legal representative gave signed and informal written consent to provide a blood sample for further studies.

Results

Study sites and populations

We performed a retrospective study on 437 participants recruited between 2003 and 2015 in three ecological areas distinguished by their malaria endemicity (Fig. 1, Table 1). Of those, 27 (6%) were from St Louis-Louga, 258 (59%) from Dakar, and 147 (33%) from Tambacounda-Kolda. In total, 153 participants had fulfilled WHO criteria for SM (Bei et al., 2017), 170 for UM, and 117 for CTR. Population CTR controls were recruited in the same ecological and endemic area as clinical cases, are free of ‘Plasmodium’ infection and represent 26.7% of the study population. The median age was 25.8 years in individuals with SM and 18.9 years in individuals with UM. Eighteen percent (18%) were children under five years of age, 15% were 5-15 years old, and 55% were above 15 years old. In our cohort, no mortality was observed in UM and CTR groups. However, 13.7% of mortality was observed in severe malaria patients, with 1.5% in children 15 years and under. 207 cases (47.36%) were male in our population, and 194 (44.39%) were female.

Table 1. General Characteristics of the study participants.

The number of patients in each group SM, UM and CTR, is shown. Age is given with median values (the number of subjects and corresponding percentages is given in parenthesis). ND the number for which, data are undetermined and/or unknown.

Severe Malaria (SM) Uncomplicated Malaria (UM) Control (CTR) Total
(N = 153) (N = 170) (N = 114) (N = 437)
Age (years) <5, n (%) 29 (18.9) 45 (26.4) 6(5.2) 80 (18.3)
[5–15], n (%) 29 (18.9) 28 (16.4) 20 (17.5) 77 (17.6)
[15–25], n (%) 37 (24.2) 15 (8.8) 6 (5.2) 58 (13.3)
>25, n (%) 52 (33.9) 48 (28.2) 76 (66.6) 176 (40.3)
ND, n (%) 6 (3.9) 31 (20) 9 (7.8) 46 (10.5)
Mean, [interval] 25.8 (25–89) 18.9 (13–77) 35.4 (32–87) 24.9 (23–89)
Survival Survivors, n (%) 117 (76.5) 169 (99.4) 117 (100) 403 (92.2)
Deaths, n (%) 21 (13.7) 0 (0) 0 (0) 21 (4.8)
ND, n (%) 15 (9.8) 1 (0.6) 0 (0) 16 (3.6)
Gender Male, n (%) 86 (56.2) 71 (41.7) 50 (42.7) 207 (47.4)
Female, n (%) 55 (36) 81 (47.6) 58 (49.6) 194 (44.4)
ND, n (%) 12 (7.8) 18 (10.7) 9 (7.7) 39 (8.9)
Locality Dakar, n (%) 82 (53.6) 67 (39.4) 114 (97.4) 263 (60.2)
Tambacounda/Kolda, n (%) 62 (40.4) 85 (50) 0 (0) 147 (33.6)
Saint-Louis/Louga, n (%) 9 (6) 18 (10.6) 0 (0) 27 (6.2)

Haemotological and parasitaemia parameters of the study population

The average values of selected haematological parameters were determined for three SM, UM and CTR groups and the association was calculated using ANOVA statistical tests. We found a difference in the level of parasitaemia between SM and UM groups, which was statistically significant (p < 0.01) with high-level parasitemia in the SM group (mean ± SD = 27220 ± 5458.3 P/µL). The haematological parameters of the three groups were compared, as shown in Table 2. Several blood cell indices, such as RBC (red blood cells) count, Hematocrit, MCV, and MCHC, differed significantly between malaria groups. We noted anaemia in the SM group with haemoglobin levels <10 g/dl of red blood cells and platelet counts were significantly lower (p-value < 0.0012). In contrast, the leucocytes count was significantly increased (p-value = 0.001) in SM compared to UM patients. In addition, we found higher levels of eosinophils, basophils and eosinophils in UM participants than SM participants (p-value < 0.01).

Table 2. Haematological and Parasitaemia parameters of the study participants.

The mean and SD (Standard Deviation) values of selected haematological parameters and parasite density were determined for three SM, UM and CTR groups. The parasitemia and haematological parameters of the three groups were compared using ANOVA test. Values were statistically significant when compared SM to UM with p < 0.001.

Parameters Severe Malaria (SM) Uncomplicated Malaria (UM) Control (CTR) P-value (SM/UM)
Mean ± SD Mean ± SD Mean ± SD
Parasitaemia (P/µL) 27220 ± 5458.3 3993 ± 1327.3 0.00 ± 0.0 <0.001
Leukocytes (×103/µL) 12.9 ± 1.0 8.8 ± 0.5 6.1 ± 0.2 <0.001
Neutrophils (%) 52.02 ± 3.7 60.29 ± 2.5 45.77 ± 1.3 <0.001
Lymphocytes (%) 22.32 ± 2.4 26.21 ± 2.1 40.09 ± 1.3 <0.001
Monocytes (%) 6.27 ± 0.9 7.56 ± 0.9 8.65 ± 0.3 0.009
Eosinophils (%) 0.49 ± 0.1 2.77 ± 0.48 4.33 ± 0.45 0.0098
Basophils (%) 0.44 ± 0.1 0.72 ± 0.1 1.14 ± 0.08 <0.001
Red Blood Cells (×106/µL) 3.71 ± 0.29 4.11 ± 0.07 4.66 ± 0.06 0.0012
Hba (g/dL) 9.49 ± 0.3 12.11 ± 0.2 12.96 ± 0.2 <0.001
Hematocrit (%) 27.88 ± 0.9 36.68 ± 0.7 38.44 ± 0.5 <0.001
MCV (fL) 79.96 ± 0.64 87.70 ± 1.37 82.98 ± 0.64 <0.001
MCHC (pg/cell) 26.48 ± 0.47 29.26 ± 0.47 28.04 ± 0.23 <0.001
Platelets (×103/µL) 124.1 ± 7.6 241.8 ± 5.2 286.8 ± 11.2 <0.001

Notes.

Hb levels and parasite density distribution are given with the mean and SD values (standard deviation). Hb* values statistically significant when compared SM to UM with p < 0.001. Parasite density* p < 0.001 when comparing SM to UM.

G6PD and HBB polymorphisms and structure in the Senegalese population

To explore the prevalence of G6PD and HBB deficiencies in our cohort, the polymorphism of G6PD (Chr.X) and HBB (Chr.11) genes and their flanking sequences were analysed in 437 Senegalese populations, including 153 SM, 170 UM and 117 CTR participants. Tables 3 and 4 summarise the frequency of each polymorphism in the SM, UM, and CTR groups living in three ecological areas. Polymorphisms were obtained both in forward and reverse sequencing, and allelic frequencies obtained in this study are like the data provided by the NCBI dbSNP database about the African population.

Table 3. Frequencies and single nucleotide polymorphism (SNP) of G6PD gene and association analysis with susceptibility to severe malaria (SM).

G6PD-SNP NCBI dbSNP number Phénotype MAF HWE
SM UM CTR Overall Population SM vs UM UM vs CTR SM vs CTR
P-value OR (95% IC) P-value P-value
+10588_A >G rs762515 0.280 0.373 0.325 0.329 0.047 0.65 (0.43–0.99) 0.344 0.415
+10707_G >A newSNP 0 0.004 0 0.001 1 1 1
+10776_G >C newSNP 0 0.008 0 0.003 0.509 0.517 1
+10869_G >A rs1050828 G6PD202A 0.022 0.032 0.018 0.026 0.568 0.536 1
+10983_G >C newSNP 0.005 0.008 0.012 0.008 1 1 0.607
+11000_C >T rs782500951 0.005 0.004 0 0.003 1 1 1

Notes.

The p values for statistical tests were performed using the linear regression model analysis for each polymorphisms. Association analysis G6PD for polymorphisms (Table 3) and HBB polymorphisms (Table 4) were performed separately, and corrections to each other were applied to reflect the real effect of each. Analysis has been carried out by comparing SM vs UM, UM vs CTR, and SM vs CTR. Borderline (0.05 ≤ p ≤ 0.1) and significant (0 −  ≤ p ≤ 0.05) p values are in bold. The OR (odds ratio) and CI (Confidence intervals) were shown when p values were significant. The phenotype attributed by from each polymorphism was mentioned (from NCBI) as: Erythrocytosis 6 familial (E6 familial), Beta-thalassemia (β-thalassemia), Fetal hemoglobin quantitative trait loci 1 (Fetal HQTL-1), Heinz body anaemia (H-B-A), METHEMOGLOBINEMIA (MTH). The SNP labelled as “benign” were ‘non-affected SNP’, then non associated with many diseases such as Thalassemia, HbSS disease, and/or any other pathology with a notorious phenotype.

SM
Severe Malaria
UM
Uncomplicated Malaria
CTR
Control group
MAF
Minor Allele Frequency
HWE
Hardy–Weinberg

Table 4. Frequencies of single nucleotide polymorphism (SNP) of HBB gene and association analysis with susceptibility to severe malaria (SM).

HBB-SNP NCBI dbSNP number Phénotype (NCBI) MAF HWE
SM UM CTR Overall population SM vs UM UM vs CTR SM vs CTR
P-value OR (95% CI) P-value OR (95% CI) P-value
+526_C >T rs7946748 β-thalassemia 0.06 0.03 0.06 0.05 0.06 0.048 0.41 (0.17–0.97) 1
+519_T >G rs7480526 Benign 0.34 0.40 0.33 0.36 0.14 0.14 1
+461_C >G rs10768683 Benign 0.09 0.06 0.09 0.08 0.26 0.14 0.75
+446_G >A rs33945777 E6 familial, Fetal HQTL-1, Hb SS disease, H-B-A, MTH β type, α- β-thalassemia. 0.008 0 0 0.002 0.19 1 0.50
+436_G >T rs35209591 E6 familial 0.39 0.41 0.42 0.41 0.61 0.93 0.56
+200_T >C rs111851677 β-thalassemia 0.02 0.03 0.02 0.02 0.79 1 1
+20_A >T rs334 HbS 0.03 0.07 0.03 0.05 0.03 0.38 (0.16-0.91) 0.26 0.45
+19_G >A rs33930165 HbC 0 0.01 0.03 0.01 0.26 0.17 0.008
+9_C >T rs713040 Hemoglobin OKAYAMA 0.11 0.07 0.12 0.10 0.10 0.03 0.51 (0.27-0.93) 0.66
−79_A >G rs34598529 E6 familial, Fetal HQTL-1, Hb SS disease, H-B-A, MTH β type, α- β-thalassemia. 0.004 0 0 0.001 0.44 1 1
−133_G >A rs72561473 β-thalassemia 0.012 0.019 0.010 0.014 0.737 0.486 1
−138_C >T rs33944208 β-thalassemia, Hemoglobinopathy 0 0 0.005 0.002 1 0.159 0.206

For the G6PD gene, PCR amplicon targeted G6PD-202 polymorphism amongst the commonly known SNPs in Western Africa and omitted polymorphisms such as G6PD +376, G6PD +542, G6PD +680 and G6PD +968. A set of 6 SNP were identified in an amplicon of 701pb (Table 3). Among them, 3 SNPs +10707 G > A, +10776 G > C, and +10983 G > C was newly identified genetic variants with low frequencies (MAF <3%). G6PD deficiency polymorphisms such as G6PD-202 G > A were 0.022, 0.032 and 0.018 in SM, UM and CTR groups, respectively. Another G6PD deficiency allele, such as +10588A/G (rs762515), was observed at a high frequency of >3% in the overall population. In specific clinical phenotype populations, the MAF were 0.28, 0.37 and 0.32 in SM, UM, and control groups. There was no evidence of genotypic deviation from HWE except +10588_A > G (rs762515).

Rather than G6PD, the HBB gene showed more polymorphisms in the Senegalese population with 12 characterised SNPs. These identified SNPs were previously registered in the NCBI database, and rs numbers were already attributed. Then our study does not reveal any novel HBB polymorphisms. Among them, 6 SNP (rs7946748, rs7480526, rs10768683, rs35209591, rs334 (HbS) and rs713040) were detected with high frequencies (with MAF > 3%), unlike 6 other SNPs (rs33945777, rs111851677, rs33930165 (HbC), rs34598529, rs72561473, rs33944208) with MAF < 3% were observed (Table 4). The MAF of the sickle cell HbS polymorphism was estimated to be 0.026, 0.069 and 0.035, and HbC polymorphism was estimated to be 0, 0.009, 0.029, in SM, UM and CTR groups, respectively.

The structure of G6PD and HBB genes in the Senegalese population was investigated. The presence of a haploblock structure gene was analysed by measuring the pairwise linkage disequilibrium (LD) in each pair of polymorphisms. Figure 2 shows the values of D’ between SNP. For G6PD rs762515 A > G, rs1050828 G > A, and for HBB rs10768683 G > C, rs7480526 T > G, rs7946748 C > T polymorphisms, the values D’ > 0.97, showing a high linkage disequilibrium, and defining a block structure in HBB genes. A genetic segment appeared to be part of a haploblock structure in the HBB loci, composed of 3 ‘successive’ SNPs, for which the confidence interval D′ was 0.9–1 (rs10768683 G > C, rs7480526 T > G, and rs7946748 C > T).

Figure 2. The HBB and G6PD polymorphisms located on chromosome chr11 and chr X, respectively.

Figure 2

The Haploview software provides the linkage disequilibrium map of HBB and G6PD genes. (A) The coding and UTR regions are indicated by black and white rectangles, respectively, for each loci. SNP positions are numbered per the initiation codon ATG considered as +1 (indicated by a black triangle). The genomic sequence used for alignment is GenBank sequence; reference genes were NG_009015.2 and NG_059281 for HBB and G6PD, respectively. The different colors are: label blue: SNP ‘promoter’ region, label ‘Black’, SNP‘intronic’ region and label ‘Red’, ‘exonic’ region. An underlined black line shows the boundaries of PCR amplicon in each gene. (B) The linkage disequilibrium map of LD plots within regions of 800 and 700 bases for HBB and G6PD genes, respectively. The LD plots show pairwise D′ values given in diamonds for each statistical comparison between the SNPs. The different shades of color represent D′ values (between 0 and 1). In white and red ‘diamond’, the values of D′ are indicated, an empty red diamond indicates that D′ > 0.97.

G6PD, HBB polymorphisms and clinical malaria

To test whether G6PD and HBB polymorphisms were associated with malaria protective effect in our overall population, a statistical association analysis was performed using logistic regression by comparing the three phenotypes groups, SM, UM and underage, sex as covariates. For G6PD, the SNP +10588 A > G (rs762515) yielded a significant association with protection from severe malaria in Senegalese populations. For SM vs UM, the SNP +10588 A > G (rs762515) showed a borderline p-value = 0.047 (OR 0.65, 95% CI [0.43–0.99]). For G6PD-202 G > A, no association with protection against severe malaria was found in our overall population (Table 3).

An analysis of HbS polymorphism showed a significant association with protection against severe malaria. For SM vs. UM comparison, the sickle cell trait HbS polymorphism (HBB +20 A > T, rs334) yielded a significant p-value = 0.033 (OR 0.38, 95% CI [0.16–0.91]) (Table 4). Surprisingly, the rs33930165 (HbC) polymorphism is not a protective variant against severe forms in our population. There is no significant difference between SM vs UM, with a p-value = 0.26 (Table 4).

G6PD, HBB polymorphisms and ecological malaria areas

We then performed comparisons of frequencies among the known malaria protective SNPs in 3 ecological areas and tried to distinguish relations with endemicity. The SNPs HbS (rs334), HbC (rs33930165) and G6PD +10588A > G (rs762515) were associated with protection against severe malaria in the three regions by comparison of SM vs UM with higher frequency in UM group (Fig. 3). However, HBB +9 C > T (rs713040) has variable protective effects depending on endemicity. In St-Louis, a low endemicity area, the SNP rs713040 and rs7946748 were associated with protecting the severe form, unlike Dakar and Tambacounda/Kolda/, which are more endemic areas. For G6PD-202 G > A, we found differences per the locality. The polymorphism was associated with protection against SM in the Dakar Capital centre but not in the other regions (St Louis-Louga and Tambacounda-Kolda). Indeed, the SNP is associated with the protection against SM form in Dakar but not in the other regions (Fig. 3).

Figure 3. Allelic frequencies of G6PD and HBB in different study areas.

Figure 3

Light grey, black grey and black areas indicate allelic frequencies from Dakar, Saint-Louis_Louga and Tamba_Kolda, respectively.

G6PD, HBB polymorphisms and biological parameters

Analyses were conducted to test whether G6PD and HBB polymorphisms were associated with biological parameters of severity, performing comparisons of SM vs UM. Then correlation of HBB polymorphisms HBB +226 C > T (rs7946748), HbS (rs334), HbC (rs33930165), HBB +9 C > T (rs713040), and G6PD-202 G > A with age, parasitaemia, and biological parameters (hemoglobin levels, blood platelets, lymphocytes, monocytes, basophils, neutrophils and eosinophils) and survival/death outcome, were performed using the Mann–Whitney test.

Our results showed a significant association between age and HBB +200_T > C (rs111851677), a variant responsible for β-thalassemia.

Then results showed significant association of rs7946748 with basophils levels p-value = 0.01 5 (beta = −0.5458, 95% CI [−0.85 to −0.11]) and lymphocytes p-value = 0.04 (beta = −8.59, 95% CI [−16.73 to −0.44]) (Fig. 4); the rs10768683 and rs713040 were associated to basophils with p-value = 0.023 (beta =−0.38, 95% CI [−0.7 to −0.06]) and p-value = 0.020 (beta = −0.37, 95% CI [−0.6 to −0.06]) respectively. HbS (rs334) triggered a significant association with eosinophils p-value = 0.024 (beta = 3.532, 95% CI [0.51–6.54]). Platelets were associated with G6PD-SNP +10707 G > A, G6PD-SNP +10983 G > C and HBB-rs72561473, with p = 0.009, 0.048, 0.006 respectively. Finally, G6PD- SNP +11000_C > T (rs782500951) was associated with parasitaemia with a p value = 0.016 (Fig. 4).

Figure 4. Forest plot of associations between HBB and G6PD polymorphisms and biological and cellular parameters.

Figure 4

Multivariate linear regression analysis with biological parameters. The blue squares represent βeta values. The horizontal lines around squares represent the 95% Confidence Intervals. SNP, single nucleotide polymorphism; MAF, minor allele frequency; beta, regression coefficient corresponding to the minor allele of each SNP; SE, Standard Error; BP, genomic location; PV, p-value; CHR, chromosome code; A1, Allele 1 (usually minor).

Discussion

In our present study, we characterised the polymorphisms of G6PD and HBB genes in the Senegalese population, two-loci previously associated with protection against severe malaria. A genetic analysis was performed to identify polymorphisms associated with clinical malaria and their prevalence, using three specific phenotypes of severity, including SM, UM and CTR groups living, from North to South, in three ecological malaria transmission areas.

We identified six polymorphisms on G6PD genes and found that the minor allele frequency of G6PD-202 G > A was 0.022, 0.032, 0.018 in SM, UM, and CTR groups, respectively. The frequency of these polymorphisms was 0.026 in the overall population. In The Gambia, a country neighbouring Senegal, it was reported that the minor allele frequency of this polymorphism in the CTR sample was 0.03, considerably lower than reports from other east African countries such as Kenya (0.18) and Malawi (0.19) (Clarke et al., 2017).

Our results confirm the weak prevalence of allele deficiency from G6PD-202 G > A in the Senegambian area compared to Kenya and Malawi populations. These results could ensure that G6PD polymorphisms prevalence is highly population-specific, and thus, attention should be paid to the study population from which future estimates may be drawn (Galatas et al., 2017). In our Senegalese population cohort, we found that G6PD-202 G > A polymorphism was not associated with protection against severe malaria at the global level. However, differences were founded in the different ecological localities. The polymorphism was associated with protection against SM in Dakar Capital centre but not in the other regions. Our results could contradict published data. Indeed the A-allele of G6PD-202 G > A polymorphism has been associated with protection against severe malaria in African populations, such as in Mali, the Gambia and Uganda (Guindo et al., 2007; Ruwende et al., 1995). However, this protective effect seems to depend on other parameters such as sex and locality. As reported by Maiga et al. (2014) the studies have shown the protective effect against severe malaria has not been consistent across extensive studies observed in females (Uyoga et al., 2015; Manjurano et al., 2015), in males (Guindo et al., 2007), in both (Ruwende et al., 1995), or no protection (Clark et al., 2009; Toure et al., 2012). This disparity between the results of the studies could be explained in part by variation in phenotype definition, choice of controls (village surveys vs hospital-based studies), age or immune status of participants, and study designs (case-control vs cohort) (Clark et al., 2009; Maiga et al., 2014). G6PD deficiency can potentially protect against uncomplicated malaria in African countries, but not severe malaria. Interestingly, this protection was mainly heterozygous, being x-linked and thus related to gender (Guindo et al., 2007; Mbanefo et al., 2017). No such protection was evident from the mosaic state of G6PD deficiency in heterozygous females. A previous study confirmed highly significant protection against severe malaria in hemizygous males but not in heterozygous females (Guindo et al., 2007). A case-control study showed a significant association between G6PD A- and risk of severe malaria, with protection against cerebral malaria but increased risk of severe anemia  (Clarke et al., 2017). However, compared to HbAS or α–thalassemia, associations between G6PD A- deficiency and risk of severe malaria have been less straightforward, with studies yielding inconsistent results (Manjurano et al., 2015; Mbanefo et al., 2017). In fact, it has been shown that the G6PD-202 polymorphism may not be a good marker of A- deficiency, and/or other polymorphisms are required to confirm the protective effect (Clark et al., 2009). Our study found that the G6PD SNP +10588_A > G (rs762515) yielded a significant association with protection from severe malaria in the Senegalese population.

In the HBB genes located on chromosome 11, we identified 12 polymorphisms and characterised their frequencies. The SNP HbS resulted in substituting a glutamate residue with a valine or lysine residue in the ß-globin chain at position 6. The Minor Allele Frequency of HbS and HbC polymorphisms was estimated to be 0.05 and 0.01, respectively, in the overall population. Several studies have shown a low frequency of the HbS allele in the control group (∼3.8%). The frequency is in keeping with other West African countries (Burkina Faso 5.2%, Cameroon 6.5%, The Gambia 7.6%, and Ghana 6.5%) and east African populations (Kenya 6.4%, Malawi 2.7%, Tanzania 7.8%) (https://malariaatlas.org/explorer/#/). A higher frequency of the HbC allele has been observed in other West African populations (Agarwal et al., 2000; Toure et al., 2012).

Our analysis confirmed the known protective effects of HbS. We found that HbS polymorphism showed a significant association with protection against severe malaria in the three ecological regions. For SM vs. UM, the sickle cell trait HbS polymorphism (HBB +20_A > T, rs334) yielded a significant p-value = 0.033. The observed frequencies for HbS (rs334) are 0.026, 0.069 and 0.035 for SM, UM, and CTR groups, respectively (Table 4). These data fit with studies and dogma on the protective effect of HbS polymorphism, corresponding to a lower frequency in ‘severe malaria form’ and an increase in the ‘uncomplicated’ malaria group and ‘healthy’ control subject, representing the general population (Table 4) (Lopera-Mesa et al., 2015; Kreuels et al., 2010; Amodu et al., 2012). In fact, the differences were not significant when comparing UMvs CTR, and SMvs CTR, probably to the fact of small numbers of controls and then, consequently, a lack of power.

Unexpectedly, the protective effect of the HbC variant has not been observed in our analysis. The HbC polymorphism (HBB +19_G > A, rs33930165) not showed significant; we found a p-value = 0.26 when comparing SM vs. UM. Previous studies showed that HbAS had been associated with 50% and 80% reduced risks of developing uncomplicated and severe malaria (Williams, 2005b); HbAC and HbCC have been associated with 30% and 93% reduced risks of developing mild and severe malaria, respectively (Modiano et al., 2001). A study conducted in Mali, a neighbouring country of Senegal, indicated that HbC is associated with protection against the SM form of P falciparum malaria in the Dogon of Bandiagara population (Clarke et al., 2017). In this ethnic group, the prevalence of HbC was significantly lower among severe malaria cases than among cases of uncomplicated malaria. The protective effect indicated an 80% reduction in the risk of severe malaria. The data also suggest that the protective effect associated with HbC may be greater than HbS in this population (Agarwal et al., 2000).

It was suggested a “balanced polymorphism” where the HbS homozygote disadvantage is repaid through the resistance of the heterozygote HbAS in regions where malaria is endemic (Lederberg, 1999).

Our results suggest differences in the protective effect of hemoglobin among the three ecological areas. The possibility that protective effects associated with different hemoglobin mutations vary among diverse human populations. This is consistent with the report’s findings. In fact, in West Africa, the Fulani ethnic group has decreased susceptibility to malaria compared to Dogon populations (Bereczky et al., 2006), and the Fulani also have a reduced prevalence of sickle hemoglobin  (Nasr et al., 2008), α-thalassemia (Modiano et al., 2001), and G6PD A- (Modiano et al., 2001; Maiga et al., 2014) compared to other groups.

We have also observed a positive association between HBB +200_T > C (rs111851677) variant and age. Moreover, this variant does not have a protective effect against severe malaria in our study. Age is an essential factor in the pathophysiology of malaria, and some reports from sub-Saharan Africa suggest that P. falciparum infections are more common among children aged 5–15 years than among younger children and adults (Cohee et al., 2020; Were et al., 2018; Toure et al., 2016). However, children are under-represented in our cohort because clinical forms were included without any preconceived notion of age categories. There are more adults because our inclusions are not oriented toward pediatric services; to fill this gap, this will be a point to consider in our future inclusions.

The results showed significant associations between several polymorphisms on HBB and G6PD genes, parasitemia, and white blood cells (basophils, eosinophils, and lymphocytes).

P. falciparum infection is associated with a profound T-cell activation with the Th1-derived cytokine interferon-γ as a major pathogenic factor (Hunt & Grau, 2003; Belnoue et al., 2008). Eosinophils and basophils seem to play important roles in malaria pathogenesis. Increases in plasma and tissue histamine derived from basophils have been associated with disease severity in human P. falciparum infections and several animal malaria (MacDonald et al., 2001). In addition, elevated plasma levels of IgE, which binds to basophils, have been associated with the severity of P. falciparum infection (Pelleau et al., 2012; Porcherie et al., 2011). Furthermore, increased eosinophil counts have been associated with recovery from infection  (Camacho et al., 1999; Pollenus et al., 2020). One of our recent studies showed that Ribonuclease 3 (RNASE 3), also known as eosinophil cationic protein (ECP), is one of the factors suggested as having a role in malaria severity. The polymorphisms: +371G/C (rs2073342), +499G/C (rs2233860) and +577A/T (rs8019343) defined a haplotype risk (G-G-T) for malaria severity (Diop et al., 2018). Despite these observations, little is known about the relationship between these cell types and disease. Our results showed that rs334 (HbS) is a variant associated with eosinophils (beta =3.532, 95% CI [0.51–6.54], p-value = 0.024) (Fig. 4). It would be interesting to study this association in-depth in future studies and compare eosinophils’ levels in patients with HBSS and HbAS variants.

The early and vigorous development of malaria immunity has been associated with a protective effect of sickle cell traits and thalassemia. Although, ethnic variations in the immunological response to P falciparum infection and rates of the clinical episodes of malaria are described in clinical studies in West Africa (Modiano et al., 2001; Agarwal et al., 2000).

In conclusion, our study gives an overview of the importance of sickle cell polymorphism on severe malaria susceptibility in the Senegalese population.

Supplemental Information

Supplemental Information 1. Raw Data.

Information about each participant: location, age, gender, malaria phenotype or control, haematological data, etc.

DOI: 10.7717/peerj.13487/supp-1
Supplemental Information 2. Biological data and SNPs for the HBB gene.

These data were used for statistical analysis.

DOI: 10.7717/peerj.13487/supp-2
Supplemental Information 3. Biological data and SNPs for the G6PD gene.

These data were used for statistical analysis.

DOI: 10.7717/peerj.13487/supp-3
Supplemental Information 4. Sequence data for HBB.
DOI: 10.7717/peerj.13487/supp-4
Supplemental Information 5. Sequences for G6PD.
DOI: 10.7717/peerj.13487/supp-5

Acknowledgments

The authors are grateful to all the patients and the medical staff who have generously collaborated in the Malaria Genomic Project Clayton Dedonder 2014 (Institut Pasteur de Dakar). The authors thank the staff of the Centre National de Recherche en Génétique humaine (CNRGH) for assistance in performing the genotyping. The CNRGH-CEA performed all sequencing of the target genes free of charge.

Funding Statement

This study was supported by grants from Dedonder Clayton 2014, Institut Pasteur Network. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Additional Information and Declarations

Competing Interests

The authors declare there are no competing interests.

Author Contributions

Fatou Thiam conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft.

Gora Diop conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft.

Cedric Coulonges analyzed the data, authored or reviewed drafts of the article, and approved the final draft.

Céline Derbois analyzed the data, authored or reviewed drafts of the article, and approved the final draft.

Babacar Mbengue conceived and designed the experiments, performed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Alassane Thiam conceived and designed the experiments, performed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Cheikh Momar Nguer conceived and designed the experiments, performed the experiments, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft.

Jean Francois Zagury analyzed the data, authored or reviewed drafts of the article, and approved the final draft.

Jean-Francois Deleuze analyzed the data, authored or reviewed drafts of the article, and approved the final draft.

Alioune Dieye conceived and designed the experiments, performed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Human Ethics

The following information was supplied relating to ethical approvals (i.e., approving body and any reference numbers):

UCAD Ethic committee (University Cheikh Anta DIOP, Senegal) (#Protocole 0344/2018/CER-UCAD#).

DNA Deposition

The following information was supplied regarding the deposition of DNA sequences:

The SNPs are available at GenBank: rs762515, rs1050828, rs782500951, rs7946748, rs7480526, rs10768683, rs33945777, rs35209591, rs111851677, rs334, rs33930165, rs713040, rs34598529, rs72561473, rs33944208.

Data Availability

The following information was supplied regarding data availability:

The raw data and sequences are available in the Supplementary Files.

References

  • Agarwal et al. (2000).Agarwal A, Guindo A, Cissoko Y, Taylor JG, Coulibaly D, Koné A, Kayentao K, Djimde A, Plowe CV, Doumbo O, Wellems TE, Diallo D. Hemoglobin C associated with protection from severe malaria in the Dogon of Mali, a West African population with a low prevalence of hemoglobin S. Blood. 2000;96(7):2358–2363. doi: 10.1182/blood.V96.7.2358. [DOI] [PubMed] [Google Scholar]
  • Agence Nationale de la Statistique et de la Démographie (ANSD) (2014).Agence Nationale de la Statistique et de la Démographie (ANSD) Senegal 2014 Population report. Dakar: ANSD 2014
  • Aldakeel et al. (2020).Aldakeel SA, Ghanem NZ, Al-Amodi AM, Osman AK, Al Asoom LI, Ahmed NR, Almandil NB, Akhtar MS, Azeez SA, Borgio JF. Identification of seven novel variants in the beta-globin gene in transfusion-dependent and normal patients. Archives of Medical Science. 2020;16(2):453–459. doi: 10.5114/aoms.2019.84825. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Amodu et al. (2012).Amodu OK, Olaniyan SA, Adeyemo AA, Troye-Blomberg M, Olumese PE, Omotade OO. Association of the sickle cell trait and the ABO blood group with clinical severity of malaria in southwest Nigeria. Acta Tropica. 2012;123(2):72–77. doi: 10.1016/j.actatropica.2012.03.013. [DOI] [PubMed] [Google Scholar]
  • Barrett et al. (2005).Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–265. doi: 10.1093/bioinformatics/bth457. [DOI] [PubMed] [Google Scholar]
  • Battistuzzi et al. (1977).Battistuzzi G, Esan GJ, Fasuan FA, Modiano G, Luzzatto L. Comparison of GdA and GdB activities in Nigerians. A study of the variation of the G6PD activity. American Journal of Human Genetics. 1977;29(1):31–36. [PMC free article] [PubMed] [Google Scholar]
  • Bei et al. (2017).Bei AK, Ahouidi AD, Dvorin JD, Miura K, Diouf A, Ndiaye D, Premji Z, Diakite M, Mboup S, Long CA, Duraisingh MT. Functional analysis reveals geographical variation in inhibitory immune responses against a polymorphic malaria antigen. The Journal of Infectious Diseases. 2017;216(2):267–275. doi: 10.1093/infdis/jix280. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Belnoue et al. (2008).Belnoue E, Potter SM, Rosa DS, Mauduit M, Grüner AC, Kayibanda M, Mitchell AJ, Hunt NH, Rénia L. Control of pathogenic CD8+ T cell migration to the brain by IFN-gamma during experimental cerebral malaria. Parasite Immunology. 2008;30(10):544–553. doi: 10.1111/j.1365-3024.2008.01053.x. [DOI] [PubMed] [Google Scholar]
  • Benito et al. (1994).Benito A, Roche J, Molina R, Amela C, Alvar J. Application and evaluation of QBC malaria diagnosis in a holoendemic area. Applied Parasitology. 1994;35(4):266–272. [PubMed] [Google Scholar]
  • Bereczky et al. (2006).Bereczky S, Dolo A, Maiga B, Hayano M, Granath F, Montgomery SM, Daou M, Arama C, Troye-Blomberg M, Doumbo OK, Färnert A. Spleen enlargement and genetic diversity of Plasmodium falciparum infection in two ethnic groups with different malaria susceptibility in Mali, West Africa. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2006;100(3):248–257. doi: 10.1016/j.trstmh.2005.03.011. [DOI] [PubMed] [Google Scholar]
  • Beutler (1989).Beutler E. Glucose-6-phosphate dehydrogenase: new perspectives. Blood. 1989;73(6):1397–1401. doi: 10.1182/blood.V73.6.1397.1397. [DOI] [PubMed] [Google Scholar]
  • Camacho et al. (1999).Camacho LH, Wilairatana P, Weiss G, Mercader MA, Brittenham GM, Looareesuwan S, Gordeuk VR. The eosinophilic response and haematological recovery after treatment for Plasmodium falciparum malaria. Tropical Medicine & International Health. 1999;4(7):471–475. doi: 10.1046/j.1365-3156.1999.00426.x. [DOI] [PubMed] [Google Scholar]
  • Carter et al. (2011).Carter N, Pamba A, Duparc S, Waitumbi JN. Frequency of glucose-6-phosphate dehydrogenase deficiency in malaria patients from six African countries enrolled in two randomized anti-malarial clinical trials. Malaria Journal. 2011;10:241. doi: 10.1186/1475-2875-10-241. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Clark et al. (2009).Clark TG, Fry AE, Auburn S, Campino S, Diakite M, Green A, Richardson A, Teo YY, Small K, Wilson J, Jallow M, Sisay-Joof F, Pinder M, Sabeti P, Kwiatkowski DP, Rockett KA. Allelic heterogeneity of G6PD deficiency in West Africa and severe malaria susceptibility. European Journal of Human Genetics. 2009;17(8):1080–1085. doi: 10.1038/ejhg.2009.8. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Clarke et al. (2017).Clarke GM, Rockett K, Kivinen K, Hubbart C, Jeffreys AE, Rowlands K, Jallow M, Conway DJ, Bojang KA, Pinder M, Usen S, Sisay-Joof F, Sirugo G, Toure O, Thera MA, Konate S, Sissoko S, Niangaly A, Poudiougou B, Mangano VD, Bougouma EC, Sirima SB, Modiano D, Amenga-Etego LN, Ghansah A, Koram KA, Wilson MD, Enimil A, Evans J, Amodu OK, Olaniyan S, Apinjoh T, Mugri R, Ndi A, Ndila CM, Uyoga S, Macharia A, Peshu N, Williams TN, Manjurano A, Sepulveda N, Clark TG, Riley E, Drakeley C, Reyburn H, Nyirongo V, Kachala D, Molyneux M, Dunstan SJ, Phu NH, Quyen NN, Thai CQ, Hien TT, Manning L, Laman M, Siba P, Karunajeewa H, Allen S, Allen A, Davis TME, Michon P, Mueller I, Molloy SF, Campino S, Kerasidou A, Cornelius VJ, Hart L, Shah SS, Band G, Spencer CCA, Agbenyega T, Achidi E, Doumbo OK, Farrar J, Marsh K, Taylor T, Kwiatkowski DP. Characterisation of the opposing effects of G6PD deficiency on cerebral malaria and severe malarial anaemia. Elife. 2017;6:e15085. doi: 10.7554/eLife.15085. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Cohee et al. (2020).Cohee LM, Opondo C, Clarke SE, Halliday KE, Cano J, Shipper AG, Barger-Kamate B, Djimde A, Diarra S, Dokras A, Kamya MR, Lutumba P, Ly AB, Nankabirwa JI, Njagi JK, Maiga H, Maiteki-Sebuguzi C, Matangila J, Okello G, Rohner F, Roschnik N, Rouhani S, Sissoko MS, Staedke SG, Thera MA, Turner EL, Van Geertruyden JP, Zimmerman MB, Jukes MCH, Brooker SJ, Allen E, Laufer MK, Chico RM. Preventive malaria treatment among school-aged children in sub-Saharan Africa: a systematic review and meta-analyses. The Lancet Global Health. 2020;8(12):e1499-e1511. doi: 10.1016/S2214-109X(20)30325-9. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Conway (2007).Conway DJ. Molecular epidemiology of malaria. Clinical Microbiology Reviews. 2007;20(1):188–204. doi: 10.1128/CMR.00021-06. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • De Araujo et al. (2006).De Araujo C, Migot-Nabias F, Guitard J, Pelleau S, Vulliamy T, Ducrocq R. The role of the G6PD AEth376G/968C allele in glucose-6-phosphate dehydrogenase deficiency in the seerer population of Senegal. Haematologica. 2006;91(2):262–263. [PubMed] [Google Scholar]
  • Delaneau et al. (2019).Delaneau O, Zagury J-F, Robinson MR, Marchini JL, Dermitzakis ET. Accurate, scalable and integrative haplotype estimation. Nature Communications. 2019;10(1):5436. doi: 10.1038/s41467-019-13225-y. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Diop et al. (2018).Diop G, Derbois C, Loucoubar C, Mbengue B, Ndao BN, Thiam F, Thiam A, Ndiaye R, Dieye Y, Olaso R, Deleuze J-F, Dieye A. Genetic variants of RNASE3 (ECP) and susceptibility to severe malaria in Senegalese population. Malaria Journal. 2018;17(1):61. doi: 10.1186/s12936-018-2205-9. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Diop et al. (2005).Diop S, Sene A, Cisse M, Toure AO, Sow O, Thiam D, Diakhate L. Prevalence and morbidity of G6PD deficiency in sickle cell disease in the homozygote. Dakar Medical. 2005;50(2):56–60. [PubMed] [Google Scholar]
  • Diouf et al. (2017).Diouf I, Rodriguez-Fonseca B, Deme A, Caminade C, Morse AP, Cisse M, Sy I, Dia I, Ermert V, Ndione J-A, Gaye AT. Comparison of malaria simulations driven by meteorological observations and reanalysis products in Senegal. International Journal of Environmental Research and Public Health. 2017;14(10):1119. doi: 10.3390/ijerph14101119. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Galatas et al. (2017).Galatas B, Mabote L, Simone W, Matambisso G, Nhamussua L, Mañú Pereira MDM, Menéndez C, Saute F, Macete E, Bassat Q, Alonso P, Aide P. Heterogeneity of G6PD deficiency prevalence in Mozambique: a school-based cross-sectional survey in three different regions. Malaria Journal. 2017;16(1):36. doi: 10.1186/s12936-016-1674-y. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Gomez-Manzo et al. (2016).Gomez-Manzo S, Marcial-Quino J, Vanoye-Carlo A, Serrano-Posada H, Ortega-Cuellar D, González-Valdez A, Castillo-Rodríguez RA, Hernández-Ochoa B, Sierra-Palacios E, Rodríguez-Bustamante E, Arreguin-Espinosa R. Glucose-6-phosphate dehydrogenase: update and analysis of new mutations around the world. International Journal of Molecular Sciences. 2016;17(12):2069. doi: 10.3390/ijms17122069. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Grant et al. (2015).Grant AV, Roussilhon C, Paul R, Sakuntabhai A. The genetic control of immunity to Plasmodium infection. BMC Immunology. 2015;16:14. doi: 10.1186/s12865-015-0078-z. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Griffiths et al. (2005).Griffiths MJ, Shafi MJ, Popper SJ, Hemingway CA, Kortok MM, Wathen A, Rockett KA, Mott R, Levin M, Newton CR, Marsh K, Relman DA, Kwiatkowski DP. Genomewide analysis of the host response to malaria in Kenyan children. The Journal of Infectious Diseases. 2005;191(10):1599–1611. doi: 10.1086/429297. [DOI] [PubMed] [Google Scholar]
  • Grignard et al. (2019).Grignard L, Mair C, Curry J, Mahey L, Bastiaens GJH, Tiono AB, Okebe J, Coulibaly SA, Goncalves BP, Affara M, Ouédraogo A, Bougouma EC, Sanou GS, Nébié I, Lanke KHW, Sirima SB, dAlessandro U, Clark TG, Campino S, Bousema T, Drakeley C. Bead-based assays to simultaneously detect multiple human inherited blood disorders associated with malaria. Malaria Journal. 2019;18(1):14. doi: 10.1186/s12936-019-2648-7. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Guindo et al. (2007).Guindo A, Fairhurst RM, Doumbo OK, Wellems TE, Diallo DA. X-linked G6PD deficiency protects hemizygous males but not heterozygous females against severe malaria. PLOS Medicine. 2007;4(3):e66. doi: 10.1371/journal.pmed.0040066. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Hedrick (2011).Hedrick PW. Population genetics of malaria resistance in humans. Heredity. 2011;107(4):283–304. doi: 10.1038/hdy.2011.16. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Hunt & Grau (2003).Hunt NH, Grau GE. Cytokines: accelerators and brakes in the pathogenesis of cerebral malaria. Trends in Immunology. 2003;24(9):491–499. doi: 10.1016/S1471-4906(03)00229-1. [DOI] [PubMed] [Google Scholar]
  • Jambou et al. (2001).Jambou R, Ranaivo L, Raharimalala L, Randrianaivo J, Rakotomanana F, Modiano D, Pietra V, Boisier P, Rabarijaona L, Rabe T, Raveloson N, De Giorgi F. Malaria in the highlands of Madagascar after five years of indoor house spraying of DDT. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2001;95(1):14–18. doi: 10.1016/S0035-9203(01)90317-7. [DOI] [PubMed] [Google Scholar]
  • Kreuels et al. (2010).Kreuels B, Kreuzberg C, Kobbe R, Ayim-Akonor M, Apiah-Thompson P, Thompson B, Ehmen C, Adjei S, Langefeld I, Adjei O, May J. Differing effects of HbS and HbC traits on uncomplicated falciparum malaria, anemia, and child growth. Blood. 2010;115(22):4551–4558. doi: 10.1182/blood-2009-09-241844. [DOI] [PubMed] [Google Scholar]
  • Lederberg (1999).Lederberg J. JBS Haldane (1949) on infectious disease and evolution. Genetics. 1999;153(1):1–3. doi: 10.1093/genetics/153.1.1. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Lopera-Mesa et al. (2015).Lopera-Mesa TM, Doumbia S, Konaté D, Anderson JM, Doumbouya M, Keita AS, Diakité SAS, Traoré K, Krause MA, Diouf A, Moretz SE, Tullo GS, Miura K, Gu W, Fay MP, Taylor SM, Long CA, Diakité M, Fairhurst RM. Effect of red blood cell variants on childhood malaria in Mali: a prospective cohort study. The Lancet Haematology. 2015;2(4):e140-e149. doi: 10.1016/S2352-3026(15)00043-5. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Lwanira et al. (2017).Lwanira CN, Kironde F, Kaddumukasa M, Swedberg G. Prevalence of polymorphisms in glucose-6-phosphate dehydrogenase, sickle haemoglobin and nitric oxide synthase genes and their relationship with incidence of uncomplicated malaria in Iganga, Uganda. Malaria Journal. 2017;16(1):322. doi: 10.1186/s12936-017-1970-1. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • MacDonald et al. (2001).MacDonald SM, Bhisutthibhan J, Shapiro TA, Rogerson SJ, Taylor TE, Tembo M, Langdon JM, Meshnick SR. Immune mimicry in malaria: Plasmodium falciparum secretes a functional histamine-releasing factor homolog in vitro and in vivo. Proceedings of the National Academy of Sciences of the United States of America. 2001;98(19):10829–10832. doi: 10.1073/pnas.201191498. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Mackinnon et al. (2005).Mackinnon MJ, Mwangi TW, Snow RW, Marsh K, Williams TN. Heritability of malaria in Africa. PLOS Medicine. 2005;2(12):e340. doi: 10.1371/journal.pmed.0020340. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Maiga et al. (2014).Maiga B, Dolo A, Campino S, Sepulveda N, Corran P, Rockett KA, Troye-Blomberg M, Doumbo OK, Clark TG. Glucose-6-phosphate dehydrogenase polymorphisms and susceptibility to mild malaria in Dogon and Fulani, Mali. Malaria Journal. 2014;13:270. doi: 10.1186/1475-2875-13-270. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Manjurano et al. (2015).Manjurano A, Sepulveda N, Nadjm B, Mtove G, Wangai H, Maxwell C, Olomi R, Reyburn H, Riley EM, Drakeley CJ, Clark TG. African glucose-6-phosphate dehydrogenase alleles associated with protection from severe malaria in heterozygous females in Tanzania. PLOS Genetics. 2015;11(2):e1004960. doi: 10.1371/journal.pgen.1004960. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Mason, Bautista & Gilsanz (2007).Mason PJ, Bautista JM, Gilsanz F. G6PD deficiency: the genotype-phenotype association. Blood Reviews. 2007;21(5):267–283. doi: 10.1016/j.blre.2007.05.002. [DOI] [PubMed] [Google Scholar]
  • Mauran-Sendrail et al. (1975).Mauran-Sendrail A, Bouloux C, Gomila J, Langaney A. Comparative study of haemoglobin types of two populations of eastern Senegal–Bedik and Niokholonko. Annals of Human Biology. 1975;2(2):129–136. doi: 10.1080/03014467500000671. [DOI] [PubMed] [Google Scholar]
  • Mbanefo et al. (2017).Mbanefo EC, Ahmed AM, Titouna A, Elmaraezy A, Trang NTH, Long NP, Anh NH, Nghi TD, Hung BT, Hieu MV, Anh NK, Huy NT, Hirayama K. Association of glucose-6-phosphate dehydrogenase deficiency and malaria: a systematic review and meta-analysis. Scientific Reports. 2017;7:45963. doi: 10.1038/srep45963. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Miller et al. (2002).Miller LH, Baruch DI, Marsh K, Doumbo OK. The pathogenic basis of malaria. Nature. 2002;415(6872):673–679. doi: 10.1038/415673a. [DOI] [PubMed] [Google Scholar]
  • Milner Jr (2018).Milner Jr DA. Malaria Pathogenesis. Cold Spring Harbor Perspectives in Medicine. 2018;8(1):a025569. doi: 10.1101/cshperspect.a025569. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Modiano et al. (2001).Modiano D, Luoni G, Sirima BS, Simporé J, Verra F, Konaté A, Rastrelli E, Olivieri A, Calissano C, Paganotti GM, D’Urbano L, Sanou I, Sawadogo A, Modiano G, Coluzzi M. Haemoglobin C protects against clinical Plasmodium falciparum malaria. Nature. 2001;414(6861):305–308. doi: 10.1038/35104556. [DOI] [PubMed] [Google Scholar]
  • Moxon et al. (2020).Moxon CA, Gibbins MP, McGuinness D, Milner Jr DA, Marti M. New insights into malaria pathogenesis. Annual Review of Pathology. 2020;15:315–343. doi: 10.1146/annurev-pathmechdis-012419-032640. [DOI] [PubMed] [Google Scholar]
  • Nasr et al. (2008).Nasr A, Elghazali G, Giha H, Troye-Blomberg M, Berzins K. Interethnic differences in carriage of haemoglobin AS and Fcgamma receptor IIa (CD32) genotypes in children living in eastern Sudan. Acta Tropica. 2008;105(2):191–195. doi: 10.1016/j.actatropica.2007.10.003. [DOI] [PubMed] [Google Scholar]
  • Ngom et al. (2007).Ngom PI, Benoist HM, Thiam F, Diagne F, Diallo PD. Influence of orthodontic anomalies on periodontal condition. Odonto-stomatologie Tropicale. 2007;30(118):9–16. [PubMed] [Google Scholar]
  • Niang et al. (2016).Niang M, Loucoubar C, Sow A, Diagne MM, Faye O, Faye O, Diallo M, Toure-Balde A, Sall AA. Genetic diversity of Plasmodium falciparum isolates from concurrent malaria and arbovirus co-infections in Kedougou, southeastern Senegal. Malaria Journal. 2016;15:155. doi: 10.1186/s12936-016-1208-7. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Onda et al. (2005).Onda M, Akaishi J, Asaka S, Okamoto J, Miyamoto S, Mizutani K, Yoshida A, Ito K, Emi M. Decreased expression of haemoglobin beta (HBB) gene in anaplastic thyroid cancer and recovery of its expression inhibits cell growth. British Journal of Cancer. 2005;92(12):2216–2224. doi: 10.1038/sj.bjc.6602634. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Papaioannou, Utzinger & Vounatsou (2019).Papaioannou I, Utzinger J, Vounatsou P. Malaria-anemia comorbidity prevalence as a measure of malaria-related deaths in sub-Saharan Africa. Scientific Reports. 2019;9(1):11323. doi: 10.1038/s41598-019-47614-6. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Pelleau et al. (2012).Pelleau S, Diop S, Badiane MD, Vitte J, Beguin P, Nato F, Diop BM, Bongrand P, Parzy D, Jambou R. Enhanced basophil reactivities during severe malaria and their relationship with the Plasmodium falciparum histamine-releasing factor translationally controlled tumor protein. Infection and Immunity. 2012;80(8):2963–2970. doi: 10.1128/IAI.00072-12. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Piel et al. (2013a).Piel FB, Howes RE, Nyangiri OA, Moyes CL, Williams TN, Weatherall DJ, Hay SI. Online biomedical resources for malaria-related red cell disorders. Human Mutation. 2013a;34(7):937–944. doi: 10.1002/humu.22330. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Piel et al. (2013b).Piel FB, Howes RE, Patil AP, Nyangiri OA, Gething PW, Bhatt S, Williams TN, Weatherall DJ, Hay SI. The distribution of haemoglobin C and its prevalence in newborns in Africa. Scientific Reports. 2013b;3:1671. doi: 10.1038/srep01671. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Pollenus et al. (2020).Pollenus E, Pham T-T, Vandermosten L, Possemiers H, Knoops S, Opdenakker G, Van den Steen PE. CCR2 is dispensable for disease resolution but required for the restoration of leukocyte homeostasis upon experimental malaria-associated acute respiratory distress syndrome. Frontiers in Immunology. 2020;11:628643. doi: 10.3389/fimmu.2020.628643. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Porcherie et al. (2011).Porcherie A, Mathieu C, Peronet R, Schneider E, Claver J, Commere P-H, Kiefer-Biasizzo H, Karasuyama H, Milon G, Dy M, Kinet J-P, Louis J, Blank U, Mécheri S. Critical role of the neutrophil-associated high-affinity receptor for IgE in the pathogenesis of experimental cerebral malaria. Journal of Experimental Medicine. 2011;208(11):2225–2236. doi: 10.1084/jem.20110845. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Prashanth (2012).Prashanth GP. Quantitative buffy coat (QBC) test for rapid diagnosis of malaria. European Journal of Pediatrics. 2012;171(1):199–200. doi: 10.1007/s00431-011-1589-0. author reply 201-202. [DOI] [PubMed] [Google Scholar]
  • Rodriguez, Gaunt & Day (2009).Rodriguez S, Gaunt TR, Day IN. Hardy-Weinberg equilibrium testing of biological ascertainment for Mendelian randomization studies. American Journal of Epidemiology. 2009;169(4):505–514. doi: 10.1093/aje/kwn359. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Ruwende et al. (1995).Ruwende C, Khoo SC, Snow RW, Yates SNR, Kwlatkowskl D, Gupta S, Warn P, Allsopp CEM, Gilbert SC, Peschu N, Newbold CI, Greenwood BM, Marsh K, Hill AVS. Natural selection of hemi- and heterozygotes for G6PD deficiency in Africa by resistance to severe malaria. Nature. 1995;376(6537):246–249. doi: 10.1038/376246a0. [DOI] [PubMed] [Google Scholar]
  • Sahondra Harisoa et al. (2001).Sahondra Harisoa LJ, Pietra V, Tombo ML, Albonico M, Ranaivo LH, De Giorgi F, Razanakolona J, D’Ancona FP, Sabatinelli G, Raveloson A, Modiano D, Rakotondramarina D. Epidemiologic surveillance system and control of malaria in the central highlands of Madagascar: results 1999–2000. Archives de l’Institut Pasteur de Madagascar. 2001;67(1–2):21–26. [PubMed] [Google Scholar]
  • Saissy, Rouvin & Koulmann (2003).Saissy JM, Rouvin B, Koulmann P. Severe malaria in intensive care units in 2003. Medecine Tropicale. 2003;63(3):258–266. [PubMed] [Google Scholar]
  • Takahashi et al. (2003).Takahashi M, Matsuda F, Margetic N, Lathrop M. Automated identification of single nucleotide polymorphisms from sequencing data. Journal of Bioinformatics and Computational Biology. 2003;1(2):253–265. doi: 10.1142/S021972000300006X. [DOI] [PubMed] [Google Scholar]
  • Taylor, Parobek & Fairhurst (2012).Taylor SM, Parobek CM, Fairhurst RM. Haemoglobinopathies and the clinical epidemiology of malaria: a systematic review and meta-analysis. The Lancet Infectious Diseases. 2012;12(6):457–468. doi: 10.1016/S1473-3099(12)70055-5. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Tetard et al. (2017).Tetard M, Milet J, Dechavanne S, Fievet N, Dorin-Semblat D, Elion J, Fairhurst RM, Deloron P, Tuikue-Ndam N, Gamain B. Heterozygous HbAC but not HbAS is associated with higher newborn birthweight among women with pregnancy-associated malaria. Scientific Reports. 2017;7(1):1414. doi: 10.1038/s41598-017-01495-9. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Thiam et al. (2021).Thiam LG, Nyarko PB, Kusi KA, Niang M, Aniweh Y, Awandare GA. Blood donor variability is a modulatory factor for P. falciparum invasion phenotyping assays. Scientific Reports. 2021;11(1):7129. doi: 10.1038/s41598-021-86438-1. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Thiam et al. (2020).Thiam L, Seck N, Diouf FN, Boiro D, Niang B, Sagna SD, Coundoul A, Faye PM, Ndiaye M, Fall AL, Ndiaye O. Clinical and paramedical features of epilepsy in children at the Ziguinchor Peace Hospital: a documentation review. The Pan African Medical Journal. 2020;37:387. doi: 10.11604/pamj.2020.37.387.21902. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Toure et al. (2012).Toure O, Konate S, Sissoko S, Niangaly A, Barry A, Sall AH, Diarra E, Poudiougou B, Sepulveda N, Campino S, Rockett KA, Clark TG, Thera MA, Doumbo O. Candidate polymorphisms and severe malaria in a Malian population. PLOS ONE. 2012;7(9):e43987. doi: 10.1371/journal.pone.0043987. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Toure et al. (2016).Toure M, Sanogo D, Dembele S, Diawara SI, Oppfeldt K, Schioler KL, Haidara DB, Traoré SF, Alifrangis M, Konradsen F, Doumbia S. Seasonality and shift in age-specific malaria prevalence and incidence in Binko and Carriere villages close to the lake in Selingue, Mali. Malaria Journal. 2016;15:219. doi: 10.1186/s12936-016-1251-4. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Travassos et al. (2015).Travassos MA, Coulibaly D, Laurens MB, Dembélé A, Tolo Y, Koné AK, Traoré K, Niangaly A, Guindo A, Wu Y, Berry AA, Jacob CG, Takala-Harrison S, Adams M, Shrestha B, Mu AZ, Kouriba B, Lyke KE, Diallo DA, Doumbo OK, Plowe CV, Thera MA. Hemoglobin C trait provides protection from clinical falciparum malaria in malian children. The Journal of Infectious Diseases. 2015;212(11):1778–1786. doi: 10.1093/infdis/jiv308. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Tregouet & Garelle (2007).Tregouet DA, Garelle V. A new JAVA interface implementation of THESIAS: testing haplotype effects in association studies. Bioinformatics. 2007;23(8):1038–1039. doi: 10.1093/bioinformatics/btm058. [DOI] [PubMed] [Google Scholar]
  • Uyoga et al. (2015).Uyoga S, Ndila CM, Macharia AW, Nyutu G, Shah S, Peshu N, Clarke GM, Kwiatkowski DP, Rockett KA, Williams TN. Glucose-6-phosphate dehydrogenase deficiency and the risk of malaria and other diseases in children in Kenya: a case-control and a cohort study. Lancet Haematol. 2015;2(10):e437-e444. doi: 10.1016/S2352-3026(15)00152-0. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Vergnes, Gherardi & Bouloux (1975).Vergnes H, Gherardi M, Bouloux C. Erythrocyte glucose-6-phosphate dehydrogenase in the Niokolonko (Malinke of the Niokolo) of Eastern Senegal. Identification of a slow variant with normal activity (Tacoma-like) Human Heredity. 1975;25(2):80–87. doi: 10.1159/000152711. [DOI] [PubMed] [Google Scholar]
  • Verra, Mangano & Modiano (2009).Verra F, Mangano VD, Modiano D. Genetics of susceptibility to Plasmodium falciparum: from classical malaria resistance genes towards genome-wide association studies. Parasite Immunology. 2009;31(5):234–253. doi: 10.1111/j.1365-3024.2009.01106.x. [DOI] [PubMed] [Google Scholar]
  • Walakira et al. (2017).Walakira A, Tukwasibwe S, Kiggundu M, Verra F, Kakeeto P, Ruhamyankaka E, Drakeley C, Dorsey G, Kamya MR, Nsobya SL, Rosenthal PJ. Marked variation in prevalence of malaria-protective human genetic polymorphisms across Uganda. Infection, Genetics and Evolution. 2017;55:281–287. doi: 10.1016/j.meegid.2017.09.021. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Weatherall & Clegg (2002).Weatherall DJ, Clegg JB. Genetic variability in response to infection: malaria and after. Genes & Immunity. 2002;3(6):331–337. doi: 10.1038/sj.gene.6363878. [DOI] [PubMed] [Google Scholar]
  • Were et al. (2018).Were V, Buff AM, Desai M, Kariuki S, Samuels A, Ter Kuile FO, Phillips-Howard PA, Patrick Kachur S, Niessen L. Socioeconomic health inequality in malaria indicators in rural western Kenya: evidence from a household malaria survey on burden and care-seeking behaviour. Malaria Journal. 2018;17(1):166. doi: 10.1186/s12936-018-2319-0. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • Williams (2015).Williams TN. Red blood cell variants and malaria: a long story not yet over. The Lancet Haematology. 2015;2(4):e130-e131. doi: 10.1016/S2352-3026(15)00049-6. [DOI] [PubMed] [Google Scholar]
  • Williams (2005a).Williams N. Africa focus raises malaria funds. Current Biology. 2005a;15(14):R528–R529. doi: 10.1016/j.cub.2005.07.012. [DOI] [PubMed] [Google Scholar]
  • Williams (2005b).Williams N. Malaria dogs African hopes. Current Biology. 2005b;15(13):R477–R478. doi: 10.1016/j.cub.2005.06.045. [DOI] [PubMed] [Google Scholar]
  • World Health Organization (2021).World Health Organization Geneva: WHOWorld malaria report 2020. 2021:1–178.

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Supplemental Information 1. Raw Data.

Information about each participant: location, age, gender, malaria phenotype or control, haematological data, etc.

DOI: 10.7717/peerj.13487/supp-1
Supplemental Information 2. Biological data and SNPs for the HBB gene.

These data were used for statistical analysis.

DOI: 10.7717/peerj.13487/supp-2
Supplemental Information 3. Biological data and SNPs for the G6PD gene.

These data were used for statistical analysis.

DOI: 10.7717/peerj.13487/supp-3
Supplemental Information 4. Sequence data for HBB.
DOI: 10.7717/peerj.13487/supp-4
Supplemental Information 5. Sequences for G6PD.
DOI: 10.7717/peerj.13487/supp-5

Data Availability Statement

The following information was supplied regarding data availability:

The raw data and sequences are available in the Supplementary Files.


Articles from PeerJ are provided here courtesy of PeerJ, Inc

RESOURCES