Genetic diversity in the C-terminus of merozoite surface protein 1 among Plasmodium knowlesi isolates from Selangor and Sabah Borneo, Malaysia

Declining Genetic Polymorphisms of the C-terminus Merozoite Surface Protein-1 Amidst Increased Plasmodium knowlesi Transmission in Thailand

Recent reports from Thailand reveal a substantial surge in Plasmodium knowlesi cases over the past decade, with a more than eightfold increase in incidence by 2023 compared to 2018. This study investigates temporal changes in genetic polymorphism associated with the escalating transmission of P. knowlesi malaria in Thailand over the past two decades. Twenty-five P. knowlesi samples collected in 2018-2023 were sequenced for the 42-kDa region of pkmsp1 and compared with 24 samples collected in 2000-2009, focusing on nucleotide diversity, natural selection, recombination rate, and population differentiation. Seven unique haplotypes were identified in recent samples, compared to 15 in earlier samples. Nucleotide and haplotype diversities were lower in recent samples (π = 0.016, Hd = 0.817) than in earlier samples (π = 0.018, Hd = 0.942). Significantly higher synonymous substitution rates were observed in both sample sets (dS - dN = 2.77 and 2.43, p < 0.05), indicating purifying selection and reduced genetic diversity over time. Additionally, 8 out of 17 mutation points were located on B-cell epitopes, suggesting an adaptive response by the parasites to evade immune recognition. Population differentiation analysis using the fixation index (Fst) revealed high genetic differentiation between parasite populations in central and southern Thailand or Malaysia. Conversely, the relatively lower Fst value between southern Thailand and Malaysia suggests a closer genetic relationship, possibly reflecting historical gene flow. In conclusion, our findings highlight a decline in genetic diversity and evidence of purifying selection associated with the recently increased incidence of P. knowlesi malaria in Thailand. The minor genetic differentiation between P. knowlesi populations from southern Thailand and Malaysia suggests a shared recent ancestry of these parasites and underscores the need for coordinated efforts between the two countries for the elimination of P. knowlesi.

Genetic diversity and in silico analysis of Plasmodium knowlesi Serine Repeat Antigen (SERA) 3 antigen 2 in Malaysia

Plasmodium knowlesi is the leading cause of malaria in Malaysia. Serine Repeat Antigens (SERAs) have an essential role in the parasite life cycle. However, genetic characterization on P. knowlesi SERA3 Ag2 (PkSERA3 Ag2) is lacking. In the present study, nucleotide diversity, natural selection, and haplotypes of PkSERA3 Ag2 in clinical samples from Peninsular Malaysia and Malaysian Borneo were investigated. A total of 50 P. knowlesi clinical samples were collected from Peninsular Malaysia and Malaysian Borneo. The PkSERA3 Ag2 gene was amplified using PCR, and subsequently cloned and sequenced. Genetic diversity, haplotype, natural selection as well as genetic structure and differentiation of PkSERA3 Ag2 were analysed. In addition, in silico analyses were performed to identify repeat motifs, B-cell epitopes, and antigenicity indices of the protein. Analysis of 114 PkSERA3 Ag2 sequences revealed high nucleotide diversity of the gene in Malaysia. A codon-based Z-test indicated that the gene underwent purifying selection. Haplotype and population structure analyses identified two distinct PkSERA3 Ag2 clusters (K = 2, ΔK = 721.14) but no clear genetic distinction between PkSERA3 Ag2 from Peninsular Malaysia and Malaysian Borneo. FST index indicated moderate differentiation of the gene. In silico analyses revealed unique repeat motifs among PkSERA3 Ag2 isolates. Moreover, the amino acid sequence of PkSERA3 Ag2 exhibited potential B-cell epitopes and possessed high antigenicity indices. These findings enhance the understanding of PkSERA3 Ag2 gene as well as its antigenic properties. Further validation is necessary to ascertain the utility of PkSERA3 Ag2 as a serological marker for P. knowlesi infection.

Genetic Diversity and Population Genetic Structure Analysis of Plasmodium knowlesi Thrombospondin-Related Apical Merozoite Protein (TRAMP) in Clinical Samples

