Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana.

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. In low transmission settings, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. We performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n = 1,409) through estimation of identity-by-descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 multi-isolate clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally observed 61 nonsynonymous substitutions that increased markedly in frequency over the study period as well as a novel pfk13 mutation (G718S). However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions, given that clones carrying drug resistance polymorphisms do not demonstrate enhanced persistence or higher abundance than clones carrying polymorphisms of comparable frequency that are unrelated to resistance. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Genotypic analysis of RTS,S/AS01E malaria vaccine efficacy against parasite infection as a function of dosage regimen and baseline malaria infection status in children aged 5-17 months in Ghana and Kenya: a longitudinal phase 2b randomised controlled trial.

BACKGROUND: The first licensed malaria vaccine, RTS,S/AS01E, confers moderate protection against symptomatic disease. Because many malaria infections are asymptomatic, we conducted a large-scale longitudinal parasite genotyping study of samples from a clinical trial exploring how vaccine dosing regimen affects vaccine efficacy. METHODS: Between Sept 28, 2017, and Sept 25, 2018, 1500 children aged 5-17 months were randomly assigned (1:1:1:1:1) to receive four different RTS,S/AS01E regimens or a rabies control vaccine in a phase 2b open-label clinical trial in Ghana and Kenya. Participants in the four RTS,S groups received two full doses at month 0 and month 1 and either full doses at month 2 and month 20 (group R012-20); full doses at month 2, month 14, month 26, and month 38 (group R012-14); fractional doses at month 2, month 14, month 26, and month 38 (group Fx012-14; early fourth dose); or fractional doses at month 7, month 20, and month 32 (group Fx017-20; delayed third dose). We evaluated the time to the first new genotypically detected infection and the total number of new infections during two follow-up periods (12 months and 20 months) in more than 36 000 dried blood spot specimens from 1500 participants. To study vaccine effects on time to the first new infection, we defined vaccine efficacy as one minus the hazard ratio (HR; RTS,S vs control) of the first new infection. We performed a post-hoc analysis of vaccine efficacy based on malaria infection status at first vaccination and force of infection by month 2. This trial (MAL-095) is registered with ClinicalTrials.gov, NCT03281291. FINDINGS: We observed significant and similar vaccine efficacy (25-43%; 95% CI union 9-53) against first new infection for all four RTS,S/AS01E regimens across both follow-up periods (12 months and 20 months). Each RTS,S/AS01E regimen significantly reduced the mean number of new infections in the 20-month follow-up period by 1·1-1·6 infections (95% CI union 0·6-2·1). Vaccine efficacy against first new infection was significantly higher in participants who were infected with malaria (68%; 95% CI 50-80) than in those who were uninfected (37%; 23-48) at the first vaccination (p=0·0053). INTERPRETATION: All tested dosing regimens blocked some infections to a similar degree. Improved vaccine efficacy in participants infected during vaccination could suggest new strategies for highly efficacious malaria vaccine development and implementation. FUNDING: GlaxoSmithKline Biologicals SA, PATH, Bill & Melinda Gates Foundation, and the German Federal Ministry of Education and Research.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana.

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. However, in low transmission settings where most mosquitoes become infected with only a single parasite clone, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. To investigate whether this clonality was potentially associated with the persistence and spatial spread of the mutation, we performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n=1,409) through estimation of identity by descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally report polymorphisms exhibiting evidence of selection for drug resistance or other phenotypes and reported a novel pfk13 mutation (G718S) as well as 61 nonsynonymous substitutions that increased markedly in frequency. However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Impact of piperaquine resistance in Plasmodium falciparum on malaria treatment effectiveness in The Guianas: a descriptive epidemiological study.

