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

Plasmodium vivax Infections in Duffy-Negative Individuals: A Paradigm Shift in Indian Malaria Epidemiology

Background

To investigate the occurrence of Plasmodium vivax infections in Duffy-negative individuals, challenging the long-held belief that P. vivax requires the Duffy antigen receptor for chemokines to infect human erythrocytes.

Materials and Methods

In the present study, 365 samples were screened using serological techniques, PCR-RFLP analysis, and DNA sequencing of the ACKR1 gene promoter region mutation to identify Duffy-negative individuals. P. vivax infection was detected using PCR targeting the 18S rRNA gene and microscopic examination of Giemsa-stained blood smears.

Results

Five individuals (1.36%) were confirmed Duffy-negative (Fy(a-b-)). Surprisingly, 3 out of these 5 Duffy-negative subjects (60%) were infected with P. vivax, as confirmed by both microscopy and PCR. Various parasite stages were observed in infected Duffy-negative samples, with parasitaemia ranging from 0.01% to 0.5%.

Discussion

Our findings provide compelling evidence that P. vivax can infect Duffy-negative individuals, suggesting the existence of alternative invasion pathways or adaptations. This has profound implications for P. vivax biology, evolution, and global distribution. The burden of vivax malaria may be underestimated, particularly in regions with a high prevalence of Duffy negativity. This study highlights the need to reevaluate P. vivax epidemiology, diagnostic approaches, and control strategies, especially in areas previously considered at low risk. Further research is needed to elucidate the mechanisms enabling P. vivax invasion of Duffy-negative erythrocytes and to assess the clinical and epidemiological consequences of these infections.

Genetic surveillance of Plasmodium-Anopheles compatibility markers during Anopheles stephensi associated malaria outbreak

Despite previous decline of malaria in Ethiopia, an outbreak in Dire Dawa in 2022 implicated invasive vector An. stephensi as responsible. The transmission of Plasmodium by invasive An. stephensi raises questions about the molecular basis of compatibility, and the origin of the Plasmodium being transmitted. The Plasmodium P47 gene is a parasite-vector interaction gene in Anopheles , and along with corresponding mosquito P47 receptor ( P47Rec ), can be critical in establishment of Plasmodium infections in anophelines. Here, we analyzed P47 and P47Rec sequences to determine the origin of Plasmodium detected in An. stephensi during the outbreak and evaluate markers of compatibility. Analysis of geographically informative SNPs in Pfs47 revealed that these P. falciparum exhibit the African haplotype. We also identified a single amino acid change in P47Rec within these An. stephensi, which could act as a marker for the propensity of An. stephensi populations to outbreaks. Together, we provide the basis for further study to deepen the understanding of invasive An. stephensi -African Plasmodium interactions to better control transmission of malaria and prevent further outbreaks.

Intra-individual variability in ancient plasmodium DNA recovery highlights need for enhanced sampling

Malaria has been a leading cause of death in human populations for centuries and remains a major public health challenge in African countries, especially affecting children. Among the five Plasmodium species infecting humans, Plasmodium falciparum is the most lethal. Ancient DNA research has provided key insights into the origins, evolution, and virulence of pathogens that affect humans. However, extensive screening of ancient skeletal remains for Plasmodium DNA has shown that such genomic material is rare, with no studies so far addressing potential intra-individual variability. Consequently, the pool of ancient mitochondrial DNA (mtDNA) or genomic sequences for P. falciparum is extremely limited, with fewer than 20 ancient sequences available for genetic analysis, and no complete P. falciparum mtDNA from Classical antiquity published to date. To investigate intra-individual diversity and genetic origins of P. falciparum from the Roman period, we generated 39 sequencing libraries from multiple teeth and two from the femur of a Roman malaria-infected individual. The results revealed considerable variability in P. falciparum recovery across different dental samples within the individual, while the femur samples showed no preservation of Plasmodium DNA. The reconstructed 43-fold P. falciparum mtDNA genome supports the hypothesis of an Indian origin for European P. falciparum and suggests mtDNA continuity in Europe over the past 2000 years.

Genomic exploration of the journey of Plasmodium vivax in Latin America

Plasmodium vivax is the predominant malaria parasite in Latin America. Its colonization history in the region is rich and complex, and is still highly debated, especially about its origin(s). Our study employed cutting-edge population genomic techniques to analyze whole genome variation from 620 P. vivax isolates, including 107 newly sequenced samples from West Africa, Middle East, and Latin America. This sampling represents nearly all potential source populations worldwide currently available. Analyses of the genetic structure, diversity, ancestry, coalescent-based inferences, including demographic scenario testing using Approximate Bayesian Computation, have revealed a more complex evolutionary history than previously envisioned. Indeed, our analyses suggest that the current American P. vivax populations predominantly stemmed from a now-extinct European lineage, with the potential contribution also from unsampled populations, most likely of West African origin. We also found evidence that P. vivax arrived in Latin America in multiple waves, initially during early European contact and later through post-colonial human migration waves in the late 19th-century. This study provides a fresh perspective on P. vivax's intricate evolutionary journey and brings insights into the possible contribution of West African P. vivax populations to the colonization history of Latin America.

Genomic analysis of global Plasmodium vivax populations reveals insights into the evolution of drug resistance

Increasing reports of chloroquine resistance (CQR) in Plasmodium vivax endemic regions have led to several countries, including Indonesia, to adopt dihydroarteminsin-piperaquine instead. However, the molecular drivers of CQR remain unclear. Using a genome-wide approach, we perform a genomic analysis of 1534 P. vivax isolates across 29 endemic countries, detailing population structure, patterns of relatedness, selection, and resistance profiling, providing insights into potential drivers of CQR. Selective sweeps in a locus proximal to pvmdr1, a putative marker for CQR, along with transcriptional regulation genes, distinguish isolates from Indonesia from those in regions where chloroquine remains highly effective. In 106 isolates from Indonesian Papua, the epicentre of CQR, we observe an increasing prevalence of novel SNPs in the candidate resistance gene pvmrp1 since the introduction of dihydroartemisinin-piperaquine. Overall, we provide novel markers for resistance surveillance, supported by evidence of regions under recent directional selection and temporal analysis in this continually evolving parasite.

Turning the needle into the haystack: Culture-independent amplification of complex microbial genomes directly from their native environment

High-throughput sequencing (HTS) has revolutionized microbiology, but many microbes exist at low abundance in their natural environment and/or are difficult, if not impossible, to culture in the laboratory. This makes it challenging to use HTS to study the genomes of many important microbes and pathogens. In this review, we discuss the development and application of selective whole genome amplification (SWGA) to allow whole or partial genomes to be sequenced for low abundance microbes directly from complex biological samples. We highlight ways in which genomic data generated by SWGA have been used to elucidate the population dynamics of important human pathogens and monitor development of antimicrobial resistance and the emergence of potential outbreaks. We also describe the limitations of this method and propose some potential innovations that could be used to improve the quality of SWGA and lower the barriers to using this method across a wider range of infectious pathogens.

Origin of the human malaria parasite Plasmodium vivax

