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.     

Plasmodium knowlesi Genome Sequences from Clinical Isolates Reveal Extensive Genomic Dimorphism

Plasmodium knowlesi is a newly described zoonosis that causes malaria in the human population that can be severe and fatal. The study of P. knowlesi parasites from human clinical isolates is relatively new and, in order to obtain maximum information from patient sample collections, we explored the possibility of generating P. knowlesi genome sequences from archived clinical isolates. Our patient sample collection consisted of frozen whole blood samples that contained excessive human DNA contamination and, in that form, were not suitable for parasite genome sequencing. We developed a method to reduce the amount of human DNA in the thawed blood samples in preparation for high throughput parasite genome sequencing using Illumina HiSeq and MiSeq sequencing platforms. Seven of fifteen samples processed had sufficiently pure P. knowlesi DNA for whole genome sequencing. The reads were mapped to the P. knowlesi H strain reference genome and an average mapping of 90% was obtained. Genes with low coverage were removed leaving 4623 genes for subsequent analyses. Previously we identified a DNA sequence dimorphism on a small fragment of the P. knowlesi normocyte binding protein xa gene on chromosome 14. We used the genome data to assemble full-length Pknbpxa sequences and discovered that the dimorphism extended along the gene. An in-house algorithm was developed to detect SNP sites co-associating with the dimorphism. More than half of the P. knowlesi genome was dimorphic, involving genes on all chromosomes and suggesting that two distinct types of P. knowlesi infect the human population in Sarawak, Malaysian Borneo. We use P. knowlesi clinical samples to demonstrate that Plasmodium DNA from archived patient samples can produce high quality genome data. We show that analyses, of even small numbers of difficult clinical malaria isolates, can generate comprehensive genomic information that will improve our understanding of malaria parasite diversity and pathobiology.     

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.     

Plasmodium knowlesi from archival blood films: Further evidence that human infections are widely distributed and not newly emergent in Malaysian Borneo

Human infections with Plasmodium knowlesi have been misdiagnosed by microscopy as Plasmodium malariae due to their morphological similarities. Although microscopy-identified P. malariae cases have been reported in the state of Sarawak (Malaysian Borno) as early as 1952, recent epidemiological studies suggest the absence of indigenous P. malariae infections. The present study aimed to determine the past incidence and distribution of P. knowlesi infections in the state of Sarawak based on archival blood films from patients diagnosed by microscopy as having P. malariae infections. Nested PCR assays were used to identify Plasmodium species in DNA extracted from 47 thick blood films collected in 1996 from patients in seven different divisions throughout the state of Sarawak. Plasmodium knowlesi DNA was detected in 35 (97.2%) of 36 blood films that were positive for Plasmodium DNA, with patients originating from all seven divisions. Only one sample was positive for P. malariae DNA. This study provides further evidence of the widespread distribution of human infections with P. knowlesi in Sarawak and its past occurrence. Taken together with data from previous studies, our findings suggest that P. knowlesi malaria is not a newly emergent disease in humans.     

The Ratio of Monocytes to Lymphocytes in Peripheral Blood Correlates with Increased Susceptibility to Clinical Malaria in Kenyan Children

Background

Plasmodium falciparum malaria remains a major cause of illness and death in sub-Saharan Africa. Young children bear the brunt of the disease and though older children and adults suffer relatively fewer clinical attacks, they remain susceptible to asymptomatic P. falciparum infection. A better understanding of the host factors associated with immunity to clinical malaria and the ability to sustain asymptomatic P. falciparum infection will aid the development of improved strategies for disease prevention.

Methods and Findings

