Role of Pfs47 in the dispersal of ancestral Plasmodium falciparum malaria through adaptation to different anopheline vectors.

Plasmodium falciparum malaria originated when Plasmodium praefalciparum, a gorilla malaria parasite transmitted by African sylvan anopheline mosquitoes, adapted to humans. Pfs47, a protein on the parasite surface mediates P. falciparum evasion of the mosquito immune system by interacting with a midgut receptor and is critical for Plasmodium adaptation to different anopheline species. Genetic analysis of 4,971 Pfs47 gene sequences from different continents revealed that Asia and Papua New Guinea harbor Pfs47 haplotypes more similar to its ortholog in P. praefalciparum at sites that determine vector compatibility, suggesting that ancestral P. falciparum readily adapted to Asian vectors. Consistent with this observation, Pfs47-receptor gene sequences from African sylvan malaria vectors, such as Anopheles moucheti and An. marshallii, were found to share greater similarity with those of Asian vectors than those of vectors of the African An. gambiae complex. Furthermore, experimental infections provide direct evidence that transformed P. falciparum parasites carrying Pfs47 orthologs of P. praefalciparum or P. reichenowi were more effective at evading the immune system of the Asian malaria vector An. dirus than An. gambiae. We propose that high compatibility of ancestral P. falciparum Pfs47 with the receptors of Asian vectors facilitated the early dispersal of human malaria to the Asian continent, without having to first adapt to sub-Saharan vectors of the An. gambiae complex.

Whole-genome analysis of Malawian Plasmodium falciparum isolates identifies possible targets of allele-specific immunity to clinical malaria.

Individuals acquire immunity to clinical malaria after repeated Plasmodium falciparum infections. Immunity to disease is thought to reflect the acquisition of a repertoire of responses to multiple alleles in diverse parasite antigens. In previous studies, we identified polymorphic sites within individual antigens that are associated with parasite immune evasion by examining antigen allele dynamics in individuals followed longitudinally. Here we expand this approach by analyzing genome-wide polymorphisms using whole genome sequence data from 140 parasite isolates representing malaria cases from a longitudinal study in Malawi and identify 25 genes that encode possible targets of naturally acquired immunity that should be validated immunologically and further characterized for their potential as vaccine candidates.

Parallelized Acquisition of Orbitrap and Astral Analyzers Enables High-Throughput Quantitative Analysis.

The growing trend toward high-throughput proteomics demands rapid liquid chromatography-mass spectrometry (LC-MS) cycles that limit the available time to gather the large numbers of MS/MS fragmentation spectra required for identification. Orbitrap analyzers scale performance with acquisition time and necessarily sacrifice sensitivity and resolving power to deliver higher acquisition rates. We developed a new mass spectrometer that combines a mass-resolving quadrupole, the Orbitrap, and the novel Asymmetric Track Lossless (Astral) analyzer. The new hybrid instrument enables faster acquisition of high-resolution accurate mass (HRAM) MS/MS spectra compared with state-of-the-art mass spectrometers. Accordingly, new proteomics methods were developed that leverage the strengths of each HRAM analyzer, whereby the Orbitrap analyzer performs full scans with a high dynamic range and resolution, synchronized with the Astral analyzer's acquisition of fast and sensitive HRAM MS/MS scans. Substantial improvements are demonstrated over previous methods using current state-of-the-art mass spectrometers.

Can incorporating genotyping data into efficacy estimators improve efficiency of early phase malaria vaccine trials?