The simian malaria parasite Plasmodium knowlesi causes a high number of zoonotic infections in Malaysia. The thrombospondin-related apical merozoite protein (TRAMP) is an essential ligand for binding to the erythrocyte cell surface, whereby it facilitates the invasion. This study is the first attempt to determine the genetic diversity, phylogeography, natural selection and population structure from 97 full-length PkTRAMP gene sequences originating from Malaysia. We found low levels of nucleotide diversity (π~0.0065) for the full-length gene despite samples originating from geographically separated regions (i.e., Peninsular Malaysia and Malaysian Borneo). The rate of synonymous substitutions was significantly higher than that of non-synonymous substitutions, indicating a purifying selection for the full-length gene within the clinical samples. The population genetic analysis revealed that the parasite population is undergoing a significant population expansion. The analysis of the amino acid sequence alignment of 97 PkTRAMP sequences identified 15 haplotypes, of which a major shared haplotype was noted Hap 1 (n = 68, Sarawak; n = 34, Sabah; n = 12, Peninsular Malaysia; n = 22). The phylogenetic analysis using DNA sequences identified two clusters that separated due to geographical distance and three mixed clusters with samples from both Peninsular Malaysia and Malaysian Borneo. Population structure analyses indicated two distinct sub-populations (K = 2). Our findings point to the potential for independent parasite evolution, which could make zoonotic malaria control and elimination even more challenging.

Plasmodium knowlesi: the game changer for malaria eradication

Plasmodium knowlesi is a zoonotic malaria parasite that has gained increasing medical interest over the past two decades. This zoonotic parasitic infection is prevalent in Southeast Asia and causes many cases with fulminant pathology. Despite several biogeographical restrictions that limit its distribution, knowlesi malaria cases have been reported in different parts of the world due to travelling and tourism activities. Here, breakthroughs and key information generated from recent (over the past five years, but not limited to) studies conducted on P. knowlesi were reviewed, and the knowledge gap in various research aspects that need to be filled was discussed. Besides, challenges and strategies required to control and eradicate human malaria with this emerging and potentially fatal zoonosis were described.

Genetic population of Plasmodium knowlesi during pre-malaria elimination in Thailand

Background

Thailand is committed to eliminating malaria by 2024. From 2013 to 2020, the total number of malaria cases have decreased, from 37,741 to 4474 (an 88.1% reduction). However, infections with Plasmodium knowlesi, a monkey malarial pathogen that can also infect humans, have been increasingly observed. This study focused on the molecular analysis of P. knowlesi parasites causing malaria in Thailand.

Methods

Under Thailand’s integrated Drug Efficacy Surveillance (iDES), which includes drug-resistance monitoring as part of routine case-based surveillance and responses, specimens were collected from malaria patients (n = 966) between 2018 and 2020. Thirty-one mono P. knowlesi infections (3.1%), most of which were from eastern and southern Thailand, were observed and confirmed by nested PCR assay and DNA sequencing. To evaluate whether these pathogens were from different lineages, cluster analysis based on seven microsatellite genotyping markers and the merozoite surface protein 1 (pkmsp1) gene was carried out. The P. knowlesi pyrimethamine resistance gene dihydrofolate reductase (pkdhfr) was sequenced and homology modelling was constructed.

Results

The results of analysing the seven microsatellite markers and pkmsp1 sequence demonstrated that P. knowlesi parasites from eastern Thailand were of the same lineage as those isolated in Cambodia, while the parasites causing malaria in southern Thailand were the same lineage as those isolated from Malaysia. The sequencing results for the pkdhfr genes indicated the presence of two mutations, Arg34Leu and a deletion at position 105. On analysis with homology modelling, the two mutations were not associated with anti-malarial drug resistance.

Conclusions

This report compared the genetic populations of P. knowlesi parasites in Thailand from 2018 to 2020 and have shown similar lineages as those isolated in Cambodia and Malaysia of P. knowlesi infection in Thailand and demonstrated that the P. knowlesi parasites were of the same lineages as those isolated in Cambodia and Malaysia. The parasites were also shown to be sensitive to pyrimethamine.

The vectors of Plasmodium knowlesi and other simian malarias Southeast Asia: challenges in malaria elimination

