Using mitochondrial genome sequences to track the origin of imported Plasmodium vivax infections diagnosed in the United States

Malaria Genomics, Vaccine Development, and Microbiome

Recent advances in malaria genetics and genomics have transformed many aspects of malaria research in areas of molecular evolution, epidemiology, transmission, host-parasite interaction, drug resistance, pathogenicity, and vaccine development. Here, in addition to introducing some background information on malaria parasite biology, parasite genetics/genomics, and genotyping methods, we discuss some applications of genetic and genomic approaches in vaccine development and in studying interactions with microbiota. Genetic and genomic data can be used to search for novel vaccine targets, design an effective vaccine strategy, identify protective antigens in a whole-organism vaccine, and evaluate the efficacy of a vaccine. Microbiota has been shown to influence disease outcomes and vaccine efficacy; studying the effects of microbiota in pathogenicity and immunity may provide information for disease control. Malaria genetics and genomics will continue to contribute greatly to many fields of malaria research.

Geographical origin of Plasmodium vivax in the Hainan Island, China: insights from mitochondrial genome

BACKGROUND

Hainan Province, China, has been an endemic region with high transmission of Plasmodium falciparum and Plasmodium vivax. Indigenous malaria caused by P. vivax was eliminated in Hainan in 2011, while imported vivax malaria remains. However, the geographical origin of P. vivax cases in Hainan remains unclear.

METHODS

Indigenous and imported P. vivax isolates (n = 45) were collected from Hainan Province, and the 6 kb mitochondrial genome was obtained. Nucleotide (π) and haplotype (h) diversity were estimated using DnaSP. The numbers of synonymous nucleotide substitutions per synonymous site (d_S) and nonsynonymous nucleotide substitutions per nonsynonymous site (d_N) were calculated using the SNAP program. Arlequin software was used to estimate the genetic diversity index and assess population differentiation. Bayesian phylogenetic analysis of P. vivax was performed using MrBayes. A haplotype network was generated using the NETWORK program.

RESULTS

In total, 983 complete mitochondrial genome sequences were collected, including 45 from this study and 938 publicly available from the NCBI. Thirty-three SNPs were identified, and 18 haplotypes were defined. The haplotype (0.834) and nucleotide (0.00061) diversity in the Hainan populations were higher than China's Anhui and Guizhou population, and the majority of pairwise F_ST values in Hainan exceeded 0.25, suggesting strong differentiation among most populations except in Southeast Asia. Most Hainan haplotypes were connected to South/East Asian and China's others haplotypes, but less connected with populations from China's Anhui and Guizhou provinces. Mitochondrial lineages of Hainan P. vivax belonged to clade 1 of four well-supported clades in a phylogenetic tree, most haplotypes of indigenous cases formed a subclade of clade 1, and the origin of seven imported cases (50%) could be inferred from the phylogenetic tree, but five imported cases (42.8%) could not be traced using the phylogenetic tree alone, necessitating epidemiological investigation.

CONCLUSIONS

Indigenous cases in Hainan display high genetic (haplotype and nucleotide) diversity. Haplotype network analysis also revealed most haplotypes in Hainan were connected to the Southeast Asian populations and divergence to a cluster of China's other populations. According to the mtDNA phylogenetic tree, some haplotypes were shared between geographic populations, and some haplotypes have formed lineages. Multiple tests are needed to further explore the origin and expansion of P. vivax populations.

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

A molecular barcode and web-based data analysis tool to identify imported Plasmodium vivax malaria

Traditionally, patient travel history has been used to distinguish imported from autochthonous malaria cases, but the dormant liver stages of Plasmodium vivax confound this approach. Molecular tools offer an alternative method to identify, and map imported cases. Using machine learning approaches incorporating hierarchical fixation index and decision tree analyses applied to 799 P. vivax genomes from 21 countries, we identified 33-SNP, 50-SNP and 55-SNP barcodes (GEO33, GEO50 and GEO55), with high capacity to predict the infection's country of origin. The Matthews correlation coefficient (MCC) for an existing, commonly applied 38-SNP barcode (BR38) exceeded 0.80 in 62% countries. The GEO panels outperformed BR38, with median MCCs > 0.80 in 90% countries at GEO33, and 95% at GEO50 and GEO55. An online, open-access, likelihood-based classifier framework was established to support data analysis (vivaxGEN-geo). The SNP selection and classifier methods can be readily amended for other use cases to support malaria control programs.

