Search

Experimental studies on the biology of malaria parasites have mostly been based on laboratory-adapted lines, but there is limited understanding of how these may differ from parasites in natural infections. Loss-of-function mutants have previously been shown to emerge during culture of some Plasmodium falciparum clinical isolates in analyses focusing on single-genotype infections. The present study included a broader array of isolates, mostly representing multiple-genotype infections, which are more typical in areas where malaria is highly endemic. Genome sequence data from multiple time points over several months of culture adaptation of 28 West African isolates were analysed, including previously available sequences along with new genome sequences from additional isolates and time points. Some genetically complex isolates eventually became fixed over time to single surviving genotypes in culture, whereas others retained diversity, although proportions of genotypes varied over time. Drug resistance allele frequencies did not show overall directional changes, suggesting that resistance-associated costs are not the main causes of fitness differences among parasites in culture. Loss-of-function mutants emerged during culture in several of the multiple-genotype isolates, affecting genes (including AP2-HS, EPAC and SRPK1) for which loss-of-function mutants were previously seen to emerge in single-genotype isolates. Parasite clones were derived by limiting dilution from six of the isolates, and sequencing identified de novo variants not detected in the bulk isolate sequences. Interestingly, several of these were nonsense mutants and frameshifts disrupting the coding sequence of EPAC, the gene with the largest number of independent nonsense mutants previously identified in laboratory-adapted lines. Analysis of genomic identity by descent to explore relatedness among clones revealed co-occurring non-identical sibling parasites, illustrative of the natural genetic structure within endemic populations.

Malaria parasites break down host haemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter, pfcrt, located in the digestive vacuole membrane, confer CQ resistance in Plasmodium falciparum, and typically also affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on a putative amino acid transporter (pfaat1) variant S258L, which increased from 0% to 97% in frequency between 1984 and 2014 in parallel with the pfcrt1 K76T variant. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containing pfaat1 that is selected by CQ treatment. Gene editing demonstrated that pfaat1 S258L potentiates CQ resistance but at a cost of reduced fitness, while pfaat1 F313S, a common southeast Asian polymorphism, reduces CQ resistance while restoring fitness. Our analyses reveal hidden complexity in CQ resistance evolution, suggesting that pfaat1 may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport.

With global progress towards malaria reduction stalling, further analysis of the epidemiology is required, particularly in countries with the highest burden. National surveys have mostly analysed infection prevalence, while large-scale data on parasite density and different developmental forms are rarely available. In Nigeria, the country with the largest burden globally, blood slide microscopy of children up to five years of age was conducted in the 2018 National Demographic and Health Survey, and infection prevalence previously reported. In the current study, malaria parasite density measurements are reported and analysed for 7783 of the children sampled across the 36 states within the six geopolitical zones of the country, with asexual and sexual stages considered separately and comparisons of infections with different malaria parasite species being performed. Across all states of Nigeria, there was a positive correlation between mean asexual parasite density in peripheral blood of infected individuals and prevalence of infection in the community (Spearman’s rho = 0.39, P = 0.02). Asexual parasite densities were highest in the northern geopolitical zones (geometric means > 2000 μL-1), extending the evidence of exceptionally high infection burden in many areas. Sexual parasite prevalence in each state was highly correlated with asexual parasite prevalence (Spearman’s rho = 0.70, P < 0.001), although sexual parasite densities were low (geometric means < 100 μL-1in all zones). Infants had lower parasite densities than children above one year of age, but there were no differences between male and female children. Most infections were ofP. falciparum, which had higher asexual densities but lower sexual parasite densities thanP. malariaeorP. ovalemono-infections. However, mixed species infections had the highest asexual parasite densities. It is recommended that future large surveys in high burden countries measure parasite densities as well as developmental stages and species, to improve the quality of malaria epidemiology and tracking of future changes.

Experimental studies on the biology of malaria parasites have been mostly based on laboratory-adapted lines, but there is limited understanding of how these may differ from parasites in natural infections. Loss-of-function mutants have previously been shown to emerge during culture of somePlasmodium falciparumclinical isolates, in analyses that focused on single-genotype infections. The present study included a broader array of isolates, mostly representing multiple-genotype infections which are more typical in areas where malaria is highly endemic. Genome sequence data from multiple time points during several months of culture adaptation of 28 West African isolates were analysed, including previously available sequences along with new genome sequences from additional isolates and timepoints. Some genetically complex isolates eventually became fixed over time to single surviving genotypes in culture, whereas others retained diversity although proportions of genotypes varied over time. Drug-resistance allele frequencies did not show overall directional changes, suggesting that resistance-associated costs are not the main causes of fitness differences among parasites in culture. Loss-of-function mutants emerged during culture in several of the multiple-genotype isolates, affecting genes (includingAP2-HS, EPACandSRPK1) for which loss-of-function mutants were previously seen to emerge in single-genotype isolates. Parasite clones were derived by limiting dilution from six of the isolates, and sequencing identifiedde novovariants not detected in the bulk isolate sequences. Interestingly, most of these were nonsense mutants and frameshifts disrupting the coding sequence ofEPAC, the gene with the largest number of independent nonsense mutants previously identified in laboratory-adapted lines. Analysis of Identity-By-Descent to explore relatedness among clones revealed co-occurring non-identical sibling parasites, illustrative of the natural genetic structure within parasite populations.

Malaria parasites break down host hemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter,pfcrt, located in the digestive vacuole membrane, confer CQ resistance inPlasmodium falciparum, but typically affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on an amino acid transporter (pfaat1)S258Lvariant which increased from 0-87% in frequency between 1984 and 2014 in parallel withpfcrt1K76T. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containingpfaat1that is selected by CQ treatment. Gene editing demonstrated thatpfaat1S258Lpotentiates CQ-resistance but at a cost of reduced fitness, whilepfaat1F313S, a common Southeast Asian polymorphism, reduces CQ-resistance while restoring fitness. Our analyses reveal hidden complexity in CQ-resistance evolution, suggesting thatpfaat1may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport.

The merozoite surface protein MSPDBL2 of Plasmodium falciparum is under strong balancing selection and is a target of naturally acquired antibodies. Remarkably, MSPDBL2 is expressed in only a minority of mature schizonts of any cultured parasite line, and mspdbl2 gene transcription increases in response to overexpression of the gametocyte development inducer GDV1, so it is important to understand its natural expression. Here, MSPDBL2 in mature schizonts was analysed in the first ex vivo culture cycle of 96 clinical isolates from four populations with varying levels of infection endemicity in different West African countries, by immunofluorescence microscopy with antibodies against a conserved region of the protein. In most isolates, less than 1% of mature schizonts were positive for MSPDBL2 (median of 0.6% overall), but the frequency distribution was highly skewed as nine isolates had more than 3% schizonts positive and one had 73% positive. To investigate whether expression of other gene loci correlated with MSPDBL2 expression, whole transcriptome sequencing was performed on schizont-enriched material from 17 of the clinical isolates with a wide range of proportions of schizonts positive. Transcripts of particular parasite genes were highly significantly positively correlated with MSPDBL2 positivity in schizonts as well as with mspdbl2 gene transcript levels, with overrepresentation of genes previously implicated as likely to be involved in gametocytogenesis, but not including the gametocytogenesis master regulator ap2g. Although MSPDBL2 is apparently not directly involved in sexual commitment, it marks a co-occurring developmental subpopulation that may be functionally distinct within blood stage infections.