Here we investigate whether full differential blood counts can predict susceptibility to clinical malaria among Kenyan children sampled at five annual cross-sectional surveys. We find that the ratio of monocytes to lymphocytes, measured in peripheral blood at the time of survey, directly correlates with risk of clinical malaria during follow-up. This association is evident among children with asymptomatic P. falciparum infection at the time the cell counts are measured (Hazard ratio (HR)  =  2.7 (95% CI 1.42, 5.01, P  =  0.002) but not in those without detectable parasitaemia (HR  =  1.0 (95% CI 0.74, 1.42, P  =  0.9).

Conclusions

We propose that the monocyte to lymphocyte ratio, which is easily derived from routine full differential blood counts, reflects an individual''s capacity to mount an effective immune response to P. falciparum infection.     

Genotyping of the Duffy Blood Group among Plasmodium knowlesi-Infected Patients in Malaysia

The Duffy blood group is of major interest in clinical medicine as it plays an important role in Plasmodium knowlesi and Plasmodium vivax infection. In the present study, the distribution of Duffy blood group genotypes and allelic frequencies among P. knowlesi infected patients as well as healthy individuals in Peninsular Malaysia were determined. The blood group of 60 healthy blood donors and 51 P. knowlesi malaria patients were genotyped using allele specific polymerase chain reaction (ASP-PCR). The data was analyzed using Fisher''s exact test in order to assess the significance of the variables. Our results show a high proportion of the FY*A/FY*A genotype (>85% for both groups) and a high frequency of the FY*A allele (>90% for both groups). The FY*A/FY*A genotype was the most predominant genotype in both infected and healthy blood samples. The genotype frequency did not differ significantly between the donor blood and the malaria patient groups. Also, there was no significant correlation between susceptibility to P. knowlesi infection with any Duffy blood genotype.     

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.     

Effective Preparation of Plasmodium vivax Field Isolates for High-Throughput Whole Genome Sequencing

Whole genome sequencing (WGS) of Plasmodium vivax is problematic due to the reliance on clinical isolates which are generally low in parasitaemia and sample volume. Furthermore, clinical isolates contain a significant contaminating background of host DNA which confounds efforts to map short read sequence of the target P. vivax DNA. Here, we discuss a methodology to significantly improve the success of P. vivax WGS on natural (non-adapted) patient isolates. Using 37 patient isolates from Indonesia, Thailand, and travellers, we assessed the application of CF11-based white blood cell filtration alone and in combination with short term ex vivo schizont maturation. Although CF11 filtration reduced human DNA contamination in 8 Indonesian isolates tested, additional short-term culture increased the P. vivax DNA yield from a median of 0.15 to 6.2 ng µl⁻¹ packed red blood cells (pRBCs) (p = 0.001) and reduced the human DNA percentage from a median of 33.9% to 6.22% (p = 0.008). Furthermore, post-CF11 and culture samples from Thailand gave a median P. vivax DNA yield of 2.34 ng µl⁻¹ pRBCs, and 2.65% human DNA. In 22 P. vivax patient isolates prepared with the 2-step method, we demonstrate high depth (median 654X coverage) and breadth (≥89%) of coverage on the Illumina GAII and HiSeq platforms. In contrast to the A+T-rich P. falciparum genome, negligible bias was observed in coverage depth between coding and non-coding regions of the P. vivax genome. This uniform coverage will greatly facilitate the detection of SNPs and copy number variants across the genome, enabling unbiased exploration of the natural diversity in P. vivax populations.     

Plasma Superoxide Dismutase-1 as a Surrogate Marker of Vivax Malaria Severity

Background

Severe outcomes have been described for both Plasmodium falciparum and P. vivax infections. The identification of sensitive and reliable markers of disease severity is fundamental to improving patient care. An intense pro-inflammatory response with oxidative stress and production of reactive oxygen species is present in malaria. Inflammatory cytokines such as tumor necrosis factor-alpha (TNF-alpha) and antioxidant agents such as superoxide dismutase-1 (SOD-1) are likely candidate biomarkers for disease severity. Here we tested whether plasma levels of SOD-1 could serve as a biomarker of severe vivax malaria.

Methodology/Principal Findings

Plasma samples were obtained from residents of the Brazilian Amazon with a high risk for P. vivax transmission. Malaria diagnosis was made by both microscopy and nested PCR. A total of 219 individuals were enrolled: non-infected volunteers (n = 90) and individuals with vivax malaria: asymptomatic (n = 60), mild (n = 50) and severe infection (n = 19). SOD-1 was directly associated with parasitaemia, plasma creatinine and alanine amino-transaminase levels, while TNF-alpha correlated only with the later enzyme. The predictive power of SOD-1 and TNF-alpha levels was compared. SOD-1 protein levels were more effective at predicting vivax malaria severity than TNF-alpha. For discrimination of mild infection, elevated SOD-1 levels showed greater sensitivity than TNF-alpha (76% vs. 30% respectively; p<0.0001), with higher specificity (100% vs. 97%; p<0.0001). In predicting severe vivax malaria, SOD-1 levels exhibited higher sensitivity than TNF-alpha (80% vs. 56%, respectively; p<0.0001; likelihood ratio: 7.45 vs. 3.14; p<0.0001). Neither SOD-1 nor TNF-alpha could discriminate P. vivax infections from those caused by P. falciparum.

Conclusion

SOD-1 is a powerful predictor of disease severity in individuals with different clinical presentations of vivax malaria.     

Transforming Growth Factor Beta 2 and Heme Oxygenase 1 Genes Are Risk Factors for the Cerebral Malaria Syndrome in Angolan Children

Background

Cerebral malaria (CM) represents a severe outcome of the Plasmodium falciparum infection. Recent genetic studies have correlated human genes with severe malaria susceptibility, but there is little data on genetic variants that increase the risk of developing specific malaria clinical complications. Nevertheless, susceptibility to experimental CM in the mouse has been linked to host genes including Transforming Growth Factor Beta 2 (TGFB2) and Heme oxygenase-1 (HMOX1). Here, we tested whether those genes were governing the risk of progressing to CM in patients with severe malaria syndromes.

Methodology/Principal Findings

We report that the clinical outcome of P. falciparum infection in a cohort of Angolan children (n = 430) correlated with nine TGFB2 SNPs that modify the risk of progression to CM as compared to other severe forms of malaria. This genetic effect was explained by two haplotypes harboring the CM-associated SNPs (Pcorrec. = 0.035 and 0.036). In addition, one HMOX1 haplotype composed of five CM-associated SNPs increased the risk of developing the CM syndrome (Pcorrec. = 0.002) and was under-transmitted to children with uncomplicated malaria (P = 0.036). Notably, the HMOX1-associated haplotype conferred increased HMOX1 mRNA expression in peripheral blood cells of CM patients (P = 0.012).

Conclusions/Significance

These results represent the first report on CM genetic risk factors in Angolan children and suggest the novel hypothesis that genetic variants of the TGFB2 and HMOX1 genes may contribute to confer a specific risk of developing the CM syndrome in patients with severe P. falciparum malaria. This work may provide motivation for future studies aiming to replicate our findings in larger populations and to confirm a role for these genes in determining the clinical course of malaria.     

The Plasmodium knowlesi MAHRP2 ortholog localizes to structures connecting Sinton Mulligan's clefts in the infected erythrocyte

During development within the host erythrocyte malaria parasites generate nascent membranous structures which serve as a pathway for parasite protein transport to modify the host cell. The molecular basis of such membranous structures is not well understood, particularly for malaria parasites other than Plasmodium falciparum. To characterize the structural basis of protein trafficking in the Plasmodium knowlesi-infected erythrocyte, we identified a P. knowlesi ortholog of MAHRP2, a marker of the tether structure that connects membranous structures in the P. falciparum-infected erythrocyte. We show that PkMAHRP2 localizes on amorphous structures that connect Sinton Mulligan's clefts (SMC) to each other and to the erythrocyte membrane. Three dimensional reconstruction of the P. knowlesi-infected erythrocyte revealed that the SMC is a plate-like structure with swollen ends, reminiscent of the morphology of the Golgi apparatus. The PkMAHRP2-localized amorphous structures are possibly functionally equivalent to P. falciparum tether structure. These findings suggest a conservation in the ultrastructure of protein trafficking between P. falciparum and P. knowlesi.     

The Plasmodium vivax Merozoite Surface Protein 3β Sequence Reveals Contrasting Parasite Populations in Southern and Northwestern Thailand

Background

Malaria control efforts have a significant impact on the epidemiology and parasite population dynamics. In countries aiming for malaria elimination, malaria transmission may be restricted to limited transmission hot spots, where parasite populations may be isolated from each other and experience different selection forces. Here we aim to examine the Plasmodium vivax population divergence in geographically isolated transmission zones in Thailand.

Methodology

We employed the P. vivax merozoite surface protein 3β (PvMSP3β) as a molecular marker for characterizing P. vivax populations based on the extensive diversity of this gene in Southeast Asian parasite populations. To examine two parasite populations with different transmission levels in Thailand, we obtained 45 P. vivax isolates from Tak Province, northwestern Thailand, where the annual parasite incidence (API) was more than 2%, and 28 isolates from Yala and Narathiwat Provinces, southern Thailand, where the API was less than 0.02%. We sequenced the PvMSP3β gene and examined its genetic diversity and molecular evolution between the parasite populations.

Principal Findings

Of 58 isolates containing single PvMSP3β alleles, 31 sequence types were identified. The overall haplotype diversity was 0.77±0.06 and nucleotide diversity 0.0877±0.0054. The northwestern vivax malaria population exhibited extensive haplotype diversity (HD) of PvMSP3β (HD = 1.0). In contrast, the southern parasite population displayed a single PvMSP3β allele (HD = 0), suggesting a clonal population expansion. This result revealed that the extent of allelic diversity in P. vivax populations in Thailand varies among endemic areas.

Conclusion

Malaria parasite populations in a given region may vary significantly in genetic diversity, which may be the result of control and influenced by the magnitude of malaria transmission intensity. This is an issue that should be taken into account for the implementation of P. vivax control measures such as drug policy and vaccine development.     

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.     

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.     

Major Burden of Severe Anemia from Non-Falciparum Malaria Species in Southern Papua: A Hospital-Based Surveillance Study

Background

The burden of anemia attributable to non-falciparum malarias in regions with Plasmodium co-endemicity is poorly documented. We compared the hematological profile of patients with and without malaria in southern Papua, Indonesia.

Methods and Findings

Clinical and laboratory data were linked for all patients presenting to a referral hospital between April 2004 and December 2012. Data were available on patient demographics, malaria diagnosis, hemoglobin concentration, and clinical outcome, but other potential causes of anemia could not be identified reliably. Of 922,120 patient episodes (837,989 as outpatients and 84,131 as inpatients), a total of 219,845 (23.8%) were associated with a hemoglobin measurement, of whom 67,696 (30.8%) had malaria. Patients with P. malariae infection had the lowest hemoglobin concentration (n = 1,608, mean = 8.93 [95% CI 8.81–9.06]), followed by those with mixed species infections (n = 8,645, mean = 9.22 [95% CI 9.16–9.28]), P. falciparum (n = 37,554, mean = 9.47 [95% CI 9.44–9.50]), and P. vivax (n = 19,858, mean = 9.53 [95% CI 9.49–9.57]); p-value for all comparisons <0.001. Severe anemia (hemoglobin <5 g/dl) was present in 8,151 (3.7%) patients. Compared to patients without malaria, those with mixed Plasmodium infection were at greatest risk of severe anemia (adjusted odds ratio [AOR] 3.25 [95% CI 2.99–3.54]); AORs for severe anaemia associated with P. falciparum, P. vivax, and P. malariae were 2.11 (95% CI 2.00–2.23), 1.87 (95% CI 1.74–2.01), and 2.18 (95% CI 1.76–2.67), respectively, p<0.001. Overall, 12.2% (95% CI 11.2%–13.3%) of severe anemia was attributable to non-falciparum infections compared with 15.1% (95% CI 13.9%–16.3%) for P. falciparum monoinfections. Patients with severe anemia had an increased risk of death (AOR = 5.80 [95% CI 5.17–6.50]; p<0.001). Not all patients had a hemoglobin measurement, thus limitations of the study include the potential for selection bias, and possible residual confounding in multivariable analyses.

Conclusions

In Papua P. vivax is the dominant cause of severe anemia in early infancy, mixed P. vivax/P. falciparum infections are associated with a greater hematological impairment than either species alone, and in adulthood P. malariae, although rare, is associated with the lowest hemoglobin concentration. These findings highlight the public health importance of integrated genus-wide malaria control strategies in areas of Plasmodium co-endemicity. Please see later in the article for the Editors'' Summary