# Potential hidden Plasmodium vivax malaria reservoirs from low parasitemia Duffy-negative Ethiopians: Molecular evidence

Abnet Abebe, Isabelle Bouyssou, Solenne Mabilotte, Sisay Dugassa, Ashenafi Assefa, Jonathan Juliano, Eugenia Lo, Didier Menard, Lemu Golassa

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

## To cite this version:

Abnet Abebe, Isabelle Bouyssou, Solenne Mabilotte, Sisay Dugassa, Ashenafi Assefa, et al.. Potential hidden Plasmodium vivax malaria reservoirs from low parasitemia Duffy-negative Ethiopians: Molecular evidence. PLoS Neglected Tropical Diseases, 2023, 17 (7), pp.e0011326. 10.1371/journal.pntd.0011326. hal-04224092

HAL Id: hal-04224092

https://hal.science/hal-04224092v1

Submitted on 19 Oct 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>
RESEARCH ARTICLE

# Potential hidden *Plasmodium vivax* malaria reservoirs from low parasitemia Duffy-negative Ethiopians: Molecular evidence

Abnet Abebe<sup>1,2*</sup>, Isabelle Bouyssou<sup>3,4</sup>, Solenne Mabilotte<sup>5</sup>, Sisay Dugassa<sup>1</sup>, Ashenafi Assefa<sup>2,6</sup>, Jonathan J. Juliano<sup>6</sup>, Eugenia Lo<sup>7</sup>, Didier Menard<sup>3,5,8,9</sup>, Lemu Golassa<sup>1*</sup>

1 Aklilu Lemma Institute of Pathobiology, Addis Ababa University, Ethiopia, 2 Ethiopian Public Health Institute, Addis Ababa, Ethiopia, 3 Institut Pasteur, Université Paris Cité, Malaria Genetics and Resistance Unit, INSERM U1201, F-75015 Paris, France, 4 Sorbonne Université, Collège Doctoral ED 515 Complexité du Vivant, Paris, France, 5 Université de Strasbourg, Institute of Parasitology and Tropical Diseases, Dynamics of Host-Pathogen Interactions, F-67000 Strasbourg, France, 6 Institute of Infectious Disease and Global Health, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America, 7 Department of Biological Sciences, Bioinformatics Research Center, University of North Carolina at Charlotte, United States of America, 8 Institut Pasteur, Université Paris Cité, Malaria Parasite Biology and Vaccines Unit, Paris, France, 9 CHU Strasbourg, Laboratory of Parasitology and Medical Mycology, Strasbourg, France

* abnetabas@gmail.com (AA); lgolassa@gmail.com (LG)

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

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

OPEN ACCESS

**Citation:** Abebe A, Bouyssou I, Mabilotte S, Dugassa S, Assefa A, Juliano JJ, et al. (2023) Potential hidden *Plasmodium vivax* malaria reservoirs from low parasitemia Duffy-negative Ethiopians: Molecular evidence. PLoS Negl Trop Dis 17(7): e0011326. https://doi.org/10.1371/journal.pntd.0011326

**Editor:** Kamala Thriemer, Menzies School of Health Research, AUSTRALIA

**Received:** April 24, 2023

**Accepted:** June 14, 2023

**Published:** July 3, 2023

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pntd.0011326

**Copyright:** This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

**Data Availability Statement:** All relevant data are within the manuscript.

## Abstract

### Background

The interaction between the *Plasmodium vivax* Duffy-binding protein and the corresponding Duffy Antigen Receptor for Chemokines (DARC) is primarily responsible for the invasion of reticulocytes by *P. vivax*. The Duffy-negative host phenotype, highly prevalent in sub-Saharan Africa, is caused by a single point mutation in the GATA-1 transcription factor binding site of the DARC gene promoter. The aim of this study was to assess the Duffy status of patients with *P. vivax* infection from different study sites in Ethiopia.

### Methods

A cross-sectional study was conducted from February 2021 to September 2022 at five varying eco-epidemiological malaria endemic sites in Ethiopia. Outpatients who were diagnosed with *P. vivax* infection (pure and mixed *P. vivax/P. falciparum*) by microscopy and Rapid Diagnostic Test (RDT) were subjected to PCR genotyping at the DARC promoter. The associations between *P. vivax* infection, host genotypes and other factors were evaluated.

### Result

In total, 361 patients with *P. vivax* infection were included in the study. Patients with pure *P. vivax* infections accounted for 89.8% (324/361), while the remaining 10.2% (37/361) had mixed *P. vivax/P. falciparum* infections. About 95.6% (345/361) of the participants were Duffy-positives (21.2% homozygous and 78.8%, heterozygous) and 4.4% (16/361) were Duffy-negatives. The mean asexual parasite density in homozygous and heterozygous
**Funding:** The author(s) received no specific funding for this work.

**Competing interests:** The authors have declared that no competing interests exist.

Duffy-positives was 12,165 p/µl (IQR25-75: 1,640–24,234 p/µl) and 11,655 p/µl (IQR25-75: 1,676–14,065 p/µl), respectively, significantly higher than that in Duffy-negatives (1,227 p/µl; IQR25-75: 539–1,732 p/µl).

## Conclusion

This study confirms that Duffy-negativity does not provide complete protection against *P. vivax* infection. The development of *P. vivax*-specific elimination strategies, including alternative antimalarial vaccines should be facilitated by a better understanding of the epidemiological landscape of *vivax* malaria in Africa. More importantly, low parasitemia associated with *P. vivax* infections in Duffy-negative patients may represent hidden reservoirs of transmission in Ethiopia.

## Author summary

*Plasmodium vivax* generally receives less attention than *P. falciparum* and was neglected in sub-Saharan Africa. However, the characteristics of *P. vivax* infection in Duffy-negative individuals, and the distribution of Duffy blood group in different eco-epidemiological zones and ethnic groups of Ethiopia are not well documented. Here, we determined the Duffy genotypes of *P. vivax* infected patients across broad regions of Ethiopia. It is clear that Duffy negative individuals in Ethiopia are not fully protected against *P. vivax* infection, and that these infections in Duffy negatives are often associated with low parasitemia. Our findings lend support to the notion that *P. vivax* may have developed a Duffy-independent erythrocyte invasion pathway and/or evolution in host immune responses.

## Introduction

In 2021, there were an estimated 247 million cases and 619,000 deaths from malaria worldwide. Of these, 2% of the cases were due to *P. vivax*, mainly in the Horn of Africa, South America, India and Asia [1]. *P. vivax* malaria is recognized as a cause of severe morbidity and mortality, with a significant negative impact on health in endemic countries [2–4]. The unique biological features of *P. vivax*, such as hypnozoite formation and early gametocytogenesis, pose a challenge to global malaria elimination efforts [5].