BACKGROUND: Early phase malaria vaccine field trials typically measure malaria infection by PCR or thick blood smear microscopy performed on serially sampled blood. Vaccine efficacy (VE) is the proportion reduction in an endpoint due to vaccination and is often calculated as VEHR = 1-hazard ratio or VERR = 1-risk ratio. Genotyping information can distinguish different clones and distinguish multiple infections over time, potentially increasing statistical power. This paper investigates two alternative VE endpoints incorporating genotyping information: VEmolFOI, the vaccine-induced proportion reduction in incidence of new clones acquired over time, and VEC, the vaccine-induced proportion reduction in mean number of infecting clones per exposure. METHODS: Power of VEmolFOI and VEC was compared to that of VEHR and VERR by simulations and analytic derivations, and the four VE methods were applied to three data sets: a Phase 3 trial of RTS,S malaria vaccine in 6912 African infants, a Phase 2 trial of PfSPZ Vaccine in 80 Burkina Faso adults, and a trial comparing Plasmodium vivax incidence in 466 Papua New Guinean children after receiving chloroquine + artemether lumefantrine with or without primaquine (as these VE methods can also quantify effects of other prevention measures). By destroying hibernating liver-stage P. vivax, primaquine reduces subsequent reactivations after treatment completion. RESULTS: In the trial of RTS,S vaccine, a significantly reduced number of clones at first infection was observed, but this was not the case in trials of PfSPZ Vaccine or primaquine, although the PfSPZ trial lacked power to show a reduction. Resampling smaller data sets from the large RTS,S trial to simulate phase 2 trials showed modest power gains from VEC compared to VEHR for data like those from RTS,S, but VEC is less powerful than VEHR for trials in which the number of clones at first infection is not reduced. VEmolFOI was most powerful in model-based simulations, but only the primaquine trial collected enough serial samples to precisely estimate VEmolFOI. The primaquine VEmolFOI estimate decreased after most control arm liver-stage infections reactivated (which mathematically resembles a waning vaccine), preventing VEmolFOI from improving power. CONCLUSIONS: The power gain from the genotyping methods depends on the context. Because input parameters for early phase power calculations are often uncertain, these estimators are not recommended as primary endpoints for small trials unless supported by targeted data analysis. TRIAL REGISTRATIONS: NCT00866619, NCT02663700, NCT02143934.

Gene Coverage Count and Classification (GC3): a locus sequence coverage assessment tool using short-read whole genome sequencing data, and its application to identify and classify histidine-rich protein 2 and 3 deletions in Plasmodium falciparum.

BACKGROUND: The ability of malaria rapid diagnostic tests (RDTs) to effectively detect active infections is being compromised by the presence of malaria strains with genomic deletions at the hrp2 and hrp3 loci, encoding the antigens most commonly targeted in diagnostics for Plasmodium falciparum detection. The presence of such deletions can be determined in publically available P. falciparum whole genome sequencing (WGS) datasets. A computational approach was developed and validated, termed Gene Coverage Count and Classification (GC3), to analyse genome-wide sequence coverage data and provide informative outputs to assess presence and coverage profile of a target locus in WGS data. GC3 was applied to detect deletions at hrp2 and hrp3 (hrp2/3) and flanking genes in different geographic regions and across time points. METHODS: GC3 uses Python and R scripts to extract locus read coverage metrics from mapped WGS data according to user-defined parameters and generates relevant tables and figures. GC3 was tested using WGS data for laboratory reference strains with known hrp2/3 genotypes, and its results compared to those of a hrp2/3-specific qPCR assay. Samples with at least 25% of coding region positions with zero coverage were classified as having a deletion. Publicly available sequence data was analysed and compared with published deletion frequency estimates. RESULTS: GC3 results matched the expected coverage of known laboratory reference strains. Agreement between GC3 and a hrp2/3-specific qPCR assay reported for 19/19 (100%) hrp2 deletions and 18/19 (94.7%) hrp3 deletions. Among Cambodian (n = 127) and Brazilian (n = 20) WGS datasets, which had not been previously analysed for hrp2/3 deletions, GC3 identified hrp2 deletions in three and four samples, and hrp3 deletions in 10 and 15 samples, respectively. Plots of hrp2/3 coding regions, grouped by year of sample collection, showed a decrease in median standardized coverage among Malawian samples (n = 150) suggesting the importance of a careful, properly controlled follow up to determine if an increase in frequency of deletions has occurred between 2007-2008 and 2014-2015. Among Malian (n = 90) samples, median standardized coverage was lower in 2002 than 2010, indicating widespread deletions present at the gene locus in 2002. CONCLUSIONS: The GC3 tool accurately classified hrp2/3 deletions and provided informative tables and figures to analyse targeted gene coverage. GC3 is an appropriate tool when performing preliminary and exploratory assessment of locus coverage data.

Association of a Novel Mutation in the Plasmodium falciparum Chloroquine Resistance Transporter With Decreased Piperaquine Sensitivity.

