Whole-genome sequencing and microarray analysis of ex vivo Plasmodium vivax reveal selective pressure on putative drug resistance genes

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.

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.

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.

No Association between the Plasmodium vivax crt-o MS334 or In9pvcrt Polymorphisms and Chloroquine Failure in a Pre-Elimination Clinical Cohort from Malaysia with a Large Clonal Expansion

Increasing reports of resistance to a frontline malaria blood-stage treatment, chloroquine (CQ), raises concerns for the elimination of Plasmodium vivax. The absence of an effective molecular marker of CQ resistance in P. vivax greatly constrains surveillance of this emerging threat. A recent genetic cross between CQ sensitive (CQS) and CQ resistant (CQR) NIH-1993 strains of P. vivax linked a moderate CQR phenotype with two candidate markers in P. vivax CQ resistance transporter gene (pvcrt-o): MS334 and In9pvcrt. Longer TGAAGH motif lengths at MS334 were associated with CQ resistance, as were shorter motifs at the In9pvcrt locus. In this study, high-grade CQR clinical isolates of P. vivax from a low endemic setting in Malaysia were used to investigate the association between the MS334 and In9pvcrt variants and treatment efficacy. Among a total of 49 independent monoclonal P. vivax isolates assessed, high-quality MS334 and In9pvcrt sequences could be derived from 30 (61%) and 23 (47%), respectively. Five MS334 and six In9pvcrt alleles were observed, with allele frequencies ranging from 2 to 76% and 3 to 71%, respectively. None of the clinical isolates had the same variant as the NIH-1993 CQR strain, and none of the variants were associated with CQ treatment failure (all P > 0.05). Multi-locus genotypes (MLGs) at 9 neutral microsatellites revealed a predominant P. vivax strain (MLG6) accounting for 52% of Day 0 infections. The MLG6 strain comprised equal proportions of CQS and CQR infections. Our study reveals complexity in the genetic basis of CQ resistance in the Malaysian P. vivax pre-elimination setting and suggests that the proposed pvcrt-o MS334 and In9pvcrt markers are not reliable markers of CQ treatment efficacy in this setting. Further studies are needed in other endemic settings, applying hypothesis-free genome-wide approaches, and functional approaches to understand the biological impact of the TGAAGH repeats linked to CQ response in a cross are warranted to comprehend and track CQR P. vivax.

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.

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.

Decrypting the complexity of the human malaria parasite biology through systems biology approaches

The human malaria parasite, Plasmodium falciparum, is a unicellular protozoan responsible for over half a million deaths annually. With a complex life cycle alternating between human and invertebrate hosts, this apicomplexan is notoriously adept at evading host immune responses and developing resistance to all clinically administered treatments. Advances in omics-based technologies, increased sensitivity of sequencing platforms and enhanced CRISPR based gene editing tools, have given researchers access to more in-depth and untapped information about this enigmatic micro-organism, a feat thought to be infeasible in the past decade. Here we discuss some of the most important scientific achievements made over the past few years with a focus on novel technologies and platforms that set the stage for subsequent discoveries. We also describe some of the systems-based methods applied to uncover gaps of knowledge left through single-omics applications with the hope that we will soon be able to overcome the spread of this life-threatening disease.

Single-genome sequencing reveals within-host evolution of human malaria parasites

Population genomics of bulk malaria infections is unable to examine intrahost evolution; therefore, most work has focused on the role of recombination in generating genetic variation. We used single-cell sequencing protocol for low-parasitaemia infections to generate 406 near-complete single Plasmodium vivax genomes from 11 patients sampled during sequential febrile episodes. Parasite genomes contain hundreds of de novo mutations, showing strong signatures of selection, which are enriched in the ApiAP2 family of transcription factors, known targets of adaptation. Comparing 315 P. falciparum single-cell genomes from 15 patients with our P. vivax data, we find broad complementary patterns of de novo mutation at the gene and pathway level, revealing the importance of within-host evolution during malaria infections.

Ten-Year Molecular Surveillance of Drug-Resistant Plasmodium spp. Isolated From the China-Myanmar Border

Antimalarial drug resistance has emerged as a major threat to global malaria control efforts, particularly in the Greater Mekong Subregion (GMS). In this study, we analyzed the polymorphism and prevalence of molecular markers associated with resistance to first-line antimalarial drugs, such as artemisinin, chloroquine, and pyrimethamine, using blood samples collected from malaria patients in the China-Myanmar border region of the GMS from 2008 to 2017, including 225 cases of Plasmodium falciparum and 194 cases of Plasmodium vivax. In artemisinin resistance, only the C580Y mutation with low frequency was detected in pfk13, and no highly frequent stable mutation was found in pvk12. In chloroquine resistance, the frequency of K76T mutation in pfcrt was always high, and the frequency of double mutations in pvmdr1 of P. vivax has been steadily increasing every year. In pyrimidine resistance, pfdhfr and pvdhfr had relatively more complex mutant types associated with drug resistance sites, and the overall mutation rate was still high. Therefore, artemisinin-based combination therapies are still suitable for use as the first choice of antimalarial strategy in the China-Myanmar border region in the future.

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

The molecular basis of antimalarial drug resistance in Plasmodium vivax

Plasmodium vivax is the most geographically widespread cause of human malaria and is responsible for the majority of cases outside of the African continent. While great progress has been made towards eliminating human malaria, drug resistant parasite strains pose a threat towards continued progress. Resistance has arisen to multiple antimalarials in P. vivax, including to chloroquine, which is currently the first line therapy for P. vivax in most regions. Despite its importance, an understanding of the molecular mechanisms of drug resistance in this species remains elusive, in large part due to the complex biology of P. vivax and the lack of in vitro culture. In this review, we will cover the extent and challenges of measuring clinical and in vitro drug resistance in P. vivax. We will consider the roles of candidate drug resistance genes. We will highlight the development of molecular approaches for studying P. vivax biology that provide the opportunity to validate the role of putative drug resistance mutations as well as identify novel mechanisms of drug resistance in this understudied parasite. Validated molecular determinants and markers of drug resistance are essential for the rapid and cost-effective monitoring of drug resistance in P. vivax, and will be useful for optimizing drug regimens and for informing drug policy in control and elimination settings.

•

Drug resistance is emerging in Plasmodium vivax, an important cause of malaria.

•

The complex biology of P. vivax and the limited range of research tools make it difficult to identify drug resistance.

•

The molecular mechanisms of drug resistance in P. vivax remain elusive.

•

This review highlights the extent of drug resistance, the putative mechanisms of resistance and new technologies for the study of P. vivax drug resistance.

Plasmodium vivax infection compromises reticulocyte stability

The structural integrity of the host red blood cell (RBC) is crucial for propagation of Plasmodium spp. during the disease-causing blood stage of malaria infection. To assess the stability of Plasmodium vivax-infected reticulocytes, we developed a flow cytometry-based assay to measure osmotic stability within characteristically heterogeneous reticulocyte and P. vivax-infected samples. We find that erythroid osmotic stability decreases during erythropoiesis and reticulocyte maturation. Of enucleated RBCs, young reticulocytes which are preferentially infected by P. vivax, are the most osmotically stable. P. vivax infection however decreases reticulocyte stability to levels close to those of RBC disorders that cause hemolytic anemia, and to a significantly greater degree than P. falciparum destabilizes normocytes. Finally, we find that P. vivax new permeability pathways contribute to the decreased osmotic stability of infected-reticulocytes. These results reveal a vulnerability of P. vivax-infected reticulocytes that could be manipulated to allow in vitro culture and develop novel therapeutics.