The *P. vivax* invasion pathway into reticulocytes was previously thought to be primarily dependent on the interaction between the Duffy binding protein (PvDBP) and its cognate receptor, the Duffy Antigen Receptor for Chemokines (DARC) [6–8]. Polymorphisms in the promoter of the DARC gene alter the expression of the Duffy antigen protein and determine an individual’s Duffy blood group [9]. A point mutation in the GATA-1 transcription factor binding site of the promoter of the gene at position–67 encodes the change of a T nucleotide to a C nucleotide and underlies the Duffy-negative (erythrocyte silent) phenotype [10–12]. Individuals with this mutation do not express the Duffy antigen on the surface of their erythrocytes [13].

Duffy-negative phenotype in the population occurs at low frequency in Caucasians (~3.5%) [14], and has a prevalence of almost 100% in West Africa, > 80% in African Americans [15], 20–36% in East Africa and 84% in southern Africa [16]. In Ethiopia, a study reported 94/416
(~23%) of the febrile patients from the six health centres/hospitals were Duffy-negatives, though only two of them were confirmed with *P. vivax* [17].

The scientific paradigm that *P. vivax* exclusively invades Duffy-positive reticulocytes through the interactions of PvDBP and *DARC* was established, and it was assumed that *vivax* malaria was rarely transmitted or even absent in regions where Duffy-negativity predominates, such as sub-Saharan Africa [18,19]. However, recent studies have shown that *P. vivax* can infect Duffy-negative reticulocytes [20], and is no longer a barrier to such infections [21]. In East Africa, the prevalence of *P. vivax* in Duffy-negative populations ranges from 3% [22] to ~12% in Ethiopia [16] and 3% [16] to ~18% in Sudan [9]. Although several studies have demonstrated the presence of *P. vivax* infection in Duffy-negative individuals, the distribution of these infections in different eco-epidemiological zones and ethnic groups are not well documented, especially in Ethiopia. This study compares the prevalence of Duffy blood group genotypes of *P. vivax* infected patients among different eco-epidemiological zones of Ethiopia.

# Materials and methods

## Ethical statement

Ethical approval was obtained from Aklilu Lema Institute of Pathobiology Institutional Review Board (IRB) (Ref No: ALIPB IRB/34/2013/20), Addis Ababa University, and ministry of education (MoE) via its national research ethical review committee (Ref No:7/2-514/m259/35), Addis Ababa, Ethiopia. An official letter was sent to the study sites and the study was initiated after obtaining permission from the respective health facilities. Written informed consent and/or assent was obtained from all study participants, using the local language where appropriate.

## Study design, period and areas

A health-facility based cross-sectional study was conducted from February 2021 to September 2022. Based on the annual parasite index (API), malaria risk stratification in Ethiopia is divided into five classes: i) malaria free, ii) very low, iii) low, iv) moderate and v) high malaria risk.

The study sites were selected from different malaria risk stratification areas. The sites include: Arbaminch, Dubti, Gambella, Metehara, and Shewarobit (Fig 1). Arbaminch hospital is located in the Southern Nation and Nationalities People Regional State in the southern part of the country (altitude 1200 m; Lat, Long: 6.02043, 37.56788); Dubti hospital is located in the Afar Regional State in the northern-east part of the country (altitude 379 m; Lat, Long: 11.72654, 41.09440); Gambella hospital is located in Gambella Regional State in the western part of the country (altitude 447 m; Lat, Long: 8.24810, 34.59071); Metehara hospital is located in Oromia Regional State in the eastern part of the country (altitude 959 m, Lat, Long: 8.89932, 39.91726); and Shewarobit hospital is located in Amhara Regional State in the northern part of the country (altitude 1268 m; Lat, Long: 10.01734, 39.91253).

## Sample collection

