Ananias A. Escalante

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.

Malaria Molecular Epidemiology: Lessons from the International Centers of Excellence for Malaria Research Network

Molecular epidemiology leverages genetic information to study the risk factors that affect the frequency and distribution of malaria cases. This article describes molecular epidemiologic investigations currently being carried out by the International Centers of Excellence for Malaria Research (ICEMR) network in a variety of malaria-endemic settings. First, we discuss various novel approaches to understand malaria incidence and gametocytemia, focusing on Plasmodium falciparum and Plasmodium vivax. Second, we describe and compare different parasite genotyping methods commonly used in malaria epidemiology and population genetics. Finally, we discuss potential applications of molecular epidemiological tools and methods toward malaria control and elimination efforts.

Genome-Wide Patterns of Genetic Polymorphism and Signatures of Selection in Plasmodium vivax

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.

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.

Multiplicity of Infection and Disease Severity in Plasmodium vivax

Background Multiplicity of infection (MOI) refers to the average number of distinct parasite genotypes concurrently infecting a patient. Although several studies have reported on MOI and the frequency of multiclonal infections in Plasmodium falciparum, there is limited data on Plasmodium vivax. Here, MOI and the frequency of multiclonal infections were studied in areas from South America where P. vivax and P. falciparum can be compared. Methodology/Principal Findings As part of a passive surveillance study, 1,328 positive malaria patients were recruited between 2011 and 2013 in low transmission areas from Colombia. Of those, there were only 38 P. vivax and 24 P. falciparum clinically complicated cases scattered throughout the time of the study. Samples from uncomplicated cases were matched in time and location with the complicated cases in order to compare the circulating genotypes for these two categories. A total of 92 P. vivax and 57 P. falciparum uncomplicated cases were randomly subsampled. All samples were genotyped by using neutral microsatellites. Plasmodium vivax showed more multiclonal infections (47.7%) than P. falciparum (14.8%). Population genetics and haplotype network analyses did not detect differences in the circulating genotypes between complicated and uncomplicated cases in each parasite. However, a Fisher exact test yielded a significant association between having multiclonal P. vivax infections and complicated malaria. No association was found for P. falciparum infections. Conclusion The association between multiclonal infections and disease severity in P. vivax is consistent with previous observations made in rodent malaria. The contrasting pattern between P. vivax and P. falciparum could be explained, at least in part, by the fact that P. vivax infections have lineages that were more distantly related among them than in the case of the P. falciparum multiclonal infections. Future research should address the possible role that acquired immunity and exposure may have on multiclonal infections and their association with disease severity.

Mode and rate of evolution of haemosporidian mitochondrial genomes: Timing the radiation of avian parasites

&NA; 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.

Population Genetics, Evolutionary Genomics, and Genome-Wide Studies of Malaria: A View Across the International Centers of Excellence for Malaria Research.

The study of the three protagonists in malaria—the Plasmodium parasite, the Anopheles mosquito, and the human host—is key to developing methods to control and eventually eliminate the disease. Genomic technologies, including the recent development of next-generation sequencing, enable interrogation of this triangle to an unprecedented level of scrutiny, and promise exciting progress toward real-time epidemiology studies and the study of evolutionary adaptation. We discuss the use of genomics by the International Centers of Excellence for Malaria Research, a network of field sites and laboratories in malaria-endemic countries that undertake cutting-edge research, training, and technology transfer in malarious countries of the world.

Profiles of low complexity regions in Apicomplexa

BackgroundLow complexity regions (LCRs) are a ubiquitous feature in genomes and yet their evolutionary history and functional roles are unclear. Previous studies have shown contrasting evidence in favor of both neutral and selective mechanisms of evolution for different sets of LCRs suggesting that modes of identification of these regions may play a role in our ability to discern their evolutionary history. To further investigate this issue, we used a multiple threshold approach to identify species-specific profiles of proteome complexity and, by comparing properties of these sets, determine the influence that starting parameters have on evolutionary inferences.ResultsWe find that, although qualitatively similar, quantitatively each species has a unique LCR profile which represents the frequency of these regions within each genome. Inferences based on these profiles are more accurate in comparative analyses of genome complexity as they allow to determine the relative complexity of multiple genomes as well as the type of repetitiveness that is most common in each. Based on the multiple threshold LCR sets obtained, we identified predominant evolutionary mechanisms at different complexity levels, which show neutral mechanisms acting on highly repetitive LCRs (e.g., homopolymers) and selective forces becoming more important as heterogeneity of the LCRs increases.ConclusionsOur results show how inferences based on LCRs are influenced by the parameters used to identify these regions. Sets of LCRs are heterogeneous aggregates of regions that include homo- and heteropolymers and, as such, evolve according to different mechanisms. LCR profiles provide a new way to investigate genome complexity across species and to determine the driving mechanism of their evolution.

Molecular Evolution of PvMSP3α Block II in Plasmodium vivax from Diverse Geographic Origins

Block II of Plasmodium vivax merozoite surface protein 3α (PvMSP3α) is conserved and has been proposed as a potential candidate for a malaria vaccine. The present study aimed to compare sequence diversity in PvMSP3a block II at a local microgeographic scale in a village as well as from larger geographic regions (countries and worldwide). Blood samples were collected from asymptomatic carriers of P. vivax in a village at the western border of Thailand and PvMSP3α was amplified and sequenced. For population genetic analysis, 237 PvMSP3α block II sequences from eleven P. vivax endemic countries were analyzed. PvMSP3α sequences from 20 village-level samples revealed two length variant types with one type containing a large deletion in block I. In contrast, block II was relatively conserved; especially, some non-synonymous mutations were extensively shared among 11 parasite populations. However, the majority of the low-frequency synonymous variations were population specific. The conserved pattern of nucleotide diversity in block II sequences was probably due to functional/structural constraints, which were further supported by the tests of neutrality. Notably, a small region in block II that encodes a predicted B cell epitope was highly polymorphic and showed signs of balancing selection, signifying that this region might be influenced by the immune selection and may serve as a starting point for designing multi-antigen/stage epitope based vaccines against this parasite.

Leucocytozoon Diversity and Possible Vectors in the Neotropical highlands of Colombia

Studies of the lowland avifauna in the Neotropical Region have shown a paucity of Leucocytozoon species. However, surveys conducted in the Colombian highlands revealed a great diversity of these parasites infecting resident birds. To further investigate the relationship between Leucocytozoon diversity, the potential vectors, and altitudinal distribution, birds from 41 families were sampled at low and high elevations in Colombia. Blood samples were screened by microscopy, and a fragment of cytochrome b was amplified from Leucocytozoon-positive samples. The complete mitochondrial genome was also obtained for each morphospecies of Leucocytozoon. Leucocytozoon species were detected in resident birds, with various degrees of host specificity, at elevations from 2,400 to 3,950 meters above sea level, where five new host-parasite associations were discovered. Phylogenetic analysis based on the cytochrome b fragment suggested that two nominal taxa, L. fringillinarum and L. majoris, are species complexes. Blood sources of Simuliidae revealed generalist-feeding habits that included avian and mammalian hosts. Molecular analysis of parasites in black flies indicated a close relationship with the parasites found in birds. Our investigation provides further evidence that the distribution and transmission of Leucocytozoon species in the Neotropics are influenced by elevation, with the highest prevalence between 2,400 and 3,200 m asl.