Subcutaneous Administration of a Monoclonal Antibody to Prevent Malaria

BACKGROUND

Subcutaneous administration of the monoclonal antibody L9LS protected adults against controlled Plasmodium falciparum infection in a phase 1 trial. Whether a monoclonal antibody administered subcutaneously can protect children from P. falciparum infection in a region where this organism is endemic is unclear.

METHODS

We conducted a phase 2 trial in Mali to assess the safety and efficacy of subcutaneous administration of L9LS in children 6 to 10 years of age over a 6-month malaria season. In part A of the trial, safety was assessed at three dose levels in adults, followed by assessment at two dose levels in children. In part B of the trial, children were randomly assigned, in a 1:1:1 ratio, to receive 150 mg of L9LS, 300 mg of L9LS, or placebo. The primary efficacy end point, assessed in a time-to-event analysis, was the first P. falciparum infection, as detected on blood smear performed at least every 2 weeks for 24 weeks. A secondary efficacy end point was the first episode of clinical malaria, as assessed in a time-to-event analysis.

RESULTS

No safety concerns were identified in the dose-escalation part of the trial (part A). In part B, 225 children underwent randomization, with 75 children assigned to each group. No safety concerns were identified in part B. P. falciparum infection occurred in 36 participants (48%) in the 150-mg group, in 30 (40%) in the 300-mg group, and in 61 (81%) in the placebo group. The efficacy of L9LS against P. falciparum infection, as compared with placebo, was 66% (adjusted confidence interval [95% CI], 45 to 79) with the 150-mg dose and 70% (adjusted 95% CI, 50 to 82) with the 300-mg dose (P<0.001 for both comparisons). Efficacy against clinical malaria was 67% (adjusted 95% CI, 39 to 82) with the 150-mg dose and 77% (adjusted 95% CI, 55 to 89) with the 300-mg dose (P<0.001 for both comparisons).

CONCLUSIONS

Subcutaneous administration of L9LS to children was protective against P. falciparum infection and clinical malaria over a period of 6 months. (Funded by the National Institute of Allergy and Infectious Diseases; ClinicalTrials.gov number, NCT05304611.).

A chromosomal reference genome sequence for the malaria mosquito, Anopheles gambiae, Giles, 1902, Ifakara strain

We present a genome assembly from an individual female Anopheles gambiae (the malaria mosquito; Arthropoda; Insecta; Diptera; Culicidae), Ifakara strain. The genome sequence is 264 megabases in span. Most of the assembly is scaffolded into three chromosomal pseudomolecules with the X sex chromosome assembled. The complete mitochondrial genome was also assembled and is 15.4 kilobases in length.

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.

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.

Limited association between Wolbachia and Plasmodium falciparum infections in natural populations of the major malaria mosquito Anopheles moucheti

Since the discovery of natural malaria vector populations infected by the endosymbiont bacterium Wolbachia, a renewed interest has arisen for using this bacterium as an alternative for malaria control. Among naturally infected mosquitoes, Anopheles moucheti, a major malaria mosquito in Central Africa, exhibits one of the highest prevalences of Wolbachia infection. To better understand whether this maternally inherited bacterium could be used for malaria control, we investigated Wolbachia influence in An. moucheti populations naturally infected by the malaria parasite Plasmodium falciparum. To this end, we collected mosquitoes in a village from Cameroon, Central Africa, where this mosquito is the main malaria vector. We found that the prevalence of Wolbachia bacterium was almost fixed in the studied mosquito population, and was higher than previously recorded. We also quantified Wolbachia in whole mosquitoes and dissected abdomens, confirming that the bacterium is also elsewhere than in the abdomen, but at lower density. Finally, we analyzed the association of Wolbachia presence and density on P. falciparum infection. Wolbachia density was slightly higher in mosquitoes infected with the malaria parasite than in uninfected mosquitoes. However, we observed no correlation between the P. falciparum and Wolbachia densities. In conclusion, our study indicates that naturally occurring Wolbachia infection is not associated to P. falciparum development within An. moucheti mosquitoes.

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 ).

A chromosomal reference genome sequence for the malaria mosquito, Anopheles moucheti, Evans, 1925

We present a genome assembly from an individual male Anopheles moucheti (the malaria mosquito; Arthropoda; Insecta; Diptera; Culicidae), from a wild population in Cameroon. The genome sequence is 271 megabases in span. The majority of the assembly is scaffolded into three chromosomal pseudomolecules with the X sex chromosome assembled. The complete mitochondrial genome was also assembled and is 15.5 kilobases in length.

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.

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 genomic epidemiology of infectious disease. AmpSeq data from infectious microbes can inform disease control in multiple ways, including 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 parasitemia 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 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.

Transcriptional correlates of malaria in RTS,S/AS01-vaccinated African children: a matched case–control study

Background:

In a phase 3 trial in African infants and children, the RTS,S/AS01 vaccine (GSK) showed moderate efficacy against clinical malaria. We sought to further understand RTS,S/AS01-induced immune responses associated with vaccine protection.

Methods:

Applying the blood transcriptional module (BTM) framework, we characterized the transcriptomic response to RTS,S/AS01 vaccination in antigen-stimulated (and vehicle control) peripheral blood mononuclear cells sampled from a subset of trial participants at baseline and month 3 (1-month post-third dose). Using a matched case–control study design, we evaluated which of these ‘RTS,S/AS01 signature BTMs’ associated with malaria case status in RTS,S/AS01 vaccinees. Antigen-specific T-cell responses were analyzed by flow cytometry. We also performed a cross-study correlates analysis where we assessed the generalizability of our findings across three controlled human malaria infection studies of healthy, malaria-naive adult RTS,S/AS01 recipients.

Results:

