Plasmodium species causing malaria in humans are not monophyletic, sharing common ancestors with nonhuman primate parasites. Plasmodium gonderi is one of the few known Plasmodium species infecting African old-world monkeys that are not found in apes. This study reports a de novo assembled P. gonderi genome with complete chromosomes. The P. gonderi genome shares codon usage, syntenic blocks, and other characteristics with the human parasites Plasmodium ovale s.l. and Plasmodium malariae, also of African origin, and the human parasite Plasmodium vivax and species found in nonhuman primates from Southeast Asia. Using phylogenetically aware methods, newly identified syntenic blocks were found enriched with conserved metabolic genes. Regions outside those blocks harbored genes encoding proteins involved in the vertebrate host-Plasmodium relationship undergoing faster evolution. Such genome architecture may have facilitated colonizing vertebrate hosts. Phylogenomic analyses estimated the common ancestor between P. vivax and an African ape parasite P. vivax-like, within the Asian nonhuman primates parasites clade. Time estimates incorporating P. gonderi placed the P. vivax and P. vivax-like common ancestor in the late Pleistocene, a time of active migration of hominids between Africa and Asia. Thus, phylogenomic and time-tree analyses are consistent with an Asian origin for P. vivax and an introduction of P. vivax-like into Africa. Unlike other studies, time estimates for the clade with Plasmodium falciparum, the most lethal human malaria parasite, coincide with their host species radiation, African hominids. Overall, the newly assembled genome presented here has the quality to support comparative genomic investigations in Plasmodium.

Multiplex bead assays (MBA) that measure IgG antibodies to the carboxy-terminal 19-kDa sub-unit of the merozoite surface protein 1 (MSP119) are currently used to determine malaria seroprevalence in human populations living in areas with both stable and unstable transmission. However, the species specificities of the IgG antibody responses to the malaria MSP119 antigens have not been extensively characterized using MBA.

The majority of malaria rapid diagnostic tests (RDTs) detect Plasmodium falciparum histidine-rich protein 2 (PfHRP2), encoded by the pfhrp2 gene. Recently, P. falciparum isolates from Peru were found to lack pfhrp2 leading to false-negative RDT results. We hypothesized that pfhrp2-deleted parasites in Peru derived from a single genetic event. We evaluated the parasite population structure and pfhrp2 haplotype of samples collected between 1998 and 2005 using seven neutral and seven chromosome 8 microsatellite markers, respectively. Five distinct pfhrp2 haplotypes, corresponding to five neutral microsatellite-based clonal lineages, were detected in 1998-2001; pfhrp2 deletions occurred within four haplotypes. In 2003-2005, outcrossing among the parasite lineages resulted in eight population clusters that inherited the five pfhrp2 haplotypes seen previously and a new haplotype; pfhrp2 deletions occurred within four of these haplotypes. These findings indicate that the genetic origin of pfhrp2 deletion in Peru was not a single event, but likely occurred multiple times.

Recent studies have demonstrated the deletion of the histidine-rich protein 2 (PfHRP2) gene (pfhrp2) in field isolates of Plasmodium falciparum, which could result in false negative test results when PfHRP2-based rapid diagnostic tests (RDTs) are used for malaria diagnosis. Although primary diagnosis of malaria in Honduras is determined based on microscopy, RDTs may be useful in remote areas. In this study, it was investigated whether there are deletions of the pfhrp2, pfhrp3 and their respective flanking genes in 68 P. falciparum parasite isolates collected from the city of Puerto Lempira, Honduras. In addition, further investigation considered the possible correlation between parasite population structure and the distribution of these gene deletions by genotyping seven neutral microsatellites.

