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

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal.

We here analyze data from the first year of an ongoing nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal. The analysis is based on 1097 samples collected at health facilities during passive malaria case detection in 2019; it provides a baseline for analyzing parasite genetic metrics as they vary over time and geographic space. The study's goal was to identify genetic metrics that were informative about transmission intensity and other aspects of transmission dynamics, focusing on measures of genetic relatedness between parasites. We found the best genetic proxy for local malaria incidence to be the proportion of polygenomic infections (those with multiple genetically distinct parasites), although this relationship broke down at low incidence. The proportion of related parasites was less correlated with incidence while local genetic diversity was uninformative. The type of relatedness could discriminate local transmission patterns: two nearby areas had similarly high fractions of relatives, but one was dominated by clones and the other by outcrossed relatives. Throughout Senegal, 58% of related parasites belonged to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci and at one novel locus, reflective of ongoing selection pressure.

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal.

Parasite genetic surveillance has the potential to play an important role in malaria control. We describe here an analysis of data from the first year of an ongoing, nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal, intended to provide actionable information for malaria control efforts. Looking for a good proxy for local malaria incidence, we found that the best predictor was the proportion of polygenomic infections (those with multiple genetically distinct parasites), although that relationship broke down in very low incidence settings (r = 0.77 overall). The proportion of closely related parasites in a site was more weakly correlated ( r = -0.44) with incidence while the local genetic diversity was uninformative. Study of related parasites indicated their potential for discriminating local transmission patterns: two nearby study areas had similarly high fractions of relatives, but one area was dominated by clones and the other by outcrossed relatives. Throughout the country, 58% of related parasites proved to belong to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci as well as at one novel locus, reflective of ongoing selection pressure.

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.

Plasmodium falciparum genomic surveillance reveals spatial and temporal trends, association of genetic and physical distance, and household clustering.

Molecular epidemiology using genomic data can help identify relationships between malaria parasite population structure, malaria transmission intensity, and ultimately help generate actionable data to assess the effectiveness of malaria control strategies. Genomic data, coupled with geographic information systems data, can further identify clusters or hotspots of malaria transmission, parasite genetic and spatial connectivity, and parasite movement by human or mosquito mobility over time and space. In this study, we performed longitudinal genomic surveillance in a cohort of 70 participants over four years from different neighborhoods and households in Thiès, Senegal-a region of exceptionally low malaria transmission (entomological inoculation rate less than 1). Genetic identity (identity by state, IBS) was established using a 24-single nucleotide polymorphism molecular barcode, identity by descent was calculated from whole genome sequence data, and a hierarchical Bayesian regression model was used to establish genetic and spatial relationships. Our results show clustering of genetically similar parasites within households and a decline in genetic similarity of parasites with increasing distance. One household showed extremely high diversity and warrants further investigation as to the source of these diverse genetic types. This study illustrates the utility of genomic data with traditional epidemiological approaches for surveillance and detection of trends and patterns in malaria transmission not only by neighborhood but also by household. This approach can be implemented regionally and countrywide to strengthen and support malaria control and elimination 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.

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.

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.