Plasmodium knowlesi, a simian malaria parasite of great public health concern has been reported from most countries in Southeast Asia and exported to various countries around the world. Currently P. knowlesi is the predominant species infecting humans in Malaysia. Besides this species, other simian malaria parasites such as P. cynomolgi and P. inui are also infecting humans in the region. The vectors of P. knowlesi and other Asian simian malarias belong to the Leucosphyrus Group of Anopheles mosquitoes which are generally forest dwelling species. Continual deforestation has resulted in these species moving into forest fringes, farms, plantations and human settlements along with their macaque hosts. Limited studies have shown that mosquito vectors are attracted to both humans and macaque hosts, preferring to bite outdoors and in the early part of the night. We here review the current status of simian malaria vectors and their parasites, knowledge of vector competence from experimental infections and discuss possible vector control measures. The challenges encountered in simian malaria elimination are also discussed. We highlight key knowledge gaps on vector distribution and ecology that may impede effective control strategies.

Knowlesi malaria: Human risk factors, clinical spectrum, and pathophysiology

Plasmodium knowlesi is endemic across Southeast Asia, and is the commonest cause of zoonotic malaria. The spectrum of clinical disease from P. knowlesi infection ranges from asymptomatic infection, through to severe malaria and death. Over 90% of clinical disease occurs in adults, mostly living in forest edge areas undergoing intensive land use change. With a 24-h asexual life cycle in humans, high parasite counts are possible, but most clinical cases of knowlesi malaria are uncomplicated with low parasitaemia. In co-endemic areas, median parasitaemia in knowlesi malaria is lower than that seen in vivax and falciparum malaria, suggesting a lower fever threshold. Severe malaria occurs in 6-9% of symptomatic adults. Manifestations of severe malaria from P. knowlesi are similar to those seen with falciparum malaria, with the notable absence of coma. Age, parasitaemia, cardiovascular comorbidities and delayed diagnosis are risk factors for severe disease and death, which are only seen in adults. Thrombocytopenia is near-universal in adults, likely related to platelet-red cell binding and clearance. Mechanisms underlying the microvascular sludging seen in fatal disease in non-natural primate hosts and the microvascular accumulation of parasites in fatal human disease are not clear. Marked reductions in deformability of both infected and uninfected red blood cells are associated with disease severity in both humans and other non-natural primate hosts, likely contributing to impaired microvascular perfusion and organ dysfunction. Endothelial activation, endothelial dysfunction, glycocalyx degradation and haemolysis are also associated with, and likely contribute to, severe disease and organ dysfunction, particularly acute kidney injury.

Identification of Plasmodium knowlesi Merozoite Surface Protein-119 (PkMSP-119 ) novel binding peptides from a phage display library

OBJECTIVE

Plasmodium knowlesi, the fifth human malaria parasite, has caused mortality in humans. We aimed to identify P. knowlesi novel binding peptides through a random linear dodecapeptide phage display targeting the 19-kDa fragment of Merozoite Surface Protein-1 protein.

METHODS

rPkMSP-1₁₉ protein was heterologously expressed using Expresso^® Solubility and Expression Screening System and competent E. cloni^® 10G cells according to protocol. Three rounds of biopanning were performed on purified rPkMSP-1₁₉ to identify binding peptides towards rPkMSP-1₁₉ using Ph.D.™-12 random phage display library. Binding sites of the identified peptides to PkMSP-1₁₉ were in silico predicted using the CABS-dock web server.

RESULTS

Four phage peptide variants that bound to PkMSP-1₁₉ were identified after three rounds of biopanning, namely Pkd1, Pkd2, Pkd3 and Pkd4. The sequences of both Pkd1 and Pkd2 consist of a large number of histidine residues. Pkd1 showed positive binding signal with 6.1× vs. BSA control. Docking results showed that Pkd1 and Pkd2 were ideal binding peptides for PkMSP-1₁₉ .

CONCLUSION

We identified two novel binding peptides of PkMSP-1₁₉ , Pkd1 (HFPFHHHKLRAH) and Pkd2 (HPMHMLHKRQHG), through phage display. They provide a valuable starting point for the development of novel therapeutics.

Plasmodium knowlesi clinical isolates from Malaysia show extensive diversity and strong differential selection pressure at the merozoite surface protein 7D (MSP7D)

BACKGROUND

The high proportion of human cases due to the simian malaria parasite Plasmodium knowlesi in Malaysia is a cause of concern, as they can be severe and even fatal. Merozoite surface protein 7 (MSP7) is a multigene family which forms a non-covalent complex with MSP-1 prior to receptor-ligand recognition in Plasmodium falciparum and thus an important antigen for vaccine development. However, no study has been done in any of the ortholog family members in P. knowlesi from clinical samples. This study investigates the level of polymorphism, haplotypes, and natural selection acting at the pkmsp-7D gene in clinical samples from Malaysia.

