# Genetic diversity of Plasmodium falciparum isolates from Baka Pygmies and their Bantu neighbours in the north of Gabon

Bertrand Mvé-Ondo, Dieudonné Nkoghe, Céline Arnathau, Virginie Rougeron, Ulrich Bisvigou, Lauriane Yacka Mouele, Larson Boundenga, Durand Patrick, Eric Elguero, Simone Lemmers, et al.

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

► To cite this version:

Bertrand Mvé-Ondo, Dieudonné Nkoghe, Céline Arnathau, Virginie Rougeron, Ulrich Bisvigou, et al. Genetic diversity of Plasmodium falciparum isolates from Baka Pygmies and their Bantu neighbours in the north of Gabon. Malaria Journal, 2015, 14 (1), 10.1186/s12936-015-0862-5. http://hal-03002433

HAL Id: hal-03002433

https://hal.umontpellier.fr/hal-03002433

Submitted on 25 May 2021

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_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# Genetic diversity of *Plasmodium falciparum* isolates from Baka Pygmies and their Bantu neighbours in the north of Gabon

Bertrand Mvé-Ondo¹,<fnref n="2" />*, Dieudonné Nkoghe¹,⁴, Céline Arnathau⁵, Virginie Rougeron¹,⁵, Ulrich Bisvigou⁴, Lauriane Yacka Mouele¹,³, Larson Boundenga¹,⁶, Patrick Durand⁵, Eric Elguero⁵, Simone Lemmers⁷, Lucrèce M. Délicat-Loembet¹, Nancy Diamella-Moukodoum¹, Christophe Paupy⁵, François Renaud⁵, Franck Prugnolle¹,⁵† and Benjamin Ollomo¹,³†

**Abstract**

**Background:** There have been many reports on the population genetic structure of *Plasmodium falciparum* from different endemic regions especially sub-Saharan Africa. However, few studies have been performed on neglected populations, such as the Pygmy populations. In this study, the population genetic structure of *P. falciparum* was investigated in the Baka Pygmies of Gabon and compared to that observed in neighboring villages composed mostly of Bantu farmers.

**Methods:** A total of 342 blood samples were collected from 170 Baka Pygmies and 172 Bantus in the north of Gabon (Woleu Ntem Province). *Plasmodium* infections were characterized by sequencing a portion of the parasite *cytochrome b* gene. Population genetic structure of *P. falciparum* in the different villages was analysed using microsatellite markers and genes coding for antigenic proteins (MSP1, MSP2, GLURP, and EBA-175).

**Results:** Overall, prevalence of *P. falciparum* was around 57 % and no significant difference of prevalence was observed between Pygmies and Bantus. No significant differences of population genetic structure of *P. falciparum* was found between Pygmy and Bantu people except for one antigen-coding gene, *glurp*, for which genetic data suggested the existence of a potentially disruptive selection acting on this gene in the two types of populations. The genetic structure of *P. falciparum* followed a pattern of isolation by distance at the scale of the study.

**Conclusion:** The prevalence and genetic diversity of *P. falciparum* observed in Baka demonstrates a significant transmission of the parasite in this population, and some exchanges of parasites with Bantu neighbours. Despite that, some antigen-coding genes seem to have had a particular evolutionary trajectory in certain Pygmy populations due to specific local human and/or mosquito characteristics.

**Keywords:** Malaria, Genetic diversity, Population structure, Microsatellite loci, Antigenic determinants, Baka Pygmies, Bantus

**Background**

Malaria is the most important parasitic disease based
on its impacts on human populations, with 207 million

cases and 627,000 deaths recorded in 2012, especially in
sub-Saharan Africa [1]. Currently, five species are known
to infect humans (*Plasmodium falciparum*, *Plasmodium
malariae*, *Plasmodium ovale*, *Plasmodium vivax*, and
*Plasmodium knowlesi*). Among them, *P. falciparum* is the
most virulent.

During the last 20 years, the genetic structure of *P. fal-
ciparum* populations has been intensely investigated all
over the world because it provides information on the

<footnote>*Correspondence: mvebertrand2000@yahoo.fr; bertrand.mveondo@cirm.fr</footnote>

<footnote>†Franck Prugnolle and Benjamin Ollomo contributed equally to the management of this study</footnote>

<footnote n="2">² Ecole Doctorale Régionale d'Afrique Centrale en Infectiologie Tropicale, BP 876, Franceville, Gabon</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
evolutionary potential of the parasite confronted by a
new drug or a vaccine, on the risk of diffusion of resist-
ance alleles throughout the world, as well as on its epi-
demiology and evolutionary history [2–11]. *Plasmodium*
*falciparum* populations on the African continent, where
most cases of infection occur every year, have largely
been studied in this context [12–16]. However, some
populations of this region have received very little atten-
tion regarding malaria infections. This is the case for
instance for the Pygmy populations.

Pygmies are indigenous hunter-gatherers of the Central African forest who diverged from other groups in Africa about 60,000 years ago [17, 18]. They traditionally lived in small, nomad groups whose livelihood strategies were based on hunting, gathering, small-scale farming, and exchange of forest products with farming neighbours. Pygmies were traditionally intimately connected to the forest. Nowadays, Pygmy communities continue to maintain a strong link with the forest whenever possible, but the large majority have now abandoned their traditional way of life, spending more time in roadside settlements, with closer contacts with neighbouring Bantu farming communities, and more reliance on farming and wage labour [19, 20]. Because Pygmy populations remained isolated in forest ecosystems for a long time and are genetically divergent from the Bantu populations (thus constituting a different environment for the parasite), it is likely that their populations of parasites also remained isolated and developed specific genetic characteristics, especially at genes directly interacting with the host immune system.

In Gabon (Central Africa), there are three distinct
Pygmy communities: the Baka, the Babongo and the
Bakoya. None of them has conserved a traditional way
of life even if some groups live in nomad forest camps
for several months for hunting at certain periods of the
year. In this country, most population genetic studies on
malaria were conducted on the Bantu populations (urban
and rural populations) [12, 21, 22]. Pygmies historically
received limited attention from the health authorities and
were largely neglected.

In this study, the population genetic structure of *P. falciparum* was characterized and compared among the Baka Pygmies and their Bantus neighbours in the same region of Gabon using both microsatellite and a set of antigen-coding genes.

**Methods**

Study area and studied populations

