Estimating Geographical Variation in the Risk of Zoonotic Plasmodium knowlesi Infection in Countries Eliminating Malaria

BackgroundInfection by the simian malaria parasite, Plasmodium knowlesi, can lead to severe and fatal disease in humans, and is the most common cause of malaria in parts of Malaysia. Despite being a serious public health concern, the geographical distribution of P. knowlesi malaria risk is poorly understood because the parasite is often misidentified as one of the human malarias. Human cases have been confirmed in at least nine Southeast Asian countries, many of which are making progress towards eliminating the human malarias. Understanding the geographical distribution of P. knowlesi is important for identifying areas where malaria transmission will continue after the human malarias have been eliminated.Conclusions/SignificanceWe have produced the first map of P. knowlesi malaria risk, at a fine-scale resolution, to identify priority areas for surveillance based on regions with sparse data and high estimated risk. Our map provides an initial evidence base to better understand the spatial distribution of this disease and its potential wider contribution to malaria incidence. Considering malaria elimination goals, areas for prioritised surveillance are identified.     

Prevalence of simian malaria parasites in macaques of Singapore

Plasmodium knowlesi is a simian malaria parasite currently recognized as the fifth causative agent of human malaria. Recently, naturally acquired P. cynomolgi infection in humans was also detected in Southeast Asia. The main reservoir of both parasites is the long-tailed and pig-tailed macaques, which are indigenous in this region. Due to increased urbanization and changes in land use, there has been greater proximity and interaction between the long-tailed macaques and the general population in Singapore. As such, this study aims to determine the prevalence of simian malaria parasites in local macaques to assess the risk of zoonosis to the general human population. Screening for the presence of malaria parasites was conducted on blood samples from 660 peridomestic macaques collected between Jan 2008 and Mar 2017, and 379 wild macaques collected between Mar 2009 and Mar 2017, using a Pan-Plasmodium-genus specific PCR. Positive samples were then screened using a simian Plasmodium species-specific nested PCR assay to identify the species of parasites (P. knowlesi, P. coatneyi, P. fieldi, P. cynomolgi, and P. inui) present. All the peridomestic macaques sampled were tested negative for malaria, while 80.5% of the 379 wild macaques were infected. All five simian Plasmodium species were detected; P. cynomolgi being the most prevalent (71.5%), followed by P. knowlesi (47.5%), P. inui (42.0%), P. fieldi (32.5%), and P. coatneyi (28.5%). Co-infection with multiple species of Plasmodium parasites was also observed. The study revealed that Singapore’s wild long-tailed macaques are natural hosts of the five simian malaria parasite species, while no malaria was detected in all peridomestic macaques tested. Therefore, the risk of simian malaria transmission to the general human population is concluded to be low. However, this can be better demonstrated with the incrimination of the vectors of simian malaria parasites in Singapore.     

Disease Progression in Plasmodium knowlesi Malaria Is Linked to Variation in Invasion Gene Family Members

