Skip to main content

Maintenance of high temporal Plasmodium falciparum genetic diversity and complexity of infection in asymptomatic and symptomatic infections in Kilifi, Kenya from 2007 to 2018

Abstract

Background

High levels of genetic diversity are common characteristics of Plasmodium falciparum parasite populations in high malaria transmission regions. There has been a decline in malaria transmission intensity over 12 years of surveillance in the community in Kilifi, Kenya. This study sought to investigate whether there was a corresponding reduction in P. falciparum genetic diversity, using msp2 as a genetic marker.

Methods

Blood samples were obtained from children (< 15 years) enrolled into a cohort with active weekly surveillance between 2007 and 2018 in Kilifi, Kenya. Asymptomatic infections were defined during the annual cross-sectional blood survey and the first-febrile malaria episode was detected during the weekly follow-up. Parasite DNA was extracted and successfully genotyped using allele-specific nested polymerase chain reactions for msp2 and capillary electrophoresis fragment analysis.

Results

Based on cross-sectional surveys conducted in 2007–2018, there was a significant reduction in malaria prevalence (16.2–5.5%: P-value < 0.001), however msp2 genetic diversity remained high. A high heterozygosity index (He) (> 0.95) was observed in both asymptomatic infections and febrile malaria over time. About 281 (68.5%) asymptomatic infections were polyclonal (> 2 variants per infection) compared to 46 (56%) polyclonal first-febrile infections. There was significant difference in complexity of infection (COI) between asymptomatic 2.3 [95% confidence interval (CI) 2.2–2.5] and febrile infections 2.0 (95% CI 1.7–2.3) (P = 0.016). Majority of asymptomatic infections (44.2%) carried mixed alleles (i.e., both FC27 and IC/3D7), while FC27 alleles were more frequent (53.3%) among the first-febrile infections.

Conclusions

Plasmodium falciparum infections in Kilifi are still highly diverse and polyclonal, despite the reduction in malaria transmission in the community.

Background

Malaria is still a major public health burden in Kenya, despite the intensification of control measures that have resulted in recent reductions in morbidity and mortality. About 70% of the population is at risk of malaria infection with the Coastal and Lake Victoria endemic regions bearing the highest community prevalence of around 4.5% and 18.9%, respectively, based on microscopy [1]. Malaria control and eventual elimination is threatened by the emergence of drug resistant parasites and insecticide resistance by mosquitoes, the perennial presence of asymptomatic Plasmodium falciparum infections and highly diverse parasite populations [2, 3]. Asymptomatic infections harbour distinct parasite sub-populations, also termed clones/variants that normally undergo recombination in the mosquito mid gut during zygote formation resulting in genetically diverse parasites [4]. An individual may thus be infected with parasites of multiple genotypes from a single mosquito bite inoculation or multiple mosquito inoculations [5]. These number of distinct parasite genotypes in an individual is referred to as complexity of infection (COI).

The genetic diversity of P. falciparum and COI are correlates of malaria transmission intensity and can be used in assessing the impact of malaria control strategies [6]. Generally, studies have shown that P. falciparum parasites have a higher within-host genetic diversity in high transmission settings than in low transmission settings [7, 8]. This has led to the notion that a reduction in transmission intensity translates to a reduction in genetic diversity due to decreased chances of recombination between genetically distinct variants [9, 10]. The extensive genetic diversity of P. falciparum vaccine targets is a major hinderance in malaria vaccine development as the host immune responses may fail to recognize all the variants of an antigen. The merozoite surface protein-2 (MSP2) has been shown to be highly polymorphic and informative in genotyping parasite populations [11]. It is a glycoprotein encoded by the msp2 gene that is located on chromosome 2. It is divided into five blocks that include a highly polymorphic central block 3 and is flanked by unique variable domains and conserved N- and C-terminal domains [12, 13]. The polymorphic block 3 contains repeats that vary in number, length and sequence that are grouped into two allelic families i.e. IC/3D7 and FC27 [14] that are associated with different malaria outcomes [15, 16].

Asymptomatic infections constitute the biggest proportion of P. falciparum infections in endemic regions [17]. They result from partial immunity developed after repeated exposure to the parasite especially in endemic areas [18]. These individuals act as a reservoir for infectious parasites. They may be associated with either increased or reduced risk of symptomatic malaria [19, 20] depending on several factors, such as age, transmission intensity, COI, parasitaemia and acquisition of new clones [21].

