# Molecular Epidemiology of G6PD Genotypes in Different Ethnic Groups Residing in Saharan and Sahelian Zones of Mauritania

Oum Kelthoum Mamadou Djigo, Mohamed Salem Ould Ahmedou Salem, Sileye Mamadou Diallo, Mohamed Abdallahi Bollahi, Boushab Mohamed Boushab, Aymeric Garre, Nasserdine Papa Mze, Leonardo Basco, Sébastien Briolant, Ali Ould Mohamed Salem Boukhary

► To cite this version:

Oum Kelthoum Mamadou Djigo, Mohamed Salem Ould Ahmedou Salem, Sileye Mamadou Diallo, Mohamed Abdallahi Bollahi, Boushab Mohamed Boushab, et al. Molecular Epidemiology of G6PD Genotypes in Different Ethnic Groups Residing in Saharan and Sahelian Zones of Mauritania. Pathogens, 2021, 10 (8), pp.931. 10.3390/pathogens10080931. hal-03429269

HAL Id: hal-03429269

https://amu.hal.science/hal-03429269

Submitted on 14 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Molecular Epidemiology of G6PD Genotypes in Different Ethnic Groups Residing in Saharan and Sahelian Zones of Mauritania

Oum Kelthoum Mamadou Djigo<fnref n="1" />, Mohamed Salem Ould Ahmedou Salem<fnref n="1" />, Sileye Mamadou Diallo<fnref n="1" />, Mohamed Abdallahi Bollahi<fnref n="2" />, Boushab Mohamed Boushab<fnref n="3" />, Aymeric Garre<fnref n="4" />,<fnref n="5" />, Nasserdine Papa Mze<fnref n="4" />,<fnref n="5" />, Leonardo Basco<fnref n="4" />,<fnref n="5" /> ID, Sébastien Briolant<fnref n="4" />,<fnref n="5" />,<fnref n="6" /> ID and Ali Ould Mohamed Salem Boukhary <fnref n="1" />,* ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Djigo, O.K.M.; Ould Ahmedou Salem, M.S.; Diallo, S.M.; Bollahi, M.A.; Boushab, B.M.; Garre, A.; Papa Mze, N.; Basco, L.; Briolant, S.; Ould Mohamed Salem Boukhary, A. Molecular Epidemiology of G6PD Genotypes in Different Ethnic Groups Residing in Saharan and Sahelian Zones of Mauritania. *Pathogens* **2021**, *10*, 931. https://doi.org/10.3390/ pathogens10080931

Academic Editor: Sébastien Besteiro

Received: 18 March 2021

Accepted: 12 April 2021

Published: 23 July 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** *Plasmodium vivax* malaria is endemic in Mauritania. Individuals with glucose-6-phosphate de-hydrogenase (G6PD) deficiency may develop acute hemolytic anemia when exposed to 8-aminoquinoline antimalarial drugs, which are indispensable for a complete cure. The prevalence of G6PD allelic variants was assessed in different ethno-linguistic groups present in Mauritania. A total of 996 blood samples (447 males and 549 females; 499 white Moors and 497 individuals of black African ancestry) were collected from febrile patients in 6 different study sites: Aleg, Atar, Kiffa, Kobeni, Nouakchott, and Rosso. The presence of the African-type G6PD A- (G202A, A376G, A542T, G680T, and T968C mutations) and the Mediterranean-type G6PD B- (C563T) variants was assessed by PCR followed by restriction fragment length polymorphism and/or DNA sequencing. The prevalence of African-type G6PD A- genotype was 3.6% (36/996), with 6.3% (28/447) of hemizygote (A-) males and 1.5% (8/549) of homozygous (A-A-) females. Forty of 549 (7.3%) women were heterozygous (AA-). The following genotypes were observed among hemizygous men and/or homozygous women: A376G/G202A (22/996; 2.2%), A376G/T968C Betica-Selma (12/996; 1.2%), and A376G/A542T Santamaria (2/996; 0.2%). The Mediterranean-type G6PD B- genotype was not observed. The prevalence rates of G6PD A- genotype in male (10/243; 4.1%) and heterozygous female (6/256; 2.3%) white Moors were lower (*p* < 0.05) than those of males (18/204; 8.8%) and heterozygous females (34/293; 11.6%) of black African ancestry. There were only a few homozygous women among both white Moors (3/256; 1.2%) and those of black African ancestry (5/293; 1.7%). The prevalence of G6PD deficiency in Mauritania was comparable to that of neighboring countries in the Maghreb. Because of the purportedly close ethnic ties between the Mauritanian white Moors and the peoples in the Maghreb, further investigations on the possible existence of the Mediterranean-type allele are required. Moreover, a surveillance system of G6PD phenotype and/or genotype screening is warranted to establish and monitor a population-based prevalence of G6PD deficiency.

**Keywords:** DNA sequencing; ethnic groups; epidemiology; G6PD deficiency; malaria elimination; Mauritania; PCR-RFLP; *Plasmodium vivax*