Background: Amplified copy number in the plasmepsin II/III genes within Plasmodium falciparum has been associated with decreased sensitivity to piperaquine. To examine this association and test whether additional loci might also contribute, we performed a genome-wide association study of ex vivo P. falciparum susceptibility to piperaquine. Methods: Plasmodium falciparum DNA from 183 samples collected primarily from Cambodia was genotyped at 33716 genome-wide single nucleotide polymorphisms (SNPs). Linear mixed models and random forests were used to estimate associations between parasite genotypes and piperaquine susceptibility. Candidate polymorphisms were evaluated for their association with dihydroartemisinin-piperaquine treatment outcomes in an independent dataset. Results: Single nucleotide polymorphisms on multiple chromosomes were associated with piperaquine 90% inhibitory concentrations (IC90) in a genome-wide analysis. Fine-mapping of genomic regions implicated in genome-wide analyses identified multiple SNPs in linkage disequilibrium with each other that were significantly associated with piperaquine IC90, including a novel mutation within the gene encoding the P. falciparum chloroquine resistance transporter, PfCRT. This mutation (F145I) was associated with dihydroartemisinin-piperaquine treatment failure after adjusting for the presence of amplified plasmepsin II/III, which was also associated with decreased piperaquine sensitivity. Conclusions: Our data suggest that, in addition to plasmepsin II/III copy number, other loci, including pfcrt, may also be involved in piperaquine resistance.

Functional analysis of Plasmodium falciparum subpopulations associated with artemisinin resistance in Cambodia.

BACKGROUND: Plasmodium falciparum malaria is one of the most widespread parasitic infections in humans and remains a leading global health concern. Malaria elimination efforts are threatened by the emergence and spread of resistance to artemisinin-based combination therapy, the first-line treatment of malaria. Promising molecular markers and pathways associated with artemisinin drug resistance have been identified, but the underlying molecular mechanisms of resistance remains unknown. The genomic data from early period of emergence of artemisinin resistance (2008-2011) was evaluated, with aim to define k13 associated genetic background in Cambodia, the country identified as epicentre of anti-malarial drug resistance, through characterization of 167 parasite isolates using a panel of 21,257 SNPs. RESULTS: Eight subpopulations were identified suggesting a process of acquisition of artemisinin resistance consistent with an emergence-selection-diffusion model, supported by the shifting balance theory. Identification of population specific mutations facilitated the characterization of a core set of 57 background genes associated with artemisinin resistance and associated pathways. The analysis indicates that the background of artemisinin resistance was not acquired after drug pressure, rather is the result of fixation followed by selection on the daughter subpopulations derived from the ancestral population. CONCLUSIONS: Functional analysis of artemisinin resistance subpopulations illustrates the strong interplay between ubiquitination and cell division or differentiation in artemisinin resistant parasites. The relationship of these pathways with the P. falciparum resistant subpopulation and presence of drug resistance markers in addition to k13, highlights the major role of admixed parasite population in the diffusion of artemisinin resistant background. The diffusion of resistant genes in the Cambodian admixed population after selection resulted from mating of gametocytes of sensitive and resistant parasite populations.

Valine/isoleucine variants drive selective pressure in the VP1 sequence of EV-A71 enteroviruses.

BACKGROUND: In 2011-2012, Northern Vietnam experienced its first large scale hand foot and mouth disease (HFMD) epidemic. In 2011, a major HFMD epidemic was also reported in South Vietnam with fatal cases. This 2011-2012 outbreak was the first one to occur in North Vietnam providing grounds to study the etiology, origin and dynamic of the disease. We report here the analysis of the VP1 gene of strains isolated throughout North Vietnam during the 2011-2012 outbreak and before. METHODS: The VP1 gene of 106 EV-A71 isolates from North Vietnam and 2 from Central Vietnam were sequenced. Sequence alignments were analyzed at the nucleic acid and protein level. Gene polymorphism was also analyzed. A Factorial Correspondence Analysis was performed to correlate amino acid mutations with clinical parameters. RESULTS: The sequences were distributed into four phylogenetic clusters. Three clusters corresponded to the subgenogroup C4 and the last one corresponded to the subgenogroup C5. Each cluster displayed different polymorphism characteristics. Proteins were highly conserved but three sites bearing only Isoleucine (I) or Valine (V) were characterized. The isoleucine/valine variability matched the clusters. Spatiotemporal analysis of the I/V variants showed that all variants which emerged in 2011 and then in 2012 were not the same but were all present in the region prior to the 2011-2012 outbreak. Some correlation was found between certain I/V variants and ethnicity and severity. CONCLUSIONS: The 2011-2012 outbreak was not caused by an exogenous strain coming from South Vietnam or elsewhere but by strains already present and circulating at low level in North Vietnam. However, what triggered the outbreak remains unclear. A selective pressure is applied on I/V variants which matches the genetic clusters. I/V variants were shown on other viruses to correlate with pathogenicity. This should be investigated in EV-A71. I/V variants are an easy and efficient way to survey and identify circulating EV-A71 strains.