Emerging pathogens undermine initiatives to control the global health impact of infectious diseases. Zoonotic malaria is no exception. Plasmodium knowlesi, a malaria parasite of Southeast Asian macaques, has entered the human population. P. knowlesi, like Plasmodium falciparum, can reach high parasitaemia in human infections, and the World Health Organization guidelines for severe malaria list hyperparasitaemia among the measures of severe malaria in both infections. Not all patients with P. knowlesi infections develop hyperparasitaemia, and it is important to determine why. Between isolate variability in erythrocyte invasion, efficiency seems key. Here we investigate the idea that particular alleles of two P. knowlesi erythrocyte invasion genes, P. knowlesi normocyte binding protein Pknbpxa and Pknbpxb, influence parasitaemia and human disease progression. Pknbpxa and Pknbpxb reference DNA sequences were generated from five geographically and temporally distinct P. knowlesi patient isolates. Polymorphic regions of each gene (approximately 800 bp) were identified by haplotyping 147 patient isolates at each locus. Parasitaemia in the study cohort was associated with markers of disease severity including liver and renal dysfunction, haemoglobin, platelets and lactate, (r = ≥0.34, p = <0.0001 for all). Seventy-five and 51 Pknbpxa and Pknbpxb haplotypes were resolved in 138 (94%) and 134 (92%) patient isolates respectively. The haplotypes formed twelve Pknbpxa and two Pknbpxb allelic groups. Patients infected with parasites with particular Pknbpxa and Pknbpxb alleles within the groups had significantly higher parasitaemia and other markers of disease severity. Our study strongly suggests that P. knowlesi invasion gene variants contribute to parasite virulence. We focused on two invasion genes, and we anticipate that additional virulent loci will be identified in pathogen genome-wide studies. The multiple sustained entries of this diverse pathogen into the human population must give cause for concern to malaria elimination strategists in the Southeast Asian region.     

Genetic Diversity and Natural Selection of the Plasmodium knowlesi Circumsporozoite Protein Nonrepeat Regions

Background

Plasmodium knowlesi is a simian malaria parasite that has been identified to cause malaria in humans. To date, several thousand cases of human knowlesi malaria have been reported around Southeast Asia. Thus far, there is no detailed study on genetic diversity and natural selection of P. knowlesi circumsporozoite protein (CSP), a prominent surface antigen on the sporozoite of the parasite. In the present study, the genetic diversity and natural selection acting on the nonrepeat regions of the gene encoding P. knowlesi CSP were investigated, focusing on the T-cell epitope regions at the C-terminal of the protein.

Methods

Blood samples from 32 knowlesi malaria patients and 2 wild monkeys (Macaca fascicularis) were used. The CSP of the P. knowlesi isolates was amplified by PCR, cloned into Escherichia coli, and sequenced. The nonrepeat regions of the CSP gene were analysed for genetic diversity, natural selection and haplotypic grouping using MEGA5 and DnaSP version 5.10.00 programmes. A haplotype network was constructed based on the C-terminal (Th2R/Th3R) T-cell epitope regions using the Median-Joining method in the NETWORK version 4.6.1.2 programme. Previously published sequences from other regions (Malaysia Borneo, Singapore) were also included in the analysis.

Results

A total of 123 P. knowlesi CSP sequences were analysed. Multiple sequence alignment revealed 58 amino acid changes, and 42 novel amino acid haplotypes were identified. Polymorphism was higher in the C-terminal Th2R/Th3R epitope (π = 0.0293, n = 123) region compared to the overall combined nonrepeat regions (π = 0.0120, n = 123). Negative natural selection was observed within the nonrepeat regions of the CSP gene. Within the C-terminal Th2R/Th3R epitope regions, there was evidence of slight positive selection. Based on haplotype network analysis of the Th2R/Th3R regions, five abundant haplotypes were identified. Sharing of haplotypes between humans and macaques were observed.

Conclusion

This study contributes to the understanding of the type and distribution of naturally occurring polymorphism in the P. knowlesi CSP gene. This study also provides a measurement of the genetic diversity of P. knowlesi and identifies the predominant haplotypes within Malaysia based on the C-terminal Th2R/Th3R regions.     

DNA replication dynamics during erythrocytic schizogony in the malaria parasites Plasmodium falciparum and Plasmodium knowlesi