<footnote n="1">¹ Unité de Recherche "Génomes et Milieux" (Jeune Equipe Associée à l'Institut de Recherche pour le Développement), Faculté des Sciences et Techniques, Université de Nouakchott Al-Aasriya, Nouakchott, Mauritania; mari1702@live.fr (O.K.M.D.); salem0606@yahoo.fr (M.S.O.A.S.); diallosiley08@gmail.com (S.M.D.)</footnote>

<footnote n="2">² Institut National de Recherche en Santé Publique (INRSP), Nouakchott, Mauritania; bollahi@yahoo.com</footnote>

<footnote n="3">³ Department of Internal Medicine and Infectious Diseases, Kiffa Regional Hospital, Assaba, Mauritania; bboushab@gmail.com</footnote>

<footnote n="4">⁴ Aix Marseille Université, Institut de Recherche pour le Développement (IRD), Assistance Publique-Hôpitaux de Marseille (AP-HM), Service de Santé des Armées (SSA), Vecteurs—Infections Tropicales et Méditerranéennes (VITROME), 13005 Marseille, France; aymericgarre@gmail.com (A.G.); npapamze@gmail.com (N.P.M.); lkbasco@yahoo.fr (L.B.); sbriolant@wanadoo.fr (S.B.)</footnote>

<footnote n="5">5 Institut Hospitalo-Universitaire (IHU)—Méditerranée Infection, 13005 Marseille, France</footnote>

<footnote n="6">⁶ Unité de Parasitologie Entomologie, Département de Microbiologie et Maladies Infectieuses, Institut de Recherche Biomédicale des Armées (IRBA), 13005 Marseille, France</footnote>

<footnote>* Correspondence: alimedsalem@gmail.com</footnote>
# 1. Introduction

Glucose-6-phosphate dehydrogenase (G6PD) is the key enzyme of pentose phosphate pathway, the unique source of the reduced form of nicotinamide adenine dinucleotide phosphate (NADPH) in mature human erythrocytes. NADPH plays a key role in enzyme-based antioxidant system to protect erythrocytes from oxidative damage. When mature erythrocytes are deficient in G6PD and are exposed to various exogenous trigger factors, including infections, a variety of drugs (e.g., 8-aminoquinoline antimalarial drugs, which include primaquine and tafenoquine), and vicine and convicine found in fava beans, the antioxidant defense system is unable to neutralize free radicals that accumulate in the erythrocytes, leading to lipid peroxidation and cell lysis, which may be manifested clinically by acute hemolytic anemia [??].

In humans, G6PD deficiency is an X-linked recessive genetic trait which occurs frequently, affecting approximately 400 million people worldwide [??]. The G6PD gene consists of 13 exons and 12 introns encoding a protein with a mass of 59 kDa [?]. The gene is highly polymorphic with more than 200 mutations reported worldwide [?]. Despite the high number of mutant G6PD variants, only a few of them result in a functionally deficient G6PD enzyme. On the African continent, the most prevalent G6PD variants of clinical importance are the African-type G6PD A- found in sub-Saharan Africa and the Mediterranean-type G6PD B- variant mainly found in North Africa. Most G6PD mutations are clinically silent unless the G6PD-deficient individual is challenged with a trigger factor. The main clinical manifestation in such individuals is potentially fatal non-immune acute hemolytic anemia.

G6PD deficiency is a major issue in areas where *Plasmodium vivax* malaria is present because 8-aminoquinolines are the only currently available drugs that can kill hypnozoites, the latent liver forms of *P. vivax* (and also *P. ovale*). The hypnozoites remain dormant for weeks or months after the primoinfection and are reactivated due to still unknown mechanisms, leading to new malarial attacks called relapse [??]. Although most acute attacks due to *P. vivax* are not fatal, each attack may worsen anemia. The therapeutic objective in *P. vivax*-infected patients is complete elimination of both blood and liver stage parasites [?]. However, primaquine and tafenoquine given at the standard doses for anti-hypnozoite therapy are contraindicated in G6PD-deficient patients. These patients are therefore subject to multiple *P. vivax* (or *P. ovale*) malaria attacks from a primoinfection and can also become a source of further transmission to other persons, rendering the objective of malaria elimination difficult to attain.

The endemicity of *P. vivax* malaria has been firmly established over the past two decades in the Saharan zone of Mauritania, including Nouakchott, the capital, and Atar, an oasis city [??]. Faced with this new reality, since 2013, the Mauritanian health authorities recommend the administration of primaquine, in addition to either artesunate-amodiaquine or artemether-lumefantrine, to prevent relapse in *P. vivax*-infected patients with normal G6PD activity [?]. However, this treatment regimen has not yet been deployed in health facilities due to the lack of point-of-care G6PD testing prior to primaquine administration.

At present, the amount of data available on G6PD genotypes (or phenotypes) in Mauritania is limited to two studies conducted in Nouakchott [?]. It is known that G6PD polymorphisms are associated with ethnicity [?]. Mauritania is home to different ethnolinguistic groups distributed into two large groups, the Moorish group comprising the white Moors (also locally called Bidhanes or Beydanes) and the black Moors (also called Haratin) and the Black African group made up of the Pular, Soninke and Wolof [?]. There are no census-based official numbers of people belonging to each ethnic group because the Mauritanian government no longer publishes the proportions of ethnic group since the 1977 census.

We hypothesized that some white Moors, who are of North African Arab-Berber descent, may have inherited the Mediterranean-type G6PD B- genotype [??]. All other ethnic groups are of black African ancestry. Based on this information, as well as on earlier studies carried
out in sub-Saharan African countries bordering Mauritania and in Mauritania, it was further hypothesized that some individuals belonging to the ethnic groups of black African descent may carry the African-type G6PD A- genotype [???????].

Given the scarcity of data on G6PD genotypes in the general Mauritanian population, the present study was conducted to determine the prevalence of African-type and Mediterranean-type G6PD genotypes in both men and women of all ethnic groups present in Mauritania and establish a baseline genetic database as the first step to the full implementation of primaquine therapy in *P. vivax*-infected Mauritanian patients in the near future.

## 2. Results

### 2.1. Characteristics of Participants

A total of 996 blood samples (447 men and 549 women) were analyzed for G6PD polymorphisms. The following numbers of samples collected from 6 study sites were genotyped: 101 from Atar, 226 from Nouakchott, 79 from Aleg, 20 from Kiffa, 323 from Kobeni, and 247 from Rosso (????). Based on ethnic groups, 499 (50.1%) blood spots were from white Moors, 316 (31.7%) from black Moors, 134 (13.5%) from Pular, 24 (2.4%) from Wolof, and 23 (2.3%) from Soninke.

<table><caption>**Table 1.** Frequency of G6PD*A and B genotypes in males in six localities in Mauritania.</caption><thead><tr><th rowspan="2">Locality</th><th rowspan="2">Ethnic Group</th><th rowspan="2">N</th><th colspan="3">Genotype</th></tr><tr><th>A-</th><th>A</th><th>B</th></tr></thead><tbody><tr><td rowspan="2">Atar</td><td>White Moors</td><td>52</td><td>4 (7.7)</td><td>2 (3.8)</td><td>46 (88.5)</td></tr><tr><td>Black Moors</td><td>2</td><td>0</td><td>0</td><td>2 (100)</td></tr><tr><td rowspan="6">Nouakchott</td><td>White Moors</td><td>71</td><td>1 (1.4)</td><td>2 (2.8)</td><td>68 (95.8)</td></tr><tr><td>Black Moors</td><td>19</td><td>1 (5.3)</td><td>0</td><td>18 (94.7)</td></tr><tr><td>Pular</td><td>9</td><td>0</td><td>1 (11.1)</td><td>8 (88.9)</td></tr><tr><td>Wolof</td><td>2</td><td>0</td><td>0</td><td>2 (100)</td></tr><tr><td>Soninke</td><td>1</td><td>0</td><td>0</td><td>1 (100)</td></tr><tr><td>White Moors</td><td>29</td><td>0</td><td>0</td><td>29 (100)</td></tr><tr><td rowspan="4">Aleg</td><td>Black Moors</td><td>12</td><td>0</td><td>2 (16.7)</td><td>10 (83.3)</td></tr><tr><td>Pular</td><td>1</td><td>0</td><td>0</td><td>1 (100)</td></tr><tr><td>Wolof</td><td>1</td><td>1 (100)</td><td>0</td><td>0</td></tr><tr><td>White Moors</td><td>9</td><td>0</td><td>2 (22.2)</td><td>7 (77.8)</td></tr><tr><td rowspan="5">Kiffa</td><td>Black Moors</td><td>2</td><td>0</td><td>1 (50.0)</td><td>1 (50.0)</td></tr><tr><td>White Moors</td><td>41</td><td>2 (4.9)</td><td>6 (14.6)</td><td>33 (80.5)</td></tr><tr><td>Black Moors</td><td>58</td><td>7 (12.1)</td><td>12 (20.7)</td><td>39 (67.2)</td></tr><tr><td>Pular</td><td>21</td><td>0</td><td>7 (33.3)</td><td>14 (66.7)</td></tr><tr><td>Soninke</td><td>4</td><td>0</td><td>1 (25.0)</td><td>3 (75.0)</td></tr><tr><td rowspan="5">Kobeni</td><td>White Moors</td><td>41</td><td>3 (7.3)</td><td>3 (7.3)</td><td>35 (85.4)</td></tr><tr><td>Black Moors</td><td>47</td><td>7 (14.9)</td><td>11 (23.4)</td><td>29 (61.7)</td></tr><tr><td>Pular</td><td>19</td><td>2 (10.5)</td><td>2 (10.5)</td><td>15 (79.0)</td></tr><tr><td>Wolof</td><td>5</td><td>0</td><td>0</td><td>5 (100)</td></tr><tr><td>Soninke</td><td>1</td><td>0</td><td>0</td><td>1 (100)</td></tr><tr><td colspan="2">Total</td><td>447</td><td>28 (6.3)</td><td>52 (11.6)</td><td>367 (82.1)</td></tr></tbody></table>
<table><caption>**Table 2.** Frequency of G6PD\*A and B genotypes in females in six localities in Mauritania.</caption><thead><tr><th rowspan="2">Locality</th><th rowspan="2">Ethnic Group</th><th rowspan="2">N</th><th colspan="6">Genotype</th></tr><tr><th>A-A-</th><th>AA-</th><th>BA-</th><th>AA</th><th>BB</th><th>BA</th></tr></thead><tbody><tr><td rowspan="2">Atar</td><td>White Moors</td><td>47</td><td>0</td><td>0</td><td>0</td><td>1 (2.1)</td><td>44 (93.6)</td><td>2 (4.3)</td></tr><tr><td>White Moors</td><td>83</td><td>1 (1.2)</td><td>0</td><td>0</td><td>1 (1.2)</td><td>70 (84.3)</td><td>11 (13.3)</td></tr><tr><td rowspan="4">Nouakchott</td><td>Black Moors</td><td>31</td><td>1 (3.2)</td><td>1 (3.2)</td><td>0</td><td>1 (3.2)</td><td>21 (67.8)</td><td>7 (22.6)</td></tr><tr><td>Pular</td><td>5</td><td>1 (20.0)</td><td>0</td><td>0</td><td>0</td><td>4 (80.0)</td><td>0</td></tr><tr><td>Wolof</td><td>3</td><td>0</td><td>0</td><td>0</td><td>1 (33.3)</td><td>2 (66.7)</td><td>0</td></tr><tr><td>Soninke</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2 (100)</td><td>0</td></tr><tr><td rowspan="3">Aleg</td><td>White Moors</td><td>23</td><td>1 (4.3)</td><td>0</td><td>0</td><td>0</td><td>21 (91.4)</td><td>1 (4.3)</td></tr><tr><td>Black Moors</td><td>11</td><td>0</td><td>0</td><td>1 (9.1)</td><td>2 (18.2)</td><td>6 (54.5)</td><td>2 (18.2)</td></tr><tr><td>Pular</td><td>2</td><td>0</td><td>0</td><td>1 (50.0)</td><td>0</td><td>1 (50.0)</td><td>0</td></tr><tr><td rowspan="3">Kiffa</td><td>White Moors</td><td>5</td><td>0</td><td>0</td><td>0</td><td>0</td><td>5 (100)</td><td>0</td></tr><tr><td>Black Moors</td><td>3</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1 (33.3)</td><td>2 (66.7)</td></tr><tr><td>Pular</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1 (100)</td></tr><tr><td rowspan="5">Kobeni</td><td>White Moors</td><td>57</td><td>1 (1.8)</td><td>1 (1.8)</td><td>3 (5.2)</td><td>1 (1.8)</td><td>48 (84.2)</td><td>3 (5.2)</td></tr><tr><td>Black Moors</td><td>77</td><td>1 (1.3)</td><td>4 (5.2)</td><td>10 (13.0)</td><td>6 (7.8)</td><td>38 (49.3)</td><td>18 (23.4)</td></tr><tr><td>Pular</td><td>48</td><td>0</td><td>2 (4.2)</td><td>2 (4.2)</td><td>0</td><td>31 (64.6)</td><td>13 (27.0)</td></tr><tr><td>Wolof</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2 (100)</td><td>0</td></tr><tr><td>Soninke</td><td>15</td><td>0</td><td>1 (6.7)</td><td>0</td><td>1 (6.7)</td><td>10 (66.6)</td><td>3 (20.0)</td></tr><tr><td rowspan="4">Rosso</td><td>White Moors</td><td>41</td><td>0</td><td>0</td><td>2 (4.9)</td><td>2 (4.9)</td><td>29 (70.7)</td><td>8 (19.5)</td></tr><tr><td>Black Moors</td><td>54</td><td>1 (1.9)</td><td>1 (1.9)</td><td>4 (7.4)</td><td>4 (7.4)</td><td>30 (55.5)</td><td>14 (25.9)</td></tr><tr><td>Pular</td><td>28</td><td>0</td><td>1 (3.6)</td><td>5 (17.9)</td><td>2 (7.1)</td><td>16 (57.1)</td><td>4 (14.3)</td></tr><tr><td>Wolof</td><td>11</td><td>1 (9.1)</td><td>0</td><td>1 (9.1)</td><td>2 (18.2)</td><td>4 (36.4)</td><td>3 (27.2)</td></tr></tbody><tfoot><tr><td colspan="2">Total</td><td>549</td><td>8 (1.5)</td><td>11 (2.0)</td><td>29 (5.3)</td><td>24 (4.4)</td><td>385 (70.1)</td><td>92 (16.7)</td></tr></tfoot></table>

## 2.2. Prevalence and Type of G6PD Genotype and Gender

Overall, 419 of 447 (93.7%, 95% confidence interval (95% CI), 91.1–95.8%) men had normal G6PD genotypes A (n = 52) or B (n = 367), and 28 (6.3%, CI95% (4.2–8.9%)) were A-(Table ??). Among women, 501 of 549 (91.3%, 95% CI (88.6–93.5%)) carried the normal G6PD genotypes AA, BB or AB, whereas 40 (7.3%, 95% CI (5.3–9.8%)) were heterozygous (AA- or BA-), and 8 (1.5%, 95% CI (0.6–2.9%)) were African-type homozygous A-A- (Table ??). The global prevalence of African-type genotypes probably associated with G6PD deficiency therefore consisted of 6.3% (28/447, 95% CI (4.2–8.9%)) hemizygous males and 1.5% (8/549, 95% CI (0.6–2.9%)) homozygous females, i.e., a total of 3.6% (36/996, 95% CI (2.5–5.0%)) in both sexes.

Among the individuals with G6PD A- genotypes, G6PD A<sup>-</sup>(202) or A376G/G202A double mutation was observed in 2.2% (22/996, 95% CI (1.4–3.3%)); 19 hemizygous males and 3 homozygous females) (????). The second most frequently observed African-type G6PD genotype associated with G6PD deficiency was Betica-Selma or A376G/T968C double mutation (12/996; 1.2%, 95% CI (0.6–2.1%)); 7 males and 5 homozygous females). Santamaria G6PD A- or A376G/A542T was observed in 2 males, i.e., with a frequency of 0.2% (2/996, 95% CI (0–0.7%)); there was no homozygous Mauritanian females with this African-type G6PD genotype. Among heterozygous females, the most common genotype was 376A/G-202G/A (n = 23/549 females, 4.2%, 95% CI (2.7–6.2%)), followed by 376A/G-T968T/C (n = 15/549 females, 2.7%, 95% CI (1.5–4.5%)) and 376A/G-542A/T (n = 2/549 females, 0.4%, 95% CI (0–1.3%)). We did not find any individuals carrying the A376G/G680T African-type
G6PD genotype. Likewise, none of the patients included in the present study had the Mediterranean-type G6PD B- genotype.