RTS,S/AS01 vaccination was associated with downregulation of B-cell and monocyte-related BTMs and upregulation of T-cell-related BTMs, as well as higher month 3 (vs. baseline) circumsporozoite protein-specific CD4⁺ T-cell responses. There were few RTS,S/AS01-associated BTMs whose month 3 levels correlated with malaria risk. In contrast, baseline levels of BTMs associated with dendritic cells and with monocytes (among others) correlated with malaria risk. The baseline dendritic cell- and monocyte-related BTM correlations with malaria risk appeared to generalize to healthy, malaria-naive adults.

Conclusions:

A prevaccination transcriptomic signature associates with malaria in RTS,S/AS01-vaccinated African children, and elements of this signature may be broadly generalizable. The consistent presence of monocyte-related modules suggests that certain monocyte subsets may inhibit protective RTS,S/AS01-induced responses.

Funding:

Funding was obtained from the NIH-NIAID (R01AI095789), NIH-NIAID (U19AI128914), PATH Malaria Vaccine Initiative (MVI), and Ministerio de Economía y Competitividad (Instituto de Salud Carlos III, PI11/00423 and PI14/01422). The RNA-seq project has been funded in whole or in part with Federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under grant number U19AI110818 to the Broad Institute. This study was also supported by the Vaccine Statistical Support (Bill and Melinda Gates Foundation award INV-008576/OPP1154739 to R.G.). C.D. was the recipient of a Ramon y Cajal Contract from the Ministerio de Economía y Competitividad (RYC-2008-02631). G.M. was the recipient of a Sara Borrell–ISCIII fellowship (CD010/00156) and work was performed with the support of Department of Health, Catalan Government grant (SLT006/17/00109). This research is part of the ISGlobal’s Program on the Molecular Mechanisms of Malaria which is partially supported by the Fundación Ramón Areces and we acknowledge support from the Spanish Ministry of Science and Innovation through the ‘Centro de Excelencia Severo Ochoa 2019–2023’ Program (CEX2018-000806-S), and support from the Generalitat de Catalunya through the CERCA Program.

Integration of whole genome sequencing and transcriptomics reveals a complex picture of the reestablishment of insecticide resistance in the major malaria vector Anopheles coluzzii

Insecticide resistance is a major threat to gains in malaria control, which have been stalling and potentially reversing since 2015. Studies into the causal mechanisms of insecticide resistance are painting an increasingly complicated picture, underlining the need to design and implement targeted studies on this phenotype. In this study, we compare three populations of the major malaria vector An. coluzzii: a susceptible and two resistant colonies with the same genetic background. The original colonised resistant population rapidly lost resistance over a 6-month period, a subset of this population was reselected with pyrethroids, and a third population of this colony that did not lose resistance was also available. The original resistant, susceptible and re-selected colonies were subject to RNAseq and whole genome sequencing, which identified a number of changes across the transcriptome and genome linked with resistance. Firstly, an increase in the expression of genes within the oxidative phosphorylation pathway were seen in both resistant populations compared to the susceptible control; this translated phenotypically through an increased respiratory rate, indicating that elevated metabolism is linked directly with resistance. Genome sequencing highlighted several blocks clearly associated with resistance, including the 2Rb inversion. Finally, changes in the microbiome profile were seen, indicating that the microbial composition may play a role in the resistance phenotype. Taken together, this study reveals a highly complicated phenotype in which multiple transcriptomic, genomic and microbiome changes combine to result in insecticide resistance.

Insecticide resistance in major malaria vectors represents the single biggest threat to malaria control programs, which are heavily reliant upon insecticide-based interventions. Studying resistance using multi-omics approaches has proven difficult due to the use of susceptible comparator populations that have been colonised in a laboratory setting for decades, leading to substantial noise in the data due to differing genetic backgrounds. Here, we utilise a resistant Anopheles coluzzii population from Burkina Faso, a derived population that rapidly lost resistance over a 6-month period, and a population re-selected after loss of resistance to explore causative mechanisms of insecticide resistance. To determine the underlying cause of this phenotype, we use RNAseq, whole genome sequencing and lab-based validation to show changes in respiratory rate, wide-ranging genomic changes and alterations in the microbiome are linked to resistance in this population. These findings demonstrate the complexity of resistance and the challenges in utilising diagnostic markers for resistance in a field setting.

Molecular surveillance of malaria scales up

Parasite and vector genetic data can guide malaria control, and technological advances are enabling more informative genetic data generation at unprecedented scales. Jacob et al. employ multiplexed amplicon sequencing to profile parasite genetic diversity from thousands of malaria samples, illuminating spatiotemporal patterns of drug resistance to inform regional drug policy change.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

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.

Plasmodium simium: population genomics reveals the origin of a reverse zoonosis

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. However, here 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.

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.

A population genomic unveiling of a new cryptic mosquito taxon within the malaria-transmitting Anopheles gambiae complex

The Anopheles gambiae complex consists of multiple morphologically indistinguishable mosquito species including the most important vectors of the malaria parasite Plasmodium falciparum in sub-Saharan Africa. Nine cryptic species have been described so far within the complex. The ecological, immunological and reproductive differences among these species will critically impact population responses to disease control strategies and environmental changes. Here, we examine whole-genome sequencing data from a longitudinal study of putative A. coluzzii in western Burkina Faso. Surprisingly, many specimens are genetically divergent from A. coluzzii and all other Anopheles species and represent a new taxon, here designated Anopheles TENGRELA (AT). Population genetic analysis suggests that the cryptic GOUNDRY subgroup, previously collected as larvae in central Burkina Faso, represents an admixed population descended from both A. coluzzii and AT. AT harbours low nucleotide diversity except for the 2La inversion polymorphism which is maintained by overdominance. It shows numerous fixed differences with A. coluzzii concentrated in several regions reflecting selective sweeps, but the two taxa are identical at standard diagnostic loci used for taxon identification, and thus, AT may often go unnoticed. We present an amplicon-based genotyping assay for identifying AT which could be usefully applied to numerous existing samples. Misidentified cryptic taxa could seriously confound ongoing studies of Anopheles ecology and evolution in western Africa, including phenotypic and genotypic surveys of insecticide resistance. Reproductive barriers between cryptic species may also complicate novel vector control efforts, for example gene drives, and hinder predictions about evolutionary dynamics of Anopheles and Plasmodium.

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.. [DOI: 10.1101/2020.11.15.20232033]

