Skip to main content

Plasmodium falciparum gametocyte carriage in symptomatic patients shows significant association with genetically diverse infections, anaemia, and asexual stage density

Abstract

Background

Multi-genotype malaria infections are frequent in endemic area, and people commonly harbour several genetically distinct Plasmodium falciparum variants. The influence of genetic multiplicity and whether some specific genetic variants are more or less likely to invest into gametocyte production is not clearly understood. This study explored host and parasite-related risk factors for gametocyte carriage, and the extent to which some specific P. falciparum genetic variants are associated with gametocyte carriage.

Methods

Gametocytes and asexual forms were detected by light microscopy on thick smears collected between 2010 and 2012 in Nanoro, Burkina Faso. Merozoite surface protein 1 and 2 were genotyped by nested PCR on clinical samples. Associations between gametocyte carriage and factors, including multiplicity of infection, parasite density, patient age, gender, haemoglobin (Hb) level, and body temperature were assessed. The relationship between the presence of a particular msp1 and msp2 genetic variants and gametocyte carriage was also explored.

Results

Of the 724 samples positive to P. falciparum and successfully genotyped, gametocytes were found in 48 samples (6.63%). There was no effect of patient gender, age and body temperature on gametocyte carriage. However, the probability of gametocyte carriage significantly increased with increasing values of multiplicity of infection (MOI). Furthermore, there was a negative association between parasite density and gametocyte carriage. MOI decreased with parasite density in gametocyte-negative patients, but increased in gametocyte carriers. The probability of gametocyte carriage decreased with Hb level. Finally, the genetic composition of the infection influenced gametocyte carriage. In particular, the presence of RO33 increased the odds of developing gametocytes by 2 while the other allelic families K1, MAD20, FC27, and 3D7 had no significant impact on the occurrence of gametocytes in infected patients.

Conclusion

This study provides insight into potential factors influencing gametocyte production in symptomatic patients. The findings contribute to enhance understanding of risk factors associated with gametocyte carriage in humans.

Trial registration NCT01232530.

Background

Malaria is an infectious disease caused by Plasmodium parasites that are transmitted through the bites of infected female Anopheles mosquitoes. It is a life-threatening disease claiming about 405,000 lives each year worldwide [1]. Burkina Faso, like several other sub-Saharan African countries, ranks among the 10 highest burden countries with over 7,875,575 cases and 4,294 malaria-related deaths every year [1]. Children under 5 years old and pregnant women remain the most vulnerable groups.

Of Plasmodium species infecting humans, Plasmodium falciparum is the most dangerous involved in about 99% of overall malaria-related deaths [1]. This is the most common Plasmodium species in tropical settings, including in Burkina Faso. Multi-genotype malaria infections are frequent in endemic area and people commonly harbour several genetically distinct P. falciparum variants [2,3,4,5]. Multi-genotype infections may result from several distinct mono-infected mosquito bites overtime and/or from a single multiple-infected mosquito bite [6, 7]. In human infections, highly polymorphic markers such as Merozoite Surface Protein 1 and 2 (msp1 and msp2) are commonly used as polymorphic markers to distinguish parasite genetic variants [5, 8, 9]. Using this approach allows differentiation of distinct parasite genotypes within a particular infection ranging from 1 (monoclonal infection) up to 10 (multiclonal infection) in endemic settings [3, 9,10,11,12].

During blood stage in humans, malaria parasites undergo either asexual or sexual differentiation leading to gametocyte production. Gametocytes are the sexual forms of the parasite assuring the continuity of the parasite life cycle in Anopheles mosquitoes. Although crucial for the understanding of the infectivity cycle and disease transmission, the ultimate and proximate mechanisms leading to gametocyte differentiation are not fully elucidated [13,14,15,16].

Contrasting predictions can be made regarding the relationship between gametocyte production and multiplicity of infection (MOI). First, the existence of within-host competition among distinct genotypes of malaria parasites [5, 9, 17] could result in an increased investment in reproduction and transmission to flee the competition [18, 19]. Second, in presence of competitors, parasites could divert resource away from transmission towards replication to maximize their competitive success and hence their within-host survival [20, 21]. These contrasting strategies illustrate the two opposite ends of the well-known general survival versus reproduction trade-off [22]. In malaria parasites, support for both of these strategies were previously provided. Some studies demonstrated that clonal multiplicity can promote both the longevity of P. falciparum infection in patients and their ability to produce gametocytes [23,24,25]. Accordingly, positive relationship between MOI and gametocyte carriage has been reported [6, 26], and can translate into a positive correlation between MOI and mosquito infection rates [6]. Other studies found that investment in gametocytes is reduced in mixed infections [27, 28]. These results are in apparent contradiction and seem difficult to reconcile but gametocyte production is a complex process. A recent modelling approach demonstrated that parasites may display different patterns of transmission investment even in controlled conditions with animal models [29].