This study investigated the temporal genetic diversity and complexity of P. falciparum infections in asymptomatic and first-febrile follow-up samples. In addition, the msp2 genetic diversity between asymptomatic and first-febrile pairs was examined. The samples were collected during the period of decline in malaria transmission in a moderate to high transmission region of Kilifi, Kenya, and msp2 gene polymorphisms assessed.

Methods

Study design

Samples from asymptomatic and febrile P. falciparum infections were collected from the Junju cohort in Kilifi, Kenya, a region of moderate to high transmission. In this cohort, 425 children are recruited at birth and followed up weekly by active clinical surveillance until the age of 15 years [22]. There are two rainy seasons per year in Kenya during which malaria transmission increases, the long rains from May to July and the short rains in October to November. Annual cross-sectional surveys were conducted in this cohort before the long rains from 2007 to 2018. During the annual cross-sectional surveys, individuals were categorized as uninfected, febrile malaria, non-malarial fever and asymptomatic P. falciparum infections based on rapid diagnostic test (RDT) and confirmed by microscopy. Asymptomatic individuals were defined as parasite positive and having: (1) an axillary temperature < 37.5 °C and no history of fever during the cross-sectional survey, (2) no recent febrile malaria episode within the month before the survey, and (3) no fever within the subsequent 7 days from the date of the survey [20]. The first-febrile episode, which is the first febrile infection detected during the weekly active surveillance after the cross-sectional survey, was defined as having ≥ 2500 parasites/µl by microscopy and a tympanic temperature > 37.5 °C based on definitions described for this cohort [22]. For this study, only microscopy positive samples were included, consequently, a total of 838 asymptomatic infections were available for genotyping at the cross-sectional survey, as well as a further 147 first-febrile infections.

Sample preparation, msp2 amplification and capillary electrophoresis

DNA was extracted from whole blood using the QIAamp® DNA mini kit (QIAGEN) according to the manufacturer’s instructions. msp2 (PF3D7_0206800) block 3 genotyping was performed using a nested PCR assay [23]. Laboratory cultured HB3 and IC/3D7 DNA were used as a positive control for FC27 and IC/3D7 alleles, respectively. The following 10 µl primary and nested PCR assay were conducted as previously described [23]. The primary PCR amplified the entire msp2 domain (forward, 5′-ATGAAGGTAATTAAAACATTGTCTATTATA-3′; reverse, 5′- CTTTGTTACCATCGGTACATTCTT-3). The nested PCR assay used fluorescently labelled oligonucleotide primers to target the msp2 allelic families: FC27 (forward, 5′- AATACTAAGAGTGTAGGTGCARATGCTCCA-3′; reverse 5′-TTTTATTTGGTGCAT TGCCAGAACTTGAAC-3′ 6-FAM) and IC/3D7 (forward, 5′- AGAAGTATGGCAGAAAGTAAKCCTYCTACT3′; reverse, 5′- GATTGTAATTCGGGGGATTCAGTTTGTTCG-3′ VIC). PCR products were visualized on 1% (w/v) agarose gels stained with RedSafe™ Nucleic Acid Staining Solution (iNtRON Biotechnology DR). Samples that failed to generate an amplicon were repeated using twice the DNA quantity. If non-amplifications persisted after the second PCR, the amplification was classified as unsuccessful. PCR fragments from each nested reaction were diluted 10 times with nuclease-free water and mixed with 9 µl of deionized formamide (Hi-Di) and 0.5 µl size standard GS-LIZ that contains 73 single-stranded DNA fragments ranging in size from 20 to 1200 bp. The solutions were transferred to 96-well Optical reaction plates and sent to the International Livestock Research Institute (ILRI) in Nairobi (Kenya) for capillary electrophoresis on the 3730xl DNA sequencer (Applied Biosystems).

msp2 data analysis

The msp2 fragment size data were analysed using GeneMapper Software version 4.0 (ThermoFisher) to determine the number of genotypes present in each sample. A fluorescent cut-off of 300 relative fluorescent units (rfu) was applied to simplify the identification of true alleles by removing the fluorescent background and non-specific low background noise [24]. Fragments were considered the same if they were within 3 bp difference in size since msp2 is a coding gene. All fragments falling within the limits of this bin were considered to belong to the same genotype. Stutter and artefact peaks were defined as peaks having a height of less than 10% the height of the true peak. Otherwise, they were considered as true peaks. COI was defined as the total number of msp2 fragment sizes in an individual infection. Samples containing both FC27 and IC/3D7 genotypes were classified as mixed infections.