METHODS

Thirty-six full-length pkmsp7D gene sequences (along with the reference H-strain: PKNH_1266000) obtained from clinical isolates of Malaysia, which were orthologous to pvmsp7H (PVX_082680) were downloaded from public databases. Population genetic, evolutionary and phylogenetic analyses were performed to determine the level of genetic diversity, polymorphism, recombination and natural selection.

RESULTS

Analysis of 36 full-length pkmsp7D sequences identified 147 SNPs (91 non-synonymous and 56 synonymous substitutions). Nucleotide diversity across the full-length gene was higher than its ortholog in Plasmodium vivax (msp7H). Region-wise analysis of the gene indicated that the nucleotide diversity at the central region was very high (π = 0.14) compared to the 5' and 3' regions. Most hyper-variable SNPs were detected at the central domain. Multiple test for natural selection indicated the central region was under strong positive natural selection however, the 5' and 3' regions were under negative/purifying selection. Evidence of intragenic recombination were detected at the central region of the gene. Phylogenetic analysis using full-length msp7D genes indicated there was no geographical clustering of parasite population.

CONCLUSIONS

High genetic diversity with hyper-variable SNPs and strong evidence of positive natural selection at the central region of MSP7D indicated exposure of the region to host immune pressure. Negative selection at the 5' and the 3' regions of MSP7D might be because of functional constraints at the unexposed regions during the merozoite invasion process of P. knowlesi. No evidence of geographical clustering among the clinical isolates from Malaysia indicated uniform selection pressure in all populations. These findings highlight the further evaluation of the regions and functional characterization of the protein as a potential blood stage vaccine candidate for P. knowlesi.

The Plasmodium knowlesi Pk41 surface protein diversity, natural selection, sub population and geographical clustering: a 6-cysteine protein family member

Introduction

The zoonotic malaria parasite Plasmodium knowlesi has currently become the most dominant form of infection in humans in Malaysia and is an emerging infectious disease in most Southeast Asian countries. The P41 is a merozoite surface protein belonging to the 6-cysteine family and is a well-characterized vaccine candidate in P. vivax and P. falciparum; however, no study has been done in the orthologous gene of P. knowlesi. This study investigates the level of polymorphism, haplotypes and natural selection of pk41 genes in clinical isolates from Malaysia.

Method

Thirty-five full-length pk41 sequences from clinical isolates of Malaysia along with four laboratory lines (along with H-strain) were downloaded from public databases. For comparative analysis between species, orthologous P41 genes from P. falciparum, P. vivax, P. coatneyi and P. cynomolgi were also downloaded. Genetic diversity, polymorphism, haplotype and natural selection were determined using DnaSP 5.10 software. Phylogenetic relationships between Pk41 genes were determined using MEGA 5.0 software.

Results

Analysis of 39 full-length pk41 sequences along with the H-strain identified 36 SNPs (20 non-synonymous and 16 synonymous substitutions) resulting in 31 haplotypes. Nucleotide diversity across the full-length gene was low and was similar to its ortholog in P. vivax; pv41. Domain-wise amino acid analysis of the two s48/45 domains indicated low level of polymorphisms for both the domains, and the glutamic acid rich region had extensive size variations. In the central domain, upstream to the glutamate rich region, a unique two to six (K-E)_n repeat region was identified within the clinical isolates. Overall, the pk41 genes were indicative of negative/purifying selection due to functional constraints. Domain-wise analysis of the s48/45 domains also indicated purifying selection. However, analysis of Tajima’s D across the genes identified non-synonymous SNPs in the s48/45 domain II with high positive values indicating possible epitope binding regions. All the 6-cysteine residues within the s48/45 domains were conserved within the clinical isolates indicating functional conservation of these regions. Phylogenetic analysis of full-length pk41 genes indicated geographical clustering and identified three subpopulations of P. knowlesi; one originating in the laboratory lines and two originating from Sarawak, Malaysian Borneo.

Conclusion