The circumsporozoite protein (CSP) is the main surface antigen of malaria sporozoites and a prime vaccine target. Responses induced by the CSP-based RTS,S vaccine towards the polymorphic C-terminal region of P.falciparum-CSP raise concerns that vaccines using single alleles may have lower efficacy against genotypic variants. We characterized the extent of C-terminal cross-reactivity of antibodies induced by RTS,S (based on the 3D7 allele) with variants representing seven circulating field isolates through a novel HTS-multiplex assay for screening closely related peptides. Reactivity to variants showed approximately 30-fold reduction in recognition relative to 3D7. The degree of reduced cross-reactivity,ranging from 21 to 69-fold, directly correlated with the number of polymorphisms between variants and 3D7. Surprisingly, protection assessed by challenge with 3D7 parasites was strongly associated with higher C-terminal antibody breadth suggesting that C-terminal specific avidity or fine-specificity may play a role in RTS,S/AS01B-mediated protection and that breadth of C-terminal CSP-specific antibody responses may be a marker of protection.

Identity-by-descent with uncertainty characterises connectivity of Plasmodium falciparum populations on the Colombian-Pacific coast

Characterising connectivity between geographically separated biological populations is a common goal in many fields. Recent approaches to understanding connectivity between malaria parasite populations, with implications for disease control efforts, have used estimates of relatedness based on identity-by-descent (IBD). However, uncertainty around estimated relatedness has not been accounted for. IBD-based relatedness estimates with uncertainty were computed for pairs of monoclonal Plasmodium falciparum samples collected from five cities on the Colombian-Pacific coast where long-term clonal propagation of P. falciparum is frequent. The cities include two official ports, Buenaventura and Tumaco, that are separated geographically but connected by frequent marine traffic. Fractions of highly-related sample pairs (whose classification using a threshold accounts for uncertainty) were greater within cities versus between. However, based on both highly-related fractions and on a threshold-free approach (Wasserstein distances between parasite populations) connectivity between Buenaventura and Tumaco was disproportionally high. Buenaventura-Tumaco connectivity was consistent with transmission events involving parasites from five clonal components (groups of statistically indistinguishable parasites identified under a graph theoretic framework). To conclude, P. falciparum population connectivity on the Colombian-Pacific coast abides by accessibility not isolation-by-distance, potentially implicating marine traffic in malaria transmission with opportunities for targeted intervention. Further investigations are required to test this hypothesis. For the first time in malaria epidemiology (and to our knowledge in ecological and epidemiological studies more generally), we account for uncertainty around estimated relatedness (an important consideration for studies that plan to use genotype versus whole genome sequence data to estimate IBD-based relatedness); we also use threshold-free methods to compare parasite populations and identify clonal components. Threshold-free methods are especially important in analyses of malaria parasites and other recombining organisms with mixed mating systems where thresholds do not have clear interpretation (e.g. due to clonal propagation) and thus undermine the cross-comparison of studies.

Population genomics reveals the expansion of highly inbred Plasmodium vivax lineages in the main malaria hotspot of Brazil

BACKGROUND

Plasmodium vivax is a neglected human malaria parasite that causes significant morbidity in the Americas, the Middle East, Asia, and the Western Pacific. Population genomic approaches remain little explored to map local and regional transmission pathways of P. vivax across the main endemic sites in the Americas, where great progress has been made towards malaria elimination over the past decades.

METHODOLOGY/PRINCIPAL FINDINGS

We analyze 38 patient-derived P. vivax genome sequences from Mâncio Lima (ML)-the Amazonian malaria hotspot next to the Brazil-Peru border-and 24 sequences from two other sites in Acre State, Brazil, a country that contributes 23% of malaria cases in the Americas. We show that the P. vivax population of ML is genetically diverse (π = 4.7 × 10-4), with a high polymorphism particularly in genes encoding proteins putatively involved in red blood cell invasion. Paradoxically, however, parasites display strong genome-wide linkage disequilibrium, being fragmented into discrete lineages that are remarkably stable across time and space, with only occasional recombination between them. Using identity-by-descent approaches, we identified a large cluster of closely related sequences that comprises 16 of 38 genomes sampled in ML over 26 months. Importantly, we found significant ancestry sharing between parasites at a large geographic distance, consistent with substantial gene flow between regional P. vivax populations.

CONCLUSIONS/SIGNIFICANCE

We have characterized the sustained expansion of highly inbred P. vivax lineages in a malaria hotspot that can seed regional transmission. Potential source populations in hotspots represent a priority target for malaria elimination in the Amazon.

The Anopheles coluzzii microbiome and its interaction with the intracellular parasite Wolbachia