Babesia microti from humans and ticks hold a genomic signature of strong population structure in the United States.

BACKGROUND: Babesia microti is an emerging tick-borne apicomplexan parasite with increasing geographic range and incidence in the United States. The rapid expansion of B. microti into its current distribution in the northeastern USA has been due to the range expansion of the tick vector, Ixodes scapularis, upon which the causative agent is dependent for transmission to humans. RESULTS: To reconstruct the history of B. microti in the continental USA and clarify the evolutionary origin of human strains, we used multiplexed hybrid capture of 25 B. microti isolates obtained from I. scapularis and human blood. Despite low genomic variation compared with other Apicomplexa, B. microti was strongly structured into three highly differentiated genetic clusters in the northeastern USA. Bayesian analyses of the apicoplast genomes suggest that the origin of the current diversity of B. microti in northeastern USA dates back 46 thousand years with a signature of recent population expansion in the last 1000 years. Human-derived samples belonged to two rarely intermixing clusters, raising the possibility of highly divergent infectious phenotypes in humans. CONCLUSIONS: Our results validate the multiplexed hybrid capture strategy for characterizing genome-wide diversity and relatedness of B. microti from ticks and humans. We find strong population structure in B. microti samples from the Northeast indicating potential barriers to gene flow.

Plasmodium falciparum parasite population structure and gene flow associated to anti-malarial drugs resistance in Cambodia.

BACKGROUND: Western Cambodia is recognized as the epicentre of emergence of Plasmodium falciparum multi-drug resistance. The emergence of artemisinin resistance has been observed in this area since 2008-2009 and molecular signatures associated to artemisinin resistance have been characterized in k13 gene. At present, one of the major threats faced, is the possible spread of Asian artemisinin resistant parasites over the world threatening millions of people and jeopardizing malaria elimination programme efforts. To anticipate the diffusion of artemisinin resistance, the identification of the P. falciparum population structure and the gene flow among the parasite population in Cambodia are essential. METHODS: To this end, a mid-throughput PCR-LDR-FMA approach based on LUMINEX technology was developed to screen for genetic barcode in 533 blood samples collected in 2010-2011 from 16 health centres in malaria endemics areas in Cambodia. RESULTS: Based on successful typing of 282 samples, subpopulations were characterized along the borders of the country. Each 11-loci barcode provides evidence supporting allele distribution gradient related to subpopulations and gene flow. The 11-loci barcode successfully identifies recently emerging parasite subpopulations in western Cambodia that are associated with the C580Y dominant allele for artemisinin resistance in k13 gene. A subpopulation was identified in northern Cambodia that was associated to artemisinin (R539T resistant allele of k13 gene) and mefloquine resistance. CONCLUSIONS: The gene flow between these subpopulations might have driven the spread of artemisinin resistance over Cambodia.

Can incorporating genotyping data into efficacy estimators improve efficiency of early phase malaria vaccine trials?

Background: Early phase malaria vaccine field trials typically measure malaria infection by PCR or thick blood smear microscopy performed on serially sampled blood. Vaccine efficacy (VE) is the proportion reduction in an endpoint due to vaccination and is often calculated as VEHR=1 - hazard ratio or VERR=1 - risk ratio. Genotyping information can distinguish different clones and distinguish multiple infections over time, potentially increasing statistical power. This paper investigates two alternative VE endpoints incorporating genotyping information: VEmolFOI, the vaccine-induced proportion reduction in incidence of new clones acquired over time, and VEC, the vaccine-induced proportion reduction in mean number of infecting clones per exposure. Methods: We used simulations and analytic derivations to compare power of these methods to VEHR and VERR and applied them to three data sets: a Phase 3 trial of RTS,S malaria vaccine in 6912 African infants, a Phase 2 trial of PfSPZ Vaccine in 80 Burkina Faso adults, and a trial comparing Plasmodium vivax incidence in 466 Papua New Guinean children after receiving chloroquine + artemether lumefantrine with or without primaquine (as these VE methods can also quantify effects of other prevention measures). By destroying hibernating liver-stage P. vivax, primaquine reduces subsequent reactivations after treatment completion. Results: The RTS,S vaccine significantly reduced the number of clones at first infection, but PfSPZ vaccine and primaquine did not. Resampling smaller data sets from the large RTS,S trial to simulate phase 2 trials showed modest power gains from VEC compared to VEHR for data like RTS,S, but VEC is less powerful than VEHR for vaccines which do not reduce the number of clones at first infection. VEmolFOI was most powerful in model-based simulations, but only the primaquine trial collected enough serial samples to precisely estimate VEmolFOI. The primaquine VEmolFOI estimate decreased after most control arm liver-stage infections reactivated (which mathematically resembles a waning vaccine), preventing VEmolFOI from improving power. Conclusions: The power gain from the genotyping methods depends on the context. Because input parameters for early phase power calculations are often uncertain, we recommend against these estimators as primary endpoints for small trials unless supported by targeted data analysis. Trial registrations: NCT00866619, NCT02663700, NCT02143934.