BACKGROUND: Plasmodium falciparum is an apicomplexan parasite responsible for lethal cases of malaria. According to WHO recommendations, P falciparum cases are treated with artemisinin-based combination therapy including dihydroartemisinin-piperaquine. However, the emergence of resistant parasites against dihydroartemisinin-piperaquine was reported in southeast Asia in 2008 and, a few years later, suspected in South America. METHODS: To characterise resistance emergence, a treatment efficacy study was performed on the reported patients infected with P falciparum and treated with dihydroartemisinin-piperaquine in French Guiana (n=6, 2016-18). Contemporary isolates collected in French Guiana were genotyped for P falciparum chloroquine resistance transporter (pfCRT; n=845) and pfpm2 and pfpm3 copy number (n=231), phenotyped using the in vitro piperaquine survival assay (n=86), and analysed through genomic studies (n=50). Additional samples from five Amazonian countries and one outside the region were genotyped (n=1440). FINDINGS: In field isolates, 40 (47%) of 86 (95% CI 35·9-57·1) were resistant to piperaquine in vitro; these phenotypes were more associated with pfCRTC350R (ie, Cys350Arg) and pfpm2 and pfpm3 amplifications (Dunn test, p<0·001). Those markers were also associated with dihydroartemisinin-piperaquine treatment failure (n=3 [50%] of 6). A high prevalence of piperaquine resistance markers was observed in Suriname in 19 (83%) of 35 isolates and in Guyana in 579 (73%) of 791 isolates. The pfCRTC350R mutation emerged before pfpm2 and pfpm3 amplification in a temporal sequence different from southeast Asia, and in the absence of artemisinin partial resistance, suggesting a geographically distinctive epistatic relationship between these genetic markers. INTERPRETATION: The high prevalence of piperaquine resistance markers in parasite populations of the Guianas, and the risk of associated therapeutic failures calls for caution on dihydroartemisinin-piperaquine use in the region. Furthermore, greater attention should be given to potential differences in genotype to phenotype mapping across genetically distinct parasite populations from different continents. FUNDING: Pan American Health Organization and WHO, French Ministry for Research, European Commission, Santé publique France, Agence Nationale de la Recherche, Fundação de Amparo à Pesquisa do Estado do Amazonas, Ministry of Health of Brazil, Oswaldo Cruz Foundation, and National Institutes of Health. TRANSLATIONS: For the French and Portuguese translations of the abstract see Supplementary Materials section.

Baseline malaria infection status and RTS,S/AS01E malaria vaccine efficacy.

Background: The only licensed malaria vaccine, RTS,S/AS01 E , confers moderate protection against symptomatic disease. Because many malaria infections are asymptomatic, we conducted a large-scale longitudinal parasite genotyping study of samples from a clinical trial exploring how vaccine dosing regimen affects vaccine efficacy (VE). Methods: 1,500 children aged 5-17 months were randomized to receive four different RTS,S/AS01 E regimens or a rabies control vaccine in a phase 2b clinical trial in Ghana and Kenya. We evaluated the time to the first new genotypically detected infection and the total number of new infections during two follow-up periods in over 36K participant specimens. We performed a post hoc analysis of VE based on malaria infection status at first vaccination and force of infection. Results: We observed significant and comparable VE (25-43%, 95% CI union 9-53%) against first new infection for all four RTS,S/AS01 E regimens across both follow-up periods (12 and 20 months). Each RTS,S/AS01 E regimen significantly reduced the number of new infections in the 20-month follow-up period (control mean 4.1 vs. RTS,S/AS01 E mean 2.6-3.0). VE against first new infection was significantly higher in participants who were malaria-infected (68%; 95% CI, 50 to 80%) versus uninfected (37%; 95% CI, 23 to 48%) at the first vaccination (P=0.0053) and in participants experiencing greater force of infection between dose 1 and 3 (P=0.059). Conclusions: All tested dosing regimens blocked some infections to a similar degree. Improved VE in participants infected during vaccination could suggest new strategies for highly efficacious malaria vaccine development and implementation. ( ClinicalTrials.gov number, NCT03276962 ).

Potent acyl-CoA synthetase 10 inhibitors kill Plasmodium falciparum by disrupting triglyceride formation.

Identifying how small molecules act to kill malaria parasites can lead to new "chemically validated" targets. By pressuring Plasmodium falciparum asexual blood stage parasites with three novel structurally-unrelated antimalarial compounds (MMV665924, MMV019719 and MMV897615), and performing whole-genome sequence analysis on resistant parasite lines, we identify multiple mutations in the P. falciparum acyl-CoA synthetase (ACS) genes PfACS10 (PF3D7_0525100, M300I, A268D/V, F427L) and PfACS11 (PF3D7_1238800, F387V, D648Y, and E668K). Allelic replacement and thermal proteome profiling validates PfACS10 as a target of these compounds. We demonstrate that this protein is essential for parasite growth by conditional knockdown and observe increased compound susceptibility upon reduced expression. Inhibition of PfACS10 leads to a reduction in triacylglycerols and a buildup of its lipid precursors, providing key insights into its function. Analysis of the PfACS11 gene and its mutations point to a role in mediating resistance via decreased protein stability.