Wolbachia, an endosymbiotic alpha-proteobacterium commonly found in insects, can inhibit the transmission of human pathogens by mosquitoes. Biocontrol programs are underway using Aedes aegypti mosquitoes trans-infected with a non-natural Wolbachia strain to reduce dengue virus transmission. Less is known about the impact of Wolbachia on the biology and vectorial capacity of Anopheles mosquitoes, the vectors of malaria parasites. A naturally occurring strain of Wolbachia, wAnga, infects populations of the major malaria vectors Anopheles gambiae and Anopheles coluzzii in Burkina Faso. Previous studies found wAnga infection was negatively correlated with Plasmodium infection in the mosquito and wAnga influenced mosquito egg-laying behavior. Here, we investigate wAnga in natural populations of An. coluzzii and its interactions with other resident microbiota using targeted 16S sequencing. Though we find no major differences in microbiota composition associated with wAnga infection, we do find several taxa that correlate with the presence or absence of wAnga in female mosquitoes following oviposition, with the caveat that we could not rule out batch effects due to the unanticipated impact of wAnga on oviposition timing. These data suggest wAnga may influence or interact with the Anopheles microbiota, which may contribute to the impact of wAnga on Anopheles biology and vectorial capacity.

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.

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.

Amplicon deep sequencing of kelch13 in Plasmodium falciparum isolates from Senegal

BACKGROUND

In 2006, the Senegalese National Malaria Control Programme recommended artemisinin-based combination therapy (ACT) with artemether-lumefantrine as the first-line treatment for uncomplicated Plasmodium falciparum malaria. To date, multiple mutations associated with artemisinin delayed parasite clearance have been described in Southeast Asia in the Pfk13 gene, such as Y493H, R539T, I543T and C580Y. Even though ACT remains clinically and parasitologically efficacious in Senegal, the spread of resistance is possible as shown by the earlier emergence of resistance to chloroquine in Southeast Asia that subsequently spread to Africa. Therefore, surveillance of artemisinin resistance in malaria endemic regions is crucial and requires the implementation of sensitive tools, such as next-generation sequencing (NGS) which can detect novel mutations at low frequency.

METHODS

Here, an amplicon sequencing approach was used to identify mutations in the Pfk13 gene in eighty-one P. falciparum isolates collected from three different regions of Senegal.

RESULTS

In total, 10 SNPs around the propeller domain were identified; one synonymous SNP and nine non-synonymous SNPs, and two insertions. Three of these SNPs (T478T, A578S and V637I) were located in the propeller domain. A578S, is the most frequent mutation observed in Africa, but has not previously been reported in Senegal. A previous study has suggested that A578S could disrupt the function of the Pfk13 propeller region.

CONCLUSION

As the genetic basis of possible artemisinin resistance may be distinct in Africa and Southeast Asia, further studies are necessary to assess the new SNPs reported in this study.

Linked-read sequencing identifies abundant microinversions and introgression in the arboviral vector Aedes aegypti

Background

Aedes aegypti is the principal mosquito vector of Zika, dengue, and yellow fever viruses. Two subspecies of Ae. aegypti exhibit phenotypic divergence with regard to habitat, host preference, and vectorial capacity. Chromosomal inversions have been shown to play a major role in adaptation and speciation in dipteran insects and would be of great utility for studies of Ae. aegypti. However, the large and highly repetitive genome of Ae. aegypti makes it difficult to detect inversions with paired-end short-read sequencing data, and polytene chromosome analysis does not provide sufficient resolution to detect chromosome banding patterns indicative of inversions.

Results

To characterize chromosomal diversity in this species, we have carried out deep Illumina sequencing of linked-read (10X Genomics) libraries in order to discover inversion loci as well as SNPs. We analyzed individuals from colonies representing the geographic limits of each subspecies, one contact zone between subspecies, and a closely related sister species. Despite genome-wide SNP divergence and abundant microinversions, we do not find any inversions occurring as fixed differences between subspecies. Many microinversions are found in regions that have introgressed and have captured genes that could impact behavior, such as a cluster of odorant-binding proteins that may play a role in host feeding preference.

Conclusions

Our study shows that inversions are abundant and widely shared among subspecies of Aedes aegypti and that introgression has occurred in regions of secondary contact. This library of 32 novel chromosomal inversions demonstrates the capacity for linked-read sequencing to identify previously intractable genomic rearrangements and provides a foundation for future population genetics studies in this species.

Resolving the cause of recurrent Plasmodium vivax malaria probabilistically

Relapses arising from dormant liver-stage Plasmodium vivax parasites (hypnozoites) are a major cause of vivax malaria. However, in endemic areas, a recurrent blood-stage infection following treatment can be hypnozoite-derived (relapse), a blood-stage treatment failure (recrudescence), or a newly acquired infection (reinfection). Each of these requires a different prevention strategy, but it was not previously possible to distinguish between them reliably. We show that individual vivax malaria recurrences can be characterised probabilistically by combined modelling of time-to-event and genetic data within a framework incorporating identity-by-descent. Analysis of pooled patient data on 1441 recurrent P. vivax infections in 1299 patients on the Thailand-Myanmar border observed over 1000 patient follow-up years shows that, without primaquine radical curative treatment, 3 in 4 patients relapse. In contrast, after supervised high-dose primaquine only 1 in 40 relapse. In this region of frequent relapsing P. vivax, failure rates after supervised high-dose primaquine are significantly lower (∼3%) than estimated previously.

Genomic architecture and introgression shape a butterfly radiation

We used 20 de novo genome assemblies to probe the speciation history and architecture of gene flow in rapidly radiating Heliconius butterflies. Our tests to distinguish incomplete lineage sorting from introgression indicate that gene flow has obscured several ancient phylogenetic relationships in this group over large swathes of the genome. Introgressed loci are underrepresented in low-recombination and gene-rich regions, consistent with the purging of foreign alleles more tightly linked to incompatibility loci. Here, we identify a hitherto unknown inversion that traps a color pattern switch locus. We infer that this inversion was transferred between lineages by introgression and is convergent with a similar rearrangement in another part of the genus. These multiple de novo genome sequences enable improved understanding of the importance of introgression and selective processes in adaptive radiation.