Statistical analysis

The student’s t-test was used to compare mean COI between asymptomatic and first-febrile infections. Mann–Whitney U test was used to compare parasitaemia between asymptomatic and first-febrile infections. Associations between categorical variables were conducted using Fisher’s exact test. The analysis of microscopy positive data trends over time was performed using Mann–Kendall trend test function in the trend package [25]. Multivariate logistic regression models were fitted to associate asymptomatic and first-febrile infections with COI after adjusting for age, parasitaemia and microscopy positivity as a categorical variable (high from 2007 to 2012 and low from 2013 to 2018). All statistical tests were conducted in R v4.0.2 [26] and all plots were generated using the R packages ggplot2 v3.3.2 [27] and ggpubr v.0.4.0 [28]. A P-value of < 0.05 was considered statistically significant. Expected heterozygosity (He) was defined as the probability that two randomly selected variants from a population will carry different alleles. He was used to estimate msp2 allelic diversity at each time-point based on the formula below.

$$H_e=\left[\text{n}/\left(\text{n}-1\right)\right]\left[\left(1-\sum\text{P}_i^2\right)\right],$$

where n is the sample size and Pi is the frequency of ith allele in the population [29].

Results

Temporal msp2 genetic diversity

A total of 410 asymptomatic and 92 first-febrile samples were amplified and successfully genotyped from 217 children between 2007 and 2018. There were no corresponding first-febrile samples in the biobank in 2007 and in 2014 PCR amplification of asymptomatic samples were unsuccessful probably due to overdiluted samples (Table 1). The children had a mean age of 8.1 years (range: 0.7–15.0) and there was an almost equal proportion of males 50.2% (109) and females 49.8 (108). The median parasitaemia was significantly lower in asymptomatic infections 800 parasites/µl (range: 1–1,320,000) compared to first-febrile infections 28,800 parasites/µl (range, 2560–910,000) (P < 0.0001). There was a significant decline in malaria positivity rate based on microscopy (P < 0.001) in this cohort of children who aged over time (Table 1). Overall, COI was stably maintained between 2 and 3 over the 12-year period and the He values were consistently high (> 0.95) in both infections over time. More IC/3D7 alleles (129 [31.4%]) were observed in the asymptomatic infections than the FC27 alleles (101 [24%]) that were predominant in the first-febrile infections (Table 2). The sizes of these genotypes ranged from 180 to 673 bp and 315–805 bp for the FC27 and IC/3D7 allelic families, respectively. There were at least 5 FC27 alleles (291 bp, 327 bp, 362 bp, 365 and 411 bp) at a relatively high frequency (dominant alleles) that persisted over the 12-year study period out of a total of 45 FC27 alleles in asymptomatic infections (Additional file 1: Table S1). Though there was a lot more genetic variation in the IC/3D7 allelic family and only three (497 bp, 548 bp, 555 bp) IC/3D7 fragments out of 78 were persistent over time (Additional file 2: Table S2). Compared with asymptomatic infections, the first febrile infections contained fewer alleles (i.e. 29 FC27 alleles and 38 IC alleles, with allele sizes ranging from 217 to 545 bp for the FC27 and 327–724 bp for the IC/3D7 allelic families (Additional file 1: Table S1). An overlap of 19 FC27 and 32 IC/3D7 alleles between asymptomatic and first-febrile infections were detected.

Table 1 Characteristics of the cohort and number of samples successfully genotyped from 2007 to 2018
Table 2 Distribution of msp2 alleles and complexity of infections over time

Complexity of infections

Asymptomatic individuals were characterized by more (281, 68.5%) polyclonal (≥ 2) infections, with a mean COI of 2.3 (1–10) (Fig. 1A). The first-febrile infections in contrast were more monoclonal (with either a single clone of FC27 or IC/3D7 allelic types) as 46 (50%) infections were observed in the wide base of the plot and they contained a maximum of 6 clones in any infection (Fig. 1A). The spread in the proportion of polyclonal asymptomatic infections over time is depicted in Fig. 1B, with the population consistently harbouring at least 5 clones every year. The mean COI for asymptomatic infections was 2.3 (95% CI 2.2–2.5). The lowest COI was 1.7 in 2008, while the highest COI, 3.0, was observed in 2016 (Table 2). While the mean COI for first-febrile infections was 2.0 (95% CI 1.7–2.3). Overall, there was a statistically significant difference in COI between asymptomatic and first-febrile infections (P = 0.015) (Fig. 1A). Further analysis revealed that the risk of being febrile reduced by 22.9% (adjusted odds ratio (AOR): 0.771; 95% CI 0.611–0.95) for every unit increase in COI.

Fig. 1
figure 1

Plasmodium falciparum genetic diversity is high in asymptomatic infections. A Distribution of COI in asymptomatic and first-febrile infections. The COI of majority of asymptomatic samples was distributed around 2 while for first-febrile infections it was distributed around 1. B Proportion of individuals with different number of P. falciparum genotypes per every 2 years in asymptomatic infections and C in first-febrile infections. In B and C, n represents the total number of successfully genotyped samples. While the various colours represent the different number of clones: dark grey (1), orange (2), light blue (3), green (4), yellow (5), blue (6), brown (7), pink (8) and red (10)

msp2 genetic diversity in paired asymptomatic and first-febrile samples

Twenty-six children had paired genotype data from their asymptomatic and corresponding follow-up first-febrile infections. Only 2 individuals maintained one allele (the prevalent FC27 alleles in the population, 327 bp and 411 bp, Additional file 1: Table S1) between their asymptomatic and first-febrile infection (Table 3). Eight FC27 and 3 IC/3D7 genotypes were common among the paired asymptomatic and first-febrile samples out of a total of 18 FC27 and 35 IC/3D7 alleles, respectively (Table 3). In contrast, about 7 (26.9%) of the asymptomatic infections did not have an FC27 genotype. Subsequently, no association was observed between asymptomatic and first-febrile infections with the allelic family types (i.e. FC27, IC/3D7 or mixed FC27 + IC/3D7 alleles). The number of asymptomatic FC27 alleles were 6, IC/3D7 alleles 7 and mixed alleles 13 which were compared to first-febrile FC27 alleles 14, IC/3D7 alleles 2 and mixed alleles 10 (P = 0.057). However, there was a significant difference (P = 0.041) when the mixed allelic infections were excluded, since the majority of FC27 and IC/3D7 genotypes were observed in first-febrile infections and asymptomatic infections, respectively.

Table 3 msp2 gene diversity in paired asymptomatic and first-febrile infections

Discussion

Despite the decline in malaria positivity prevalence over 12 years in the community cohort, malaria is still characterized by highly genetically diverse P. falciparum infections, the stability of msp2 alleles and a high complexity of infection. This corresponds with the sustained moderate-high transmission in the study area. There was no temporal change in msp2 genetic diversity or COI, suggesting that in this moderate to high transmission area though malaria positivity rate significantly declined between 2007 and 2018, it was not substantial enough to result in a change in the parasite genetic profile. In the wider study area, Kilifi County, a significant decline in county referral hospital malaria admissions was described between 2002 and 2009 [30, 31]. The decline was not sustained and thereafter from 2009 there was an increase in hospital admission malaria positivity in older children [30]. The decline in the localized community population observed in this study, during a period of an overall rise in malaria hospital admissions [30], highlights the differences in the surveillance populations. The hospital surveillance data provides a better representation of the population since it covers a wider catchment area of the county, compared to the local community cohort analysis that is a subset of the wider county population and includes asymptomatic infections in the counts. It is possible that a genetically diverse parasite population was maintained by the sustained transmission in the county despite the local decline in the Junju area. This hypothesis is consistent with findings at a household level, where serological surveys showed evidence of diverse populations in homesteads at low malaria risk where the surrounding area was at high transmission, and vice versa evidence of less diverse populations in homesteads at high malaria risk where surrounding areas were at low transmission [32].

Thus, the extensive parasite genetic diversity is maintained. In great contrast, a dramatic reduction in malaria transmission as observed in Grande Comore Island, Union of Comoros, from 108,260 cases in 2006 to 1072 in 2015, was followed by a commensurable decline in MOI based on msp2 genotyping from 2.75 to 1.35 in healthcare facility samples obtained from 2006/2007 and 2013–2016 [9]. Furthermore, there was a significant reduction in msp2 alleles between the two time-points [9]. Altogether the msp2 genetic profile corresponded to the decline in malaria transmission, indicating COI as a marker of assessing the changes in transmission. Similarly, intensification of malaria control interventions in Senegal between 2006 and 2011 resulted into a reduction in genetic diversity of parasite populations [33]. On the contrary, reduction in malaria transmission in the Kingdom of Eswatini did not result into low parasite genetic diversity mainly due to malaria importation from neighbouring countries with high malaria transmission intensity [34]. Thus, inferring malaria transmission intensity from parasite genetic data ought to consider the impact of external factors affecting the parasite population genetics.

The apparent preference for the msp2 FC27 alleles was a significant feature of first-febrile infections in the asymptomatic-first-febrile paired analysis. This observation that has been made before in Congo and Tanzania, FC27 alleles were associated with severity of disease and were more predominant in children who had two or more febrile malaria episodes [16, 35]. Interestingly, in a case-control study conducted in Papua New-Guinea, the FC27 genotypes were twice as likely to be found in symptomatic than asymptomatic individuals [36]). The FC27 allelic family is potentially an important set of genetic variation to interrogate further to determine their impact on immunity. The IC/3D7 family has been associated with asymptomatic infections and is thought to protect against clinical malaria [16, 37, 38]. However, contradictory findings have reported that parasites carrying FC27 like alleles are more prevalent among asymptomatic carriers [15, 39]. There is no clear consensus on whether the two msp2 allelic families are likely to be found in asymptomatic or symptomatic infections. Larger studies in regions with different transmission intensities are needed to gain more insights into the effect of each allelic family on clinical outcome.