During Plasmodium intra-erythrocytic developmental, parasites compromise the structural integrity of host red-blood cells. Here, Clark et al. develop a flow cytometric osmotic stability assay to show that P. vivax infection destabilizes host reticulocytes, which are less stable than P. falciparum-infected normocytes.

Whole genome sequencing of Plasmodium vivax isolates reveals frequent sequence and structural polymorphisms in erythrocyte binding genes

Plasmodium vivax malaria is much less common in Africa than the rest of the world because the parasite relies primarily on the Duffy antigen/chemokine receptor (DARC) to invade human erythrocytes, and the majority of Africans are Duffy negative. Recently, there has been a dramatic increase in the reporting of P. vivax cases in Africa, with a high number of them being in Duffy negative individuals, potentially indicating P. vivax has evolved an alternative invasion mechanism that can overcome Duffy negativity. Here, we analyzed single nucleotide polymorphism (SNP) and copy number variation (CNV) in Whole Genome Sequence (WGS) data from 44 P. vivax samples isolated from symptomatic malaria patients in southwestern Ethiopia, where both Duffy positive and Duffy negative individuals are found. A total of 123,711 SNPs were detected, of which 22.7% were nonsynonymous and 77.3% were synonymous mutations. The largest number of SNPs were detected on chromosomes 9 (24,007 SNPs; 19.4% of total) and 10 (16,852 SNPs, 13.6% of total). There were particularly high levels of polymorphism in erythrocyte binding gene candidates including merozoite surface protein 1 (MSP1) and merozoite surface protein 3 (MSP3.5, MSP3.85 and MSP3.9). Two genes, MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. Variation in gene copy number was also concentrated in genes involved in host-parasite interactions, including the expansion of the Duffy binding protein gene (PvDBP) on chromosome 6 and MSP3.11 on chromosome 10. Based on the phylogeny constructed from the whole genome sequences, the expansion of these genes was an independent process among the P. vivax lineages in Ethiopia. We further inferred transmission patterns of P. vivax infections among study sites and showed various levels of gene flow at a small geographical scale. The genomic features of P. vivax provided baseline data for future comparison with those in Duffy-negative individuals and allowed us to develop a panel of informative Single Nucleotide Polymorphic markers diagnostic at a micro-geographical scale.

Plasmodium vivax is the most geographically widespread parasite species that causes malaria in humans. Although it occurs in Africa as a member of a mix of Plasmodium species, P. vivax is dominant in other parts of the world outside of Africa (e.g., Brazil). It was previously thought that most African populations were immune to P. vivax infections due to the absence of Duffy antigen chemokine receptor (DARC) gene expression required for erythrocyte invasion. However, several recent reports have indicated the emergence and potential spread of P. vivax across human populations in Africa. Compared to Southeast Asia and South America where P. vivax is highly endemic, data on polymorphisms in erythrocyte binding gene candidates of P. vivax from Africa is limited. Filling this knowlege gap is critical for identifying functional genes in erythrocyte invasion, biomarkers for tracking the P. vivax isolates from Africa, as well as potential gene targets for vaccine development. This paper examined the level of genetic polymorphisms in a panel of 43 potential erythrocyte binding protein genes based on whole genome sequences and described transmission patterns of P. vivax infections from different study sites in Ethiopia based on the genetic variants. Our analyses showed that chromosomes 9 and 10 of the P. vivax genomes isolated in Ethiopia had the most high-quality genetic polymorphisms. Among all erythrocyte binding protein gene candidates, the merozoite surface proteins 1 and merozoite surface protein 3 showed high levels of polymorphism. MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. The expansion of the Duffy binding protein and merozoite surface protein 3 gene copies was an independent process among the P. vivax lineages in Ethiopia. Various levels of gene flow were observed even at a smaller geographical scale. Our study provided baseline data for future comparison with P. vivax in Duffy negative individuals and help develop a panel of genetic markers that are informative at a micro-geographical scale.

Putative circumsporozoite protein (CSP) of Plasmodium vivax is considerably distinct from the well-known CSP and plays a role in the protein ubiquitination pathway

Amidst technical challenges which limit successful culture and genetic manipulation of P. vivax parasites, we used a computational approach to identify a critical target with evolutionary significance. The putative circumsporozoite protein on chromosome 13 of P. vivax (PvpuCSP)is distinct from the well-known vaccine candidate PfCSP. The aim of this study was to understand the role of PvpuCSP and its relatedness to the well-known CSP. The study revealed PvpuCSP as a membrane bound E3 ubiquitin ligase involved in ubiquitination. It has a species-specific tetra-peptide unit which is differentially repeated in various P. vivax strains. The PvpuCSP is different from CSP in terms of stage-specific expression and function. Since E3 ubiquitin ligases are known antimalarial drug targets targeting the proteasome pathway, PvpuCSP, with evolutionary connotation and a key role in orchestrating protein degradation in P. vivax, can be explored as a potential drug target.

•

PvpuCSP is predicted as E3 ubiquitin ligase, a part of ubiquitination pathway.

•

Tetra-peptide tandem repeat at C terminal of PvpuCSP is exclusive to P. vivax.

•

Moderately expressed during all parasitic stages in host and vector

•

Partially disordered protein with both structured domains and two distinct IDRs

•

A transmembrane protein with highly conserved functional domain across Apicomplexa

Plasmodium vivax Malaria Viewed through the Lens of an Eradicated European Strain

The protozoan Plasmodium vivax is responsible for 42% of all cases of malaria outside Africa. The parasite is currently largely restricted to tropical and subtropical latitudes in Asia, Oceania, and the Americas. Though, it was historically present in most of Europe before being finally eradicated during the second half of the 20th century. The lack of genomic information on the extinct European lineage has prevented a clear understanding of historical population structuring and past migrations of P. vivax. We used medical microscope slides prepared in 1944 from malaria-affected patients from the Ebro Delta in Spain, one of the last footholds of malaria in Europe, to generate a genome of a European P. vivax strain. Population genetics and phylogenetic analyses placed this strain basal to a cluster including samples from the Americas. This genome allowed us to calibrate a genomic mutation rate for P. vivax, and to estimate the mean age of the last common ancestor between European and American strains to the 15th century. This date points to an introduction of the parasite during the European colonization of the Americas. In addition, we found that some known variants for resistance to antimalarial drugs, including Chloroquine and Sulfadoxine, were already present in this European strain, predating their use. Our results shed light on the evolution of an important human pathogen and illustrate the value of antique medical collections as a resource for retrieving genomic information on pathogens from the past.

Plasmodium vivax chloroquine resistance links to pvcrt transcription in a genetic cross