This is the first study to report on the polymorphism and natural selection of pk41 genes from clinical isolates of Malaysia. The results reveal that there is low level of polymorphism in both s48/45 domains, indicating that this antigen could be a potential vaccine target. However, genetic and molecular immunology studies involving higher number of samples from various parts of Malaysia would be necessary to validate this antigen’s candidacy as a vaccine target for P. knowlesi.

Genetic polymorphism and natural selection in the C-terminal 42 kDa region of merozoite surface protein-1 (MSP-1) among Plasmodium knowlesi samples from Malaysia

Background

The merozoite surface protein-1 (MSP-1) gene encodes for a leading malaria vaccine candidate antigen. However, its extensive polymorphic nature represents a major obstacle to the development of a protective vaccine. Previously, a pilot study was carried out to explore the sequence variation of the C-terminal 42 kDa fragment within P. knowlesi MSP-1 gene (PkMSP-1₄₂) based on 12 clinical samples; however, further study on an adequate sample size is vital in estimating the genetic diversity of the parasite population.

Methods

In the present study, we included a larger sample size of P. knowlesi (83 samples) covering eight states of Malaysia to determine the genetic polymorphism, natural selection and haplotype groups of the gene fragment coding PkMSP-1₄₂. The region flanking PkMSP-1₄₂ was amplified by PCR and directly sequenced. Genetic diversity, haplotype diversity, population genetic differentiation and natural selection were determined in order to study the polymorphic characteristic of PkMSP-1₄₂.

Results

A high level of genetic diversity (Hd = 0.970 ± 0.007; л = 0.01079 ± 0.00033) was observed among the 83 P. knowlesi samples, confirming the extensive genetic polymorphism exhibited among the P. knowlesi population found in Malaysia. A total of 18 distinct haplotypes with 17 amino acid changes were identified, whereby 15 were new haplotypes. High population differentiation values were observed within samples from Peninsular Malaysia and Malaysian Borneo. The 42 kDa fragments of P. knowlesi from Malaysian Borneo were found to be acting on balancing selection whilst purifying selection was suggested to act on isolates from Peninsular Malaysia. The separation of PkMSP-1₄₂ haplotypes into two main groups based on geographical separation has further supported the existence of two distinct P. knowlesi lineages.

Conclusions

A high level of genetic diversity was observed among PkMSP-1₄₂ in Malaysia, whereby most of the polymorphisms were found within the 33 kDa region. Taken together, these data will be useful in order to understand the nature of P. knowlesi population in Malaysia as well as the design and development of a MSP-1₄₂ based knowlesi malaria vaccine.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3234-5) contains supplementary material, which is available to authorized users.

Within-population genetic diversity and population structure of Plasmodium knowlesi merozoite surface protein 1 gene from geographically distinct regions of Malaysia and Thailand

BACKGROUND

The C-terminal 42 kDa domain of Plasmodium knowlesi merozoite surface protein 1 (PkMSP1) is a potential asexual blood-stage vaccine candidate, however, only a limited number of clinical isolates have been analysed from Malaysia and no inter-country comparative diversity study has been conducted. In the present study, nucleotide diversity, haplotypes and natural selection levels of pkmsp1 in clinical samples from geographically distinct regions of Malaysia and Thailand were investigated. The overall population structure of the parasite from the region was determined.

METHODS

Eleven full-length pkmsp1 sequences obtained from clinical isolates of Malaysia along with the H-strain were downloaded from the database for domain wise characterization of pkmsp1 gene. Additionally, 76 pkmsp-1₄₂ sequences from Thailand and Malaysia were downloaded from the database for intra and inter-population analysis. DnaSP 5.10 and MEGA 5.0 software were used to determine genetic diversity, polymorphism, haplotypes and natural selection. Genealogical relationships were determined using haplotype network tree in NETWORK software v5.0. Population genetic differentiation index (F_ST) of parasites were analysed using Arlequin v3.5.

RESULTS