A PCR-RFLP Technique to Assess the Geographic Origin of Plasmodium falciparum Strains in Central America

The elimination of malaria requires strengthening diagnosis and offering adequate and timely treatment. Imported cases of falciparum malaria represent a major challenge for pre-elimination areas, such as Central America, where chloroquine and primaquine continue to be used as first-line treatment. The pfs47 gene has been previously described as a precise molecular marker to track the geographic origin of the parasite. The aim of this study was to design a simple and low-cost technique using the polymorphic region of pfs47 to assess the geographic origin of P. falciparum strains. A PCR-RFLP technique was developed and evaluated using the MseI enzyme that proved capable of discriminating, with reasonable precision, the geographical origin of the parasites. This method could be used by national surveillance laboratories and malaria elimination programs in countries such as Honduras and Nicaragua in cases of malaria where an origin outside the Central American isthmus is suspected.

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.

Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region

The Asia-Pacific region faces formidable challenges in achieving malaria elimination by the proposed target in 2030. Molecular surveillance of Plasmodium parasites can provide important information on malaria transmission and adaptation, which can inform national malaria control programmes (NMCPs) in decision-making processes. In November 2019 a parasite genotyping workshop was held in Jakarta, Indonesia, to review molecular approaches for parasite surveillance and explore ways in which these tools can be integrated into public health systems and inform policy. The meeting was attended by 70 participants from 8 malaria-endemic countries and partners of the Asia Pacific Malaria Elimination Network. The participants acknowledged the utility of multiple use cases for parasite genotyping including: quantifying the prevalence of drug resistant parasites, predicting risks of treatment failure, identifying major routes and reservoirs of infection, monitoring imported malaria and its contribution to local transmission, characterizing the origins and dynamics of malaria outbreaks, and estimating the frequency of Plasmodium vivax relapses. However, the priority of each use case varies with different endemic settings. Although a one-size-fits-all approach to molecular surveillance is unlikely to be applicable across the Asia-Pacific region, consensus on the spectrum of added-value activities will help support data sharing across national boundaries. Knowledge exchange is needed to establish local expertise in different laboratory-based methodologies and bioinformatics processes. Collaborative research involving local and international teams will help maximize the impact of analytical outputs on the operational needs of NMCPs. Research is also needed to explore the cost-effectiveness of genetic epidemiology for different use cases to help to leverage funding for wide-scale implementation. Engagement between NMCPs and local researchers will be critical throughout this process.

Dynamics of Plasmodium vivax populations in border areas of the Greater Mekong sub-region during malaria elimination

Background

Countries within the Greater Mekong Sub-region (GMS) of Southeast Asia have committed to eliminating malaria by 2030. Although the malaria situation has greatly improved, malaria transmission remains at international border regions. In some areas, Plasmodium vivax has become the predominant parasite. To gain a better understanding of transmission dynamics, knowledge on the changes of P. vivax populations after the scale-up of control interventions will guide more effective targeted control efforts.

Methods

This study investigated genetic diversity and population structures in 206 P. vivax clinical samples collected at two time points in two international border areas: the China-Myanmar border (CMB) (n = 50 in 2004 and n = 52 in 2016) and Thailand-Myanmar border (TMB) (n = 50 in 2012 and n = 54 in 2015). Parasites were genotyped using 10 microsatellite markers.

Results

Despite intensified control efforts, genetic diversity remained high (H_E = 0.66–0.86) and was not significantly different among the four populations (P > 0.05). Specifically, H_E slightly decreased from 0.76 in 2004 to 0.66 in 2016 at the CMB and increased from 0.80 in 2012 to 0.86 in 2015 at the TMB. The proportions of polyclonal infections varied significantly among the four populations (P < 0.05), and showed substantial decreases from 48.0% in 2004 to 23.7 at the CMB and from 40.0% in 2012 to 30.7% in 2015 at the TMB, with corresponding decreases in the multiplicity of infection. Consistent with the continuous decline of malaria incidence in the GMS over time, there were also increases in multilocus linkage disequilibrium, suggesting more fragmented and increasingly inbred parasite populations. There were considerable genetic differentiation and sub-division among the four tested populations. Temporal genetic differentiation was observed at each site (F_ST = 0.081 at the CMB and F_ST = 0.133 at the TMB). Various degrees of clustering were evident between the older parasite samples collected in 2004 at the CMB and the 2016 CMB and 2012 TMB populations, suggesting some of these parasites had shared ancestry. In contrast, the 2015 TMB population was genetically distinctive, which may reflect a process of population replacement. Whereas the effective population size (N_e) at the CMB showed a decrease from 4979 in 2004 to 3052 in 2016 with the infinite allele model, the N_e at the TMB experienced an increase from 6289 to 10,259.