Asexual blood stage malaria parasites must produce sexual progeny in order to infect mosquitoes. It is therefore important to understand the scope of variation in sexual commitment rates of the most important human parasite Plasmodium falciparum, and to assess responses to conditions that may affect this commitment. First, using a genetically modified line with inducible elevated sexual commitment we compared two different methods of measuring the rates, either staining the Pfs16 marker in very early gametocytes or counting gametocytes after further development. The methods showed correlated estimates in evaluation across multiple independent replicate assays, with higher sensitivity and precision being achieved by Pfs16 marker detection in early gametocytes, so this method was used to survey a wide range of P. falciparum lines. A panel of six recently culture-established clinical isolates from Ghana showed mean sexual commitment rates per cycle ranging from 3.3% to 12.2%, with some significant differences among isolates as well as variation among biological replicates of each. Comparing 13 long-term laboratory-adapted lines of diverse origins, eight had sexual commitment rates similar to those of the recent clinical isolates, with means across multiple biological replicates for each line ranging from 4.7% to 14.5%, while four lines had significantly lower rates with means ranging from 0.3 to 1.6%, and one line with a non-functional ap2-g gene always showed zero sexual commitment. Among a subset of parasite lines, adding choline to suppress commitment had quantitatively variable effects, although these were significant in most assays of lines that had relatively high rates. This study indicates the importance of multiple assay replicates and comparisons of diverse isolates to understand natural and culture-induced variation in this key reproductive trait, relevant to investigating potential effects of parasite adaptation on transmission.Author SummaryOnly sexual malaria parasites are transmitted from humans to mosquitoes, so it is vital to understand variation in sexual commitment rates of blood stage malaria parasites, and responses to conditions that affect this. We compared two different methods of measuring the rates in a Plasmodium falciparum line with engineered sexual commitment variation, demonstrating higher sensitivity and precision by detection of an early differentiation marker, and this method was then used to survey diverse P. falciparum lines in culture with extensive biological replicate measurements. Recent clinical isolates from Ghana showed mean sexual commitment rates per cycle ranging from 3% to 12%, with significant differences among isolates as well as variation among biological replicates. Long-term laboratory-adapted lines of diverse origins had a wide range, most having rates similar to those of the clinical isolates, while a minority consistently had much lower or zero rates. There was quantitative variation among lines in the effects of adding choline to suppress commitment, although most assays of lines that had relatively high rates showed significant effects. Performing multiple biological replicates and comparisons of recent parasite isolates is vital to understand intra-specific variation in this important reproductive trait, and move towards investigating direct effects on disease transmission.

Background In areas where Plasmodium falciparum malaria is seasonal, a dry season reservoir of blood-stage infection is essential for initiating transmission during the following wet season. Methods In The Gambia, a cohort of 42 individuals with quantitative polymerase chain reaction-positive P falciparum infections at the end of the transmission season (December) were followed monthly until the end of the dry season (May) to evaluate infection persistence. The influence of human host and parasitological factors was investigated. Results A large proportion of individuals infected at the end of the wet season had detectable infections until the end of the dry season (40.0%; 16 of 40). At the start of the dry season, the majority of these persistent infections (82%) had parasite densities >10 p/µL compared to only 5.9% of short-lived infections. Persistent infections (59%) were also more likely to be multiclonal than short-lived infections (5.9%) and were associated with individuals having higher levels of P falciparum-specific antibodies (P = .02). Conclusions Asymptomatic persistent infections were multiclonal with higher parasite densities at the beginning of the dry season. Screening and treating asymptomatic infections during the dry season may reduce the human reservoir of malaria responsible for initiating transmission in the wet season.

Most malaria in Malaysia is caused by Plasmodium knowlesi parasites through zoonotic infection from macaque reservoir hosts. We obtained genome sequences from 28 clinical infections in Peninsular Malaysia to clarify the emerging parasite population structure and test for evidence of recent adaptation. The parasites all belonged to a major genetic population of P. knowlesi (cluster 3) with high genomewide divergence from populations occurring in Borneo (clusters 1 and 2). We also observed unexpected local genetic subdivision; most parasites belonged to 2 subpopulations sharing a high level of diversity except at particular genomic regions, the largest being a region of chromosome 12, which showed evidence of recent directional selection. Surprisingly, we observed a third subpopulation comprising P. knowlesi infections that were almost identical to each other throughout much of the genome, indicating separately maintained transmission and recent genetic isolation. Each subpopulation could evolve and present a broader health challenge in Asia.

BackgroundIn areas where Plasmodium falciparum malaria is highly seasonal, a dry season reservoir of blood-stage infection is essential for initiating transmission during the following wet season, bridging transmission seasons several months apart. Understanding infections during the dry season could thus inform approaches for malaria control.MethodsIn The Gambia, a cohort of 42 individuals with qPCR positive P. falciparum infections at the end of the transmission season (December) were followed monthly until the end of the dry season (May) to evaluate the duration of detectable infections. The influence of human host (age, sex, haemoglobin concentration and genotype, and P. falciparum-specific antibodies), and parasitological (parasite density, gametocyte density and genotypic multiplicity of infection) factors was investigated.ResultsA large proportion of individuals infected at the end of the wet season had detectable infections until the end of the dry season (40.0%; 16/40), with the majority of these infections also harbouring gametocytes (81.3%; 13/16). 22 infections were classified as persistent (detectable for at least 3 months), 17 were classified as short-lived (undetectable within 2 months), and 3 were treated (due to symptoms). At the start of the dry season, the majority of persistent infections (82%; 18/22) had parasite densities >10 p/µL compared to only 5.9% (1/17) of short-lived infections. Persistent infections (59%; 13/22) were also more likely to be multi-clonal than short-lived infections (5.9%; 1/17), they were most common in 5 to 15 year old children (63%; 12/19), and were associated with individuals having higher levels of P. falciparum-specific antibodies (p = 0.058).ConclusionsAsymptomatic persistent dry season infections in The Gambia were multiclonal with higher parasite densities at the beginning of the dry season, mostly occurring in school age children and adults with higher P. falciparum-specific antibodies. Screening and treating asymptomatic, malaria-infected individuals during the dry season may reduce the human reservoir of malaria responsible initiating transmission in the wet-season.

Molecular epidemiology has been central to uncovering P. knowlesi as an important cause of human malaria in Southeast Asia, and to understanding the complex nature of this zoonosis. Species-specific parasite detection and characterization of sequences were vital to show that P. knowlesi was distinct from the human parasite species that had been presumed to cause all malaria. With established sensitive and specific molecular detection tools, surveys subsequently indicated the distribution of P. knowlesi infections in humans, wild primate reservoir host species, and mosquito vector species. The importance of studying P. knowlesi genetic polymorphism was indicated initially by analysing a few nuclear gene loci as well as the mitochondrial genome, and subsequently by multi-locus microsatellite analyses and whole-genome sequencing. Different human infections generally have unrelated P. knowlesi genotypes, acquired from the diverse local parasite reservoirs in macaques. However, individual human infections are usually less genetically complex than those of wild macaques which experience more frequent superinfection with different P. knowlesi genotypes. Multi-locus analyses have revealed deep population subdivisions within P. knowlesi, which are structured both geographically and in relation to different macaque reservoir host species. Simplified genotypic discrimination assays now enable efficient large-scale surveillance of the sympatric P. knowlesi subpopulations within Malaysian Borneo. The whole-genome sequence analyses have also identified loci under recent positive natural selection in the P. knowlesi genome, with evidence that different loci are affected in different populations. These provide a foundation to understand recent adaptation of the zoonotic parasite populations, and to track and interpret future changes as they emerge.