Estimating Relatedness Between Malaria Parasites

Understanding the relatedness of individuals within or between populations is a common goal in biology. Increasingly, relatedness features in genetic epidemiology studies of pathogens. These studies are relatively new compared to those in humans and other organisms, but are important for designing interventions and understanding pathogen transmission. Only recently have researchers begun to routinely apply relatedness to apicomplexan eukaryotic malaria parasites, and to date have used a range of different approaches on an ad hoc basis. Therefore, it remains unclear how to compare different studies and which measures to use. Here, we systematically compare measures based on identity-by-state (IBS) and identity-by-descent (IBD) using a globally diverse data set of malaria parasites, Plasmodium falciparum and P. vivax, and provide marker requirements for estimates based on IBD. We formally show that the informativeness of polyallelic markers for relatedness inference is maximized when alleles are equifrequent. Estimates based on IBS are sensitive to allele frequencies, which vary across populations and by experimental design. For portability across studies, we thus recommend estimates based on IBD. To generate estimates with errors below an arbitrary threshold of 0.1, we recommend ∼100 polyallelic or 200 biallelic markers. Marker requirements are immediately applicable to haploid malaria parasites and other haploid eukaryotes. C.I.s facilitate comparison when different marker sets are used. This is the first attempt to provide rigorous analysis of the reliability of, and requirements for, relatedness inference in malaria genetic epidemiology. We hope it will provide a basis for statistically informed prospective study design and surveillance strategies.

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.

De Novo Mutations Resolve Disease Transmission Pathways in Clonal Malaria

Detecting de novo mutations in viral and bacterial pathogens enables researchers to reconstruct detailed networks of disease transmission and is a key technique in genomic epidemiology. However, these techniques have not yet been applied to the malaria parasite, Plasmodium falciparum, in which a larger genome, slower generation times, and a complex life cycle make them difficult to implement. Here, we demonstrate the viability of de novo mutation studies in P. falciparum for the first time. Using a combination of sequencing, library preparation, and genotyping methods that have been optimized for accuracy in low-complexity genomic regions, we have detected de novo mutations that distinguish nominally identical parasites from clonal lineages. Despite its slower evolutionary rate compared with bacterial or viral species, de novo mutation can be detected in P. falciparum across timescales of just 1–2 years and evolutionary rates in low-complexity regions of the genome can be up to twice that detected in the rest of the genome. The increased mutation rate allows the identification of separate clade expansions that cannot be found using previous genomic epidemiology approaches and could be a crucial tool for mapping residual transmission patterns in disease elimination campaigns and reintroduction scenarios.

Correction to: Mapping malaria by combining parasite genomic and epidemiologic data

The original article [1] contained an error in the presentation of Figure 1; this error has now been rectified and Figure 1 is now presented correctly.

Improved reference genome of Aedes aegypti informs arbovirus vector control

Female Aedes aegypti mosquitoes infect more than 400 million people each year with dangerous viral pathogens including dengue, yellow fever, Zika and chikungunya. Progress in understanding the biology of mosquitoes and developing the tools to fight them has been slowed by the lack of a high-quality genome assembly. Here we combine diverse technologies to produce the markedly improved, fully re-annotated AaegL5 genome assembly, and demonstrate how it accelerates mosquito science. We anchored physical and cytogenetic maps, doubled the number of known chemosensory ionotropic receptors that guide mosquitoes to human hosts and egg-laying sites, provided further insight into the size and composition of the sex-determining M locus, and revealed copy-number variation among glutathione S-transferase genes that are important for insecticide resistance. Using high-resolution quantitative trait locus and population genomic analyses, we mapped new candidates for dengue vector competence and insecticide resistance. AaegL5 will catalyse new biological insights and intervention strategies to fight this deadly disease vector.

An improved, fully re-annotated Aedes aegypti genome assembly (AaegL5) provides insights into the sex-determining M locus, chemosensory systems that help mosquitoes to hunt humans and loci involved in insecticide resistance and will help to generate intervention strategies to fight this deadly disease vector.

Mapping malaria by combining parasite genomic and epidemiologic data

BACKGROUND

Recent global progress in scaling up malaria control interventions has revived the goal of complete elimination in many countries. Decreasing transmission intensity generally leads to increasingly patchy spatial patterns of malaria transmission in elimination settings, with control programs having to accurately identify remaining foci in order to efficiently target interventions.

FINDINGS

The role of connectivity between different pockets of local transmission is of increasing importance as programs near elimination since humans are able to transfer parasites beyond the limits of mosquito dispersal, thus re-introducing parasites to previously malaria-free regions. Here, we discuss recent advances in the quantification of spatial epidemiology of malaria, particularly Plasmodium falciparum, in the context of transmission reduction interventions. Further, we highlight the challenges and promising directions for the development of integrated mapping, modeling, and genomic approaches that leverage disparate datasets to measure both connectivity and transmission.

CONCLUSION

A more comprehensive understanding of the spatial transmission of malaria can be gained using a combination of parasite genetics and epidemiological modeling and mapping. However, additional molecular and quantitative methods are necessary to answer these public health-related questions.

hmmIBD: software to infer pairwise identity by descent between haploid genotypes

Background

A number of recent malaria studies have used identity by descent (IBD) to study epidemiological processes relevant to malaria control. In this paper, a software package, hmmIBD, is introduced for estimating pairwise IBD between haploid genomes, such as those of the malaria parasite, sampled from one or two populations. Source code is freely available.

Methods

The performance of hmmIBD was verified using simulated data and benchmarked against an existing method for detecting IBD within populations. Code for all tests is freely available. The utility of hmmIBD for detecting IBD across populations was demonstrated using Plasmodium falciparum data from Cambodia and Ghana.

Results