The geographic origin of Plasmodium vivax, a leading cause of human malaria, has been the subject of much speculation. Here we review the evolutionary history of P. vivax and P. vivax-like parasites in humans and non-human primates on three continents, providing overwhelming evidence for an African origin. This conclusion is consistent with recent reports showing that Duffy-negative humans in Africa are, in fact, susceptible to P. vivax, with parasites invading Duffy-antigen-expressing erythroid precursors. Thus, the African origin of P. vivax not only explains the distribution of the Duffy-negative genotype but also provides new insight into the history and status of P. vivax malaria in Africa and efforts geared toward its eradication.

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

Plasmodium vivax genomic surveillance in the Peruvian Amazon with Pv AmpliSeq assay

BACKGROUND

Plasmodium vivax is the most predominant malaria species in Latin America, constituting 71.5% of malaria cases in 2021. With several countries aiming for malaria elimination, it is crucial to prioritize effectiveness of national control programs by optimizing the utilization of available resources and strategically implementing necessary changes. To support this, there is a need for innovative approaches such as genomic surveillance tools that can investigate changes in transmission intensity, imported cases and sources of reintroduction, and can detect molecular markers associated with drug resistance.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we apply a modified highly-multiplexed deep sequencing assay: Pv AmpliSeq v2 Peru. The tool targets a newly developed 41-SNP Peru barcode for parasite population analysis within Peru, the 33-SNP vivaxGEN-geo panel for country-level classification, and 11 putative drug resistance genes. It was applied to 230 samples from the Peruvian Amazon (2007-2020), generating baseline surveillance data. We observed a heterogenous P. vivax population with high diversity and gene flow in peri-urban areas of Maynas province (Loreto region) with a temporal drift using all SNPs detected by the assay (nSNP = 2909). In comparison, in an indigenous isolated area, the parasite population was genetically differentiated (FST = 0.07-0.09) with moderate diversity and high relatedness between isolates in the community. In a remote border community, a clonal P. vivax cluster was identified, with distinct haplotypes in drug resistant genes and ama1, more similar to Brazilian isolates, likely representing an introduction of P. vivax from Brazil at that time. To test its applicability for Latin America, we evaluated the SNP Peru barcode in P. vivax genomes from the region and demonstrated the capacity to capture local population clustering at within-country level.

CONCLUSIONS/SIGNIFICANCE

Together this data shows that P. vivax transmission is heterogeneous in different settings within the Peruvian Amazon. Genetic analysis is a key component for regional malaria control, offering valuable insights that should be incorporated into routine surveillance.

Analytic Approaches in Genomic Epidemiological Studies of Parasitic Protozoa

Whole genome sequencing (WGS) plays an important role in the advanced characterization of pathogen transmission and is widely used in studies of major bacterial and viral diseases. Although protozoan parasites cause serious diseases in humans and animals, WGS data on them are relatively scarce due to the large genomes and lack of cultivation techniques for some. In this review, we have illustrated bioinformatic analyses of WGS data and their applications in studies of the genomic epidemiology of apicomplexan parasites. WGS has been used in outbreak detection and investigation, studies of pathogen transmission and evolution, and drug resistance surveillance and tracking. However, comparative analysis of parasite WGS data is still in its infancy, and available WGS data are mainly from a few genera of major public health importance, such as Plasmodium, Toxoplasma, and Cryptosporidium. In addition, the utility of third-generation sequencing technology for complete genome assembly at the chromosome level, studies of the biological significance of structural genomic variation, and molecular surveillance of pathogens has not been fully exploited. These issues require large-scale WGS of various protozoan parasites of public health and veterinary importance using both second- and third-generation sequencing technologies.

Population genomic evidence of structured and connected Plasmodium vivax populations under host selection in Latin America

Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub-patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high-quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub-populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility - crucial for the parasite's survival in diverse environments. Understanding these population-level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

Vivax malaria: a possible stumbling block for malaria elimination in India

Plasmodium vivax is geographically the most widely dispersed human malaria parasite species. It has shown resilience and a great deal of adaptability. Genomic studies suggest that P. vivax originated from Asia or Africa and moved to the rest of the world. Although P. vivax is evolutionarily an older species than Plasmodium falciparum, its biology, transmission, pathology, and control still require better elucidation. P. vivax poses problems for malaria elimination because of the ability of a single primary infection to produce multiple relapses over months and years. P. vivax malaria elimination program needs early diagnosis, and prompt and complete radical treatment, which is challenging, to simultaneously exterminate the circulating parasites and dormant hypnozoites lodged in the hepatocytes of the host liver. As prompt surveillance and effective treatments are rolled out, preventing primaquine toxicity in the patients having glucose-6-phosphate dehydrogenase (G6PD) deficiency should be a priority for the vivax elimination program. This review sheds light on the burden of P. vivax, changing epidemiological patterns, the hurdles in elimination efforts, and the essential tools needed not just in India but globally. These tools encompass innovative treatments for eliminating dormant parasites, coping with evolving drug resistance, and the development of potential vaccines against the parasite.

Comparative transcriptomics reveal differential gene expression among Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms

The documentation of Plasmodium vivax malaria across Africa especially in regions where Duffy negatives are dominant suggests possibly alternative erythrocyte invasion mechanisms. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa is unclear. In this study, we examined the expression of 4,404 gene transcripts belong to 12 functional groups and 43 erythrocyte binding gene candidates in Ethiopian isolates and compared them with the Cambodian and Brazilian P. vivax transcriptomes. Overall, there were 10-26% differences in the gene expression profile amongst geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian isolates. Members of the reticulocyte binding protein PvRBP2a and PvRBP3 expressed six-fold higher than Duffy binding protein PvDBP1 and 60-fold higher than PvEBP/DBP2 in the Ethiopian isolates. Other genes including PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression patterns were observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, gametocyte genes including PvAP2-G, PvGAP (female gametocytes), and Pvs47 (male gametocytes) were highly expressed across geographical samples.

Genomic analysis of Plasmodium vivax describes patterns of connectivity and putative drivers of adaptation in Ethiopia

Ethiopia has the greatest burden of Plasmodium vivax in Africa, but little is known about the epidemiological landscape of parasites across the country. We analysed the genomic diversity of 137 P. vivax isolates collected nine Ethiopian districts from 2012 to 2016. Signatures of selection were detected by cross-country comparisons with isolates from Thailand (n = 104) and Indonesia (n = 111), representing regions with low and high chloroquine resistance respectively. 26% (35/137) of Ethiopian infections were polyclonal, and 48.5% (17/35) of these comprised highly related clones (within-host identity-by-descent > 25%), indicating frequent co-transmission and superinfection. Parasite gene flow between districts could not be explained entirely by geographic distance, with economic and cultural factors hypothesised to have an impact on connectivity. Amplification of the duffy binding protein gene (pvdbp1) was prevalent across all districts (16-75%). Cross-population haplotype homozygosity revealed positive selection in a region proximal to the putative chloroquine resistance transporter gene (pvcrt-o). An S25P variant in amino acid transporter 1 (pvaat1), whose homologue has recently been implicated in P. falciparum chloroquine resistance evolution, was prevalent in Ethiopia (96%) but not Thailand or Indonesia (35-53%). The genomic architecture in Ethiopia highlights circulating variants of potential public health concern in an endemic setting with evidence of stable transmission.