The high COI and large proportion of polyclonal asymptomatic infections is a result of the frequent and repeated exposure to genetically distinct malaria parasites in endemic areas, as described in previous studies [40]. This leads to the development of partial immunity that results in a reduction in clinical symptoms and carriage of low-level parasitaemia [41, 42]. The paired samples revealed the rapid turnover of alleles between asymptomatic and first-febrile infections, which is expected given ongoing malaria transmission in the study area. Asymptomatic P. falciparum infections can act as precursors to symptomatic malaria [43]. Genotyping of msp2 has previously been used to assess whether the development of symptoms is due to persistence of an existing clone or due to infection with a new clone [44]. In this study, the febrile infections were characterized by more monoclonal infections, an overall lower COI and new alleles unobserved in the prior asymptomatic infection. The new alleles likely escape immune responses, rapidly increasing parasitaemia thereby causing massive tissue damage that manifests as symptoms. Similar findings have been reported in other studies, implicating the lack of protective immune responses against the new clones [44,45,46,47]. Although the study used the more sensitive capillary electrophoresis to determine fragment sizes, a strict inclusion criterion was used to define true peaks during data analysis, which may have underestimated the fragment numbers impacting the estimation of COI. The presence of stutter peaks in the capillary electrophoresis data also presented technical challenges in the definition of true peaks. Future studies should consider using more sensitive methods like targeted amplicon deep sequencing (TADS) to define COI.