Besides genetic complexity, other factors generally associated with stressful conditions and pointing to deterioration in host conditions can drive gametocyte differentiation [13, 15]. These drivers include drug treatment, immune response, host anaemia, and nutritional status [15,16,17,18, 30,31,32,33,34]. Male gender was also previously reported as an independent risk factor associated with gametocyte carriage [35]. Another level of complexity is whether some specific genetic variants are more or less likely to invest into gametocyte production. In animal models, the existence of clones of Plasmodium berghei with different capability to produce gametocytes were previously described [36]. Genetically determined capability of some clones to produce gametocytes after a prolonged asexual stage was also reported in vitro [36, 37] and in vivo [38]. In natural human infections, a previous study suggested that some particular P. falciparum clones have a relatively higher capacity to produce gametocytes [39]. Causal relationship underlying sexual forms differentiation is difficult to establish in human infections. However, associated parasite and host-related risk factors could be identified through correlation between these factors and gametocyte carriage. Genes located in chromosome 9 of the parasite and variation in DNA repeats were both listed as factors involved in gametocytogenesis mechanism [40, 41]. Variation in DNA repeats represents the main molecular basis of the differentiation among msp1 (chromosome 9) and msp2 allelic families. Therefore, there is a good reason to postulate that allelic variability in msp1 and msp2 genes could influence gametocytogenesis in humans. This study explored parasite- and host-related factors associated with gametocyte carriage and the extent to which some specific P. falciparum variants are associated with gametocyte carriage in symptomatic patients in Burkina Faso.

Methods

Source of samples and microscopic analyses

Samples for this report were collected between 2010 and 2012 during a pharmacovigilance for artemisinin-based combination therapy study carried out in Nanoro, Burkina Faso. The latter was described in detail elsewhere [42, 43]. Malaria slides were prepared from peripheral blood obtained from finger prick and were stained with 3% Giemsa for 30 min. Slides were double read using Olympus CX21 microscope (Olympus Corporation, Tokyo, Japan) for the detection of asexual forms as first intention against 200 white blood cells and a negative result were declared after examination of 100 microscopic fields. During that process, when sexual forms were met, sexual form density was determined against 500 white blood cells count. Final parasite density represented arithmetic mean of the two readings. A third reader was appealed in case of discrepancies (huge difference between the two readers) defined as follow: (i) difference in Plasmodium species identification; (ii) positive versus negative results; and, (iii) if the higher count divided by the lower count was ≥ 2. In that case, the two closest readings among the three were considered.

Haemoglobin and body temperature

Haemoglobin (Hb) level was measured using Hemocue® 301+ (HemoCue AB, Ängelholm, Sweden). Undiluted blood obtained from finger prick was drawn into the microcuvette, which was inserted in the analyzer and Hb value in g/dL was read immediately. Electronic digital thermometers (Omron, Dalian, China) were used for the measurement of axillary temperature. The tip of the thermometer was inserted under the armpit and numeric value of axillary temperature in degree Celsius was recorded.

Molecular analyses

Parasite genomic DNA was extracted from dried blood spots collected during the screening process of patients, i.e., before treatment administration using QIamp DNA Kit (Qiagen, Hilden, Germany) following manufacturer’s instruction. DNA extract was used for merozoite surface proteins (msp) 1 block 2 gene (K1, MAD20 and RO33) and msp2 block 3 gene (3D7 and FC27) genotyping by nested PCR as previously described. Briefly, 5 µL of DNA extract was used to initiate the first PCR round using a Mastercycler® Gradient, (Eppendorf, Hamburg, Germany) and a Biometra thermal cycler (Analytik Jena, Jena, Germany) PCR machines. For the primary multiplex PCR round, the cycling conditions consisted of an initial denaturation step of 5 min at 94 °C, followed by 36 cycles of 1 min at 94 °C, 2 min at 58 °C and 2 min at 72 °C, and a final extension step of 10 min at 72 °C. One µl of this PCR product was used as DNA template to launch the nested round. For the nested msp1 round, the cycling conditions consisted of an initial denaturation step of 5 min at 94 °C, followed by 30 cycles of 1 min at 94 °C, 2 min at 59 °C and 2 min at 72 °C, and a final extension step of 10 min at 72 °C. For the nested msp2 round, the cycling conditions consisted of an initial denaturation step of 2 min at 94 °C, followed by 40 cycles of 30 s at 94 °C, 45 s at 58 °C and 1.5 min at 72 °C, and a final extension step of 10 min at 72 °C. PCR amplicons were detected under UV light transillumination after ethidium bromide-stained agarose gel electrophoresis. DNA fragments sizes were calculated using Photo CaptMW (version 11.01) regarding fragments sizes of the reference which was a 100 Pb molecular weight marker from Thermo Scientific Fermentas® (Fermentas UAB, subsidiary of Thermo Fisher Scientific Inc. Vilnius, Lithuania). Details about the genotyping procedure were previously reported [5, 42]. MOI referred to the number of different parasite genotypes co-existing within a given infection. MOI was calculated separately for msp1 and msp2 and the final MOI value for each clinical isolate represented the maximum MOI value from both msp1 and msp2 loci as previously described in detail [5, 9]. A sample was classified as belonging to a given allelic family (K1, MAD20 or RO33 for msp1 and 3D7 or FC27 for msp2) on the basis of an occurrence of at least one band after DNA amplification using the family specific primers [5, 9, 44].