Host genetic factors can confer resistance against malaria, raising the question of whether this has led to evolutionary adaptation of parasite populations. In this study we investigated the correlation between host and parasite genetic variation in 4,171 Gambian and Kenya children ascertained with severe malaria due to Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and variation in three regions of the parasite genome, including nonsynonymous variants in the acyl-CoA synthetase family member PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The HbS-associated parasite alleles are in strong linkage disequilibrium and have frequencies which covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome, and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

Periodic fever is the most characteristic clinical feature of human malaria1-3, but how parasites survive febrile episodes is not known. While Plasmodium spp. genomes encode a full complement of chaperones4, they lack an ortholog of the conserved transcription factor HSF1, which in most eukaryotes activates the expression of key chaperones upon heat shock (HS)5-8. Here we identified PfAP2-HS, a transcription factor of the ApiAP2 family9-11, as the key regulator of the P. falciparum protective HS response. The PfAP2-HS-dependent HS response is largely restricted to rapid activation of hsp70-1, the predominant direct target of PfAP2-HS, and hsp90. Deletion of PfAP2-HS dramatically reduced HS survival and also resulted in severe growth defects at 37°C, but not at 35°C, and increased sensitivity to imbalances in protein homeostasis (proteostasis) produced by artemisinin, the current frontline antimalarial drug12,13. These results demonstrate that PfAP2-HS contributes to general maintenance of proteostasis and drives a rapid chaperone-based protective response against febrile temperatures. While several ApiAP2 transcription factors regulate life cycle transitions in malaria parasites11,14,15, PfAP2-HS is the first identified Plasmodium transcription factor that controls a protective response to a within-host environmental challenge.

Periodic fever is a characteristic clinical feature of human malaria, but how parasites survive febrile episodes is not known. Although the genomes of Plasmodium species encode a full set of chaperones, they lack the conserved eukaryotic transcription factor HSF1, which activates the expression of chaperones following heat shock. Here, we show that PfAP2-HS, a transcription factor in the ApiAP2 family, regulates the protective heat-shock response in Plasmodium falciparum. PfAP2-HS activates the transcription of hsp70-1 and hsp90 at elevated temperatures. The main binding site of PfAP2-HS in the entire genome coincides with a tandem G-box DNA motif in the hsp70-1 promoter. Engineered parasites lacking PfAP2-HS have reduced heat-shock survival and severe growth defects at 37 °C but not at 35 °C. Parasites lacking PfAP2-HS also have increased sensitivity to imbalances in protein homeostasis (proteostasis) produced by artemisinin, the frontline antimalarial drug, or the proteasome inhibitor epoxomicin. We propose that PfAP2-HS contributes to the maintenance of proteostasis under basal conditions and upregulates specific chaperone-encoding genes at febrile temperatures to protect the parasite against protein damage.

BackgroundGlobal progress in reducing malaria has stalled since 2015. Analysis of the situation is particularly needed in Nigeria, the country with by far the largest share of the burden, where approximately a quarter of all cases in the world are estimated to occur.MethodsWe analysed data from three nationwide surveys (Malaria Indicator Surveys in 2010 and 2015, and a National Demographic and Health Survey in 2018), with malaria parasite prevalence in children under five years of age determined by sampling from all 36 states of Nigeria, and blood slide microscopy performed in the same accredited laboratory for all samples. Changes over time were evaluated by calculating prevalence ratio (PR) values with 95% confidence intervals (CI) for each state, together with Mantel-Haenszel adjusted prevalence ratios (PRadj) for each of the six major geopolitical zones of the country.ResultsBetween 2010 and 2018 there were significant reductions in parasite prevalence in 25 states, but not in the remaining 11 states. Prevalence decreased most in southern zones of the country (South West PRadj = 0.53; South East PRadj = 0.59; South South PRadj = 0.51) and the North Central zone (PRadj = 0.36). Changes in the north were less marked, but were significant and indicated overall reductions by more than 20% (North West PRadj = 0.74; North East PRadj = 0.76). Changes in the south occurred mostly between 2010 and 2015, whereas those in the north were more gradual and most continued after 2015. Recent changes were not correlated with survey-reported variation in use of preventive measures.ConclusionReductions in malaria infection in children under five have occurred in most individual states in Nigeria since 2010, but substantial geographical variation in the timing and extent indicate challenges to be overcome to enable global malaria reduction.

Molecular epidemiology has been central to uncovering P. knowlesi as an important cause of human malaria in Southeast Asia, and to understanding the complex nature of this zoonosis. Species-specific parasite detection and characterization of sequences were vital to show that P. knowlesi was distinct from the human parasite species that had been presumed to cause all malaria. With established sensitive and specific molecular detection tools, surveys subsequently indicated the distribution of P. knowlesi infections in humans, wild primate reservoir host species, and mosquito vector species. The importance of studying P. knowlesi genetic polymorphism was indicated initially by analysing a few nuclear gene loci as well as the mitochondrial genome, and subsequently by multi-locus microsatellite analyses and whole-genome sequencing. Different human infections generally have unrelated P. knowlesi genotypes, acquired from the diverse local parasite reservoirs in macaques. However, individual human infections are usually less genetically complex than those of wild macaques which experience more frequent superinfection with different P. knowlesi genotypes. Multi-locus analyses have revealed deep population subdivisions within P. knowlesi, which are structured both geographically and in relation to different macaque reservoir host species. Simplified genotypic discrimination assays now enable efficient large-scale surveillance of the sympatric P. knowlesi subpopulations within Malaysian Borneo. The whole-genome sequence analyses have also identified loci under recent positive natural selection in the P. knowlesi genome, with evidence that different loci are affected in different populations. These provide a foundation to understand recent adaptation of the zoonotic parasite populations, and to track and interpret future changes as they emerge.

Population genetic analysis revealed that Plasmodium knowlesi infections in Malaysian Borneo are caused by 2 divergent parasites associated with long-tailed (cluster 1) and pig-tailed (cluster 2) macaques. Because the transmission ecology is likely to differ for each macaque species, we developed a simple genotyping PCR to efficiently distinguish between and survey the 2 parasite subpopulations. This assay confirmed differences in the relative proportions in areas of Kapit division of Sarawak state, consistent with multilocus microsatellite analyses. Analyses of 1,204 human infections at Kapit Hospital showed that cluster 1 caused approximately two thirds of cases with no significant temporal changes from 2000 to 2018. We observed an apparent increase in overall numbers in the most recent 2 years studied, driven mainly by increased cluster 1 parasite infections. Continued monitoring of the frequency of different parasite subpopulations and correlation with environmental alterations are necessary to determine whether the epidemiology will change substantially.

Pathogen multiplication rate is theoretically an important determinant of virulence, although often poorly understood and difficult to measure accurately. We show intrinsic asexual blood stage multiplication rate variation of the major human malaria parasite Plasmodium falciparum to be associated with blood-stage infection intensity in patients. A panel of clinical isolates from a highly endemic West African population was analysed repeatedly during five months of continuous laboratory culture, showing a range of exponential multiplication rates at all timepoints tested, mean rates increasing over time. All isolates had different genome sequences, many containing within-isolate diversity that decreased over time in culture, but increases in multiplication rates were not primarily attributable to genomic selection. New mutants, including premature stop codons emerging in a few isolates, did not attain sufficiently high frequencies to substantially affect overall multiplication rates. Significantly, multiplication rate variation among the isolates at each of the assayed culture timepoints robustly correlated with parasite levels seen in patients at clinical presentation, indicating innate parasite control of multiplication rate that contributes to virulence., Lindsay Stewart et al. analyze clinical isolates of the human malaria parasite Plasmodium falciparum from a highly endemic West African population and show that intrinsic multiplication rate variation is associated with blood-stage infection intensity. Their results indicate that parasite control of multiplication contributes to virulence.

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.