Resolving drug selection and migration in an inbred South American Plasmodium falciparum population with identity-by-descent analysis.

The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination.

Declines in prevalence alter the optimal level of sexual investment for the malaria parasite Plasmodium falciparum.

Successful infectious disease interventions can result in large reductions in parasite prevalence. Such demographic change has fitness implications for individual parasites and may shift the parasite's optimal life history strategy. Here, we explore whether declining infection rates can alter Plasmodium falciparum's investment in sexual versus asexual growth. Using a multiscale mathematical model, we demonstrate how the proportion of polyclonal infections, which decreases as parasite prevalence declines, affects the optimal sexual development strategy: Within-host competition in multiclone infections favors a greater investment in asexual growth whereas single-clone infections benefit from higher conversion to sexual forms. At the same time, drug treatment also imposes selection pressure on sexual development by shortening infection length and reducing within-host competition. We assess these models using 148 P. falciparum parasite genomes sampled in French Guiana over an 18-y period of intensive intervention (1998 to 2015). During this time frame, multiple public health measures, including the introduction of new drugs and expanded rapid diagnostic testing, were implemented, reducing P. falciparum malaria cases by an order of magnitude. Consistent with this prevalence decline, we see an increase in the relatedness among parasites, but no single clonal background grew to dominate the population. Analyzing individual allele frequency trajectories, we identify genes that likely experienced selective sweeps. Supporting our model predictions, genes showing the strongest signatures of selection include transcription factors involved in the development of P. falciparum's sexual gametocyte form. These results highlight how public health interventions impose wide-ranging selection pressures that affect basic parasite life history traits.

Design and implementation of multiplexed amplicon sequencing panels to serve genomic epidemiology of infectious disease: A malaria case study.

Multiplexed PCR amplicon sequencing (AmpSeq) is an increasingly popular application for cost-effective monitoring of threatened species and managed wildlife populations, and shows strong potential for the genomic epidemiology of infectious disease. AmpSeq data from infectious microbes can inform disease control in multiple ways, such as by measuring drug resistance marker prevalence, distinguishing imported from local cases, and determining the effectiveness of therapeutics. We describe the design and comparative evaluation of two new AmpSeq assays for Plasmodium falciparum malaria parasites: a four-locus panel ("4CAST") composed of highly diverse antigens, and a 129-locus panel ("AMPLseq") composed of drug resistance markers, highly diverse loci for inferring relatedness, and a locus to detect Plasmodium vivax co-infection. We explore the performance of each panel in various public health use cases with in silico simulations as well as empirical experiments. The 4CAST panel appears highly suitable for evaluating the number of distinct parasite strains within samples (complexity of infection), showing strong performance across a wide range of parasitaemia levels without a DNA pre-amplification step. For relatedness inference, the larger AMPLseq panel performs similarly to two existing panels of comparable size, despite differences in the data and approach used for designing each panel. Finally, we describe an R package (paneljudge) that facilitates the design and comparative evaluation of genetic panels for relatedness estimation, and we provide general guidance on the design and implementation of AmpSeq panels for the genomic epidemiology of infectious disease.

A novel CSP C-terminal epitope targeted by an antibody with protective activity against Plasmodium falciparum.