Statistical analysis

Clinical data were entered in an ACCESS database while molecular data were captured in a separate EXCEL file. All statistical analyses were performed with the R (version 3.5.1) software [45]. A generalized linear model (GLM) with binomial errors was used to investigate the effect of MOI, parasite density, patient age, gender, Hb level, body temperature, and two-way interactions on gametocyte carriage. GLMs with Poisson, Gaussian and negative binomials errors were used to compare, respectively, (i) MOI levels; (ii) Hb levels; and, (iii) the asexual stage density, between gametocyte-positive and -negative samples. The effect of the presence of a particular allelic family (msp1: K1, MAD20, RO33, and msp2: 3D7 and FC27) on gametocyte carriage was tested using a GLM with binomial errors. Finally, a GLM with negative binomial errors was used to investigate the effect of MOI, parasite density, patient age, gender, Hb level, body temperature, and two-way interactions on gametocytaemia. Model simplification used stepwise removal of terms, followed by likelihood ratio tests (LRT). Term removals that significantly reduced explanatory power (p < 0.05) were retained in the minimal adequate model. Odds ratios and 95% confidence interval (CI) of the binomial models were also computed. Collinearity between explanatory variables was checked using the “vif” function of the “car” package.

Results

Of the 724 samples positive to P. falciparum and successfully genotyped, gametocytes were found along with asexual forms in 48 samples, i.e., a gametocyte carriage prevalence was 6.63% (48/724) with a mean gametocyte density of 186.9 ± 35.1 gametocytes/µl (min = 15, max = 1,240). This represents a minimum gametocyte prevalence since the sensitivity of microscopy is relatively low. All 724 positives samples were P. falciparum mono-infection.

Effect of MOI on gametocyte carriage

The probability of gametocyte carriage significantly increased with increasing values of MOI (GLM binomial, LRT X21 = 5, P = 0.025, Fig. 1a, Additional file 1: Table S1). In particular, one additional unit of MOI enhanced the odds of being gametocyte-positive by 1.26 (95% CI: 1.04–1.5). Accordingly, the number of parasite genotypes identified in blood samples from gametocyte carriers (median = 3.5) was significantly higher than in samples from gametocyte-negative patients (median = 2) (GLM Poisson: LRT X21 = 4.78, P = 0.029, Fig. 1b).

Fig. 1
figure 1

Effects of MOI and Hb level (g/dL) on gametocyte carriage. a Estimated probability of gametocyte carriage (± 95% confidence interval) as a function of MOI level. Gametocyte stages are more likely to occur in parasite isolates harbouring high number of P. falciparum genotypes. b MOI levels in gametocyte-negative (0) and –positive (1) patients. Each point represents a blood sample from a P. falciparum-infected patient (n = 676 gametocyte negative patients and 48 gametocyte carriers). Dark blue points indicate samples with low values of MOI while light blue points show samples with high values of MOI. The horizontal bold black line represents the median value of MOI for each of the two group, and the upper and lower boundaries of the box indicate the 75th and 25th percentile, respectively. Note that the 25th percentile of the gametocyte-negative group (0) was equal to its median (i.e., 2). The presence of gametocytes in patient blood was determined by microscopic observation and MOI was defined as the number of different parasite genotypes (based on msp1 and msp2 genotyping) co-existing within a particular infection. c Estimated probability of gametocyte carriage (± 95% confidence interval) as a function of Hb level. Gametocyte stages are more likely to occur in anaemic patients. d Hb levels in gametocyte-negative (0) and –positive (1) patients. Each point represents a blood sample from a P. falciparum-infected patient. There were 3 Hb missing values giving a total of 673 gametocyte-negative patients and 48 gametocyte carriers. The presence of gametocytes in patient blood was determined by microscopic observation and Hb level was measured using Hemocue® 301 + 

Effect of Hb on gametocyte carriage

The probability of gametocyte carriage decreased with Hb level (GLM binomial, LRT X21 = 5.01, P = 0.0245, Fig. 1c, Additional file 1: Table S1). In particular, for every 1-point increase in Hb level, the odds of being gametocyte-positive decreased by 19% (OR = 0.81, 95% CI: 0.70–0.94). Accordingly, Hb level was lower in gametocyte carriers than in patients infected with asexual stages only (GLM Gaussian, F1,719 = 6.94, P = 0.009, Fig. 1d). There was no collinearity between Hb level and MOI (vif = 1.003) such that the effect of anaemia on gametocyte carriage was independent of the effect of MOI.

Effect of parasite density on gametocyte carriage and interaction with MOI and body temperature