Malaria parasites are unusual, early-diverging protozoans with non-canonical cell cycles. They do not undergo binary fission, but divide primarily by schizogony. This involves the asynchronous production of multiple nuclei within the same cytoplasm, culminating in a single mass cytokinesis event. The rate and efficiency of parasite reproduction is fundamentally important to malarial disease, which tends to be severe in hosts with high parasite loads. Here, we have studied for the first time the dynamics of schizogony in two human malaria parasite species, Plasmodium falciparum and Plasmodium knowlesi. These differ in their cell-cycle length, the number of progeny produced and the genome composition, among other factors. Comparing them could therefore yield new information about the parameters and limitations of schizogony. We report that the dynamics of schizogony differ significantly between these two species, most strikingly in the gap phases between successive nuclear multiplications, which are longer in P. falciparum and shorter, but more heterogenous, in P. knowlesi. In both species, gaps become longer as schizogony progresses, whereas each period of active DNA replication grows shorter. In both species there is also extreme variability between individual cells, with some schizonts producing many more nuclei than others, and some individual nuclei arresting their DNA replication for many hours while adjacent nuclei continue to replicate. The efficiency of schizogony is probably influenced by a complex set of factors in both the parasite and its host cell.     

Plasmodium knowlesi infecting humans in Southeast Asia: What’s next?

Plasmodium knowlesi, a simian malaria parasite, has been in the limelight since a large focus of human P. knowlesi infection was reported from Sarawak (Malaysian Borneo) in 2004. Although this infection is transmitted across Southeast Asia, the largest number of cases has been reported from Malaysia. The increasing number of knowlesi malaria cases has been attributed to the use of molecular tools for detection, but environmental changes including deforestation likely play a major role by increasing human exposure to vector mosquitoes, which coexist with the macaque host. In addition, with the reduction in human malaria transmission in Southeast Asia, it is possible that human populations are at a greater risk of P. knowlesi infection due to diminishing cross-species immunity. Furthermore, the possibility of increasing exposure of humans to other simian Plasmodium parasites such as Plasmodium cynomolgi and Plasmodium inui should not be ignored. We here review the current status of these parasites in humans, macaques, and mosquitoes to support necessary reorientation of malaria control and elimination in the affected areas.     

Admixture in Humans of Two Divergent Plasmodium knowlesi Populations Associated with Different Macaque Host Species

Human malaria parasite species were originally acquired from other primate hosts and subsequently became endemic, then spread throughout large parts of the world. A major zoonosis is now occurring with Plasmodium knowlesi from macaques in Southeast Asia, with a recent acceleration in numbers of reported cases particularly in Malaysia. To investigate the parasite population genetics, we developed sensitive and species-specific microsatellite genotyping protocols and applied these to analysis of samples from 10 sites covering a range of >1,600 km within which most cases have occurred. Genotypic analyses of 599 P. knowlesi infections (552 in humans and 47 in wild macaques) at 10 highly polymorphic loci provide radical new insights on the emergence. Parasites from sympatric long-tailed macaques (Macaca fascicularis) and pig-tailed macaques (M. nemestrina) were very highly differentiated (FST = 0.22, and K-means clustering confirmed two host-associated subpopulations). Approximately two thirds of human P. knowlesi infections were of the long-tailed macaque type (Cluster 1), and one third were of the pig-tailed-macaque type (Cluster 2), with relative proportions varying across the different sites. Among the samples from humans, there was significant indication of genetic isolation by geographical distance overall and within Cluster 1 alone. Across the different sites, the level of multi-locus linkage disequilibrium correlated with the degree of local admixture of the two different clusters. The widespread occurrence of both types of P. knowlesi in humans enhances the potential for parasite adaptation in this zoonotic system.     

Low Levels of Polymorphisms and No Evidence for Diversifying Selection on the Plasmodium knowlesi Apical Membrane Antigen 1 Gene