The high msp2 genetic diversity maintained across the study period was expected as Kilifi is a region of moderate to high malaria transmission. The 291 bp, 327 and 411 bp FC27 and 555 bp IC/3D7 fragment sizes were common in both asymptomatic and first-febrile infections. Strikingly, some of these genotypes have been reported in other countries, such as Mali [48], as the most common genotypes, suggesting that they can be selected as candidates for malaria vaccine development. However, identical fragment lengths may not always represent identical sequence lengths and sequencing is required for confirmation.

Conclusions

Malaria surveillance should also focus on asymptomatic infections, in addition to symptomatic infections, given the extensive genetic diversity and the impact they have on sustaining malaria transmission. Similar studies should be conducted to monitor the trends in parasite genetic diversity to associate this with changes in malaria transmission.

Availability of data and materials

The datasets supporting the conclusion of this article are available in the Harvard Dataverse repository: https://doi.org/10.7910/DVN/3UE1NB.

Abbreviations

msp2 :

Merozoite surface protein 2

COI:

Complexity of infections

bp:

Base pairs

PCR:

Polymerase chain reaction

He:

Expected heterozygosity

References

  1. Kenya National Bureau of Statistics. Kenya malaria indicator survey 2020. Kenya: Nairobi; 2020.

    Google Scholar 

  2. Phillips MA, Burrows JN, Manyando C, Van Huijsduijnen RH, Van Voorhis WC, Wells TNC. Malaria. Nat Rev Dis Prim. 2017;3:1–24.

    Google Scholar 

  3. Alonso PL, Brown G, Arevalo-Herrera M, Binka F, Chitnis C, Collins F, et al. A research agenda to underpin malaria eradication. PLoS Med Med. 2011;8:e1000406.

    Article  Google Scholar 

  4. Babiker HA, Ranford-Cartwright LC, Currie D, Charlwood JD, Billingsley P, Teuscher T, et al. Random mating in a natural population of the malaria parasite Plasmodium falciparum. Parasitology. 1994;109:413.

    Article  PubMed  Google Scholar 

  5. Nkhoma SC, Trevino SG, Gorena KM, Nair S, Khoswe S, Jett C, et al. Co-transmission of related malaria parasite lineages shapes within-host parasite diversity. Cell Host Microbe. 2020;27:93-103.e4.

    Article  CAS  PubMed  Google Scholar 

  6. Apinjoh TO, Ouattara A, Titanji VPK, Djimde A, Amambua-Ngwa A. Genetic diversity and drug resistance surveillance of Plasmodium falciparum for malaria elimination: is there an ideal tool for resource-limited sub-Saharan Africa? Malar J. 2019;18:217.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Auburn S, Campino S, Miotto O, Djimde AA, Zongo I, Manske M, et al. Characterization of within-host Plasmodium falciparum diversity using next-generation sequence data. PLoS ONE. 2012;7:e32891.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  9. Huang B, Tuo F, Liang Y, Wu W, Wu G, Huang S, et al. Temporal changes in genetic diversity of msp-1, msp-2, and msp-3 in Plasmodium falciparum isolates from Grande Comore Island after introduction of ACT. Malar J. 2018;17:83.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Anderson TJC, Haubold B, Williams JT, Estrada-Franco JG, Richardson L, Mollinedo R, et al. Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum. Mol Biol Evol. 2000;17:1467–82.

    Article  CAS  PubMed  Google Scholar 

  11. Owuso-Afyei S, Poku Asante K, Adjuik M, Adjei G, Awini E, Chandramohan D. Epidemiology of malaria in the forest-savanna transitional zone of Ghana. Malar J. 2009;8:220.

    Article  Google Scholar 

  12. Takala S, Escalante A, Branch O, Kariuki S, Biswas S, Chaiyaroj S, et al. Genetic diversity in the Block 2 region of the merozoite surface protein 1 (MSP-1) of Plasmodium falciparum: additional complexity and selection and convergence in fragment size polymorphism. Infect Genet Evol. 2006;6:417–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Ferreira MU, Hartl DL. Plasmodium falciparum: worldwide sequence diversity and evolution of the malaria vaccine candidate merozoite surface protein-2 (MSP-2). Exp Parasitol. 2007;115:32–40.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Gnagne AP, Konate A, Bedia-Tanoh AV, Amiah-Droh M, Menan HIE, N’Guetta ASP, et al. Dynamics of Plasmodium falciparum genetic diversity among asymptomatic and symptomatic children in three epidemiological areas in Cote d’Ivoire. Pathog Glob Health. 2019;113:133–42.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Kidima W, Nkwengulila G. Plasmodium falciparum msp2 genotypes and multiplicity of infections among children under five years with uncomplicated malaria in Kibaha, Tanzania. J Parasitol Res. 2015;2015:721201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  18. Bull PC, Marsh K. The role of antibodies to Plasmodium falciparum-infected-erythrocyte surface antigens in naturally acquired immunity to malaria. Trends Microbiol. 2002;10:55–8.

    Article  CAS  PubMed  Google Scholar 

  19. Rono J, Osier FHA, Olsson D, Montgomery S, Mhoja L, Rooth I, et al. Breadth of anti-merozoite antibody responses is associated with the genetic diversity of asymptomatic Plasmodium falciparum infections and protection against clinical malaria. Clin Infect Dis. 2013;57:1409–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Wamae K, Wambua J, Nyangweso G, Mwambingu G, Osier F, Ndung’u F, et al. Transmission and age impact the risk of developing febrile malaria in children with asymptomatic Plasmodium falciparum parasitemia. J Infect Dis. 2018;219:936–44.

    Article  PubMed Central  Google Scholar 

  21. Kimenyi KM, Wamae K, Ochola-Oyier LI. Understanding P. falciparum asymptomatic infections: a proposition for a transcriptomic approach. Front Immunol. 2019;10:2398.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Mwangi TW, Ross A, Snow RW, Marsh K. Case definitions of clinical malaria under different transmission conditions in Kilifi District, Kenya. J Infect Dis. 2005;191:1932–9.

    Article  PubMed  Google Scholar 

  23. Liljander A, Wiklund L, Falk N, Kweku M, Martensson A, Felger I, et al. Optimization and validation of multi-coloured capillary electrophoresis for genotyping of Plasmodium falciparum merozoite surface proteins (msp1 and 2). Malar J. 2009;8:78.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Falk N, Maire N, Sama W, Owusu-Agyei S, Smith T, Beck HP, et al. Comparison of PCR-RFLP and genescan-based genotyping for analyzing infection dynamics of Plasmodium falciparum. Am J Trop Med Hyg. 2006;74:944–50.

    Article  CAS  PubMed  Google Scholar 

  25. Pohlert T. Non-parametric trend tests and change-point detection. R package version 1.1.4. 2020.

  26. R Core Team. R: a language and environment for statistical computing. Vienna, Austria; 2020.

  27. Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2016.

    Book  Google Scholar 

  28. Kassambara A. ggpubr: “ggplot2” based publication ready plots. R package version 0.4.0; 2020.

  29. Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978;89:583–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Mogeni P, Williams TN, Fegan G, Nyundo C, Bauni E, Mwai K, et al. Age, spatial, and temporal variations in hospital admissions with malaria in Kilifi County, Kenya: a 25-year longitudinal observational study. PLoS Med. 2016;13:e1002047.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Njuguna P, Maitland K, Nyaguara A, Mwanga D, Mogeni P, Mturi N, et al. Observational study: 27 years of severe malaria surveillance in Kilifi, Kenya. BMC Med. 2019;17:124.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Bejon P, Turner L, Lavstsen T, Cham G, Olotu A, Drakeley CJ, et al. Serological evidence of discrete spatial clusters of Plasmodium falciparum parasites. PLoS ONE. 2011;6:e21711.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Daniels R, Chang HH, Séne PD, Park DC, Neafsey DE, Schaffner SF, et al. Genetic surveillance detects both clonal and epidemic transmission of malaria following enhanced intervention in Senegal. PLoS ONE. 2013;8:e60780.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Roh ME, Tessema SK, Murphy M, Nhlabathi N, Mkhonta N, Vilakati S, et al. High genetic diversity of Plasmodium falciparum in the low-transmission setting of the Kingdom of Eswatini. J Infect Dis. 2019;220:1346–54.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Ibara-Okabande R, Koukouikila-Koussounda F, Ndounga M, Vouvoungui J, Malonga V, Casimiro PN, et al. Reduction of multiplicity of infections but no change in msp2 genetic diversity in Plasmodium falciparum isolates from Congolese children after introduction of artemisinin-combination therapy. Malar J. 2012;11:410.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Engelbrecht F, Felger I, Genton B, Alpers M, Beck FP. Plasmodium falciparum: malaria morbidity is associated with specific merozoite surface antigen 2 genotypes. Exp Parasitol. 1995;81:90–6.

    Article  CAS  PubMed  Google Scholar 

  37. Abukari Z, Okonu R, Nyarko SB, Lo AC, Dieng CC, Salifu SP, et al. The diversity, multiplicity of infection and population structure of P. falciparum parasites circulating in asymptomatic carriers living in high and low malaria transmission settings of Ghana. Genes (Basel). 2019;10:434.

    Article  CAS  Google Scholar 

  38. Botwe AK, Asante KP, Adjei G, Assafuah S, Dosoo D, Owusu-Agyei S. Dynamics in multiplicity of Plasmodium falciparum infection among children with asymptomatic malaria in central Ghana. BMC Genet. 2017;18:67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Amodu OK, Oyedeji SI, Ntoumi F, Orimadegun AE, Gbadegesin RA, Olumese PE, et al. Complexity of the msp2 locus and the severity of childhood malaria, in south-western Nigeria. Ann Trop Med Parasitol. 2008;102:95–102.

    Article  CAS  PubMed  Google Scholar 

  40. Touray AO, Mobegi VA, Wamunyokoli F, Herren JK. Diversity and multiplicity of P. falciparum infections among asymptomatic school children in Mbita, Western Kenya. Sci Rep. 2020;10:5924.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Bereczky S, Liljander A, Rooth I, Faraja L, Granath F, Montgomery SM, et al. Multiclonal asymptomatic Plasmodium falciparum infections predict a reduced risk of malaria disease in a Tanzanian population. Microbes Infect. 2007;9:103–10.

    Article  PubMed  Google Scholar 

  42. Smith T, Felger I, Tanner M, Beck HP. Premunition in Plasmodium falciparum infection: insights from the epidemiology of multiple infections. Trans R Soc Trop Med Hyg. 1999;93(Suppl 1):59–64.

    Article  PubMed  Google Scholar 

  43. Galatas B, Bassat Q, Mayor A. Malaria parasites in the asymptomatic: looking for the hay in the haystack. Trends Parasitol. 2016;32:296–308.

    Article  PubMed  Google Scholar 

  44. Kun JFJ, Missinou MA, Lell B, Sovric M, Knoop H, Bojowald B, et al. New emerging Plasmodium falciparum genotypes in children during the transition phase from asymptomatic parasitemia to malaria. Am J Trop Med Hyg. 2002;66:653–8.

    Article  PubMed  Google Scholar 

  45. Contamin H, Fandeur T, Rogier C, Bonnefoy S, Konate L, Trape J-F, et al. Different genetic characteristics of Plasmodium falciparum isolates collected during successive clinical malaria episodes in Senegalese children. Am J Trop Med Hyg. 1996;54:632–43.

    Article  CAS  PubMed  Google Scholar 

  46. Nsobya SL, Parikh S, Kironde F, Lubega G, Kamya MR, Rosenthal PJ, et al. Molecular evaluation of the natural history of asymptomatic parasitemia in Ugandan children. J Infect Dis. 2004;189:2220–6.

    Article  PubMed  Google Scholar 

  47. Roper C, Richardson W, Elhassan IM, Giha H, Hviid L, Satti GMH, et al. Seasonal changes in the Plasmodium falciparum population in individuals and their relationship to clinical malaria: a longitudinal study in a Sudanese village. Parasitology. 1998;116:501–10.

    Article  CAS  PubMed  Google Scholar 

  48. Sondén K, Doumbo S, Hammar U, Vafa Homann M, Ongoiba A, Traoré B, et al. Asymptomatic multiclonal Plasmodium falciparum infections carried through the dry season predict protection against subsequent clinical malaria. J Infect Dis. 2015;212:608–16.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