Mainstay treatment for Plasmodium vivax malaria has long relied on chloroquine (CQ) against blood-stage parasites plus primaquine against dormant liver-stage forms (hypnozoites), however drug resistance confronts this regimen and threatens malaria control programs. Understanding the basis of P. vivax chloroquine resistance (CQR) will inform drug discovery and malaria control. Here we investigate the genetics of P. vivax CQR by a cross of parasites differing in drug response. Gametocytogenesis, mosquito infection, and progeny production are performed with mixed parasite populations in nonhuman primates, as methods for P. vivax cloning and in vitro cultivation remain unavailable. Linkage mapping of progeny surviving >15 mg/kg CQ identifies a 76 kb region in chromosome 1 including pvcrt, an ortholog of the Plasmodium falciparum CQR transporter gene. Transcriptional analysis supports upregulated pvcrt expression as a mechanism of CQR.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

ApiAP2 Transcription Factors in Apicomplexan Parasites

Apicomplexan parasites are protozoan organisms that are characterised by complex life cycles and they include medically important species, such as the malaria parasite Plasmodium and the causative agents of toxoplasmosis (Toxoplasma gondii) and cryptosporidiosis (Cryptosporidium spp.). Apicomplexan parasites can infect one or more hosts, in which they differentiate into several morphologically and metabolically distinct life cycle stages. These developmental transitions rely on changes in gene expression. In the last few years, the important roles of different members of the ApiAP2 transcription factor family in regulating life cycle transitions and other aspects of parasite biology have become apparent. Here, we review recent progress in our understanding of the different members of the ApiAP2 transcription factor family in apicomplexan parasites.

Molecular and cellular interactions defining the tropism of Plasmodium vivax for reticulocytes

Plasmodium vivax is uniquely restricted to invading reticulocytes, the youngest of red blood cells. Parasite invasion relies on the sequential deployment of multiple parasite invasion ligands. Correct targeting of the host reticulocyte is mediated by two families of invasion ligands: the reticulocyte binding proteins (RBPs) and erythrocyte binding proteins (EBPs). The Duffy receptor has long been established as a key determinant for P. vivax invasion. However, recently, the RBP protein PvRBP2b has been shown to bind to transferrin receptor, which is expressed on reticulocytes but lost on normocytes, implicating the ligand-receptor in the reticulocyte tropism of P. vivax. Furthermore there is increasing evidence for P. vivax growth and sexual development in reticulocyte-enriched tissues such as the bone marrow.

In silico analysis of putative dormancy genes in Plasmodium vivax

Plasmodium vivax is the most widely spread species causing human malaria. The control of malaria caused by P. vivax has been largely hampered by its ability to develop a dormant liver stage that can generate a new blood infection at different periods of time. Unfortunately, the mechanisms of dormancy in P. vivax have not been thoroughly elucidated to date. In this study, the putative dormancy genes were analyzed to select genes with less genetic variability to maintain the function of relapsing. Expression data concerning these genes were searched to support the selection. Protein interactions among selected gene products were identified based on known and predicted protein-protein interaction using String database. Potentially interacting proteins (n = 15) were used to propose a mechanism involved in dormancy based on the differential vesicular transport due to the iron available in the hepatocyte.

Genome-wide SNP and InDel mutations in Mycobacterium tuberculosis associated with rifampicin and isoniazid resistance

OBJECTIVE

Multiple resistances to isoniazid and rifampicin lead to the majority of death associated with M. tuberculosis infection. This study aimed to characterize the single nucleotide polymorphisms (SNPs) and insertion and deletion (InDel) mutations associated with isoniazid and rifampicin resistance.

METHODS

The M. tuberculosis strain H37Rv was cultured and treated with isoniazid or rifampicin for generations. Total DNA samples from different generations were extracted for construction of DNA library, and the SNP and InDel mutation in different samples were detected by whole genome sequencing. Bioinformatics analysis such as phylogenetic tree and heap map were also performed.

RESULTS

Totally 58 nonsynonymous SNP mutations, 64 synonymous SNP mutations, and 99 SNP mutations in intergenic regions were detected in M. tuberculosis strains treated with rifampicin or isoniazid. Seven InDel mutations were found in the intergenic regions, and also six frameshift InDel mutation and three non-frameshift InDel mutations were also characterized. The phylogenetic tree showed clustering of all samples into three main subgroups. A great number of known and newly identified genes associated with drug resistance were detected in M. tuberculosis, showing distinct mutation patterns.

CONCLUSION

By whole genome sequencing, many genetic mutations in both known and new genes associated with isoniazid and rifampicin resistance were characterized in M. tuberculosis.

Immunomics analysis of Babesia microti protein markers by high-throughput screening assay

Babesia microti is a protozoan considered to be a major etiological agent of emerging human babesiosis. It imposes an increasing public-health threat and can be overlooked because of low parasitemia or mixed infection with other pathogens. More sensitive and specific antigens are needed to improve the diagnosis of babesiosis. To screen the immune diagnostic antigens of B. microti, 204 sequences from homologue proteins between B. microti and B. bovis genome sequences in PiroplasmaDB were selected. The high throughput cloned and expressed B. microti proteins were screened with the sera from the BALB/c mice infected by B. microti using protein arrays. Ten (5.9%, 10/169) highly immunoreactive proteins were identified, and most (80%, 8/10) of these highly immunoreactive proteins had not been characterized before, making them potentially useful as candidate antigens for the development of diagnostic tools for babesiosis.

Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics

During the last decade, the vast omics field has revolutionized biological research, especially the genomics, transcriptomics and proteomics branches, as technological tools become available to the field researcher and allow difficult question-driven studies to be addressed. Parasitology has greatly benefited from next generation sequencing (NGS) projects, which have resulted in a broadened comprehension of basic parasite molecular biology, ecology and epidemiology. Malariology is one example where application of this technology has greatly contributed to a better understanding of Plasmodium spp. biology and host-parasite interactions. Among the several parasite species that cause human malaria, the neglected Plasmodium vivax presents great research challenges, as in vitro culturing is not yet feasible and functional assays are heavily limited. Therefore, there are gaps in our P. vivax biology knowledge that affect decisions for control policies aiming to eradicate vivax malaria in the near future. In this review, we provide a snapshot of key discoveries already achieved in P. vivax sequencing projects, focusing on developments, hurdles, and limitations currently faced by the research community, as well as perspectives on future vivax malaria research.

Malaria Epidemiology at the Clone Level

Genotyping to distinguish between parasite clones is nowadays a standard in many molecular epidemiological studies of malaria. It has become crucial in drug trials and to follow individual clones in epidemiological studies, and to understand how drug resistance emerges and spreads. Here, we review the applications of the increasingly available genotyping tools and whole-genome sequencing data, and argue for a better integration of population genetics findings into malaria-control strategies.

Genome-wide diversity and differentiation in New World populations of the human malaria parasite Plasmodium vivax

Background

The Americas were the last continent colonized by humans carrying malaria parasites. Plasmodium falciparum from the New World shows very little genetic diversity and greater linkage disequilibrium, compared with its African counterparts, and is clearly subdivided into local, highly divergent populations. However, limited available data have revealed extensive genetic diversity in American populations of another major human malaria parasite, P. vivax.