Blood samples diagnosed positive for *P. vivax* (pure or mixed *P. vivax*/*P. falciparum*) by qPCR from patients seen at health facilities were included to determine their Duffy status. The sample size was calculated based on the proportion of Duffy-negative (20%) cases previously reported in Ethiopia [23].
<figure><img src="image_5.png" /><figcaption>Fig 1. A map of the study sites with malaria risk stratification background, Ethiopia.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0011326.g001

Malaria diagnosis by Rapid Diagnostic Tests (RDT) and microscopy

Finger prick blood samples were collected from self-reported febrile patients seeking malaria diagnosis at a participating health centre. The malaria laboratory diagnosis was carried out for each suspected malaria patient using RDT and microscopy. SD BIOLINE Malaria Ag Pf/Pv (05FK80) RDT was used for the detection of HRP-II (Histidine-rich protein II) specific for *P. falciparum* and *Plasmodium* lactate dehydrogenase (pLDH) specific for *P. vivax*. Based on WHO public report, the panel detection score of the kit was 96.0% at 200 parasites/µl, and 95.0% at 200 parasites/µl for *P. falciparum*, and *P. vivax*, respectively [24]. Five microliters (5µl) of sample were obtained for the test and the result was interpreted after 15 minutes (up to 30 minutes) according to the manufacturer's guidelines.

Thick and thin blood films (stained with 10% Giemsa staining solution—pH 7.2—for 10 minutes) were prepared from finger-prick blood samples. Parasite densities were determined on thick films by expert malaria microscopists for both sexual and asexual stages of the parasites, considering a minimum of 200 WBCs, the reported number of parasites per microliter of
blood, assuming a total white blood cell count of 8000/µL. *P. vivax* sexual stages/gametocytes quantification was performed on the blood film with *P. vivax* (pure) and a mixed species (*P. vivax*/*P. falciparum*). Parasitemia was then classified as low (<1000 parasites per microliter of blood), intermediate (1000–4999 parasites per microliter of blood) and high (≥5000 parasites per microliter of blood) [25]. Slides were declared negative after examination of at least 100 high-power microscope fields [26]. Blood film slides were re-examined and quantified by WHO certified malaria microscopists. Patients who were positive for *P. vivax* (pure) and mixed *P. vivax*/*P. falciparum* infection by both RDT and microscopy were included in the study. Venous blood samples were collected and spotted onto Whatman 903 Protein Saver Card, US (Dried Blood Spot, DBS), for molecular analysis of the Duffy genotypes.

## DNA extraction and molecular diagnosis of *Plasmodium* species

Genomic DNA (human and parasite) was extracted from a 6 mm diameter spot of a DBS sample and eluted in 100µl of TE Buffer (Tris-EDTA) using the QIAamp DNA Extraction Kit (Cat.No: 79216, Lot: 172018338, Germany) according to the manufacturer's instructions for the molecular identification of malaria parasites and Duffy genotypes. The extracted DNA was stored at -20°C for further use.

Primers previously designed to target the mitochondrial species-specific *cytochrome b* (*cytb*) gene were used for the detection and identification of *P. vivax* and *P. falciparum* (*Table 1*). Given the multiple copies of the mitochondria, this assay allows the detection of low density *Plasmodium* infections in malaria pre-elimination settings [27].

The PCR amplifications were carried out using 10 µL of Advanced Mix Sybgreen (2x) (Catalog # 1725271, Bio-Rad laboratories Inc. United State), 4 µL of water, 0.5 µL of each primer, and 5 µL of DNA template with the final volume of 20 µL in a Bio-Rad CFX96 Real-Time PCR detection system.

## Molecular identification of the Duffy genotypes

Molecular identification of the Duffy genotype was based on a novel sequence-specific PCR primers (Amplification Refractory Mutation System, ARMS) that allows amplifications of the GATA-1 transcription factor-binding site of the *DARC* gene promoter only when the target allele (TCT for Duffy-positive and CCT for Duffy-negative genotypes) is detected in the sample. Different sets of primers were designed and tested using Duffy-positive and Duffy-negative DNA controls which were validated by PCR and Sanger sequencing (see S1 Appendix) [28] to detect the Duffy genotype by qPCR. The forward primers contain a nucleotide mismatch at the 3' end (T for Duffy-positive and C for Duffy-negative genotypes) in addition to a mismatch base near to the SNP at the 3' end (C for Duffy-positive and A for Duffy-negative genotypes), to avoid false-amplification and prevent binding to non-complementary sequence and terminate amplification. The reverse primer are the same (Table 2). All Duffy negative samples found by PCR ARMS have been confirmed by Sanger sequencing.

<table><caption>Table 1. Species-specific primers used for the detection and the species identification of *P. vivax* and *P. falciparum*.</caption><thead><tr><th>Species</th><th>Forward or Reverse</th><th>Primer</th></tr></thead><tbody><tr><td rowspan="2">P. vivax</td><td>Forward</td><td>5'-TGCTACAGGTGCATCTCTTGTATTC -3'</td></tr><tr><td>Reverse</td><td>5'-ATTTGTCCCCAAGGTAAAACG -3'</td></tr><tr><td rowspan="2">P. falciparum</td><td>Forward</td><td>5'-ATGGATATCTGGATTGATTTTATTTATGA -3'</td></tr><tr><td>Reverse</td><td>5'-TCCTCCACATATCCAAATTACTGC -3'</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0011326.t001
<table><caption>Table 2. AMRS PCR primers used for the detection of Duffy-positive and -negative genotypes.</caption><thead><tr><th>Duffy genotype</th><th>Forward or Reverse</th><th>Primer</th></tr></thead><tbody><tr><td rowspan="2">positive</td><td>Forward</td><td>5'- CCCTCATTAGTCCTTGGCTCTCAT -3'</td></tr><tr><td>Reverse</td><td>5'- CACCCTGTGCAGACAGTTCC -3'</td></tr><tr><td rowspan="2">negative</td><td>Forward</td><td>5'- CCCTCATTAGTCCTTGGCTCTAAC -3'</td></tr><tr><td>Reverse</td><td>5'- CACCCTGTGCAGACAGTTCC -3'</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0011326.t002

Briefly, polymorphisms in the GATA-1 were detected in two separate real-time PCR reactions. Real-time PCR reactions were performed with 10 µL of ssAdvanced mix sybgreen 2x, 10µM of primers and 5 µL of DNA extract, and amplified on the CFX96 PCR instrument (Bio-Rad, Marnes-la-Coquette, France). PCR amplification was carried out under the following conditions: heating at 98°C for 3 min, followed by 35 cycles of heating at 95°C for 20 s, and of annealing/extension at 62°C for 30 s and a final extension step at 72°C for 10 min. After completion of the amplification reaction, a melting curve (65-90°C) was generated to detect non-specific PCR amplification. For each run, homozygous (TCT/TCT), heterozygous (TCT/CCT) Duffy-positive, and homozygous Duffy-negative (CCT/CCT) samples were included as positive controls and a no template (water) as negative control (see S2 Appendix).

## Statistical analyses

Data were analysed using Statistical Package for Social Sciences (SPSS, version 25). Chi-squared, or Fisher's exact tests were used for frequency data (expressed as percentages). Associations between Duffy genotypes and different variables were assessed and two-sided p-values less than 0.05 were considered as statistically significant.

## Results

### Characteristics of study participants

A total of 361 PCR-confirmed *P. vivax* (mono or mixed) infected samples were included for the molecular analysis of Duffy status. Among them, 126 (34.9%) were female, the mean age was 25 years old (ranging from 1-80 years old), 108 (29.9%) had no occupation, and 152 (42.1%) participants live in rural areas (Table 3).

### Distribution of the Duffy genotypes by study sites, study participants, ethnicity, gender, and age

From the 361 samples tested, 345 (95.6%) of the *P. vivax* patients were Duffy-positives (272 (78.8%) were heterozygous and 73 (21.2%) were homozygous) and 16 (4.4%) were Duffy-negatives.

Based on Duffy genotyping, 73 (20.2%) were TCT/TCT, 272 (75.4%) were TCT/CCT and 16 (4.4%) were CCT/CCT. The Duffy-negative patients (CCT/CCT genotype) from Arbaminch, Gambella, Metehara and Shewa robit was 2.7%, 8.1%, 4.0%, and 7.1%, respectively (Fig 2). No Duffy-negatives were detected from Dubti. The highest proportion of Duffy-negatives was in Gambella (8.1%, 6/74) and the lowest in Arbaminch (2.7%, 2/74). There was no significant association between Duffy genotype and study site (p = 0.12). The highest proportion of Duffy-negatives was observed in the Nuer ethnic group (12.5%, 2/16) followed by the Agnuwak group (7.7%, 1/13). No Duffy-negatives was identified among the Somale, Afar, Tigre, Welayita and Hadiya ethnic groups. No significant association was found between Duffy genotype and ethnic group (p = 0.50). About 5.5% (13/235) male and 2.4% (3/126)
<table><caption>Table 3. Sociodemographic characteristics of study participants in Ethiopia (n = 361).</caption><thead><tr><th>Variables</th><th>Category</th><th>Number</th><th>Percent (%)</th></tr></thead><tbody><tr><td rowspan="2">Gender</td><td>Female</td><td>126</td><td>34.9</td></tr><tr><td>Male</td><td>235</td><td>65.1</td></tr><tr><td rowspan="4">Age stratification in years</td><td>&le;14</td><td>40</td><td>11.1</td></tr><tr><td>15-24</td><td>135</td><td>37.4</td></tr><tr><td>25-34</td><td>122</td><td>33.8</td></tr><tr><td>&ge;35</td><td>64</td><td>17.7</td></tr><tr><td rowspan="10">Ethnic</td><td>Amhara</td><td>106</td><td>29.4</td></tr><tr><td>Anuak</td><td>13</td><td>3.6</td></tr><tr><td>Somale</td><td>6</td><td>1.7</td></tr><tr><td>Oromo</td><td>76</td><td>21.1</td></tr><tr><td>Afar</td><td>52</td><td>14.4</td></tr><tr><td>Tigre</td><td>11</td><td>3.0</td></tr><tr><td>Welayita</td><td>30</td><td>8.3</td></tr><tr><td>Gamo</td><td>37</td><td>10.2</td></tr><tr><td>Sidama</td><td>11</td><td>3.0</td></tr><tr><td>Hadiya</td><td>3</td><td>0.8</td></tr><tr><td>Nuer</td><td>16</td><td>4.4</td></tr><tr><td rowspan="5">Study Sites</td><td>Arbaminch</td><td>74</td><td>20.5</td></tr><tr><td>Dubti</td><td>68</td><td>18.8</td></tr><tr><td>Gambella</td><td>74</td><td>20.5</td></tr><tr><td>Metehara</td><td>75</td><td>20.8</td></tr><tr><td>Shewarobit</td><td>70</td><td>19.4</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0011326.t003

female participants were Duffy-negatives, with no significant association between Duffy genotype and gender (p = 0.16). From different age stratification [22], the lowest proportion of Duffy-negatives was identified from the age groups ≤14 (2.5%, 1/40) and 25-34 (2.5%, 3/122), while the highest proportion was identified from the age group ≥35 (7.8%, 5/64), no significant association between Duffy genotype and age groups (p = 0.337). (Table 4). Similarly, among the Duffy-positive participants, the relative proportions of heterozygous and homozygous were not associated with either study site (p = 0.47), ethnic group (p = 0.47) or gender (p = 0.17) (Table 5).

## Distribution of the Duffy genotype by P. vivax infections

For the 16 *P. vivax* infections in Duffy negatives, 12 (75%) were mono-infection and four (25%) were mixed *P. vivax/P. falciparum* infections. For the 345 *P. vivax* infections in Duffy positives, 312 (90%) were mono-infections and 33 (~10%) were mixed *P. vivax/P. falciparum* infections. The proportion of mixed *P. vivax/P. falciparum* infections was two-fold higher in Duffy-negative than in Duffy positive individuals (p = 0.047), despite small samples for the Duffy negatives (Table 6). Of the 73 homozygous Duffy-positives, 68 (93.2%) were pure *P. vivax* infections and 5 (6.8%) were mixed *P. vivax/P. falciparum* infections, slightly lower than that observed in heterozygous Duffy-positives of which 244 (78.2%) were pure *P. vivax* infections and 28 (10.3%) were mixed *P. vivax/P. falciparum* infections (p = 0.37) (Table 6).

## Asexual and sexual parasitemia comparison by Duffy groups

Of the samples with an asexual stage density of <1,000 p/µL and ranging from 1,000 to 4,999 p/µL, 9/60 (15.0%) and 7/131 (5.3%) were obtained from Duffy-negative individuals,
<figure><img src="image_6.png" /><figcaption>Fig 2. Map showing the distribution of the Duffy genotypes at five different malaria endemic sites in Ethiopia.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0011326.g002

respectively. Interestingly, no Duffy-negatives were detected with asexual parasite density ≥5000 p/µL but only Duffy positive patients (n = 170; Table 7). The mean asexual parasite density in homozygous and heterozygous Duffy-positives was 12,165 p/µl (IQR25-75: 1,640–24,234 p/µl) and 11,655 p/µl (IQR25-75: 1,676–14,065 p/µl), respectively, significantly higher than that in Duffy-negatives (1,227p/µl; IQR25-75: 539–1,732p/µl) (*p*<0.001; Fig 3).

Of the samples with a sexual gametocytemia of <1,000 p/µL, 8/180 (4.4%) were detected in Duffy-negative individuals (Table 7). The other eight *P. vivax* infections in Duffy negatives were detected with no gametocytes. By contrast, 172/180 (95.6%) of the Duffy-positives were detected with sexual gametocytemia of <1,000 p/µL, 31 were detected with gametocytemia of 1,000–4999 p/µL, and 4 with gametocytemia of ≥5,000 p/µL.

Discussion

We report here the association between 361 pure *P. vivax* and mixed *P. vivax*/*P. falciparum* infections collected from five eco-epidemiological study sites in Ethiopia and the Duffy genotype of the patients. We assessed malaria infection by qPCR and polymorphisms in the
<table><caption>Table 4. Analysis of Duffy status by study sites (altitude), study participants ethnicity and gender (n = 361).</caption><thead><tr><th rowspan="2">Variables</th><th rowspan="2"></th><th colspan="3">Duffy Status</th><th rowspan="2">P-value</th></tr><tr><th>Positive n (%)</th><th>Negative n (%)</th><th>Total</th></tr></thead><tbody><tr><td rowspan="5">Study sites (altitude)</td><td>Arbaminch(1200 m)</td><td>72 (97.3)</td><td>2 (2.7)</td><td>74</td><td rowspan="5">0.089</td></tr><tr><td>Dubti (379 m)</td><td>68 (100)</td><td>0</td><td>68</td></tr><tr><td>Gambella (447 m)</td><td>68 (91.9)</td><td>6 (8.1)</td><td>74</td></tr><tr><td>Metehara (959 m)</td><td>72 (96.0)</td><td>3 (4.0)</td><td>75</td></tr><tr><td>Shewarobit(1268 m)</td><td>65 (92.9)</td><td>5 (7.1)</td><td>70</td></tr><tr><td rowspan="12">Ethnic of Study participant</td><td>Amhara</td><td>100 (94.3)</td><td>6 (5.7)</td><td>106</td><td rowspan="12">0.41</td></tr><tr><td>Agnuwak</td><td>12 (92.3)</td><td>1 (7.7)</td><td>13</td></tr><tr><td>Somale</td><td>6 (100)</td><td>0</td><td>6</td></tr><tr><td>Oromo</td><td>71 (93.4)</td><td>5 (6.6)</td><td>76</td></tr><tr><td>Afar</td><td>52 (100)</td><td>0</td><td>52</td></tr><tr><td>Tigre</td><td>11 (100)</td><td>0</td><td>11</td></tr><tr><td>welayita</td><td>30 (100)</td><td>0</td><td>30</td></tr><tr><td>Gamo</td><td>35 (94.6)</td><td>2 (5.4)</td><td>37</td></tr><tr><td>Sidama</td><td>11 (100)</td><td>0 (0)</td><td>11</td></tr><tr><td>Hadiya</td><td>3 (100%)</td><td>0</td><td>3</td></tr><tr><td>Nuer</td><td>14 (87.5)</td><td>2 (12.5)</td><td>16</td></tr><tr><td rowspan="2">Gender</td><td>Female</td><td>123 (97.6)</td><td>3 (2.4)</td><td>126</td><td rowspan="2">0.19</td></tr><tr><td>Male</td><td>222 (94.5)</td><td>13 (5.5)</td><td>235</td></tr><tr><td rowspan="4">Age in years</td><td>≤14</td><td>39 (97.5%)</td><td>1 (2.5%)</td><td>40</td><td rowspan="4">0.34</td></tr><tr><td>15–24</td><td>128 (94.8%)</td><td>7 (5.2%)</td><td>135</td></tr><tr><td>25–34</td><td>119 (97.5%)</td><td>3 (2.5%)</td><td>122</td></tr><tr><td>≥35</td><td>59 (92.2%)</td><td>5 (7.8%)</td><td>64</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0011326.t004