The study was carried out in Woleu Ntem Province,
located in the north of Gabon (Fig. 1). Three coun-
tries surround this province: Cameroon, Equatorial
Guinea and the Republic of Congo. The project focused

on Baka Pygmies and Bantus living in rural and forest
environments.

The study was conducted following two stages: (1) a sensitization stage in May 2013 during which persons were informed of the aim of the study, and, (2) blood sampling from June to July 2013. The sampling was conducted just after the rainy season (occurring from Febru- ary to May) when malaria transmission is the highest.

Participants were recruited from ten villages (Fig. 1;
Table 1): three in the Woleu region and seven in the region
of Ntem Minvoul. Two villages (Bitouga and Etho) corre-
sponded to what will be hereafter be referred as 'isolated
Pygmy villages' (IPV). These villages were located in iso-
lated forest zones, with no road access and were only
composed by Pygmies. Five villages were 'mixed villages'
(MV). Located along roads, these villages are mostly popu-
lated by Bantu peoples and a minority of Pygmy peoples.
In these villages, both Bantus and Pygmies were sampled
when possible (Table 1). Finally, three villages were villages
built by forestry companies for their workers [hereafter
named 'forestry villages' (FV)]. Only Bantu peoples lived
in these villages, accessible by road. Distances between vil-
lages ranged from 1 to 136 km (Fig. 1). In mixed villages
and other sites, participants were subjected to a question-
naire in which they would indicate their ethnic origin and
the origin of their parents. Participants born from inter-
ethnic marriages were not included in this study.

In each village, all volunteers were included in the study.
Volunteers were recruited with their informed consent
and in the case of children the consent of their parents.
The translation of consent for participants was done by
members of the NGO AGAFI (Gabonese Association of
Assistance to Indigenous Women and Indigent). Before
blood collection, a medical doctor examined each volun-
teer. Age and sex were retrieved for each patient. Blood
was collected in EDTA tubes using venipuncture. Each
sample was aliquoted and stored at -20 °C until used.

A total of 342 blood samples were obtained: 170 from Baka Pygmies (121 from MV and 49 from IPV) and 172 from Bantus (38 in MV and 134 in FV). This study was authorized by the Ministry of Health and the agreement of the National Ethics Committee on Research (authorization number: PROT No 0030/2013/SG/CNE).

Determination of *Plasmodium falciparum* infections

For each blood sample, total DNA was extracted from
200 µl of blood using the DNeasy blood and tissue kit
Qiagen (QIAamp® DNA Blood Mini Kit) according to
manufacturer's instructions.

For each sample, *Plasmodium* infections were detected
using the method based on the amplification of a por-
tion of the *Plasmodium* cytochrome b (*cyt-b*) gene by
a nested PCR as described in Additional file 1 [23]. The
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Geographical location of the study sites in the departments of Woleu and High Ntem in Woleu Ntem Province. Ten sites were sampled: (1) green isolated Pygmy villages; (2) red mixed villages (Bantus and Pygmies) and (3) blue Bantu populations in forestry villages</figcaption></figure>

PCR amplified products were sequenced by Eurofins (MWG, France). Then, a logistic regression (binomial family) was used to test for variations of *P. falciparum* prevalence between Pygmies and Bantus (fixed effect) or between village types (fixed effect). In both cases, the factor ‘village of collection’ was considered as a random effect. Tests were realized using the R software [24] and the library lme4 [25, 26].

## Molecular genotyping of *Plasmodium falciparum* isolates

Seven microsatellite markers as well as four candidate genes coding for antigenic proteins were genotyped.

The microsatellite markers amplified by semi-nested PCR are shown in Additional file 2. The PCR primer sets and amplification conditions were those described by [9, 27]. Fluorescence-labelled PCR products were sized on ABI Prism 310 genetic analyzer (Applied Biosystems), with a Genescan 500 LIZ internal size standard.

Regarding the candidate genes, merozoite surface protein 1 (*msp1*—block 2), merozoite surface protein 2 (*msp2*—block 3), glutamate-rich protein (*glurp*-region 2), and erythrocyte binding protein-175 (*eba-175*—region III) genes were genotyped as described in [21, 22, 28]. *Msp1* is a single copy gene located on chromosome 9. It encodes a merozoite surface antigen of approximately 190 kDa. The C-terminal part (*MSP1*<sub>19</sub>) is immunogenic and contains Epidermal Growth Factor (EGF)-like domains. *MSP1*<sub>19</sub> is very rich in cysteine and is involved in the invasion of red blood cells by merozoites [29, 30]. *MSP1* is composed of 17 blocks. Block 2, which is the first variable block, located at the N-terminal region is the most polymorphic. Block 2 of *msp1* was genotyped in this study. *MSP-2* is located on chromosome 2 and encodes a 46–53 kDa antigenic protein located at the surface of the merozoite [31]. It consists of three regions: two semi-conserved (at the N-and C-terminal regions)
<table><caption>**Table 1** Number of blood samples collected per village and ethnic group (Pygmies and Bantu) as well as prevalence (%) of *Plasmodium falciparum*</caption><thead><tr><th rowspan="2">Collection sites and GPS positions</th><th colspan="2">Pygmies</th><th colspan="2">Bantu</th></tr><tr><th>Number of samples</th><th>Number of P. falciparum infections (prevalence %)</th><th>Number of samples</th><th>Number of P. falciparum infections (prevalence %)</th></tr></thead><tbody><tr><td colspan="5"><strong>Isolated Pygmy villages</strong></td></tr><tr><td>Bitouga (2°07.133'N, 12°07.770'E)</td><td>21</td><td>13 (62)</td><td>0</td><td>–</td></tr><tr><td>Etho (2°05.728'N, 12°09.460'E)</td><td>28</td><td>9 (32)</td><td>0</td><td>–</td></tr><tr><td colspan="5"><strong>Mixed villages</strong></td></tr><tr><td>Doumassi (2°09.212'N, 12°07.781'E)</td><td>32</td><td>13 (40.6)</td><td>5</td><td>1 (20)</td></tr><tr><td>Zangaville (2°04.474'N, 12°02.476'E)</td><td>13</td><td>5 (38.4)</td><td>10</td><td>6 (60)</td></tr><tr><td>Esseng (2°09.726'N, 12°08.321'E)</td><td>54</td><td>30 (55.5)</td><td>0</td><td>–</td></tr><tr><td>NkokAkom (2°10.900'N, 12°07.826'E)</td><td>11</td><td>8 (72.7)</td><td>12</td><td>5 (41.6)</td></tr><tr><td>Mimbang (2°11.593'N, 12°07.620'E)</td><td>11</td><td>7 (63.6)</td><td>11</td><td>8 (72.7)</td></tr><tr><td colspan="5"><strong>Forestry villages</strong></td></tr><tr><td>TTIB (1°36. 201'N, 11°59. 986'E)</td><td>0</td><td>–</td><td>56</td><td>39 (69.6)</td></tr><tr><td>Bordamure (1°02.305'N, 11°41.629'E)</td><td>0</td><td>–</td><td>39</td><td>21 (53.8)</td></tr><tr><td>STIBG (1°00.472'N, 11°55.789'E)</td><td>0</td><td>–</td><td>39</td><td>31 (79.4)</td></tr><tr><td><strong>Total</strong></td><td>170</td><td>85 (50)</td><td>172</td><td>111 (64.5)</td></tr></tbody></table>