Methods

We used an improved sample preparation strategy and next-generation sequencing to characterize 9 high-quality P. vivax genome sequences from northwestern Brazil. These new data were compared with publicly available sequences from recently sampled clinical P. vivax isolates from Brazil (BRA, total n = 11 sequences), Peru (PER, n = 23), Colombia (COL, n = 31), and Mexico (MEX, n = 19).

Principal findings/Conclusions

We found that New World populations of P. vivax are as diverse (nucleotide diversity π between 5.2 × 10⁻⁴ and 6.2 × 10⁻⁴) as P. vivax populations from Southeast Asia, where malaria transmission is substantially more intense. They display several non-synonymous nucleotide substitutions (some of them previously undescribed) in genes known or suspected to be involved in antimalarial drug resistance, such as dhfr, dhps, mdr1, mrp1, and mrp-2, but not in the chloroquine resistance transporter ortholog (crt-o) gene. Moreover, P. vivax in the Americas is much less geographically substructured than local P. falciparum populations, with relatively little between-population genome-wide differentiation (pairwise F_ST values ranging between 0.025 and 0.092). Finally, P. vivax populations show a rapid decline in linkage disequilibrium with increasing distance between pairs of polymorphic sites, consistent with very frequent outcrossing. We hypothesize that the high diversity of present-day P. vivax lineages in the Americas originated from successive migratory waves and subsequent admixture between parasite lineages from geographically diverse sites. Further genome-wide analyses are required to test the demographic scenario suggested by our data.

Plasmodium vivax is the most common human malaria parasite in the Americas, but how and when this species arrived in the New World remains unclear. Here we describe high-quality whole-genome sequence data for nine P. vivax isolates from Brazil, a country that accounts for 37% of the malaria burden in this continent, and compare these data with additional publicly available P. vivax genomes from Brazil, Peru, Colombia, and Mexico. P. vivax populations from the New World were found to be as diverse as their counterparts from areas with substantially higher malaria transmission, such as Southeast Asia, and to carry several non-synonymous substitutions in candidate drug-resistance genes. Moreover, genome-wide patterns of linkage disequilibrium between pairs of polymorphic sites are consistent with very frequent outcrossing in these populations. Interestingly, local P. vivax is more polymorphic, with less between-population differentiation, than sympatric populations of P. falciparum, possibly as a result of different demographic histories of these two species in the Americas. We hypothesize that local P. vivax lineages originated from successive migratory waves and subsequent admixture between parasites from geographically diverse sites.

Further evaluation of the NWF filter for the purification of Plasmodium vivax-infected erythrocytes

BACKGROUND

Isolation of Plasmodium-infected red blood cells (iRBCs) from clinical blood samples is often required for experiments, such as ex vivo drug assays, in vitro invasion assays and genome sequencing. Current methods for removing white blood cells (WBCs) from malaria-infected blood are time-consuming or costly. A prototype non-woven fabric (NWF) filter was developed for the purification of iRBCs, which showed great efficiency for removing WBCs in a pilot study. Previous work was performed with prototype filters optimized for processing 5-10 mL of blood. With the commercialization of the filters, this study aims to evaluate the efficiency and suitability of the commercial NWF filter for the purification of Plasmodium vivax-infected RBCs in smaller volumes of blood and to compare its performance with that of Plasmodipur^® filters.

METHODS

Forty-three clinical P. vivax blood samples taken from symptomatic patients attending malaria clinics at the China-Myanmar border were processed using the NWF filters in a nearby field laboratory. The numbers of WBCs and iRBCs and morphology of P. vivax parasites in the blood samples before and after NWF filtration were compared. The viability of P. vivax parasites after filtration from 27 blood samples was examined by in vitro short-term culture. In addition, the effectiveness of the NWF filter for removing WBCs was compared with that of the Plasmodipur^® filter in six P. vivax blood samples.

RESULTS

Filtration of 1-2 mL of P. vivax-infected blood with the NWF filter removed 99.68% WBCs. The densities of total iRBCs, ring and trophozoite stages before and after filtration were not significantly different (P > 0.05). However, the recovery rates of schizont- and gametocyte-infected RBCs, which were minor parasite stages in the clinical samples, were relatively low. After filtration, the P. vivax parasites did not show apparent morphological changes. Culture of 27 P. vivax-infected blood samples after filtration showed that parasites successfully matured into the schizont stage. The WBC removal rates and iRBC recovery rates were not significantly different between the NWF and Plasmodipur^® filters (P > 0.05).

CONCLUSIONS

When tested with 1-2 mL of P. vivax-infected blood, the NWF filter could effectively remove WBCs and the recovery rates for ring- and trophozoite-iRBCs were high. P. vivax parasites after filtration could be successfully cultured in vitro to reach maturity. The performance of the NWF and Plasmodipur^® filters for removing WBCs and recovering iRBCs was comparable.

Genomic variation in Plasmodium vivax malaria reveals regions under selective pressure

BACKGROUND

Although Plasmodium vivax contributes to almost half of all malaria cases outside Africa, it has been relatively neglected compared to the more deadly P. falciparum. It is known that P. vivax populations possess high genetic diversity, differing geographically potentially due to different vector species, host genetics and environmental factors.

RESULTS

We analysed the high-quality genomic data for 46 P. vivax isolates spanning 10 countries across 4 continents. Using population genetic methods we identified hotspots of selection pressure, including the previously reported MRP1 and DHPS genes, both putative drug resistance loci. Extra copies and deletions in the promoter region of another drug resistance candidate, MDR1 gene, and duplications in the Duffy binding protein gene (PvDBP) potentially involved in erythrocyte invasion, were also identified. For surveillance applications, continental-informative markers were found in putative drug resistance loci, and we show that organellar polymorphisms could classify P. vivax populations across continents and differentiate between Plasmodia spp.

CONCLUSIONS

This study has shown that genomic diversity that lies within and between P. vivax populations can be used to elucidate potential drug resistance and invasion mechanisms, as well as facilitate the molecular barcoding of the parasite for surveillance applications.

A systematic review on malaria sero-epidemiology studies in the Brazilian Amazon: insights into immunological markers for exposure and protection

Background

Considerable success in reducing malaria incidence and mortality has been achieved in Brazil, leading to discussions over the possibility of moving towards elimination. However, more than reporting and counting clinical cases, elimination will require the use of efficient tools and strategies for measuring transmission dynamics and detecting the infectious reservoir as the primary indicators of interest for surveillance and evaluation. Because acquisition and maintenance of anti-malarial antibodies depend on parasite exposure, seroprevalence rates could be used as a reliable tool for assessing malaria endemicity and an adjunct measure for monitoring transmission in a rapid and cost-effective manner.

Methods

This systematic review synthesizes the existing literature on seroprevalence of malaria in the Brazilian Amazon Basin. Different study designs (cross-sectional surveys and longitudinal studies) with reported serological results in well-defined Brazilian populations were considered. Medline (via PubMed), EMBASE and LILACS databases were screened and the articles were included per established selection criteria. Data extraction was performed by two authors and a modified critical appraisal tool was applied to assess the quality and completeness of cross-sectional studies regarding defined variables of interest.