Plasmodium vivax is the most prevalent human malaria parasite outside of Africa. Yet, studies aimed to identify genes with signatures consistent with natural selection are rare. Here, we present a comparative analysis of the pattern of genetic variation of five sequenced isolates of P. vivax and its divergence with two closely related species, Plasmodium cynomolgi and Plasmodium knowlesi, using a set of orthologous genes. In contrast to Plasmodium falciparum, the parasite that causes the most lethal form of human malaria, we did not find significant constraints on the evolution of synonymous sites genome wide in P. vivax. The comparative analysis of polymorphism and divergence across loci allowed us to identify 87 genes with patterns consistent with positive selection, including genes involved in the "exportome" of P. vivax, which are potentially involved in evasion of the host immune system. Nevertheless, we have found a pattern of polymorphism genome wide that is consistent with a significant amount of constraint on the replacement changes and prevalent negative selection. Our analyses also show that silent polymorphism tends to be larger toward the ends of the chromosomes, where many genes involved in antigenicity are located, suggesting that natural selection acts not only by shaping the patterns of variation within the genes but it also affects genome organization.

Detection of histidine-rich protein 2 (HRP2) from the malaria parasite Plasmodium falciparum provides evidence for active or recent infection, and is utilized for both diagnostic and surveillance purposes, but current laboratory immunoassays for HRP2 are hindered by low sensitivities and high costs. Here we present a new HRP2 immunoassay based on antigen capture through a bead-based system capable of detecting HRP2 at sub-picogram levels. The assay is highly specific and cost-effective, allowing fast processing and screening of large numbers of samples. We utilized the assay to assess results of HRP2-based rapid diagnostic tests (RDTs) in different P. falciparum transmission settings, generating estimates for true performance in the field. Through this method of external validation, HRP2 RDTs were found to perform well in the high-endemic areas of Mozambique and Angola with 86.4% and 73.9% of persons with HRP2 in their blood testing positive by RDTs, respectively, and false-positive rates of 4.3% and 0.5%. However, in the low-endemic setting of Haiti, only 14.5% of persons found to be HRP2 positive by the bead assay were RDT positive. Additionally, 62.5% of Haitians showing a positive RDT test had no detectable HRP2 by the bead assay, likely indicating that these were false positive tests. In addition to RDT validation, HRP2 biomass was assessed for the populations in these different settings, and may provide an additional metric by which to estimate P. falciparum transmission intensity and measure the impact of interventions.

Social network structure is a fundamental determinant of human health, from infectious to chronic diseases. However, quantitative and unbiased approaches to measuring social network structure are lacking. We hypothesized that genetic relatedness of oral commensal bacteria could be used to infer social contact between humans, just as genetic relatedness of pathogens can be used to determine transmission chains of pathogens. We used a traditional, questionnaire survey-based method to characterize the contact network of the School of Public Health at a large research university. We then collected saliva from a subset of individuals to analyze their oral microflora using a modified deep sequencing multilocus sequence typing (MLST) procedure. We examined micro-evolutionary changes in the S. viridans group to uncover transmission patterns reflecting social network structure. We amplified seven housekeeping gene loci from the Streptococcus viridans group, a group of ubiquitous commensal bacteria, and sequenced the PCR products using next-generation sequencing. By comparing the generated S. viridans reads between pairs of individuals, we reconstructed the social network of the sampled individuals and compared it to the network derived from the questionnaire survey-based method. The genetic relatedness significantly (p-value < 0.001) correlated with social distance in the questionnaire-based network, and the reconstructed network closely matched the network derived from the questionnaire survey-based method. Oral commensal bacterial are thus likely transmitted through routine physical contact or shared environment. Their genetic relatedness can be used to represent a combination of social contact and shared physical space, therefore reconstructing networks of contact. This study provides the first step in developing a method to measure direct social contact based on commensal organism genotyping, potentially capable of unmasking hidden social networks that contribute to pathogen transmission.

Recent anti-malarial resistance monitoring in Angola has shown efficacy of artemether-lumefantrine (AL) in certain sites approaching the key 90% lower limit of efficacy recommended for artemisinin-based combination therapy. In addition, a controversial case of malaria unresponsive to artemisinins was reported in a patient infected in Lunda Sul Province in 2013.