Potent and durable vaccine responses will be required for control of malaria caused by Plasmodium falciparum (Pf). RTS,S/AS01 is the first, and to date, the only vaccine that has demonstrated significant reduction of clinical and severe malaria in endemic cohorts in Phase 3 trials. Although the vaccine is protective, efficacy declines over time with kinetics paralleling the decline in antibody responses to the Pf circumsporozoite protein (PfCSP). Although most attention has focused on antibodies to repeat motifs on PfCSP, antibodies to other regions may play a role in protection. Here, we expressed and characterized seven monoclonal antibodies to the C-terminal domain of CSP (ctCSP) from volunteers immunized with RTS,S/AS01. Competition and crystal structure studies indicated that the antibodies target two different sites on opposite faces of ctCSP. One site contains a polymorphic region (denoted α-ctCSP) and has been previously characterized, whereas the second is a previously undescribed site on the conserved ß-sheet face of the ctCSP (denoted ß-ctCSP). Antibodies to the ß-ctCSP site exhibited broad reactivity with a diverse panel of ctCSP peptides whose sequences were derived from field isolates of P. falciparum whereas antibodies to the α-ctCSP site showed very limited cross reactivity. Importantly, an antibody to the ß-site demonstrated inhibition activity against malaria infection in a murine model. This study identifies a previously unidentified conserved epitope on CSP that could be targeted by prophylactic antibodies and exploited in structure-based vaccine design.

Multiplex serological assay for establishing serological profiles of polymorphic, closely related peptide antigens.

Profiling of serological responses to establish the landscape of antibody specificities in individuals exposed to pathogens or vaccines is crucial for (a) revealing humoral immune correlates of protection; (b) uncovering markers of pathogen exposure; and (c) identifying antigens and epitopes associated with disease vs. protection. Establishing the antigenic profile of serological responses requires either expensive microarrays or labor- and time-intensive ELISA assays. Multiplex assay platforms are increasingly being evaluated for their usefulness for high-throughput testing of sera or plasma. The methodology described here utilizes a plate-based assay that allows the simultaneous detection of up to ten antigens per well in a 96 well format using an electrochemiluminescence immunoassay (ECLIA).•The newly developed protocol outlines high-throughput profiling of serological responses using a multiplex testing platform with subsequent computational analysis.•The protocol is a modification of the basic assay development manual from the manufacturer of the MESO QuickPlex SQ 120 instrument (MSD, Gaithersburg, MD) and can be used for synthetic peptides as well as full length proteins.•The protocol can be applied to map serological responses to pathogens or pathogen-derived antigens to establish serological profiles in search for biomarkers or immune correlates.

A newly characterized malaria antigen on erythrocyte and merozoite surfaces induces parasite inhibitory antibodies.

We previously identified a Plasmodium falciparum (Pf) protein of unknown function encoded by a single-copy gene, PF3D7_1134300, as a target of antibodies in plasma of Tanzanian children in a whole-proteome differential screen. Here we characterize this protein as a blood-stage antigen that localizes to the surface membranes of both parasitized erythrocytes and merozoites, hence its designation as Pf erythrocyte membrane and merozoite antigen 1 (PfEMMA1). Mouse anti-PfEMMA1 antisera and affinity-purified human anti-PfEMMA1 antibodies inhibited growth of P. falciparum strains by up to 68% in growth inhibition assays. Following challenge with uniformly fatal Plasmodium berghei (Pb) ANKA, up to 40% of mice immunized with recombinant PbEMMA1 self-cured, and median survival of lethally infected mice was up to 2.6-fold longer than controls (21 vs. 8 d, P = 0.005). Furthermore, high levels of naturally acquired human anti-PfEMMA1 antibodies were associated with a 46% decrease in parasitemia over 2.5 yr of follow-up of Tanzanian children. Together, these findings suggest that antibodies to PfEMMA1 mediate protection against malaria.

Plasmodium simium: Population Genomics Reveals the Origin of a Reverse Zoonosis.

BACKGROUND: The population history of Plasmodium simium, which causes malaria in sylvatic Neotropical monkeys and humans along the Atlantic Coast of Brazil, remains disputed. Genetically diverse P vivax populations from various sources, including the lineages that founded the species P simium, are thought to have arrived in the Americas in separate migratory waves. METHODS: We use population genomic approaches to investigate the origin and evolution of P simium. RESULTS: We find a minimal genome-level differentiation between P simium and present-day New World P vivax isolates, consistent with their common geographic origin and subsequent divergence on this continent. The meagre genetic diversity in P simium samples from humans and monkeys implies a recent transfer from humans to non-human primates - a unique example of malaria as a reverse zoonosis of public health significance. Likely genomic signatures of P simium adaptation to new hosts include the deletion of >40% of a key erythrocyte invasion ligand, PvRBP2a, which may have favored more efficient simian host cell infection. CONCLUSIONS: New World P vivax lineages that switched from humans to platyrrhine monkeys founded the P simium population that infects nonhuman primates and feeds sustained human malaria transmission in the outskirts of major cities.

