Skip to main content

Age-structured gametocyte allocation links immunity to epidemiology in malaria parasites



Despite a long history of attempts to model malaria epidemiology, the over-riding conclusion is that a detailed understanding of host-parasite interactions leading to immunity is required. It is still not known what governs the duration of an infection and how within-human parasite dynamics relate to malaria epidemiology.

Presentation of the hypothesis

Immunity to Plasmodium falciparum develops slowly and requires repeated exposure to the parasite, which thus generates age-structure in the host-parasite interaction. An age-structured degree of immunity would present the parasite with humans of highly variable quality. Evolutionary theory suggests that natural selection will mould adaptive phenotypes that are more precise (less variant) in "high quality" habitats, where lifetime reproductive success is best. Variability in malaria parasite gametocyte density is predicted to be less variable in those age groups who best infect mosquitoes. Thus, the extent to which variation in gametocyte density is a simple parasite phenotype reflecting the complex within-host parasite dynamics is addressed.

Testing the hypothesis

Gametocyte densities and corresponding infectiousness to mosquitoes from published data sets and studies in both rural and urban Cameroon are analysed. The mean and variation in gametocyte density according to age group are considered and compared with transmission success (proportion of mosquitoes infected). Across a wide range of settings endemic for malaria, the age group that infected most mosquitoes had the least variation in gametocyte density, i.e. there was a significant relationship between the variance rather than the mean gametocyte density and age-specific parasite transmission success. In these settings, the acquisition of immunity over time was evident as a decrease in asexual parasite densities with age. By contrast, in an urban setting, there were no such age-structured relationships either with variation in gametocyte density or asexual parasite density.

Implications of the hypothesis

Gametocyte production is seemingly predicted by evolutionary theory, insofar as a reproductive phenotype (gametocyte density) is most precisely expressed (i.e. is most invariant) in the most infectious human age group. This human age group would thus be expected to be the habitat most suitable for the parasite. Comprehension of the immuno-epidemiology of malaria, a requisite for any vaccine strategies, remains poor. Immunological characterization of the human population stratified by parasite gametocyte allocation would be a step forward in identifying the salient immunological pathways of what makes a human a good habitat.


Mathematical models of infectious diseases and most especially those based on the calculation of R0, the basic reproductive number, have proved to be very powerful and robust in interpreting epidemiological trends, most notably for viral and macro-parasitic diseases [1]. However, despite a long history of attempts to model malaria epidemiology, and most especially that of lethal human malaria Plasmodium falciparum, the over-riding conclusion is that a detailed understanding of host-parasite interactions leading to immunity is required. The parasite's high degree of genetic diversity and mechanisms of clonal antigenic variation highlight the complex nature of the human-parasite interaction. However, it is still not known what governs the duration of an infection and how within-human parasite dynamics relate to malaria epidemiology. Indeed, controversy still rages over the rate of acquisition of immunity to malaria [2, 3].

The classical Ross-Macdonald model ( R 0 = m a 2 b c μ γ e x p μ T MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGsbGudaWgaaWcbaGaeGimaadabeaakiabg2da9maalaaabaGaemyBa0Maemyyae2aaWbaaSqabeaacqaIYaGmaaGccqWGIbGycqWGJbWyaeaaiiGacqWF8oqBcqWFZoWzaaacbiGae4xzauMae4hEaGNae4hCaa3aaWbaaSqabeaacqGHsislcqWF8oqBcqWGubavaaaaaa@4212@ ) captures the essence of malaria epidemiology and consists of mosquito components that determine, to a large extent, the force of infection [mosquito density (m), biting rate (a) and mortality rate (μ) and the duration of sporogonic development in the vector (T)] and human components that directly relate to the human-parasite interaction. Two of these parameters are concerned with within-host parasite dynamics: the human recovery rate from infection (γ) and the infectiousness of an infected individual to mosquitoes (c). Parameter (b) is the proportion of infective bites that lead to an infection in man.

Immunity to P. falciparum develops slowly and requires repeated exposure to the parasite. This generates age-structure in the host-parasite interaction and most notably in the duration of episodes of infection [4]. In areas endemic and stable for malaria, recovery rate is lowest in the young/intermediate age groups [5, 6]. Recent studies [2, 7], however, suggest that the rate of absolute recovery from infection may not increase with age, even though infecting parasite density may be reduced. Epidemiological models incorporating age-dependency, whether through increased absolute recovery from infection or increased parasite death rate, improve markedly the fit with reality [1]. Measures of human infectiousness to mosquitoes have also shown strong age-specificity such that infections in particular age groups most successfully infect mosquitoes [811].

In this work, an attempt is made to identify key measurable parameters linking infection and immunity to epidemiology.

Presentation of the hypothesis

Age-dependency in the development of an immune response against the asexual parasite stages would present the parasite with habitats of highly variable quality. Evolutionary theory predicts that, in a spatially heterogeneous environment, natural selection will mould adaptive phenotypes such that they are more precise (less variant) in the better and the more frequent habitats and less so in infrequent habitats and those where reproductive performance is poor [12], i.e. natural selection will act on both the mean and the variance of the phenotype in question. Heterogeneity in habitat quality can severely alter R0 and it is trivial to note that optimising transmission (c) in the age groups where recovery rate (γ) is the lowest will significantly increase R0. If the parasite is indeed optimising transmission according to the immune-dependent quality of individuals, this should be apparent from the within-host parasite dynamics and notably that pertaining to gametocyte production.