GATA-1 promoter affecting Duffy status by a novel ARMS-PCR approach. Here, we observed that the proportion of Duffy-negative (4.4%) among *P. vivax* infected patients was lower than previous studies conducted in different countries; Madagascar (9.3%) [28], different parts of Ethiopia (10.5–11.9%) [16], and Sudan 17.9% [9], but relatively higher than study reports such as Ethiopia 2.9% [22], south-eastern Iran 1.3% [29], Brazilian Amazon 0.6% [18]. Duffy blood groups remain highly polymorphic between different populations [30]. Historically, individuals with the Duffy-negative phenotype were thought to be resistant to *P. vivax* infection [18,19,31–33], but recent reports have shown *P. vivax* infection in Duffy-negative individuals [20,34–38] in many parts of Africa.

This study identified 4.4% (16/361) *P. vivax* infections in Duffy-negative patients. This is different from Central or West Africa, where the Duffy-negative genotype is much more common and almost all *P. vivax* infections occur among Duffy-negative individuals [37,39,40], despite relatively low number of *P. vivax* cases compared to East Africa. Several hypotheses have been put forward to explain the change in the pattern of infection in Duffy-negative individuals. One is a change in the host susceptibility, which in turn may drastically affect the distribution of *P. vivax* [9] and another is that *P. vivax* may have evolved to adapt to a wider range of hosts, environments, and highly competent transmission vectors in Africa [34]. Alternatively, *P. vivax* parasites may have developed alternative mechanisms to bind and invade reticulocytes as other host cell receptors have been shown to be involved in erythrocyte invasion [41].