Alongside an existing method, hmmIBD was highly accurate, sensitive and specific. It is fast, requiring only 70 s on average to analyse 50 whole genome sequences on a laptop computer, and scales linearly in the number of pairwise comparisons. Treatment of different populations under hmmIBD improves detection of IBD across populations.

Conclusion

Fast and accurate software for detecting IBD in malaria parasite genetic data sampled from one or two populations is presented. The latter will likely be a useful feature for malaria elimination efforts, since it could facilitate identification of imported malaria cases. Software is robust to possible misspecification of the genotyping error and the recombination rate. However, exclusion of data in regions whose rates vary greatly from their genome-wide average is recommended.

Electronic supplementary material

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

Transmission of molecularly undetectable circulating parasite clones leads to high infection complexity in mosquitoes post feeding

•

Additional parasite alleles were consistently identified in mosquitoes compared with the human blood sample they had fed on.

•

Assessments of Plasmodium falciparum complexity relying on single time-point collections miss transmissible clones.

•

Low-density gametocyte – producing clones are capable of successfully establishing infections in mosquitoes.

Plasmodium falciparum malaria infections often comprise multiple distinct parasite clones. Few datasets have directly assessed infection complexity in humans and mosquitoes they infect. Examining parasites using molecular tools may provide insights into the selective transmissibility of isolates. Using capillary electrophoresis genotyping and next generation amplicon sequencing, we analysed complexity of parasite infections in human blood and in the midguts of mosquitoes that became infected in membrane feeding experiments using the same blood material in two West African settings. Median numbers of clones in humans and mosquitoes were higher in samples from Burkina Faso (4.5, interquartile range 2–8 for humans; and 2, interquartile range 1–3 for mosquitoes) than in The Gambia (2, interquartile range 1–3 and 1, interquartile range 1–3, for humans and mosquitoes, respectively). Whilst the median number of clones was commonly higher in human blood samples, not all transmitted alleles were detectable in the human peripheral blood. In both study sample sets, additional parasite alleles were identified in mosquitoes compared with the matched human samples (10–88.9% of all clones/feeding assay, n = 73 feeding assays). The results are likely due to preferential amplification of the most abundant clones in peripheral blood but confirm the presence of low density clones that produce transmissible sexual stage parasites.

Host-mediated selection impacts the diversity of Plasmodium falciparum antigens within infections

Host immunity exerts strong selective pressure on pathogens. Population-level genetic analysis can identify signatures of this selection, but these signatures reflect the net selective effect of all hosts and vectors in a population. In contrast, analysis of pathogen diversity within hosts provides information on individual, host-specific selection pressures. Here, we combine these complementary approaches in an analysis of the malaria parasite Plasmodium falciparum using haplotype sequences from thousands of natural infections in sub-Saharan Africa. We find that parasite genotypes show preferential clustering within multi-strain infections in young children, and identify individual amino acid positions that may contribute to strain-specific immunity. Our results demonstrate that natural host defenses to P. falciparum act in an allele-specific manner to block specific parasite haplotypes from establishing blood-stage infections. This selection partially explains the extreme amino acid diversity of many parasite antigens and suggests that vaccines targeting such proteins should account for allele-specific immunity.

Host immune responses exert selective pressure on Plasmodium falciparum. Here, the authors show that allele-specific immunity impacts the antigenic diversity of individual malaria infections. This process partially explains the extreme amino acid diversity of many parasite antigens and suggests that vaccines should account for allele-specific immunity.

Probing Plasmodium falciparum sexual commitment at the single-cell level

Background: Malaria parasites go through major transitions during their complex life cycle, yet the underlying differentiation pathways remain obscure. Here we apply single cell transcriptomics to unravel the program inducing sexual differentiation in Plasmodium falciparum. Parasites have to make this essential life-cycle decision in preparation for human-to-mosquito transmission. Methods: By combining transcriptional profiling with quantitative imaging and genetics, we defined a transcriptional signature in sexually committed cells. Results: We found this transcriptional signature to be distinct from general changes in parasite metabolism that can be observed in response to commitment-inducing conditions. Conclusions: This proof-of-concept study provides a template to capture transcriptional diversity in parasite populations containing complex mixtures of different life-cycle stages and developmental programs, with important implications for our understanding of parasite biology and the ongoing malaria elimination campaign.

Probing Plasmodium falciparum sexual commitment at the single-cell level

Background: Malaria parasites go through major transitions during their complex life cycle, yet the underlying differentiation pathways remain obscure. Here we apply single cell transcriptomics to unravel the program inducing sexual differentiation in Plasmodium falciparum. Parasites have to make this essential life-cycle decision in preparation for human-to-mosquito transmission. Methods: By combining transcriptional profiling with quantitative imaging and genetics, we defined a transcriptional signature in sexually committed cells. Results: We found this transcriptional signature to be distinct from general changes in parasite metabolism that can be observed in response to commitment-inducing conditions. Conclusions: This proof-of-concept study provides a template to capture transcriptional diversity in parasite populations containing complex mixtures of different life-cycle stages and developmental programs, with important implications for our understanding of parasite biology and the ongoing malaria elimination campaign.

Probing Plasmodium falciparum sexual commitment at the single-cell level

Background: Malaria parasites go through major transitions during their complex life cycle, yet the underlying differentiation pathways remain obscure. Here we apply single cell transcriptomics to unravel the program inducing sexual differentiation in Plasmodium falciparum. Parasites have to make this essential life-cycle decision in preparation for human-to-mosquito transmission. Methods: By combining transcriptional profiling with quantitative imaging and genetics, we defined a transcriptional signature in sexually committed cells. Results: We found this transcriptional signature to be distinct from general changes in parasite metabolism that can be observed in response to commitment-inducing conditions. Conclusions: This proof-of-concept study provides a template to capture transcriptional diversity in parasite populations containing complex mixtures of different life-cycle stages and developmental programs, with important implications for our understanding of parasite biology and the ongoing malaria elimination campaign.