Testing the hypothesis

Gametocyte densities and corresponding infectiousness to mosquitoes from published data sets [9, 10, 13] and studies in both rural [11, 14] and urban [15] settings in Cameroon are examined. The mean and variation in gametocyte density (as proxies of the parasite's transmission phenotype) according to age group are considered and compared with actual measures of transmission success (the proportion of mosquitoes that become infected). Maturation of gametocytes occurs slowly and infections with apparent gametocytes may not be infectious at that exact time but may be so hours later. Therefore, all observed gametocyte densities are included. The age groups chosen in the published and the rural Cameroon data sets reflect the classical age distribution as defined by malaria prevalence rates that provide a gross measure of the force of infection and gradual development of immunity. The decrease in asexual parasite densities with age is a more sensitive indication of the slow acquisition of immunity that generates age structure and is thus presented. By contrast, malaria transmission in urban settings is low and unstable, resulting in a poor acquisition of immunity. Therefore, age-stratification of the urban Cameroon data set would not be expected to reflect age-specific development of immunity or be apparent from asexual parasite densities.

Cameroon study sites and data collection

The study sites and methods have been previously published [11, 14, 15]. Briefly, in the rural setting, individuals, symptomatic or not, were recruited in two adjacent villages in the district of Mengang, South Cameroon, in an area of degraded forest with seasonal transmission (intensity of 170 infected bites/person per year) and hyperendemic malaria (P. falciparum prevalence rates range from 49% to 82%, depending on the season). Mosquitoes were gorged directly on individuals whether positive or not for gametocytes. In the urban setting, patients were recruited at a Public Health Centre in central Yaounde, the capital town of Cameroon. From those individuals positive for gametocytes by thick blood smear, venous blood was immediately placed into a membrane feeder upon which mosquitoes gorged. Previous comparative studies have shown that although direct feeding of mosquitoes on infectious individuals yields higher infections in mosquitoes than feeding through a membrane, there is good concordance between the two tests for both mosquito infection percentages and oocyst loads [16, 17]. In both studies, only experiments where at least 20 mosquitoes survived to day of dissection for oocyst counts (day 7 post-gorging) were used in analyses. The mean number of mosquitoes dissected per experiment was 35 and 36 in the rural and urban settings respectively. Blood parasite counts were established by thick smear per 1,000 leucocytes for asexual parasite density and 1,000 and 3,000 leucocytes for gametocytes respectively in urban and rural settings, assuming an average number of 8,000 leucocytes per microlitre of blood. The studies used the same laboratory mosquito strain reared in the same insectary (Anopheles gambiae s.s. Yaounde strain; OCEAC insectary).

Statistical analyses

The effects of age group and gametocyte density per infection on transmission success in the Cameroon data sets were analysed in Genstat version 7 using logistic regression specifying a binomial error structure. The effect of age group on asexual parasite density was analysed by loglinear regression fitting a GLM, specifying a Poisson error structure. Dispersion parameters were estimated and thus the analyses generate an F statistic in the analysis of deviance.

Using data from the Cameroon studies as well as the two published studies, how the mean and the variability in gametocyte density per age group affected transmission success were then examined. This combined data set was analysed by two alternative methods: a REML meta-analysis and by fitting a GLMM (Mixed model) with study specified in the random model and the mean proportion of mosquitoes infected by age group as the response variate. A REML meta-analysis takes into account differences between sites by incorporating the residual variance of each experiment. GLMMs account for variation due to study site and calculate the variance component due to this random factor independently of that due to the factor of interest – i.e. the response variate. Because the measure of variability in gametocyte density per age group was previously given as confidence intervals (CI) in the published studies, this measure was used for the Cameroon data sets. The lower CI subtracted from the upper CI was calculated to generate a value of dispersion in gametocyte density by age group. Both this value and mean gametocyte density per age group were normalized for each study by dividing by the lowest value, thus giving the lower value of 1 per study. Such normalization enables comparison of the general effect of gametocyte density and variability in density by age group without study specific effects, most evident of which are the very different overall gametocyte densities. Both the normalized gametocyte density mean and variability were then fitted as explanatory variates. Because the data were over-dispersed, a dispersion parameter was estimated. Statistical significance is presented as Wald statistics, which approximate to a χ2 distribution. For the Cameroon data sets, variance/mean ratios and the Standardized Morisita index (Ip) were additionally calculated. Sample size influences both confidence intervals and variances. The larger the sample size the smaller these measures will be for a data set of the same intrinsic variability. The Standardized Morisita index (Ip), however, is considered the best measure of dispersion being independent of sample size [18]. To calculate this index, the Morisita's index of dispersion (Id) and two critical values, the Uniform index (Mu) and the clumped index (Mc) are calculated as follows:

I d = n [ x 2 x ( x ) 2 x ] MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqsdaWgaaWcbaGaemizaqgabeaakiabg2da9iabd6gaUnaadmaabaWaaSaaaeaadaaeabqaaiabdIha4naaCaaaleqabaGaeGOmaidaaOGaeyOeI0YaaabqaeaacqWG4baEaSqabeqaniabggHiLdaaleqabeqdcqGHris5aaGcbaWaaeWaaeaadaaeabqaaiabdIha4bWcbeqab0GaeyyeIuoaaOGaayjkaiaawMcaamaaCaaaleqabaGaeGOmaidaaOGaeyOeI0YaaabqaeaacqWG4baEaSqabeqaniabggHiLdaaaaGccaGLBbGaayzxaaaaaa@47C2@
M u = χ .975 2 n + x ( x ) 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGnbqtdaWgaaWcbaGaemyDauhabeaakiabg2da9maalaaabaacciGae83Xdm2aa0baaSqaaiabc6caUiabiMda5iabiEda3iabiwda1aqaaiabikdaYaaakiabgkHiTiabd6gaUjabgUcaRmaaqaeabaGaemiEaGhaleqabeqdcqGHris5aaGcbaWaaeWaaeaadaaeabqaaiabdIha4bWcbeqab0GaeyyeIuoaaOGaayjkaiaawMcaaiabgkHiTiabigdaXaaaaaa@4515@
M c = χ .025 2 n + x ( x ) 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGnbqtdaWgaaWcbaGaem4yamgabeaakiabg2da9maalaaabaacciGae83Xdm2aa0baaSqaaiabc6caUiabicdaWiabikdaYiabiwda1aqaaiabikdaYaaakiabgkHiTiabd6gaUjabgUcaRmaaqaeabaGaemiEaGhaleqabeqdcqGHris5aaGcbaWaaeWaaeaadaaeabqaaiabdIha4bWcbeqab0GaeyyeIuoaaOGaayjkaiaawMcaaiabgkHiTiabigdaXaaaaaa@44D5@

where n is the sample size, x is the number of gametocytes in each individual, χ .975 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFhpWydaqhaaWcbaGaeiOla4IaeGyoaKJaeG4naCJaeGynaudabaGaeGOmaidaaaaa@3361@ and χ .025 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFhpWydaqhaaWcbaGaeiOla4IaeGimaaJaeGOmaiJaeGynaudabaGaeGOmaidaaaaa@3345@ are the values of the chi-squared with (n-1) degrees of freedom that have 97.5% or 2.5% of the area to the right. The Standardized Morisita index (Ip) is then calculated by one of the four following formulae:

(a) when Id ≥ Mc > 1, I p = 0.5 + 0.5 × ( I d M c n M c ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqsdaWgaaWcbaGaemiCaahabeaakiabg2da9iabicdaWiabc6caUiabiwda1iabgUcaRiabicdaWiabc6caUiabiwda1iabgEna0oaabmaabaWaaSaaaeaacqWGjbqsdaWgaaWcbaGaemizaqgabeaakiabgkHiTiabd2eannaaBaaaleaacqWGJbWyaeqaaaGcbaGaemOBa4MaeyOeI0Iaemyta00aaSbaaSqaaiabdogaJbqabaaaaaGccaGLOaGaayzkaaaaaa@45C3@ , (b) when Mc > Id ≥ 1, I p = 0.5 × ( I d 1 M u 1 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqsdaWgaaWcbaGaemiCaahabeaakiabg2da9iabicdaWiabc6caUiabiwda1iabgEna0oaabmaabaWaaSaaaeaacqWGjbqsdaWgaaWcbaGaemizaqgabeaakiabgkHiTiabigdaXaqaaiabd2eannaaBaaaleaacqWG1bqDaeqaaOGaeyOeI0IaeGymaedaaaGaayjkaiaawMcaaaaa@400E@ ,

(c) when 1 > Id > Mu, I p = 0.5 × ( I d 1 M u 1 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqsdaWgaaWcbaGaemiCaahabeaakiabg2da9iabgkHiTiabicdaWiabc6caUiabiwda1iabgEna0oaabmaabaWaaSaaaeaacqWGjbqsdaWgaaWcbaGaemizaqgabeaakiabgkHiTiabigdaXaqaaiabd2eannaaBaaaleaacqWG1bqDaeqaaOGaeyOeI0IaeGymaedaaaGaayjkaiaawMcaaaaa@40FB@ and (d) when 1 > Mu > Id, I p = 0.5 + 0.5 × ( I d M u M u ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqsdaWgaaWcbaGaemiCaahabeaakiabg2da9iabgkHiTiabicdaWiabc6caUiabiwda1iabgUcaRiabicdaWiabc6caUiabiwda1iabgEna0oaabmaabaWaaSaaaeaacqWGjbqsdaWgaaWcbaGaemizaqgabeaakiabgkHiTiabd2eannaaBaaaleaacqWG1bqDaeqaaaGcbaGaemyta00aaSbaaSqaaiabdwha1bqabaaaaaGccaGLOaGaayzkaaaaaa@44A6@ .

The Standardized Morisita index (Ip) ranges from -1 to +1. Random patterns give a value of zero, clumped patterns above zero and uniform patterns below zero. With respect to gametocyte data, the higher the value above zero the more over-dispersed, and thus more variable, are the gametocyte densities.


Mosquito infection rates have been classically shown to correlate positively with individual gametocyte density. This positive correlation was confirmed here for the Cameroon rural data set, including only those individuals who infected at least one mosquito (F1,45 = 13.9 P < 0.001 N = 50). However, gametocyte density per se explained only 16.6% of the variation in transmission success and there was a significant effect of age group on mosquito infection rates (F3,45 = 3.43, P = 0.025 N = 50). Within any of the published data sets in Kenya and Sri Lanka and the rural Cameroon data set, the mean gametocyte densities were relatively constant across age groups (Table 1) and typical of those observed in rural endemic situations across the globe. It is notable that the most infectious age groups in any of these settings were not those with the highest mean gametocyte densities (Table 1). Indeed, mean gametocyte density per age group did not correlate with mosquito infection rates (GLMM χ21 = 1.03, P = 0.31). By contrast, it is striking that the age groups that had the lowest variability in gametocyte density infected the highest proportion of mosquitoes (Figure 1) (REML meta-analysis: χ21 = 7.95, P = 0.005; GLMM χ21 = 6.72, P = 0.01). Variation in gametocyte density (i.e. gametocyte precision) explained 17.1% of variation in age group transmission success. Using the rural Cameroon data set, this relationship was confirmed when measuring gametocyte variability using several different dispersion indices (Table 1), including the Standardized Morisita index (Ip), which is considered independent of sample size [18]; all dispersion indices show the same pattern in the Cameroon data set. Moreover, Table 1 shows that, with the exception of P. vivax gametocyte data from Sri Lanka, gametocyte variability is lower in the best transmitting age group rather than the age group with the largest sample size. Larger sample sizes would be expected to yield lower values of dispersion if intrinsic variability were the same among groups and could therefore explain degrees of dispersion. Here this is not the case.

Figure 1

Age group variation in gametocyte density with corresponding infectiousness to mosquitoes. Data shown in Table 1. The higher minus the lower 95% confidence intervals are taken as a measure of gametocyte density variation and are normalized within each study area by dividing by the lowest value, thus giving the lower value of 1. Circles – P. vivax Sri Lanka [9]; Triangles – P. falciparum Sri Lanka [9]; Squares – P. falciparum Kenya [10]; Crosses – P. falciparum Cameroon [11].

Table 1 Age-specific gametocyte densities with measures of transmission success in rural Cameroon [11] and published data sets [9,10].

Immune-generated age-structure is apparent in rural Cameroon, where asexual parasite densities (and variation in density) decreased significantly with age from a mean of 3965 parasites/μ l (SE ± 690) in the youngest age group to 829 parasites/μ l (SE ± 168) in the oldest age group (F3,410 = 11.51 P < 0.001 N = 413)(Table 2). Likewise, the prevalence of infection and the proportion of infections bearing gametocytes also decreased significantly in the eldest age group [11]. By contrast in the urban setting there was no effect of age group on asexual parasite densities (F4,90 = 0.68 P = 0.68, N = 95), with a high overall mean of 7626 parasites/μ l (SE ± 1013) (Table 2), as would be expected if acquisition of immunity were low. This lack of age-structure is mirrored in mosquito infection rates where there was no effect of age (F4,80 = 1.3, P = 0.27, N = 86), nor was there any relationship between gametocyte variability and transmission in the urban sample (Table 3) (F1,3 = 0.12 P = 0.752). However, there was a strong effect of gametocyte density per se (F1,84 = 23.1 P < 0.001 N = 86), which accounted for 28.1% of variation in mosquito infection rates.

Table 2 Age-specific asexual Plasmodium spp. prevalence rates and densities.
Table 3 Age-specific gametocyte densities with measures of transmission success in urban Cameroon [15].

Implications of the hypothesis

These analyses reveal a significant relationship between precision in gametocyte density and the transmission success of age groups in endemic settings of malaria transmission. Importantly, the actual age group most infectious to mosquitoes varied according to the transmission intensity, becoming younger with increasing intensity. Mean gametocyte densities by age group were not, however, related to transmission success except in urban Cameroon. As expected in endemic conditions, an age-dependent degree of acquired immunity was evident from the gradual decrease in mean asexual parasite densities. By contrast, the relationship between gametocyte variability and transmission success did not occur in an urban setting where individuals of all ages did not differ in their overall degree of immunity, as apparent from their similar asexual parasite densities.

Studies to date have examined human infectiousness to mosquitoes at three distinct levels: (1) Infectiousness of the individual according to gametocyte density [e.g. [19]], (2) the age group reservoir of infection within a site [e.g. [8, 11, 20]], and most recently (3) across populations of differing intensity [21]. Here, a gametocyte phenotype (low variance of gametocyte density) measured at the level of age group, but which incorporates individual infection data, is consistently related to transmission success across populations from differing transmission intensities. This relationship is consistent with the thesis that natural selection should mould phenotypes (gametocyte density) such that they are less variant in the most productive habitat. Why certain age groups are the most productive (i.e. infectious) is less clear. Gametocyte epidemiology is complex but does display some age-structure (For an excellent review see [22]). In semi-immune populations, gametocyte prevalence does not always decrease with age [20] and, when it does, tends to reflect age-asexual parasite prevalence profiles [8, 10, 23]. This suggests that low gametocyte prevalence reflects the development of immunity against the asexual stages [24] rather than any effective anti-gametocyte immune response. However, a recent reappraisal suggests that gametocyte prevalence may decrease more rapidly with age than asexual parasite prevalence rates and that there might be additional anti-gametocyte immune responses [22]. Moreover, when gametocytes are present in older age groups, their densities, relative to the asexual parasite density whence they arise, are generally increased [22]. This is consistent with the known influence of both specific and non-specific anti-asexual parasite immune mechanisms on the rate of conversion from asexual to gametocyte stage parasites [2527].

How such variation in gametocyte prevalence rates and densities per human age group relate to transmission success is unclear. Although increasing gametocyte density tends to result in greater infectiousness to mosquitoes [19], it has been repeatedly demonstrated that high gametocyte densities do not guarantee high mosquito infection rates [15, 19, 20]. Recent work within the context of developing a transmission-blocking vaccine has highlighted concurrent gametocyte density, but not age, as being positively correlated with transmission-blocking immunity [28]. Not unsurprisingly, there exists a minimum gametocyte density necessary for effective transmission, but one which is variable according to study site [29]. Low asexual parasite densities, such as those found in older semi-immune individuals, might therefore generate insufficient gametocytes for transmission. In the rural Cameroon data set here, however, mean and even mode gametocyte densities did not differ among age groups. Conversely, high asexual parasite densities, such as those seen in infants, generate an excess of gametocytes [30] and thus may induce transmission blocking anti-gametocyte specific immunity. Moreover, concurrent fever, induced by high asexual parasitaemia, reduces infectivity [31, 32]. An age-dependent pyrogenic threshold of P. falciparum parasitaemia has been demonstrated [33]; notably the individuals most tolerant of parasite density were neither the very young nor the older age groups. Such a threshold is expected to vary according to the transmission intensity. This might be an explanation for our observed differences in age of the best-transmitting groups in areas of differing transmission intensity.

Although age-specific transmission success could be explained with reference to age-specific pyrogenic thresholds, it is not clear why variation in gametocyte density should correlate similarly with age. Variation in asexual parasitaemia was found to decrease with age as might be expected with the development of immunity and yet variation in gametocyte density did not. Although variation in gametocyte density in the youngest age groups may thus be explained as a consequence of the variation in asexual parasite densities, the same argument can not simultaneously explain the increased variability in the older age groups. Specific anti-asexual parasite immunity may additionally be involved in older age groups acting indirectly in conflicting ways on gametocyte density: both reducing the asexual parasite source and increasing the rate of conversion to gametocytes. This dual effect on gametocyte density could generate increased variability in both gametocyte density and transmission success. Potentially there exists a hitherto unidentified anti-gametocyte specific response that increases with the acquisition of anti-parasite (asexual) immunity.

Hitherto unmentioned, is the documented relationship between gametocyte production and haematological factors [3436]. Although oxidative stress by cytokine response to infection can lead to alteration of uninfected RBC membranes and thus enhance their removal by the spleen [37], to date no age-specific anti-parasite immune response affecting haemolysis and hence erythropoiesis has been described [38]. Indeed, in regions of high intensity, anaemia occurs in infants and in certain risk groups (pregnant women); older age groups do not apparently suffer from malaria-induced anaemia though haematopoietic activity has not been surveyed as such across a population.

Malaria parasite gametocyte allocation has proven highly informative, revealing adaptive strategies that optimize transmission according to parasite population structure [39, 40] and to ensure transmission success [4143]. Here, examination of variation in gametocyte density reveals a measure that captures the essence of the parasite's behaviour with respect to transmission. The parasite has a very precise gametocyte allocation in individuals of an age from whom transmission to mosquitoes is best. This suggests that the parasite is best adapted to specific age groups, i.e. the parasite is being more precise in the better habitat.

The human immune response characterizes habitat quality of humans for malaria parasites. Infection longevity would be expected to be beneficial for the parasite and thus those individuals with the longest duration of infection be predicted to be the highest quality habitats. The duration of an infection has proved, however, a very difficult parameter to measure and recently the subject of renewed discussion [2, 3, 7]. Nevertheless, estimates made in areas of high transmission intensity suggest that the longest duration occurs in intermediate age groups [5, 6]; it is not known, however, how this relates to transmission success or what are the principle immunological mechanisms involved. Immuno-tolerance is an accepted feature of malaria parasite sustained infection and a role for the innate immune response has been recently proposed [44]. The age groups most infectious to mosquitoes are those who seemingly have the lowest recovery rate from infection [5, 6], when extrapolating from studies in regions of similar transmission intensity to the Kenyan and Cameroon study examined here [10, 11]. The longer duration of infection in individuals apparently tolerant to increased parasite densities may provide sufficient time, in a relatively stable environment, with which to produce gametocytes. This would provide the link between immune-dependent habitat quality and transmission success.

Malaria parasite epidemiology may be best understood with reference to reproductive allocation. A deeper understanding of the immunological basis to the host-parasite interaction and what defines a "better host" will be more readily achieved by inspection of host anti-parasite immune responses in age groups defined by the reproductive precision of their infecting parasites. Importantly, the best human habitats were those who had had previous exposure to the parasite. Vaccine strategies designed to reduce morbidity in infants may be at risk of generating of better quality individuals for parasite transmission. In the light of the current widespread vaccine trials this must be carefully monitored to avoid undesirable side-effects.


  1. 1.

    Anderson RM, May RM: Infectious diseases of humans: Dynamics and Control. 1991, Oxford: Oxford University Press

    Google Scholar 

  2. 2.

    Smith DL, Dushoff J, Snow RW, Hay SI: The entomological inoculation rate and Plasmodium falciparum in African children. Nature. 2005, 438: 492-495. 10.1038/nature04024.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  3. 3.

    Filion GJ, Paul RE, Robert V: Transmission and immunity: the importance of heterogeneity in the fight against malaria. Trends Parasitol. 2006, 22: 345-348. 10.1016/

    Article  PubMed  Google Scholar 

  4. 4.

    Bruce MC, Donnelly CA, Packer M, Lagog M, Gibson N, Narara A, Walliker D, Alpers MP, Day KP: Age- and species-specific duration of infection in asymptomatic malaria infections in Papua New Guinea. Parasitology. 2000, 121: 247-256. 10.1017/S0031182099006344.

    Article  PubMed  Google Scholar 

  5. 5.

    Bekessy A, Molineaux L, Storey J: Estimation of incidence and recovery rates of Plasmodium falciparum parasitaemia from longitudinal data. Bull World Health Organ. 1976, 54: 685-691.

    PubMed Central  CAS  PubMed  Google Scholar 

  6. 6.

    Smith T, Beck HP, Kitua A, Mwankusye S, Felger I, Fraser-Hurt N, Irion A, Alonso P, Teuscher T, Tanner M: Age dependence of the multiplicity of Plasmodium falciparum infections and of other malariological indices in an area of high endemicity. Trans R Soc Trop Med Hyg. 1999, 93 (Suppl 1): 15-20. 10.1016/S0035-9203(99)90322-X.

    Article  PubMed  Google Scholar 

  7. 7.

    Sama S, Owusu-Agyei S, Felger I, Dietz K, Smith T: Age and seasonal variation in the transition rates and detectability of Plasmodium falciparum malaria. Parasitology. 2006, 132: 13-21. 10.1017/S0031182005008607.

    Article  CAS  PubMed  Google Scholar 

  8. 8.

    Muirhead-Thomson RC: The malarial infectivity of an African village population to mosquitoes (Anopheles gambiae). Am J Trop Med Hyg. 1957, 6: 971-979.

    CAS  PubMed  Google Scholar 

  9. 9.

    Gamage-Mendis AC, Rajakaruna J, Carter R, Mendis KN: Infectious reservoir of Plasmodium vivax and Plasmodium falciparum malaria in an endemic region of Sri Lanka. Am J Trop Med Hyg. 1991, 45: 479-487.

    CAS  PubMed  Google Scholar 

  10. 10.

    Githeko AK, Brandling-Bennett AD, Beier M, Atieli F, Owaga M, Collins FH: The reservoir of Plasmodium falciparum malaria in a holoendemic area of western Kenya. Trans R Soc Trop Med Hyg. 1992, 86: 355-358. 10.1016/0035-9203(92)90216-Y.

    Article  CAS  PubMed  Google Scholar 

  11. 11.

    Bonnet S, Paul REL, Gouagna C, Safeukui I, Meunier J-Y, Gounoue R, Boudin C: Level and dynamics of malaria transmission and morbidity in an equatorial area of South Cameroon: entomological and parasitological data in a population exposed to low and seasonal transmission. Trop Med Int Health. 2002, 7: 1-8. 10.1046/j.1365-3156.2002.00861.x.

    Article  Google Scholar 

  12. 12.

    Kawecki TJ, Stearns SC: The evolution of life histories in spatially heterogeneous environments: optimal reaction norms revisited. Evol Ecol. 1993, 7: 155-174. 10.1007/BF01239386.

    Article  Google Scholar 

  13. 13.

    Mendis C, Gamage-Mendis AC, De Zoysa AP, Abhayawardena TA, Carter R, Herath PR, Mendis KN: Characteristics of malaria transmission in Kataragama, Sri Lanka: a focus for immuno-epidemiological studies. Am J Trop Med Hyg. 1990, 42: 298-308.

    CAS  PubMed  Google Scholar 

  14. 14.

    Bonnet S, Gouagna CL, Paul RE, Safeukui I, Meunier J-Y, Boudin C: Estimation of malaria transmission from humans to mosquitoes in two neighbouring villages in south Cameroon: evaluation and comparison of several indices. Trans R Soc Trop Med Hyg. 2003, 97: 53-59. 10.1016/S0035-9203(03)90022-8.

    Article  CAS  PubMed  Google Scholar 

  15. 15.

    Tchuinkam T, Mulder B, Deschering K, Stoffels H, Verhave JP, Cot M, Carnevale P, Meuwissen JH, Robert V: Experimental infections of Anopheles gambiae with Plasmodium falciparum of naturally infected gametocyte carriers in Cameroon: factors influencing the infectivity to mosquitoes. Trop Med Parasitol. 1993, 44: 271-276.

    CAS  PubMed  Google Scholar 

  16. 16.

    Bonnet S, Gouagna C, Safeukui I, Meunier J-Y, Boudin C: Comparison of artificial membrane feeding with direct skin feeding to estimate infectiousness of Plasmodium falciparum gametocyte carriers to mosquitoes. Trans R Soc Trop Med Hyg. 2000, 94: 103-106. 10.1016/S0035-9203(00)90456-5.

    Article  CAS  PubMed  Google Scholar 

  17. 17.

    Awono-Ambene HP, Diawara L, Robert V: Comparison of direct and membrane feeding methods to infect Anopheles arabiensis with Plasmodium falciparum. Am J Trop Med Hyg. 2001, 64: 32-34.

    CAS  PubMed  Google Scholar 

  18. 18.

    Krebs CJ: Ecological Methodology. 1999, Harlow: Addison Wesley Longman, 2

    Google Scholar 

  19. 19.

    Jeffery GM, Eyles DE: Infectivity to mosquitoes of Plasmodium falciparum as related to gametocyte density and duration of infection. Am J Trop Med Hyg. 1955, 4: 781-789.

    CAS  PubMed  Google Scholar 

  20. 20.

    Boudin C, Olivier M, Molez JF, Chiron JP, Ambroise-Thomas P: High human malarial infectivity to laboratory-bred Anopheles gambiae in a village in Burkina Faso. Am J Trop Med Hyg. 1993, 48: 700-706.

    CAS  PubMed  Google Scholar 

  21. 21.

    Killeen GF, Ross A, Smith T: Infectiousness of malaria-endemic human populations to vectors. Am J Trop Med Hyg. 2006, 75: 38-45.

    PubMed  Google Scholar 

  22. 22.

    Drakeley C, Sutherland C, Bousema JT, Sauerwein RW, Targett GA: The epidemiology of Plasmodium falciparum gametocytes: weapons of mass dispersion. Trends Parasitol. 2006, 22: 424-30. 10.1016/

    Article  PubMed  Google Scholar 

  23. 23.

    Jones TR, McElroy PD, Oster CN, Beier JC, Oloo AJ, Onyango FK, Chumo DK, Sherwood JA, Hoffman SL: Plasmodium falciparum gametocytaemia in Kenyan children: associations among age, intensity of exposure to transmission, and prevalence and density of subsequent gametocytaemia. Am J Trop Med Hyg. 1997, 56: 133-136.

    CAS  PubMed  Google Scholar 

  24. 24.

    Carter R, Mendis K: Transmission immunity in malaria: reflections on the underlying immune mechanisms during natural infections and following artificial immunization. Mem Inst Oswaldo Cruz Rio. 1992, 87 (Suppl III): 169-173.

    Article  PubMed  Google Scholar 

  25. 25.

    Smalley ME, Brown J: Plasmodium falciparum gametocytogenesis stimulated by lymphocytes and serum from infected Gambian children. Trans R Soc Trop Med Hyg. 1981, 75: 316-317. 10.1016/0035-9203(81)90348-5.

    Article  CAS  PubMed  Google Scholar 

  26. 26.

    Kwiatkowski D, Sambou I, Twumasi P, Greenwood BM, Hill AVS, Manogue KR, Cerami A, Castracane J, Brewster DR: TNF concentration in fatal cerebral, non-fatal cerebral, and uncomplicated Plasmodium falciparum malaria. Lancet. 1990, 336: 1201-1204. 10.1016/0140-6736(90)92827-5.

    Article  CAS  PubMed  Google Scholar 

  27. 27.

    Dyer M, Day KP: Commitment to gametocytogenesis in Plasmodium falciparum. Parasitol Today. 2000, 16: 102-107. 10.1016/S0169-4758(99)01608-7.

    Article  CAS  PubMed  Google Scholar 

  28. 28.

    Boudin C, Van Der Kolk M, Tchuinkam T, Gouagna C, Bonnet S, Safeukui I, Mulder B, Meunier JY, Verhave JP: Plasmodium falciparum transmission blocking immunity under conditions of low and high endemicity in Cameroon. Parasite Immunol. 2004, 26: 105-110. 10.1111/j.0141-9838.2004.00689.x.

    Article  CAS  PubMed  Google Scholar 

  29. 29.

    Paul REL, Bonnet S, Boudin C, Tchuinkam T, Robert V: Aggregation in malaria parasites places limits on mosquito infection rates. Infect Genet Evol.

  30. 30.

    Davidson G, Draper CC: Field studies of some of the basic factors concerned in the transmission of malaria. Trans R Soc Trop Med Hyg. 1953, 47: 522-535. 10.1016/S0035-9203(53)80005-2.

    Article  CAS  PubMed  Google Scholar 

  31. 31.

    Naotunne TS, Karunaweera ND, Del Giudice G, Kularatne MU, Grau GE, Carter R, Mendis KN: Cytokines kill malaria parasites during infection crisis: extracellular complementary factors are essential. J Exp Med. 1991, 173: 523-529. 10.1084/jem.173.3.523.

    Article  CAS  PubMed  Google Scholar 

  32. 32.

    Gouagna LC, Ferguson HM, Okech BA, Killeen GF, Kabiru EW, Beier JC, Githure JI, Yan G: Plasmodium falciparum malaria disease manifestations in humans and transmission to Anopheles gambiae: a field study in Western Kenya. Parasitology. 2004, 128: 235-243. 10.1017/S003118200300444X.

    Article  CAS  PubMed  Google Scholar 

  33. 33.

    Rogier C, Commenges D, Trape JF: Evidence for an age-dependent pyrogenic threshold of Plasmodium falciparum parasitemia in highly endemic populations. Am J Trop Med Hyg. 1996, 54: 613-619.

    CAS  PubMed  Google Scholar 

  34. 34.

    Trager W, Gill GS, Lawrence C, Nagel RL: Plasmodium falciparum: enhanced gametocyte formation in vitro in reticulocyte-rich blood. Exp Parasitol. 1999, 91: 115-118. 10.1006/expr.1998.4347.

    Article  CAS  PubMed  Google Scholar 

  35. 35.

    Price R, Nosten F, Simpson JA, Luxemburger C, Phaipun L, ter Kuile F, van Vugt M, Chongsuphajaisiddhi T, White NJ: Risk factors for gametocyte carriage in uncomplicated falciparum malaria. Am J Trop Med Hyg. 1999, 60: 1019-1023.

    CAS  PubMed  Google Scholar 

  36. 36.

    Nacher M, Singhasivanon P, Silachamroon U, Treeprasertsuk S, Tosukhowong T, Vannaphan S, Gay F, Mazier D, Looareesuwan S: 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 

  37. 37.

    Ho M, White NJ, Looareesuwan S, Wattanagoon Y, Lee SH, Walport MJ, Bunnag D, Harinasuta T: Splenic Fc receptor function in host defense and anemia in acute Plasmodium falciparum malaria. J Infect Dis. 1990, 161: 555-561.

    Article  CAS  PubMed  Google Scholar 

  38. 38.

    Paul REL, Brey PT: Malaria parasites and red blood cells: from anaemia to transmission. Mol Cells. 2003, 15: 139-149.

    CAS  PubMed  Google Scholar 

  39. 39.

    Read AF, Narara A, Nee S, Keymer AE, Day KP: Gametocyte sex ratios as indirect measures of outcrossing rates in malaria. Parasitology. 1992, 104: 387-395.

    Article  PubMed  Google Scholar 

  40. 40.

    Paul REL, Packer M, Walmsley M, Ranford-Cartwright LC, Paru R, Lagog M, Day K: Mating patterns in malaria parasite populations of Papua New Guinea. Science. 1995, 269: 1709-1711. 10.1126/science.7569897.

    Article  CAS  PubMed  Google Scholar 

  41. 41.

    Paul REL, Coulson TN, Raibaud A, Brey PT: Sex determination in malaria parasites. Science. 2000, 287: 128-131. 10.1126/science.287.5450.128.

    Article  CAS  PubMed  Google Scholar 

  42. 42.

    West SA, Smith TG, Nee S, Read AF: Fertility insurance and the sex ratios of malaria and related hemospororin blood parasites. J Parasitol. 2002, 88: 258-263.

    Article  PubMed  Google Scholar 

  43. 43.

    Reece SE, Duncan AB, West SA, Read AF: Host cell preference and variable transmission strategies in malaria parasites. Proc R Soc B-Biol Sci. 2005, 272: 511-517. 10.1098/rspb.2004.2972.

    Article  Google Scholar 

  44. 44.

    Boutlis CS, Yeo TW, Anstey NM: Malaria tolerance – for whom the cell tolls?. Trends Parasitol. 2006, 22: 371-377. 10.1016/

    PubMed Central  Article  CAS  PubMed  Google Scholar 

Download references


This work was supported by the Strategic Anopheles Horizontal Research Programme, Institut Pasteur.

Author information



Corresponding author

Correspondence to Richard E Paul.

Additional information

Authors' contributions

RELP conceived the idea, analyzed the data and wrote the paper. SB & CB performed all the data collection in rural Cameroon, made comments on and edited the paper. TT & VR performed all the data collection in Yaounde, TT made comments on and edited the paper and VR co-wrote sections of the paper.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Paul, R.E., Bonnet, S., Boudin, C. et al. Age-structured gametocyte allocation links immunity to epidemiology in malaria parasites. Malar J 6, 123 (2007).

Download citation


  • Malaria
  • Transmission Success
  • Gametocyte Density
  • Malaria Epidemiology
  • Asexual Parasitaemia