Mutations in the dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) genes of Plasmodium falciparum are associated with resistance to anti-folate drugs, most notably sulphadoxine-pyrimethamine (SP). Molecular studies document the prevalence of these mutations in parasite populations across the African continent. However, there is no systematic review examining the collective epidemiological significance of these studies. This meta-analysis attempts to: 1) summarize genotype frequency data that are critical for molecular surveillance of anti-folate resistance and 2) identify the specific challenges facing the development of future molecular databases.

RTS,S is the most advanced malaria vaccine candidate, currently under phase-III clinical trials in Africa. This Plasmodium falciparum vaccine contains part of the central repeat region and the complete C-terminal T cell epitope region (Th2R and Th3R) of the circumsporozoite protein (CSP). Since naturally occurring polymorphisms at the vaccine candidate loci are critical determinants of the protective efficacy of the vaccines, it is imperative to investigate these polymorphisms in field isolates. In this study we have investigated the genetic diversity at the central repeat, C-terminal T cell epitope (Th2R and Th3R) and N-terminal T cell epitope regions of the CSP, in P. falciparum isolates from Madhya Pradesh state of India. These isolates were collected through a 5-year prospective study aimed to develop a well-characterized field-site for the future evaluation of malaria vaccine in India. Our results revealed that the central repeat (63 haplotypes, n = 161) and C-terminal Th2R/Th3R epitope (24 haplotypes, n = 179) regions were highly polymorphic, whereas N-terminal non-repeat region was less polymorphic (5 haplotypes, n = 161) in this population. We did not find any evidence of the role of positive natural selection in maintaining the genetic diversity at the Th2R/Th3R regions of CSP. Comparative analysis of the Th2R/Th3R sequences from this study to the global isolates (n = 1160) retrieved from the GenBank database revealed two important points. First, the majority of the sequences (~61%, n = 179) from this study were identical to the Dd2/Indochina type, which is also the predominant Th2R/Th3R haplotype in Asia (~59%, n = 974). Second, the Th2R/Th3R sequences in Asia, South America and Africa are geographically distinct with little allele sharing between continents. In conclusion, this study provides an insight on the existing polymorphisms in the CSP in a parasite population from India that could potentially influence the efficacy of RTS,S vaccine in this region.

Molecular methods elucidate evolutionary and ecological processes in parasites, where interaction between hosts and parasites enlighten the evolution of parasite lifestyles and host defenses. Population genetics of Plasmodium vivax parasites accurately describe transmission dynamics of the parasites and evaluation of malaria control measures. As a first generation vaccine candidate against malaria, the Circumsporozoite Protein (CSP) has demonstrated significant potential in P. falciparum. Extensive polymorphism hinders the development of a potent malaria vaccine. Hence, the genetic diversity of Pvcsp was investigated for the first time in 60 Sri Lankan clinical isolates by obtaining the nucleotide sequence of the central repeat (CR) domain and examining the polymorphism of the peptide repeat motifs (PRMs), the genetic diversity indices and phylogenetic relationships. PCR amplicons determined size polymorphism of 610, 700 and 710 bp in Pvcsp of Sri Lanka where all amino acid sequences obtained were of the VK210 variant, consisting variable repeats of 4 different PRMs. The two most abundant PRMs of the CR domain, GDRADGQPA and GDRAAGQPA consisted ~2-4 repeats, while GNRAAGQPA was unique to the island. Though, different nucleotide sequences termed repeat allotypes (RATs) were observed for each PRM, these were synonymous contributing to a less polymorphic CR domain. The genetic diversity of Pvcsp in Sri Lanka was due to the number of repetitive peptide repeat motifs, point mutations, and intragenic recombination. The 19 amino acid haplotypes defined were exclusive to Sri Lanka, whereas the 194 Pvcsp sequences of global isolates generated 57 more distinct a.a. haplotypes of the VK210 variant. Strikingly, the CR domain of both VK210 and VK247 variants was under purifying selection interpreting the scarcity of CSP non-synonymous polymorphisms. Insights to the distribution of RATs in the CR region with geographic clustering of the P. vivax VK210 variant were revealed. The cladogram reiterated this unique geographic clustering of local (VK210) and global isolates (VK210 and VK247), which was further validated by the elevated fixation index values of the VK210 variant.