Genomics of Plasmodium vivax in Colombia reveals evidence of local bottle-necking and inter-country connectivity in the Americas

Colombia aims to eliminate malaria by 2030 but remains one of the highest burden countries in the Americas. Plasmodium vivax contributes half of all malaria cases, with its control challenged by relapsing parasitaemia, drug resistance and cross-border spread. Using 64 Colombian P. vivax genomes collected between 2013 and 2017, we explored diversity and selection in two major foci of transmission: Chocó and Córdoba. Open-access data from other countries were used for comparative assessment of drug resistance candidates and to assess cross-border spread. Across Colombia, polyclonal infections were infrequent (12%), and infection connectivity was relatively high (median IBD = 5%), consistent with low endemicity. Chocó exhibited a higher frequency of polyclonal infections (23%) than Córdoba (7%), although the difference was not significant (P = 0.300). Most Colombian infections carried double pvdhfr (95%) and single pvdhps (71%) mutants, but other drug resistance mutations were less prevalent (< 10%). There was no evidence of selection at the pvaat1 gene, whose P. falciparum orthologue has recently been implicated in chloroquine resistance. Global population comparisons identified other putative adaptations. Within the Americas, low-level connectivity was observed between Colombia and Peru, highlighting potential for cross-border spread. Our findings demonstrate the potential of molecular data to inform on infection spread and adaptation.

A new Plasmodium vivax reference genome for South American isolates

BACKGROUND

Plasmodium vivax is the second most important cause of human malaria worldwide, and accounts for the majority of malaria cases in South America. A high-quality reference genome exists for Papua Indonesia (PvP01) and Thailand (PvW1), but is lacking for South America. A reference genome specifically for South America would be beneficial though, as P. vivax is a genetically diverse parasite with geographical clustering.

RESULTS

This study presents a new high-quality assembly of a South American P. vivax isolate, referred to as PvPAM (P. vivax Peruvian AMazon). The genome was obtained from a low input patient sample from the Peruvian Amazon and sequenced using PacBio technology, resulting in a highly complete assembly with 6497 functional genes. Telomeric ends were present in 17 out of 28 chromosomal ends, and additional (sub)telomeric regions are present in 12 unassigned contigs. A comparison of multigene families between PvPAM and the PvP01 genome revealed remarkable variation in vir genes, and the presence of merozoite surface proteins (MSP) 3.6 and 3.7. Three dhfr and dhps drug resistance associated mutations are present in PvPAM, similar to those found in other Peruvian isolates. Mapping of publicly available South American whole genome sequencing (WGS) data to PvPAM resulted in significantly fewer variants and truncated reads compared to the use of PvP01 or PvW1 as reference genomes. To minimize the number of core genome variants in non-South American samples, PvW1 is most suited for Southeast Asian isolates, both PvPAM and PvW1 are suited for South Asian isolates, and PvPAM is recommended for African isolates. Interestingly, non-South American samples still contained the least subtelomeric variants when mapped to PvPAM, indicating high quality of the PvPAM subtelomeric regions.

CONCLUSIONS

Our findings show that the PvPAM reference genome more accurately represents South American P. vivax isolates in comparison to PvP01 and PvW1. In addition, PvPAM has a high level of completeness, and contains a similar number of annotated genes as PvP01 or PvW1. The PvPAM genome therefore will be a valuable resource to improve future genomic analyses on P. vivax isolates from the South American continent.

Genetic diversity and molecular evolution of Plasmodium vivax Duffy Binding Protein and Merozoite Surface Protein-1 in northwestern Thailand

The local diversity and population structure of malaria parasites vary across different regions of the world, reflecting variations in transmission intensity, host immunity, and vector species. This study aimed to use amplicon sequencing to investigate the genotypic patterns and population structure of P. vivax isolates from a highly endemic province of Thailand in recent years. Amplicon deep sequencing was performed on 70 samples for the 42-kDa region of pvmsp1 and domain II of pvdbp. Unique haplotypes were identified and a network constructed to illustrate genetic relatedness in northwestern Thailand. Based on this dataset of 70 samples collected between 2015 and 2021, 16 and 40 unique haplotypes were identified in pvdbpII and pvmsp142kDa, respectively. Nucleotide diversity was higher in pvmsp142kDa than in pvdbpII (π = 0.027 and 0.012), as was haplotype diversity (Hd = 0.962 and 0.849). pvmsp142kDa also showed a higher recombination rate and higher levels of genetic differentiation (Fst) in northwestern Thailand versus other regions (0.2761-0.4881). These data together suggested that the genetic diversity of P. vivax in northwestern Thailand at these two studied loci evolved under a balancing selection, most likely host immunity. The lower genetic diversity of pvdbpII may reflect its stronger functional constrain. In addition, despite the balancing selection, a decrease in genetic diversity was observed. Hd of pvdbpII decreased from 0.874 in 2015-2016 to 0.778 in 2018-2021; π of pvmsp142kDa decreased from 0.030 to 0.022 over the same period. Thus, the control activities must have had a strong impact on the parasite population size. The findings from this study provide an understanding of P. vivax population structure and the evolutionary force on vaccine candidates. They also established a new baseline for tracking future changes in P. vivax diversity in the most malarious area of Thailand.

How can the complex epidemiology of malaria in India impact its elimination?

Malaria is a human health hazard in the tropical and subtropical zones of the globe and is poised to be eliminated by the year 2030. Despite a decrease in incidence in the past two decades, many endemic countries, including India, report cases regularly. The epidemiology of malaria in India is unique owing to several features of the Plasmodium parasites, Anopheles vectors, ecoepidemiological situations conducive to disease transmission, and susceptible humans living in rural and forested areas. Limitations in public health reach, and poor health-seeking behaviour of vulnerable populations living in hard-to-reach areas, add to the problem. We bring all of these factors together in a comprehensive framework and opine that, in spite of complexities, targeted elimination of malaria in India is achievable with planned programmatic approaches.

Etiology of acute febrile illness in the peruvian amazon as determined by modular formatted quantitative PCR: a protocol for RIVERA, a health facility-based case-control study

BACKGROUND

The study of the etiology of acute febrile illness (AFI) has historically been designed as a prevalence of pathogens detected from a case series. This strategy has an inherent unrealistic assumption that all pathogen detection allows for causal attribution, despite known asymptomatic carriage of the principal causes of acute febrile illness in most low- and middle-income countries (LMICs). We designed a semi-quantitative PCR in a modular format to detect bloodborne agents of acute febrile illness that encompassed common etiologies of AFI in the region, etiologies of recent epidemics, etiologies that require an immediate public health response and additional pathogens of unknown endemicity. We then designed a study that would delineate background levels of transmission in the community in the absence of symptoms to provide corrected estimates of attribution for the principal determinants of AFI.

METHODS