Performance of the malaria Rapid Diagnostic Test (RDT) OptiMal-IT® was evaluated in Mauritania where malaria is low and dependent on a short transmission season. Slide microscopy was considered as the reference method of diagnosis. Febrile patients with suspected malaria were recruited from six health facilities, 3 urban and 3 rural, during two periods (December 2011 to February 2012, and August 2012 to March 2013). Overall, 780 patients were sampled, with RDT and thick blood film microscopy results being obtained for 759 of them. Out of 774 slides examined, of which 200 were positive, P. falciparum and P. vivax mono-infections were detected in 63.5% (127) and 29.5% (59), while P. falciparum/P. vivax coinfections were detected in 7% (14). Both species were observed in all study sites, although in significantly different proportions. The proportions of thick blood film and OptiMal-IT® RDT positive individuals was 26.3% and 30.3% respectively. Sensitivity and specificity of OptiMal-IT® RDT were 89% [95% CI, 84.7-93.3] and 91.1% [88.6-93.4]. Positives and negative predictive values were 78.1% [72.2-83.7] and 95.9% [94.1-97.5]. These diagnostic values are similar to those generally reported elsewhere, and support the use of RDTs as the main diagnostic tool for malaria in Mauritanian health facilities. In the future, choice of RDTs to be used must take account of thermostability in a hot, dry environment and their ability to detect P. falciparum and P. vivax.

We conducted a genome-wide association study of host resistance to severePlasmodium falciparummalaria in over 17,000 individuals from 11 malaria-endemic countries, undertaking a wide ranging analysis which identifies five replicable associations with genome-wide levels of evidence. Our findings include a newly implicated variant on chromosome 6 associated with risk of cerebral malaria, and the discovery of an erythroid-specific transcription start site underlying the association inATP2B4. Previously reported HLA associations cannot be replicated in this dataset. We estimate substantial heritability of severe malaria (h2~ 23%), of which around 10% is explained by the currently identified associations. Our dataset will provide a major building block for future research on the genetic determinants of disease in these diverse human populations.

Background Malaria parasites are genetically polymorphic and phenotypically plastic. In studying transcriptome variation among parasites from different infections, it is challenging to overcome potentially confounding technical and biological variation between samples. We investigate variation in the major human parasite Plasmodium falciparum, generating RNA-seq data on multiple independent replicate sample preparations of merozoite-containing intra-erythrocytic schizonts from a panel of clinical isolates and from long-term laboratory-adapted clones, with a goal of robustly identifying differentially expressed genes. Results Analysis of biological sample replicates shows that increased numbers improve the true discovery rate of differentially expressed genes, and that six independent replicates of each parasite line allowed identification of most differences that could be detected with larger numbers. For highly expressed genes, focusing on the top quartile at schizont stages, there was more power to detect differences. Comparing cultured clinical isolates and laboratory-adapted clones, genes more highly expressed in the laboratory-adapted clones include those encoding an AP2 transcription factor (PF3D7_0420300), a ubiquitin-binding protein and two putative methyl transferases. In contrast, higher expression in clinical isolates was seen for the merozoite surface protein gene dblmsp2, proposed to be a marker of schizonts forming merozoites committed to sexual differentiation. Variable expression was extremely strongly, but not exclusively, associated with genes known to be targeted by Heterochromatin Protein 1. Clinical isolates show variable expression of several known merozoite invasion ligands, as well as other genes for which new RT-qPCR assays validate the quantitation and allow characterisation in samples with more limited material. Expression levels of these genes vary among schizont preparations of different clinical isolates in the first ex vivo cycle in patient erythrocytes, but mean levels are similar to those in continuously cultured clinical isolates. Conclusions Analysis of multiple biological sample replicates greatly improves identification of genes variably expressed between different cultured parasite lines. Clinical isolates recently established in culture show differences from long-term adapted clones in transcript levels of particular genes, and are suitable for analyses requiring biological replicates to understand parasite phenotypes and variable expression likely to be relevant in nature. Electronic supplementary material The online version of this article (10.1186/s12864-018-5257-x) contains supplementary material, which is available to authorized users.

Background: Within Plasmodium falciparum merozoite surface protein 1 (MSP1), the N-terminal block 2 region is a highly polymorphic target of naturally acquired antibody responses. The antigenic diversity is determined by complex repeat sequences as well as non-repeat sequences, grouping into three major allelic types that appear to be maintained within populations by natural selection. Within these major types, many distinct allelic sequences have been described in different studies, but the extent and significance of the diversity remains unresolved.Methods: To survey the diversity more extensively, block 2 allelic sequences in the msp1 gene were characterized in 2400 P. falciparum infection isolates with whole genome short read sequence data available from the Pf3K project, and compared with the data from previous studies.Results: Mapping the short read sequence data in the 2400 isolates to a reference library of msp1 block 2 allelic sequences yielded 3815 allele scores at the level of major allelic family types, with 46% of isolates containing two or more of these major types. Overall frequencies were similar to those previously reported in other samples with different methods, the K1-like allelic type being most common in Africa, MAD20-like most common in Southeast Asia, and RO33-like being the third most abundant type in each continent. The rare MR type, formed by recombination between MAD20-like and RO33-like alleles, was only seen in Africa and very rarely in the Indian subcontinent but not in Southeast Asia. A combination of mapped short read assembly approaches enabled 1522 complete msp1 block 2 sequences to be determined, among which there were 363 different allele sequences, of which 246 have not been described previously. In these data, the K1-like msp1 block 2 alleles are most diverse and encode 225 distinct amino acid sequences, compared with 123 different MAD20-like, 9 RO33-like and 6 MR type sequences. Within each of the major types, the different allelic sequences show highly skewed geographical distributions, with most of the more common sequences being detected in either Africa or Asia, but not in both.Conclusions: Allelic sequences of this extremely polymorphic locus have been derived from whole genome short read sequence data by mapping to a reference library followed by assembly of mapped reads. The catalogue of sequence variation has been greatly expanded, so that there are now more than 500 different msp1 block 2 allelic sequences described. This provides an extensive reference for molecular epidemiological genotyping and sequencing studies, and potentially for design of a multi-allelic vaccine.

IntroductionGlobal progress in reducing malaria has stalled since 2015. Analysis of the situation is particularly needed in Nigeria, the country with by far the largest share of the burden, where approximately a quarter of all cases in the world are estimated to occur.MethodsWe analysed data from three nationwide surveys (Malaria Indicator Surveys in 2010 and 2015 and a National Demographic and Health Survey in 2018), with malaria parasite prevalence in children under 5 years of age determined by sampling from all 36 states of Nigeria, and blood slide microscopy performed in the same accredited laboratory for all samples. Changes over time were evaluated by calculating prevalence ratio (PR) values with 95% CIs for each state, together with Mantel-Haenszel-adjusted PRs (PRadj) for each of the six major geopolitical zones of the country.ResultsBetween 2010 and 2018, there were significant reductions in parasite prevalence in 25 states, but not in the remaining 11 states. Prevalence decreased most in southern zones of the country (South West PRadj=0.53; South East PRadj=0.59; South South PRadj=0.51) and the North Central zone (PRadj=0.36). Changes in the north were less marked, but were significant and indicated overall reductions by more than 20% (North-West PRadj=0.74; North East PRadj=0.70). Changes in the south occurred mostly between 2010 and 2015, whereas those in the north were more gradual and most continued after 2015. Recent changes were not correlated with survey-reported variation in use of preventive measures.ConclusionReductions in malaria infection in children under 5 have occurred in most individual states in Nigeria since 2010, but substantial geographical variation in the timing and extent indicate challenges to be overcome to enable global malaria reduction.

