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.

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.

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed. Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Background Plasmodium vivax has been recently discovered as a significant cause of malaria in Mauritania, although very rare elsewhere in West Africa. It has not been known if this is a recently introduced or locally remnant parasite population, nor whether the genetic structure reflects epidemic or endemic transmission. Methodology/Principal findings To investigate the P. vivax population genetic structure in Mauritania and compare with populations previously analysed elsewhere, multi-locus genotyping was undertaken on 100 clinical isolates, using a genome-wide panel of 38 single nucleotide polymorphisms (SNPs), plus seven SNPs in drug resistance genes. The Mauritanian P. vivax population is shown to be genetically diverse and divergent from populations elsewhere, indicated consistently by genetic distance matrix analysis, principal components analyses, and fixation indices. Only one isolate had a genotype clearly indicating recent importation, from a southeast Asian source. There was no linkage disequilibrium in the local parasite population, and only a small number of infections appeared to be closely genetically related, indicating that there is ongoing genetic recombination consistent with endemic transmission. The P. vivax diversity in a remote mining town was similar to that in the capital Nouakchott, with no indication of local substructure or of epidemic population structure. Drug resistance alleles were virtually absent in Mauritania, in contrast with P. vivax in other areas of the world. Conclusions/Significance The molecular epidemiology indicates that there is long-standing endemic transmission that will be very challenging to eliminate. The virtual absence of drug resistance alleles suggests that most infections have been untreated, and that this endemic infection has been more neglected in comparison to P. vivax elsewhere., Author summary Plasmodium vivax is a widespread cause of malaria in Mauritania, in contrast to its rarity elsewhere throughout West Africa. To investigate whether the parasite may be recently introduced or epidemic, multi-locus genotyping was performed on 100 Mauritanian P. vivax malaria cases. Analysis of a genome-wide panel of single nucleotide polymorphisms showed the P. vivax population to be genetically diverse and divergent from populations elsewhere, indicating that there has been long-standing endemic transmission. Almost all infections appear to be locally acquired, with the exception of one that was presumably imported with a genotype similar to infections seen in Southeast Asia. The Mauritanian P. vivax population shows no linkage disequilibrium, and very few infections have closely related genotypes, indicating ongoing recombination. The parasite showed no indication of local substructure or epidemic population structure. Drug resistance alleles were virtually absent, suggesting that most infections have been untreated historically. The molecular epidemiology indicates that there has been long-standing endemic transmission of this neglected parasite that requires special attention for control.

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.

Parasites infect hosts in widely varying environments, encountering diverse challenges for adaptation. To identify malaria parasite genes under locally divergent selection across a large endemic region with a wide spectrum of transmission intensity, genome sequences were obtained from 284 clinical Plasmodium falciparum infections from four newly sampled locations in Senegal, The Gambia, Mali and Guinea. Combining these with previous data from seven other sites in West Africa enabled a multi-population analysis to identify discrete loci under varying local selection. A genome-wide scan showed the most exceptional geographical divergence to be at the early gametocyte gene locus gdv1 which is essential for parasite sexual development and transmission. We identified a major structural dimorphism with alternative 1.5 kb and 1.0 kb sequence deletions at different positions of the 3′-intergenic region, in tight linkage disequilibrium with the most highly differentiated single nucleotide polymorphism, one of the alleles being very frequent in Senegal and The Gambia but rare in the other locations. Long non-coding RNA transcripts were previously shown to include the entire antisense of the gdv1 coding sequence and the portion of the intergenic region with allelic deletions, suggesting adaptive regulation of parasite sexual development and transmission in response to local conditions.

Genome sequences of 247 Plasmodium falciparum isolates collected in The Gambia in 2008 and 2014 were analysed to identify changes possibly related to the scale-up of antimalarial interventions that occurred during this period. Overall, there were 15 regions across the genomes with signatures of positive selection. Five of these were sweeps around known drug resistance and antigenic loci. Signatures at antigenic loci such as thrombospodin related adhesive protein (Pftrap) were most frequent in eastern Gambia, where parasite prevalence and transmission remain high. There was a strong temporal differentiation at a non-synonymous SNP in a cysteine desulfarase (Pfnfs) involved in iron-sulphur complex biogenesis. During the 7-year period, the frequency of the lysine variant at codon 65 (Pfnfs-Q65K) increased by 22% (10% to 32%) in the Greater Banjul area. Between 2014 and 2015, the frequency of this variant increased by 6% (20% to 26%) in eastern Gambia. IC50 for lumefantrine was significantly higher in Pfnfs-65K isolates. This is probably the first evidence of directional selection on Pfnfs or linked loci by lumefantrine. Given the declining malaria transmission, the consequent loss of population immunity, and sustained drug pressure, it is important to monitor Gambian P. falciparum populations for further signs of adaptation.