Infection with Plasmodium knowlesi, a zoonotic primate malaria, is a growing human health problem in Southeast Asia. P. knowlesi is being used in malaria vaccine studies, and a number of proteins are being considered as candidate malaria vaccine antigens, including the Apical Membrane Antigen 1 (AMA1). In order to determine genetic diversity of the ama1 gene and to identify epitopes of AMA1 under strongest immune selection, the ama1 gene of 52 P. knowlesi isolates derived from human infections was sequenced. Sequence analysis of isolates from two geographically isolated regions in Sarawak showed that polymorphism in the protein is low compared to that of AMA1 of the major human malaria parasites, P. falciparum and P. vivax. Although the number of haplotypes was 27, the frequency of mutations at the majority of the polymorphic positions was low, and only six positions had a variance frequency higher than 10%. Only two positions had more than one alternative amino acid. Interestingly, three of the high-frequency polymorphic sites correspond to invariant sites in PfAMA1 or PvAMA1. Statistically significant differences in the quantity of three of the six high frequency mutations were observed between the two regions. These analyses suggest that the pkama1 gene is not under balancing selection, as observed for pfama1 and pvama1, and that the PkAMA1 protein is not a primary target for protective humoral immune responses in their reservoir macaque hosts, unlike PfAMA1 and PvAMA1 in humans. The low level of polymorphism justifies the development of a single allele PkAMA1-based vaccine.     

Recognition of Human Erythrocyte Receptors by the Tryptophan-Rich Antigens of Monkey Malaria Parasite Plasmodium knowlesi

Background

The monkey malaria parasite Plasmodium knowlesi also infect humans. There is a lack of information on the molecular mechanisms that take place between this simian parasite and its heterologous human host erythrocytes leading to this zoonotic disease. Therefore, we investigated here the binding ability of P. knowlesi tryptophan-rich antigens (PkTRAgs) to the human erythrocytes and sharing of the erythrocyte receptors between them as well as with other commonly occurring human malaria parasites.

Methods

Six PkTRAgs were cloned and expressed in E.coli as well as in mammalian CHO-K1 cell to determine their human erythrocyte binding activity by cell-ELISA, and in-vitro rosetting assay, respectively.

Results

Three of six PkTRAgs (PkTRAg38.3, PkTRAg40.1, and PkTRAg67.1) showed binding to human erythrocytes. Two of them (PkTRAg40.1 and PkTRAg38.3) showed cross-competition with each other as well as with the previously described P.vivax tryptophan-rich antigens (PvTRAgs) for human erythrocyte receptors. However, the third protein (PkTRAg67.1) utilized the additional but different human erythrocyte receptor(s) as it did not cross-compete for erythrocyte binding with either of these two PkTRAgs as well as with any of the PvTRAgs. These three PkTRAgs also inhibited the P.falciparum parasite growth in in-vitro culture, further indicating the sharing of human erythrocyte receptors by these parasite species and the biological significance of this receptor-ligand interaction between heterologous host and simian parasite.

Conclusions

Recognition and sharing of human erythrocyte receptor(s) by PkTRAgs with human parasite ligands could be part of the strategy adopted by the monkey malaria parasite to establish inside the heterologous human host.     

Congenital Plasmodium vivax malaria mimicking neonatal sepsis: a case report

A recently published comment on a report of Plasmodium knowlesi infections in Vietnam states that this may not accurately represent the situation in the study area because the PCR primers used may cross-hybridize with Plasmodium vivax. Nevertheless, P. knowlesi infections have been confirmed by sequencing. In addition, a neighbour-joining tree based on the 18S S-Type SSUrRNA gene shows that the Vietnamese samples clearly cluster with the P. knowlesi isolates identified in Malaysia and are distinct from the corresponding P. vivax sequences. All samples came from asymptomatic individuals who did not consult for fever during the months preceding or following the survey, indicating that asymptomatic P. knowlesi infections occur in this population, although this does not exclude the occurrence of symptomatic cases. Large-scale studies to determine the extent and the epidemiology of P. knowlesi malaria in Vietnam are further needed.     

Population genetic analysis of Plasmodium falciparum parasites using a customized Illumina GoldenGate genotyping assay