Plasmodium cynomolgi is a malaria parasite that typically infects Asian macaque monkeys, and humans on rare occasions. P. cynomolgi serves as a model system for the human malaria parasite Plasmodium vivax, with which it shares such important biological characteristics as formation of a dormant liver stage and a preference to invade reticulocytes. While genomes of three P. cynomolgi strains have been sequenced, genetic diversity of P. cynomolgi has not been widely investigated. To address this we developed the first panel of P. cynomolgi microsatellite markers to genotype eleven P. cynomolgi laboratory strains and 18 field isolates from Sarawak, Malaysian Borneo. We found diverse genotypes among most of the laboratory strains, though two nominally different strains were found to be genetically identical. We also investigated sequence polymorphism in two erythrocyte invasion gene families, the reticulocyte binding protein and Duffy binding protein genes, in these strains. We also observed copy number variation in rbp genes.

Erythrocyte invasion by Plasmodium falciparum merozoites is an essential step for parasite survival and proliferation. Invasion is mediated by multiple ligands, which could be promising vaccine targets. The usage and sequence of these ligands differs between parasites, yet most studies of them have been carried out in only a few laboratory-adapted lines. To understand the true extent of natural variation in invasion phenotypes and prioritize vaccine candidates on a relevant evidence base, we need to develop and apply standardized assays to large numbers of field isolates. The West African Merozoite Invasion Network (WAMIN) has been formed to meet these goals, expand training in Plasmodium phenotyping, and perform large-scale field phenotyping studies in order to prioritize blood stage vaccine candidates.

The pathogenesis of infectious diseases depends on the interaction of host and pathogen. In Plasmodium falciparum malaria, host and parasite processes can be assessed by dual RNA-sequencing of blood from infected patients. Here we performed dual transcriptome analyses on samples from 46 malaria-infected Gambian children to reveal mechanisms driving the systemic pathophysiology of severe malaria. Integrating these transcriptomic data with estimates of parasite load and detailed clinical information allowed consideration of potentially confounding effects due to differing leukocyte proportions in blood, parasite developmental stage, and whole-body pathogen load. We report hundreds of human and parasite genes differentially expressed between severe and uncomplicated malaria, with distinct profiles associated with coma, hyperlactatemia, and thrombocytopenia. High expression of neutrophil granule-related genes was consistently associated with all severe malaria phenotypes. We observed severity-associated variation in the expression of parasite genes which determine cytoadhesion to vascular endothelium, rigidity of infected erythrocytes, and parasite growth rate. Up to 99% of human differential gene expression in severe malaria was driven by differences in parasite load, whereas parasite gene expression showed little association with parasite load. Co-expression analyses revealed interactions between human and P. falciparum, with prominent co-regulation of translation genes in severe malaria between host and parasite. Multivariate analyses suggested that increased expression of granulopoiesis and interferon-γ related genes, together with inadequate suppression of type-1 interferon signalling, best explained severity of infection. These findings provide a framework for understanding the contributions of host and parasite to the pathogenesis of severe malaria and identifying targets for adjunctive therapy.

Malaria parasite genes exhibit variation in both sequence and expression level. There is much information on sequence polymorphism, but less resolution on natural variation in transcriptomes of parasites at specific developmental stages. This is largely because it is challenging to obtain highly replicated sampling of transcriptomes to overcome potentially confounding technical and biological variation. We address the issue in the major human parasite Plasmodium falciparum by obtaining RNA-seq profiles of multiple independent replicate preparations of mature schizont-stage parasites from a panel of clinical isolates recently established in culture and from long-term laboratory-adapted clones. With a goal of robustly identifying variably expressed genes, we show that increasing the numbers of biological sample replicates greatly improves the discovery rate. Generally, six independent replicates of each parasite culture is recommendable as being significantly to lower numbers, although for highly expressed genes variable expression can be detected when fewer replicates are available. A broad comparison identifies genes differing in relative expression between cultured clinical isolates and laboratory-adapted clones. Genes more highly expressed in the laboratory-adapted clones include an AP2 transcription factor gene Pf3D7_0420300 and putative methyl transferase genes. The variable expression of several known merozoite invasion ligands is confirmed, and previously uncharacterised genes are shown to be differentially expressed among clinical isolates. New RT-qPCR assays validate the variation in transcript levels of these genes, and allow quantitation of expression to be extended to a wider panel of clinical isolate samples. These variably expressed genes are new candidates for investigation as potential determinants of alternative parasite developmental pathways or targets of immunity.Author summaryUnderstanding parasite diversity and adaptation may require characterisation of gene expression variation, and is vital if chemotherapeutic or vaccine development is to consider new candidate targets, but it is technically challenging to generate precise data on clinical isolates. Here, we analyse the transcriptomes of mature Plasmodium falciparum schizonts using RNA-sequencing, using large numbers of biological replicate samples to minimise the impact of inter-replicate variation on observed patterns of differential expression. This identifies genes that are differentially expressed in long term laboratory-adapted parasites and recently cultured clinical isolates, as well as among different clinical isolates. In additional samples of schizonts grown in the first cycle ex vivo prior to any erythrocyte invasion, expression levels of a selected panel of these genes vary among isolates, but mean levels are similar to those in the continuously cultured clinical isolates, indicating that the latter are useful for experimental studies requiring biological replication.

SUMMARY Transcriptomics, the analysis of genome-wide RNA expression, is a common approach to investigate host and pathogen processes in infectious diseases. Technical and bioinformatic advances have permitted increasingly thorough analyses of the association of RNA expression with fundamental biology, immunity, pathogenesis, diagnosis, and prognosis. Transcriptomic approaches can now be used to realize a previously unattainable goal, the simultaneous study of RNA expression in host and pathogen, in order to better understand their interactions. This exciting prospect is not without challenges, especially as focus moves from interactions in vitro under tightly controlled conditions to tissue- and systems-level interactions in animal models and natural and experimental infections in humans. Here we review the contribution of transcriptomic studies to the understanding of malaria, a parasitic disease which has exerted a major influence on human evolution and continues to cause a huge global burden of disease. We consider malaria a paradigm for the transcriptomic assessment of systemic host-pathogen interactions in humans, because much of the direct host-pathogen interaction occurs within the blood, a readily sampled compartment of the body. We illustrate lessons learned from transcriptomic studies of malaria and how these lessons may guide studies of host-pathogen interactions in other infectious diseases. We propose that the potential of transcriptomic studies to improve the understanding of malaria as a disease remains partly untapped because of limitations in study design rather than as a consequence of technological constraints. Further advances will require the integration of transcriptomic data with analytical approaches from other scientific disciplines, including epidemiology and mathematical modeling.

Plasmodium falciparum merozoites use diverse alternative erythrocyte receptors for invasion and variably express cognate ligands encoded by the erythrocyte binding antigen ( eba ) and reticulocyte binding-like homologue ( Rh ) gene families. Previous analyses conducted on parasites from single populations in areas of endemicity revealed a wide spectrum of invasion phenotypes and expression profiles, although comparisons across studies have been limited by the use of different protocols. For direct comparisons within and among populations, clinical isolates from three different West African sites of endemicity (in Ghana, Guinea, and Senegal) were cryopreserved and cultured ex vivo after thawing in a single laboratory to assay invasion of target erythrocytes pretreated with enzymes affecting receptor subsets. Complete invasion assay data from 67 isolates showed no differences among the populations in the broad range of phenotypes measured by neuraminidase treatment (overall mean, 40.6% inhibition) or trypsin treatment (overall mean, 83.3% inhibition). The effects of chymotrypsin treatment (overall mean, 79.2% inhibition) showed heterogeneity across populations (Kruskall-Wallis P = 0.023), although the full phenotypic range was seen in each. Schizont-stage transcript data for a panel of 8 invasion ligand genes ( eba175 , eba140 , eba181 , Rh1 , Rh2a , Rh2b , Rh4 , and Rh5 ) were obtained for 37 isolates, showing similar ranges of variation in each population except that eba175 levels tended to be higher in parasites from Ghana than in those from Senegal (whereas levels of eba181 and Rh2b were lower in parasites from Ghana). The broad diversity in invasion phenotypes and gene expression seen within each local population, with minimal differences among them, is consistent with a hypothesis of immune selection maintaining parasite variation.