<table><caption>**Table 3.** Frequency of G6PD*A- allelic variants in males with G6PD*A- genotypes in Mauritania.</caption><thead><tr><th rowspan="2">Ethnic Group</th><th rowspan="2">n/N (%)</th><th colspan="4">G6PD*A- Allelic Variants (n; % Among Men with African-Type G6PD Genotypes)</th></tr><tr><th>202A</th><th>542T</th><th>680T</th><th>968C</th></tr></thead><tbody><tr><td>White Moors</td><td>10/243 (4.1)</td><td>6 (21.4)</td><td>2 (7.1)</td><td>0</td><td>2 (7.1)</td></tr><tr><td>Black Moors</td><td>15/140 (10.7)</td><td>12 (42.8)</td><td>0</td><td>0</td><td>3 (10.7)</td></tr><tr><td>Pular</td><td>2/50 (4.0)</td><td>0</td><td>0</td><td>0</td><td>2 (7.1)</td></tr><tr><td>Wolof</td><td>1/8 (12.5)</td><td>1 (3.6)</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Soninke</td><td>0/6 (0)</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Total</td><td>28/447 (6.3)</td><td>19</td><td>2</td><td>0</td><td>7</td></tr></tbody></table>

<table><caption>**Table 4.** Frequency of G6PD*A- allelic variants in females with G6PD*A- genotypes in Mauritania.</caption><thead><tr><th rowspan="3">Ethnic Group</th><th rowspan="3">n/N (%)</th><th colspan="6">G6PD*A- Allelic Variants (n; % Among Females with the African-Type G6PD Genotypes)</th></tr><tr><th colspan="3">Heterozygote</th><th colspan="3">Homozygote</th></tr><tr><th>202G/A</th><th>542A/T</th><th>968T/C</th><th>202A/A</th><th>968C/C</th></tr></thead><tbody><tr><td>White Moors</td><td>9/256 (3.5)</td><td>5 (10.4)</td><td>1 (2.1)</td><td>0</td><td>2 (4.2)</td><td>1 (2.1)</td></tr><tr><td>Black Moors</td><td>24/176 (13.6)</td><td>14 (29.2)</td><td>1 (2.1)</td><td>6 (12.5)</td><td>1 (2.1)</td><td>2 (4.2)</td></tr><tr><td>Pular</td><td>12/84 (14.3)</td><td>3 (6.2)</td><td>0</td><td>8 (16.7)</td><td>0</td><td>1 (2.1)</td></tr><tr><td>Wolof</td><td>2/16 (12.5)</td><td>1 (2.1)</td><td>0</td><td>0</td><td>0</td><td>1 (2.1)</td></tr><tr><td>Soninke</td><td>1/17 (5.9)</td><td>0</td><td>0</td><td>1 (2.1)</td><td>0</td><td>0</td></tr><tr><td>Total</td><td>48/549 (8.7)</td><td>23 (47.9)</td><td>2 (4.2)</td><td>15 (31.2)</td><td>3 (6.2)</td><td>5 (10.4)</td></tr></tbody></table>

Data are expressed as the number of individual males belonging to one of the ethnic groups (n) and percentage in each locality. G6PD genotype: African-type genotype in hemizygous male, A-; African-type normal male with a single A376G mutation, A; normal male, B. There were no male Pulars in Kiffa and no male Wolofs in Kobeni.

Data are expressed as the number of individual females belonging to one of the ethnic groups (n) and percentage in each locality. G6PD genotype: African-type genotype in homozygous female, A-A-; African-type genotype in female heterozygous carrier, AA- or BA-; normal female with a single African-type A376G mutation with either AA or BA and normal female without any mutation, BB. There were no female Black Moors in Atar, no female Wolofs in Aleg, and no female Soninkes in Rosso in the present study.

## 2.3. Prevalence and Type of G6PD Genotype and Ethnic Groups

???? summarize the data on individuals with African-type G6PD genotype in relation to ethnic groups. Among white Moors with African-type G6PD variant (19/499, 3.8%, 95% CI (2.3–5.9%)), 10 males were hemizygous, 3 females were homozygous, and 6 females were heterozygous. Among black Moors with African-type G6PD variant (39/316, 12.3%, 95% CI (8.9–16.5%)), 15 males were hemizygous, 3 females were homozygous, and 21 females were heterozygous. Among the individuals belonging to the minority ethnic groups composed of Pulars, Wolofs, and Soninkes (19/181, 10.5%, 95% CI (6.4–15.9%)), 3 males were hemizygous, 2 females were homozygous, and 13 females were heterozygous. Taken
together, there were no statistically significant difference between the number of males of black African ancestry (18/204; 8.8%, 95% CI (5.3–13.6%)) and the male white Moors (10/243; 4.1%, 95% CI (2.0–7.4%)), who carried African-type G6PD genotype associated with G6PD deficiency (p = 0.0772). Likewise, there were significantly more heterozygous females among those of black African origin (34/293; 11.6%, 95% CI (7.9–15.3%)) than heterozygous white Moor women (6/256; 2.3%, 95% CI (0.9–5.0%)) (p < 0.0001). There was no statistically significant difference between the number of homozygous females of black African ancestry (5/293; 1.7%, 95% CI (0.6–3.9%)) and that of homozygous female white Moors (3/256; 1.2%, 95% CI (0.2–3.4%)) (p = 0.73), but the numbers of homozygous females were too low for a meaningful statistical comparison.

Data expressed as the number (n) of affected males with African-type genotype among the total number (N) of males belonging to each ethnic group and percentage of affected males in each ethnic group in parentheses. The proportion of A- allelic variants refers to the type of A- variants in relation to 28 males (i.e., denominator) with the African-type G6PD variant. There were a total of 447 men in the study.