and one polymorphic region composed of tandemly repetitive sequences [32–34]. The central polymorphic region (block 3) of *msp2* was genotyped in this study. The EBA-175 antigen is a molecule expressed at the surface of the parasite, which binds to red blood cells. It is involved in the invasion of the red blood cells [35–37]. It is a dimorphic molecule composed of two alleles: CAMP (C) and FCR3 (F). Dimorphism is related to two segments of region III. This is this region that has been genotyped in this study. Finally, GLURP is a 220 kDa protein expressed in both pre-erythrocytic and erythrocytic stages of *P. falciparum*, as well as on the surface of newly released merozoites in human host [38]. It consists of a N-terminal region of limited diversity named R0 followed by two polymorphic repeat regions named R1 and R2 [39]. The region 2 (R2) of GLURP has been genotyped in this study.

The specific primers used for genotyping each gene are given in Additional file 1. For *msp1* gene, three allelic families (K1, MAD20 and RO33) were determined. The two allelic families 3D7 and FC27 of the *msp2* gene were considered and the F (FCR3) and C (CAMP) alleles for the gene *eba-175*. For *glurp*, alleles were defined based on the amplicon size. PCR products were all visualized on a 1.5 % agarose gel electrophoresis.

### Plasmodium falciparum population genetic analyses

Allele frequency was computed as the proportion of the total of all alleles detected among the isolates examined (see Additional file 3).

All population genetic analyses were performed using FSTAT version 2.9.4 [40]. Genetic diversity within populations was assessed by the unbiased expected heterozygosity (*He*) [41]. Comparison of genetic diversity was performed either between Pygmy and Bantu populations or between IPs against the others village types using a permutation test (1000 permutations).

Population genetic differentiation was estimated using the estimator ($\theta$) of the Wright's *F*-statistics $F_{ST}$ [42, 43]. This index measures the genetic differentiation (i.e., the differences of allelic frequencies) between populations of interest. Deviation of $\theta$ from 0 was tested using a permutation test (1000 permutations). Mantel tests were used to test for patterns of isolation by distance (IBD). A partial Mantel test was performed to test for an effect of the ethnic group on *P. falciparum* genetic differentiation while controlling for geographic distance (1000 permutations). For this test, the first matrix was the matrix of $F_{ST}$ computed between all pairs of populations. The second matrix was the matrix of pairwise geographic distances computed between populations and the third one was a matrix containing 1 when the populations belonged to the same ethnicity (Pygmy/Pygmy and Bantu/Bantu) and 2 when they belonged to different ethnic groups (Pygmy/Bantu).

An additional partial Mantel test was performed with two new groups of populations: one including only the IPs that are strictly inhabited by Pygmies and the other including the other villages that are mostly inhabited by
Bantus (MV + FV). MV isolates were indeed considered more likely to be representative of Bantus because, inside MVs, pygmies represent only a minority of individuals (even if they were sometimes sampled in majority) and also because no genetic differentiation was observed between *P. falciparum* populations collected from Pygmies and Bantus within these villages (see "Results"). For this test, the first matrix was the matrix of $F_{ST}$ between populations, the second matrix corresponded to the geographic distances between populations and the third one contained 1 when the populations belonged to the same village type (IPV or MV + FV) and 2 when they belonged to two village types (IPV/MV + FV).

## Results

### Malaria prevalence and Plasmodium species

A total of 342 blood samples were analysed: 172 from Bantus and 170 from Pygmies. Distribution and origin of blood samples per village, village types (IPV, MV and FV) and ethnic groups are given in Table 1. A total of 196 samples were detected positive to *P. falciparum*. Prevalence varied from one village to another and is given in Table 1. Logistic regressions revealed no effect of the ethnic group (Pygmy vs Bantu, *P-value* = 0.154) on the prevalence of infection (while accounting for variations between villages) but differences between village types. Thus, although no significant differences were observed between IPVs and MVs (*P-value* = 0.49), a significant difference was observed between FVs

and IPVs (*P-value* = 0.029) and between FVs and MVs (*P-value* = 0.041). Average infection prevalence was significantly higher in the villages from the forest concessions (FV) than in any other village types (IPV + MV).

### Microsatellite population structure

Seven microsatellite markers were genotyped in a total of 196 *P. falciparum* isolates from two IPVs, five Bantu/ Pygmy: MVs, and from three villages of FVs representing 111 isolates from Bantus and 85 from Pygmies. All markers were highly polymorphic with a total number of alleles per locus observed varying from sixto 17. For each village type or ethnic group (Pygmy/Bantu), mean expected heterozygosity (*He*) is given in Table 2. No significant difference was observed between the estimates of genetic diversity (*He*) compared either between ethnic groups (Pygmies vs Bantu, *P-value* = 0.24) or village types (IPV vs MV + FV, *P-value* = 0.25).

Regarding genetic differentiation, the average $F_{ST}$ ($\theta$) computed over all loci and populations was low ($\theta$ = 0.005) and not significantly different from 0 (*P-value* = 0.102). Pairwise $F_{ST}$ ($\theta$) are given in Table 3. Within MVs with isolates from Pygmies and Bantus, no significant genetic differentiation was observed between the *P. falciparum* isolates collected from the different ethnic groups (Table 3). At a regional scale, a significant isolation by distance was observed among populations (Fig. 2, correlation coefficient r = 0.387; *P-value* = 0.002).