Polymorphisms in Duffy blood group among different ethnic groups may contribute to the changing epidemiology of *vivax* malaria in Africa [42]. The Ethiopian population has a very
<table><caption>Table 5. Analysis of heterozygous and homozygous Duffy positive patients with study site (altitude), ethnicity, and gender (n = 312).</caption><thead><tr><th rowspan="2">Variables</th><th rowspan="2"></th><th colspan="3">Duffy Positive</th><th rowspan="2">P-value</th></tr><tr><th>Homozygous n (%)</th><th>Heterozygous n (%)</th><th>Total</th></tr></thead><tbody><tr><td rowspan="5">Study Sites (altitude)</td><td>Arbaminch</td><td>20 (27.8)</td><td>52 (72.2)</td><td>72</td><td rowspan="5">0.47</td></tr><tr><td>Dubti</td><td>12 (17.6)</td><td>56 (82.4)</td><td>68</td></tr><tr><td>Gambella</td><td>16 (23.5)</td><td>52 (76.5)</td><td>68</td></tr><tr><td>Metehara</td><td>12 (16.7)</td><td>60 (83.3)</td><td>72</td></tr><tr><td>Shewarobit</td><td>13 (20.0)</td><td>52 (80.0)</td><td>65</td></tr><tr><td rowspan="10">Ethnicity of study participants</td><td>Amhara</td><td>20 (20.0)</td><td>80 (80.0)</td><td>100</td><td rowspan="10">0.47</td></tr><tr><td>Agnuwak</td><td>4 (33.3)</td><td>8 (66.7)</td><td>12</td></tr><tr><td>Somale</td><td>0 (0.0)</td><td>6 (100)</td><td>6</td></tr><tr><td>Oromo</td><td>13 (18.3)</td><td>58 (81.7)</td><td>71</td></tr><tr><td>Afar</td><td>10 (19.2)</td><td>42 (80.8)</td><td>52</td></tr><tr><td>Tigre</td><td>2 (18.2)</td><td>9 (81.8)</td><td>11</td></tr><tr><td>welayita</td><td>7 (23.3)</td><td>23 (76.7)</td><td>30</td></tr><tr><td>Gamo</td><td>10 (28.6)</td><td>25 (71.4)</td><td>35</td></tr><tr><td>Sidama</td><td>5 (45.5)</td><td>6 (54.5)</td><td>11</td></tr><tr><td>Hadiya</td><td>0 (0)</td><td>3 (100)</td><td>3</td></tr><tr><td>Nuer</td><td>2 (14.3)</td><td>12 (85.7)</td><td>14</td></tr><tr><td rowspan="2">Gender</td><td>Female</td><td>31 (25.2)</td><td>92 (74.8)</td><td>123</td><td rowspan="2">0.17</td></tr><tr><td>Male</td><td>42 (18.9)</td><td>180 (81.1)</td><td>222</td></tr><tr><td rowspan="4">Age in years</td><td>≤14</td><td>4 (10.3%)</td><td>35 (89.7%)</td><td>39</td><td rowspan="4">0.284</td></tr><tr><td>15–24</td><td>31 (24.2%)</td><td>97 (75.8%)</td><td>128</td></tr><tr><td>25–34</td><td>24 (20.2%)</td><td>95 (79.8%)</td><td>119</td></tr><tr><td>≥35</td><td>14 (23.7%)</td><td>45 (76.3%)</td><td>59</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0011326.t005

heterogeneous ethnic composition. Eleven different ethnic groups were included in this study and the proportion of Duffy-negative and Duffy-positive varied between ethnic groups and study sites. A high proportion of Duffy-negative samples were observed among the Nuer followed by the Agnuwak ethnic groups, but no Duffy-negative was identified among the Somale, Afar, Tigre, Welayita and Hadiya ethnic groups. Other studies have also showed variation of Duffy genotype between ethnic groups. For example, in Sudan the Duffy-negative phenotype was identified in 83% of Shagia and 77% of Manasir tribes [43]; in Mauritania 54% of Moors were Duffy-positives while only 2% of black ethnic groups were Duffy-positives [42]; in Colombia Duffy genotypes were significantly associated with ethnicity [44]. These findings suggest that ethnic groups may be genetically distinct, resulting in different levels of susceptibility to *P. vivax* infection, in addition to different hypothesis such as *P. vivax* trains with a