Conclusions

With enhanced control efforts on malaria, P. vivax at the TMB and CMB showed considerable spatial and temporal differentiation, but the presence of large P. vivax reservoirs still sustained genetic diversity and transmission. These findings provide new insights into P. vivax transmission dynamics and population structure in these border areas of the GMS. Coordinated and integrated control efforts on both sides of international borders are essential to reach the goal of regional malaria elimination.

A molecular barcode to inform the geographical origin and transmission dynamics of Plasmodium vivax malaria

Although Plasmodium vivax parasites are the predominant cause of malaria outside of sub-Saharan Africa, they not always prioritised by elimination programmes. P. vivax is resilient and poses challenges through its ability to re-emerge from dormancy in the human liver. With observed growing drug-resistance and the increasing reports of life-threatening infections, new tools to inform elimination efforts are needed. In order to halt transmission, we need to better understand the dynamics of transmission, the movement of parasites, and the reservoirs of infection in order to design targeted interventions. The use of molecular genetics and epidemiology for tracking and studying malaria parasite populations has been applied successfully in P. falciparum species and here we sought to develop a molecular genetic tool for P. vivax. By assembling the largest set of P. vivax whole genome sequences (n = 433) spanning 17 countries, and applying a machine learning approach, we created a 71 SNP barcode with high predictive ability to identify geographic origin (91.4%). Further, due to the inclusion of markers for within population variability, the barcode may also distinguish local transmission networks. By using P. vivax data from a low-transmission setting in Malaysia, we demonstrate the potential ability to infer outbreak events. By characterising the barcoding SNP genotypes in P. vivax DNA sourced from UK travellers (n = 132) to ten malaria endemic countries predominantly not used in the barcode construction, we correctly predicted the geographic region of infection origin. Overall, the 71 SNP barcode outperforms previously published genotyping methods and when rolled-out within new portable platforms, is likely to be an invaluable tool for informing targeted interventions towards elimination of this resilient human malaria.

Using the Plasmodium mitochondrial genome for classifying mixed-species infections and inferring the geographical origin of P. falciparum parasites imported to the U.S

The ability to identify mixed-species infections and track the origin of Plasmodium parasites can further enhance the development of treatment and prevention recommendations as well as outbreak investigations. Here, we explore the utility of using the full Plasmodium mitochondrial genome to classify Plasmodium species, detect mixed infections, and infer the geographical origin of imported P. falciparum parasites to the United States (U.S.). Using the recently developed standardized, high-throughput Malaria Resistance Surveillance (MaRS) protocol, the full Plasmodium mitochondrial genomes of 265 malaria cases imported to the U.S. from 2014-2017 were sequenced and analyzed. P. falciparum infections were found in 94.7% (251/265) of samples. Five percent (14/265) of samples were identified as mixed- Plasmodium species or non-P. falciparum, including P. vivax, P. malariae, P. ovale curtisi, and P. ovale wallikeri. P. falciparum mitochondrial haplotypes analysis revealed greater than eighteen percent of samples to have at least two P. falciparum mitochondrial genome haplotypes, indicating either heteroplasmy or multi-clonal infections. Maximum-likelihood phylogenies of 912 P. falciparum mitochondrial genomes with known country origin were used to infer the geographical origin of thirteen samples from persons with unknown travel histories as: Africa (country unspecified) (n = 10), Ghana (n = 1), Southeast Asia (n = 1), and the Philippines (n = 1). We demonstrate the utility and current limitations of using the Plasmodium mitochondrial genome to classify samples with mixed-infections and infer the geographical origin of imported P. falciparum malaria cases to the U.S. with unknown travel history.

Mitochondrial genome sequence variation as a useful marker for assessing genetic heterogeneity among Cyclospora cayetanensis isolates and source-tracking

Background

Cyclospora cayetanensis is an important enteric pathogen, causing diarrhea and food-borne cyclosporiasis outbreaks. For effective outbreak identification and investigation, it is essential to rapidly assess the genetic heterogeneity of C. cayetanensis specimens from cluster cases and identify the likely occurrence of outbreaks.