<table><caption>**Table 2** Mean genetic diversity (He) of Plasmodium falciparum populations estimated from the different village types and the different ethnic groups (All Pygmies and All Bantu) for microsatellite markers and antigen-coding genes</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Isolated pygmyvillages (IPV)</th><th colspan="2">Mixed villages (MV)</th><th rowspan="2">Forestry villages (FV)</th><th rowspan="2">All pygmies</th><th rowspan="2">All bantu</th></tr><tr><th>Pygmies</th><th>Bantu</th></tr></thead><tbody><tr><td colspan="7"><strong>Microsatellites</strong></td></tr><tr><td>Ta60</td><td>0.835</td><td>0.817</td><td>0.828</td><td>0.789</td><td>0.826</td><td>0.808</td></tr><tr><td>Ara2</td><td>0.833</td><td>0.877</td><td>0.889</td><td>0.853</td><td>0.855</td><td>0.871</td></tr><tr><td>Taa87</td><td>0.819</td><td>0.847</td><td>0.916</td><td>0.832</td><td>0.833</td><td>0.874</td></tr><tr><td>Pfpk2</td><td>0.930</td><td>0.808</td><td>0.889</td><td>0.890</td><td>0.869</td><td>0.890</td></tr><tr><td>Polyα</td><td>0.944</td><td>0.899</td><td>0.932</td><td>0.891</td><td>0.922</td><td>0.912</td></tr><tr><td>Ta1</td><td>0.944</td><td>0.929</td><td>0.941</td><td>0.848</td><td>0.937</td><td>0.894</td></tr><tr><td>Taa81</td><td>0.814</td><td>0.834</td><td>0.792</td><td>0.835</td><td>0.824</td><td>0.814</td></tr><tr><td>Mean</td><td>0.874</td><td>0.859</td><td>0.884</td><td>0.848</td><td>0.866</td><td>0.866</td></tr><tr><td colspan="7"><strong>Candidate genes</strong></td></tr><tr><td>MSP2</td><td>0.711</td><td>0.836</td><td>0.603</td><td>0.741</td><td>0.773</td><td>0.672</td></tr><tr><td>MSP1</td><td>0.821</td><td>0.718</td><td>0.829</td><td>0.724</td><td>0.769</td><td>0.776</td></tr><tr><td>EBA-175</td><td>0.220</td><td>0.361</td><td>0.423</td><td>0.384</td><td>0.291</td><td>0.404</td></tr><tr><td>GLURP</td><td>0.359</td><td>0.741</td><td>0.612</td><td>0.644</td><td>0.550</td><td>0.628</td></tr><tr><td>Mean</td><td>0.528</td><td>0.664</td><td>0.617</td><td>0.623</td><td>0.596</td><td>0.620</td></tr></tbody></table>

Values within each category correspond to the average of genetic diversity computed within each population
<table><caption>Table 3 Pairwise FST estimates (θ) measured between each pair of populations</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Group</th><th rowspan="2">Village</th><th colspan="2">Isolated Pygmy villages</th><th colspan="8">Mixed villages</th><th colspan="5">Forestry Villages</th></tr><tr><th>Pygmies</th><th>Etho</th><th>Pygmies</th><th>Doumassi</th><th>Esseng</th><th>Pygmies</th><th>Mimbang</th><th>Bantus</th><th>Pygmies</th><th>NkokAkom</th><th>Bantus</th><th>Pygmies</th><th>NkokAkom</th><th>Bantus</th><th>Bordaville</th><th>Bordamure</th><th>STIBG</th><th>Bantus</th><th>TTIB</th></tr></thead><tbody><tr><td rowspan="10">Microsatellites</td><td rowspan="2">Isolated Pygmy villages</td><td>Pygmies</td><td>Bitouga</td><td rowspan="2">-0.0256</td><td rowspan="2">-0.0101</td><td>-0.0033</td><td>-0.0068</td><td>-0.0016</td><td>-0.0005</td><td>-0.039</td><td>0.0075</td><td>-0.0217</td><td>0.0075</td><td>0.0025</td><td>-0.0017</td></tr><tr><td>Pygmies</td><td>Etho</td><td>-0.0138</td><td>-0.0168</td><td>0.0162</td><td>-0.0135</td><td>-0.0175</td><td>-0.033</td><td>-0.0117</td><td>-0.0212</td><td>0.008</td><td>0.0045</td><td>-0.0073</td></tr><tr><td rowspan="2">Mixed villages</td><td>Pygmies</td><td>Doumassi</td><td rowspan="2">-0.0069</td><td rowspan="2">-0.0175</td><td>0.0116</td><td>-0.0035</td><td>-0.0188</td><td>0.0392</td><td>-0.0301</td><td>0.0277</td><td>0.0058</td><td>0.02</td><td>0.0113</td><td>0.0141</td></tr><tr><td>Pygmies</td><td>Esseng</td><td>-0.0022</td><td>0.0146</td><td>0.0164</td><td>-0.0166</td><td>-0.0202</td><td>-0.0138</td><td>0.0151</td><td>0.0228</td><td>0.0142</td><td>0.0319</td></tr><tr><td></td><td>Pygmies</td><td>Mimbang</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Bantus</td><td>Mimbang</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Pygmies</td><td>NkokAkom</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Bantus</td><td>NkokAkom</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Pygmies</td><td>Zangaville</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Bantus</td><td>Zangaville</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Forestry villages</td><td>Bantus</td><td>Bordamure</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Bantus</td><td>STIBG</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="10">Candidate genes</td><td rowspan="2">Isolated Pygmy villages</td><td>Pygmies</td><td>Bitouga</td><td rowspan="2">-0.0417</td><td rowspan="2">0.0421</td><td>0.1063</td><td>0.1041</td><td>0.1398</td><td>0.0372</td><td>0.2664*</td><td>0.149</td><td>0.1301</td><td>0.1632</td><td>0.1509*</td><td>0.1428</td></tr><tr><td>Pygmies</td><td>Etho</td><td>0.0795</td><td>0.0244</td><td>0.1629</td><td>0.0272</td><td>0.2818</td><td>0.0324</td><td>0.0543</td><td>0.1335</td><td>0.1627*</td><td>0.1245</td></tr><tr><td rowspan="2">Mixed villages</td><td>Pygmies</td><td>Doumassi</td><td rowspan="2">0.0469</td><td rowspan="2">-0.0436</td><td>0.0255</td><td>-0.064</td><td>0.1047</td><td>-0.0423</td><td>-0.0001</td><td>0.0128</td><td>0.114*</td><td>0.0194</td><td>0.01047*</td></tr><tr><td>Pygmies</td><td>Esseng</td><td>0.1123</td><td>0.1308</td><td>0.1079</td><td>0.2229*</td><td>-0.0417</td><td>0.0472</td><td>0.065</td><td>0.0752</td><td>0.1047*</td></tr><tr><td></td><td>Pygmies</td><td>Mimbang</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Bantus</td><td>Mimbang</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Pygmies</td><td>NkokAkom</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Bantus</td><td>NkokAkom</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Pygmies</td><td>Zangaville</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Bantus</td><td>Zangaville</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">Forestry villages</td><td>Bantus</td><td>Bordamure</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Bantus</td><td>STIBG</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Values in italics are FST values obtained within mixed villages between P. falciparum from Pygmies and Bantus