<table><caption>Table 6. Analysis of *P. vivax* pure and mixed *P. vivax*/*P. falciparum* infections with *Duffy* genotype (n = 361).</caption><thead><tr><th rowspan="2">Variables</th><th rowspan="2"></th><th colspan="3">P. vivax infections</th><th rowspan="2">Chi-Square</th><th rowspan="2">P-value</th></tr><tr><th>Pure n (%)</th><th>Mixed n (%)</th><th>Total</th></tr></thead><tbody><tr><td rowspan="3">Duffy genotype</td><td>positive</td><td>312 (90.4)</td><td>33 (9.6)</td><td>345</td><td rowspan="3">3.96</td><td rowspan="3">0.047</td></tr><tr><td>negative</td><td>12 (75)</td><td>4 (25)</td><td>16</td></tr><tr><td>Total</td><td>324 (89.8)</td><td>37 (10.2)</td><td>361</td></tr><tr><td rowspan="3">DuffyPositive</td><td>homozygous</td><td>68 (93.2)</td><td>5 (6.8)</td><td>73</td><td rowspan="3">0.79</td><td rowspan="3">0.37</td></tr><tr><td>heterozygous</td><td>244 (89.7)</td><td>28 (10.3)</td><td>272</td></tr><tr><td>Total</td><td>312 (90.4)</td><td>33 (9.6)</td><td>345</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0011326.t006
<table><caption>Table 7. Association of Duffy status with asexual parasite density (n = 361).</caption><thead><tr><th rowspan="2">Variables</th><th rowspan="2"></th><th colspan="3">Duffy Status</th><th rowspan="2">Chi-Square</th><th rowspan="2">p-value</th></tr><tr><th>Positive n (%)</th><th>Negative n (%)</th><th>Total</th></tr></thead><tbody><tr><td rowspan="3">Asexual stage density</td><td>&lt;1000</td><td>51 (85.0)</td><td>9 (15.0)</td><td>60</td><td rowspan="3">23.96</td><td rowspan="3">&lt;0.001</td></tr><tr><td>1000-4999</td><td>124 (94.7)</td><td>7 (5.3)</td><td>131</td></tr><tr><td>&ge;5000</td><td>170 (100)</td><td>0 (0)</td><td>170</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0011326.t007

new capacity for erythrocyte invasion [28], or use several receptor-ligand interactions to tightly bind erythrocytes in the absence of a Duffy receptor [21].

Duffy-negatives with pure *P. vivax* and mixed *P. vivax/P. falciparum* infections was identified in the current study. The proportion of mixed *P. vivax/P. falciparum* infections in Duffy-negative was 25% in the present study, which was similar with the findings from Nigeria (25%).

<figure><img src="image_7.png" /><figcaption>Fig 3. Box plot showing the distribution of asexual parasite density among identified Duffy status. Key: Fy+/+; homozygous Duffy-positive, Fy+/-; Heterozygous Duffy-positive, Fy-/-; Duffy-negative.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0011326.g003
[45], and Cameroon (25%) [46], while it was lower than the study conducted in Madagascar (26.7%) [47], and Democratic Republic of the Congo (64%) [39]. The results provide evidence that Duffy-negative individuals can be infected by pure *P. vivax* alone or together with *P. falciparum*.

The mean asexual parasite density among Duffy-negatives (CCT/CCT) was 1,227 p/µl, around tenfold lower compared to homozygous Duffy-positives (TCT/TCT; 12,165 p/µl) and heterozygous Duffy-positive (TCT/CCT; 11,655 p/µl). Several prior studies have shown that *P. vivax* infections observed in Duffy-negatives have a significantly lower parasitemia than in Duffy-positive infections and are asymptomatic [13,16]. The low parasite densities found in Duffy-negative malaria patients may indicate that *P. vivax* uses an alternative pathway to invade Duffy-negative reticulocytes, which can be much less efficient than the classical pathway involving the interaction between PvDBP and DARC [22]. Alternatively, it has been hypothesized that the reticulocytes of Duffy-negative individuals may express the DARC protein in a transient and weak manner, thereby mediating a low invasion capacity [16]. Among the 16 Duffy negative samples detected in this study, half were detected with low-level gametocytes. The presence of gametocytes among Duffy negatives may further imply that these infections have the potential contributing to onward transmission.

To conclude, this study reports heterogeneous proportions of Duffy-negative individuals among *P. vivax* malaria patients across Ethiopia, confirms that Duffy-negativity does not provide complete protection against *P. vivax* infection, and these infections are often associated with low parasitemia, which may represent hidden reservoirs that can contribute to transmission. Therefore, the development of *P. vivax*-specific elimination strategies, including alternative antimalarial vaccines other than a Duffy-binding protein-based vaccine, should be facilitated by a better understanding of the relationship between Duffy-negativity and the invasion mechanisms of *P. vivax*.

**Supporting information**

S1 Appendix. Figures summarizing Sanger sequencing data.
(TIF)

S2 Appendix. Figures summarizing AMRS PCR-based allelic discrimination amplification curves in a control samples validated by PCR/Sanger sequencing.
(TIF)

**Acknowledgments**

We would like to thank i) SIDA project, ii) Human, Heredity and Health in Africa (H3Africa) [H3A-18-002]. H3Africa is managed by the Science for Africa Foundation (SFA Foundation) in partnership with Wellcome, NIH and AfSHG and iii) the Institute Pasteur, Paris, the University of Strasbourg, and the French Government (Agence Nationale de la Recherche; ANR-18-CE15-0026, ANR-21-CE35-0006 and the Laboratoire d'Excellence 'French Parasitology Alliance for Health Care', ANR-11-15 LABX-0024-PARAFRAP), and iv) EPHI for the materials support. We would also like to thank all the individuals who contributed technically to this study: the laboratory staff at the study sites, the staff at ALIPB (parasitology department), the staff at the Institute of Parasitology and Tropical Diseases, University of Strasbourg (Cécile Doderer-Lang), and those who were provided support on Arc Map and statistical analysis (Dereje Dilu, Abriham Keraleme and Chalie). We must express our sincere gratitude to all the regional health offices and the study sites for their cooperation. Last but not least, our sincere thanks go to all the participants who took part in this study.
**Author Contributions**

**Conceptualization:** Abnet Abebe.

**Data curation:** Abnet Abebe, Sisay Dugassa, Ashenafi Assefa, Jonathan J. Juliano, Eugenia Lo, Didier Menard, Lemu Golassa.

**Formal analysis:** Abnet Abebe, Isabelle Bouyssou, Solenne Mabilotte, Sisay Dugassa, Ashenafi Assefa, Jonathan J. Juliano, Eugenia Lo, Didier Menard, Lemu Golassa.

**Funding acquisition:** Abnet Abebe.

**Investigation:** Abnet Abebe, Didier Menard, Lemu Golassa.

**Methodology:** Abnet Abebe, Eugenia Lo, Didier Menard, Lemu Golassa.

**Project administration:** Abnet Abebe.

**Resources:** Abnet Abebe, Didier Menard, Lemu Golassa.

**Software:** Abnet Abebe.

**Supervision:** Abnet Abebe, Sisay Dugassa, Ashenafi Assefa, Jonathan J. Juliano, Eugenia Lo, Didier Menard, Lemu Golassa.

**Validation:** Abnet Abebe, Eugenia Lo, Didier Menard, Lemu Golassa.

**Visualization:** Abnet Abebe, Lemu Golassa.

**Writing – original draft:** Abnet Abebe.

**Writing – review & editing:** Abnet Abebe, Sisay Dugassa, Ashenafi Assefa, Jonathan J. Juliano, Eugenia Lo, Didier Menard, Lemu Golassa.

**References**

1. World Health Organization (WHO). World Malaria Report 2022. https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2022.

2. Poespoprodjo JR, Fobia W, Kenangalem E, Lampah DA, Hasanuddin A, Warikar N, et al. Vivax malaria: A major cause of morbidity in early infancy. Clinical Infectious Diseases. 2009; 48(12):1704–12. https://doi.org/10.1086/599041 PMID: 19438395

3. Kevin Baird J. Evidence and implications of mortality associated with acute plasmodium vivax malaria. Clinical Microbiology Reviews. 2013; 26(1):36–57. https://doi.org/10.1128/CMR.00074-12 PMID: 23297258

4. Moreira CM, Abo-Shehada M, Price RN, Drakeley CJ. A systematic review of sub-microscopic Plasmodium vivax infection. Malaria Journal. 2015; 14(1):1–10. https://doi.org/10.1186/s12936-015-0884-z PMID: 26390924

5. Ford A, Kepple D, Williams J, Kolesar G, Ford CT, Abebe A, et al. Gene Polymorphisms Among Plasmodium vivax Geographical Isolates and the Potential as New Biomarkers for Gametocyte Detection. Front. Cell. Infect. Microbiol. 2022; 11:789417.

6. Caudhurit Asok, Polyakovat Julia, Zbrzenat Valerea, William Kennett, Subhash Gulati Pogo OA. Cloning of glycoprotein D cDNA, which encodes the major subunit of the Duffy blood group system and the receptor for the Plasmodium vivax malaria parasite. Proc Natl Acad Sci USA. 1993; 90:10793–10797.

7. Tournamille C, Colin Y, Cartron JP, Le Van Kim C. "Disruption of a GATA motif in the Duffy gene promoter abolishes erythroid gene expression in Duffy-negative individuals." Nat Genet. 1995; 10(2):224–8. https://doi.org/10.1038/ng0695-224 PMID: 7663520

8. Howes RE, Patil AP, Piel FB, Nyangiri OA, Kabaria CW, Gething PW, et al. The global distribution of the Duffy blood group. Nat Commun. 2011; 2(4):266. https://doi.org/10.1038/ncomms1265 PMID: 21468018

9. Musab MA, Albsheer KP, Safaa A, Mohammed E, Eiman G, Salma MA E al. Distribution of Duffy Phenotypes among Plasmodium vivax Infections in Sudan. Genes. 2019; 10:417.
10. Popovici J, Roesch C, Rougeron V. The enigmatic mechanisms by which Plasmodium vivax infects Duffy-negative individuals. PLoS Pathog. 2020; 16(2):e1008258. https://doi.org/10.1371/journal.ppat.1008258 PMID: 32078643

11. Kepple D, Pestana K, Tomida J, Abebe A, Golassa L, Lo E. Alternative invasion mechanisms and host immune response to Plasmodium vivax malaria: Trends and future directions. Microorganisms. 2021; 9(1):1-17.

12. Daniels G. The molecular genetics of blood group polymorphism. Transplant Immunology. 2005; 14(3-4 SPEC. ISS.):143-53. https://doi.org/10.1016/j.trim.2005.03.003 PMID: 15982556

13. Gunalan K, Niangaly A, Thera MA, Doumbo OK, Miller LH. Plasmodium vivax Infections of Duffy-Negative Erythrocytes: Historically Undetected or a Recent Adaptation? Trends in Parasitology. 2018; 34(5):420-9. Available from: http://dx.doi.org/10.1016/j.pt.2018.02.006

14. Yazdanbakhsh K, Rios M, Storry JR, Kosower N, Parasol N, Chaudhuri A, et al. Molecular mechanisms that lead to reduced expression of duffy antigens. Transfusion. 2000; 40(3):310-20. https://doi.org/10.1046/j.1537-2995.2000.40030310.x PMID: 10738032

15. Schmid P, Ravenell KR, Sheldon SL FW. DARC alleles and Duffy phenotypes in African Americans. Transfusion. 2012; 52:1260-7. https://doi.org/10.1111/j.1537-2995.2011.03431.x PMID: 22082243

16. Loa E, Gianluca Russob, Pestanaa K, Kepplea D, Abageroc BR, Ghyslaine Bruna, et A. Contrasting epidemiology and genetic variation of Plasmodium vivax infecting Duffy-negative individuals across Africa. International Journal of Infectious Diseases. 2021; 108:63-71.

17. Lo E, Yewhalaw D, Zhong D, Zemene E, Degefa T, Tushune K. Molecular epidemiology of Plasmodium vivax and Plasmodium falciparum malaria among Duffy-positive and Duffy-negative populations in Ethiopia. Malaria Journal. 2015; 14(84). https://doi.org/10.1186/s12936-015-0596-4 PMID: 25884875

18. Carlos EC, Luiz CM, Álvaro AR, D'Almeida C, Vanja SC, Yuri G, et A. Duffy blood group gene polymorphisms among malaria vivax patients in four areas of the Brazilian Amazon region. Malaria Journal. 2007; 6:167. https://doi.org/10.1186/1475-2875-6-167 PMID: 18093292

19. Roesch C, Popovici J, Bin S, Run V, Kim S, Ramboarina S, et al. Genetic diversity in two Plasmodium vivax protein ligands for reticulocyte invasion. PLoS Negl Trop Dis. 2018;12(10):e0006555.

20. Mendes C, Dias F, Figueiredo J, Mora VG, Cano J, de Sousa B, et A. Duffy Negative Antigen Is No Longer a Barrier to Plasmodium vivax-Molecular Evidences from the African West Coast (Angola and Equatorial Guinea). PLoS Negl Trop Dis. 2011; 5(6):e1192.

21. Wilairatana P, Masangkay FR, Kotepui KU, Milanez GDJ, Kotepui M. Prevalence and risk of Plasmodium vivax infection among Duffy-negative individuals: a systematic review and meta-analysis. Scientific Reports. 2022;1-13. Available from: https://doi.org/10.1038/s41598-022-07711-5

22. Abate A, Bouyssou I, Mabilotte S, Lang CD, Dembele L, Menard D, et al. Vivax malaria in Duffy-negative patients shows invariably low asexual parasitaemia: implication towards malaria control in Ethiopia. Malaria Journal. 2022;1-10. Available from: https://doi.org/10.1186/s12936-022-04250-2

23. Tamirat Gebru Woldearegai Peter G. Kremsner, Jürgen F J. Kun BM. Plasmodium vivax malaria in Duffy-negative individuals from Ethiopia. Transactions of The Royal Society of Tropical Medicine and Hygiene. 2013; 107(5):328-31.

24. WHO. WHO Prequalification of Diagnostics Programme PUBLIC REPORT Product: SD Bioline Malaria Ag Pf/Pv Number: PQDx 0031-012-00. 2022. Available from: https://extranet.who.int/pqweb/sites/default/files/PQDx_0125-012-00_MalariaAgP.f-P.v_8.0.pdf

25. Debash H, Bisetegn H, Ebrahim H, Feleke DG, Gedefie A, Tilahun M, et al. Prevalence and associated risk factors of malaria among febrile under-five children visiting health facilities in Ziqala district, North-east Ethiopia: A multicenter cross-sectional study. PLoS ONE. 2022; 17(10 October):1-13. Available from: http://dx.doi.org/10.1371/journal.pone.0276899 PMID: 36301956

26. Nicholas D, Julie S, Aung P, Hadjar S, Armedy H, Enny K, et A. Gametocyte Dynamics and the Role of Drugs in Reducing the Transmission Potential of Plasmodium vivax. JID. 2013; 208:802-12.

27. Canier L, Khim N, Kim S, Sluydts V, Heng S, Dourng D, et al. An innovative tool for moving malaria PCR detection of parasite reservoir into the field. Malaria Journal. 2013; 12(1):1-12. https://doi.org/10.1186/1475-2875-12-405 PMID: 24206649

28. Ménard D, Barnadas C, Bouchier C, Henry-Halldin C, Gray LR, Ratsimbasoa A, et A. Plasmodium vivax clinical malaria is commonly observed in Duffy-negative Malagasy people. Proc Natl Acad Sci USA. 2010; 107(13):5967-71. https://doi.org/10.1073/pnas.0912496107 PMID: 20231434

29. Miri-moghaddam E, Bameri Z, Mohamadi M. Duffy blood group genotypes among malaria Plasmodium vivax patients of Baoulch population in southeastern Iran. Asian Pacific Journal of Tropical Medicine. 2014; 7(3):206-7. Available from: https://doi.org/10.1016/S1995-7645(14)60021-3 PMID: 24507640

30. Łukasik E, Nowak I, Czerwinski M, Waśniowska K. Duffy blood group system - the frequency of Duffy antigen polymorphisms and novel mutations in the Polish population. Transfus Apher Sci. 2019; 58 (2):156-161.
31. Guerra CA, Howes RE, Patil AP, Gething PW, van Boeckel TP, Temperley WH, et al. The international limits and population at risk of Plasmodium vivax transmission in 2009. PLoS Neglected Tropical Diseases. 2010; 4(8). https://doi.org/10.1371/journal.pntd.0000774 PMID: 20689816

32. Miller LH, Mason SJ, Clyde DF MM. The resistance factor to Plasmodium vivax in blacks: the Duffy-blood-group genotype, FyFy. N Engl J Med. 1976; 295:302-4.

33. Rondônia D, Brasileira AO, Cavasini CE, José F, Pereira T, Ribeiro WL, et al. Duffy blood group genotypes among malaria patients in Rondônia, Western Brazilian Amazon Genótipos do sistema sanguíneo Duffy em pacientes maláricos. 2001; 34(6):591-5.

34. Abdelraheem MH, Albsheer MMA, Mohamed HS, Amin M, Hamid MMA. Transmission of Plasmodium vivax in Duffy-negative individuals in central Sudan. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2016; 110(4):258-60. https://doi.org/10.1093/trstmh/trw014 PMID: 27076512

35. Zimmerman PA. Plasmodium vivax infection in Duffy-Negative People in Africa. Am J Trop Med Hyg. 2017;97. https://doi.org/10.4269/ajtmh.17-0461 PMID: 28990906

36. Dongho GBD, Gunalan K, L'Episcopia M, Paganotti GM, Menegon M, Sangong RE, et al. Plasmodium vivax Infections Detected in a Large Number of Febrile Duffy-Negative Africans in Dschang, Cameroon. American Journal of Tropical Medicine and Hygiene. 2021; 104(3):987-92. https://doi.org/10.4269/ajtmh.20-1255 PMID: 33534776

37. Oboh MA, Badiane AS, Ntadom G, Ndiaye YD, Diongue K, Diallo MA, et al. Molecular identification of Plasmodium species responsible for malaria reveals Plasmodium vivax isolates in Duffy negative individuals from southwestern Nigeria. Malaria Journal. 2018; 17(1):1-12. Available from: https://doi.org/10.1186/s12936-018-2588-7

38. Ryan JR, Stoute JA, Amon J, Dunton RF, Mtalib R, Koros J, et al. Evidence for transmission of Plasmodium vivax among a Duffy antigen negative population in Western Kenya. American Journal of Tropical Medicine and Hygiene. 2006; 75(4):575-81. PMID: 17038676

39. Brazeau NF, Whitesell AN, Doctor SM, Keeler C, Mwandagalirwa MK, Tshefu AK, et al. Plasmodium vivax infections in duffy-negative individuals in the democratic republic of the Congo. American Journal of Tropical Medicine and Hygiene. 2018; 99(5):1128-33. https://doi.org/10.4269/ajtmh.18-0277 PMID: 30203741

40. Poirier P, Doderer-Lang C, Atchade PS, Lemoine JP, De L'Isle MLC, Abou-Bacar A, et al. The hide and seek of Plasmodium vivax in West Africa: Report from a large-scale study in Beninese asymptomatic subjects. Malaria Journal. 2016; 15(1):1-9. https://doi.org/10.1186/s12936-016-1620-z PMID: 27887647

41. Gruszczyk J, Kanjee U, Chan LJ, Menant S, Malleret B, Lim NTY, et al. Transferrin receptor 1 is a reticulocyte-specific receptor for Plasmodium vivax. Science. 2018; 359(6371):48-55. https://doi.org/10.1126/science.aan1078 PMID: 29302006

42. Wurtz N, Mint Lekweiry K, Bogreau H, Pradines B, Rogier C, Ould Mohamed Salem Boukhary A, et al. Vivax malaria in Mauritania includes infection of a Duffy-negative individual. Malaria Journal. 2011; 10:1-8.

43. Kempi A, Knap O, Drozd A, Kaczmarczyk M, Parafiniuk M. Analysis for genotyping Duffy blood group in inhabitants of Sudan, the Fourth Cataract of the Nile. 2012;1-6.

44. Gonzalez L, Veg J, Ramirez Jose L, Bedoya G C-F. Relationship between Duffy blood groups genotypes and malaria infection in different ethnic groups of Choco- Colombia. Colomb Med [online] 2012, vol43, n3, pp189-195. 2012; 43(3):189-95.

45. Oboh MA, Singh US, Ndiaye D, Badiane AS, Ali NA, Bharti PK, et al. Presence of additional Plasmodium vivax malaria in Duffy negative individuals from Southwestern Nigeria. Malaria Journal. 2020;1-8. Available from: https://doi.org/10.1186/s12936-020-03301-w

46. Gaelle H, Mbenda N, Das A. Molecular Evidence of Plasmodium vivax Mono and Mixed Malaria Parasite Infections in Duffy-Negative Native Cameroonians. 2014; 9(8):1-9.

47. Ménard D, Barnadas C, Bouchier C, Henry-Halldin C, Gray LR, Ratsimbasoa A, et al. Plasmodium vivax clinical malaria is commonly observed in Duffy-negative Malagasy people. Proceedings of the National Academy of Sciences of the United States of America. 2010; 107(13):5967-71. https://doi.org/10.1073/pnas.0912496107 PMID: 20231434