Summary Background Although many malaria control programmes in sub-Saharan Africa use indoor residual spraying with long-lasting insecticidal nets (LLINs), the two studies assessing the benefit of the combination of these two interventions gave conflicting results. We aimed to assess whether the addition of indoor residual spraying to LLINs provided a significantly different level of protection against clinical malaria in children or against house entry by vector mosquitoes. Methods In this two-arm cluster, randomised, controlled efficacy trial we randomly allocated clusters of Gambian villages using a computerised algorithm to LLINs alone (n=35) or indoor residual spraying with dichlorodiphenyltrichloroethane plus LLINs (n=35). In each cluster, 65–213 children, aged 6 months to 14 years, were surveyed at the start of the 2010 transmission season and followed in 2010 and 2011 by passive case detection for clinical malaria. Exposure to parasite transmission was assessed by collection of vector mosquitoes with both light and exit traps indoors. Primary endpoints were the incidence of clinical malaria assessed by passive case detection and number of Anopheles gambiae sensu lato mosquitoes collected per light trap per night. Intervention teams had no role in data collection and the data collection teams were not informed of the spray status of villages. The trial is registered at the ISRCTN registry, number ISRCTN01738840. Findings LLIN coverage in 2011 was 3510 (93%) of 3777 children in the indoor residual spraying plus LLIN group and 3622 (95·5%) of 3791 in the LLIN group. In 2010, 7845 children were enrolled, 7829 completed passive case detection, and 7697 (98%) had complete clinical and covariate data. In 2011, 7009 children remained in the study, 648 more were enrolled, 7657 completed passive case detection, and 7545 (98·5%) had complete data. Indoor residual spraying coverage per cluster was more than 80% for both years in the indoor residual spraying plus LLIN group. Incidence of clinical malaria was 0·047 per child-month at risk in the LLIN group and 0·044 per child-month at risk in the indoor residual spraying plus LLIN group in 2010, and 0·032 per child-month at risk in the LLIN group and 0·034 per child-month at risk in the indoor residual spraying plus LLIN group in 2011. The incident rate ratio was 1·08 (95% CI 0·80–1·46) controlling for confounders and cluster by mixed-effect negative binomial regression on all malaria attacks for both years. No significant difference was recorded in the density of vector mosquitoes caught in light traps in houses over the two transmission seasons; the mean number of A gambiae sensu lato mosquitoes per trap per night was 6·7 (4·0–10·1) in the LLIN group and 4·5 (2·4–7·4) in the indoor residual spraying plus LLIN group (p=0·281 in the random-effects linear regression model). Interpretation We identified no significant difference in clinical malaria or vector density between study groups. In this area with high LLIN coverage, moderate seasonal transmission, and susceptible vectors, indoor residual spraying did not provide additional benefit. Funding UK Medical Research Council.

Highlights • Discovery of vaccine candidate antigens by parasite genome sequence analyses. • Genetic crosses, linkage group selection, and functional studies on parasites. • Characterizing developmental and epigenetic variation alongside allelic polymorphism. • Selection by naturally acquired immune responses helps to focus vaccine design., More human death and disease is caused by malaria parasites than by all other eukaryotic pathogens combined. As early as the sequencing of the first human genome, malaria parasite genomics was prioritized to fuel the discovery of vaccine candidate antigens. This stimulated increased research on malaria, generating new understanding of the cellular and molecular mechanisms of infection and immunity. This review of recent developments illustrates how new approaches in parasite genomics, and increasingly large amounts of data from population studies, are helping to identify antigens that are promising lead targets. Although these results have been encouraging, effective discovery and characterization need to be coupled with more innovation and funding to translate findings into newly designed vaccine products for clinical trials.

Multilocus microsatellite genotyping of Plasmodium knowlesi isolates previously indicated 2 divergent parasite subpopulations in humans on the island of Borneo, each associated with a different macaque reservoir host species. Geographic divergence was also apparent, and independent sequence data have indicated particularly deep divergence between parasites from mainland Southeast Asia and Borneo. To resolve the overall population structure, multilocus microsatellite genotyping was conducted on a new sample of 182 P. knowlesi infections (obtained from 134 humans and 48 wild macaques) from diverse areas of Malaysia, first analyzed separately and then in combination with previous data. All analyses confirmed 2 divergent clusters of human cases in Malaysian Borneo, associated with long-tailed macaques and pig-tailed macaques, and a third cluster in humans and most macaques in peninsular Malaysia. High levels of pairwise divergence between each of these sympatric and allopatric subpopulations have implications for the epidemiology and control of this zoonotic species.

