Skip to main content

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

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.

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. falciparum populations has been intensely investigated all over the world because it provides information on the evolutionary potential of the parasite confronted by a new drug or a vaccine, on the risk of diffusion of resistance alleles throughout the world, as well as on its epidemiology and evolutionary history [211]. Plasmodium falciparum populations on the African continent, where most cases of infection occur every year, have largely been studied in this context [1216]. However, some populations of this region have received very little attention 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 countries 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.

Fig. 1
figure 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

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 February 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) corresponded to what will be hereafter be referred as ‘isolated Pygmy villages’ (IPV). These villages were located in isolated 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 populated 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 villages ranged from 1 to 136 km (Fig. 1). In mixed villages and other sites, participants were subjected to a questionnaire 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.

Table 1 Number of blood samples collected per village and ethnic group (Pygmies and Bantus) as well as prevalence (%) of Plasmodium falciparum

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 volunteer. 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 portion of the Plasmodium cytochrome b (cyt-b) gene by a nested PCR as described in Additional file 1 [23]. The 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 lme 4 [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 (MSP119) is immunogenic and contains Epidermal Growth Factor (EGF)-like domains. MSP119 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) and one polymorphic region composed of tandemly repetitive sequences [3234]. 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 [3537]. 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 IPVs against the others village types using a permutation test (1000 permutations).

Population genetic differentiation was estimated using the estimator (θ) 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 θ 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 IPVs 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 Bantus, P-value = 0.24) or village types (IPV vs MV + FV, P-value = 0.25).

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 Bantus) for microsatellite markers and antigen-coding genes

Regarding genetic differentiation, the average F ST (θ) computed over all loci and populations was low (θ = 0.005) and not significantly different from 0 (P-value = 0.102). Pairwise F ST (θ) 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). Partial Mantel tests revealed no significant effect of the ethnic group or of the village type (Table 4).

Table 3 Pairwise F ST estimates (θ) measured between each pair of populations
Fig. 2
figure 2

Relationship between pairwise geographic distances (in km) and pairwise F ST estimates (θ) measured between each pair of populations, for microsatellite markers

Table 4 Partial correlation coefficients of partial Mantel tests realized for microsatellites (all loci) and each of the four candidate genes

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 ST measured between populations over all loci was low (F ST = 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.

Fig. 3
figure 3

Relationship between pairwise geographic distances (in km) and pairwise F ST estimates (θ) measured between each pair of populations, for antigenic markers

Discussion

In this study, the characterization of the genetic polymorphism of P. falciparum was under taken in isolated Baka 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.

References

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

    Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  PubMed Central  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Prugnolle F, Durand P, Jacob K, Razakandrainibe 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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  PubMed Central  PubMed  Google Scholar 

  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.

    Article  PubMed Central  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Mobegi VA, Loua KM, Ahouidi 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.

    Article  PubMed Central  PubMed  Google Scholar 

  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.

    Article  PubMed Central  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  PubMed Central  PubMed  Google Scholar 

  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.

    PubMed Central  CAS  PubMed  Google Scholar 

  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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  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.

    Article  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    CAS  PubMed  Google Scholar 

  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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  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.

    Article  Google Scholar 

  26. Bates D, Maechler M, Bolker B. lme4: Linear Mixed-Effects Models Using S4 Classes. R package version 0.999375-38. http://CRAN.R-projectorg/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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  30. Holder AA, Guevara Patino JA, Uthaipibull 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.

    CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

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

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

    PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  Google Scholar 

  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.

    Google Scholar 

  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é: Cameroun; 1955.

    Google Scholar 

  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.

    Google Scholar 

  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.

    Google Scholar 

  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.

    CAS  PubMed  Google Scholar 

  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.

    Article  CAS  Google Scholar 

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

    Article  Google Scholar 

  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.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

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.

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, 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.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Bertrand Mvé-Ondo.

Additional information

Franck Prugnolle and Benjamin Ollomo contributed equally to the management of this study

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 GLUPR by Population Group.

Additional file 4:

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

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mvé-Ondo, B., Nkoghe, D., Arnathau, C. et al. Genetic diversity of Plasmodium falciparum isolates from Baka Pygmies and their Bantu neighbours in the north of Gabon. Malar J 14, 395 (2015). https://doi.org/10.1186/s12936-015-0862-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-015-0862-5

Keywords