A case-control study of acute febrile illness in patients ten years or older seeking health care in Iquitos, Loreto, Peru, was planned. Upon enrollment, we will obtain blood, saliva, and mid-turbinate nasal swabs at enrollment with a follow-up visit on day 21-28 following enrollment to attain vital status and convalescent saliva and blood samples, as well as a questionnaire including clinical, socio-demographic, occupational, travel, and animal contact information for each participant. Whole blood samples are to be simultaneously tested for 32 pathogens using TaqMan array cards. Mid-turbinate samples will be tested for SARS-CoV-2, Influenza A and Influenza B. Conditional logistic regression models will be fitted treating case/control status as the outcome and with pathogen-specific sample positivity as predictors to attain estimates of attributable pathogen fractions for AFI.

DISCUSSION

The modular PCR platforms will allow for reporting of all primary results of respiratory samples within 72 h and blood samples within one week, allowing for results to influence local medical practice and enable timely public health responses. The inclusion of controls will allow for a more accurate estimate of the importance of specific prevalent pathogens as a cause of acute illness.

STUDY REGISTRATION

Project 1791, Registro de Proyectos de Investigación en Salud Pública (PRISA), Instituto Nacional de Salud, Perú.

Etiology of Acute Febrile Illness in the Peruvian Amazon as determined by modular formatted quantitative PCR: A Protocol for RIVERA, a Health Facility-Based Case-Control Study

Background

The study of the etiology of acute febrile illness (AFI) has historically been designed as a prevalence of pathogens detected from a case series. This strategy has an inherent unrealistic assumption that all pathogen detection allows for causal attribution, despite known asymptomatic carriage of the principal causes of acute febrile illness in most low- and middle-income countries (LMICs). We designed a semi-quantitative PCR in a modular format to detect bloodborne agents of acute febrile illness that encompassed common etiologies of AFI in the region, etiologies of recent epidemics, etiologies that require an immediate public health response and additional pathogens of unknown endemicity. We then designed a study that would delineate background levels of transmission in the community in the absence of symptoms to provide corrected estimates of attribution for the principal determinants of AFI.

Methods

A case-control study of acute febrile illness in patients ten years or older seeking health care in Iquitos, Loreto, Peru, was planned. Upon enrollment, we will obtain blood, saliva, and mid-turbinate nasal swabs at enrollment with a follow-up visit on day 21-28 following enrollment to attain vital status and convalescent saliva and blood samples, as well as a questionnaire including clinical, socio-demographic, occupational, travel, and animal contact information for each participant. Whole blood samples are to be simultaneously tested for 32 pathogens using TaqMan array cards. Mid-turbinate samples will be tested for SARS-CoV-2, Influenza A and Influenza B. Conditional logistic regression models will be fitted treating case/control status as the outcome and with pathogen-specific sample positivity as predictors to attain estimates of attributable pathogen fractions for AFI.

Discussion

The modular PCR platforms will allow for reporting of all primary results of respiratory samples within 72 hours and blood samples within one week, allowing for results to influence local medical practice and enable timely public health responses. The inclusion of controls will allow for a more accurate estimate of the importance of specific, prevalent pathogens as a cause of acute illness.

Study Registration

Project 1791, Registro de Proyectos de Investigación en Salud Pública (PRISA), Instituto Nacional de Salud, Perú.

Lineage-informative microhaplotypes for spatio-temporal surveillance of Plasmodium vivax malaria parasites