Results

From 220 single records identified, 23 studies were included in this systematic review for the qualitative synthesis. Five studies reported serology results on Plasmodium falciparum, 14 papers assessed Plasmodium vivax and four articles reported results on both Plasmodium species. Considerable heterogeneity among the evaluated malarial antigens, including sporozoite and blood stage antigens, was observed. The majority of recent studies analysed IgG responses against P. vivax antigens reflecting the species distribution pattern in Brazil over the last decades. Most of the published papers were cross-sectional surveys (73.9%) and only six cohort studies were included in this review. Three studies pointed to an association between antibodies against circumsporozoite protein of both P. falciparum and P. vivax and malaria exposure. Furthermore, five out 13 cross-sectional studies evidenced a positive association between IgG antibodies to the conserved 19-kDa C-terminal region of the merozoite surface protein 1 of P. vivax (PvMSP1₁₉) and malaria exposure.

Conclusions

This systematic review identifies potential biomarkers of P. falciparum and P. vivax exposure in areas with variable and unstable malaria transmission in Brazil. However, this study highlights the need for standardization of further studies to provide an ideal monitoring tool to evaluate trends in malaria transmission and the effectiveness of malaria intervention programmes in Brazil. Moreover, the score-based weighted tool developed and used in this study still requires further validation.

Whole-genome sequencing of a Plasmodium vivax clinical isolate exhibits geographical characteristics and high genetic variation in China-Myanmar border area

Background

Currently in China, the trend of Plasmodium vivax cases imported from Southeast Asia was increased especially in the China-Myanmar border area. Driven by the increase in P. vivax cases and stronger need for vaccine and drug development, several P. vivax isolates genome sequencing projects are underway. However, little is known about the genetic variability in this area until now.

Results

The sequencing of the first P. vivax isolate from China-Myanmar border area (CMB-1) generated 120 million paired-end reads. A percentage of 10.6 of the quality-evaluated reads were aligned onto 99.9% of the reference strain Sal I genome in 62-fold coverage with an average of 4.8 SNPs per kb. We present a 539-SNP marker data set for P. vivax that can identify different parasites from different geographic origins with high sensitivity. We also identified exceptionally high levels of genetic variability in members of multigene families such as RBP, SERA, vir, MSP3 and AP2. The de-novo assembly yielded a database composed of 8,409 contigs with N50 lengths of 6.6 kb and revealed 661 novel predicted genes including 78 vir genes, suggesting a greater functional variation in P. vivax from this area.

Conclusion

Our result contributes to a better understanding of P. vivax genetic variation, and provides a fundamental basis for the geographic differentiation of vivax malaria from China-Myanmar border area using a direct sequencing approach without leukocyte depletion. This novel sequencing method can be used as an essential tool for the genomic research of P. vivax in the near future.

Electronic supplementary material

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

Design, construction and validation of a Plasmodium vivax microarray for the transcriptome profiling of clinical isolates

High density oligonucleotide microarrays have been used on Plasmodium vivax field isolates to estimate whole genome expression. However, no microarray platform has been experimentally optimized for studying the transcriptome of field isolates. In the present study, we adopted both bioinformatics and experimental testing approaches to select best optimized probes suitable for detecting parasite transcripts from field samples and included them in designing a custom 15K P. vivax microarray. This microarray has long oligonucleotide probes (60mer) that were in-situ synthesized onto glass slides using Agilent SurePrint technology and has been developed into an 8X15K format (8 identical arrays on a single slide). Probes in this array were experimentally validated and represents 4180 P. vivax genes in sense orientation, of which 1219 genes have also probes in antisense orientation. Validation of the 15K array by using field samples (n=14) has shown 99% of parasite transcript detection from any of the samples. Correlation analysis between duplicate probes (n=85) present in the arrays showed perfect correlation (r²=0.98) indicating the reproducibility. Multiple probes representing the same gene exhibited similar kind of expression pattern across the samples (positive correlation, r≥0.6). Comparison of hybridization data with the previous studies and quantitative real-time PCR experiments were performed to highlight the microarray validation procedure. This array is unique in its design, and results indicate that the array is sensitive and reproducible. Hence, this microarray could be a valuable functional genomics tool to generate reliable expression data from P. vivax field isolates.

Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia

Cambodia, in which both Plasmodium vivax and Plasmodium falciparum are endemic, has been the focus of numerous malaria-control interventions, resulting in a marked decline in overall malaria incidence. Despite this decline, the number of P vivax cases has actually increased. To understand better the factors underlying this resilience, we compared the genetic responses of the two species to recent selective pressures. We sequenced and studied the genomes of 70 P vivax and 80 P falciparum isolates collected between 2009 and 2013. We found that although P falciparum has undergone population fracturing, the coendemic P vivax population has grown undisrupted, resulting in a larger effective population size, no discernable population structure, and frequent multiclonal infections. Signatures of selection suggest recent, species-specific evolutionary differences. Particularly, in contrast to P falciparum, P vivax transcription factors, chromatin modifiers, and histone deacetylases have undergone strong directional selection, including a particularly strong selective sweep at an AP2 transcription factor. Together, our findings point to different population-level adaptive mechanisms used by P vivax and P falciparum parasites. Although population substructuring in P falciparum has resulted in clonal outgrowths of resistant parasites, P vivax may use a nuanced transcriptional regulatory approach to population maintenance, enabling it to preserve a larger, more diverse population better suited to facing selective threats. We conclude that transcriptional control may underlie P vivax's resilience to malaria control measures. Novel strategies to target such processes are likely required to eradicate P vivax and achieve malaria elimination.

A new Plasmodium vivax reference sequence with improved assembly of the subtelomeres reveals an abundance of pir genes

Plasmodium vivax is now the predominant cause of malaria in the Asia-Pacific, South America and Horn of Africa. Laboratory studies of this species are constrained by the inability to maintain the parasite in continuous ex vivo culture, but genomic approaches provide an alternative and complementary avenue to investigate the parasite’s biology and epidemiology. To date, molecular studies of P. vivax have relied on the Salvador-I reference genome sequence, derived from a monkey-adapted strain from South America. However, the Salvador-I reference remains highly fragmented with over 2500 unassembled scaffolds.  Using high-depth Illumina sequence data, we assembled and annotated a new reference sequence, PvP01, sourced directly from a patient from Papua Indonesia. Draft assemblies of isolates from China (PvC01) and Thailand (PvT01) were also prepared for comparative purposes. The quality of the PvP01 assembly is improved greatly over Salvador-I, with fragmentation reduced to 226 scaffolds. Detailed manual curation has ensured highly comprehensive annotation, with functions attributed to 58% core genes in PvP01 versus 38% in Salvador-I. The assemblies of PvP01, PvC01 and PvT01 are larger than that of Salvador-I (28-30 versus 27 Mb), owing to improved assembly of the subtelomeres.  An extensive repertoire of over 1200 Plasmodium interspersed repeat ( pir) genes were identified in PvP01 compared to 346 in Salvador-I, suggesting a vital role in parasite survival or development. The manually curated PvP01 reference and PvC01 and PvT01 draft assemblies are important new resources to study vivax malaria. PvP01 is maintained at GeneDB and ongoing curation will ensure continual improvements in assembly and annotation quality.