Data expressed as the number (n) of affected females with African-type genotype among the total number (N) of females belonging to each ethnic group and percentage of affected females in each ethnic group in parentheses. The proportion of A- allelic variants refers to the type of A- variants in relation to 48 heterozygous or homozygous females (i.e., denominator) with the African-type G6PD variant. There were a total of 549 females in the study. Heterozygous 680G/T and homozygous 542T/T and 680T/T genotypes were not observed.

## 2.4. Distribution of G6PD Genotypes in Mauritania

The distribution of G6PD genotypes in six study sites in Mauritania is depicted in Figure 2. Hemizygous or homozygous genotypes associated with the African-type G6PD deficiency were identified in all six study sites, except for Kiffa. However, the situation in Kiffa was not representative of the general population in the city due to the small number of samples available from this study site. If data from Kiffa are excluded, the proportions of hemizygous and homozygous genotypes associated with G6PD deficiency ranged from 2.2% (in Aleg and Nouakchott) to 5.7% (in Rosso). Heterozygous females tended to be more present in southern Sahelian regions of the country where relatively high numbers of Mauritanians of black African ancestry reside, as compared to the desertic regions in the north.
<figure><img src="image_4.png" /><figcaption>Distribution of *G6PD* genotypes in Mauritania.</figcaption></figure>

**Figure 1.** The proportions of *G6PD* normal genotypes, genotypes associated with G6PD deficiency, and heterozygous females in six study sites in Mauritania are presented. Blue (AA, BB, and AB normal females and A and B normal hemizygous males); red, A- hemizygous males and A-A-homozygous females; green, AA- and BA- heterozygous females.

**3. Discussion**

One of the hypotheses of the present study that some white Moors may carry the Mediterranean-type G6PD deficiency seems to be invalidated. The results of the present study are in contradiction with those of our earlier study based on multiplex allele-specific PCR [? ]. In that study, of 44 phenotypically G6PD deficient samples according to the rapid diagnostic test of G6PD, 13 (8 white Moors, 3 black Moors, and 2 black Africans) were reported to have the Mediterranean variant. These 13 samples were re-analyzed by the same multiplex allele-specific PCR kit, PCR-RFLP protocol described in the present study, and DNA sequencing, and they were all characterized to be of normal *G6PD* genotype by three methods (data not shown). We can therefore tentatively conclude that, based on more than 500 samples (499 from the present study + 8 from the earlier study) from white Moors in six different study sites in Mauritania, the Mediterranean-type G6PD B- variant was not found in our study.

It can be hypothesized that Mediterranean-type *G6PD* genotype is rare (<0.2%, i.e., less than 1 in 500 based on the present study) or possibly absent in white Moors despite the supposed ethnic and historical ties with peoples from North Africa. It has been estimated that the prevalence rates of G6PD deficiency in the general population in Algeria and Morocco, two neighboring Maghreb countries north of Mauritania, are 2.8% and 0.6%, respectively [? ]. In Algeria, an early study had shown that 3.2% of males living in Algiers, the capital city, have G6PD deficiency [? ]. A later study conducted in Algerian males with clinical manifestations associated with G6PD deficiency (favism, neonatal jaundice, acute hemolytic anemia) demonstrated that 23% (23/100) were due to the Mediterranean-type *G6PD* genotype, while 46% had the African-type A- genotype [? ]. In a study performed in Morocco, 15% (54/360) of anemic patients were found to be G6PD-deficient, of whom 50% were characterized as the African-type G6PD A- and 7.4% had the Mediterranean-type variant [? ]. By contrast, our study showed that all genotypes associated with G6PD
deficiency were of African type among Mauritanians. Another hypothesis that may explain the absence of the Mediterranean-type G6PD deficiency in Mauritanians in the present study is the possibility that white Moors may have had in reality less ties with peoples in North Africa and remained geographically isolated from them by the Sahara Desert for centuries. The white Moors are composed of several social strata (warriors or "hassani," religious leaders or "zawaya," vassals or "zenaga," artisans, and bards or "griots"). Some of these social strata are known to have intermarried with black African populations over the past centuries [? ], which could have diluted the prevalence of G6PD B- disorder, if ever it was present to begin with. These hypotheses on the genetic relatedness of Mauritanian white Moors and peoples dwelling in North Africa today could be further strengthened by performing genetic analysis of other codons that have been reported to undergo unique mutations in some Algerian and Moroccan subjects with G6PD deficiency, such as the Kabyle and Aures G6PD A- variants [???]. These considerations need further historical and comparative human genome analysis in order to fully understand why Mediterranean-type G6PD genotype was not found in Mauritanians in our study.

The data presented here support the other hypothesis of this study that some Mauritanians of black African origin may be affected by the African-type G6PD deficiency. Among men belonging to black African ancestry (i.e., black Moors and one of three minority ethnic groups), 18 men (18/204; 8.8%, 95% CI (5.3–13.6%)) had the A- genotype. This prevalence rate in hemizygous Mauritanian men of black African origin is comparable to that reported for Senegalese males belonging to Sereer ethnic group (10%) and Malian males belonging to the Bambara ethnic group (10.5%), but higher than that of Dogon (2.2%) and Malinké (4.0%) ethnic groups in Mali [???].

Although we found a much higher prevalence of G6PD deficiency in black Moors and three minority ethnic groups than in white Moors in Mauritania, in practice, it may be argued that many of these individuals of black African ancestry may not be confronted with the potential problem of anti-hypnozoite therapy with 8-aminoquinolines because they tend to be Duffy-negative, as in the majority of black Africans elsewhere in sub-Saharan Africa, and may be naturally "resistant" to *P. vivax* infections [???]. However, cases of *P. vivax* infections have been reported in recent years in Duffy-negative individuals in Africa, including in Mauritania [???]. At present, *P. vivax* infection in Duffy-negative individuals is generally thought to be more of an exception rather than the rule, and the underlying biochemical and molecular mechanisms are poorly understood [??]. Furthermore, many of the minority ethnic groups of black African ancestry reside in southern Mauritania, where the chances of being exposed to the risk of *P. vivax* malaria are slim [????]. Nonetheless, many Mauritanians living in the south often travel to the north, where *P. vivax* is endemic, and some migrate to the northern Saharan zone, often for economic reasons, especially to Nouakchott, which has become a melting pot of various ethnic groups. Mauritanians of all ethnic backgrounds are therefore exposed to the risk of *P. vivax* infection, albeit to a lesser extent in those with black African ancestry due to the relative protection provided by the absence of Duffy antigen on the surface of red blood cells which serves as the principal known receptor for entry of *P. vivax* merozoites into the cell [?].

G6PD A⁻ is generally associated with mild to moderate deficiency (class III, i.e., 10–60% residual activity, except for Santamaria, class II, i.e., less than 10% residual activity) and G6PD B⁻ with severe deficiency (class II) [???]. Based on the WHO classification, the data obtained in the present study suggest that Mauritanians with mutant G6PD genotype are likely to present phenotypically mild to moderate deficiency. However, this assumption is not supported by solid evidence, and WHO experts have recommended against the use of this old classification to make clinical decisions [?]. The results of the present study are in agreement with earlier studies conducted in other West African countries which showed that, among different mutations resulting in African-type G6PD deficiency, G6PD⁻²⁰² is most frequently observed [?]. In several previous studies, A542T, G680T, and T968C mutations were not observed in West Africa (Burkina Faso, Ghana, Mali, and Nigeria) [??].
[]. In Mauritania, A376G/A542T Santamaria and A376G/T968C Betica-Selma were found, but not A376G/G680T.

Regardless of the A- variant and the ethnic group to which an individual *P. vivax*-infected patient belongs, G6PD screening test or standard biochemical test is mandatory before treatment with a standard 14-day regimen of primaquine can be considered. If the person is G6PD normal male, the standard dose of primaquine can be given if a qualitative or quantitative G6PD test displays normal enzymatic activity (> 60% of normal). If the person is G6PD normal female, the standard dose of primaquine can be given only if a reliable quantitative G6PD test is available and confirms a normal G6PD activity [? ]. If a Mauritanian individual is G6PD-deficient, it is likely to be due to the mutant African-type *G6PD* genotype. The clinician would have to decide on a case-by-case basis whether a modified scheme of anti-hypnozoite therapy with primaquine (once weekly dose for 8 weeks) would be beneficial and outweigh possible risks [? ]. In some cases, as in heterozygous females and even G6PD normal females with normal qualitative G6PD screening test but without quantitative test confirmation, a precise measurement of G6PD enzyme activity is required before making a decision.