Challenges in understanding the origin of recurrent Plasmodium vivax infections constrains the surveillance of antimalarial efficacy and transmission of this neglected parasite. Recurrent infections within an individual may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or new inoculations (reinfection). Molecular inference of familial relatedness (identity-by-descent or IBD) based on whole genome sequence data, together with analysis of the intervals between parasitaemic episodes ("time-to-event" analysis), can help resolve the probable origin of recurrences. Whole genome sequencing of predominantly low-density P. vivax infections is challenging, so an accurate and scalable genotyping method to determine the origins of recurrent parasitaemia would be of significant benefit. We have developed a P. vivax genome-wide informatics pipeline to select specific microhaplotype panels that can capture IBD within small, amplifiable segments of the genome. Using a global set of 615 P. vivax genomes, we derived a panel of 100 microhaplotypes, each comprising 3-10 high frequency SNPs within <200 bp sequence windows. This panel exhibits high diversity in regions of the Asia-Pacific, Latin America and the horn of Africa (median HE = 0.70-0.81) and it captured 89% (273/307) of the polyclonal infections detected with genome-wide datasets. Using data simulations, we demonstrate lower error in estimating pairwise IBD using microhaplotypes, relative to traditional biallelic SNP barcodes. Our panel exhibited high accuracy in predicting the country of origin (median Matthew's correlation coefficient >0.9 in 90% countries tested) and it also captured local infection outbreak and bottlenecking events. The informatics pipeline is available open-source and yields microhaplotypes that can be readily transferred to high-throughput amplicon sequencing assays for surveillance in malaria-endemic regions.

Plasmodium vivax MSP1-42 kD Variant Proteins Detected Naturally Induced IgG Antibodies in Patients Regardless of the Infecting Parasite Phenotype in Mesoamerica

Background: The serological tests using blood stage antigens might be helpful for detecting recent exposure to Plasmodium parasites, and seroepidemiological studies would aid in the elimination of malaria. This work produced recombinant proteins of PvMSP142 variants and evaluated their capacity to detect IgG antibodies in symptomatic patients from Mesoamerica. Methods: Three variant Pvmsp142 genes were cloned in the pHL-sec plasmid, expressed in the Expi293F™ eukaryotic system, and the recombinant proteins were purified by affinity chromatography. Using an ELISA, 174 plasma or eluted samples from patients infected with different P. vivax haplotypes were evaluated against PvMSP142 proteins and to a native blood stage antigen (NBSA). Results: The antibody IgG OD values toward PvMSP142 variants (v88, v21, and v274) were heterogeneous (n = 178; median = 0.84 IQR 0.28–1.64). The correlation of IgG levels among all proteins was very high (spearman’s rho = 0.96–0.98; p < 0.0001), but was lower between them and the NBSA (rho = 0.771; p < 0.0001). In only a few samples, higher reactivity to the homologous protein was evident. Patients with a past infection who were seropositive had higher IgG levels and lower parasitemia levels than those who did not (p < 0.0001). Conclusions: The PvMSP142 variants were similarly efficient in detecting specific IgG antibodies in P. vivax patients from Mesoamerica, regardless of the infecting parasite’s haplotype, and might be good candidates for malaria surveillance and epidemiological studies in the region.

Retrospective analysis of Plasmodium vivax genomes from a pre-elimination China inland population in the 2010s

Introduction

In malaria-free countries, imported cases are challenging because interconnections with neighboring countries with higher transmission rates increase the risk of parasite reintroduction. Establishing a genetic database for rapidly identifying malaria importation or reintroduction is crucial in addressing these challenges. This study aimed to examine genomic epidemiology during the pre-elimination stage by retrospectively reporting whole-genome sequence variation of 10 Plasmodium vivax isolates from inland China.

Methods

The samples were collected during the last few inland outbreaks from 2011 to 2012 when China implemented a malaria control plan. After next-generation sequencing, we completed a genetic analysis of the population, explored the geographic specificity of the samples, and examined clustering of selection pressures. We also scanned genes for signals of positive selection.

Results

China's inland populations were highly structured compared to the surrounding area, with a single potential ancestor. Additionally, we identified genes under selection and evaluated the selection pressure on drug-resistance genes. In the inland population, positive selection was detected in some critical gene families, including sera, msp3, and vir. Meanwhile, we identified selection signatures in drug resistance, such as ugt, krs1, and crt, and noticed that the ratio of wild-type dhps and dhfr-ts increased after China banned sulfadoxine-pyrimethamine (SP) for decades.

Discussion

Our data provides an opportunity to investigate the molecular epidemiology of pre-elimination inland malaria populations, which exhibited lower selection pressure on invasion and immune evasion genes than neighbouring areas, but increased drug resistance in low transmission settings. Our results revealed that the inland population was severely fragmented with low relatedness among infections, despite a higher incidence of multiclonal infections, suggesting that superinfection or co-transmission events are rare in low-endemic circumstances. We identified selective signatures of resistance and found that the proportion of susceptible isolates fluctuated in response to the prohibition of specific drugs. This finding is consistent with the alterations in medication strategies during the malaria elimination campaign in inland China. Such findings could provide a genetic basis for future population studies, assessing changes in other pre-elimination countries.

Population-based genomic study of Plasmodium vivax malaria in seven Brazilian states and across South America

Background

Brazil is a unique and understudied setting for malaria, with complex foci of transmission associated with human and environmental conditions. An understanding of the population genomic diversity of P. vivax parasites across Brazil can support malaria control strategies.

Methods

Through whole genome sequencing of P. vivax isolates across 7 Brazilian states, we use population genomic approaches to compare genetic diversity within country (n = 123), continent (6 countries, n = 315) and globally (26 countries, n = 885).

Findings

We confirm that South American isolates are distinct, have more ancestral populations than the other global regions, with differentiating mutations in genes under selective pressure linked to antimalarial drugs (pvmdr1, pvdhfr-ts) and mosquito vectors (pvcrmp3, pvP45/48, pvP47). We demonstrate Brazil as a distinct parasite population, with signals of selection including ABC transporter (PvABCI3) and PHIST exported proteins.

Interpretation

Brazil has a complex population structure, with evidence of P. simium infections and Amazonian parasites separating into multiple clusters. Overall, our work provides the first Brazil-wide analysis of P. vivax population structure and identifies important mutations, which can inform future research and control measures.

Funding

AI is funded by an MRC LiD PhD studentship. TGC is funded by the Medical Research Council (Grant no. MR/M01360X/1, MR/N010469/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1). SC is funded by Medical Research Council UK grants (MR/M01360X/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1) and Bloomsbury SET (ref. CCF17-7779). FN is funded by The Shloklo Malaria Research Unit - part of the Mahidol Oxford Research Unit, supported by the Wellcome Trust (Grant no. 220211). ARSB is funded by São Paulo Research Foundation - FAPESP (Grant no. 2002/09546-1). RLDM is funded by Brazilian National Council for Scientific and Technological Development - CNPq (Grant no. 302353/2003-8 and 471605/2011-5); CRFM is funded by FAPESP (Grant no. 2020/06747-4) and CNPq (Grant no. 302917/2019-5 and 408636/2018-1); JGD is funded by FAPESP fellowships (2016/13465-0 and 2019/12068-5) and CNPq (Grant no. 409216/2018-6).

Genome-wide functional screening of drug-resistance genes in Plasmodium falciparum

The global spread of drug resistance is a major obstacle to the treatment of Plasmodium falciparum malaria. The identification of drug-resistance genes is an essential step toward solving the problem of drug resistance. Here, we report functional screening as a new approach with which to identify drug-resistance genes in P. falciparum. Specifically, a high-coverage genomic library of a drug-resistant strain is directly generated in a drug-sensitive strain, and the resistance gene is then identified from this library using drug screening. In a pilot experiment using the strain Dd2, the known chloroquine-resistant gene pfcrt is identified using the developed approach, which proves our experimental concept. Furthermore, we identify multidrug-resistant transporter 7 (pfmdr7) as a novel candidate for a mefloquine-resistance gene from a field-isolated parasite; we suggest that its upregulation possibly confers the mefloquine resistance. These results show the usefulness of functional screening as means by which to identify drug-resistance genes.

The Impact, Emerging Needs, and New Research Questions Arising from 12 Years of the Center for the Study of Complex Malaria in India

The Center for the Study of Complex Malaria in India (CSCMi) was launched in 2010 with the overall goal of addressing major gaps in our understanding of "complex malaria" in India through projects on the epidemiology, transmission, and pathogenesis of the disease. The Center was mandated to adopt an integrated approach to malaria research, including building capacity, developing infrastructure, and nurturing future malaria leaders while conducting relevant and impactful studies to assist India as it moves from control to elimination. Here, we will outline some of the interactions and impacts the Center has had with malaria policy and control counterparts in India, as well as describe emerging needs and new research questions that have become apparent over the past 12 years.

Malaria Resilience in South America: Epidemiology, Vector Biology, and Immunology Insights from the Amazonian International Center of Excellence in Malaria Research Network in Peru and Brazil

The 1990s saw the rapid reemergence of malaria in Amazonia, where it remains an important public health priority in South America. The Amazonian International Center of Excellence in Malaria Research (ICEMR) was designed to take a multidisciplinary approach toward identifying novel malaria control and elimination strategies. Based on geographically and epidemiologically distinct sites in the Northeastern Peruvian and Western Brazilian Amazon regions, synergistic projects integrate malaria epidemiology, vector biology, and immunology. The Amazonian ICEMR's overarching goal is to understand how human behavior and other sociodemographic features of human reservoirs of transmission-predominantly asymptomatically parasitemic people-interact with the major Amazonian malaria vector, Nyssorhynchus (formerly Anopheles) darlingi, and with human immune responses to maintain malaria resilience and continued endemicity in a hypoendemic setting. Here, we will review Amazonian ICEMR's achievements on the synergies among malaria epidemiology, Plasmodium-vector interactions, and immune response, and how those provide a roadmap for further research, and, most importantly, point toward how to achieve malaria control and elimination in the Americas.

Drug resistance and population structure of Plasmodium falciparum and Plasmodium vivax in the Peruvian Amazon

Malaria is a major health problem in Peru despite substantial progress achieved by the ongoing malaria elimination program. This study explored the population genetics of 63 Plasmodium falciparum and 170 P. vivax cases collected in the Peruvian Amazon Basin between 2015 and 2019. Microscopy and PCR were used for malaria detection and positive samples were genotyped at neutral and drug resistance-associated regions. The P. falciparum population exhibited a low nucleotide diversity (π = 0.02) whereas the P. vivax population presented a higher genetic diversity (π = 0.34). All P. falciparum samples (n = 63) carried chloroquine (CQ) resistant mutations on Pfcrt. Most P. falciparum samples (53 out of 54) carried sulfadoxine (SD) resistant mutations on Pfdhfr and Pfdhps. No evidence was found of artemisinin resistance mutations on kelch13. Population structure showed that a single cluster accounted for 93.4% of the P. falciparum samples whereas three clusters were found for P. vivax. Our study shows a low genetic diversity for both species with significant differences in genetic sub-structuring. The high prevalence of CQ-resistance mutations could be a result of indirect selection pressures driven by the P. vivax treatment scheme. These results could be useful for public health authorities to safeguard the progress that Peru has achieved towards malaria elimination.

Prospects for genomic surveillance for selection in schistosome parasites

Schistosomiasis is a neglected tropical disease caused by multiple parasitic Schistosoma species, and which impacts over 200 million people globally, mainly in low- and middle-income countries. Genomic surveillance to detect evidence for natural selection in schistosome populations represents an emerging and promising approach to identify and interpret schistosome responses to ongoing control efforts or other environmental factors. Here we review how genomic variation is used to detect selection, how these approaches have been applied to schistosomes, and how future studies to detect selection may be improved. We discuss the theory of genomic analyses to detect selection, identify experimental designs for such analyses, and review studies that have applied these approaches to schistosomes. We then consider the biological characteristics of schistosomes that are expected to respond to selection, particularly those that may be impacted by control programs. Examples include drug resistance, host specificity, and life history traits, and we review our current understanding of specific genes that underlie them in schistosomes. We also discuss how inherent features of schistosome reproduction and demography pose substantial challenges for effective identification of these traits and their genomic bases. We conclude by discussing how genomic surveillance for selection should be designed to improve understanding of schistosome biology, and how the parasite changes in response to selection.

Population genomics in neglected malaria parasites

Malaria elimination includes neglected human malaria parasites Plasmodium vivax, Plasmodium ovale spp., and Plasmodium malariae. Biological features such as association with low-density infection and the formation of hypnozoites responsible for relapse make their elimination challenging. Studies on these parasites rely primarily on clinical samples due to the lack of long-term culture techniques. With improved methods to enrich parasite DNA from clinical samples, whole-genome sequencing of the neglected malaria parasites has gained increasing popularity. Population genomics of more than 2200 P. vivax global isolates has improved our knowledge of parasite biology and host-parasite interactions, identified vaccine targets and potential drug resistance markers, and provided a new way to track parasite migration and introduction and monitor the evolutionary response of local populations to elimination efforts. Here, we review advances in population genomics for neglected malaria parasites, discuss how the rich genomic information is being used to understand parasite biology and epidemiology, and explore opportunities for the applications of malaria genomic data in malaria elimination practice.

Novel highly-multiplexed AmpliSeq targeted assay for Plasmodium vivax genetic surveillance use cases at multiple geographical scales

Although the power of genetic surveillance tools has been acknowledged widely, there is an urgent need in malaria endemic countries for feasible and cost-effective tools to implement in national malaria control programs (NMCPs) that can generate evidence to guide malaria control and elimination strategies, especially in the case of Plasmodium vivax. Several genetic surveillance applications ('use cases') have been identified to align research, technology development, and public health efforts, requiring different types of molecular markers. Here we present a new highly-multiplexed deep sequencing assay (Pv AmpliSeq). The assay targets the 33-SNP vivaxGEN-geo panel for country-level classification, and a newly designed 42-SNP within-country barcode for analysis of parasite dynamics in Vietnam and 11 putative drug resistance genes in a highly multiplexed NGS protocol with easy workflow, applicable for many different genetic surveillance use cases. The Pv AmpliSeq assay was validated using: 1) isolates from travelers and migrants in Belgium, and 2) routine collections of the national malaria control program at sentinel sites in Vietnam. The assay targets 229 amplicons and achieved a high depth of coverage (mean 595.7 ± 481) and high accuracy (mean error-rate of 0.013 ± 0.007). P. vivax parasites could be characterized from dried blood spots with a minimum of 5 parasites/µL and 10% of minority-clones. The assay achieved good spatial specificity for between-country prediction of origin using the 33-SNP vivaxGEN-geo panel that targets rare alleles specific for certain countries and regions. A high resolution for within-country diversity in Vietnam was achieved using the designed 42-SNP within-country barcode that targets common alleles (median MAF 0.34, range 0.01-0.49. Many variants were detected in (putative) drug resistance genes, with different predominant haplotypes in the pvmdr1 and pvcrt genes in different provinces in Vietnam. The capacity of the assay for high resolution identity-by-descent (IBD) analysis was demonstrated and identified a high rate of shared ancestry within Gia Lai Province in the Central Highlands of Vietnam, as well as between the coastal province of Binh Thuan and Lam Dong. Our approach performed well in geographically differentiating isolates at multiple spatial scales, detecting variants in putative resistance genes, and can be easily adjusted to suit the needs in other settings in a country or region. We prioritize making this tool available to researchers and NMCPs in endemic countries to increase ownership and ensure data usage for decision-making and malaria policy.

Complete Genomes of Theileria orientalis Chitose and Buffeli Genotypes Reveal within Species Translocations and Differences in ABC Transporter Content

Theileria orientalis causes losses to cattle producers in Eastern Asia, Oceania and, more recently, North America. One pathogenic genotype (Ikeda) has been sequenced to the chromosomal level, while only draft genomes exist for globally distributed Chitose and Buffeli genotypes. To provide an accurate comparative gene-level analysis and help further understand their pathogenicity, we sequenced isolates of the Chitose and Buffeli genotypes of T. orientalis using long-read sequencing technology. A combination of several long-read assembly methods and short reads produced chromosomal-level assemblies for both Fish Creek (Chitose) and Goon Nure (Buffeli) isolates, including the first complete and circular apicoplast genomes generated for T. orientalis. Comparison with the Shintoku (Ikeda) reference sequence showed both large and small translocations in T. orientalis Buffeli, between chromosomes 2 and 3 and chromosomes 1 and 4, respectively. Ortholog clustering showed expansion of ABC transporter genes in Chitose and Buffeli. However, differences in several genes of unknown function, including DUF529/FAINT-domain-containing proteins, were also identified and these genes were more prevalent in Ikeda and Chitose genotypes. Phylogenetics and similarity measures were consistent with previous short-read genomic analysis. The generation of chromosomal sequences for these highly prevalent T. orientalis genotypes will also support future studies of population genetics and mixed genotype infections.

Genetic diversity of Plasmodium vivax reticulocyte binding protein 2b in global parasite populations

Background

Plasmodium vivax reticulocyte binding protein 2b (PvRBP2b) plays a critical role in parasite invasion of reticulocytes by binding the transferrin receptor 1. PvRBP2b is a vaccine candidate based on the negative correlation between antibody titers against PvRBP2b recombinant proteins and parasitemia and risk of vivax malaria. The aim of this study was to analyze the genetic diversity of the PvRBP2b gene in the global P. vivax populations.

Methods

Near full-length PvRBP2b nucleotide sequences (190–8349 bp) were obtained from 88 P. vivax isolates collected from the China–Myanmar border (n = 44) and Thailand (n = 44). An additional 224 PvRBP2b sequences were retrieved from genome sequences from parasite populations worldwide. The genetic diversity, neutral selection, haplotype distribution and genetic differentiation of PvRBP2b were examined.

Results

The genetic diversity of PvRBP2b was distributed unevenly, with peak diversity found in the reticulocyte binding region in the N-terminus. Neutrality analysis suggested that this region is subjected to balancing selection or population bottlenecks. Several amino acid variants were found in all or nearly all P. vivax endemic regions. However, the critical residues responsible for reticulocyte binding were highly conserved. There was substantial population differentiation according to the geographical separation. The distribution of haplotypes in the reticulocyte binding region varied among regions; even the two major haplotypes Hap_6 and Hap_8 were found in only five populations.

Conclusions

Our data show considerable genetic variations of PvRBPb in global parasite populations. The geographic divergence may pose a challenge to PvRBP2b-based vaccine development.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05296-6.

Evolutionary history of Plasmodium vivax and Plasmodium simium in the Americas

Malaria is a vector-borne disease caused by protozoan parasites of the genus Plasmodium. Plasmodium vivax is the most prevalent human-infecting species in the Americas. However, the origins of this parasite in this continent are still debated. Similarly, it is now accepted that the existence of Plasmodium simium is explained by a P. vivax transfer from humans to monkey in America. However, many uncertainties still exist concerning the origin of the transfer and whether several transfers occurred. In this review, the most recent studies that addressed these questions using genetic and genomic approaches are presented.

An open dataset of Plasmodium vivax genome variation in 1,895 worldwide samples

This report describes the MalariaGEN Pv4 dataset, a new release of curated genome variation data on 1,895 samples of Plasmodium vivax collected at 88 worldwide locations between 2001 and 2017. It includes 1,370 new samples contributed by MalariaGEN and VivaxGEN partner studies in addition to previously published samples from these and other sources. We provide genotype calls at over 4.5 million variable positions including over 3 million single nucleotide polymorphisms (SNPs), as well as short indels and tandem duplications. This enlarged dataset highlights major compartments of parasite population structure, with clear differentiation between Africa, Latin America, Oceania, Western Asia and different parts of Southeast Asia. Each sample has been classified for drug resistance to sulfadoxine, pyrimethamine and mefloquine based on known markers at the dhfr, dhps and mdr1 loci. The prevalence of all of these resistance markers was much higher in Southeast Asia and Oceania than elsewhere. This open resource of analysis-ready genome variation data from the MalariaGEN and VivaxGEN networks is driven by our collective goal to advance research into the complex biology of P. vivax and to accelerate genomic surveillance for malaria control and elimination.

Whole Genome Sequencing Contributions and Challenges in Disease Reduction Focused on Malaria

Malaria elimination remains an important goal that requires the adoption of sophisticated science and management strategies in the era of the COVID-19 pandemic. The advent of next generation sequencing (NGS) is making whole genome sequencing (WGS) a standard today in the field of life sciences, as PCR genotyping and targeted sequencing provide insufficient information compared to the whole genome. Thus, adapting WGS approaches to malaria parasites is pertinent to studying the epidemiology of the disease, as different regions are at different phases in their malaria elimination agenda. Therefore, this review highlights the applications of WGS in disease management, challenges of WGS in controlling malaria parasites, and in furtherance, provides the roles of WGS in pursuit of malaria reduction and elimination. WGS has invaluable impacts in malaria research and has helped countries to reach elimination phase rapidly by providing required information needed to thwart transmission, pathology, and drug resistance. However, to eliminate malaria in sub-Saharan Africa (SSA), with high malaria transmission, we recommend that WGS machines should be readily available and affordable in the region.

Primate malarias as a model for cross-species parasite transmission

Parasites regularly switch into new host species, representing a disease burden and conservation risk to the hosts. The distribution of these parasites also gives insight into characteristics of ecological networks and genetic mechanisms of host-parasite interactions. Some parasites are shared across many species, whereas others tend to be restricted to hosts from a single species. Understanding the mechanisms producing this distribution of host specificity can enable more effective interventions and potentially identify genetic targets for vaccines or therapies. As ecological connections between human and local animal populations increase, the risk to human and wildlife health from novel parasites also increases. Which of these parasites will fizzle out and which have the potential to become widespread in humans? We consider the case of primate malarias, caused by Plasmodium parasites, to investigate the interacting ecological and evolutionary mechanisms that put human and nonhuman primates at risk for infection. Plasmodium host switching from nonhuman primates to humans led to ancient introductions of the most common malaria-causing agents in humans today, and new parasite switching is a growing threat, especially in Asia and South America. Based on a wild host-Plasmodium occurrence database, we highlight geographic areas of concern and potential areas to target further sampling. We also discuss methodological developments that will facilitate clinical and field-based interventions to improve human and wildlife health based on this eco-evolutionary perspective.

Spatiotemporal Changes in Plasmodium vivax msp142 Haplotypes in Southern Mexico: From the Control to the Pre-Elimination Phase

For 20 years, Plasmodium vivax has been the only prevalent malaria species in Mexico, and cases have declined significantly and continuously. Spatiotemporal genetic studies can be helpful for understanding parasite dynamics and developing strategies to weaken malaria transmission, thus facilitating the elimination of the parasite. The aim of the current contribution was to analyze P. vivax-infected blood samples from patients in southern Mexico during the control (1993–2007) and pre-elimination phases (2008–2011). Nucleotide and haplotype changes in the pvmsp142 fragment were evaluated over time. The majority of multiple genotype infections occurred in the 1990s, when the 198 single nucleotide sequences exhibited 57 segregating sites, 64 mutations, and 17 haplotypes. Nucleotide and genetic diversity parameters showed subtle fluctuations from across time, in contrast to the reduced haplotype diversity and the increase in the R2 index and Tajima’s D value from 2008 to 2011. The haplotype network consisted of four haplogroups, the geographical distribution of which varied slightly over time. Haplogroup-specific B-cell epitopes were predicted. Since only high-frequency and divergent haplotypes persisted, there was a contraction of the parasite population. Given that 84% of haplotypes were exclusive to Mesoamerica, P. vivax flow is likely circumscribed to this region, representing important information for parasite surveillance.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

Integrating genomic and epidemiologic data to accelerate progress toward schistosomiasis elimination

The global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

Plasmodium vivax From Duffy-Negative and Duffy-Positive Individuals Share Similar Gene Pools in East Africa

Plasmodium vivax malaria was thought to be rare in Africa, but an increasing number of P. vivax cases reported across Africa and in Duffy-negative individuals challenges this dogma. The genetic characteristics of P. vivax in Duffy-negative infections, the transmission of P. vivax in East Africa, and the impact of environments on transmission remain largely unknown. This study examined genetic and transmission features of P. vivax from 107 Duffy-negative and 305 Duffy-positive individuals in Ethiopia and Sudan. No clear genetic differentiation was found in P. vivax between the 2 Duffy groups, indicating between-host transmission. P. vivax from Ethiopia and Sudan showed similar genetic clusters, except samples from Khartoum, possibly due to distance and road density that inhibited parasite gene flow. This study is the first to show that P. vivax can transmit to and from Duffy-negative individuals and provides critical insights into the spread of P. vivax in sub-Saharan Africa.

The genome of the zoonotic malaria parasite Plasmodium simium reveals adaptations to host switching

BACKGROUND

Plasmodium simium, a malaria parasite of non-human primates (NHP), was recently shown to cause zoonotic infections in humans in Brazil. We sequenced the P. simium genome to investigate its evolutionary history and to identify any genetic adaptions that may underlie the ability of this parasite to switch between host species.

RESULTS

Phylogenetic analyses based on whole genome sequences of P. simium from humans and NHPs reveals that P. simium is monophyletic within the broader diversity of South American Plasmodium vivax, suggesting P. simium first infected NHPs as a result of a host switch of P. vivax from humans. The P. simium isolates show the closest relationship to Mexican P. vivax isolates. Analysis of erythrocyte invasion genes reveals differences between P. vivax and P. simium, including large deletions in the Duffy-binding protein 1 (DBP1) and reticulocyte-binding protein 2a genes of P. simium. Analysis of P. simium isolated from NHPs and humans revealed a deletion of 38 amino acids in DBP1 present in all human-derived isolates, whereas NHP isolates were multi-allelic.

CONCLUSIONS

Analysis of the P. simium genome confirmed a close phylogenetic relationship between P. simium and P. vivax, and suggests a very recent American origin for P. simium. The presence of the DBP1 deletion in all human-derived isolates tested suggests that this deletion, in combination with other genetic changes in P. simium, may facilitate the invasion of human red blood cells and may explain, at least in part, the basis of the recent zoonotic infections.

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.

A population genetic perspective on the origin, spread and adaptation of the human malaria agents Plasmodium falciparum and Plasmodium vivax

Malaria is considered one of the most important scourges that humanity has faced during its history, being responsible every year for numerous deaths worldwide. The disease is caused by protozoan parasites, among which two species are responsible of the majority of the burden, Plasmodium falciparum and Plasmodium vivax. For these two parasite species, the questions of their origin (how and when they appeared in humans), of their spread throughout the world, as well as how they have adapted to humans have long been of interest to the scientific community. Here, we review the current knowledge that has accumulated on these different questions, thanks in particular to the analysis of the genetic and genomic variability of these parasites and comparison with related Plasmodium species infecting other host species (like non-human primates). In this paper we review the existing body of knowledge, including current research dealing with these questions, focusing particularly on genetic analysis and genomic variability of these parasites and comparison with related Plasmodium species infecting other species of host (such as non-human primates).

Genetic diversity in two leading Plasmodium vivax malaria vaccine candidates AMA1 and MSP119 at three sites in India

Plasmodium vivax, a major contributor to the malaria burden in India, has the broadest geographic distribution and shows higher genetic diversity than P. falciparum. Here, we investigated the genetic diversity of two leading P. vivax vaccine candidate antigens, at three geographically diverse malaria-endemic regions in India. Pvama1 and Pvmsp119 partial coding sequences were generated from one hundred P. vivax isolates in India (Chennai n = 28, Nadiad n = 50 and Rourkela n = 22) and ~1100 published sequences from Asia, South America, North America, and Oceania regions included. These data were used to assess the genetic diversity and potential for vaccine candidacy of both antigens on a global scale. A total of 44 single nucleotide polymorphism (SNPs) were identified among 100 Indian Pvama1 sequences, including 10 synonymous and 34 nonsynonymous mutations. Nucleotide diversity was higher in Rourkela and Nadiad as compared to Chennai. Nucleotide diversity measures showed a strong balancing selection in Indian and global population for domain I of Pvama1, which suggests that it is a dominant target of the protective immune response. In contrast, the Pvmsp119 region showed highly conserved sequences in India and across the Oceania, South America, North America and Asia, demonstrating low genetic diversity in the global population when compared to Pvama1. Results suggest the possibility of including Pvmsp119 in a multivalent vaccine formulation against P. vivax infections. However, the high genetic diversity seen in Pvama1 would be more challenging for vaccine development.

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.

The Potential Role of Clinical Metagenomics in Infectious Diseases: Therapeutic Perspectives

Clinical metagenomics (CMg) is the process of sequencing nucleic acid of clinical samples to obtain clinically relevant information such as the identification of microorganisms and their susceptibility to antimicrobials. Over the last decades, sequencing and bioinformatic solutions supporting CMg have much evolved and an increasing number of case reports and series covering various infectious diseases have been published. Metagenomics is a new approach to infectious disease diagnosis that is currently being developed and is certainly one of the most promising for the coming years. However, most CMg studies are retrospective, and few address the potential impact CMg could have on patient management, including initiation, adaptation, or cessation of antimicrobials. In this narrative review, we have discussed the potential role of CMg in bacteriology, virology, mycology, and parasitology. Several reports and case-series confirm that CMg is an innovative tool with which one can (i) identify more microorganisms than with conventional methods in a single test, (ii) obtain results within hours, and (iii) tailor the antimicrobial regimen of patients. However, the cost-efficiency of CMg and its real impact on patient management are still to be determined.

Ape Origins of Human Malaria

African apes harbor at least twelve Plasmodium species, some of which have been a source of human infection. It is now well established that Plasmodium falciparum emerged following the transmission of a gorilla parasite, perhaps within the last 10,000 years, while Plasmodium vivax emerged earlier from a parasite lineage that infected humans and apes in Africa before the Duffy-negative mutation eliminated the parasite from humans there. Compared to their ape relatives, both human parasites have greatly reduced genetic diversity and an excess of nonsynonymous mutations, consistent with severe genetic bottlenecks followed by rapid population expansion. A putative new Plasmodium species widespread in chimpanzees, gorillas, and bonobos places the origin of Plasmodium malariae in Africa. Here, we review what is known about the origins and evolutionary history of all human-infective Plasmodium species, the time and circumstances of their emergence, and the diversity, host specificity, and zoonotic potential of their ape counterparts.

PvP01-DB: computational structural and functional characterization of soluble proteome of PvP01 strain of Plasmodium vivax

Despite Plasmodium vivax being the main offender in the majority of malarial infections, very little information is available about its adaptation and development in humans. Its capability for activating relapsing infections through its dormant liver stage and resistance to antimalarial drugs makes it as one of the major challenges in eradicating malaria. Noting the immediate necessity for the availability of a comprehensive and reliable structural and functional repository for P. vivax proteome, here we developed a web resource for the new reference genome, PvP01, furnishing information on sequence, structure, functions, active sites and metabolic pathways compiled and predicted using some of the state-of-the-art methods in respective fields. The PvP01 web resource comprises organized data on the soluble proteome consisting of 3664 proteins in blood and liver stages of malarial cycle. The current public resources represent only 163 proteins of soluble proteome of PvP01, with complete information about their molecular function, biological process and cellular components. Also, only 46 proteins of P. vivax have experimentally determined structures. In this milieu of extreme scarcity of structural and functional information, PvP01 web resource offers meticulously validated structures of 3664 soluble proteins. The sequence and structure-based functional characterization led to a quantum leap from 163 proteins available presently to whole soluble proteome offered through PvP01 web resource. We believe PvP01 web resource will serve the researchers in identifying novel protein drug targets and in accelerating the development of structure-based new drug candidates to combat malaria.

Database Availability: http://www.scfbio-iitd.res.in/PvP01