The diversity in the Plasmodium falciparum genome can be used to explore parasite population dynamics, with practical applications to malaria control. The ability to identify the geographic origin and trace the migratory patterns of parasites with clinically important phenotypes such as drug resistance is particularly relevant. With increasing single-nucleotide polymorphism (SNP) discovery from ongoing Plasmodium genome sequencing projects, a demand for high SNP and sample throughput genotyping platforms for large-scale population genetic studies is required. Low parasitaemias and multiple clone infections present a number of challenges to genotyping P. falciparum. We addressed some of these issues using a custom 384-SNP Illumina GoldenGate assay on P. falciparum DNA from laboratory clones (long-term cultured adapted parasite clones), short-term cultured parasite isolates and clinical (non-cultured isolates) samples from East and West Africa, Southeast Asia and Oceania. Eighty percent of the SNPs (n = 306) produced reliable genotype calls on samples containing as little as 2 ng of total genomic DNA and on whole genome amplified DNA. Analysis of artificial mixtures of laboratory clones demonstrated high genotype calling specificity and moderate sensitivity to call minor frequency alleles. Clear resolution of geographically distinct populations was demonstrated using Principal Components Analysis (PCA), and global patterns of population genetic diversity were consistent with previous reports. These results validate the utility of the platform in performing population genetic studies of P. falciparum.     

Crystal Structure of Plasmodium knowlesi Apical Membrane Antigen 1 and Its Complex with an Invasion-Inhibitory Monoclonal Antibody

The malaria parasite Plasmodium knowlesi, previously associated only with infection of macaques, is now known to infect humans as well and has become a significant public health problem in Southeast Asia. This species should therefore be targeted in vaccine and therapeutic strategies against human malaria. Apical Membrane Antigen 1 (AMA1), which plays a role in Plasmodium merozoite invasion of the erythrocyte, is currently being pursued in human vaccine trials against P. falciparum. Recent vaccine trials in macaques using the P. knowlesi orthologue PkAMA1 have shown that it protects against infection by this parasite species and thus should be developed for human vaccination as well. Here, we present the crystal structure of Domains 1 and 2 of the PkAMA1 ectodomain, and of its complex with the invasion-inhibitory monoclonal antibody R31C2. The Domain 2 (D2) loop, which is displaced upon binding the Rhoptry Neck Protein 2 (RON2) receptor, makes significant contacts with the antibody. R31C2 inhibits binding of the Rhoptry Neck Protein 2 (RON2) receptor by steric blocking of the hydrophobic groove and by preventing the displacement of the D2 loop which is essential for exposing the complete binding site on AMA1. R31C2 recognizes a non-polymorphic epitope and should thus be cross-strain reactive. PkAMA1 is much less polymorphic than the P. falciparum and P. vivax orthologues. Unlike these two latter species, there are no polymorphic sites close to the RON2-binding site of PkAMA1, suggesting that P. knowlesi has not developed a mechanism of immune escape from the host’s humoral response to AMA1.     

Seasonal and Spatial Dynamics of the Primary Vector of Plasmodium knowlesi within a Major Transmission Focus in Sabah,Malaysia

Background

The simian malaria parasite Plasmodium knowlesi is emerging as a public health problem in Southeast Asia, particularly in Malaysian Borneo where it now accounts for the greatest burden of malaria cases and deaths. Control is hindered by limited understanding of the ecology of potential vector species.

Methodology/Principal Findings

We conducted a one year longitudinal study of P. knowlesi vectors in three sites within an endemic area of Sabah, Malaysia. All mosquitoes were captured using human landing catch. Anopheles mosquitoes were dissected to determine, oocyst, sporozoites and parous rate. Anopheles balabacensis is confirmed as the primary vector of. P. knowlesi (using nested PCR) in Sabah for the first time. Vector densities were significantly higher and more seasonally variable in the village than forest or small scale farming site. However An. balabacensis survival and P. knowlesi infection rates were highest in forest and small scale farm sites. Anopheles balabacensis mostly bites humans outdoors in the early evening between 1800 to 2000hrs.