There was a negative association between the presence of gametocyte and parasite density (GLM Binomial, LRT X21 = 15, P < 0.001, Fig. 2a, Additional file 1: Table S1) (OR = 0.99, 0.99–1). Accordingly, gametocyte carriers exhibited lower parasite densities compared to gametocyte-negative patients (GLM negative binomial: LRT X21 = 23.8, P < 0.001, Fig. 2b). There was a significant interaction between parasite density and MOI (LRT X21 = 5.4, P = 0.02, Additional file 1: Table S1): MOI decreased with parasite density in gametocyte-negative patients, but increased in gametocyte carriers (Fig. 2c). There was also an interaction between parasite density and body temperature on gametocyte carriage (LRT X21 = 6, P = 0.015, Additional file 1: Table S1): body temperature increased with parasite density in gametocyte-negative patients but tended to decrease in gametocyte carriers (Fig. 2d).

Fig. 2
figure 2

Relationship between parasite density and gametocyte carriage. a Estimated probability of gametocyte carriage as a function of parasite density. Gametocyte stages are more likely to occur in infections with low parasite density. b Parasite density (number of parasites / µl of blood) in gametocyte-negative (0) and –positive (1) patients. Each point represents a blood sample from a P. falciparum-infected patient (n = 676 gametocyte negative patients and 48 gametocyte carriers). The horizontal bold black line represents the median value of parasite density for each of the two groups, and the upper and lower boundaries of the box indicate the 75th and 25th percentile, respectively. c MOI level as a function of parasite density for both gametocyte-negative (0) and –positive (1) patients. Each colour line represents a linear relationship (± se) fitted to the MOI values for each group (red = non-carriers and blue = gametocyte carriers). The presence of gametocytes in patient blood was determined by microscopic observation and MOI was defined as the number of different parasite genotypes (based on msp1 and msp2 genotyping) co-existing within a particular infection. d Body temperature (°C) as a function of parasite density for both gametocyte-negative (in red) and –positive (in blue) patients. The lines represent a linear relationship (± se) fitted to the temperature values for each group. Because there were 52 missing values in the temperature dataset temperature, the sample size in d is different from that of (a) and (c) (n = 625 gametocyte negative patients and 47 gametocyte carriers)

Effect of age, gender and body temperature on gametocyte carriage

There was no main effect of patient gender, age, body temperature and other two-way interactions on gametocyte carriage prevalence (Additional file 1: Table S1).

Effect of msp1 and msp2 allelic families on gametocyte carriage

The genetic composition of the infection influenced gametocyte prevalence. In particular, the presence of RO33 increased the odds of developing gametocytes by 2 (95% CI: 1.09–4.02, Fig. 3). The other allelic families K1, MAD20, FC27, and 3D7 had no significant impact on the occurrence of gametocytes in infected patients (Fig. 3).

Fig. 3
figure 3

Relationship between gametocyte carriage and parasite genetic composition. Gametocyte prevalence in samples containing (“yes”) or lacking (“no”) parasites genotypes from the a K1 allelic family, b MAD20 allelic family, c RO33 allelic family, d 3D7 allelic family, and e FC27 allelic family. The horizontal dashed line indicates the overall prevalence in the studied population (0.066). Statistical influence of the presence/absence of a given allelic family is showed using the odds ratio (OR) from the binomial model and the associated 95% confidence interval and P-value

Effects of MOI, Hb, parasite density, age, gender, and body temperature on gametocyte density

Although gametocyte density tended to increase with MOI (Fig. 4a) and decrease with Hb (Fig. 4b), similar to the observation made for gametocyte carriage. This was not statistically significant, likely due to a lack of power (i.e., only 48 gametocyte-positive samples) (MOI effect: LRT X21 = 418, P = 0.12; Hb effect: LRT X21 = 419, P = 0.12, Additional file 2: Table S2). The only significant predictor of gametocyte density was asexual stage density (LRT X21 = 2090, P < 0.001, Additional file 2: Table S2, Fig. 4c). However, there was no effect of gender, age, body temperature, and two-way interactions (Additional file 2: Table S2).

Fig. 4
figure 4

Relationship between gametocytaemia and a MOI, b Hb level and c asexual stage density. Blue lines represent a linear relationship (± se) fitted to the MOI, Hb and parasite density values, respectively. On each panel, gametocytaemia on the y-axis is on a log10 scale. The x-axis of c is on a log10 scale

Discussion

In this study, the probability of gametocyte carriage significantly increased with increasing values of MOI. Given the within-host competition among different parasite strains [5], the presence of multiple, genetically distinct, parasite strains may represent unfavourable environment for the less competitive strains. The stress resulting from this unfavourable environment was described as a potential mechanism that triggers the sexual differentiation pathway [13, 46]. Previous reports indicated that MOI can promote either longer persistence or continuous production of gametocytes [23, 25, 26]. Beside supporting the terminal investment strategy, the positive association between MOI and gametocyte carriage may suggest that infections with multiple clones have simply more chance to contain some clones that will evade the host immune response, persist in the host and result in gametocyte development [24, 26].