Understanding the spatial clustering of Plasmodium falciparum populations can assist efforts to contain drug-resistant parasites and maintain the efficacy of future drugs. We sequenced single nucleotide polymorphisms (SNPs) in the dihydropteroate synthase gene (dhps) associated with sulfadoxine resistance and 5 microsatellite loci flanking dhps in order to investigate the genetic backgrounds, genetic relatedness, and geographic clustering of falciparum parasites in the Democratic Republic of the Congo (DRC). Resistant haplotypes were clustered into subpopulations: one in the northeast DRC, and the other in the balance of the DRC. Network and clonal lineage analyses of the flanking microsatellites indicate that geographically-distinct mutant dhps haplotypes derive from separate lineages. The DRC is therefore a watershed for haplotypes associated with sulfadoxine resistance. Given the importance of central Africa as a corridor for the spread of antimalarial resistance, the identification of the mechanisms of this transit can inform future policies to contain drug-resistant parasite strains.

Determining the source of malaria outbreaks in Ecuador and identifying remaining transmission foci will help in malaria elimination efforts. In this study, the genetic signatures of Plasmodium falciparum isolates, obtained from an outbreak that occurred in northwest Ecuador from 2012 to 2013, were characterized.

In vivo efficacy assessments of antimalarials are essential for ensuring effective case management. In Ethiopia, chloroquine (CQ) without primaquine is the first-line treatment for Plasmodium vivax in malarious areas, but artemether-lumefantrine (AL) is also commonly used.

False-negative results for Plasmodium falciparum histidine-rich protein (HRP) 2-based rapid diagnostic tests (RDTs) are increasing in Eritrea. We investigated HRP gene 2/3 (pfhrp2/pfhrp3) status in 50 infected patients at 2 hospitals. We showed that 80.8% (21/26) of patients at Ghindae Hospital and 41.7% (10/24) at Massawa Hospital were infected with pfhrp2-negative parasites and 92.3% (24/26) of patients at Ghindae Hospital and 70.8% (17/24) at Massawa Hospital were infected with pfhrp3-negative parasites. Parasite densities between pfhrp2-positive and pfhrp2-negative patients were comparable. All pfhrp2-negative samples had no detectable HRP2/3 antigen and showed negative results for HRP2-based RDTs. pfhrp2-negative parasites were genetically less diverse and formed 2 clusters with no close relationships to parasites from Peru. These parasites probably emerged independently by selection in Eritrea. High prevalence of pfhrp2-negative parasites caused a high rate of false-negative results for RDTs. Determining prevalence of pfhrp2-negative parasites is urgently needed in neighboring countries to assist case management policies.

Universal coverage with long-lasting insecticidal nets (LLINs) is a primary control strategy against Plasmodium falciparum malaria. However, its impact on the three other main species of human malaria and lymphatic filariasis (LF), which share the same vectors in many co-endemic areas, is not as well characterized. The recent development of multiplex antibody detection provides the opportunity for simultaneous evaluation of the impact of control measures on the burden of multiple diseases.

Haemosporidians are a diverse group of vector-borne parasitic protozoa that includes the agents of human malaria; however, most of the described species are found in birds and reptiles. Although our understanding of these parasites' diversity has expanded by analyses of their mitochondrial genes, there is limited information on these genes' evolutionary rates. Here, 114 mitochondrial genomes (mtDNA) were studied from species belonging to four genera: Leucocytozoon, Haemoproteus, Hepatocystis, and Plasmodium. Contrary to previous assertions, the mtDNA is phylogenetically informative. The inferred phylogeny showed that, like the genus Plasmodium, the Leucocytozoon and Haemoproteus genera are not monophyletic groups. Although sensitive to the assumptions of the molecular dating method used, the estimated times indicate that the diversification of the avian haemosporidian subgenera/genera took place after the Cretaceous-Paleogene boundary following the radiation of modern birds. Furthermore, parasite clade differences in mtDNA substitution rates and strength of negative selection were detected. These differences may affect the biological interpretation of mtDNA gene lineages used as a proxy to species in ecological and parasitological investigations. Given that the mitochondria are critically important in the parasite life cycle stages that take place in the vector and that the transmission of parasites belonging to particular clades has been linked to specific insect families/subfamilies, this study suggests that differences in vectors have affected the mode of evolution of haemosporidian mtDNA genes. The observed patterns also suggest that the radiation of haemosporidian parasites may be the result of community-level evolutionary processes between their vertebrate and invertebrate hosts.