Conclusions/Significance

This study indicates transmission is unlikely to be prevented by bednets. This combined with its high vectorial capacity poses a threat to malaria elimination programmes within the region.     

Genetic Analysis and Species Specific Amplification of the Artemisinin Resistance-Associated Kelch Propeller Domain in P. falciparum and P. vivax

Plasmodium falciparum resistance to artemisinin has emerged in the Greater Mekong Subregion and now poses a threat to malaria control and prevention. Recent work has identified mutations in the kelch propeller domain of the P. falciparum K13 gene to be associated artemisinin resistance as defined by delayed parasite clearance and ex vivo ring stage survival assays. Species specific primers for the two most prevalent human malaria species, P. falciparum and P. vivax, were designed and tested on multiple parasite isolates including human, rodent, and non- humans primate Plasmodium species. The new protocol described here using the species specific primers only amplified their respective species, P. falciparum and P. vivax, and did not cross react with any of the other human malaria Plasmodium species. We provide an improved species specific PCR and sequencing protocol that could be effectively used in areas where both P. falciparum and P. vivax are circulating. To design this improved protocol, the kelch gene was analyzed and compared among different species of Plasmodium. The kelch propeller domain was found to be highly conserved across the mammalian Plasmodium species.     

A case of severe Plasmodium knowlesi in a splenectomized patient

Plasmodium knowlesi, a zoonotic malaria, is now considered the fifth species of Plasmodium causing malaria in humans. With its 24-hour erythrocytic stage of development, it has raised concern regarding its high potential in replicating and leading to severe illness. Spleen is an important site for removal of parasitized red blood cells and generating immunity. We reported a case of knowlesi malaria in a non-immune, splenectomized patient. We observed the delay in parasite clearance, high parasitic counts, and severe illness at presentation. A thorough search through literature revealed several case reports on falciparum and vivax malaria in splenectomized patients. However, literature available for knowlesi malaria in splenectomized patient is limited. Further studies need to be carried out to clarify the role of spleen in host defense against human malaria especially P. knowlesi.     

Hybrid selection for sequencing pathogen genomes from clinical samples

We have adapted a solution hybrid selection protocol to enrich pathogen DNA in clinical samples dominated by human genetic material. Using mock mixtures of human and Plasmodium falciparum malaria parasite DNA as well as clinical samples from infected patients, we demonstrate an average of approximately 40-fold enrichment of parasite DNA after hybrid selection. This approach will enable efficient genome sequencing of pathogens from clinical samples, as well as sequencing of endosymbiotic organisms such as Wolbachia that live inside diverse metazoan phyla.     

Dihydrofolate-Reductase Mutations in Plasmodium knowlesi Appear Unrelated to Selective Drug Pressure from Putative Human-To-Human Transmission in Sabah,Malaysia

Background

Malaria caused by zoonotic Plasmodium knowlesi is an emerging threat in Eastern Malaysia. Despite demonstrated vector competency, it is unknown whether human-to-human (H-H) transmission is occurring naturally. We sought evidence of drug selection pressure from the antimalarial sulfadoxine-pyrimethamine (SP) as a potential marker of H-H transmission.

Methods

The P. knowlesi dihdyrofolate-reductase (pkdhfr) gene was sequenced from 449 P. knowlesi malaria cases from Sabah (Malaysian Borneo) and genotypes evaluated for association with clinical and epidemiological factors. Homology modelling using the pvdhfr template was used to assess the effect of pkdhfr mutations on the pyrimethamine binding pocket.

Results