Results of this study support two predictions made by a model of transmission strategies in the rodent malaria parasite Plasmodium chabaudi, that is: investment in gametocytogenesis should increase: (i) at low parasite densities; and, (ii) with host anaemia [47]. This study showed that gametocyte carriers exhibited lower parasite densities compared to gametocyte-negative patients. This negative association between parasite density and gametocyte carriage was previously reported highlighting lower parasitaemia as a risk factor of gametocytaemia. As lower parasitaemia are commonly observed in asymptomatic population, this finding raises the question about the contribution of asymptomatic population in sustainability of malaria transmission. Thus, restricting this investigation to symptomatic patients represents a limit of the study. Indeed, higher prevalence of gametocyte carriers was previously reported in asymptomatic population than in symptomatic population [48, 49]. This could mean that the principal infectious reservoir could be people outside the vulnerable groups, i.e., children and pregnant women who fall regularly clinically ill from the disease. In such context, interventions targeting specifically the reservoir need to be implemented in order to achieve best impact in malaria control and elimination. However, conflicting findings pointing to hyper-parasitaemia as risk factor of gametocytes carriage was also reported [50]. The differences are mainly attributable to population profile (symptomatic versus asymptomatic including sub-patent parasite carriers) and parasite detection method, i.e., light microscopy against qPCR [51]. Gametocyte detection based on light microscopy, rather than sensitive methods such as qPCR, represent one limitation of this study, underestimating gametocyte carrier prevalence as well as the possible contribution of other sub-patent Plasmodium species. In addition, reading malaria slides against 500 leucocytes rather than 1,000 leucocytes for gametocyte detection could have underestimated detection power.

In gametocyte-negative patients, MOI was lower in patients with heavy infections than in patients with mild infections, thus supporting the existence of within host competition with suppressive effect for disadvantaged strains [5]. However, in gametocyte-positive patients, MOI increased with increasing parasite density suggesting that the suppressive effect resulting from within-host competition and strain-specific immunity (observed in gametocyte-negative population) is mitigated and the disadvantaged strains commit to sexual differentiation instead of being eliminated.

Furthermore, as pointed out by several other authors, this study identified host anaemia as a risk factor for gametocyte carriage, such that for every 1-point increase in Hb level, the odds of being gametocyte-positive decreased by 19% [31, 48, 52]. Possible mechanism of this correlation between anaemia and gametocyte carriage would be the haemolysis of infected erythrocytes which was previously identified as a factor inducing gametocytogenesis [53]. These results suggest that combining malaria and malnutrition programmes would enhance the effectiveness of control interventions.

There was no main effect of patient gender, age, body temperature, and other two-way interactions on gametocyte carriage, yet some of these host factors such as gender and body temperature were previously described as associated risk factors of gametocytaemia [35, 48]. This may also be attributable to population profile because of the lack of asymptomatic population in this study.

This study used msp1 and msp2 as polymorphic markers which are under immune pressure due to their involvement in the merozoite invasion [54]. This represents one limitation of this study in terms of discriminatory power compared to neutral markers such as microsatellites. New methods with much higher genetic resolution such as Next Generation Sequencing would also provide better characterization of the parasite population profile compared to msp1 and msp2 genotyping. Nevertheless, the genetic composition of the infection influenced gametocyte prevalence. Previous reports suggested that different P. falciparum parasite variants may influence gametocyte carriage differently. A study reported that individuals carrying both 3D7 and FC27 allelic types of parasite had a higher risk of harbouring gametocytes compared to single 3D7 or FC27 allelic types [26]. In addition, an association between gametocytaemia and particular dhfr allele was previously reported [55]. In this study, the RO33 allelic family was mostly associated with gametocyte carriage than others. Nevertheless, these findings, although highlighting an association with gametocyte carriage do not establish causal relationship, i.e., these allelic families are mostly engaged in sexual differentiation than others. For example, if this was true, the RO33 allelic family should have been the commonest allelic family (mostly transmitted) yet it is the less represented in the study area [9]. Therefore, further investigations exploring the effect of allelic family variability on gametocytaemia with special attention towards the infectivity of RO33 allelic family for anopheline mosquitoes are needed.

Conclusion

This study provides insight into potential factors influencing gametocyte production in symptomatic patients. The findings contribute to enhance understanding of risk factors associated with gametocytes carriage in humans.

Availability of data and materials

Data supporting the conclusions of this article are included within the article and its additional files.

Abbreviations

dhfr:

Dihydrofolate reductase

DNA:

Deoxyribonucleic acid

GLM:

Generalized linear model

GPM:

Global precipitation measurement

LRT:

Likelihood ratio test

MOI:

Multiplicity of infection

MSP1:

Merozoite surface protein 1

MSP2:

Merozoite surface protein 2

PCR:

Polymerase chain reaction