Dense time-course gene expression profiling of the Drosophila melanogaster innate immune response.

BACKGROUND: Immune responses need to be initiated rapidly, and maintained as needed, to prevent establishment and growth of infections. At the same time, resources need to be balanced with other physiological processes. On the level of transcription, studies have shown that this balancing act is reflected in tight control of the initiation kinetics and shutdown dynamics of specific immune genes. RESULTS: To investigate genome-wide expression dynamics and trade-offs after infection at a high temporal resolution, we performed an RNA-seq time course on D. melanogaster with 20 time points post Imd stimulation. A combination of methods, including spline fitting, cluster analysis, and Granger causality inference, allowed detailed dissection of expression profiles, lead-lag interactions, and functional annotation of genes through guilt-by-association. We identified Imd-responsive genes and co-expressed, less well characterized genes, with an immediate-early response and sustained up-regulation up to 5 days after stimulation. In contrast, stress response and Toll-responsive genes, among which were Bomanins, demonstrated early and transient responses. We further observed a strong trade-off with metabolic genes, which strikingly recovered to pre-infection levels before the immune response was fully resolved. CONCLUSIONS: This high-dimensional dataset enabled the comprehensive study of immune response dynamics through the parallel application of multiple temporal data analysis methods. The well annotated data set should also serve as a useful resource for further investigation of the D. melanogaster innate immune response, and for the development of methods for analysis of a post-stress transcriptional response time-series at whole-genome scale.

Breadth of humoral immune responses to the C-terminus of the circumsporozoite protein is associated with protective efficacy induced by the RTS,S malaria vaccine.

The circumsporozoite protein (CSP) is the main surface antigen of malaria sporozoites, a prime vaccine target, and is known to have polymorphisms in the C-terminal region. Vaccines using a single allele may have lower efficacy against genotypic variants. Recent studies have found evidence suggesting the efficacy of the CSP-based RTS,S malaria vaccine may be limited against P. falciparum CSP alleles that diverge from the 3D7 vaccine allele, particularly in this polymorphic C-terminal region. In order to assess the breadth of the RTS,S-induced antibody responses against CSP C-terminal antigenic variants, we used a novel multiplex assay to measure reactivity of serum samples from a recent RTS,S study against C-terminal peptides from 3D7 and seven additional CSP alleles that broadly represent the genetic diversity found in circulating P. falciparum field isolates. We found that responses to the variants showed, on average, a ~ 30-fold reduction in reactivity relative to the vaccine-matched 3D7 allele. The extent of this reduction, ranging from 21 to 69-fold, correlated with the number of polymorphisms between the variants and 3D7. We calculated antibody breadth of each sample as the median relative reactivity to the seven CSP variants compared to 3D7. Surprisingly, protection from 3D7 challenge in the RTS,S study was associated with higher C-terminal antibody breadth. These findings suggest CSP C-terminal-specific avidity or fine-specificity may play a role in RTS,S-mediated protection and that breadth of C-terminal CSP-specific antibody responses may be a marker of protection.

Population genomics of Plasmodium vivax in Panama to assess the risk of case importation on malaria elimination.

Malaria incidence in Panama has plateaued in recent years in spite of elimination efforts, with almost all cases caused by Plasmodium vivax. Notwithstanding, overall malaria prevalence remains low (fewer than 1 case per 1000 persons). We used selective whole genome amplification to sequence 59 P. vivax samples from Panama. The P. vivax samples were collected from two periods (2007-2009 and 2017-2019) to study the population structure and transmission dynamics of the parasite. Imported cases resulting from increased levels of human migration could threaten malaria elimination prospects, and four of the samples evaluated came from individuals with travel history. We explored patterns of recent common ancestry among the samples and observed that a highly genetically related lineage (termed CL1) was dominant among the samples (47 out of 59 samples with good sequencing coverage), spanning the entire period of the collection (2007-2019) and all regions of the country. We also found a second, smaller clonal lineage (termed CL2) of four parasites collected between 2017 and 2019. To explore the regional context of Panamanian P. vivax we conducted principal components analysis and constructed a neighbor-joining tree using these samples and samples collected worldwide from a previous study. Three of the four samples with travel history clustered with samples collected from their suspected country of origin (consistent with importation), while one appears to have been a result of local transmission. The small number of Panamanian P. vivax samples not belonging to either CL1 or CL2 clustered with samples collected from Colombia, suggesting they represent the genetically similar ancestral P. vivax population in Panama or were recently imported from Colombia. The low diversity we observe in Panama indicates that this parasite population has been previously subject to a severe bottleneck and may be eligible for elimination. Additionally, while we confirmed that P. vivax is imported to Panama from diverse geographic locations, the lack of impact from imported cases on the overall parasite population genomic profile suggests that onward transmission from such cases is limited and that imported cases may not presently pose a major barrier to elimination.