Fourteen non-synonymous mutations were detected, with the most common being at codon T91P (10.2%) and R34L (10.0%), resulting in 21 different genotypes, including the wild-type, 14 single mutants, and six double mutants. One third of the P. knowlesi infections were with pkdhfr mutants; 145 (32%) patients had single mutants and 14 (3%) had double-mutants. In contrast, among the 47 P. falciparum isolates sequenced, three pfdhfr genotypes were found, with the double mutant 108N+59R being fixed and the triple mutants 108N+59R+51I and 108N+59R+164L occurring with frequencies of 4% and 8%, respectively. Two non-random spatio-temporal clusters were identified with pkdhfr genotypes. There was no association between pkdhfr mutations and hyperparasitaemia or malaria severity, both hypothesized to be indicators of H-H transmission. The orthologous loci associated with resistance in P. falciparum were not mutated in pkdhfr. Subsequent homology modelling of pkdhfr revealed gene loci 13, 53, 120, and 173 as being critical for pyrimethamine binding, however, there were no mutations at these sites among the 449 P. knowlesi isolates.

Conclusion

Although moderate diversity was observed in pkdhfr in Sabah, there was no evidence this reflected selective antifolate drug pressure in humans.     

Genome‐wide mosaicism in divergence between zoonotic malaria parasite subpopulations with separate sympatric transmission cycles

Plasmodium knowlesi is a significant cause of human malaria transmitted as a zoonosis from macaque reservoir hosts in South‐East Asia. Microsatellite genotyping has indicated that human infections in Malaysian Borneo are an admixture of two highly divergent sympatric parasite subpopulations that are, respectively, associated with long‐tailed macaques (Cluster 1) and pig‐tailed macaques (Cluster 2). Whole‐genome sequences of clinical isolates subsequently confirmed the separate clusters, although fewer of the less common Cluster 2 type were sequenced. Here, to analyse population structure and genomic divergence in subpopulation samples of comparable depth, genome sequences were generated from 21 new clinical infections identified as Cluster 2 by microsatellite analysis, yielding a cumulative sample size for this subpopulation similar to that for Cluster 1. Profound heterogeneity in the level of intercluster divergence was distributed across the genome, with long contiguous chromosomal blocks having high or low divergence. Different mitochondrial genome clades were associated with the two major subpopulations, but limited exchange of haplotypes from one to the other was evident, as was also the case for the maternally inherited apicoplast genome. These findings indicate deep divergence of the two sympatric P. knowlesi subpopulations, with introgression likely to have occurred recently. There is no evidence yet of specific adaptation at any introgressed locus, but the recombinant mosaic types offer enhanced diversity on which selection may operate in a currently changing landscape and human environment. Loci responsible for maintaining genetic isolation of the sympatric subpopulations need to be identified in the chromosomal regions showing fixed differences.     

A microarray platform and novel SNP calling algorithm to evaluate Plasmodium falciparum field samples of low DNA quantity

Background

Analysis of single nucleotide polymorphisms (SNPs) derived from whole-genome studies allows for rapid evaluation of genome-wide diversity, and genomic epidemiology studies of Plasmodium falciparum provide insights into parasite population structure, gene flow, drug resistance and vaccine development. In areas with adequate cold chain facilities, large volumes of leukocyte-depleted patient blood can be frozen for use in parasite genomic analyses. In more remote endemic areas smaller volumes of infected blood are taken by finger prick, and dried and stored on filter paper. These dried blood spots do not generally yield enough concentrated parasite DNA for whole-genome sequencing.

Results

A DNA microarray was designed for use on field samples to type a genome-wide set of SNPs which prior sequencing had shown to be variable in Africa, Southeast Asia, and Papua New Guinea. An algorithm was designed to call SNPs in samples with low parasite DNA. With this new algorithm SNP-calling accuracy of 98% was measured by hybridizing purified DNA from malaria lab strains and comparing calls with SNPs called from full genome sequences. An average accuracy of >98% was likewise obtained for DNA extracted from malaria field samples collected in studies in Southeast Asia, with an average call rate of > 82%.

Conclusion

This new high-density microarray provided high quality SNP calls from a wide range of parasite DNA quantities, and represents a robust tool for genome-wide analysis of malaria parasites in diverse settings.