References

  1. WHO. World malaria report 2019. Geneva: World Health Organization; 2019.

    Google Scholar 

  2. Arnot D. Clone multiplicity of Plasmodium falciparum infections in individuals exposed to variable levels of disease transmission. Trans R Soc Trop Med Hyg. 1998;92:580–5.

    Article  CAS  PubMed  Google Scholar 

  3. Babiker HA, Ranford-Cartwright LC, Walliker D. Genetic structure and dynamics of Plasmodium falciparum infections in the Kilombero region of Tanzania. Trans R Soc Trop Med Hyg. 1999;93:11–4.

    Article  PubMed  Google Scholar 

  4. Soulama I, Nébié I, Ouédraogo 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  PubMed Central  CAS  Google Scholar 

  5. Sondo P, Derra K, Lefevre T, Diallo-Nakanabo S, Tarnagda Z, Zampa O, et al. Genetically diverse Plasmodium falciparum infections, within-host competition and symptomatic malaria in humans. Sci Rep. 2019;9:127.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Nwakanma D, Kheir A, Sowa M, Dunyo S, Jawara M, Pinder M, et al. High gametocyte complexity and mosquito infectivity of Plasmodium falciparum in the Gambia. Int J Parasitol. 2008;38:219–27.

    Article  PubMed  Google Scholar 

  7. Walliker D, Quakyi IA, Wellems TE, McCutchan TF, Szarfman A, London WT, et al. Genetic analysis of the human malaria parasite Plasmodium falciparum. Science. 1987;236:1661–6.

    Article  CAS  PubMed  Google Scholar 

  8. Tanabe K, Sakihama N, Kaneko O, Saito-Ito A, Kimura M. A PCR method for molecular epidemiology of Plasmodium falciparum msp-1. Tokai J Exp Clin Med. 1998;23:375–82.

    CAS  PubMed  Google Scholar 

  9. Sondo P, Derra K, Rouamba T, Diallo SN, Taconet P, Kazienga A, et al. Determinants of Plasmodium falciparum multiplicity of infection and genetic diversity in Burkina Faso. Parasit Vectors. 2020;13:427.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Konaté L, Zwetyenga J, Rogier C, Bischoff E, Fontenille D, Tall A, et al. Variation of Plasmodium falciparum msp1 block 2 and msp2 allele prevalence and of infection complexity in two neighbouring Senegalese villages with different transmission conditions. Trans R Soc Trop Med Hyg. 1999;93:21–8.

    Article  PubMed  Google Scholar 

  11. Beck S, Mockenhaupt FP, Bienzle U, Eggelte TA, Thompson W, Stark K. Multiplicity of Plasmodium falciparum infection in pregnancy. Am J Trop Med Hyg. 2001;65:631–6.

    Article  CAS  PubMed  Google Scholar 

  12. Magesa S, Mdira K, Babiker H, Alifrangis M, Färnert A, Simonsen P, et al. Diversity of Plasmodium falciparum clones infecting children living in a holoendemic area in north-eastern Tanzania. Acta Trop. 2002;84:83–92.

    Article  CAS  PubMed  Google Scholar 

  13. Dyer M, Day K. Commitment to gametocytogenesis in Plasmodium falciparum. Parasitol Today. 2000;16:102–7.

    Article  CAS  PubMed  Google Scholar 

  14. Josling GA, Llinás M. Sexual development in Plasmodium parasites: knowing when it’s time to commit. Nat Rev Microbiol. 2015;13:573–87.

    Article  CAS  PubMed  Google Scholar 

  15. Carter LM, Kafsack BF, Llinás M, Mideo N, Pollitt LC, Reece SE. Stress and sex in malaria parasites: why does commitment vary? Evol Med Public Health. 2013;2013:135–47.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Taylor LH, Read AF. Why so few transmission stages? Reproductive restraint by malaria parasites. Parasitol Today. 1997;13:135–40.

    Article  CAS  PubMed  Google Scholar 

  17. De Roode JC, Helinski ME, Anwar MA, Read AF. Dynamics of multiple infection and within-host competition in genetically diverse malaria infections. Am Nat. 2005;166:531–42.

    Article  PubMed  Google Scholar 

  18. Koella JC, Antia R. Optimal pattern of replication and transmission for parasites with two stages in their life cycle. Theor Popul Biol. 1995;47:277–91.

    Article  Google Scholar 

  19. Antia R, Nowak MA, Anderson RM. Antigenic variation and the within-host dynamics of parasites. Proc Natl Acad Sci USA. 1996;93:985–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Mideo N. Parasite adaptations to within-host competition. Trends Parasitol. 2009;25:261–8.

    Article  PubMed  Google Scholar 

  21. Mideo N, Day T. On the evolution of reproductive restraint in malaria. Proc Biol Sci. 2008;275:1217–24.

    PubMed  PubMed Central  Google Scholar 

  22. Williams G. Adaptation and natural selection. Reprint. New Jersey. Princeton Science Library: Princeton; 1966.

    Google Scholar 

  23. Wampfler R, Timinao L, Beck H-P, Soulama I, Tiono AB, Siba P, et al. Novel genotyping tools for investigating transmission dynamics of Plasmodium falciparum. J Infect Dis. 2014;210:1188–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Nassir E, Abdel-Muhsin A-MA, Suliaman S, Kenyon F, Kheir A, Geha H, et al. Impact of genetic complexity on longevity and gametocytogenesis of Plasmodium falciparum during the dry and transmission-free season of eastern Sudan. Int J Parasitol. 2005;35:49–55.

    Article  CAS  PubMed  Google Scholar 

  25. Vardo-Zalik AM, Schall JJ. Clonal diversity alters the infection dynamics of a malaria parasite (Plasmodium mexicanum) in its vertebrate host. Ecology. 2009;90:529–36.

    Article  PubMed  Google Scholar 

  26. Lamptey H, Ofori MF, Kusi KA, Adu B, Owusu-Yeboa E, Kyei-Baafour E, et al. The prevalence of submicroscopic Plasmodium falciparum gametocyte carriage and multiplicity of infection in children, pregnant women and adults in a low malaria transmission area in Southern Ghana. Malar J. 2018;17:331.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Pollitt LC, Mideo N, Drew DR, Schneider P, Colegrave N, Reece SE. Competition and the evolution of reproductive restraint in malaria parasites. Am Nat. 2011;177:358–67.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Huijben S, Sim DG, Nelson W, Read AF. The fitness of drug-resistant malaria parasites in a rodent model: multiplicity of infection. J Evol Biol. 2011;24:2410–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Greischar MA, Mideo N, Read AF, Bjørnstad ON. Predicting optimal transmission investment in malaria parasites. Evolution. 2016;70:1542–58.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Buckling A, Ranford-Cartwright L, Miles A, Read A. Chloroquine increases Plasmodium falciparum gametocytogenesis in vitro. Parasitology. 1999;118:339–46.

    Article  CAS  PubMed  Google Scholar 

  31. Nacher M, Singhasivanon P, Silachamroon U, Treeprasertsuk S, Tosukhowong T, Vannaphan S, et al. Decreased hemoglobin concentrations, hyperparasitemia, and severe malaria are associated with increased Plasmodium falciparum gametocyte carriage. J Parasitol. 2002;88:97–101.

    Article  PubMed  Google Scholar 

  32. Paul R, Lafond T, MĂĽller-Graf C, Nithiuthai S, Brey P, Koella J. Experimental evaluation of the relationship between lethal or non-lethal virulence and transmission success in malaria parasite infections. BMC Evol Biol. 2004;4:30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Reece SE, Duncan AB, West SA, Read AF. Host cell preference and variable transmission strategies in malaria parasites. Proc Biol Sci. 2005;272:511–7.

    PubMed  PubMed Central  Google Scholar 

  34. Ali E, Mackinnon MJ, Abdel-Muhsin A-MA, Ahmed S, Walliker D, Babiker HA. Increased density but not prevalence of gametocytes following drug treatment of Plasmodium falciparum. Trans R Soc Trop Med Hyg. 2006;100:176–83.

    Article  CAS  PubMed  Google Scholar 

  35. Sowunmi A, Fateye B, Adedeji A, Fehintola F, Happi T. Risk factors for gametocyte carriage in uncomplicated falciparum malaria in children. Parasitology. 2004;129:255.

    Article  CAS  PubMed  Google Scholar 

  36. Janse C, Boorsma E, Ramesar J, Van Vianen P, Van der Meer R, Zenobi P, et al. Plasmodium berghei: gametocyte production, DNA content, and chromosome-size polymorphisms during asexual multiplication in vivo. Exp Parasitol. 1989;68:274–82.

    Article  CAS  PubMed  Google Scholar 

  37. Graves P, Carter R, McNeill K. Gametocyte production in cloned lines of Plasmodium falciparum. Am J Trop Med Hyg. 1984;33:1045–50.

    Article  CAS  PubMed  Google Scholar 

  38. Diebner HH, Eichner M, Molineaux L, Collins WE, Jeffery GM, Dietz K. Modelling the transition of asexual blood stages of Plasmodium falciparum to gametocytes. J Theor Biol. 2000;202:113–27.

    Article  CAS  PubMed  Google Scholar 

  39. Abdel-Wahab A, Abdel-Muhsin A-MA, Ali E, Suleiman S, Ahmed S, Walliker D, et al. Dynamics of gametocytes among Plasmodium falciparum clones in natural infections in an area of highly seasonal transmission. J Infect Dis. 2002;185:1838–42.

    Article  PubMed  Google Scholar 

  40. Casaglia O, Dore E, Frontali C, Zenobi P, Walliker D. Re-examination of earlier work on repetitive DNA and mosquito infectivity in rodent malaria. Mol Biochem Parasitol. 1985;16:35–42.

    Article  CAS  PubMed  Google Scholar 

  41. Bousema T, Drakeley C. Epidemiology and infectivity of Plasmodium falciparum and Plasmodium vivax gametocytes in relation to malaria control and elimination. Clin Microbiol Rev. 2011;24:377–410.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Sondo P, Derra K, Diallo-Nakanabo S, Tarnagda Z, Zampa O, Kazienga A, et al. Effectiveness and safety of artemether–lumefantrine versus artesunate–amodiaquine for unsupervised treatment of uncomplicated falciparum malaria in patients of all age groups in Nanoro, Burkina Faso: a randomized open label trial. Malar J. 2015;14:325.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Derra K, Rouamba E, Kazienga A, Ouedraogo S, Tahita MC, Sorgho H, et al. Profile: Nanoro health and demographic surveillance system. Int J Epidemiol. 2012;41:1293–301.

    Article  PubMed  Google Scholar 

  44. Snounou G, Zhu X, Siripoon N, Jarra W, Thaithong S, Brown KN, et al. Biased distribution of msp1 and msp2 allelic variants in Plasmodium falciparum populations in Thailand. Trans R Soc Trop Med Hyg. 1999;93:369–74.

    Article  CAS  PubMed  Google Scholar 

  45. Crawley M. The R Book. 2nd ed. Chichester: Wiley; 2007.

    Book  Google Scholar 

  46. Baker DA. Malaria gametocytogenesis. Mol Biochem Parasitol. 2010;172:57–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Cameron A, Reece SE, Drew DR, Haydon DT, Yates AJ. Plasticity in transmission strategies of the malaria parasite, Plasmodium chabaudi: environmental and genetic effects. Evol Appl. 2013;6:365–76.

    Article  PubMed  Google Scholar 

  48. von Seidlein L, Drakeley C, Greenwood B, Walraven G, Targett G. Risk factors for gametocyte carriage in Gambian children. Am J Trop Med Hyg. 2001;65:523–7.

    Article  Google Scholar 

  49. Lindblade KA, Steinhardt L, Samuels A, Kachur SP, Slutsker L. The silent threat: asymptomatic parasitemia and malaria transmission. Expert Rev Anti Infect Ther. 2013;11:623–39.

    Article  CAS  PubMed  Google Scholar 

  50. Koepfli C, Robinson LJ, Rarau P, Salib M, Sambale N, Wampfler R, et al. Blood-stage parasitaemia and age determine Plasmodium falciparum and P. vivaxgametocytaemia in Papua New Guinea. PLoS One. 2015;10:e0126747.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  51. Ayanful-Torgby R, Oppong A, Abankwa J, Acquah F, Williamson KC, Amoah LE. Plasmodium falciparum genotype and gametocyte prevalence in children with uncomplicated malaria in coastal Ghana. Malar J. 2016;15:592.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Price R, Nosten F, Simpson JA, Luxemburger C, Phaipun L, Ter Kuile F, et al. Risk factors for gametocyte carriage in uncomplicated falciparum malaria. Am J Trop Med Hyg. 1999;60:1019–23.

    Article  CAS  PubMed  Google Scholar 

  53. Schneweis S, Maier W, Seitz H. Haemolysis of infected erythrocytes—a trigger for formation of Plasmodium falciparum gametocytes? Parasitol Res. 1991;77:458–60.

    Article  CAS  PubMed  Google Scholar 

  54. Blackman MJ, Heidrich H-G, Donachie S, McBride J, Holder A. 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 

  55. Mockenhaupt FP, Teun Bousema J, Eggelte TA, Schreiber J, Ehrhardt S, Wassilew N, et al. Plasmodium falciparum dhfr but not dhps mutations associated with sulphadoxine-pyrimethamine treatment failure and gametocyte carriage in northern Ghana. Trop Med Int Health. 2005;10:901–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Special thanks to the local community and health workers of the Nanoro health district. We thank SANOFI Aventis, WHO-TDR and the Institute of Tropical Medicine of Antwerp, Belgium for the financial support allocated to the pharmacovigilance project and the molecular analyses.