The Thailand-Cambodia border is the epicenter for drug-resistant falciparum malaria. Previous studies have shown that chloroquine (CQ) and pyrimethamine resistance originated in this region and eventually spread to other Asian countries and Africa. However, there is a dearth in understanding the origin and evolution of dhps alleles associated with sulfadoxine resistance. The present study was designed to reveal the origin(s) of sulfadoxine resistance in Cambodia and its evolutionary relationship to African and South American dhps alleles. We sequenced 234 Cambodian Plasmodium falciparum isolates for the dhps codons S436A/F, A437G, K540E, A581G and A613S/T implicated in sulfadoxine resistance. We also genotyped 10 microsatellite loci around dhps to determine the genetic backgrounds of various alleles and compared them with the backgrounds of alleles prevalent in Africa and South America. In addition to previously known highly-resistant triple mutant dhps alleles SGEGA and AGEAA (codons 436, 437, 540, 581, 613 are sequentially indicated), a large proportion of the isolates (19.3%) contained a 540N mutation in association with 437G/581G yielding a previously unreported triple mutant allele, SGNGA. Microsatellite data strongly suggest the strength of selection was greater on triple mutant dhps alleles followed by the double and single mutants. We provide evidence for at least three independent origins for the double mutants, one each for the SGKGA, AGKAA and SGEAA alleles. Our data suggest that the triple mutant allele SGEGA and the novel allele SGNGA have common origin on the SGKGA background, whereas the AGEAA triple mutant was derived from AGKAA on multiple, albeit limited, genetic backgrounds. The SGEAA did not share haplotypes with any of the triple mutants. Comparative analysis of the microsatellite haplotypes flanking dhps alleles from Cambodia, Kenya, Cameroon and Venezuela revealed an independent origin of sulfadoxine resistant alleles in each of these regions.

There is an increased interest in potential zoonotic malarias. To date, Plasmodium malariae that infects humans remains indistinguishable from Plasmodium brasilianum, which is widespread among New World primates. Distributed throughout tropical Central and South America, the Callitrichidae are small arboreal primates in which detection of natural Plasmodium infection has been extremely rare. Most prior screening efforts have been limited to small samples, the use of low-probability detection methods, or both. Rarely have screening efforts implemented a longitudinal sampling design. Through an annual mark-recapture program of two sympatric callitrichids, the emperor (Saguinus imperator) and saddleback (Saguinus fuscicollis) tamarins, whole blood samples were screened for Plasmodium by microscopy and nested PCR of the cytochrome b gene across four consecutive years (2012-2015). Following the first field season, approximately 50% of the samples collected each subsequent year were from recaptured individuals. In particular, out of 245 samples from 129 individuals, 11 samples from 6 individuals were positive for Plasmodium, and all but one of these infections was found in S. imperator. Importantly, the cytochrome b sequences were 100% identical to former isolates of P. malariae from humans and P. brasilianum from Saimiri sp. Chronic infections were detected as evidenced by repeated infections (7) from two individuals across the 4-year study period. Furthermore, 4 of the 5 infected emperor tamarins were part of a single group spanning the entire study period. Overall, the low prevalence reported here is consistent with previous findings. This study identifies two new natural hosts for P. brasilianum and provides evidence in support of chronic infections in wildlife populations. Given that callitrichids are often found in mixed-species associations with other primates and can be resilient to human-disturbed environments, they could contribute to the maintenance of P. malariae populations if future work provides entomological and epidemiological evidence indicating human zoonotic infections.

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.