Probing Plasmodium falciparum sexual commitment at the single-cell level

Background: Malaria parasites go through major transitions during their complex life cycle, yet the underlying differentiation pathways remain obscure. Here we apply single cell transcriptomics to unravel the program inducing sexual differentiation in Plasmodium falciparum. Parasites have to make this essential life-cycle decision in preparation for human-to-mosquito transmission. Methods: By combining transcriptional profiling with quantitative imaging and genetics, we defined a transcriptional signature in sexually committed cells. Results: We found this transcriptional signature to be distinct from general changes in parasite metabolism that can be observed in response to commitment-inducing conditions. Conclusions: This proof-of-concept study provides a template to capture transcriptional diversity in parasite populations containing complex mixtures of different life-cycle stages and developmental programs, with important implications for our understanding of parasite biology and the ongoing malaria elimination campaign.

Quantitative Proteomic Profiling Reveals Novel Plasmodium falciparum Surface Antigens and Possible Vaccine Candidates

Despite recent efforts toward control and elimination, malaria remains a major public health problem worldwide. Plasmodium falciparum resistance against artemisinin, used in front line combination drugs, is on the rise, and the only approved vaccine shows limited efficacy. Combinations of novel and tailored drug and vaccine interventions are required to maintain the momentum of the current malaria elimination program. Current evidence suggests that strain-transcendent protection against malaria infection can be achieved using whole organism vaccination or with a polyvalent vaccine covering multiple antigens or epitopes. These approaches have been successfully applied to the human-infective sporozoite stage. Both systemic and tissue-specific pathology during infection with the human malaria parasite P. falciparum is caused by asexual blood stages. Tissue tropism and vascular sequestration are the result of specific binding interactions between antigens on the parasite-infected red blood cell (pRBC) surface and endothelial receptors. The major surface antigen and parasite ligand binding to endothelial receptors, PfEMP1 is encoded by about 60 variants per genome and shows high sequence diversity across strains. Apart from PfEMP1 and three additional variant surface antigen families RIFIN, STEVOR, and SURFIN, systematic analysis of the infected red blood cell surface is lacking. Here we present the most comprehensive proteomic investigation of the parasitized red blood cell surface so far. Apart from the known variant surface antigens, we identified a set of putative single copy surface antigens with low sequence diversity, several of which are validated in a series of complementary experiments. Further functional and immunological investigation is underway to test these novel P. falciparum blood stage proteins as possible vaccine candidates.

Quantifying connectivity between local Plasmodium falciparum malaria parasite populations using identity by descent

With the rapidly increasing abundance and accessibility of genomic data, there is a growing interest in using population genetic approaches to characterize fine-scale dispersal of organisms, providing insight into biological processes across a broad range of fields including ecology, evolution and epidemiology. For sexually recombining haploid organisms such as the human malaria parasite P. falciparum, however, there have been no systematic assessments of the type of data and methods required to resolve fine scale connectivity. This analytical gap hinders the use of genomics for understanding local transmission patterns, a crucial goal for policy makers charged with eliminating this important human pathogen. Here we use data collected from four clinics with a catchment area spanning approximately 120 km of the Thai-Myanmar border to compare the ability of divergence (FST) and relatedness based on identity by descent (IBD) to resolve spatial connectivity between malaria parasites collected from proximal clinics. We found no relationship between inter-clinic distance and FST, likely due to sampling of highly related parasites within clinics, but a significant decline in IBD-based relatedness with increasing inter-clinic distance. This association was contingent upon the data set type and size. We estimated that approximately 147 single-infection whole genome sequenced parasite samples or 222 single-infection parasite samples genotyped at 93 single nucleotide polymorphisms (SNPs) were sufficient to recover a robust spatial trend estimate at this scale. In summary, surveillance efforts cannot rely on classical measures of genetic divergence to measure P. falciparum transmission on a local scale. Given adequate sampling, IBD-based relatedness provides a useful alternative, and robust trends can be obtained from parasite samples genotyped at approximately 100 SNPs.

Malaria Genomics in the Era of Eradication

The first reference genome assembly for the Plasmodium falciparum malaria parasite was completed over a decade ago, and the impact of this and other genomic resources on malaria research has been significant. Genomic resources for other malaria parasites are being established, even as P. falciparum continues to be the focus of development of new genomic methods and applications. Here we review the impact and applications of genomic data on malaria research, and discuss future needs and directions as genomic data generation becomes less expensive and more decentralized. Specifically, we focus on how population genomic strategies can be utilized to advance the malaria eradication agenda.

Longitudinal genomic surveillance of Plasmodium falciparum malaria parasites reveals complex genomic architecture of emerging artemisinin resistance

Background

Artemisinin-based combination therapies are the first line of treatment for Plasmodium falciparum infections worldwide, but artemisinin resistance has risen rapidly in Southeast Asia over the past decade. Mutations in the kelch13 gene have been implicated in this resistance. We used longitudinal genomic surveillance to detect signals in kelch13 and other loci that contribute to artemisinin or partner drug resistance. We retrospectively sequenced the genomes of 194 P. falciparum isolates from five sites in Northwest Thailand, over the period of a rapid increase in the emergence of artemisinin resistance (2001–2014).

Results

We evaluate statistical metrics for temporal change in the frequency of individual SNPs, assuming that SNPs associated with resistance increase in frequency over this period. After Kelch13-C580Y, the strongest temporal change is seen at a SNP in phosphatidylinositol 4-kinase, which is involved in a pathway recently implicated in artemisinin resistance. Furthermore, other loci exhibit strong temporal signatures which warrant further investigation for involvement in artemisinin resistance evolution. Through genome-wide association analysis we identify a variant in a kelch domain-containing gene on chromosome 10 that may epistatically modulate artemisinin resistance.