Methods

In this study, we developed a quantitative PCR (qPCR) targeting the polymorphic link region between copies of the mitochondrial genome of C. cayetanensis, and evaluated the genetic heterogeneity among 36 specimens from six countries using melt curve, gel electrophoresis, and sequence analyses of the qPCR products.

Results

All specimens were amplified successfully in the qPCR and produced melt peaks with different Tm values in the melt curve analysis. In gel electrophoresis of the qPCR products, the specimens yielded bands of variable sizes. Nine genotypes were identified by DNA sequencing of the qPCR products. Geographical segregation of genotypes was observed among specimens analyzed, which could be useful in geographical source-tracking.

Conclusions

The length and nucleotide sequence variations in the mitochondrial genome marker allow rapid assessment of the genetic heterogeneity among C. cayetanensis specimens by melt curve, gel electrophoresis, or DNA sequence analysis of qPCR products. The sequence data generated could be helpful in the initial source-tracking of the pathogen.

Paralogs vs. genotypes? Variability of Babesia canis assessed by 18S rDNA and two mitochondrial markers

Canine babesiosis caused by Babesia canis sensu stricto became an emerging disease of dogs across Europe calling for attention also in countries where it was an only rare imported disease. An easy accessibility of molecular methods and the growing amount of sequencing data led to the description of intraspecific variability in 18S rDNA sequences designated as "genotypes". Using material from a homogenous cohort of dogs with microscopically confirmed canine babesiosis caused by B. canis, we evaluated Babesia intraspecific variability and amplification sensitivity of three different genes (18S rDNA, COI, Cytb) to assess their potential as diagnostic or phylogenetic markers. In raw sequencing data obtained, we observed at least 3 ambiguous positions in up to 86% of chromatograms within the ∼560 bp fragment of 18S rDNA suggesting the existence of several, not identical copies of this gene. Our COI haplotype analysis resulted in a star-like pattern indicating a recent origin of most haplotypes, but not supporting the existence of two dominant haplotypes. Similarly, the Cytb sequences obtained from samples with all variants of 18S rDNA were identical. We corroborate previous observations from three other European countries and bring the evidence of the existence of 18S rDNA paralogs in B. canis genome replacing currently used "genotype" theory.

Population Genetics and Molecular Epidemiology of Eukaryotes

Molecular epidemiology uses the distribution and organization of a pathogen's DNA to understand the distribution and determinants of disease. Since the biology of DNA for eukaryotic pathogens differs substantially from that of bacteria, the analytic approach to their molecular epidemiology can also differ. While many of the genotyping techniques presented earlier in this series, "Advances in Molecular Epidemiology of Infectious Diseases," can be applied to eukaryotes, the output must be interpreted in the light of how DNA is distributed from one generation to the next. In some cases, parasite populations can be evaluated in ways reminiscent of bacteria. They differ, however, when analyzed as sexually reproducing organisms, where all individuals are unique but the genetic composition of the population does not change unless a limited set of events occurs. It is these events (migration, mutation, nonrandom mating, selection, and genetic drift) that are of interest. At a given time, not all of them are likely to be equally important, so the list can easily be narrowed down to understand the driving forces behind the population as it is now and even what it will look like in the future. The main population characteristics measured to assess these events are differentiation and diversity, interpreted in the light of what is known about the population from observation. The population genetics of eukaryotes is important for planning and evaluation of control measures, surveillance, outbreak investigation, and monitoring of the development and spread of drug resistance. *This article is part of a curated collection.

Nationwide genetic surveillance of Plasmodium vivax in Papua New Guinea reveals heterogeneous transmission dynamics and routes of migration amongst subdivided populations