<footnote>* P-value <0.05 (after Bonferroni correction)</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Relationship between pairwise geographic distances (in km) and pairwise *F*<sub>ST</sub> estimates (θ) measured between each pair of populations, for microsatellite markers</figcaption></figure>

<table><caption>**Table 4** Partial correlation coefficients of partial Mantel tests realized for microsatellites (all loci) and each of the four candidate genes</caption><thead><tr><th></th><th>Distance</th><th>Within/between ethnies</th><th>Within/between village types</th></tr></thead><tbody><tr><td>Microsatellites</td><td>0.387**</td><td>-0.145 ns</td><td>-0.15 ns</td></tr><tr><td>MSP2</td><td>0.258 ns</td><td>0.0156 ns</td><td>0.15 ns</td></tr><tr><td>MSP1</td><td>0.08 ns</td><td>-0.06 ns</td><td>-0.23 ns</td></tr><tr><td>EBA-175</td><td>0.025 ns</td><td>-0.05 ns</td><td>0.02 ns</td></tr><tr><td>GLURP</td><td>0.025 ns</td><td>0.09 ns</td><td>0.54**</td></tr><tr><td>All candidate genes</td><td>0.197 ns</td><td>0.11 ns</td><td>0.39**</td></tr></tbody></table>

ns non-significant

** *P*-value ≤0.01. *P*-values are given after Bonferroni corrections for multiple testing

Partial Mantel tests revealed no significant effect of the ethnic group or of the village type (Table 4).

## Candidate gene population structure

The genotyping of the different genes (*msp1*, *msp2*, *eba-175*, and *glurp*) was also performed using 196 *P. falciparum* isolates. For these genes, alleles were classified according to the size of the amplified fragment. The different alleles and their distribution in the different collection sites are presented in Additional file 3. Unbiased genetic diversity (*He*) estimated for each gene, each type of village and each population group is presented in Table 2. Comparison among ethnic groups revealed no difference between Pygmies and Bantus for any of these genes. The only significant difference was observed for *glurp* between IPV and the other villages (MV + FV, *P*-value = 0.018).

Regarding genetic differentiation, the average *F*<sub>ST</sub> measured between populations over all loci was low (*F*<sub>ST</sub> = 0.09) and significantly different from 0 (*P*-value ≤0.001).

Within MVs, no significant genetic differentiation was observed between the isolates collected in Pygmies and Bantus (Table 3; Additional file 4). Although non-significant, the Mantel test revealed a similar regional pattern as the one observed for microsatellites, namely a global isolation by distance (r = 0.197; *P*-value = 0.093) (Fig. 3). Partial Mantel tests were realized for each gene. Results are given in Table 4. As shown, correlation coefficients for the ethnic group effect were non-significant. Regarding the village type effect (IPVs vs the others), all correlation coefficients were not significantly different from 0 except for GLURP for which a significant strong negative correlation coefficient was observed meaning that, for this gene, genetic differentiation was higher between village types (IPVs vs the others) than within.

## Discussion

In this study, the characterization of the genetic polymorphism of *P. falciparum* was under taken in isolated Baka
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Relationship between pairwise geographic distances (in km) and pairwise FST estimates (θ) measured between each pair of populations, for antigenic markers</figcaption></figure>

Pygmy villages in Gabon and compared to neighbouring villages composed of mixed Bantu/Pygmy populations or villages of logging companies (only Bantus).

**Plasmodium falciparum infection prevalence**

Over all samples (n = 342), around 57 % of the individuals were infected by *P. falciparum* (n = 196). No significant difference was observed in the prevalence of *P. falciparum* between Pygmies and Bantus. This result is at odds with various studies realized in the middle of the 20th Century in Central Africa, showing a global malaria infection rate lower in Pygmy camps than in surrounding Bantu villages. Explanations for these differences varied from one study to another. Schwetz attributed the differences to ethnic characteristics because they considered that the different human communities were submitted to a similar malaria risk [44]. However, some authors have noted a far lower number of malaria vectors in Pygmy camps (even in those temporarily sedentarized) than in neighbouring Bantu villages as well as differences in the composition of the vector communities [45, 46]. Thus, *Anopheles gambiae s.s.*, the main vector of *Plasmodium* in Equatorial regions, was mostly absent from the undergrowth and so the Pygmy camps because of the ecological preferences of their larvae for temporarily sunned breeding sites [47].

The absence of differences observed in this study between the IPVs and the neighbouring Bantu villages (especially the MVs) can mainly be explained by the drastic cultural changes Pygmies have experienced in the last 50 years. Most Pygmies are no longer nomadic and now live in settlements very similar to those of the Bantus. Pygmy communities spend more time outside the forest in fixed settlements, within-village population densities have increased as well as their contact with the Bantus

[20]. As a consequence, their settlements are now favourable to the development of *An. gambiae*, the main vector of *P. falciparum* and to an efficient malaria transmission. All these factors may explain the increase of malaria prevalence in present Pygmy communities and the similarity in the intensity of transmission with the neighbouring Bantu communities.

**Plasmodium falciparum population genetic structure: microsatellites**

Because microsatellite markers are expected to be neutral, population genetic structure estimated using these markers is considered to be mainly the consequence of demographic events and not selection acting on these populations. Regarding *P. falciparum* population genetic diversity and its distribution, no substantial difference between village types or ethnic groups (Pygmies vs. Bantus) was observed. Genetic diversity as measured by expected heterozygosity was high in the different population groups and equivalent to levels of diversity measured in other African areas with a hyper endemic transmission [27, 48] (as is the case for Gabon).