BACKGROUND: Malaria control interventions have led to a decline in transmission intensity in many endemic areas, and resulted in elimination in some areas. This decline, however, will lead to delayed acquisition of protective immunity and thus impact disease manifestation and outcomes. Therefore, the variation in clinical and haematological parameters in children with malaria was assessed across three areas in Ghana with varying transmission intensities. METHODS: A total of 568 children between the ages of 2 and 14 years with confirmed malaria were recruited in hospitals in three areas with varying transmission intensities (Kintampo > Navrongo > Accra) and a comprehensive analysis of parasitological, clinical, haematological and socio-economic parameters was performed. RESULTS: Areas of lower malaria transmission tended to have lower disease severity in children with malaria, characterized by lower parasitaemias and higher haemoglobin levels. In addition, total white cell counts and percent lymphocytes decreased with decreasing transmission intensity. The heterozygous sickle haemoglobin genotype was protective against disease severity in Kintampo (P = 0.016), although this was not significant in Accra and Navrongo. Parasitaemia levels were not a significant predictor of haemoglobin level after controlling for age and gender. However, higher haemoglobin levels in children were associated with certain socioeconomic factors, such as having fathers who had any type of employment (P

Cultured human pathogens may differ significantly from source populations. To investigate the genetic basis of laboratory adaptation in malaria parasites, clinical Plasmodium falciparum isolates were sampled from patients and cultured in vitro for up to three months. Genome sequence analysis was performed on multiple culture time point samples from six monoclonal isolates, and single nucleotide polymorphism (SNP) variants emerging over time were detected. Out of a total of five positively selected SNPs, four represented nonsense mutations resulting in stop codons, three of these in a single ApiAP2 transcription factor gene, and one in SRPK1. To survey further for nonsense mutants associated with culture, genome sequences of eleven long-term laboratory-adapted parasite strains were examined, revealing four independently acquired nonsense mutations in two other ApiAP2 genes, and five in Epac. No mutants of these genes exist in a large database of parasite sequences from uncultured clinical samples. This implicates putative master regulator genes in which multiple independent stop codon mutations have convergently led to culture adaptation, affecting most laboratory lines of P. falciparum. Understanding the adaptive processes should guide development of experimental models, which could include targeted gene disruption to adapt fastidious malaria parasite species to culture.

Many human genetic associations with resistance to malaria have been reported, but few have been reliably replicated. We collected data on 11,890 cases of severe malaria due to Plasmodium falciparum and 17,441 controls from 12 locations in Africa, Asia and Oceania. We tested 55 SNPs in 27 loci previously reported to associate with severe malaria. There was evidence of association at P < 1 × 10(-4) with the HBB, ABO, ATP2B4, G6PD and CD40LG loci, but previously reported associations at 22 other loci did not replicate in the multicenter analysis. The large sample size made it possible to identify authentic genetic effects that are heterogeneous across populations or phenotypes, with a striking example being the main African form of G6PD deficiency, which reduced the risk of cerebral malaria but increased the risk of severe malarial anemia. The finding that G6PD deficiency has opposing effects on different fatal complications of P. falciparum infection indicates that the evolutionary origins of this common human genetic disorder are more complex than previously supposed.

The Anopheles gambiae complex of mosquitoes includes malaria vectors at different stages of speciation, whose study enables a better understanding of how adaptation to divergent environmental conditions leads to evolution of reproductive isolation. We investigated the population genetic structure of closely related sympatric taxa that have recently been proposed as separate species (An. coluzzii and An. gambiae), sampled from diverse habitats along the Gambia river in West Africa. We characterized putatively neutral microsatellite loci as well as chromosomal inversion polymorphisms known to be associated with ecological adaptation. The results revealed strong ecologically associated population subdivisions within both species. Microsatellite loci on chromosome-3L revealed clear differentiation between coastal and inland populations, which in An. coluzzii is reinforced by a unusual inversion polymorphism pattern, supporting the hypothesis of genetic divergence driven by adaptation to the coastal habitat. A strong reduction of gene flow was observed between An. gambiae populations west and east of an extensively rice-cultivated region apparently colonized exclusively by An. coluzzii. Notably, this 'intraspecific' differentiation is higher than that observed between the two species and involves also the centromeric region of chromosome-X which has previously been considered a marker of speciation within this complex, possibly suggesting that the two populations may be at an advanced stage of differentiation triggered by human-made habitat fragmentation. These results confirm ongoing ecological speciation within these most important Afro-tropical malaria vectors and raise new questions on the possible effect of this process in malaria transmission.

Killer-cell immunoglobulin-like receptors (KIRs) are a group of natural killer cell receptors (NKRs) that regulate NK-cell-mediated production of interferon gamma (IFN-γ) in response to infection. These receptors have recently been suggested to influence the severity of clinical Plasmodium falciparum malaria infection. We examined the KIR locus in relation to malaria in children from southwest Nigeria. Sequence specific priming (SSP)-PCR was used to detect the KIR genes. The presence or absence of fifteen different KIR genes was determined in each individual and the proportions compared across 3 clinical groups; asymptomatic malaria, uncomplicated clinical malaria and severe clinical malaria. The genes KIR2DL5, KIR2DS3 and KIR2DS5 were present in a significantly higher proportion of individuals in the asymptomatic control group than in the malaria cases. Furthermore, KIR2DS3 and KIR2DS5 were present in a higher proportion of uncomplicated malaria cases than severe malaria cases. Carriage c-AB2 genotype (which comprises all centromeric KIR genes including KIR2DL5, KIR2DS3 and KIR2DS5) decreases with severity of the disease suggesting that the KIR AB profile might be associated with protection from severe malaria infection in this population in Nigeria.

INTRODUCTION Malaria parasites cause human disease by invading and replicating inside red blood cells. In the case of Plasmodium falciparum , this can lead to severe forms of malaria that are a major cause of childhood mortality in Africa. This species of parasite enters the red blood cell through interactions with surface proteins including the glycophorins GYPA and GYPB, which determine the polymorphic MNS blood group system. In a recent genome-wide association study, we identified alleles associated with protection against severe malaria near the cluster of genes encoding these invasion receptors. RATIONALE Investigation of genetic variants at this locus and their relation to severe malaria is challenging because of the high sequence similarity between the neighboring glycophorin genes and the relative lack of available sequence data capturing the genetic diversity of sub-Saharan Africa. To better assess whether variation in the glycophorin genes could explain the signal of association, we generated additional sequence data from sub-Saharan African populations and developed an analytical approach to characterize structural variation at this complex locus. RESULTS Using 765 newly sequenced human genomes from 10 African ethnic groups along with data from the 1000 Genomes Project, we generated a reference panel of haplotypes across the glycophorin region. In addition to single-nucleotide polymorphisms and short indels, we assayed large copy number variants (CNVs) using sequencing read depth and uncovered extensive structural diversity. By imputing from this reference panel into 4579 severe malaria cases and 5310 controls from three African populations, we found that a complex CNV, here called DUP4, is associated with resistance to severe malaria and fully explains the previously reported signal of association. In our sample, DUP4 is present only in east Africa, and this localization, as well as the extent of similarity between DUP4 haplotypes, suggests that it has recently increased in frequency, presumably under natural selection due to malaria. To evaluate the potential functional consequences of this structural variant, we analyzed high-coverage sequence-read data from multiple individuals to generate a model of the DUP4 chromosome structure. The DUP4 haplotype contains five glycophorin genes, including two hybrid genes that juxtapose the extracellular domain of GYPB with the transmembrane and intracellular domains of GYPA. Noting that these predicted hybrids are characteristic of the Dantu antigen in the MNS blood group system, we sequenced a Dantu positive individual and confirmed that DUP4 is the molecular basis of the Dantu NE blood group variant. CONCLUSION Although a role for GYPA and GYPB in parasite invasion is well known, a direct link between glycophorin polymorphisms and clinical susceptibility to malaria has been elusive. Here we have provided a systematic catalog of CNVs, describing structural diversity that may have functional importance at this locus. Our results identify a specific variant that encodes hybrid glycophorin proteins and is associated with protection against severe malaria. This discovery calls for further work to determine how this particular molecular rearrangement affects parasite invasion and the red blood cell response and may lead us toward new parasite vulnerabilities that can be utilized in future interventions against this deadly disease.

Prospective studies continue to identify malaria parasite genes with particular patterns of polymorphism which indicate they may be under immune selection, and the encoded proteins require investigation. Sixteen new recombinant protein reagents were designed to characterize three such polymorphic proteins expressed in Plasmodium falciparum schizonts and merozoites: MSPDBL1 (also termed MSP3.4) and MSPDBL2 (MSP3.8), which possess Duffy binding-like (DBL) domains, and SURFIN4.2, encoded by a member of the surface-associated interspersed ( surf ) multigene family. After testing the antigenicities of these reagents by murine immunization and parasite immunofluorescence, we analyzed naturally acquired antibody responses to the antigens in two cohorts in coastal Kenya in which the parasite was endemic (Chonyi [ n = 497] and Ngerenya [ n = 461]). As expected, the prevalence and levels of serum antibodies increased with age. We then investigated correlations with subsequent risk of clinical malaria among children

Understanding genetic causes and effects of speciation in sympatric populations of sexually reproducing eukaryotes is challenging, controversial, and of practical importance for controlling rapidly evolving pests and pathogens. The major African malaria vector mosquito Anopheles gambiae sensu stricto (s.s.) is considered to contain two incipient species with strong reproductive isolation, hybrids between the M and S molecular forms being very rare. Following recent observations of higher proportions of hybrid forms at a few sites in West Africa, we conducted new surveys of 12 sites in four contiguous countries (The Gambia, Senegal, Guinea-Bissau, and Republic of Guinea). Identification and genotyping of 3499 A. gambiae s.s. revealed high frequencies of M/S hybrid forms at each site, ranging from 5 to 42%, and a large spectrum of inbreeding coefficient values from 0.11 to 0.76, spanning most of the range expected between the alternative extremes of panmixia and assortative mating. Year-round sampling over 2 years at one of the sites in The Gambia showed that M/S hybrid forms had similar relative frequencies throughout periods of marked seasonal variation in mosquito breeding and abundance. Genome-wide scans with an Affymetrix high-density single-nucleotide polymorphism (SNP) microarray enabled replicate comparisons of pools of different molecular forms, in three separate populations. These showed strong differentiation between M and S forms only in the pericentromeric region of the X chromosome that contains the molecular form-specific marker locus, with only a few other loci showing minor differences. In the X chromosome, the M/S hybrid forms were more differentiated from M than from S forms, supporting a hypothesis of asymmetric introgression and backcrossing.

Chronic infections with the human malaria parasite Plasmodium falciparum depend on antigenic variation. P. falciparum erythrocyte membrane protein 1 (PfEMP1), the major erythrocyte surface antigen mediating parasite sequestration in the microvasculature, is encoded in parasites by a highly diverse family of var genes. Antigenic switching is mediated by clonal variation in var expression, and recent in vitro studies have demonstrated a role for epigenetic processes in var regulation. Expression of particular PfEMP1 variants may result in parasite enrichment in different tissues, a factor in the development of severe disease. Here, we study in vivo human infections and provide evidence that infection-induced stress responses in the host can modify PfEMP1 expression via the perturbation of epigenetic mechanisms. Our work suggests that severe disease may not be the direct result of an adaptive virulence strategy to maximize parasite survival but that it may indicate a loss of control of the carefully regulated process of antigenic switching that maintains chronic infections.

With the paradigm shift from the reduction of morbidity and mortality to the interruption of transmission, the focus of malaria control broadens from symptomatic infections in children ≤ 5 years of age to include asymptomatic infections in older children and adults. In addition, as control efforts intensify and the number of interventions increases, there will be decreases in prevalence, incidence and transmission with additional decreases in morbidity and mortality. Expected secondary consequences of these changes include upward shifts in the peak ages for infection (parasitemia) and disease, increases in the ages for acquisition of antiparasite humoral and cellular immune responses and increases in false-negative blood smears and rapid diagnostic tests. Strategies to monitor these changes must include: 1] studies of the entire population (that are not restricted to children ≤ 5 or ≤ 10 years of age), 2] study sites in both cities and rural areas (because of increasing urbanization across sub-Saharan Africa) and 3] innovative strategies for surveillance as the prevalence of infection decreases and the frequency of false-negative smears and rapid diagnostic tests increases.