Genomic analysis of local variation and recent evolution in Plasmodium vivax

The widespread distribution and relapsing nature of Plasmodium vivax infection present major challenges for the elimination of malaria. To characterize the genetic diversity of this parasite in individual infections and across the population, we performed deep genome sequencing of >200 clinical samples collected across the Asia-Pacific region and analyzed data on >300,000 SNPs and nine regions of the genome with large copy number variations. Individual infections showed complex patterns of genetic structure, with variation not only in the number of dominant clones but also in their level of relatedness and inbreeding. At the population level, we observed strong signals of recent evolutionary selection both in known drug resistance genes and at new loci, and these varied markedly between geographical locations. These findings demonstrate a dynamic landscape of local evolutionary adaptation in the parasite population and provide a foundation for genomic surveillance to guide effective strategies for control and elimination of P. vivax.

An immunomics approach for the analysis of natural antibody responses to Plasmodium vivax infection

High throughput immunomics is a powerful platform to discover potential targets of host immunity and develop diagnostic tests for infectious diseases. We screened the sera of Plasmodium vivax-exposed individuals to profile the antibody response to blood-stage antigens of P. vivax using a P. vivax protein microarray. A total of 1936 genes encoding the P. vivax proteins were expressed, printed and screened with sera from P. vivax-exposed individuals and normal subjects. Total of 151 (7.8% of the 1936 targets) highly immunoreactive antigens were identified, including five well-characterized antigens of P. vivax (ETRAMP11.2, Pv34, SUB1, RAP2 and MSP4). Among the highly immunoreactive antigens, 5 antigens were predicted as adhesins by MAAP, and 11 antigens were predicted as merozoite invasion-related proteins based on homology with P. falciparum proteins. There are 40 proteins that have serodiagnostic potential for antibody surveillance. These novel Plasmodium antigens identified provide the clues for understanding host immune response to P. vivax infection and the development of antibody surveillance tools.

Advances in genetics and genomics: use and limitations in achieving malaria elimination goals

Success of the global research agenda towards eradication of malaria will depend on the development of new tools, including drugs, vaccines, insecticides and diagnostics. Genetic and genomic information now available for the malaria parasites, their mosquito vectors and human host, can be harnessed to both develop these tools and monitor their effectiveness. Here we review and provide specific examples of current technological advances and how these genetic and genomic tools have increased our knowledge of host, parasite and vector biology in relation to malaria elimination and in turn enhanced the potential to reach that goal. We then discuss limitations of these tools and future prospects for the successful achievement of global malaria elimination goals.

New insights into the Plasmodium vivax transcriptome using RNA-Seq

Historically seen as a benign disease, it is now becoming clear that Plasmodium vivax can cause significant morbidity. Effective control strategies targeting P. vivax malaria is hindered by our limited understanding of vivax biology. Here we established the P. vivax transcriptome of the Intraerythrocytic Developmental Cycle (IDC) of two clinical isolates in high resolution by Illumina HiSeq platform. The detailed map of transcriptome generates new insights into regulatory mechanisms of individual genes and reveals their intimate relationship with specific biological functions. A transcriptional hotspot of vir genes observed on chromosome 2 suggests a potential active site modulating immune evasion of the Plasmodium parasite across patients. Compared to other eukaryotes, P. vivax genes tend to have unusually long 5′ untranslated regions and also present multiple transcription start sites. In contrast, alternative splicing is rare in P. vivax but its association with the late schizont stage suggests some of its significance for gene function. The newly identified transcripts, including up to 179 vir like genes and 3018 noncoding RNAs suggest an important role of these gene/transcript classes in strain specific transcriptional regulation.

Whole Genome Sequencing of Field Isolates Reveals Extensive Genetic Diversity in Plasmodium vivax from Colombia

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.

Expression of Plasmodium vivax crt-o Is Related to Parasite Stage but Not Ex Vivo Chloroquine Susceptibility

Chloroquine (CQ)-resistant Plasmodium vivax is present in most countries where P. vivax infection is endemic, but the underlying molecular mechanisms responsible remain unknown. Increased expression of P. vivaxcrt-o (pvcrt-o) has been correlated with in vivo CQ resistance in an area with low-grade resistance. We assessed pvcrt-o expression in isolates from Papua (Indonesia), where P. vivax is highly CQ resistant. Ex vivo drug susceptibilities to CQ, amodiaquine, piperaquine, mefloquine, and artesunate were determined using a modified schizont maturation assay. Expression levels of pvcrt-o were measured using a novel real-time quantitative reverse transcription-PCR method. Large variations in pvcrt-o expression were observed across the 51 isolates evaluated, with the fold change in expression level ranging from 0.01 to 59 relative to that seen with the P. vivax β-tubulin gene and from 0.01 to 24 relative to that seen with the P. vivax aldolase gene. Expression was significantly higher in isolates with the majority of parasites at the ring stage of development (median fold change, 1.7) compared to those at the trophozoite stage (median fold change, 0.5; P < 0.001). Twenty-nine isolates fulfilled the criteria for ex vivo drug susceptibility testing and showed high variability in CQ responses (median, 107.9 [range, 6.5 to 345.7] nM). After controlling for the parasite stage, we found that pvcrt-o expression levels did not correlate with the ex vivo response to CQ or with that to any of the other antimalarials tested. Our results highlight the importance of development-stage composition for measuring pvcrt-o expression and suggest that pvcrt-o transcription is not a primary determinant of ex vivo drug susceptibility. A comprehensive transcriptomic approach is warranted for an in-depth investigation of the role of gene expression levels and P. vivax drug resistance.

Emerging Plasmodium vivax resistance to chloroquine in South America: an overview

The global emergence of Plasmodium vivax strains resistant to chloroquine (CQ) since the late 1980s is complicating the current international efforts for malaria control and elimination. Furthermore, CQ-resistant vivax malaria has already reached an alarming prevalence in Indonesia, East Timor and Papua New Guinea. More recently, in vivo studies have documented CQ-resistant P. vivax infections in Guyana, Peru and Brazil. Here, we summarise the available data on CQ resistance across P. vivax-endemic areas of Latin America by combining published in vivo and in vitro studies. We also review the current knowledge regarding the molecular mechanisms of CQ resistance in P. vivax and the prospects for developing and standardising reliable molecular markers of drug resistance. Finally, we discuss how the Worldwide Antimalarial Resistance Network, an international collaborative effort involving malaria experts from all continents, might contribute to the current regional efforts to map CQ-resistant vivax malaria in South America.

Next-Generation Sequencing of Plasmodium vivax Patient Samples Shows Evidence of Direct Evolution in Drug-Resistance Genes