Levels of genetic differentiation measured between villages were small and not significantly higher between populations belonging to different village types (IPV, MV and FV) or ethnic groups. No significant genetic differentiation was observed between Pygmies and Bantus within MVs. Over all the studied area, a pattern of isolation by distance was observed (genetic differentiation is on average lower for populations that are geographically close) thus suggesting that exchanges of parasites among villages occur mostly at small geographic distances and that no particular barrier to migration exists between Pygmy and Bantu communities (Fig. 2).
## Plasmodium falciparum population genetic structure: candidate genes

The distribution of genetic diversity at candidate genes, if under selection, is the consequence of the demographic history of the populations (as for microsatellites) but also of the action of selection. Comparing their population structure to that of microsatellites is thus a way to detect phenomenon of selection and in particular phenomenon of local adaptation for certain genes.

Overall, similar patterns of population structure were observed at the candidate genes, namely no differences in the levels of genetic diversity among village types or ethnic groups and a global isolation by distance at the scale of the study (although only marginally significant). The only substantial difference was observed for the gene GLURP. Not only the level of genetic diversity of this gene was significantly lower in the IPVs (Table 2) than in the other village types but also the level of genetic differentiation was significantly more important between village types (IPV against MV and FV) than inside, everything else being equal (geographical distance). This suggests therefore the existence of a potential phenomenon of disruptive selection acting on this gene between the IPVs and the other villages.

GLURP is expressed at different stages in the parasite including mosquito and vertebrate stages [49]. If disruptive selection indeed occurs at this gene, this could then be due to differential effects of the Pygmy and Bantu immune systems or to the mosquito vectors that are involved in transmission in the different village types. If ethnicity is the selective factor, one could however wonder why no significant genetic differences were observed within mixed villages between isolates collected from Pygmies or Bantus. This could be due to a lack of statistical power (due to small sample sizes) or to the strong parasite migration occurring between ethnic groups within village, a phenomenon known to prevent the evolution of local adaptation [50].

The findings of this are compatible with results of another study comparing the level of population structure among *P. falciparum* populations collected in different regions of the world at genes encoding ten leading vaccine antigens [51], including the four candidate genes studied here. In this study, the strongest level of genetic differentiation among different populations was also observed at the *glurp* gene. Authors attributed this observation to several possible factors, including variations in human genetic resistance to malaria or adaptation to different anopheline species transmitting *P. falciparum* worldwide [51].

## Conclusion

Very few studies have been conducted on the population genetic structure of *P. falciparum* in Pygmies.

In this study, the prevalence and genetic diversity of *P. falciparum* observed in Baka demonstrates a significant transmission of the parasite in this population, and some exchanges of parasites with the Bantu neighbours. Despite that, some antigen-coding genes seem to have had a particular evolutionary trajectory in certain Pygmy populations due to specific local human and/or mosquito characteristics. It is necessary to extend this research to other sites, including the Bakoya Pygmies in the northeast and Babongo Pygmies in southern Gabon, as well as to other sites in Central Africa. It would in particular be interesting to analyse the prevalence and population genetic structure in the last Pygmy groups that have conserved a traditional nomadic and forest-dependent way of life, such as in the Central African Republic.

## Additional files

**Additional file 1:** Sequences of primers used for genotyping *P. falciparum* genes. This table shows the primers sequences (cytochrome b, MSP1, MSP2, EBA 175 and GLURP) used.

**Additional file 2:** Microsatellite loci characteristics for *P. falciparum*. Presentation of the seven microsatellite sequences used and their repeat motif.

**Additional file 3:** Allelic diversity of MSP-1, MSP-2, EBA-175 and Glurp of *P. falciparum* in Bantus and Baka Pygmies. Distribution of alleles of genes MSP1, MSP2, EBA 175 and GLURP by Population Group.

**Additional file 4:** Pairwise *F*<sub>ST</sub> estimates (θ) measured between each pair of populations for each of the four candidate genes (MSP2, MSP1, EBA-175, GLURP). Values in bold are *F*<sub>ST</sub> values obtained within mixed villages between *P. falciparum* from Pygmies and Bantus. *P-value* < 0.05 (after Bonferroni correction). NA: not available. Pairwise *F*<sub>ST</sub> estimates (θ). Values in bold are *F*<sub>ST</sub> values obtained within mixed villages between *P. falciparum* from Pygmies and Bantus. *P-value* < 0.05 (after Bonferroni correction). NA: not available.

## Authors' contributions

BMO, DN, UB, and BO contributed to the acquisition of samples in fieldwork; BMO, FP, LYM, CA, EE, VR, BO, LB, NDM, LMDL, SL, PD, and FR analysed and interpreted the data; BMO, BO, FP, DN, PD, CP, VR, and FR conducted and supervised this work; BMO, BO, FP, and FR wrote this paper. All authors read and approved the final manuscript.

## Author details

## Acknowledgements

The Authors thank the two reviewers for their constructive comments. The study was funded by Centre International de Recherches Médicales de Franceville (CIRMF, Gabon), Agence Nationale de la Recherche (ANR, France,

<footnote n="1">¹ Centre International de Recherches Médicales de Franceville (CIRMF), BP 769, Franceville, Gabon. ² Ecole Doctorale Régionale d'Afrique Centrale en Infectiologie Tropicale, BP 876, Franceville, Gabon. ³ Fac de Médecine, UMR-BIOMED, BP 8507, Libreville, Gabon. ⁴ Ministry of Health, BP 50, Libreville, Gabon.</footnote>

<footnote n="5">⁵ Laboratoire Maladies Infectieuses et Vecteurs, Ecologie, Génétique, Evolution et Contrôle (MIVEGEC), UMR CNRS 5290/IRD 224, Université Montpellier, CHRU de Montpellier, 39 Avenue Charles Flahault, 34295 Montpellier, France.</footnote>

<footnote n="6">⁶ Département de Biologie Animale, Faculté des Sciences et Techniques, Laboratoire d'Écologie et Biologie évolutives, Université Cheikh AntaDiop de Dakar, BP 5005, Dakar, Senegal. ⁷ Department of Anthropology, Durham University, Dawson Building, South Road, Durham DH1 3LE, England, UK.</footnote>
grant ORIGIN JCJC 2012),Office des Ports et Rades du Gabon (OPRAG) and Institut de Recherche pour le développement (IRD, France). We thank NGO AGAFI for the sensitizing and interpreter role in participants, the nurses Fabrice Mounguegui and Sheila Minkoue for participation in the collection of blood and the technician Thélesfort Mbang-Mboro for sample processing. We also thank Aye Lydie for her help.