Merozoite surface protein 1 (MSP1) of Plasmodium falciparum has been implicated as an important target of acquired immunity, and candidate components for a vaccine include polymorphic epitopes in the N-terminal polymorphic block 2 region. We designed a polyvalent hybrid recombinant protein incorporating sequences of the three major allelic types of block 2 together with a composite repeat sequence of one of the types and N-terminal flanking T cell epitopes, and compared this with a series of recombinant proteins containing modular sub-components and similarly expressed in Escherichia coli. Immunogenicity of the full polyvalent hybrid protein was tested in both mice and rabbits, and comparative immunogenicity studies of the sub-component modules were performed in mice. The full hybrid protein induced high titre antibodies against each of the major block 2 allelic types expressed as separate recombinant proteins and against a wide range of allelic types naturally expressed by a panel of diverse P. falciparum isolates, while the sub-component modules had partial antigenic coverage as expected. This encourages further development and evaluation of the full MSP1 block 2 polyvalent hybrid protein as a candidate blood-stage component of a malaria vaccine.

Parasite antigen genes might evolve under frequency-dependent immune selection. The distinctive patterns of polymorphism that result can be detected using population genetic methods that test for signatures of balancing selection, allowing genes encoding important targets of immunity to be identified. Analyses can be complicated by population structures, histories and features of a parasite's genome. However, new sequencing technologies facilitate scans of polymorphism throughout parasite genomes to identify the most exceptional gene specific signatures. We focus on malaria parasites to illustrate challenges and opportunities for detecting targets of frequency-dependent immune selection to discover new potential vaccine candidates.

Signatures of balancing selection operating on specific gene loci in endemic pathogens can identify candidate targets of naturally acquired immunity. In malaria parasites, several leading vaccine candidates convincingly show such signatures when subjected to several tests of neutrality, but the discovery of new targets affected by selection to a similar extent has been slow. A small minority of all genes are under such selection, as indicated by a recent study of 26 Plasmodium falciparum merozoite-stage genes that were not previously prioritized as vaccine candidates, of which only one (locus PF10_0348) showed a strong signature. Therefore, to focus discovery efforts on genes that are polymorphic, we scanned all available shotgun genome sequence data from laboratory lines of P. falciparum and chose six loci with more than five single nucleotide polymorphisms per kilobase (including PF10_0348) for in-depth frequency-based analyses in a Kenyan population (allele sample sizes >50 for each locus) and comparison of Hudson-Kreitman-Aguade (HKA) ratios of population diversity (π) to interspecific divergence (K) from the chimpanzee parasite Plasmodium reichenowi. Three of these (the msp3/6-like genes PF10_0348 and PF10_0355 and the surf(4.1) gene PFD1160w) showed exceptionally high positive values of Tajima's D and Fu and Li's F indices and have the highest HKA ratios, indicating that they are under balancing selection and should be prioritized for studies of their protein products as candidate targets of immunity. Combined with earlier results, there is now strong evidence that high HKA ratio (as well as the frequency-independent ratio of Watterson's /K) is predictive of high values of Tajima's D. Thus, the former offers value for use in genome-wide screening when numbers of genome sequences within a species are low or in combination with Tajima's D as a 2D test on large population genomic samples.

SUMMARY Longitudinal cohort studies are important to describe the dynamics of naturally acquired antibody response profiles to defined Plasmodium falciparum malaria antigens relative to clinical malaria episodes. In children under 7 years of age in The Gambia, serum IgG responses were measured to P. falciparum merozoite antigens AMA1, EBA175, MSP119, MSP2 and crude schizont extract, over a 10-month period. Persistence of antibody responses was measured in 152 children during the dry season when there was virtually no malaria transmission, and 103 children were monitored for new episodes of clinical malaria during the subsequent wet season when transmission occurred. Children who experienced clinical malaria had lower antibody levels at the start of the study than those who remained free from malaria. Associations between dry season antibody persistence and subsequent wet season antibody levels suggested robust immunological memory responses. Mean antibody levels to all antigens were elevated by the end of the wet season in children who experienced clinical malaria; each of these children had a boosted antibody response to at least one antigen. In all children, antibody avidities were lower against MSP2 than other antigens, a difference that did not change throughout the study period or in relation to clinical malaria episodes.

Larviciding to control malaria was assessed in rural areas with extensive seasonal flooding. Larval and adult mosquitoes and malaria incidence were surveyed routinely in four 100-km(2) areas either side of the Gambia River. Baseline data were collected in 2005. Microbial larvicide was applied to all water bodies by hand application with water-dispersible granular formulations and corn granules weekly from May to November in two areas in 2006 and in the other two areas in 2007 in a cross-over design. The intervention was associated with a reduction in habitats with late stage anopheline larvae and an 88% reduction in larval densities (P < 0.001). The effect of the intervention on mosquito densities was not pronounced and was confounded by the distance of villages to the major breeding sites and year (P = 0.002). There was no reduction in clinical malaria or anemia. Ground applications of non-residual larvicides with simple equipment are not effective in riverine areas with extensive flooding, where many habitats are poorly demarcated, highly mobile, and inaccessible on foot.