Immunogenomic profile at baseline predicts host susceptibility to clinical malaria.

Introduction: Host gene and protein expression impact susceptibility to clinical malaria, but the balance of immune cell populations, cytokines and genes that contributes to protection, remains incompletely understood. Little is known about the determinants of host susceptibility to clinical malaria at a time when acquired immunity is developing. Methods: We analyzed peripheral blood mononuclear cells (PBMCs) collected from children who differed in susceptibility to clinical malaria, all from a small town in Mali. PBMCs were collected from children aged 4-6 years at the start, peak and end of the malaria season. We characterized the immune cell composition and cytokine secretion for a subset of 20 children per timepoint (10 children with no symptomatic malaria age-matched to 10 children with >2 symptomatic malarial illnesses), and gene expression patterns for six children (three per cohort) per timepoint. Results: We observed differences between the two groups of children in the expression of genes related to cell death and inflammation; in particular, inflammatory genes such as CXCL10 and STAT1 and apoptotic genes such as XAF1 were upregulated in susceptible children before the transmission season began. We also noted higher frequency of HLA-DR+ CD4 T cells in protected children during the peak of the malaria season and comparable levels cytokine secretion after stimulation with malaria schizonts across all three time points. Conclusion: This study highlights the importance of baseline immune signatures in determining disease outcome. Our data suggests that differences in apoptotic and inflammatory gene expression patterns can serve as predictive markers of susceptibility to clinical malaria.

An In Silico Analysis of Malaria Pre-Erythrocytic-Stage Antigens Interpreting Worldwide Genetic Data to Suggest Vaccine Candidate Variants and Epitopes.

Failure to account for genetic diversity of antigens during vaccine design may lead to vaccine escape. To evaluate the vaccine escape potential of antigens used in vaccines currently in development or clinical testing, we surveyed the genetic diversity, measured population differentiation, and performed in silico prediction and analysis of T-cell epitopes of ten such Plasmodium falciparum pre-erythrocytic-stage antigens using whole-genome sequence data from 1010 field isolates. Of these, 699 were collected in Africa (Burkina Faso, Cameroon, Guinea, Kenya, Malawi, Mali, and Tanzania), 69 in South America (Brazil, Colombia, French Guiana, and Peru), 59 in Oceania (Papua New Guinea), and 183 in Asia (Cambodia, Myanmar, and Thailand). Antigens surveyed include cell-traversal protein for ookinetes and sporozoites, circumsporozoite protein, liver-stage antigens 1 and 3, sporozoite surface proteins P36 and P52, sporozoite asparagine-rich protein-1, sporozoite microneme protein essential for cell traversal-2, and upregulated-in-infectious-sporozoite 3 and 4 proteins. The analyses showed that a limited number of these protein variants, when combined, would be representative of worldwide parasite populations. Moreover, predicted T-cell epitopes were identified that could be further explored for immunogenicity and protective efficacy. Findings can inform the rational design of a multivalent malaria vaccine.

A genetically modified Plasmodium berghei parasite as a surrogate for whole-sporozoite vaccination against P. vivax malaria.