The team acknowledges the study participants and sincerely appreciates to their parents.

Funding

This work was supported through the DELTAS Africa Initiative (Grant No. DEL-15-003), an independent funding scheme of the African Academy of Sciences (AAS) Alliance for Accelerating Excellence in Science in Africa (AESA) and supported by the New Partnership for Africa’s Development Planning and Coordinating Agency (NEPAD Agency) with funding from the Wellcome Trust and the UK government, DELGEME Grant 107740/Z/15/Z to support KK and IDeAL Grant No. 107769/Z/10/Z to support KW. The views expressed in this publication are those of the author(s) and not necessarily those of AAS, NEPAD Agency, Wellcome Trust or the UK government. There was also support from Wellcome Trust Intermediate Fellowships, Grant Numbers 107568/Z/15/Z and 209289/Z/17/Z awarded to LIO-O and AA, respectively.

Author information

Authors and Affiliations

Authors

Contributions

KMK wrote the manuscript, generated and analysed the data, and interpreted the results. LIO-O conceptualized the work, interpreted results and wrote the manuscript. KW analysed data, interpreted results and wrote the manuscript. LM, ZL, VO and JN generated data. AA and GO were involved in the conceptualization of the work and reviewed the manuscript. PB was involved in the conceptualization of the work, interpreted the results and wrote the manuscript. All authors reviewed the final manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kelvin M. Kimenyi.

Ethics declarations

Ethical approval and consent to partcipate

For this study was obtained from the ethics review committee of the Kenya Medical Research Institute under protocol number SERU 3149. Informed consent was obtained from parents/guardians of all study participants before sample collection.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1.

Temporal changes in FC27 allele frequencies across asymptomatic and febrile infections in Kilifi, Kenya from 2007 to 2018. n refers to the number of successfully genotyped individuals per year. m denotes the total number of genotypes per year.

Additional file 2: Table S2.

Temporal changes in IC/3D7 allele frequencies across asymptomatic and febrile infections in Kilifi, Kenya from 2007 to 2018. n refers to the number of successfully genotyped individuals per year. m denotes the total number of genotypes per year.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kimenyi, K.M., Wamae, K., Ngoi, J.M. et al. Maintenance of high temporal Plasmodium falciparum genetic diversity and complexity of infection in asymptomatic and symptomatic infections in Kilifi, Kenya from 2007 to 2018. Malar J 21, 192 (2022). https://doi.org/10.1186/s12936-022-04213-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-022-04213-7

Keywords