BackgroundPlasmodium falciparum uses a repertoire of merozoite-stage proteins for invasion of erythrocytes. Antibodies against some of these proteins halt the replication cycle of the parasite by preventing erythrocyte invasion and are implicated as contributors to protective immunity against malaria.MethodsWe assayed antibody reactivity against a panel of 9 recombinant antigens based on erythrocyte-binding antigen (EBA) and reticulocyte-like homolog (Rh) proteins in plasma from children with malaria and healthy adults residing in 3 endemic areas in Ghana using enzyme-linked immunosorbent assay. Purified immunoglobulin (Ig)G from adult plasma samples was also tested for invasion inhibition against 7 different P falciparum culture lines, including clinical isolates.ResultsAntibodies against the antigens increased in an age-dependent manner in children. Breadth of reactivity to the different antigens was strongly associated with in vitro parasite growth inhibitory activity of IgG purified from the adults. The strongest predictors of breadth of antibody reactivity were age and transmission intensity, and a combination of reactivities to Rh2, Rh4, and Rh5 correlated strongly with invasion inhibition.ConclusionsGrowth inhibitory activity was significantly associated with breadth of antibody reactivity to merozoite antigens, encouraging the prospect of a multicomponent blood-stage vaccine.

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.

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.

Background Plasmodium falciparum may evade complement-mediated host defense by hijacking complement Factor H (FH), a negative regulator of the alternative complement pathway. Plasma levels of FH vary between individuals and may therefore influence malaria susceptibility and severity. Methods We measured convalescent FH plasma levels in 149 Gambian children who had recovered from uncomplicated or severe P. falciparum malaria and in 173 healthy control children. We compared FH plasma levels between children with malaria and healthy controls, and between children with severe (n = 82) and uncomplicated malaria (n = 67). We determined associations between FH plasma levels and laboratory features of severity and used multivariate analyses to examine associations with FH when accounting for other determinants of severity. Results FH plasma levels differed significantly between controls, uncomplicated malaria cases, and severe malaria cases (mean [95% confidence interval], 257 [250 to 264], 288 [268 to 309], and 328 [313 to 344] µg/mL, respectively; analysis of variance P < .0001). FH plasma levels correlated with severity biomarkers, including lactate, parasitemia, and parasite density, but did not correlate with levels of PfHRP2, which represent the total body parasite load. Associations with severity and lactate remained significant when adjusting for age and parasite load. Conclusions Natural variation in FH plasma levels is associated with malaria susceptibility and severity. A prospective study will be needed to strengthen evidence for causation, but our findings suggest that interfering with FH binding by P. falciparum might be useful for malaria prevention or treatment.

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.

Improved methods are needed to identify host mechanisms which directly protect against human infectious diseases in order to develop better vaccines and therapeutics1,2. Pathogen load determines the outcome of many infections3, and is a consequence of pathogen multiplication rate, duration of the infection, and inhibition or killing of pathogen by the host (resistance). If these determinants of pathogen load could be quantified then their mechanistic correlates might be determined. In humans the timing of infection is rarely known and treatment cannot usually be withheld to monitor serial changes in pathogen load and host response. Here we present an approach to overcome this and identify potential mechanisms of resistance which control parasite load inPlasmodium falciparummalaria. Using a mathematical model of longitudinal infection dynamics for orientation, we made individualized estimates of parasite multiplication and growth inhibition in Gambian children at presentation with acute malaria and used whole blood RNA-sequencing to identify their correlates. We identified novel roles for secreted proteases cathepsin G and matrix metallopeptidase 9 (MMP9) as direct effector molecules which inhibitP. falciparumgrowth. Cathepsin G acts on the erythrocyte membrane, cleaving surface receptors required for parasite invasion, whilst MMP9 acts on the parasite. In contrast, the type 1 interferon response and expression ofCXCL10(IFN-γ-inducible protein of 10 kDa, IP-10) were detrimental to control of parasite growth. Natural variation in iron status and plasma levels of complement factor H were determinants of parasite multiplication rate. Our findings demonstrate the importance of accounting for the dynamic interaction between host and pathogen when seeking to identify correlates of protection, and reveal novel mechanisms controlling parasite growth in humans. This approach could be extended to identify additional mechanistic correlates of natural- and vaccine-induced immunity to malaria and other infections.

Background: Although thousands of clinical isolates of Plasmodium falciparum are being sequenced and analysed by short read technology, the data do not resolve the highly variable subtelomeric regions of the genomes that contain polymorphic gene families involved in immune evasion and pathogenesis. There is also no current standard definition of the boundaries of these variable subtelomeric regions. Methods: Using long-read sequence data (Pacific Biosciences SMRT technology), we assembled and annotated the genomes of 15 P. falciparum isolates, ten of which are newly cultured clinical isolates. We performed comparative analysis of the entire genome with particular emphasis on the subtelomeric regions and the internal var genes clusters. Results: The nearly complete sequence of these 15 isolates has enabled us to define a highly conserved core genome, to delineate the boundaries of the subtelomeric regions, and to compare these across isolates. We found highly structured variable regions in the genome. Some exported gene families purportedly involved in release of merozoites show copy number variation. As an example of ongoing genome evolution, we found a novel CLAG gene in six isolates. We also found a novel gene that was relatively enriched in the South East Asian isolates compared to those from Africa. Conclusions: These 15 manually curated new reference genome sequences with their nearly complete subtelomeric regions and fully assembled genes are an important new resource for the malaria research community. We report the overall conserved structure and pattern of important gene families and the more clearly defined subtelomeric regions.