Two malaria parasite species, Plasmodium falciparum (Pf) and P. vivax (Pv) are responsible for most of the disease burden caused by malaria. Vaccine development against this disease has focused mainly on Pf. Whole-sporozoite (WSp) vaccination, targeting pre-erythrocytic (PE) parasite stages, is a promising strategy for immunization against malaria and several PfWSp-based vaccine candidates are currently undergoing clinical evaluation. In contrast, no WSp candidates have been developed for Pv, mainly due to constraints in the production of Pv sporozoites in the laboratory. Recently, we developed a novel approach for WSp vaccination against Pf based on the use of transgenic rodent P. berghei (Pb) sporozoites expressing immunogens of this human-infective parasite. We showed that this platform can be used to deliver PE Pf antigens, eliciting both targeted humoral responses and cross-species cellular immune responses against Pf. Here we explored this WSp platform for the delivery of Pv antigens. As the Pv circumsporozoite protein (CSP) is a leading vaccine candidate antigen, we generated a transgenic Pb parasite, PbviVac, that, in addition to its endogenous PbCSP, expresses PvCSP under the control of a strictly PE promoter. Immunofluorescence microscopy analyses confirmed that both the PbCSP and the PvCSP antigens are expressed in PbviVac sporozoites and liver stages and that PbviVac sporozoite infectivity of hepatic cells is similar to that of its wild-type Pb counterpart. Immunization of mice with PbviVac sporozoites elicits the production of anti-PvCSP antibodies that efficiently recognize and bind to Pv sporozoites. Our results warrant further development and evaluation of PbviVac as a surrogate for WSp vaccination against Pv malaria.

Plasmodium falciparum 7G8 challenge provides conservative prediction of efficacy of PfNF54-based PfSPZ Vaccine in Africa.

Controlled human malaria infection (CHMI) has supported Plasmodium falciparum (Pf) malaria vaccine development by providing preliminary estimates of vaccine efficacy (VE). Because CHMIs generally use Pf strains similar to vaccine strains, VE against antigenically heterogeneous Pf in the field has been required to establish VE. We increased the stringency of CHMI by selecting a Brazilian isolate, Pf7G8, which is genetically distant from the West African parasite (PfNF54) in our PfSPZ vaccines. Using two regimens to identically immunize US and Malian adults, VE over 24 weeks in the field was as good as or better than VE against CHMI at 24 weeks in the US. To explain this finding, here we quantify differences in the genome, proteome, and predicted CD8 T cell epitopes of PfNF54 relative to 704 Pf isolates from Africa and Pf7G8. We show that Pf7G8 is more distant from PfNF54 than any African isolates tested. We propose VE against Pf7G8 CHMI for providing pivotal data for malaria vaccine licensure for travelers to Africa, and potentially for endemic populations, because the genetic distance of Pf7G8 from the Pf vaccine strain makes it a stringent surrogate for Pf parasites in Africa.

Integration of population and functional genomics to understand mechanisms of artemisinin resistance in Plasmodium falciparum.

Resistance to antimalarial drugs, and in particular to the artemisinin derivatives and their partner drugs, threatens recent progress toward regional malaria elimination and eventual global malaria eradication. Population-level studies utilizing whole-genome sequencing approaches have facilitated the identification of regions of the parasite genome associated with both clinical and in vitro drug-resistance phenotypes. However, the biological relevance of genes identified in these analyses and the establishment of a causal relationship between genotype and phenotype requires functional characterization. Here we examined data from population genomic and transcriptomic studies in the context of data generated from recent functional studies, using a new population genetic approach designed to identify potential favored mutations within the region of a selective sweep (iSAFE). We identified several genes functioning in pathways now known to be associated with artemisinin resistance that were supported in early population genomic studies, as well as potential new drug targets/pathways for further validation and consideration for treatment of artemisinin-resistant Plasmodium falciparum. In addition, we establish the utility of iSAFE in identifying positively-selected mutations in population genomic studies, potentially accelerating the time to functional validation of candidate genes.

Oral Selective TLR8 Agonist Selgantolimod Induces Multiple Immune Cell Responses in Humans.