## Compliance with ethical guidelines

### Competing interests

The authors have declared that they have no competing interests.

Received: 6 June 2015 Accepted: 23 August 2015

Published online: 09 October 2015

## References

1. WHO. World Malaria Report 2013. Geneva: World Health Organization; 2013.

2. Day KP, Marsh K. Naturally acquired immunity to *Plasmodium falciparum*. *Immunol Today*. 1991;12:A68–71.

3. Mwingira F, Nkwengulila G, Schoepflin S, Sumari D, Beck HP, Snounou G, et al. *Plasmodium falciparum* msp1, msp2 and glurp allele frequency and diversity in sub-Saharan Africa. *Malar J*. 2011;10:79.

4. Healer J, Murphy V, Hodder AN, Masciantonio R, Gemmill AW, Anders RF, et al. Allelic polymorphisms in apical membrane antigen-1 are responsible for evasion of antibody-mediated inhibition in *Plasmodium falciparum*. *Mol Microbiol*. 2004;52:159–68.

5. Manske M, Miotto O, Campino S, Auburn S, Almagro-Garcia J, Maslen G, et al. Analysis of *Plasmodium falciparum* diversity in natural infections by deep sequencing. *Nature*. 2012;487:375–9.

6. Prugnolle F, Durand P, Jacob K, Razakrainibe F, Arnathau C, Villarreal D, et al. A comparison of *Anopheles gambiae* and *Plasmodium falciparum* genetic structure over space and time. *Microbes Infect*. 2008;10:269–75.

7. Larranaga N, Mejia RE, Hormaza JI, Montoya A, Soto A, Fontecha GA. Genetic structure of *Plasmodium falciparum* populations across the Honduras-Nicaragua border. *Malar J*. 2013;12:354.

8. Annan Z, Durand P, Ayala FJ, Arnathau C, Awono-Ambene P, Simard F, et al. Population genetic structure of *Plasmodium falciparum* in the two main African vectors, Anopheles gambiae and Anopheles funestus. *Proc Natl Acad Sci USA*. 2007;104:7987–92.

9. Anderson TJ, Su XZ, Bockarie M, Lagog M, Day KP. Twelve microsatellite markers for characterization of *Plasmodium falciparum* from finger-prick blood samples. *Parasitology*. 1999;119:113–25.

10. Yalcindag E, Elguero E, Arnathau C, Durand P, Akiana J, Anderson TJ, et al. Multiple independent introductions of *Plasmodium falciparum* in South America. *Proc Natl Acad Sci USA*. 2012;109:511–6.

11. Mobegi VA, Loua KM, Ahoudi AD, Satoguina J, Nwakanma DC, Amambua-Ngwa A, et al. Population genetic structure of *Plasmodium falciparum* across a region of diverse endemicity in West Africa. *Malar J*. 2012;11:223.

12. Aubouy A, Migot-Nabias F, Deloron P. Polymorphism in two merozoite surface proteins of *Plasmodium falciparum* isolates from Gabon. *Malar J*. 2003;2:12.

13. Aubouy A, Migot-Nabias F, Deloron P. Correlations between treatment outcome and both anti-MSP119 antibody response and erythrocyte-related genetic factors in *Plasmodium falciparum* malaria. *Infect Genet Evol*. 2007;7:147–54.

14. Peyerl-Hoffmann G, Jelinek T, Kilian A, Kabagambe G, Metzger WG, von Sonnenburg F. Genetic diversity of *Plasmodium falciparum* and its relationship to parasite density in an area with different malaria endemicities in West Uganda. *Trop Med Int Health*. 2001;6:607–13.

15. Soulama I, Nebie I, Ouedraogo A, Gansane A, Diarra A, Tiono AB, et al. *Plasmodium falciparum* genotypes diversity in symptomatic malaria of children living in an urban and a rural setting in Burkina Faso. *Malar J*. 2009;8:135.

16. Olasehinde GI, Yah CS, Singh R, Ojuronbge OO, Ajayi AA, Valecha N, et al. Genetic diversity of *Plasmodium falciparum* field isolates from south western Nigeria. *Afr Health Sci*. 2012;12:355–61.

17. Quintana-Murci L, Quach H, Harmant C, Luca F, Massonnet B, Patin E, et al. Maternal traces of deep common ancestry and asymmetric gene flow between Pygmy hunter-gatherers and Bantu-speaking farmers. *Proc Natl Acad Sci USA*. 2008;105:1596–601.

18. Verdu P, Austerlitz F, Estoup A, Vitalis R, Georges M, Thery S, et al. Origins and genetic diversity of pygmy hunter-gatherers from Western Central Africa. *Curr Biol*. 2009;19:312–8.

19. Ohenjo N, Willis R, Jackson D, Nettleton C, Good K, Mugarura B. Health of Indigenous people in Africa. *Lancet*. 2006;367:1937–46.

20. Stephens C, Porter J, Nettleton C, Willis R. Disappearing, displaced, and undervalued: a call to action for Indigenous health worldwide. *Lancet*. 2006;367:2019–28.

21. Ranford-Cartwright LC, Taylor J, Umasunthar T, Taylor LH, Babiker HA, Lell B, et al. Molecular analysis of recrudescent parasites in a *Plasmodium falciparum* drug efficacy trial in Gabon. *Trans R Soc Trop Med Hyg*. 1997;91:719–24.

22. Ntoumi F, Ngoundou-Landji J, Lekoulou F, Luty A, Deloron P, Ringwald P. Site-based study on polymorphism of *Plasmodium falciparum* MSP-1 and MSP-2 genes in isolates from two villages in Central Africa. *Parassitologia*. 2000;42:197–203.

23. Prugnolle F, Durand P, Neel C, Ollomo B, Ayala FJ, Arnathau C, et al. African great apes are natural hosts of multiple related malaria species, including *Plasmodium falciparum*. *Proc Natl Acad Sci USA*. 2010;107:1458–63.

24. R-Development-Core-Team. R: a language and environment for statistical computing. In *R Foundation for Statistical Computing*, Vienna, Austria http://www.R-projetorg/, 2014.

25. Doran H, Bates D, Bliese P, Dowling M. Estimating the multilevel Rasch Model: with the lme4 package. *J Stat Softw*. 2007;20:1–18.