It has been estimated that the worldwide prevalence of G6PD deficiency is 5% and that the prevalence is highest in Africa (8.5%) [???]. In the present study, the prevalence of G6PD deficiency based on genotyping was 6.3% in males and 1.5% in homozygous females, or 3.6% if the prevalence of both genders is combined. Data obtained in the present work, in particular samples from Kobeni city, were slightly skewed towards ethnic groups of black African ancestry and were not representative of the proportions of different ethnic groups present in Kobeni. Moreover, the samples used in the present study were collected from febrile patients in health centers and hospitals, where persons with an adequate financial resource and access to health tend to visit more often. Data obtained in health structures do not reflect the situation in the general population. A population-based study will be required for a more accurate assessment of the distribution of *G6PD* genotypes in the country. In addition, in some study sites, especially in the oasis city of Atar, none of the black African ethnic groups were seen at the health center. These considerations probably led to sampling bias in our study, with white Moors (50.1%, including men and women) over-representing the general population. The actual proportions of the ethnic groups of black African ancestry in Mauritania today are not known. Nevertheless, their proportion (18.2% composed of 13.5% Pular, 2.4% Wolof, and 2.3% Soninke) observed in the present study remains most likely less than what is expected. These confounding factors may have considerably modified the representativity of different ethnic groups in the samples, which may, at least in part, explain the lower-than-expected prevalence of African-type G6PD deficiency observed in the present study. This limitation of the study is counterbalanced by the fact that samples were collected from febrile patients spontaneously presenting themselves to health centers and hospitals. A subset of these patients may have been infected by *P. vivax* and may have required primaquine therapy.

Our study has other limitations. The present study was focused on the well-known G6PD 202, 376, 542, 680, and 968 polymorphisms associated with African-type G6PD A⁻ and G6PD 563 associated with the Mediterranean-type G6PD B⁻. Other less common mutations were not analyzed. Although one of our initial hypotheses about the possible presence of the Mediterranean variant associated with moderate G6PD deficiency among white Moors turned out to be unconfirmed among the study participants, preconceived ideas should not preclude investigators and clinicians to rule out G6PD deficiency. An unexpected variant with severe phenotypic deficiency may sometimes be detected in Africa [? ]. Moreover, due to the lack of sophisticated equipment and pre-existing technical know-how, G6PD phenotype was not characterized using the gold standard method, which is spectrophotometry. A study correlating G6PD genotype and G6PD phenotype is being planned to address this gap in the current knowledge of G6PD epidemiology in the country.
# 4. Materials and Methods

## 4.1. Study Sites