Funding

World Health Organization—the Special Program for Research and Training in Tropical Diseases and SANOFI AVENTIS supported the Pharmacovigilance trial (NCT01232530) from which clinical samples were collected. The Institute for Tropical Medicine (ITM), Belgium supported through the FA3-DGCD program supported the molecular analyses.

Author information

Authors and Affiliations

Authors

Contributions

SP, IV, DK, ZT, HS, TH designed the study, TH, SP, DK, ZT, SND, HS, BB, TMC, TH implemented the study and supervised field work, SP, TMC supervised laboratory work analyses, SP, DK, AK, contributed in data management, IH, BB, AK, LT contributed in statistical analysis, SP, LT, IH, TMC, BB, TH, contributed in drafting the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Paul Sondo.

Ethics declarations

Ethics approval and consent to participate

The study was carried out in full compliance with fundamental ethical principles and relevant regulations in Burkina Faso. The pharmacovigilance study protocol was approved by the Ethical Committee for health research in Burkina Faso (Deliberation N°2010-27, April 14th; 2010). A signed/thumb printed inform consent was obtained from all participants before enrolment.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Determinants of the frequency of the K1 allelic family.

Additional file 2: Table S2.

Determinants of the frequency of the MAD20 allelic family.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sondo, P., Bihoun, B., Tahita, M.C. et al. Plasmodium falciparum gametocyte carriage in symptomatic patients shows significant association with genetically diverse infections, anaemia, and asexual stage density. Malar J 20, 31 (2021). https://doi.org/10.1186/s12936-020-03559-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-020-03559-0

Keywords