26. Bates D, Maechler M, Bolker B. lme4: Linear Mixed-Effects Models Using S4 Classes. R package version 0.99375-38. http://CRAN.R-project.org/package=lme4, 2011.

27. Leclerc MC, Durand P, de Meeus T, Robert V, Renaud F. Genetic diversity and population structure of *Plasmodium falciparum* isolates from Dakar, Senegal, investigated from microsatellite and antigen determinant loci. *Microbes Infect*. 2002;4:685–92.

28. Toure FS, Mavoungou E, Ndong JM, Tshipamba P, Deloron P. Erythrocyte binding antigen (EBA-175) of *Plasmodium falciparum*: improved genotype determination by nested polymerase chain reaction. *Trop Med Int Health*. 2001;6:767–9.

29. Blackman MJ, Heidrich HG, Donachie S, McBride JS, Holder AA. A single fragment of a malaria merozoite surface protein remains on the parasite during red cell invasion and is the target of invasion-inhibiting antibodies. J Exp Med. 1990;172:379–82.

30. Holder AA, Guevara Patino JA, Uthaipibul C, Syed SE, Ling IT, Scott-Finnigan T, et al. Merozoite surface protein 1, immune evasion, and vaccines against asexual blood stage malaria. *Parassitologia*. 1999;41:409–14.

31. Clark JT, Donachie S, Anand R, Wilson CF, Heidrich HG, McBride JS. 46-53 kilodalton glycoprotein from the surface of *Plasmodium falciparum* merozoites. Mol Biochem Parasitol. 1989;32:15–24.

32. Smythe JA, Coppel RL, Day KP, Martin RK, Oduola AM, Kemp DJ, et al. Structural diversity in the *Plasmodium falciparum* merozoite surface antigen 2. *Proc Natl Acad Sci USA*. 1991;88:1751–5.

33. Fenton B, Clark JT, Khan CM, Robinson JV, Walliker D, Ridley R, et al. Structural and antigenic polymorphism of the 35- to 48-kilodalton merozoite surface antigen (MSA-2) of the malaria parasite *Plasmodium falciparum*. Mol Cell Biol. 1991;11:963–71.

34. Irion A, Beck HP, Felger I. New repeat unit and hot spot of recombination in FC27-type alleles of the gene coding for *Plasmodium falciparum* merozoite surface protein 2. Mol Biochem Parasitol. 1997;90:367–70.

35. Sim BK, Chitnis CE, Wasniowska K, Hadley TJ, Miller LH. Receptor and ligand domains for invasion of erythrocytes by *Plasmodium falciparum*. Science. 1994;264:1941–4.

36. Klotz FW, Orlandi PA, Reuter G, Cohen SJ, Haynes JD, Schauer R, et al. Binding of *Plasmodium falciparum* 175-kilodalton erythrocyte binding antigen and invasion of murine erythrocytes requires N-acetylneuraminic acid but not its O-acetylated form. Mol Biochem Parasitol. 1992;51:49–54.

37. Orlandi PA, Klotz FW, Haynes JD. A malaria invasion receptor, the 175-kilodalton erythrocyte binding antigen of *Plasmodium falciparum* recognizes the terminal Neu5Ac(alpha 2-3)Gal- sequences of glycophorin A. J Cell Biol. 1992;116:901–9.
38. Borre MB, Dziegiel M, Hogh B, Petersen E, Rieneck K, Riley E, et al. Primary structure and localization of a conserved immunogenic *Plasmodium falciparum* glutamate rich protein (GLURP) expressed in both the preerythrocytic and erythrocytic stages of the vertebrate life cycle. Mol Biochem Parasitol. 1991;49:119-31.

39. de Stricker K, Vuust J, Jepsen S, Oeuvray C, Theisen M. Conservation and heterogeneity of the glutamate-rich protein (GLURP) among field isolates and laboratory lines of *Plasmodium falciparum*. Mol Biochem Parasitol. 2000;111:123-30.

40. Goudet J. FSTAT version 2.9.4: a program to estimation genetics parameters. Updated from Goudet (1995). Available from http://www.unilch/izea/software/fstathtm 2003.

41. Nei M. Modification of linkage intensity by natural selection. Genetics. 1967;57:625-41.

42. Wright S. The genetical structure of populations. Ann Eugen. 1951;15:323-54.

43. Weir BS, Cockerham CC. Estimating F-statistic for the analysis of population structure. Evolution. 1984;38:1358-70.

44. Schwetz J, Baumann H, Peel A, Droeshaut F. Etude comparative de la malaria chez les Pygmées et les indigènes ordinaires de la forêt de l'Ituri (Congo Belge). Bull Soc Path Exot. 1933;26:639-51.

45. Languillon J, Mouchet J. Résultats des enquêtes épidémiologiques et entomologiques sur le paludisme dans les régions de Lom et Kadei et du

Boumba-Ngoko. Rapport ronéotypé IRCAM (ORSTOM) Yaoundé: Cameroon; 1955.

46. Carnevale P, Rickenbach A, Frezil JL. Les maladies transmises et leurs vecteurs chez les Pygmées Aka de Basse Lobaye. Paris: RCA. DGRST; 1978.

47. Mouchet J, Carnevale P, Coosemans M, Fontenille D, Ravaon-Jaahary C, Richard A, et al. Typologie du paludisme en Afrique. Cahiers de santé. 1993;3:220-30.

48. Bogreau H, Renaud F, Bouchiba H, Durand P, Assi SB, Henry MC, et al. Genetic diversity and structure of African *Plasmodium falciparum* populations in urban and rural areas. Am J Trop Med Hyg. 2006;74:953-9.

49. Lopez-Barragan MJ, Lemieux J, Quinones M, Williamson KC, Molina-Cruz A, Cui K, et al. Directional gene expression and antisense transcripts in sexual and asexual stages of *Plasmodium falciparum*. BMC Genom. 2011;12:587.

50. Lenormand T. Gene flow and the limits to natural selection. Trends Ecol Evol. 2002;17:183-9.

51. Barry AE, Schultz L, Buckee CO, Reeder JC. Contrasting population structures of the genes encoding ten leading vaccine-candidate antigens of the human malaria parasite, *Plasmodium falciparum*. PLoS One. 2009;4:e8497.

Submit your next manuscript to BioMed Central
and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

( ) **BioMed** Central