Field studies were conducted in 2015–2018 in the following cities in Mauritania, from the north to the south, located in six of 12 provinces (locally called “Wilaya”) of Mauritania (Figure ??): Atar (20°31′00″ N 13°03′00″ W; 38,877 inhabitants, according to the latest available data from the National Statistic Office [? ], the regional capital of Adrar, located in northern Saharan region; Nouakchott (18°06′ N 15°57′ W; 958,399 inhabitants), the capital city of Mauritania where one-fourth of the population reside; Aleg (17°3′ N 13°55′ W; 111,512 inhabitants), the regional capital of Brakna in central Mauritania; Kiffa (16°37′ N 11°24′ W; 110,714 inhabitants), the regional capital of Assaba, southern Mauritania, located about 600 km from Nouakchott; Kobeni city (15°49′02″ N 9°24′39″ W; 92,690 inhabitants) in the province of Hodh Elgharbi, southeastern region of the country near the Mauritania-Malian border; and Rosso (16°30′46″ N 15°48′18″ W; 57,726 inhabitants), the regional capital of Trarza, located along the Senegal River in southwestern Mauritania along the border between Mauritania and Senegal. In total, approximately 63% of the total Mauritanian population reside in these cities [? ].

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Study sites in Mauritania. Nouakchott and Atar are situated in the Saharan zone. Kobeni, Kiffa, and Aleg are in the Sahelo-Saharan transition zone. Rosso is located in the Sahelian zone along the Senegal River.</figcaption></figure>

Blood samples were collected in the following public health structures: health center of Atar; Teyarett health center and Maternity and Infant Hospital Center (Centre Hospitalier Mère et Enfant) in Nouakchott; health center of Aleg; Department of Internal Medicine and Infectious Diseases, Kiffa Regional Hospital, Kiffa; health center of Kobeni; and
health center of Rosso. Except for Kiffa and Aleg (which is located midway between Nouakchott and Kiffa), these study sites were selected on the basis of previous series of field surveys carried out by our research teams to assess malaria epidemiology in three different settings defined by geoclimatic conditions, namely, the northern Saharan zone (Atar and Nouakchott), southern and southeastern Sahelo-Saharan transition zone (Kiffa and Kobeni; also Aleg), and southwestern Sahelian zone (Rosso) [?????]. In addition, Atar, Nouakchott, Rosso, and Kobeni are sentinel sites of the Mauritanian National Malaria Control Programme [??].

## 4.2. Collection of Blood Samples

Fingerprick capillary blood samples were collected as part of the series of epidemiological studies conducted in the country between 2015 and 2018 [????]. Patients were recruited if they were febrile or had a recent history of fever within 48 h prior to consultation, regardless of age or sex. After explaining the purpose of the study and obtaining an informed consent, two to three drops (100-150 µL) of fingerprick capillary blood were imbibed onto Whatman grade 3MM filter paper (GE Healthcare UK Ltd., Little Chalfont, Buckinghamshire, UK) or Whatman FTA card (GE Healthcare), air dried, and stored in a sealed plastic sachet with a desiccant at -20 °C until DNA extraction.

## 4.3. Genotyping

A 1 mm-diameter disc of filter paper imbibed with dried blood was punched out and placed in 96-well plates. Human DNA was extracted using an automated MagMAX™-Express system (Thermo Fisher Scientific, Montigny-le-Bretonneux, France) according to the manufacturer's instructions. At the end of the DNA extraction protocol, approximately 120 µL of supernatant containing DNA were obtained from each sample.

In the present study, DNA fragments including the following G6PD nucleotide positions, where mutations have been reported to occur most frequently in Africa, giving rise to G6PD A- variant, were amplified by PCR and analyzed by restriction fragment length polymorphism (RFLP) and/or DNA sequencing in a stepwise manner: 376, 202, 542, 680, and 968 [? ] (Figure ??). In addition, nucleotide 563 was characterized in all samples to determine the presence of C563T mutation, which gives rise to the Mediterranean variant (G6PD B-).

Genomic human DNA was first amplified using 2.5-5 µl of eluted DNA template and primer pairs A376G_F (5'-CCCAGGCCACCCCAGAGGAGA-3') and A376G_R (5'-CGGCCCCGGACACGCTCATAG-3'; 0.5 µM each) in a PCR master mix containing Dream-Taq DNA polymerase, buffer, 2 mM MgCl₂, and deoxyribose nucleoside triphosphate (dNTP) (Thermo Fisher DreamTaq™ Green PCR Master mix; ThermoFisher Scientific, Illkirch, France) and the following thermal cycling program: initial step of 95 °C for 5 min, 35 cycles of 95 °C for 1 min, 53 °C (annealing temperature) for 1 min, and 72 °C for 1 min, then an extension at 72 °C for 10 min. PCR products (308 base pairs (bp)) were incubated at 37 °C overnight (at least 18 h) with *Fok* I (also called *Bse*G I or *Bts*C I) endonuclease (ThermoFisher Scientific, Illkirch, France).

PCR amplification was also performed for all samples to determine nucleotide 563 using the primer pairs C563T_F (5'-ACTCCCCAAGAGGGGTTCAAGG-3') and C563T_R (5'-CCAGCCTCCCAGGAGAGAGGAAG-3') (0.5 µM each) [? ]. The PCR program was as follows: initial step of 95 °C for 5 min, 35 cycles of 95 °C for 1 min, 63 °C (annealing temperature) for 1 min, and 72 °C for 1 min, then an extension at 72 °C for 10 min. PCR products (547 bp) were incubated at 37 °C overnight (at least 18 h) with *Mbo* II endonuclease (ThermoFisher Scientific, Illkirch, France).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Flow diagram depicting PCR-sequencing strategy to analyze G6PD genotypes. * If a sample was characterized as wild-type in nucleotides 376 and 563, it was considered as wild-type G6PD genotype in this study. ** The PCR-amplified fragment of exons 6-7 was sequenced in samples with A376G mutation if digestion did not display G202A or T968C mutation. *** For samples with A376G, but without G202A or T968C mutation, sequencing further confirmed the result of PCR-RFLP. † African-type G6PD A denotes the presence of A376G mutation, without a second mutation in nucleotide 202, 542, 680, or 968. The single mutation in 376 does not change the G6PD phenotype (i.e., the individual has a normal G6PD activity).</figcaption></figure>

Samples with A376G nucleotide substitution were further analyzed to detect the presence of an additional mutation in 202, 542, 680, or 968. The 216-bp fragment spanning nucleotide 202 was amplified using the primer pairs G202A_F (5'-CCACCACTGCCCCTGTGA CCT-3') and G202A_R (5'-GGCCCTGACACCACCCA CCTT-3') in DreamTaq™ Green PCR Master mix [?]. The thermal cycler was programmed as follows: initial denaturation at 95 °C for 5 min, 35 cycles of 95 °C for 1 min, 65 °C (annealing temperature) for 1 min, and 72 °C for 1 min, and an extension at 72 °C for 10 min. The PCR product was incubated at 37 °C overnight with *Nla* III (also called *Hin* 1 II) (ThermoFisher Scientific, Illkirch, France).

If nucleotide 202 was wild-type, PCR was performed to amplify a 283-bp fragment spanning nucleotide 968 using the primer pairs T968C_F (5'-TCCCTGCACCCCAACTCAAC-3') and T968C_R (5'-CCAGTTCTGCCTTGCTGGGC-3') [?]. PCR program was as for other fragments, with the anneal temperature of 65 °C. The amplicons were digested by incubat-
ing at 37 °C overnight with Nci I (ThermoFisher Scientific, Illkirch, France). After enzymatic digestion, digested amplificons were analyzed by 2% agarose gel electrophoresis. The PCR conditions and primers are summarized in Table ??.

<table><caption>**Table 5.** PCR primers used to amplify fragments of G6PD gene.</caption><thead><tr><th>Primers</th><th>Primer Sequence</th><th>Tm (°C)</th><th>Size (bp)</th><th>Enzyme</th></tr></thead><tbody><tr><td colspan="5"><strong>Exon 5</strong></td></tr><tr><td>A376G_F</td><td>5'-CCCAGGCCACCCCAGAGGAGA-3'</td><td rowspan="2">58</td><td rowspan="2">308</td><td rowspan="2">Fok I</td></tr><tr><td>A376G_R</td><td>5'-CGGCCCCGGACACGCTCATAG-3'</td></tr><tr><td colspan="5"><strong>Exons 6-7</strong></td></tr><tr><td>C563T_F</td><td>5'-ACTCCCCGAAGAGGGGTCAAGG-3'</td><td rowspan="2">62</td><td rowspan="2">547</td><td rowspan="2">Mbo II (C563T)<br>Acc III (G542T)<br>Bst NI (G680T)</td></tr><tr><td>C563T_R</td><td>5'-CCAGCCTCCCAGGAGAGAGGAAG-3'</td></tr><tr><td colspan="5"><strong>Exon 4</strong></td></tr><tr><td>G202A_F</td><td>5'-CCACCACTGCCCCTGTGACCT-3'</td><td rowspan="2">65</td><td rowspan="2">216</td><td rowspan="2">Nla III</td></tr><tr><td>G202A_R</td><td>5'-GGCCCTGACACCACCCACCTT-3'</td></tr><tr><td colspan="5"><strong>Exon 9</strong></td></tr><tr><td>T968C_F</td><td>5'-TCCCTGCACCCCAACTCAAC-3'</td><td rowspan="2">65</td><td rowspan="2">283</td><td rowspan="2">Nci I</td></tr><tr><td>T968C_R</td><td>5'-CCAGTTCTGCCTTGCTGGGC-3'</td></tr></tbody></table>

F for forward primer; R for reverse primer; Tm, annealing temperature. The size in base pairs (bp) refers to the expected fragment size of the PCR product. Three restriction sites were analyzed in the fragment spanning exons 6 and 7. The primer sequences of the fragments amplified from exons 4, 5, and 9 (G202A, A376G, T968C) were published by Carter et al. [?], and those in exons 6 and 7 (A542T, C563T, G680T) were published by Ezz El-Deen et al. [?]. The published primer sequence of C563T_F has one missing C. This was corrected based on the complete human G6PD gene sequence (GenBank accession number X55448.1) [?].

The 547-bp fragment spanning nucleotides 542, 563, and 680, which was initially amplified for RFLP using Mbo II to determine the nucleotide 563 in all samples, was sequenced if samples with A376G mutation had neither G202A nor T968C mutation. Sequencing was outsourced (Macrogen Europe; Amsterdam, The Netherlands).

F for forward primer; R for reverse primer; Tm, annealing temperature. The size in base pairs (bp) refers to the expected fragment size of the PCR product. Three restriction sites were analyzed in the fragment spanning exons 6 and 7. The primer sequences of the fragments amplified from exons 4, 5, and 9 (G202A, A376G, T968C) were published by Carter et al. [?], and those in exons 6 and 7 (A542T, C563T, G680T) were published by Ezz El-Deen et al. [?]. The published primer sequence of C563T_F has one missing C. This was corrected based on the complete human G6PD gene sequence (GenBank accession number X55448.1) [?].

### 4.4. Data Analysis

In all study sites, samples were randomly selected for analysis, with the exception of patients enrolled in Kobeni where, due to a large number of available samples (more than 2300) [?], 327 samples were selected based on the following arbitrary criteria: absence of malaria parasites and preference for patients belonging to one of the minority groups of black African origin. These latter groups are therefore over-represented in Kobeni in the present study.

The wild-type phenotype of normal enzyme is designated G6PD B. The corresponding genotype in hemizygous male is “B” for haplotype GAACGT (based on nucleotides G202/A376/A542/C563/G680/T968) and “BB” for homozygous normal female. The African variant genotype with a single A376G nucleotide change is associated with the phenotype G6PD A, which is not associated with any change in enzyme activity compared to the normal enzyme. The corresponding genotypes are denoted “A” in hemizygous male (haplotype GGACGT) and AA or AB in female. The genotypes denoted A and B in males and AA, AB, and BB in females are considered G6PD normal. The African-type deficient
enzyme is designated G6PD A-; there are two mutations, one of which is A376G. The corresponding genotypes are “A-” in hemizygous male (haplotypes **AGACGT**, **GGTCGT**, **GGACTT**, or **GGACGC**) and “A-A-” in homozygous deficient female. The genotypes of heterozygous females are either “AA-” or “BA-.” The Mediterranean variant with C563T (haplotype **GAATGT**) is designated Mediterranean G6PD B-.

If A376G nucleotide change (AAT to GAT; amino acid substitution Asn126Asp, exon 5) occurs, the genotype is either G6PD A or G6PD A-. Fok I yields two bands (116 bp + 192 bp) in the presence of the mutant 376G allele, but it has no restriction site in the PCR fragment with the wild-type A376 allele (a single 308-bp band). A heterozygote female with G6PD A or A- genotype is characterized by the presence of three bands (116 bp + 192 bp + 308 bp).

Enzymatic digestion of the fragment spanning nucleotide 202 with *Nla* III yields two fragments (81 bp + 135 bp) if it carries the mutant allele 202A, one uncut 216-bp fragment if it has the wild-type allele G202, or three fragments (81 bp + 135 bp + 216 bp) in heterozygous females. The genotype with double mutations A376G and G202A is designated “G6PD A<sup>-</sup>(202).”

*Nci* I restriction enzyme digests the 283-bp fragment into two (121 bp + 162 bp) fragments in the presence of mutation T968C or into three fragments (121 bp + 162 bp + 283 bp) in heterozygous females. The genotype characterized by double mutations A376G and T968C is called “G6PD A<sup>-</sup> Betica-Selma.”

Sequencing of some PCR products of samples with A376G mutation, but without G202A or T968C, allowed not only to determine A542T and G680T mutations but also to confirm PCR-RFLP results using *Mbo* II to determine nucleotide 563. The genotypes with double mutations A376G and A542T or A376G and G680T are called “G6PD A<sup>-</sup> Santamaria” and “G6PD A<sup>-</sup>(680),” respectively [?].

*Mbo* II enzyme digests the 547-bp fragment into 4 fragments (25 + 26 + 199 + 377 bp) if there is a wild-type C563 allele and 5 fragments (25 + 26 + 100 + 199 + 277 bp) in the presence of mutant 563T (Mediterranean-type G6PD B-).

In the present paper, the term “normal G6PD genotype” is synonymous with the wild-type G6PD B and denotes the absence of mutation in one of the following nucleotide positions: 202, 376, 542, 563, 680, and/or 968.

Proportions were compared using a two-tailed Fisher's exact test. The *p*-value of <0.05 was used to test for significance.

**5. Conclusions**

This is the first nationwide survey of G6PD genotyping that addresses the relation between ethnic origin and G6PD deficiency in Mauritania. The African-type G6PD mutations were found in all ethnic groups, more frequently in those of black African ancestry than in white Moors. From a practical viewpoint, every *P. vivax*-infected Mauritanian patient would require a reliable G6PD screening test and, for females, a reliable quantitative G6PD measurement, before a combination of blood schizonticidal and anti-hypnozoite antimalarial drugs can be administered to ensure that radical cure is achieved to completely eliminate the parasites from the patient. A phenotypic and/or genotypic surveillance of G6PD deficiency is warranted as the country attempts to move on to implement anti-hypnozoite therapy as part of the malaria elimination program.

**Author Contributions:** Conceptualization, M.S.O.A.S. and A.O.M.S.B.; field investigation, O.K.M.D., S.M.D., B.M.B., and M.A.B.; molecular methods, O.K.M.D., A.G., N.P.M., S.B., and L.B.; validation and formal analysis, S.B.; data curation, S.B.; writing—original draft preparation, A.O.M.S.B. and S.B.; writing—review and editing, L.B.; field supervision, M.S.O.A.S.; laboratory supervision, N.P.M. and S.B.; project administration, A.O.M.S.B. and L.B.; funding acquisition, A.O.M.S.B. and L.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Institut de Recherche pour le Développement (IRD), Jeune Equipe Associée à l'IRD (RI3M), Institut de Recherche Biomédicale des Armées (IRBA), and Direction Générale de l'Armement (DGA). O.K.M.D. received a scholarship from the French government. The
APC was funded by Institut de Recherche Biomédicale des Armées (IRBA) and Direction Générale de l'Armement (DGA).

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki and was reviewed and approved by the institutional ethics committee of the Université de Nouakchott Al-Aasriya, Nouakchott, Mauritania (approval no. 112/12-09-2014/USTM, 003/2020/CE/UNA) and the Institutional Review Board of the Institut de Recherche pour le Développement (IRD), Marseille, France (Comité consultative de déontologie et d'éthique approval no. 15/12/2012).

**Informed Consent Statement:** The purpose of the study was explained in local dialect to all adult patients or to the parents or legal guardians of pediatric patients. Informed consent was obtained from all adult subjects or either the parents or legal guardians of children involved in the study.

**Data Availability Statement:** The authors confirm that the data supporting the findings of this study are available within the article.

**Acknowledgments:** We thank the physicians, technicians, and nurses at health structures in the six study sites for facilitating the study.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Beutler, E. Glucose-6-phosphate dehydrogenase deficiency. *Br. J. Haematol.* **1970**, *18*, 117-121. [CrossRef]

2. Vulliamy, T.; Mason, P.; Luzzatto, L. The molecular basis of glucose-6-phosphate dehydrogenase deficiency. *Trends Genet.* **1992**, *8*, 138-143. [CrossRef]

3. Nkhoma, E.T.; Poole, C.; Vannappagari, V.; Hall, S.A.; Beutler, E. The global prevalence of glucose-6-phosphate dehydrogenase deficiency: A systematic review and meta-analysis. *Blood Cells Mol. Dis.* **2009**, *42*, 267-278. [CrossRef] [PubMed]

4. Howes, R.E.; Battle, K.E.; Satyagraha, A.W.; Baird, J.K.; Hay, S.I. G6PD deficiency: Global distribution, genetic variants and primaquine therapy. *Adv. Parasitol.* **2013**, *81*, 133-201. [CrossRef] [PubMed]

5. Luzzatto, L.; Arese, P. Favism and glucose-6-phosphate dehydrogenase deficiency. *N. Engl. J. Med.* **2018**, *378*, 60-71. [CrossRef] [PubMed]

6. Lippi, G.; Mattiuzzi, C. Updated worldwide epidemiology of inherited erythrocyte disorders. *Acta Haematol.* **2020**, *143*, 196-203. [CrossRef] [PubMed]

7. Persico, M.G.; Viglietto, G.; Martini, G.; Toniolo, D.; Paonessa, G.; Moscatelli, C.; Dono, R.; Vulliamy, T.; Luzzatto, L.; D'Urso, M. Isolation of human glucose-6-phosphate dehydrogenase (G6PD) cDNA clones: Primary structure of the protein and unusual 5' non-coding region. *Nucleic Acids Res.* **1986**, *14*, 2511-2522. [CrossRef] [PubMed]

8. Minucci, A.; Moradkhani, K.; Hwang, M.J.; Zuppi, C.; Giardina, B.; Capoluongo, E. Glucose-6-phosphate dehydrogenase (G6PD) mutations database: Review of the "old" and update of the new mutations. *Blood Cells Mol. Dis.* **2012**, *48*, 154-165. [CrossRef]

9. Gómez-Manzo, S.; Marcial-Quino, J.; Vanoye-Carlo, A.; Serrano-Posada, H.; Ortega-Cuellar, D.; González-Valdez, A.; Castillo-Rodríguez, R.A.; Hernández-Ochoa, B.; Sierra-Palacios, E.; Rodríguez-Bustamante, E.; et al. Glucose-6-phosphate dehydrogenase: Update and analysis of new mutations around the world. *Int. J. Mol. Sci.* **2016**, *17*, 2069. [CrossRef]

10. Collins, W.E.; Jeffery, G.M. *Plasmodium ovale*: Parasite and disease. *Clin. Microbiol. Rev.* **2005**, *18*, 570-581. [CrossRef]

11. White, N.J.; Imwong, M. Relapse. *Adv. Parasitol.* **2012**, *80*, 113-150. [CrossRef]

12. WHO. *Guidelines for the Treatment of Malaria*, 3rd ed.; World Health Organization: Geneva, Switzerland, 2015; pp. 1-313.

13. Lekweiry, K.M.; Salem, M.S.; Basco, L.K.; Briolant, S.; Hafid, J.; Boukhary, A.O. Malaria in Mauritania: Retrospective and prospective overview. *Malar. J.* **2015**, *14*, 100. [CrossRef]

14. Ould Ahmedou Salem, M.S.; Mint Lekweiry, K.; Mint Deida, J.; Ould Emouh, A.; Ould Weddady, M.; Ould Mohamed Salem Boukhary, A.; Basco, L.K. Increasing prevalence of *Plasmodium vivax* among febrile patients in Nouakchott, Mauritania. *Am. J. Trop. Med. Hyg.* **2015**, *92*, 537-540. [CrossRef]

15. Mint Deida, J.; Ould Khalef, Y.; Mint Semane, E.; Ould Ahmedou Salem, M.S.; Bogreau, H.; Basco, L.; Ould Mohamed Salem Boukhary, A.; Tahar, R. Assessment of drug resistance associated genetic diversity in Mauritanian isolates of *Plasmodium vivax* reveals limited polymorphism. *Malar. J.* **2018**, *17*, 416. [CrossRef] [PubMed]

16. Deida, J.; Tahar, R.; Khalef, Y.O.; Lekweiry, K.M.; Hmeyade, A.; Khairy, M.L.O.; Simard, F.; Bogreau, H.; Basco, L.; Boukhary, A.O.M.S. Oasis malaria, northern Mauritania. *Emerg. Infect. Dis.* **2019**, *25*, 273-280. [CrossRef]

17. Mauritanian Ministry of Health. *Guide Clinique et Thérapeutique à L'usage du Personnel des Centres de Santé de la Mauritanie*, 2nd ed.; Agencia Española de Cooperación Internacional para el Desarrollo (AECID): Madrid, Spain, 2013; pp. 1-244. Available online: http://www.sante.gov.mr/?wpfb_dl=140 (accessed on 22 January 2021).
18. Djigo, O.K.M.; Bollahi, M.A.; Hasni Ebou, M.; Ould Ahmedou Salem, M.S.; Tahar, R.; Bogreau, H.; Basco, L.; Ould Mohamed Salem Boukhary, A. Assessment of glucose-6-phosphate dehydrogenase activity using CareStart G6PD rapid diagnostic test and associated genetic variants in *Plasmodium vivax* malaria endemic setting in Mauritania. *PLoS ONE* **2019**, *14*, e0220977. [CrossRef] [PubMed]

19. Mohamed, G.S.; Lemine, S.M.; Cheibetta, S.; Mohamed, A. Neonatal screening for glucose-6-phosphate dehydrogenase (G6PD) deficiency in Mauritania. *Pan Afr. Med. J.* **2018**, *30*, 224. [CrossRef]

20. WHO Working Group. Glucose-6-phosphate dehydrogenase deficiency. *Bull. World Health Organ.* **1989**, *67*, 601–611.

21. Handloff, R.E. (Ed.) *Mauritania: A Country Study*, 2nd ed.; Federal Research Division; U.S. Government Printing Office for the Library of Congress: Washington, DC, USA, 1990; pp. 1–218.

22. Nafa, K.; Reghis, A.; Osmani, N.; Baghli, L.; Aït-Abbes, H.; Benabadji, M.; Kaplan, J.C.; Vulliamy, T.; Luzzatto, L. At least five polymorphic mutants account for the prevalence of glucose-6-phosphate dehydrogenase deficiency in Algeria. *Hum. Genet.* **1994**, *94*, 513–517. [CrossRef] [PubMed]

23. Doss, C.G.; Alasmar, D.R.; Bux, R.I.; Sneha, P.; Bakhsh, F.D.; Al-Azwani, I.; Bekay, R.E.; Zayed, H. Genetic epidemiology of glucose-6-phosphate dehydrogenase deficiency in the Arab world. *Sci. Rep.* **2016**, *6*, 37284. [CrossRef]

24. Migot-Nabias, F.; Pelleau, S.; Watier, L.; Guitard, J.; Toly, C.; De Araujo, C.; Ngom, M.I.; Chevillard, C.; Gaye, O.; Garcia, A. Red blood cell polymorphisms in relation to *Plasmodium falciparum* asymptomatic parasite densities and morbidity in Senegal. *Microbes Infect.* **2006**, *8*, 2352–2358. [CrossRef] [PubMed]

25. Carter, N.; Pamba, A.; Duparc, S.; Waitumbi, J.N. Frequency of glucose-6-phosphate dehydrogenase deficiency in malaria patients from six African countries enrolled in two randomized anti-malarial clinical trials. *Malar. J.* **2011**, *10*, 241. [CrossRef] [PubMed]

26. Tine, R.C.K.; Ndiaye, M.; Hansson, H.H.; Ndour, C.T.; Faye, B.; Alifrangis, M.; Sylla, K.; Ndiaye, J.L.; Magnussen, P.; Bygbjerg, I.C.; et al. The association between malaria parasitaemia, erythrocyte polymorphisms, malnutrition and anaemia in children less than 10 years in Senegal: A case control study. *BMC Res. Notes* **2012**, *5*, 565. [CrossRef] [PubMed]

27. Maiga, B.; Dolo, A.; Campino, S.; Sepulveda, N.; Corran, P.; Rockett, K.A.; Troye-Blomberg, M.; Doumbo, O.K.; Clark, T.G. Glucose-6-phosphate dehydrogenase polymorphisms and susceptibility to mild malaria in Dogon and Fulani, Mali. *Malar. J.* **2014**, *13*, 270. [CrossRef]

28. Benyachou, B. L'anémie par Déficit en G6PD Chez L'enfant. MD Thesis, Faculté de Médecine et de Pharmacie de Fès, Fès, Morocco, 2008. Available online: http://scolarite.fmp-usmba.ac.ma/cdim/mediatheque/e_theses/4-08.pdf (accessed on 30 January 2021).

29. Benabadji, M.; Merad, F.; Benmoussa, M.; Trabuchet, G.; Junien, C.; Dreyfus, J.C.; Kaplan, J.C. Heterogeneity of glucose-6-phosphate dehydrogenase deficiency in Algeria. Study in Northern Algeria with description of five new variants. *Hum. Genet.* **1978**, *40*, 177–184. [CrossRef]

30. Saoud, M.Z.; Benkirane, S.; Bennis, F.Z.; Mechal, Y.; Mamad, H.; Dahmani, F.; Berchane, Z.; Woumki, A.; Masrar, A. Déficit en glucose-6-phosphate déshydrogénase: À propos de 54 cas. *Hématologie* **2019**, *25*, 16. [CrossRef]

31. 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*, 262–263.

32. Dolo, A.; Maiga, B.; Guindo, A.; Diakité, S.A.; Diakite, M.; Tapily, A.; Traoré, M.; Sangaré, B.; Arama, C.; Daou, M.; et al. Frequency of glucose-6-phosphate dehydrogenase deficiency (A-376/202) in three Malian ethnic groups. *Bull. Soc. Pathol. Exot. 2014*, **107**, 165–170. [CrossRef]

33. Hadley, T.J. Invasion of erythrocytes by malaria parasites: A cellular and molecular overview. *Annu. Rev. Microbiol.* **1986**, *40*, 451–477. [CrossRef]

34. Langhi, D.M., Jr.; Bordin, J.O. Duffy blood group and malaria. *Hematology* **2006**, *11*, 389–398. [CrossRef] [PubMed]

35. Rosenberg, R. *Plasmodium vivax* in Africa: Hidden in plain sight? *Trends Parasitol.* **2007**, *23*, 193–196. [CrossRef] [PubMed]

36. Wurtz, N.; Mint Lekweiry, K.; Bogreau, H.; Pradines, B.; Rogier, C.; Ould Mohamed Salem Boukhary, A.; Hafid, J.E.; Ould Ahmedou Salem, M.S.; Trape, J.F.; Basco, L.K.; et al. Vivax malaria in Mauritania includes infection of a Duffy-negative individual. *Malar. J.* **2011**, *10*, 336. [CrossRef]

37. Culleton, R.; Carter, R. African *Plasmodium vivax*: Distribution and origins. *Int. J. Parasitol.* **2012**, *42*, 1091–1097. [CrossRef] [PubMed]

38. Golassa, L.; Amenga-Etego, L.; Lo, E.; Amambua-Ngwa, A. The biology of unconventional invasion of Duffy-negative reticulo-cytes by *Plasmodium vivax* and its implication in malaria epidemiology and public health. *Malar. J.* **2020**, *19*, 299. [CrossRef]

39. Popovici, J.; Roesch, C.; Rougeron, V. The enigmatic mechanisms by which *Plasmodium vivax* infects Duffy-negative individuals. *PLoS Pathog.* **2020**, *16*, e1008258. [CrossRef]

40. Ouldabdallahi, M.; Ouldbezeid, M.; Dieye, M.; Yacine, B.; Faye, O. Study of the incidence of malaria in febrile patients and among schoolchildren in the region of Trarza, Islamic Republic of Mauritania. *Bull. Soc. Pathol. Exot. 2011*, *104*, 288–290. [CrossRef]

41. Ould Ahmedou Salem, M.S.; Basco, L.K.; Ouldabdallahi, M.; Mint Lekweiry, K.; Konaté, L.; Faye, O.; Ould Mohamed Salem Boukhary, A. Malaria-associated morbidity during the rainy season in Saharan and Sahelian zones in Mauritania. *Acta Trop.* **2015**, *152*, 1–7. [CrossRef] [PubMed]

42. Ouldabdallahi Moukah, M.; Ba, O.; Ba, H.; Ould Khairy, M.L.; Faye, O.; Bogreau, H.; Simard, F.; Basco, L.K. Malaria in three epidemiological strata in Mauritania. *Malar. J.* **2016**, *15*, 204. [CrossRef]
43. Ba, H.; Duffy, C.W.; Ahouidi, A.D.; Deh, Y.B.; Diallo, M.Y.; Tandia, A.; Conway, D.J. Widespread distribution of *Plasmodium vivax* malaria in Mauritania on the interface of the Maghreb and West Africa. *Malar. J.* **2016**, *15*, 80. [CrossRef]

44. Ba, H.; Auburn, S.; Jacob, C.G.; Goncalves, S.; Duffy, C.W.; Stewart, L.B.; Price, R.N.; Deh, Y.B.; Diallo, M.Y.; Tandia, A.; et al. Multi-locus genotyping reveals established endemicity of a geographically distinct *Plasmodium vivax* population in Mauritania, West Africa. *PLoS Negl. Trop. Dis.* **2020**, *14*, e0008945. [CrossRef]

45. Diallo, S.M.; Bogreau, H.; Papa Mze, N.; Ould Ahmedou Salem, M.S.; Ould Khairy, M.L.; Parola, P.; Basco, L.; Ould Mohamed Salem Boukhary, A. Malaria epidemiology in Kobeni department, southeastern Mauritania from 2015 to 2017. *Infect. Dis. Poverty* **2020**, *9*, 21. [CrossRef]

46. Yoshida, A.; Beutler, E.; Motulsky, A.G. Human glucose-6-phosphate dehydrogenase variants. *Bull. World Health Organ.* **1971**, *45*, 243–253. [PubMed]

47. World Health Organization. Malaria Policy Advisory Committee Meeting. Point-of-care G6PD testing to support safe use of primaquine for the treatment of vivax malaria. In *WHO Evidence Review Group Meeting Report, 8–9 October 2014*; WHO/HTM/GMP/MPAC/2015.6; World Health Organization: Geneva, Switzerland, 2015.

48. Howes, R.E.; Chan, E.R.; Rakotomanga, T.A.; Schulte, S.; Gibson, J.; Zikursh, M.; Franchard, T.; Ramiranirina, B.; Ratsimbasoa, A.; Zimmerman, P.A. Prevalence and genetic variants of G6PD deficiency among two Malagasy populations living in *Plasmodium vivax*-endemic areas. *Malar. J.* **2017**, *16*, 139. [CrossRef] [PubMed]

49. Office National de la Statistique. Annuaire statistique 2019. *Ministère des Affaires Economiques et de Promotion des Secteurs Productifs. November 2020*. Available online: http://www.ons.mr/images/Archive/doc/publication/Annuaire_Statistique_2019.pdf (accessed on 13 March 2021).

50. Boushab, M.B.; Ould Ahmedou Salem, M.S.; Ould Mohamed Salem Boukhary, A.; Parola, P.; Basco, L. Clinical features and mortality associated with severe malaria in adults in southern Mauritania. *Trop. Med. Infect. Dis.* **2021**, *22*, 1. [CrossRef]

51. Ouldabdallahi, M.; Alew, I.; Salem, M.S.; Dit Dialaw Ba, M.; Boukhary, A.O.; Khairy, M.L.; Aziz, M.B.; Ringwald, P.; Basco, L.K.; Niang, S.D.; et al. Efficacy of artesunate-amodiaquine for the treatment of acute uncomplicated falciparum malaria in southern Mauritania. *Malar. J.* **2014**, *13*, 496. [CrossRef]

52. Ould Ahmedou Salem, M.S.; Mohamed Lemine, Y.O.; Deida, J.M.; Lemrabott, M.A.; Ouldabdallahi, M.; Ba, M.D.; Boukhary, A.O.; Khairy, M.L.; Abdel Aziz, M.B.; Ringwald, P.; et al. Efficacy of chloroquine for the treatment of *Plasmodium vivax* in the Saharan zone in Mauritania. *Malar. J.* **2015**, *14*, 39. [CrossRef]

53. Ezz El-Deen, Z.M.; Hussin, N.F.; Abdel Hamid, T.A.; Abdel Migeed, O.R.; Samy, R.M. G6PD deficiency and G6PD (Mediterranean and silent) polymorphisms in Egyptian infants with neonatal hyperbilirubinemia. *Lab. Med.* **2013**, *44*, 228–234. [CrossRef]

54. Chen, E.Y.; Cheng, A.; Lee, A.; Kuang, W.J.; Hillier, L.; Green, P.; Schlessinger, D.; Ciccodicola, A.; D'Urso, M. Sequence of human glucose-6-phosphate dehydrogenase cloned in plasmids and a yeast artificial chromosome. *Genomics* **1991**, *10*, 792–800. [CrossRef]