Sequence analysis of 11 full-length pkmsp1 sequences along with the H-strain identified 477 (8.4%) polymorphic sites, of which 107 were singleton sites. The overall diversity observed in the full-length genes were high in comparison to its ortholog pvmsp1 and the 4 variable domains showed extensive size variations. The nucleotide diversity was low towards the pkmsp1-42 compared to the conserved domains. The 19 kDa domain was less diverse and completely conserved among isolates from Malaysian Borneo. The nucleotide diversity of isolates from Peninsular Malaysia and Thailand were higher than Malaysian Borneo. Network analysis of pkmsp1-42 haplotypes showed geographical clustering of the isolates from Malaysian Borneo and grouping of isolates from Peninsular Malaysia and Thailand. Population differentiation analysis indicated high F_ST values between parasite populations originating from Malaysian Borneo, Peninsular Malaysia and Thailand attributing to geographical distance. Moderate genetic differentiation was observed for parasite populations from Thailand and Peninsular Malaysia. Evidence of population expansion and purifying selection were observed in all conserved domains with strongest selection within the pkmsp1-42 domain.

CONCLUSIONS

This study is the first to report on inter country genetic diversity and population structure of P. knowlesi based on msp1. Strong evidence of negative selection was observed in the 42 kDa domain, indicating functional constrains. Geographical clustering of P. knowlesi and moderate to high genetic differentiation values between populations identified in this study highlights the importance of further evaluation using larger number of clinical samples from Southeast Asian countries.

Plasmodium knowlesi malaria: current research perspectives

Originally known to cause simian malaria, Plasmodium knowlesi is now known as the fifth human malaria species. Since the publishing of a report that largely focused on human knowlesi cases in Sarawak in 2004, many more human cases have been reported in nearly all of the countries in Southeast Asia and in travelers returning from these countries. The zoonotic nature of this infection hinders malaria elimination efforts. In order to grasp the current perspective of knowlesi malaria, this literature review explores the different aspects of the disease including risk factors, diagnosis, treatment, and molecular and functional studies. Current studies do not provide sufficient data for an effective control program. Therefore, future direction for knowlesi research is highlighted here with a final aim of controlling, if not eliminating, the parasite.

Diversity and natural selection on the thrombospondin-related adhesive protein (TRAP) gene of Plasmodium knowlesi in Malaysia

Background

Plasmodium knowlesi a parasite of the macaques is currently the most common cause of human malaria in Malaysia. The thrombospondin-related adhesive protein (TRAP) gene is pre-erythrocytic stage antigen. It is a well-characterized vaccine candidate in Plasmodium vivax and Plasmodium falciparum, however, no study has been done in the orthologous gene of P. knowlesi. This study investigates nucleotide diversity, haplotypes, natural selection and population differentiation of full-length pktrap genes in clinical samples from Malaysia.

Methods

Forty full-length pktrap sequences from clinical isolates of Malaysia along with the reference H-strain were downloaded from published databases. Genetic diversity, polymorphism, haplotype and natural selection were determined using DnaSP 5.10 software. McDonald–Kreitman test was conducted using P. vivax and Plasmodium coatneyi as ortholog sequence in DnaSP 5.10 software. Population genetic differentiation index (F_ST) of parasite populations was determined using Arlequin v3.5. Phylogenetic relationships between trap ortholog genes were determined using MEGA 5.0 software.

Results

Comparison of 40 full-length pktrap sequences along with the H-strain identified 74 SNPs (53 non-synonymous and 21 synonymous substitutions) resulting in 29 haplotypes. Analysis of the full-length gene showed that the nucleotide diversity was lower compared to its nearest ortholog pvtrap. Domain-wise analysis indicated that the proline/asparagine rich region had higher nucleotide diversity compared to the von Willebrand factor domain and the thrombospondin-type-1 domain. McDonald–Kreitman test identified that the ratio of the number of nonsynonymous to synonymous polymorphic sites within P. knowlesi was significantly higher than that of the number of nonsynonymous to synonymous fixed sites between P. knowlesi and P. vivax. The von Willebrand factor domain also indicated balancing selection using MK test, however, it did not give significant results when tested with P. coatneyi as an outgroup. Phylogenetic analysis of full-length genes identified three distinct sub-clusters of P. knowlesi, one originating from Peninsular Malaysia and two originating from Malaysian Borneo. High population differentiation values was observed within samples from Peninsular Malaysia and Malaysian Borneo.

Conclusions

This study is the first to report on the genetic diversity and natural selection of full-length pktrap. Low level of genetic diversity was found across the full-length gene of pktrap. Balancing selection of the von Willebrand factor domain indicated that TRAP could be a target in inducing immune response against P. knowlesi infections. However, higher number of samples would be necessary to further confirm the findings.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2423-1) contains supplementary material, which is available to authorized users.