TLR8 agonists have the potential for use as immunomodulatory components in therapeutic modalities for viral infections such as chronic HBV (CHB) and HIV. In this study, using peripheral blood samples from a phase 1a clinical trial, we examined the acute effects of a single oral administration of a selective TLR8 agonist on immune cell phenotypes. Administration of the TLR8 agonist selgantolimod (SLGN) in healthy individuals resulted in alteration in frequencies of peripheral blood monocytes, pDCs, mDCs and MAIT cells. Frequencies of mDCs and lymphoid cells significantly reduced after 8 h of SLGN administration, whereas pDC frequencies significantly increased, with changes possibly reflecting migration of different cell types between peripheral and tissue compartments in response to the agonist. Myeloid cell activation was evident by an upregulated expression of co-stimulatory molecules CD40 and CD86 accompanied by the production of IL-6 and IL-18 from these cells. Concomitantly, there was induction of the early activation marker CD69 on innate and adaptive lymphoid cells, including MAIT and NK cell subsets. Further, these activated lymphoid cells had enhanced expression of the effector molecules granzyme B and perforin. Microarray analysis of isolated lymphocytes and monocytes from baseline and post-SLGN treatment revealed changes in expression of genes involved in cellular response to cytokine stimulus, innate immune response, myeloid cell differentiation and antigen receptor-mediated signaling pathway. In a preliminary analysis of samples from CHB patients treated with selgantolimod, activation of innate and adaptive lymphocytes was evident. In conclusion, this first in-human study shows that selgantolimod administration in humans results in activation of multiple immune cell responses with antiviral potential.

A genotyping assay to determine geographic origin and transmission potential of Plasmodium falciparum malaria cases.

As countries work towards malaria elimination, it is important to monitor imported cases to prevent reestablishment of local transmission. The Plasmodium falciparum Pfs47 gene has strong geographic population structure, because only those parasites with Pfs47 haplotypes compatible with the mosquito vector species in a given continent are efficiently transmitted. Analysis of 4,971 world-wide Pfs47 sequences identified two SNPs (at 707 and 725 bp) as sufficient to establish the likely continent of origin of P. falciparum isolates. Pfs47 sequences from Africa, Asia, and the New World presented more that 99% frequency of distinct combinations of the SNPs 707 and 725 genotypes. Interestingly, Papua New Guinea Pfs47 sequences have the highest diversity in SNPs 707 and 725. Accurate and reproducible High-Resolution Melting (HRM) assays were developed to genotype Pfs47 SNPs 707 and 725 in laboratory and field samples, to assess the geographic origin and risk of local transmission of imported P. falciparum malaria cases.

Successful Profiling of Plasmodium falciparum var Gene Expression in Clinical Samples via a Custom Capture Array.

var genes encode Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) antigens. These highly diverse antigens are displayed on the surface of infected erythrocytes and play a critical role in immune evasion and sequestration of infected erythrocytes. Studies of var expression using non-leukocyte-depleted blood are challenging because of the predominance of host genetic material and lack of conserved var segments. Our goal was to enrich for parasite RNA, allowing de novo assembly of var genes and detection of expressed novel variants. We used two overall approaches: (i) enriching for total mRNA in the sequencing library preparations and (ii) enriching for parasite RNA with a custom capture array based on Roche's SeqCap EZ enrichment system. The capture array was designed with probes based on the whole 3D7 reference genome and an additional >4,000 full-length var gene sequences from other P. falciparum strains. We tested each method on RNA samples from Malian children with severe or uncomplicated malaria infections. All reads mapping to the human genome were removed, the remaining reads were assembled de novo into transcripts, and from these, var-like transcripts were identified and annotated. The capture array produced the longest maximum length and largest numbers of var gene transcripts in each sample, particularly in samples with low parasitemia. Identifying the most-expressed var gene sequences in whole-blood clinical samples without the need for extensive processing or generating sample-specific reference genome data is critical for understanding the role of PfEMP1s in malaria pathogenesis. IMPORTANCE Malaria parasites display antigens on the surface of infected red blood cells in the human host that facilitate attachment to blood vessels, contributing to the severity of infection. These antigens are highly variable, allowing the parasite to evade the immune system. Identifying these expressed antigens is critical to understanding the development of severe malarial disease. However, clinical samples contain limited amounts of parasite genetic material, a challenge for sequencing efforts further compounded by the extreme diversity of the parasite surface antigens. We present a method that enriches for these antigen sequences in clinical samples using a custom capture array, requiring minimal processing in the field. While our results are focused on the malaria parasite Plasmodium falciparum, this approach has broad applicability to other highly diverse antigens from other parasites and pathogens such as those that cause giardiasis and leishmaniasis.

Strains used in whole organism Plasmodium falciparum vaccine trials differ in genome structure, sequence, and immunogenic potential.

BACKGROUND: Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions. METHODS: Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia. RESULTS: While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8+ T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region. CONCLUSIONS: These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.