- Open Access
Plasmodium falciparum gametocyte dynamics in areas of different malaria endemicity
Malaria Journal volume 7, Article number: 249 (2008)
The aim of this study was to identify and compare factors associated with Plasmodium falciparum gametocyte carriage in three regions of differing malaria endemicity.
Retrospective data from Thailand, The Gambia and Tanzania were used. The data came from large prospective field-based clinical trials, which investigated gametocyte carriage after different anti-malarial drug treatments.
Gametocytaemia was detected during the observation period in 12% of patients (931 out of 7548) in Thailand, 34% (683 out of 2020) in The Gambia, and 31% (430 out of 1400) in Tanzania (p < 0.001). Approximately one third (33%, 680/2044) of the patients with gametocytaemia during the observation period, already had patent gametocytaemia at enrolment (day 0 or day 1): 35% (318/931) in Thailand, 37% (250/683) in The Gambia, 26% (112/430) in Tanzania. Maximum gametocytaemia was usually observed on or before the seventh day after starting treatment (93% in Thailand, 70% in Tanzania and 78% in The Gambia). Lowest gametocyte carriage rates were observed following treatment with artemisinin derivatives, while sulphadoxine-pyrimethamine (SP) was associated with significantly greater development of gametocytaemia than other drug treatments (p < 0.001). The duration of gametocyte carriage was shorter in Thailand by 86% and Tanzania by 65% than in The Gambia. Gametocyte carriage was 27% longer among people presenting with anaemia, and was shorter in duration among patients who received artemisinin derivatives, by 27% in Thailand and by 71% in Tanzania and The Gambia.
This study confirms the independent association of gametocytaemia with anaemia, and the significantly lower prevalence and duration of gametocyte carriage following treatment with an artemisinin derivative. The large differences in gametocyte carriage rates between regions with different levels of malaria transmission suggest that drug interventions to prevent transmission will have different effects in different places.
Malaria control rests traditionally on two strategies; vector control (reducing the numbers of anopheline vectors and reducing the probability of being bitten), and drug treatment. Effective anti-malarials reduce morbidity, prevent mortality, and reduce asexual parasite biomass. This, in turn, reduces the numbers of sexual (transmissible) parasites (gametocytes). In Plasmodium falciparum infections, it remains unclear whether gametocyte production is programmed early on after hepatic schizogony or is a response to stimuli acting upon the parasite population. The proportion of parasites committed to sexual stage development may change during the course of an infection. The developing sexual stages (stages I to IV) remain sequestered in the microvasculature for approximately 10 days before appearing as morphologically distinct male and female stage V gametocytes in the peripheral blood. One male (containing eight microgametes) and one female (macro-gamete) are required per mosquito blood meal (approx 2 μL) for infection to occur. Thus gametocyte densities of 1 per μL are theoretically sufficient to infect mosquitoes, a density beneath the limit of detection for most routine microscopy. This explains malaria transmission from subjects without apparent gametocytaemia.
In areas of low or unstable transmission most malaria infections are transmitted by people who are ill, or recovering from symptomatic malaria. In such areas asymptomatic infections are unusual, so treatment-seeking behaviour and the pharmacokinetic and pharmacodynamic properties of the anti-malarial drugs used are important determinants of transmission. But low transmission settings often contain small foci of higher transmission and the few asymptomatic individuals in these areas are important in sustaining malaria through the dry season. In higher transmission settings the situation is more complex. Anti-disease controlling immunity is acquired which results in an increasing proportion of infections stabilizing at relatively high parasite densities. These may be either asymptomatic or tolerated in older patients who are less likely to seek treatment. These infections may still be transmissible . In addition to immunity against asexual stages there is the development of a specific immunity against sexual stage parasites which can further reduce the transmission probability per infection [2, 3]. Several investigators have assessed the contributions of different age and patient groups to overall malaria transmission in endemic areas by direct measurement of infectivity to vector mosquitoes. These studies, which have been conducted in different geographic regions with differing patterns of malaria epidemiology, emphasize the importance of older often asymptomatic individuals in sustaining transmission of both falciparum and vivax malaria despite usually lower parasite densities [4–10].
Recrudescence (treatment failure), and the subsequent extended duration of infection is an important source of transmission particularly of drug-resistant parasite genotypes . The factors associated with gametocytogenesis include the parasite density itself [12, 13], anaemia , the duration of infection, stresses on the parasite population such as host immunity (related to age) or anti-malarial treatment, and the stage specificity of the anti-malarial drugs used. The relationship between gametocyte density in blood and the transmission probability is generally sigmoid [15–17], although it varies between individuals and is affected by the factors described above [18, 19].
The aim of this study was to identify factors associated with gametocyte carriage at different levels of malaria transmission and thus clinical epidemiology. Retrospective data from study sites in three different endemic areas were used. The data come from large prospective field-based clinical trials which investigated gametocyte carriage after different anti-malarial drug treatments.
Description of study sites
The Tanzanian studies were carried out near Ifakara, in the Kilombero district in the southeast of the country between 1997 and 1999. At the time falciparum malaria was the foremost health problem in the district. In a one-year period, 37% of children (aged 0–15) admitted to hospital in Ifakara were diagnosed as having malaria with a case fatality rate of 3% . Most of these were infants. At the time these data were collected, the parasite prevalence in rural areas was approximately 70%, rising to 90% in children under five years, and remained high all year-round. Nearly all these cases were P. falciparum, although Plasmodium malariae and Plasmodium ovale infections also occurred occasionally. Entomological inoculation rates (EIRs) were between 200 and 300 infective bites per person per year, with no seasonal variation. In this region, there were a large number of government dispensaries and chloroquine was widely available in local shops. This study on gametocyte dynamics was part of a larger study investigating natural transmission-blocking immunity in this area of intense perennial transmission .
The Gambian study catchment area was roughly 275 km2, located at two sites, one to the east and the other to the west of Farafenni, a rural town 170 km to the east of the Atlantic coast. During cross-sectional malaria surveys, between 7% to 9% of children under five years of age were found to be parasitaemic. The overall the malaria attributable mortality in this region in the 1980s was estimated to be 6.3 per 1,000 in children under the age of one year and 10.7 per 1000 in children aged 1–4 years . Transmission was highly seasonal following the rains in July to September with an EIR of approximately 4  with a peak of clinical attacks at the end of the wet season. Symptomatic malaria was rare between February and July. Data used in our study comes from three studies [16, 24, 25].
The studies in Thailand were carried out in camps for displaced persons of the Karen ethnic minority on the western border of Thailand. Transmission of malaria was unstable and seasonal, occurring through out the year with peaks in May-July and December-January . The estimated EIR and corresponding incidence of malaria was low (approximately 0.5 to 1.5 cases/person/year) with prevalence rates of 1–4% for P. falciparum . Overall, P. falciparum accounted for 37% of malaria infections, with the remainder due to Plasmodium vivax. All P. falciparum infections and approximately 90% of P. vivax infections were symptomatic, but the case fatality rate was low (~0.1%). Severe disease was more common in young children and pregnant women . Between 1989 and 2000, more than 10,000 people were enrolled in more than 25 drug studies. The data come from an amalgamation of 18 studies (mostly comparative trials, but some simply monitoring therapy) conducted between 1990 and 1996 [28, 29]. Data from the artemether-lumefantrine studies were gathered between 1996 and 1999 [30, 31].
Design of studies
In all studies symptomatic patients attending outpatient clinics were considered for enrolment. In Thailand, patients of all ages were enrolled providing that they weighed more than five kilograms, whereas in Tanzania, all patients were older than one year (range 1–67 years). In The Gambia all patients were under ten years of age, weighed more than five kilograms and had a parasitaemia greater than 500 per microlitre in the presence of a febrile illness. Fully informed consent was obtained from all patients or their parents or guardians. Pregnant women and patients with severe disease were excluded at all three sites.
Patients were monitored following drug treatment in Thailand weekly until day 28 in early studies, and until day 63 in later studies. In Tanzania post-treatment follow up was weekly until day 28, while in The Gambia follow up was on day 7, day 14 and day 28.
Parasite counts were calculated from Giemsa-stained thick and thin blood films. In Thailand, parasitaemia was expressed as the number of parasitized erythrocytes per 1,000 red cells or the number of parasites seen on the thick film per 500 white blood cells; in The Gambia and in Tanzania it was expressed per 200 white blood cells. Gametocyte density was determined on the basis on the number of gametocytes per 500 white blood cells in Thailand and Tanzania, and per 1,000 white cells in The Gambia. In Tanzania slides were declared negative if no parasite was seen in 100 high power fields.
Description of the data
In all studies patient characteristics were recorded at enrolment including: age, sex, haematocrit, temperature, parasitaemia, length of the illness before enrolment, size of liver, size of spleen and the patient's symptoms. The size of the liver, spleen, the haematocrit and parasitaemia were also measured at the follow-up examinations. In majority of patients in Tanzania parasitaemia was only recorded as a categorical variable with five levels: < 800/μL, 800–1,599/μL, 1,600–3,999/μL, 4,000–19,999/μL and ≥ 20,000/μL. Patients with no recorded parasitaemia on enrolment were excluded from the analyses.
There was a considerable variability in malaria treatments between the studies. In Thailand chloroquine and sulphadoxine-pyrimethamine were not evaluated because of the high level resistance. Artemisinin derivatives either alone or in combination were the only drugs evaluated at all three locations. Eight treatment groups were defined (Table 1).
Four principal outcomes were evaluated:
presence of gametocytes on enrolment (binary);
appearance of gametocytes after enrolment up to day 14 (binary);
length of gametocyte carriage within 28 days follow-up (continuous);
level of maximum gametocytaemia (continuous).
Presence of gametocytes on enrolment was defined as a gametocytaemia detected by microscopy on day 0. Only patients with a recorded count available on day 0 were included in the analysis of this endpoint. Appearance of gametocytes after enrolment up to day 14 was defined in patients with no detectable gametocytaemia on enrolment and, for the purpose of comparison between sites, was assessed only from patients with blood slide assessments on days 7 and 14. For the analysis of factors associated with gametocytaemia within each site, in addition to mandatory counts on days 7 and 14, counts between days 1 and 6, if available, were used to evaluate carriage. Thus for this analysis, if there was a positive count at anytime in the first two weeks, the patient was classified as gametocytaemic regardless of availability of counts on day 7 or 14.
The duration of gametocyte carriage was defined as the time from the first recording of gametocytaemia until half way between the last positive record and the subsequent negative record. For patients with missing records after a positive record of gametocytaemia, half way between the last positive record and the next scheduled assessment (7, 14, 21 or 28 days) was used to calculate duration of carriage, and these observations were treated as censored. For patients with intermittent gametocytaemia, duration of carriage was defined as time from the first positive count to the last positive count recorded. Missing assessments between positive counts were treated as positive. Maximum gametocytaemia was defined as the maximum density of gametocytaemia recorded over the whole follow-up period.
Infectivity was not assessed in these studies. Whilst accepting that the relationship between infectivity and parasite densities is dependent on many factors and varies considerably between individuals, there is nevertheless a quantitative relationship between the two. To illustrate this, gametocytaemia-infectivity relationships from two data sets were used to construct putative infectivity. For each patient gametocyte densities on days 0, 7, 14, 28 were used to interpolate (using the straight line interpolation between two subsequent points) levels for all days between day 0 and day 28. These levels were then converted into the probabilities pi of infecting a feeding anopheline mosquito on the i-th day using a sigmoid relationship (Figure 1) obtained from fitting a non-symmetric Gompertz model to the estimated proportion of infected mosquitoes:
Method A: 1.08 exp(- exp(-0.86*(log10(gametocyte count) * – 1.48)))
This is a data set derived from non-immune volunteers .
Method B: 0.07 + 0.38 exp(-exp(-4.09 *(log10(gametocyte count) – 2.09)))
This is a data set derived from semi-immune children in The Gambia .
The model fits were satisfactory (both R2 ≥ 0.99) with the estimated values between the observed values close to linear interpolation. The same relationship was used for all treatments.
As a measure reflecting overall infectivity of an individual, the area under the curve (AUIC) of the probability of infecting a biting female anopheline mosquito and time was estimated  and expressed in infectivity days. This was calculated using trapezoid rule as any short segment of the infectivity curve could be approximated by a straight line. This measure reflects the potential gametocyte exposure to feeding anopheline mosquitoes, but does not take into account any individual drug effect on infectivity which is independent of gametocyte density (such as sporontocidal activity). As a measure of overall infectivity of a population the individual AUICs from all patients (carrying gametocytes or not) were added up and presented per 1,000 patients.
Failure rates were estimated by survival analysis using Kaplan-Meier method. Recurrence of parasitaemia was treated as failure and patients lost to follow-up were censored at the last follow-up visit.
The association between the area under the gametocyte density time curve (AUC) and gametocyte count on day 7 was examined using Spearman's correlation. AUC was calculated using trapezoid method over 28 day follow-up period and was restricted to patients who had measurements taken on day 0,7,14 and 28. To calculate AUC and AUIC, it was assumed that missing counts which follow gametocytaemia clearance (i.e. at least one negative count after positive counts) remained negative. To avoid confounding with recrudescence or reinfection in assessing the asexual to sexual stage transition, the ratio between maximum gametocytaemia and maximum asexual parasitaemia was calculated only in patients with maximum gametocytaemia recorded within 14 days after starting treatment. Parasitaemia on enrolment was used for the maximum parasitaemia.
For the purpose of analyses the following variables were dichotomized: temperature (oral > = 38C or axillary > 37.5C), history of illness (> = 3 days), liver size (> 0 cm), spleen size (> 0 cm). Anaemia and severe anaemia were defined on the basis of haematocrit using cutoffs of 30% and 20% respectively. Parasitaemia on enrolment was examined as a categorical variable. The effect of continuous parasitaemia (after log 10 transformation) was also examined in studies with available data.
The duration of gametocyte carriage was analysed by survival methods. Kaplan-Meier estimates (K-M) of gametocyte carriage rates were calculated at day 7 and 14 and compared between groups using the log-rank test stratified by study site. Heterogeneity of results between studies within each country was examined through the test of interaction in the Cox proportional hazards model. As the proportional hazards assumption was not satisfied between sites, effects of covariates in the combined dataset and their interaction with site were examined in a lognormal model with accelerated failure-time parameterization and with gamma frailty (to account for study and site effects) .
Associations between maximum gametocytaemia and explanatory variables were examined using a negative binomial regression model [33, 34] with random effects to accommodate study variation in dispersion. Results were summarized in groups of patients by ratios of mean densities (exponent of the regression coefficients). Only patients who carried gametocytes at some time during the follow-up contributed to this analysis.
Relationships between binary outcomes and explanatory variables were examined using the Mantel-Haenzel method, stratified by study site. The effects of covariates and their interaction with country in the combined data set were explored using logistic regression models with random effects. In all analyses, unless stated otherwise, age, haematocrit and enrolment parasitaemia were treated as continuous covariates. Enrolment parasitaemia was log-transformed and haematocrit was adjusted for each site population median.
In multivariate analyses, all covariates examined in corresponding univariate analyses were investigated. Final models were selected by stepwise forward variable selection procedure and only covariates significant at 5% were included in the final model. For continuous variables categorical and continuous representations were investigated and the one which maximized the log-likelihood function was selected. If the difference in the log-likelihood was not significant, the continuous representation was selected.
All analyses were performed on the maximum possible data set, i.e. on all patients with non-missing required data and since different analysis used different endpoints and covariates different subsets of patients were analysed throughout the paper. Numbers of patients used in each analysis are given.
In total, there were 10,968 patients with confirmed falciparum malaria on enrolment and at least one gametocyte assessment available on enrolment or during the follow-up. Treatments administered and their efficacies are summarized in Table 1.
The demographic characteristics of patients in the three populations studied are shown in Table 2. There was a significant difference in the age distribution between sites. After adjustment for age the patients in Thailand were less anaemic (p < 0.001) and had lower parasite counts (p < 0.001) than in the two African study sites. Figure 2 shows the distributions of enrolment parasitaemia, age and haematocrit for the three sites.
Overall, 19% (2044/10,968) of patients had patent gametocytaemia during the observation period: 12% (931/7548) in Thailand, 34% (683/2020) in The Gambia and 31% (430/11400) in Tanzania (p < 0.001, chi-square test). Approximately one third (33%, 680/2044) of the patients with gametocytaemia at some stage of the infection, already had patent gametocytaemia at enrolment (day 0 or day 1): 35% (318/931) in Thailand, 37% (250/683) in The Gambia, 26% (112/430) in Tanzania. Maximum gametocyte densities were different between the sites, being lowest in Tanzania (median = 55/μL), highest in The Gambia (median = 256/μL) and with a median maximum gametocytaemia of 180/μL in Thailand (p < 0.001, Kruskal-Wallis test). In patients presenting with gametocytes the maximum recorded gametocytaemia occurred on or before day 7 in 99% of Thai patients, 90% of those in Tanzania and 95% in patients from The Gambia. In patients who developed gametocytaemia after enrolment, the maximum gametocytaemia was usually recorded on day 7 but in 30% patients in Gambia, 40% in Tanzania, and 11% in Thailand it occurred after day 7 of follow-up. The area under the gametocyte count-time curve was correlated with gametocytaemia level on day 7; Spearman correlation coefficient of 0.86 (The Gambia; n = 416), 0.57 (Tanzania, n = 198) and 0.42 (Thailand; n = 729) (all p < 0.001). The median (IQR) number of measurements per patients was 4 (4–4), 5(4–5) and 8 (7–11), respectively.
The ratio between maximum gametocytaemia and maximum asexual parasitaemia (a measure of the asexual to sexual stage transition rate) was the lowest (median, 90% range) in Tanzania 0.0008 (0.0001–0.04) n = 16, compared with The Gambia 0.01 (0.0002–1.10) n = 614, and Thailand 0.04 (0.0005–1.77) n = 725 (p < 0.001, Kruskal-Wallis test). The ratio was significantly lower (p < 0.001) in the youngest age group in The Gambia and there was no difference between age groups in Thailand. The ratio was also significantly lower in patients who developed gametocytes after enrolment (median [90% range] = 0.01 [0.0002–0.80]) as compared to those with gametocytes at enrolment (0.05 [0.001–2.67]); p < 0.001. In patients who had any parasitaemia assessments between day 0 and day 7, 79% had maximum parasitaemia recorded on day 0. In the others, the ratio between maximum parasitaemia and enrolment parasitaemia was median (90% range) = 2.6 (1–50).
Gametocytaemia on enrolment
Overall, 12% (250/2018) of The Gambian patients had gametocytaemia recorded on enrolment compared to 8% (112/1400) in Tanzania and 4% (318/7502) in Thailand (Table 3); p < 0.001.
Risk factors for admission gametocyte carriage included: lower parasitaemia and anaemia in all sites; younger age in Thailand and Tanzania; lower body temperature in Thailand and The Gambia; and longer history of illness, palpable liver, palpable spleen, and pure P. falciparum infection in Thailand (Table 3).
A multivariate analysis could be conducted in three sites separately (Table 4) and for combined data from Thailand and The Gambia only because of different characteristics being collected in Tanzania. For The Gambia, Thailand and their combined data, only enrolment parasitaemia, haematocrit and temperature were selected in the final model. In the combined model, after adjusting for these covariates, gametocyte carriage rates were still much lower in Thailand (OR [95%CI] = 0.048 [0.029–0.075], p < 0.001) then in The Gambia. Although duration of illness and presence of mixed infection was a risk factor in Thailand they could not be tested in the overall model as they were not recorded in The Gambia. The odds ratios confirm the univariate results, and for parasitaemia (OR [95%CI] = 0.580 [0.492–0.683], p < 0.001) and temperature (OR [95%CI] = 0.593 [0.452–0.779], p < 0.001) were similar for the two sites, whereas there were significantly different values for enrolment haematocrit (OR [95%CI] = 0.902 [0.876–0.930] in The Gambia and 0.834 [0.805–0.864] in Thailand, p < 0.001, likelihood ratio test). In Tanzania only age and parasitaemia were selected in the final model.
Gametocytaemia after enrolment
In those patients who presented without gametocytes on enrolment, only 3% (122/4,083) overall developed gametocytes carriage during follow up after treatment with artemisinin derivatives. This proportion varies across sites. It was highest in The Gambia; 15% (96/645) compared with 0% (0/100) and 1% (26/3338) in Tanzania and Thailand. However, it should be noted that 10 patients in Tanzania and 109 in Thailand had gametocytaemia recorded between days 1 and 6 but as their readings were negative on day 7 and 14, they were included in these comparisons as not gametocytaemic. Sulphadoxine-pyrimethamine (SP) was associated with the greatest carriage rates (64% 221/348) and this effect was similar in Tanzania and The Gambia (p = 0.215, chi square test). The corresponding overall rates for developing patent gametocytaemia on day 7 or day 14 after the other treatments were 29% (114/390) following chloroquine, 3% (10/349) following halofantrine, 6% (53/878) after mefloquine, and 22% (34/153) after quinine.
In the univariate analyses (Table 5), statistically significant associations were also found between the development of gametocytaemia and young age, long history of illness, anaemia and high parasitaemia on enrolment. In Thailand also palpable spleen, palpable liver and pure P. falciparum infection were associated with increased risk of gametocytaemia. In the multivariate analyses (Table 6), in Thailand and The Gambia, low haematocrit on enrolment (OR [95%CI] = 0.875 [0.852–0.898] in Thailand and 0.965 [0.938–0.994] in The Gambia) was associated with an increased risk of the appearance of gametocytaemia. In addition, Thai patients with a palpable spleen (OR [95%CI] = 1.433 [1.050–1.955]) and long history of illness (OR [95%CI] = 2.909 [2.126–3.981]) and pure P.falciparum infection (OR [95%CI] = 2.762 [1.559–4.901]) and Gambian patients with high parasitaemia on enrolment (OR [95%CI] = 1.503 [1.177–1.919]) were at increased risk of developing gametocytaemia In all countries, treatment with artemisinin derivatives alone or combined with other drugs resulted in the lowest risk of gametocytaemia. In Tanzania, apart from the treatment, only young age was associated with gametocyte carriage after enrolment. Different sampling schedules prohibited us from combining data from different sites.
Length of gametocyte carriage
Of the 2027 patients with gametocytes before day 28, 12 (0.6%) had three and 164 (8%) had two separated episodes of gametocytes carriage, with a negative count between positive counts. Among a total of 2213 episodes, only 58 (3%) had one missing measurement and 7 (0.3%) episodes had 2 or more missing measurements between the first and last counts. Quality of the follow-up information varied between sites, with 68% censored observations in The Gambia, 61% in Tanzania and 19% in Thailand, although one third of these in The Gambia and one quarter in Tanzania were censored after day 14. The proportion of patients who had missing measurements before the first positive gametocyte count was 7%, 6% and 2% in Tanzania, The Gambia and Thailand.
Gametocyte carriage rates (95% CI) after 7 and 14 days were smallest in Thailand (Table 7): 34(31–37)% and 19(16–22)%. In The Gambia and Tanzania rates were very similar (p = 0.99, logrank test) and at least twice the rates in Thailand. In each site, gametocyte carriage rates after 7 days in patients who were treated with artemisinin derivatives alone or with other drugs were significantly lower than rates of carriage in patients on other treatments: 33(30 – 36)% versus 62% (58–65)% for patients treated with monotherapies (groups C-F) and 85 (71–88)% for patients treated with sulphadoxine-pyrimethamine (groups G-H) in all sites combined (p < 0.001, logrank test stratified by study). The corresponding gametocyte carriage rates (95%CI) after 14 days were 18 (15–21)%, 47 (43–51)% and 80 (75–84)%.
In Thailand other risk factors for prolonged carriage included high enrolment parasitaemia (p = 0.005), a prolonged history of illness (p = 0.038) and anaemia (p = 0.034). These associations were not apparent in the other sites.
Site, presence of gametocytaemia on enrolment, haematocrit and treatment were all independent predictors of the duration of gametocyte carriage in the multivariate analysis of the combined data from Thailand and The Gambia (Table 8). In Thailand carriage was 86% shorter than in The Gambia. In patients who received artemisinin derivatives carriage was shortened in duration by 27% in Thailand and by 71% in the Gambia. In The Gambia, gametocytaemia present before treatment had a 72% longer clearance time than gametocytaemia which appeared after treatment. In Thailand the difference in clearance time was not significant. Carriage was longer in people presenting with anaemia by 27% in both countries.
When a common model was fitted for all three data sets using common covariates (anaemia could not be included since it was missing for majority of patients in Tanzania), gametocyte carriage in Tanzania was estimated to be shorter by 65% than gametocyte carriage in The Gambia. Effects of gametocytaemia on admission and treatment with artemisinin derivatives in Tanzania were estimated to be the same as in The Gambia.
In Thailand, the maximum gametocyte density was increased in patients with gametocytaemia on enrolment (Incidence Rate Ratio (IRR) [95%CI] = 1.261 [1.081–1.430], p = 0.001, n = 925), severe anaemia (IRR [95%CI] = 1.744 [1.136–2.677], p = 0.011, n = 403) and a prolonged history of illness (> 3 days) (IRR [95%CI] = 1.234 [1.047–1.455], p = 0.019, n = 491). The two children age groups (< 5 years and 5–15 years) had similar levels of maximum gametocytaemia which were significantly higher than in adults (IRR [95%CI] = 1.135 [1.002–1.284], p = 0.046, n = 931). The effects of these covariates were the same for patients presenting with gametocytaemia on enrolment and patients who developed gametocytaemia later, as tested by the interaction term in the model. No covariates were significant in the multivariate analysis.
In The Gambia, the only determinant of maximum gametocyte density was drug treatment, but this effect was different for patients presenting with gametocytaemia on enrolment and patients who developed gametocytaemia later (p = 0.009, n = 651). In both cases, patients treated with SP produced higher number of gametocytes than patients treated with artemisinin derivatives (IRR [95%CI] = 1.880 [1.533–2.306] for patients presenting without gametocytaemia, 1.404 [1.106–1.7] for patients presenting with gametocytaemia).
In Tanzania, higher gametocyte densities were found in patients with high parasitaemia on enrolment (IRR [95%CI] = 1.461 [1.301–1.640] for increase in parasitaemia up to the next category, p < 0.001 test for trend, n = 430), females (IRR [95%CI] = 1.595 [1.244–2.045], p < 0.001, n = 391), anaemic patients (IRR [95%CI] = 2.427 [1.128–5.222], p = 0.012, n = 33), and symptomatic patients (IRR [95%CI] = 1.730 [1.186–2.524], p = 0.004, n = 362). Children aged 5–15 years (IRR [95%CI] = 0.584 [0.431–0.790], p < 0.001) and adults (IRR [95%CI] = 0.253 [0.176–0.363], p < 0.001) had significantly lower gametocyte densities compared to children less than 5 years of age (n = 427). In the multivariate model only age (IRR [95%CI] = 0.583 [0.439–0.795] for 5–15 years, 0.317 [0.207–0.424] for > 15 years, as compared to children 0–5 years) and parasitaemia (IRR [95%CI] = 2.201 [1.542–3.142] for parasitaemia increase up to the next category, p < 0.001, test for trend) remained independent determinants of gametocyte density. Anaemia was not investigated in the multivariate analysis due to the large number of missing values. Since multivariate analyses in each site resulted in models with no common covariates no analysis was performed on the combined data set.
Predictors of potential infectivity
The putative patient infectivity was estimated from the serial gametocyte counts. Among patients with gametocytes measurements available on days 0, 7, 14 and 28, all four blood film examinations were negative for gametocytes in 92% (3,643/3,970) of Thai patients, compared to 63% (177/281) of Tanzanian patients and 54% (373/685) of patients from the Gambia; p < 0.001. These patients were excluded from the analysis of infectivity. In nearly all cases (1012/1017) the area under the infectivity curve (AUIC) calculated using Method A  was median (range) 102 (4 – 505)% higher than that calculated using Method B . Among those who carried gametocytes at some stage of the infection, the area under the infectivity curve (median [90% range]) was highest in The Gambia (Method A: 11.2 [2.0–24.2] infectivity days and Method B: 4.4 [1.9–12.04] infectivity days, n = 426), markedly lower in Thailand (Method A: 6.3 [2.6–18.3] and Method B: 2.9 [1.9–11.6] infectivity days, n = 421) and lowest in Tanzania (Method A: 5.6 [2.3–14.3] and Method B: 1.9 [1.9–5.6] infectivity days, n = 170) (Figure 6). Figures 7 and 8 show distribution of AUIC by treatment group in Thailand and The Gambia. In Tanzania majority of patients (78%) received the same treatment. Assuming that Method A is more appropriate for non-immune patients in Thailand while Method B is more appropriate for patients in The Gambia and Tanzania and that 8% of symptomatic patients in Thailand, 37% in Tanzania and 46% in The Gambia will carry gametocytes at some stage of the infection; the estimated overall infectivity of a population of 1000 patients during the 28 days follow-up are 609 person infectivity days in Thailand, 901 in Tanzania and 2688 in The Gambia. Assuming that everybody in this population received one treatment, then these estimates become 486 person infectivity days for artemisinin derivatives, 530 for halofantrine and 771 for mefloquine in Thailand; 1812 person infectivity days for SP+AS and 3474 for SP alone in The Gambia; 892 person infectivity days for chloroquine in Tanzania. The proportion of patients presenting with gametocytes at enrolment was the same for the two treatments in Thailand (5%) and The Gambia (12%).
This comparison between three malariaous regions with very different epidemiological characteristics detected considerable heterogeneity in P. falciparum gametocyte carriage. Some of this heterogeneity is explained by known factors such as the differences in age (a proxy for cumulative exposure and thus immunity), levels of anaemia, and type of anti-malarial drug used . On enrolment, before drug treatment, gametocyte carriage in Thailand was 4%. This was three-fold lower than in The Gambia, and two-fold lower than in Tanzania. Within 14 days from starting treatment a further 40% of patients developed gametocytes in The Gambia and Tanzania, while only 7% more patients developed gametocytaemia in Thailand. Duration of carriage was also significantly longer in the African countries than in Thailand. In contrast, gametocytes densities were lowest in Tanzania. All this is reflected in the AUIC estimates with The Gambia having the largest individual estimates and Tanzania the smallest. Nevertheless, after adjusting for the number of gametocytes carriers, even if different relationships appropriate to the immunity status in the three regions were applied (i.e. Method A in Thailand and Method B in The Gambia and Tanzania), the Thai population of symptomatic patients was still the least infectious, with Tanzanian patients being 50% more infectious and patients in The Gambia being nearly four times as infectious.
Admission gametocytaemia was associated independently with lower haematocrit, lower parasitaemia, and lower temperature. Slowing the expansion of asexual parasitaemia may have less effect on the persistence of gametocytaemia, as gametocytes have greater longevity than the asexual stages. In Thailand, where duration of illness was recorded, gametocyte prevalence was associated with a longer period of illness before presentation. In low transmission settings the majority of gametocytaemic individuals will be symptomatic and seek treatment during the acute phase of the illness, and the majority of gametocyte carriage will occur after starting treatment. In the high transmission settings contribution of asymptomatic carriers is significant and it has been demonstrated that carriers with submicroscopic densities of gametocytes are capable of infecting mosquitoes [1, 9]. The dynamics of infection in a high transmission setting where multiple infections coexist within the human host may also be more conducive to the production of gametocytes . Whichever the explanation effective anti-malarial treatment will have a much greater impact on overall malaria transmission in areas of low transmission intensity than in high transmission areas, where a smaller fraction of parasitaemic individuals seeks treatment.
However, transmission intensity often varies greatly over small geographic distances, and low transmission settings often contain microfoci of much higher transmission intensity. For example even in Thailand, there is evidence of an untreated reservoir of infectious individuals in some areas , although these are now few. Given marked heterogeneity in gametocyte dynamics among sites, it is remarkable that several parameters were consistently associated with gametocyte carriage in all three sites; all of which are indicative of a longer period of infection, which could reflect an effective host response and the development of premunition, poor access to effective treatment in impoverished rural areas, exacerbated by drug resistance or a combination of these . Further, the substantial impact of ACT in reducing post-treatment gametocyte carriage was seen across all studies in all three countries .
In nearly all patients with gametocytaemia, densities peaked in the first week after starting treatment (day 0–7). The day 7 gametocyte density proved a good surrogate for the area under the gametocytaemia-time curve. The factors associated with post-treatment gametocytaemia were similar to those associated with enrolment gametocytaemia. There were large differences between treatment regimens; those containing an artemisinin derivative were associated with lower and shorter carriage, and treatment with SP being associated with higher carriage rates and longer duration. Although there is a relationship between the administration of anti-malarial drugs and gametocyte density, these data cannot be translated directly into transmission potential because of the different effects of the drugs on gametocyte viability [1, 16, 38–40].
As drug resistant parasites become more prevalent, the duration of malaria infections lengthens and the proportion of recrudescent infections increases. Recrudescent infections are cumulatively of longer duration than primary infections, and are associated with higher rates of anaemia and gametocyte carriage . Drug resistant parasites have been associated with higher gametocyte carriage and enhanced mosquito infectivity in the absence of appreciable clinical or parasitological treatment failure [41–43]. Therefore, in high transmission areas, where the impact of treatment on overall malaria transmission is likely to be less than in areas of low endemicity, ACT can provide substantial benefit in reducing the spread of drug resistant parasite genotypes .
This study has several limitations. Data sets from several studies were combined for the purpose of this analysis. This was a retrospective analysis so inclusion criteria for the trials were not standardized which resulted in the differences in the populations participating in the three sites. For example, in The Gambia only children were enrolled while in the other two sites a broader age spectrum was enrolled. Different covariates were assessed and laboratory methods were not standardized between the sites. Different detection limits for gametocytes affected the estimated densities and rates of carriage. It is also possible that the apparent differences in gametocyte carriage rates, which are dependent on the frequency of observation, are influenced by these differences in trial conduct. Gametocytes counts were collected at slightly different sampling schedules and there was a considerable number of missing values (20–40%) at visits scheduled for all three sites. The estimates of infectivity are illustrative, being based on two data-sets which characterized the sigmoid relationship between gametocyte densities and infectivity to an anopheline vector. This hides considerable inter-individual variability and the confounding effects of immunity and other factors. Much more data are needed to characterize and quantify the sources of variance in this assessment and produce more valid assessments necessary to model the impact of interventions on malaria transmission.
This study confirms the independent association of increased rates of gametocyte carriage with anaemia [29, 44], and with treatment with sulphadoxine-pyrimethamine [45, 46]. Gametocyte carriage is significantly reduced by combination treatment with the artemisinin derivatives (ACT), which were associated with the lowest rates of gametocyte carriage [11, 29]. The large differences in gametocyte carriage rates between regions with different levels of malaria transmission and access to treatment suggest that drug interventions to prevent transmission will have different effects in different places. In areas of high stable transmission, less reduction in transmission can be expected from effective treatment, with the major benefit of ACT deployment being reducing the spread of resistance to non-artemisinin drugs. At low levels of transmission, prompt treatment with gametocytocidal drugs can have a major effect on the transmission and thus incidence of falciparum malaria.
Artemisinin combination therapy
Area under the infectivity curve
Area under the gametocyte density time curve
Entomological inoculation rate
Gametocyte carriage rate
Incidence Rate Ratio
- P. falciparum :
- P. malariae :
- P. ovale :
- P. vivax :
Schneider P, Bousema JT, Gouagna LC, Otieno S, Vegte-Bolmer van de M, Omar SA, Sauerwein RW: Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection. Am J Trop Med Hyg. 2007, 76: 470-474.
Saeed M, Roeffen W, Alexander N, Drakeley CJ, Targett GA, Sutherland CJ: Plasmodium falciparum antigens on the surface of the gametocyte-infected erythrocyte. PLoS ONE. 2008, 3: e2280-10.1371/journal.pone.0002280.
Bousema JT, Drakeley CJ, Sauerwein RW: Sexual-stage antibody responses to P. falciparum in endemic populations. Curr Mol Med. 2006, 6: 223-229. 10.2174/156652406776055140.
Bonnet S, Gouagna LC, Paul RE, Safeukui I, Meunier JY, 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.
Drakeley CJ, Akim NI, Sauerwein RW, Greenwood BM, Targett GA: Estimates of the infectious reservoir of Plasmodium falciparum malaria in The Gambia and in Tanzania. Trans R Soc Trop Med Hyg. 2000, 94: 472-476. 10.1016/S0035-9203(00)90056-7.
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.
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.
Paul RE, Bonnet S, Boudin C, Tchuinkam T, Robert V: Age-structured gametocyte allocation links immunity to epidemiology in malaria parasites. Malar J. 2007, 6: 123-10.1186/1475-2875-6-123.
Paul RE, Diallo M, Brey PT: Mosquitoes and transmission of malaria parasites – not just vectors. Malar J. 2004, 3: 39-10.1186/1475-2875-3-39.
Paul RE, Bonnet S, Boudin C, Tchuinkam T, Robert V: Aggregation in malaria parasites places limits on mosquito infection rates. Infect Genet Evol. 2007, 7: 577-586. 10.1016/j.meegid.2007.04.004.
Price RN, Nosten F, Luxemburger C, ter Kuile FO, Paiphun L, Chongsuphajaisiddhi T, White NJ: Effects of artemisinin derivatives on malaria transmissibility. Lancet. 1996, 347: 1654-1658. 10.1016/S0140-6736(96)91488-9.
Carter R, Miller LH: Evidence for environmental modulation of gametocytogenesis in Plasmodium falciparum in continuous culture. Bull World Health Organ. 1979, 57 (Suppl 1): 37-52.
Bruce MC, Alano P, Duthie S, Carter R: Commitment of the malaria parasite Plasmodium falciparum to sexual and asexual development. Parasitology. 1990, 100 (Pt 2): 191-200.
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.
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.
Drakeley CJ, Secka I, Correa S, Greenwood BM, Targett GA: Host haematological factors influencing the transmission of Plasmodium falciparum gametocytes to Anopheles gambiae s.s. mosquitoes. Trop Med Int Health. 1999, 4: 131-138. 10.1046/j.1365-3156.1999.00361.x.
Barnes KI, White NJ: Population biology and antimalarial resistance: The transmission of antimalarial drug resistance in Plasmodium falciparum. Acta Trop. 2005, 94: 230-240.
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-430. 10.1016/j.pt.2006.07.001.
Baird JK, Jones TR, Purnomo , Masbar S, Ratiwayanto S, Leksana B: Evidence for specific suppression of gametocytemia by Plasmodium falciparum in residents of hyperendemic Irian Jaya. Am J Trop Med Hyg. 1991, 44: 183-190.
Schellenberg D, Menendez C, Kahigwa E, Font F, Galindo C, Acosta C, Schellenberg JA, Aponte JJ, Kimario J, Urassa H, Mshinda H, Tanner M, Alonso P: African children with malaria in an area of intense Plasmodium falciparum transmission: features on admission to the hospital and risk factors for death. Am J Trop Med Hyg. 1999, 61: 431-438.
Akim NI, Drakeley C, Kingo T, Simon B, Senkoro K, Sauerwein RW: Dynamics of P. falciparum gametocytemia in symptomatic patients in an area of intense perennial transmission in Tanzania. Am J Trop Med Hyg. 2000, 63: 199-203.
Greenwood BM, Pickering H: A malaria control trial using insecticide-treated bed nets and targeted chemoprophylaxis in a rural area of The Gambia, west Africa. 1. A review of the epidemiology and control of malaria in The Gambia, west Africa. Trans R Soc Trop Med Hyg. 1993, 87 (Suppl 2): 3-11. 10.1016/0035-9203(93)90169-Q.
Thomson MC, Adiamah JH, Connor SJ, Jawara M, Bennett S, D'Alessandro U, Quinones M, Langerock P, Greenwood BM: Entomological evaluation of the Gambia's National Impregnated Bednet Programme. Ann Trop Med Parasitol. 1995, 89: 229-241.
von Seidlein L, Bojang K, Jones P, Jaffar S, Pinder M, Obaro S, Doherty T, Haywood M, Snounou G, Gemperli B, Gathmann I, Royce C, McAdam K, Greenwood B: A randomized controlled trial of artemether/benflumetol, a new antimalarial and pyrimethamine/sulfadoxine in the treatment of uncomplicated falciparum malaria in African children. Am J Trop Med Hyg. 1998, 58: 638-644.
von Seidlein L, Jawara M, Coleman R, Doherty T, Walraven G, Targett G: Parasitaemia and gametocytaemia after treatment with chloroquine, pyrimethamine/sulfadoxine, and pyrimethamine/sulfadoxine combined with artesunate in young Gambians with uncomplicated malaria. Trop Med Int Health. 2001, 6: 92-98. 10.1046/j.1365-3156.2001.00683.x.
Luxemburger C, Thwai KL, White NJ, Webster HK, Kyle DE, Maelankirri L, Chongsuphajaisiddhi T, Nosten F: The epidemiology of malaria in a Karen population on the western border of Thailand. Trans R Soc Trop Med Hyg. 1996, 90: 105-111. 10.1016/S0035-9203(96)90102-9.
Luxemburger C, Ricci F, Nosten F, Raimond D, Bathet S, White NJ: The epidemiology of severe malaria in an area of low transmission in Thailand. Trans R Soc Trop Med Hyg. 1997, 91: 256-262. 10.1016/S0035-9203(97)90066-3.
Nosten F, van Vugt M, Price R, Luxemburger C, Thway KL, Brockman A, McGready R, ter Kuile F, Looareesuwan S, White NJ: Effects of artesunate-mefloquine combination on incidence of Plasmodium falciparum malaria and mefloquine resistance in western Thailand: a prospective study. Lancet. 2000, 356: 297-302. 10.1016/S0140-6736(00)02505-8.
Price RN, Nosten F, Luxemburger C, van Vugt M, Phaipun L, Chongsuphajaisiddhi T, White NJ: Artesunate/mefloquine treatment of multi-drug resistant falciparum malaria. Trans R Soc Trop Med Hyg. 1997, 91: 574-577. 10.1016/S0035-9203(97)90032-8.
van Vugt M, Brockman A, Gemperli B, Luxemburger C, Gathmann I, Royce C, Slight T, Looareesuwan S, White NJ, Nosten F: Randomized comparison of artemether-benflumetol and artesunate-mefloquine in treatment of multidrug-resistant falciparum malaria. Antimicrob Agents Chemother. 1998, 42: 135-139.
van Vugt M, Looareesuwan S, Wilairatana P, McGready R, Villegas L, Gathmann I, Mull R, Brockman A, White NJ, Nosten F: Artemether-lumefantrine for the treatment of multidrug-resistant falciparum malaria. Trans R Soc Trop Med Hyg. 2000, 94: 545-548. 10.1016/S0035-9203(00)90082-8.
Glidden DV, Vittinghoff E: Modelling clustered survival data from multicentre clinical trials. Stat Med. 2004, 23: 369-388. 10.1002/sim.1599.
Barnes KI, Durrheim DN, Little F, Jackson A, Mehta U, Allen E, Dlamini SS, Tsoka J, Bredenkamp B, Mthembu DJ, White NJ, Sharp BL: Effect of artemether-lumefantrine policy and improved vector control on malaria burden in KwaZulu-Natal, South Africa. PLoS Med. 2005, 2: e330-10.1371/journal.pmed.0020330.
Sutherland CJ, Ord R, Dunyo S, Jawara M, Drakeley CJ, Alexander N, Coleman R, Pinder M, Walraven G, Targett GA: Reduction of malaria transmission to Anopheles mosquitoes with a six-dose regimen of co-artemether. PLoS Med. 2005, 2: e92-10.1371/journal.pmed.0020092.
Babiker HA, Schneider P, Reece SE: Gametocytes: insights gained during a decade of molecular monitoring. Trends Parasitol. 2008
Pethleart A, Prajakwong S, Suwonkerd W, Corthong B, Webber R, Curtis C: Infectious reservoir of Plasmodium infection in Mae Hong Son Province, north-west Thailand. Malar J. 2004, 3: 34-10.1186/1475-2875-3-34.
Meerman L, Ord R, Bousema JT, van Niekerk M, Osman E, Hallett R, Pinder M, Walraven G, Sutherland CJ: Carriage of chloroquine-resistant parasites and delay of effective treatment increase the risk of severe malaria in Gambian children. J Infect Dis. 2005, 192: 1651-1657. 10.1086/496887.
Tchuinkam T, Mulder B, Dechering 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.
Hallett RL, Sutherland CJ, Alexander N, Ord R, Jawara M, Drakeley CJ, Pinder M, Walraven G, Targett GA, Alloueche A: Combination therapy counteracts the enhanced transmission of drug-resistant malaria parasites to mosquitoes. Antimicrob Agents Chemother. 2004, 48: 3940-3943. 10.1128/AAC.48.10.3940-3943.2004.
Govere JM, Durrheim DN, Mngomezulu NM, Barnes K, Sharp B: Infectivity of Plasmodium falciparum gametocytes to Anopheles arabiensis after treatment with sulfadoxine-pyrimethamine. Trans R Soc Trop Med Hyg. 2003, 97: 707-708. 10.1016/S0035-9203(03)80108-6.
Barnes KI, Little F, Mabuza A, Mngomezulu N, Govere J, Durrheim D, Roper C, Watkins B, White NJ: Increased gametocytemia after treatment: an early parasitological indicator of emerging sulfadoxine-pyrimethamine resistance in falciparum malaria. J Infect Dis. 2008, 197: 1605-1613. 10.1086/587645.
Hallett RL, Dunyo S, Ord R, Jawara M, Pinder M, Randall A, Alloueche A, Walraven G, Targett GA, Alexander N: Chloroquine/sulphadoxine-pyrimethamine for gambian children with malaria: transmission to mosquitoes of multidrug-resistant Plasmodium falciparum. PLoS Clin Trials. 2006, 1: e15-10.1371/journal.pctr.0010015.
Mendez F, Herrera S, Murrain B, Gutierrez A, Moreno LA, Manzano M, Munoz A, Plowe CV: Selection of antifolate-resistant Plasmodium falciparum by sulfadoxine-pyrimethamine treatment and infectivity to Anopheles mosquitoes. Am J Trop Med Hyg. 2007, 77: 438-443.
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.
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-527.
Govere JM, la Grange JJ, Durrheim DN, Freese JA, Sharp BL, Mabuza A, Mngomezulu N, Bredenkamp BL: Sulfadoxine-pyrimethamine effectiveness against Plasmodium falciparum malaria in Mpumalanga Province, South Africa. Trans R Soc Trop Med Hyg. 1999, 93: 644-10.1016/S0035-9203(99)90082-2.
We thank Geoffrey Targett for providing trial data from The Gambia, and for helpful discussions. CJS is supported by the UK Health Protection Agency. The work in Tanzania was conducted at the Ifakara Health Research and Development Center and supported by funds from the Directorate General of Development Cooperation DGIS (NL002701) of the Dutch Ministry of Foreign Affairs: it could not have been completed without the contribution of I. Akim, T. Kingo, B. Simon, H. Mshinda and R. Sauerwein. We also thank the staff of the Shoklo Malaria Research Unit for their work and all the patients who participated in the studies.
The authors declare that they have no competing interests.
KS analysed the data and drafted the manuscript; RNP, CJS, CJD, LS, FN provided the data and contributed to the manuscript writing, NJW contributed to the manuscript writing and conceptually to data analysis.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.