Understanding the mechanisms of drug resistance in Plasmodium vivax, the parasite that causes the most widespread form of human malaria, is complicated by the lack of a suitable long-term cell culture system for this parasite. In contrast to P. falciparum, which can be more readily manipulated in the laboratory, insights about parasite biology need to be inferred from human studies. Here we analyze the genomes of parasites within 10 human P. vivax infections from the Peruvian Amazon. Using next-generation sequencing we show that some P. vivax infections analyzed from the region are likely polyclonal. Despite their polyclonality we observe limited parasite genetic diversity by showing that three or fewer haplotypes comprise 94% of the examined genomes, suggesting the recent introduction of parasites into this geographic region. In contrast we find more than three haplotypes in putative drug-resistance genes, including the gene encoding dihydrofolate reductase-thymidylate synthase and the P. vivax multidrug resistance associated transporter, suggesting that resistance mutations have arisen independently. Additionally, several drug-resistance genes are located in genomic regions with evidence of increased copy number. Our data suggest that whole genome sequencing of malaria parasites from patients may provide more insight about the evolution of drug resistance than genetic linkage or association studies, especially in geographical regions with limited parasite genetic diversity.

One Size Doesn't Fit All - RefEditor: Building Personalized Diploid Reference Genome to Improve Read Mapping and Genotype Calling in Next Generation Sequencing Studies

With rapid decline of the sequencing cost, researchers today rush to embrace whole genome sequencing (WGS), or whole exome sequencing (WES) approach as the next powerful tool for relating genetic variants to human diseases and phenotypes. A fundamental step in analyzing WGS and WES data is mapping short sequencing reads back to the reference genome. This is an important issue because incorrectly mapped reads affect the downstream variant discovery, genotype calling and association analysis. Although many read mapping algorithms have been developed, the majority of them uses the universal reference genome and do not take sequence variants into consideration. Given that genetic variants are ubiquitous, it is highly desirable if they can be factored into the read mapping procedure. In this work, we developed a novel strategy that utilizes genotypes obtained a priori to customize the universal haploid reference genome into a personalized diploid reference genome. The new strategy is implemented in a program named RefEditor. When applying RefEditor to real data, we achieved encouraging improvements in read mapping, variant discovery and genotype calling. Compared to standard approaches, RefEditor can significantly increase genotype calling consistency (from 43% to 61% at 4X coverage; from 82% to 92% at 20X coverage) and reduce Mendelian inconsistency across various sequencing depths. Because many WGS and WES studies are conducted on cohorts that have been genotyped using array-based genotyping platforms previously or concurrently, we believe the proposed strategy will be of high value in practice, which can also be applied to the scenario where multiple NGS experiments are conducted on the same cohort. The RefEditor sources are available at https://github.com/superyuan/refeditor.

Predicting and exploring network components involved in pathogenesis in the malaria parasite via novel subnetwork alignments

BACKGROUND

Malaria is a major health threat, affecting over 40% of the world's population. The latest report released by the World Health Organization estimated about 207 million cases of malaria infection, and about 627,000 deaths in 2012 alone. During the past decade, new therapeutic targets have been identified and are at various stages of characterization, thanks to the emerging omics-based technologies. However, the mechanism of malaria pathogenesis remains largely unknown. In this paper, we employ a novel neighborhood subnetwork alignment approach to identify network components that are potentially involved in pathogenesis.

RESULTS

Our module-based subnetwork alignment approach identified 24 functional homologs of pathogenesis-related proteins in the malaria parasite P. falciparum, using the protein-protein interaction networks in Escherichia coli as references. Eighteen out of these 24 proteins are associated with 418 other proteins that are related to DNA replication, transcriptional regulation, translation, signaling, metabolism, cell cycle regulation, as well as cytoadherence and entry to the host.

CONCLUSIONS

The subnetwork alignments and subsequent protein-protein association network mining predicted a group of malarial proteins that may be involved in parasite development and parasite-host interaction, opening a new systems-level view of parasite pathogenesis and virulence.

Uncovering the transmission dynamics of Plasmodium vivax using population genetics

Population genetic analysis of malaria parasites has the power to reveal key insights into malaria epidemiology and transmission dynamics with the potential to deliver tools to support control and elimination efforts. Analyses of parasite genetic diversity have suggested that Plasmodium vivax populations are more genetically diverse and less structured than those of Plasmodium falciparum indicating that P. vivax may be a more ancient parasite of humans and/or less susceptible to population bottlenecks, as well as more efficient at disseminating its genes. These population genetic insights into P. vivax transmission dynamics provide an explanation for its relative resilience to control efforts. Here, we describe current knowledge on P. vivax population genetic structure, its relevance to understanding transmission patterns and relapse and how this information can inform malaria control and elimination programmes.

Heavy path mining of protein-protein associations in the malaria parasite

Annotating and understanding the function of proteins and other elements in a genome can be difficult in the absence of a well-studied and evolutionarily close relative. The causative agent of malaria, one of the oldest and most deadly global infectious diseases, is a good example of this problem. The burden of malaria is huge and there is a pressing need for new, more effective antimalarial strategies. However, techniques such as homology-dependent annotation transfer are severely impaired in this parasite because there are no well-understood close relatives. To circumvent this approach we developed a network-based method that uses a heavy path network-mining algorithm. We uncovered the protein-protein associations that are implicated in important cellular processes including genome integrity, DNA repair, transcriptional regulation, invasion, and pathogenesis, thus demonstrating the utility of this method. The URL of the source code for super-sequence mining method is http://www.cs.utsa.edu/~korkmaz/research/heavy-path-mining/.

Development of a single nucleotide polymorphism barcode to genotype Plasmodium vivax infections

Plasmodium vivax, one of the five species of Plasmodium parasites that cause human malaria, is responsible for 25–40% of malaria cases worldwide. Malaria global elimination efforts will benefit from accurate and effective genotyping tools that will provide insight into the population genetics and diversity of this parasite. The recent sequencing of P. vivax isolates from South America, Africa, and Asia presents a new opportunity by uncovering thousands of novel single nucleotide polymorphisms (SNPs). Genotyping a selection of these SNPs provides a robust, low-cost method of identifying parasite infections through their unique genetic signature or barcode. Based on our experience in generating a SNP barcode for P. falciparum using High Resolution Melting (HRM), we have developed a similar tool for P. vivax. We selected globally polymorphic SNPs from available P. vivax genome sequence data that were located in putatively selectively neutral sites (i.e., intergenic, intronic, or 4-fold degenerate coding). From these candidate SNPs we defined a barcode consisting of 42 SNPs. We analyzed the performance of the 42-SNP barcode on 87 P. vivax clinical samples from parasite populations in South America (Brazil, French Guiana), Africa (Ethiopia) and Asia (Sri Lanka). We found that the P. vivax barcode is robust, as it requires only a small quantity of DNA (limit of detection 0.3 ng/μl) to yield reproducible genotype calls, and detects polymorphic genotypes with high sensitivity. The markers are informative across all clinical samples evaluated (average minor allele frequency > 0.1). Population genetic and statistical analyses show the barcode captures high degrees of population diversity and differentiates geographically distinct populations. Our 42-SNP barcode provides a robust, informative, and standardized genetic marker set that accurately identifies a genomic signature for P. vivax infections.