Local emergence in Amazonia of Plasmodium falciparum k13 C580Y mutants associated with in vitro artemisinin resistance.

Antimalarial drug resistance has historically arisen through convergent de novo mutations in Plasmodium falciparum parasite populations in Southeast Asia and South America. For the past decade in Southeast Asia, artemisinins, the core component of first-line antimalarial therapies, have experienced delayed parasite clearance associated with several pfk13 mutations, primarily C580Y. We report that mutant pfk13 has emerged independently in Guyana, with genome analysis indicating an evolutionary origin distinct from Southeast Asia. Pfk13 C580Y parasites were observed in 1.6% (14/854) of samples collected in Guyana in 2016-2017. Introducing pfk13 C580Y or R539T mutations by gene editing into local parasites conferred high levels of in vitro artemisinin resistance. In vitro growth competition assays revealed a fitness cost associated with these pfk13 variants, potentially explaining why these resistance alleles have not increased in frequency more quickly in South America. These data place local malaria control efforts at risk in the Guiana Shield.

All recommended treatments against malaria include a drug called artemisinin or some of its derivatives. However, there are concerns that Plasmodium falciparum, the parasite that causes most cases of malaria, will eventually develop widespread resistance to the drug. A strain of P. falciparum partially resistant to artemisinin was seen in Cambodia in 2008, and it has since spread across Southeast Asia. The resistance appears to be frequently linked to a mutation known as pfk13 C580Y. Southeast Asia and Amazonia are considered to be hotspots for antimalarial drug resistance, and the pfk13 C580Y mutation was detected in the South American country of Guyana in 2010. To examine whether the mutation was still circulating in this part of the world, Mathieu et al. collected and analyzed 854 samples across Guyana between 2016 and 2017. Overall, 1.6% of the samples had the pfk13 C580Y mutation, but this number was as high as 8.8% in one region. Further analyses revealed that the mutation in Guyana had not spread from Southeast Asia, but that it had occurred in Amazonia independently. To better understand the impact of the pfk13 C580Y mutation, Mathieu et al. introduced this genetic change into non-resistant parasites from a country neighbouring Guyana. As expected, the mutation made P. falciparum highly resistant to artemisinin, but it also slowed the growth rate of the parasite. This disadvantage may explain why the mutation has not spread more rapidly through Guyana in recent years. Artemisinin and its derivatives are always associated with other antimalarial drugs to slow the development of resistance; there are concerns that reduced susceptibility to artemisinin leads to the parasites becoming resistant to the partner drugs. Further research is needed to evaluate how the pfk13 C580Y mutation affects the parasite's response to the typical combination of drugs that are given to patients.

Detection of low-density Plasmodium falciparum infections using amplicon deep sequencing.