The Asia Pacific Leaders in Malaria Alliance (APLMA) have committed to eliminate malaria from the region by 2030. Papua New Guinea (PNG) has the highest malaria burden in the Asia-Pacific region but with the intensification of control efforts since 2005, transmission has been dramatically reduced and Plasmodium vivax is now the dominant malaria infection in some parts of the country. To gain a better understanding of the transmission dynamics and migration patterns of P. vivax in PNG, here we investigate population structure in eight geographically and ecologically distinct regions of the country. A total of 219 P. vivax isolates (16-30 per population) were successfully haplotyped using 10 microsatellite markers. A wide range of genetic diversity (H_e=0.37-0.87, R_s=3.60-7.58) and significant multilocus linkage disequilibrium (LD) was observed in six of the eight populations (I_A^S=0.08-0.15 p-value<0.05) reflecting a spectrum of transmission intensities across the country. Genetic differentiation between regions was evident (Jost's D=0.07-0.72), with increasing divergence of populations with geographic distance. Overall, P. vivax isolates clustered into three major genetic populations subdividing the Mainland lowland and coastal regions, the Islands and the Highlands. P. vivax gene flow follows major human migration routes, and there was higher gene flow amongst Mainland parasite populations than among Island populations. The Central Province (samples collected in villages close to the capital city, Port Moresby), acts as a sink for imported infections from the three major endemic areas. These insights into P. vivax transmission dynamics and population networks will inform targeted strategies to contain malaria infections and to prevent the spread of drug resistance in PNG.

Mitochondrial genome of Plasmodium vivax/simium detected in an endemic region for malaria in the Atlantic Forest of Espírito Santo state, Brazil: do mosquitoes, simians and humans harbour the same parasite?

BACKGROUND

The transmission of malaria in the extra-Amazonian regions of Brazil, although interrupted in the 1960s, has persisted to the present time in some areas of dense Atlantic Forest, with reports of cases characterized by particular transmission cycles and clinical presentations. Bromeliad-malaria, as it is named, is particularly frequent in the state of Espírito Santo, with Plasmodium vivax being the parasite commonly recognized as the aetiologic agent of human infections. With regard to the spatial and temporal distances between cases reported in this region, the transmission cycle does not fit the traditional malaria cycle. The existence of a zoonosis, with infected simians participating in the epidemiology, is therefore hypothesized. In the present study, transmission of bromeliad-malaria in Espírito Santo is investigated, based on the complete mitochondrial genome of DNA extracted from isolates of Plasmodium species, which had infected humans, a simian from the genus Allouata, and Anopheles mosquitoes. Plasmodium vivax/simium was identified in the samples by both nested PCR and real-time PCR. After amplification, the mitochondrial genome was completely sequenced and compared with a haplotype network which included all sequences of P. vivax/simium mitochondrial genomes sampled from humans and simians from all regions in Brazil.

RESULTS

The haplotype network indicates that humans and simians from the Atlantic Forest become infected by the same haplotype, but some isolates from humans are not identical to the simian isolate. In addition, the plasmodial DNA extracted from mosquitoes revealed sequences different from those obtained from simians, but similar to two isolates from humans.

CONCLUSIONS

These findings strengthen support for the hypothesis that in the Atlantic Forest, and especially in the state with the highest frequency of bromeliad-malaria in Brazil, parasites with similar molecular backgrounds are shared by humans and simians. The recognized identity between P. vivax and P. simium at the species level, the sharing of haplotypes, and the participation of the same vector in transmitting the infection to both host species indicate interspecies transference of the parasites. However, the intensity, frequency and direction of this transfer remain to be clarified.

Outbreak of human malaria caused by Plasmodium simium in the Atlantic Forest in Rio de Janeiro: a molecular epidemiological investigation

BACKGROUND

Malaria was eliminated from southern and southeastern Brazil over 50 years ago. However, an increasing number of autochthonous episodes attributed to Plasmodium vivax have recently been reported from the Atlantic Forest region of Rio de Janeiro state. As the P vivax-like non-human primate malaria parasite species Plasmodium simium is locally enzootic, we performed a molecular epidemiological investigation to determine whether zoonotic malaria transmission is occurring.

METHODS

We examined blood samples from patients presenting with signs or symptoms suggestive of malaria as well as from local howler monkeys by microscopy and PCR. Samples were included from individuals if they had a history of travel to or resided in areas within the Rio de Janeiro Atlantic Forest, but not if they had malaria prophylaxis, blood transfusion or tissue or organ transplantation, or had travelled to known malaria endemic areas in the preceding year. Additionally, we developed a molecular assay based on sequencing of the parasite mitochondrial genome to distinguish between P vivax and P simium, and applied this assay to 33 cases from outbreaks that occurred in 2015, and 2016.

FINDINGS

A total of 49 autochthonous malaria cases were reported in 2015-16. Most patients were male, with a mean age of 44 years (SD 14·6), and 82% lived in urban areas of Rio de Janeiro state and had visited the Atlantic Forest for leisure or work-related activities. 33 cases were used for mitochondrial DNA sequencing. The assay was successfully performed for 28 samples, and all were shown to be P simium, indicative of zoonotic transmission of this species to human beings in this region. Sequencing of the whole mitochondrial genome of three of these cases showed that P simium is most closely related to P vivax parasites from South America. The malaria outbreaks in this region were caused by P simium, previously considered to be a monkey-specific malaria parasite, related to but distinct from P vivax, and which has never conclusively been shown to infect people before.

INTERPRETATION

This unequivocal demonstration of zoonotic transmission, 50 years after the only previous report of P simium in people, leads to the possibility that this parasite has always infected people in this region, but that it has been consistently misdiagnosed as P vivax because of an absence of molecular typing techniques. Thorough screening of local non-human primates and mosquitoes (Anopheline) is required to evaluate the extent of this newly recognised zoonotic threat to public health and malaria elimination in Brazil.

FUNDING

Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado de Rio de Janeiro, The Brazilian National Council for Scientific and Technological Development (CNPq), JSPS Grant-in-Aid for scientific research, Secretary for Health Surveillance of the Brazilian Ministry of Health, Global Fund, Fundaçao de amparo à pesquisa do estado de Minas Gerais (Fapemig), and PRONEX Program of the CNPq.

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.

Comparative sequence analysis of Cyclospora cayetanensis apicoplast genomes originating from diverse geographical regions

BACKGROUND

Cyclospora cayetanensis is an emerging coccidian parasite that causes endemic and epidemic diarrheal disease called cyclosporiasis, and this infection is associated with consumption of contaminated produce or water in developed and developing regions. Food-borne outbreaks of cyclosporiasis have occurred almost every year in the USA since the 1990s. Investigations of these outbreaks are currently hampered due to lack of molecular epidemiological tools for trace back analysis. The apicoplast of C. cayetanensis, a relict non-photosynthetic plastid with an independent genome, provides an attractive target to discover sequence polymorphisms useful as genetic markers for detection and trace back analysis of the parasite. Distinct differences in the apicoplast genomes of C. cayetanensis could be useful in designing advanced molecular methods for rapid detection and, subtyping and geographical source attribution, which would aid outbreak investigations and surveillance studies.

METHODS

To obtain the genome sequence of the C. cayetanensis apicoplast, we sequenced the C. cayetanensis genomic DNA extracted from clinical stool samples, assembled and annotated a 34,146 bp-long circular sequence, and used this sequence as a reference genome in this study. We compared the genome and the predicted proteome to the data available from other apicomplexan parasites. To initialize the search for genetic markers, we mapped the raw sequence reads from an additional 11 distinct clinical stool samples originating from Nepal, New York, Texas, and Indonesia to the apicoplast reference genome.

RESULTS

We identified several high quality single nucleotide polymorphisms (SNPs) and small insertion/deletions spanning the apicoplast genome supported by extensive sequencing reads data, and a 30 bp sequence repeat at the terminal spacer region in a Nepalese sample. The predicted proteome consists of 29 core apicomplexan peptides found in most of the apicomplexans. Cluster analysis of these C. cayetanensis apicoplast genomes revealed a familiar pattern of tight grouping with Eimeria and Toxoplasma, separated from distant species such as Plasmodium and Babesia.

CONCLUSIONS

SNPs and sequence repeats identified in this study may be useful as genetic markers for identification and differentiation of C. cayetanensis isolates found and could facilitate outbreak investigations.

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

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

Further Evidence of Increasing Diversity of Plasmodium vivax in the Republic of Korea in Recent Years

BACKGROUND

Vivax malaria was successfully eliminated from the Republic of Korea (ROK) in the late 1970s but re-emerged in 1993. Two decades later as the ROK enters the final stages of malaria elimination, dedicated surveillance of the local P. vivax population is critical. We apply a population genetic approach to gauge P. vivax transmission dynamics in the ROK between 2010 and 2012.

METHODOLOGY/PRINCIPAL FINDINGS

P. vivax positive blood samples from 98 autochthonous cases were collected from patients attending health centers in the ROK in 2010 (n = 27), 2011 (n = 48) and 2012 (n = 23). Parasite genotyping was undertaken at 9 tandem repeat markers. Although not reaching significance, a trend of increasing population diversity was observed from 2010 (HE = 0.50 ± 0.11) to 2011 (HE = 0.56 ± 0.08) and 2012 (HE = 0.60 ± 0.06). Conversely, linkage disequilibrium declined during the same period: IAS = 0.15 in 2010 (P = 0.010), 0.09 in 2011 (P = 0.010) and 0.05 in 2012 (P = 0.010). In combination with data from other ROK studies undertaken between 1994 and 2007, our results are consistent with increasing parasite divergence since re-emergence. Polyclonal infections were rare (3% infections) suggesting that local out-crossing alone was unlikely to explain the increased divergence. Cases introduced from an external reservoir may therefore have contributed to the increased diversity. Aside from one isolate, all infections carried a short MS20 allele (142 or 149 bp), not observed in other studies in tropical endemic countries despite high diversity, inferring that these regions are unlikely reservoirs.

CONCLUSIONS

Whilst a number of factors may explain the observed population genetic trends, the available evidence suggests that an external geographic reservoir with moderate diversity sustains the majority of P. vivax infection in the ROK, with important implications for malaria elimination.

Targeting vivax malaria in the Asia Pacific: The Asia Pacific Malaria Elimination Network Vivax Working Group

The Asia Pacific Malaria Elimination Network (APMEN) is a collaboration of 18 country partners committed to eliminating malaria from within their borders. Over the past 5 years, APMEN has helped to build the knowledge, tools and in-country technical expertise required to attain this goal. At its inaugural meeting in Brisbane in 2009, Plasmodium vivax infections were identified across the region as a common threat to this ambitious programme; the APMEN Vivax Working Group was established to tackle specifically this issue. The Working Group developed a four-stage strategy to identify knowledge gaps, build regional consensus on shared priorities, generate evidence and change practice to optimize malaria elimination activities. This case study describes the issues faced and the solutions found in developing this robust strategic partnership between national programmes and research partners within the Working Group. The success of the approach adopted by the group may facilitate similar applications in other regions seeking to deploy evidence-based policy and practice.

Genetic similarities between Cyclospora cayetanensis and cecum-infecting avian Eimeria spp. in apicoplast and mitochondrial genomes

Background

Cyclospora cayetanensis is an important cause for diarrhea in children in developing countries and foodborne outbreaks of cyclosporiasis in industrialized nations. To improve understanding of the basic biology of Cyclospora spp. and development of molecular diagnostic tools and therapeutics, we sequenced the complete apicoplast and mitochondrial genomes of C. cayetanensis.

Methods

The genome of one Chinese C. cayetanensis isolate was sequenced using Roche 454 and Illumina technologies. The assembled genomes of the apicoplast and mitochondrion were retrieved, annotated, and compared with reference genomes for other apicomplexans to infer genome organizations and phylogenetic relationships. Sequence variations in the mitochondrial genome were identified by comparison of two C. cayetanensis nucleotide sequences from this study and a recent publication.

Results

The apicoplast and mitochondrial genomes of C. cayetanensis are 34,155 and 6,229 bp in size and code for 65 and 5 genes, respectively. Comparative genomic analysis showed high similarities between C. cayetanensis and Eimeria tenella in both genomes; they have 85.6 % and 90.4 % nucleotide sequence similarities, respectively, and complete synteny in gene organization. Phylogenetic analysis of the genomic sequences confirmed the genetic similarities between cecum-infecting avian Eimeria spp. and C. cayetanensis. Like in other coccidia, both genomes of C. cayetanensis are transcribed bi-directionally. The apicoplast genome is circular, codes for the complete machinery for protein biosynthesis, and contains two inverted repeats that differ slightly in LSU rRNA gene sequences. In contrast, the mitochondrial genome has a linear concatemer or circular mapping topology. Eight single-nucleotide and one 7-bp multiple-nucleotide variants were detected between the mitochondrial genomes of C. cayetanensis from this and recent studies.

Conclusions

The apicoplast and mitochondrial genomes of C. cayetanensis are highly similar to those of cecum-infecting avian Eimeria spp. in both genome organization and sequences. The availability of sequence data beyond rRNA and heat shock protein genes could facilitate studies of C. cayetanensis biology and development of genotyping tools for investigations of cyclosporiasis outbreaks.

Plasmodium vivax Diversity and Population Structure across Four Continents

Plasmodium vivax is the geographically most widespread human malaria parasite. To analyze patterns of microsatellite diversity and population structure across countries of different transmission intensity, genotyping data from 11 microsatellite markers was either generated or compiled from 841 isolates from four continents collected in 1999–2008. Diversity was highest in South-East Asia (mean allelic richness 10.0–12.8), intermediate in the South Pacific (8.1–9.9) Madagascar and Sudan (7.9–8.4), and lowest in South America and Central Asia (5.5–7.2). A reduced panel of only 3 markers was sufficient to identify approx. 90% of all haplotypes in South Pacific, African and SE-Asian populations, but only 60–80% in Latin American populations, suggesting that typing of 2–6 markers, depending on the level of endemicity, is sufficient for epidemiological studies. Clustering analysis showed distinct clusters in Peru and Brazil, but little sub-structuring was observed within Africa, SE-Asia or the South Pacific. Isolates from Uzbekistan were exceptional, as a near-clonal parasite population was observed that was clearly separated from all other populations (F_ST>0.2). Outside Central Asia F_ST values were highest (0.11–0.16) between South American and all other populations, and lowest (0.04–0.07) between populations from South-East Asia and the South Pacific. These comparisons between P. vivax populations from four continents indicated that not only transmission intensity, but also geographical isolation affect diversity and population structure. However, the high effective population size results in slow changes of these parameters. This persistency must be taken into account when assessing the impact of control programs on the genetic structure of parasite populations.

Plasmodium vivax is the predominant malaria parasite in Latin America, Asia and the South Pacific. Different factors are expected to shape diversity and population structure across continents, e.g. transmission intensity which is much lower in South America as compared to Southeast-Asia and the South Pacific, or geographical isolation of P. vivax populations in the South Pacific. We have compiled data from 841 isolates from South and Central America, Africa, Central Asia, Southeast-Asia and the South Pacific typed with a panel of 11 microsatellite markers. Diversity was highest in Southeast-Asia, where transmission is intermediate-high and migration of infected hosts is high, and lowest in South America and Central Asia where malaria transmission is low and focal. Reducing the panel of microsatellites showed that 2–6 markers are sufficient for genotyping for most drug trials and epidemiological studies, as these markers can identify >90% of all haplotypes. Parasites clustered according to continental origin, with high population differentiation between South American and Central Asian populations and the other populations, and lowest differences between Southeast-Asia and the South Pacific. Current attempts to reduce malaria transmission might change this pattern, but only after transmission is reduced for an extended period of time.

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.

Tracking malaria parasites in the eradication era

As more endemic countries enter the elimination phase, the detection of sporadic malaria infections or outbreaks elicits many questions: (i) are the infections locally acquired or imported? (ii) If imported, where do they come from? (iii) Do outbreak strains have a single or multiple geographic origins? New molecular barcoding methods provide ways to analyze clinical malaria parasite samples and answer these and other crucial questions.

Genetic diversity and population structure of Plasmodium vivax in Central China

BACKGROUND

In Central China the declining incidence of Plasmodium vivax has been interrupted by epidemic expansions and imported cases. The impact of these changes on the local parasite population, and concurrent risks of future resurgence, was assessed.

METHODS

Plasmodium vivax isolates collected from Anhui and Jiangsu provinces, Central China between 2007 and 2010 were genotyped using capillary electrophoresis at seven polymorphic short tandem repeat markers. Spatial and temporal analyses of within-host and population diversity, population structure, and relatedness were conducted on these isolates.

RESULTS

Polyclonal infections were infrequent in the 94 isolates from Anhui (4%) and 25 from Jiangsu (12%), with a trend for increasing frequency from 2008 to 2010 (2 to 19%) when combined. Population diversity was high in both provinces and across the years tested (H(E) = 0.8 - 0.85). Differentiation between Anhui and Jiangsu was modest (F'(ST) = 0.1). Several clusters of isolates with identical multi-locus haplotypes were observed across both Anhui and Jiangsu. Linkage disequilibrium was strong in both populations and in each year tested (I(A)(S) = 0.2 - 0.4), but declined two- to four-fold when identical haplotypes were accounted for, indicative of occasional epidemic transmission dynamics. None of five imported isolates shared identical haplotypes to any of the central Chinese isolates.

CONCLUSIONS

The population genetic structure of P. vivax in Central China highlights unstable transmission, with limited barriers to gene flow between the central provinces. Despite low endemicity, population diversity remained high, but the reservoirs sustaining this diversity remain unclear. The challenge of imported cases and risks of resurgence emphasize the need for continued surveillance to detect early warning signals. Although parasite genotyping has potential to inform the management of outbreaks, further studies are required to identify suitable marker panels for resolving local from imported P. vivax isolates.