Plasmodium vivax malaria is a major global public health problem, with nearly 2.5 billion people at risk for infection and approximately 132–391 million clinical infections annually. It has a wide geographical range, with a high disease burden in Asia, Central and South America, the Middle East, Oceania, and East Africa. Advances in sequencing technology and sample processing have made it possible to characterize the genetic diversity of P. vivax populations. This genetic variation provides a means to identify parasites by unique genetic signatures, or “barcodes.” We developed such a genetic barcode for P. vivax, composed of 42 robust and informative variants. Here we report its development and validation based on 87 clinical samples identified by microscopy to contain P. vivax from geographically diverse parasite populations from South America (Brazil, French Guiana), Africa (Ethiopia) and Asia (Sri Lanka). We show that the SNP barcode provides a genotyping tool that can be performed at low cost, providing a means to uniquely identify parasite infections and distinguish geographic origins, and that barcode data may offer new insights into P. vivax population structure and diversity.

Comparative analysis of field-isolate and monkey-adapted Plasmodium vivax genomes

Significant insights into the biology of Plasmodium vivax have been gained from the ability to successfully adapt human infections to non-human primates. P. vivax strains grown in monkeys serve as a renewable source of parasites for in vitro and ex vivo experimental studies and functional assays, or for studying in vivo the relapse characteristics, mosquito species compatibilities, drug susceptibility profiles or immune responses towards potential vaccine candidates. Despite the importance of these studies, little is known as to how adaptation to a different host species may influence the genome of P. vivax. In addition, it is unclear whether these monkey-adapted strains consist of a single clonal population of parasites or if they retain the multiclonal complexity commonly observed in field isolates. Here we compare the genome sequences of seven P. vivax strains adapted to New World monkeys with those of six human clinical isolates collected directly in the field. We show that the adaptation of P. vivax parasites to monkey hosts, and their subsequent propagation, did not result in significant modifications of their genome sequence and that these monkey-adapted strains recapitulate the genomic diversity of field isolates. Our analyses also reveal that these strains are not always genetically homogeneous and should be analyzed cautiously. Overall, our study provides a framework to better leverage this important research material and fully utilize this resource for improving our understanding of P. vivax biology.

In this study we compare the genome sequences of Plasmodium vivax collected directly from patients with those of parasites propagated in laboratory monkeys. We show that the adaptation and continuous propagation of Plasmodium vivax in monkeys does not induce systematic changes in the genome and, therefore, that these parasites constitute an unbiased resource for studying this important pathogen. Our analyses also reveal that some monkey-adapted Plasmodium vivax strains are not genetically homogenous and retain multiple genetically different parasites present in the original patient infection. Overall, our study confirms the utility of monkey-adapted Plasmodium vivax strains for malaria research but also shows that this resource should be analyzed cautiously as different samples of the same strain might provide different biological material.

A novel sterol regulatory element-binding protein gene (sreA) identified in penicillium digitatum is required for prochloraz resistance, full virulence and erg11 (cyp51) regulation

Penicilliumdigitatum is the most destructive postharvest pathogen of citrus fruits, causing fruit decay and economic loss. Additionally, control of the disease is further complicated by the emergence of drug-resistant strains due to the extensive use of triazole antifungal drugs. In this work, an orthologus gene encoding a putative sterol regulatory element-binding protein (SREBP) was identified in the genome of P. digitatum and named sreA. The putative SreA protein contains a conserved domain of unknown function (DUF2014) at its carboxyl terminus and a helix-loop-helix (HLH) leucine zipper DNA binding domain at its amino terminus, domains that are functionally associated with SREBP transcription factors. The deletion of sreA (ΔsreA) in a prochloraz-resistant strain (PdHS-F6) by Agrobacteriumtumefaciens-mediated transformation led to increased susceptibility to prochloraz and a significantly lower EC₅₀ value compared with the HS-F6 wild-type or complementation strain (COsreA). A virulence assay showed that the ΔsreA strain was defective in virulence towards citrus fruits, while the complementation of sreA could restore the virulence to a large extent. Further analysis by quantitative real-time PCR demonstrated that prochloraz-induced expression of cyp51A and cyp51B in PdHS-F6 was completely abolished in the ΔsreA strain. These results demonstrate that sreA is a critical transcription factor gene required for prochloraz resistance and full virulence in P. digitatum and is involved in the regulation of cyp51 expression.

The biology of Plasmodium vivax explored through genomics

Malaria is a mosquito-borne disease caused by the Plasmodium parasite. Of the four Plasmodium species that routinely cause human malaria, Plasmodium vivax is the most widespread species outside Africa, causing ∼18.9 million cases in 2012. P. vivax cannot be cultured continuously in vitro, which severely hampers research in nonendemic and endemic countries alike. Consequently, whole-genome sequencing has become an effective means to interrogate the biology of the P. vivax parasite. Our comparative genomic analysis of five P. vivax reference genomes and several whole-genome sequences of the closely related monkey malaria species P. cynomolgi has revealed an extraordinary level of genetic diversity and enabled characterization of novel multigene families and important single-copy genes. The generation of whole-genome sequences from multiple clinical isolates is also driving forward knowledge concerning the biology and evolution of the species. Understanding the biology of P. vivax is crucial to develop potential antimalarial drugs and vaccines and to achieve the goal of eliminating malaria.

A high resolution case study of a patient with recurrent Plasmodium vivax infections shows that relapses were caused by meiotic siblings

Plasmodium vivax infects a hundred million people annually and endangers 40% of the world's population. Unlike Plasmodium falciparum, P. vivax parasites can persist as a dormant stage in the liver, known as the hypnozoite, and these dormant forms can cause malaria relapses months or years after the initial mosquito bite. Here we analyze whole genome sequencing data from parasites in the blood of a patient who experienced consecutive P. vivax relapses over 33 months in a non-endemic country. By analyzing patterns of identity, read coverage, and the presence or absence of minor alleles in the initial polyclonal and subsequent monoclonal infections, we show that the parasites in the three infections are likely meiotic siblings. We infer that these siblings are descended from a single tetrad-like form that developed in the infecting mosquito midgut shortly after fertilization. In this natural cross we find the recombination rate for P. vivax to be 10 kb per centimorgan and we further observe areas of disequilibrium surrounding major drug resistance genes. Our data provide new strategies for studying multiclonal infections, which are common in all types of infectious diseases, and for distinguishing P. vivax relapses from reinfections in malaria endemic regions. This work provides a theoretical foundation for studies that aim to determine if new or existing drugs can provide a radical cure of P. vivax malaria.

Plasmodium vivax is capable of remaining dormant in the human liver for months to years after an initial infection, creating an asymptomatic human reservoir. This unique aspect of parasite biology makes eliminating P. vivax distinctly different from P. falciparum elimination, and yet very little is known about this dormant parasite stage. Lack of knowledge about the dormant liver stage prevents the creation of new drugs and public health interventions directed at P. vivax. In order to better understand this particular parasite stage, we used whole genome sequencing to analyze three sequential P. vivax infections, two of which could be definitively categorized as having arisen from dormant liver stages. Our whole genome sequencing data demonstrates that dormant liver stage parasites are closely related yet not, as had previously been postulated, identical. These data highlight the need for a new paradigm to investigate P. vivax dormant liver stages in order to design the next generation of P. vivax drugs and effective global health interventions.