BACKGROUND: Deep sequencing of targeted genomic regions is becoming a common tool for understanding the dynamics and complexity of Plasmodium infections, but its lower limit of detection is currently unknown. Here, a new amplicon analysis tool, the Parallel Amplicon Sequencing Error Correction (PASEC) pipeline, is used to evaluate the performance of amplicon sequencing on low-density Plasmodium DNA samples. Illumina-based sequencing of two Plasmodium falciparum genomic regions (CSP and SERA2) was performed on two types of samples: in vitro DNA mixtures mimicking low-density infections (1-200 genomes/µl) and extracted blood spots from a combination of symptomatic and asymptomatic individuals (44-653,080 parasites/µl). Three additional analysis tools-DADA2, HaplotypR, and SeekDeep-were applied to both datasets and the precision and sensitivity of each tool were evaluated. RESULTS: Amplicon sequencing can contend with low-density samples, showing reasonable detection accuracy down to a concentration of 5 Plasmodium genomes/µl. Due to increased stochasticity and background noise, however, all four tools showed reduced sensitivity and precision on samples with very low parasitaemia (< 5 copies/µl) or low read count (< 100 reads per amplicon). PASEC could distinguish major from minor haplotypes with an accuracy of 90% in samples with at least 30 Plasmodium genomes/µl, but only 61% at low Plasmodium concentrations (< 5 genomes/µl) and 46% at very low read counts (< 25 reads per amplicon). The four tools were additionally used on a panel of extracted parasite-positive blood spots from natural malaria infections. While all four identified concordant patterns of complexity of infection (COI) across four sub-Saharan African countries, COI values obtained for individual samples differed in some cases. CONCLUSIONS: Amplicon deep sequencing can be used to determine the complexity and diversity of low-density Plasmodium infections. Despite differences in their approach, four state-of-the-art tools resolved known haplotype mixtures with similar sensitivity and precision. Researchers can therefore choose from multiple robust approaches for analysing amplicon data, however, error filtration approaches should not be uniformly applied across samples of varying parasitaemia. Samples with very low parasitaemia and very low read count have higher false positive rates and call for read count thresholds that are higher than current default recommendations.

Naturally acquired immunity against immature Plasmodium falciparum gametocytes.

The recent decline in global malaria burden has stimulated efforts toward Plasmodium falciparum elimination. Understanding the biology of malaria transmission stages may provide opportunities to reduce or prevent onward transmission to mosquitoes. Immature P. falciparum transmission stages, termed stages I to IV gametocytes, sequester in human bone marrow before release into the circulation as mature stage V gametocytes. This process likely involves interactions between host receptors and potentially immunogenic adhesins on the infected red blood cell (iRBC) surface. Here, we developed a flow cytometry assay to examine immune recognition of live gametocytes of different developmental stages by naturally exposed Malawians. We identified strong antibody recognition of the earliest immature gametocyte-iRBCs (giRBCs) but not mature stage V giRBCs. Candidate surface antigens (n = 30), most of them shared between asexual- and gametocyte-iRBCs, were identified by mass spectrometry and mouse immunizations, as well as correlations between responses by protein microarray and flow cytometry. Naturally acquired responses to a subset of candidate antigens were associated with reduced asexual and gametocyte density, and plasma samples from malaria-infected individuals were able to induce immune clearance of giRBCs in vitro. Infected RBC surface expression of select candidate antigens was validated using specific antibodies, and genetic analysis revealed a subset with minimal variation across strains. Our data demonstrate that humoral immune responses to immature giRBCs and shared iRBC antigens are naturally acquired after malaria exposure. These humoral immune responses may have consequences for malaria transmission potential by clearing developing gametocytes, which could be leveraged for malaria intervention.

Mutations in Plasmodium falciparum actin-binding protein coronin confer reduced artemisinin susceptibility.

Drug resistance is an obstacle to global malaria control, as evidenced by the recent emergence and rapid spread of delayed artemisinin (ART) clearance by mutant forms of the PfKelch13 protein in Southeast Asia. Identifying genetic determinants of ART resistance in African-derived parasites is important for surveillance and for understanding the mechanism of resistance. In this study, we carried out long-term in vitro selection of two recently isolated West African parasites (from Pikine and Thiès, Senegal) with increasing concentrations of dihydroartemisinin (DHA), the biologically active form of ART, over a 4-y period. We isolated two parasite clones, one from each original isolate, that exhibited enhanced survival to DHA in the ring-stage survival assay. Whole-genome sequence analysis identified 10 mutations in seven different genes. We chose to focus on the gene encoding PfCoronin, a member of the WD40-propeller domain protein family, because mutations in this gene occurred in both independent selections, and the protein shares the ß-propeller motif with PfKelch13 protein. For functional validation, when pfcoronin mutations were introduced into the parental parasites by CRISPR/Cas9-mediated gene editing, these mutations were sufficient to reduce ART susceptibility in the parental lines. The discovery of a second gene for ART resistance may yield insights into the molecular mechanisms of resistance. It also suggests that pfcoronin mutants could emerge as a nonkelch13 type of resistance to ART in natural settings.