Conclusions

This analysis demonstrates the potential of a longitudinal genomic surveillance approach to detect resistance-associated gene loci to improve our mechanistic understanding of how resistance develops. Evidence for additional genomic regions outside of the kelch13 locus associated with artemisinin-resistant parasites may yield new molecular markers for resistance surveillance, which may be useful in efforts to reduce the emergence or spread of artemisinin resistance in African parasite populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1204-4) contains supplementary material, which is available to authorized users.

Genetic relatedness analysis reveals the cotransmission of genetically related Plasmodium falciparum parasites in Thiès, Senegal

Background

As public health interventions drive parasite populations to elimination, genetic epidemiology models that incorporate population genomics can be powerful tools for evaluating the effectiveness of continued intervention. However, current genetic epidemiology models may not accurately simulate the population genetic profile of parasite populations, particularly with regard to polygenomic (multi-strain) infections. Current epidemiology models simulate polygenomic infections via superinfection (multiple mosquito bites), despite growing evidence that cotransmission (a single mosquito bite) may contribute to polygenomic infections.

Methods

Here, we quantified the relatedness of strains within 31 polygenomic infections collected from patients in Thiès, Senegal using a hidden Markov model to measure the proportion of the genome that is inferred to be identical by descent.

Results

We found that polygenomic infections can be composed of highly related parasites and that superinfection models drastically underestimate the relatedness of strains within polygenomic infections.

Conclusions

Our findings suggest that cotransmission is a major contributor to polygenomic infections in Thiès, Senegal. The incorporation of cotransmission into existing genetic epidemiology models may enhance our ability to characterize and predict changes in population structure associated with reduced transmission intensities and the emergence of important phenotypes like drug resistance that threaten to undermine malaria elimination activities.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0398-0) contains supplementary material, which is available to authorized users.

THE REAL McCOIL: A method for the concurrent estimation of the complexity of infection and SNP allele frequency for malaria parasites

As many malaria-endemic countries move towards elimination of Plasmodium falciparum, the most virulent human malaria parasite, effective tools for monitoring malaria epidemiology are urgent priorities. P. falciparum population genetic approaches offer promising tools for understanding transmission and spread of the disease, but a high prevalence of multi-clone or polygenomic infections can render estimation of even the most basic parameters, such as allele frequencies, challenging. A previous method, COIL, was developed to estimate complexity of infection (COI) from single nucleotide polymorphism (SNP) data, but relies on monogenomic infections to estimate allele frequencies or requires external allele frequency data which may not available. Estimates limited to monogenomic infections may not be representative, however, and when the average COI is high, they can be difficult or impossible to obtain. Therefore, we developed THE REAL McCOIL, Turning HEterozygous SNP data into Robust Estimates of ALelle frequency, via Markov chain Monte Carlo, and Complexity Of Infection using Likelihood, to incorporate polygenomic samples and simultaneously estimate allele frequency and COI. This approach was tested via simulations then applied to SNP data from cross-sectional surveys performed in three Ugandan sites with varying malaria transmission. We show that THE REAL McCOIL consistently outperforms COIL on simulated data, particularly when most infections are polygenomic. Using field data we show that, unlike with COIL, we can distinguish epidemiologically relevant differences in COI between and within these sites. Surprisingly, for example, we estimated high average COI in a peri-urban subregion with lower transmission intensity, suggesting that many of these cases were imported from surrounding regions with higher transmission intensity. THE REAL McCOIL therefore provides a robust tool for understanding the molecular epidemiology of malaria across transmission settings.

Mosquito Vectors and the Globalization of Plasmodium falciparum Malaria

Plasmodium falciparum malaria remains a devastating public health problem. Recent discoveries have shed light on the origin and evolution of Plasmodium parasites and their interactions with their vertebrate and mosquito hosts. P. falciparum malaria originated in Africa from a single horizontal transfer between an infected gorilla and a human, and became global as the result of human migration. Today, P. falciparum malaria is transmitted worldwide by more than 70 different anopheline mosquito species. Recent studies indicate that the mosquito immune system can be a barrier to malaria transmission and that the P. falciparum Pfs47 gene allows the parasite to evade mosquito immune detection. Here, we review the origin and globalization of P. falciparum and integrate this history with analysis of the biology, evolution, and dispersal of the main mosquito vectors. This new perspective broadens our understanding of P. falciparum population structure and the dispersal of important parasite genetic traits.

Population genomics studies identify signatures of global dispersal and drug resistance in Plasmodium vivax

Plasmodium vivax is a major public health burden, responsible for the majority of malaria infections outside Africa. We explored the impact of demographic history and selective pressures on the P. vivax genome by sequencing 182 clinical isolates sampled from 11 countries across the globe, using hybrid selection to overcome human DNA contamination. We confirmed previous reports of high genomic diversity in P. vivax relative to the more virulent Plasmodium falciparum species; regional populations of P. vivax exhibited greater diversity than the global P. falciparum population, indicating a large and/or stable population. Signals of natural selection suggest that P. vivax is evolving in response to antimalarial drugs and is adapting to regional differences in the human host and the mosquito vector. These findings underline the variable epidemiology of this parasite species and highlight the breadth of approaches that may be required to eliminate P. vivax globally.

Genomic Dark Matter Illuminated: Anopheles Y Chromosomes

Hall et al. have strategically used long-read sequencing technology to characterize the structure and highly repetitive content of the Y chromosome in Anopheles malaria mosquitoes. Their work confirms that this important but elusive heterochromatic sex chromosome is evolving extremely rapidly and harbors a remarkably small number of genes.