Early origin and recent expansion of Plasmodium falciparum

Evolutionary structure of Plasmodium falciparum major variant surface antigen genes in South America: Implications for epidemic transmission and surveillance

Strong founder effects resulting from human migration out of Africa have led to geographic variation in single nucleotide polymorphisms (SNPs) and microsatellites (MS) of the malaria parasite, Plasmodium falciparum. This is particularly striking in South America where two major founder populations of P. falciparum have been identified that are presumed to have arisen from the transatlantic slave trade. Given the importance of the major variant surface antigen of the blood stages of P. falciparum as both a virulence factor and target of immunity, we decided to investigate the population genetics of the genes encoding “Plasmodium falciparum Erythrocyte Membrane Protein 1” (PfEMP1) among several countries in South America, in order to evaluate the transmission patterns of malaria in this continent. Deep sequencing of the DBLα domain of var genes from 128 P. falciparum isolates from five locations in South America was completed using a 454 high throughput sequencing protocol. Striking geographic variation in var DBLα sequences, similar to that seen for SNPs and MS markers, was observed. Colombia and French Guiana had distinct var DBLα sequences, whereas Peru and Venezuela showed an admixture. The importance of such geographic variation to herd immunity and malaria vaccination is discussed.

Extreme mutation bias and high AT content in Plasmodium falciparum

For reasons that remain unknown, the Plasmodium falciparum genome has an exceptionally high AT content compared to other Plasmodium species and eukaryotes in general - nearly 80% in coding regions and approaching 90% in non-coding regions. Here, we examine how this phenomenon relates to genome-wide patterns of de novo mutation. Mutation accumulation experiments were performed by sequential cloning of six P. falciparum isolates growing in human erythrocytes in vitro for 4 years, with 279 clones sampled for whole genome sequencing at different time points. Genome sequence analysis of these samples revealed a significant excess of G:C to A:T transitions compared to other types of nucleotide substitution, which would naturally cause AT content to equilibrate close to the level seen across the P. falciparum reference genome (80.6% AT). These data also uncover an extremely high rate of small indel mutation relative to other species, primarily associated with repetitive AT-rich sequences, in addition to larger-scale structural rearrangements focused in antigen-coding var genes. In conclusion, high AT content in P. falciparum is driven by a systematic mutational bias and ultimately leads to an unusual level of microstructural plasticity, raising the question of whether this contributes to adaptive evolution.

Population Genomics of Paramecium Species

Population-genomic analyses are essential to understanding factors shaping genomic variation and lineage-specific sequence constraints. The dearth of such analyses for unicellular eukaryotes prompted us to assess genomic variation in Paramecium, one of the most well-studied ciliate genera. The Paramecium aurelia complex consists of ∼15 morphologically indistinguishable species that diverged subsequent to two rounds of whole-genome duplications (WGDs, as long as 320 MYA) and possess extremely streamlined genomes. We examine patterns of both nuclear and mitochondrial polymorphism, by sequencing whole genomes of 10-13 worldwide isolates of each of three species belonging to the P. aurelia complex: P. tetraurelia, P. biaurelia, P. sexaurelia, as well as two outgroup species that do not share the WGDs: P. caudatum and P. multimicronucleatum. An apparent absence of global geographic population structure suggests continuous or recent dispersal of Paramecium over long distances. Intergenic regions are highly constrained relative to coding sequences, especially in P. caudatum and P. multimicronucleatum that have shorter intergenic distances. Sequence diversity and divergence are reduced up to ∼100-150 bp both upstream and downstream of genes, suggesting strong constraints imposed by the presence of densely packed regulatory modules. In addition, comparison of sequence variation at non-synonymous and synonymous sites suggests similar recent selective pressures on paralogs within and orthologs across the deeply diverging species. This study presents the first genome-wide population-genomic analysis in ciliates and provides a valuable resource for future studies in evolutionary and functional genetics in Paramecium.

Experimental evidence for hybridization of closely related lineages in Plasmodium relictum

Over 50 avian Plasmodium species have been described. However, PCR-based information shows much broader diversity of genetic lineages in these parasites. This discrepancy indicates insufficient knowledge about taxonomic diversity and boundaries of a single species in avian Plasmodium species. In recent taxonomy, most of genetically closely related lineages that share the same morphology and development patterns are attributed to the same biological species, but there is no information if these lineages are able to cross. This information is crucial to understand if these lineages form single or multiple evolutionary units. Due to presence of sexual process and sporogonic development of Plasmodium parasites in mosquitoes, self and cross-fertilization can occur and be identified during the oocyst stage. We initiated in vivo hybridization experiments of two widespread Plasmodium relictum lineages (pSGS1 and pGRW11) in experimentally infected Culex pipiens pipiens form molestus mosquitoes. To study putative hybrid oocysts, we used a laser microdissection technique together with PCR-based analyses of mitochondrial and nuclear genes. We demonstrate that both pSGS1 and pGRW11 lineages develop in infected mosquitoes in parallel, but also form hybrid oocysts of these two lineages. Our results are in accord to a recent global phylogeographic study of P. relictum that suggested that cross-fertilization between pSGS1 and pGRW11 might occur. This information helps to understand population structure, gene flow and the evolutionary process of haemosporidian parasites.

Genome content analysis yields new insights into the relationship between the human malaria parasite Plasmodium falciparum and its anopheline vectors

Background

The persistent and growing gap between the availability of sequenced genomes and the ability to assign functions to sequenced genes led us to explore ways to maximize the information content of automated annotation for studies of anopheline mosquitos. Specifically, we use genome content analysis of a large number of previously sequenced anopheline mosquitos to follow the loss and gain of protein families over the evolutionary history of this group.

The importance of this endeavor lies in the potential for comparative genomic studies between Anopheles and closely related non-vector species to reveal ancestral genome content dynamics involved in vector competence. In addition, comparisons within Anopheles could identify genome content changes responsible for variation in the vectorial capacity of this family of important parasite vectors.

Results

The competence and capacity of P. falciparum vectors do not appear to be phylogenetically constrained within the Anophelinae. Instead, using ancestral reconstruction methods, we suggest that a previously unexamined component of vector biology, anopheline nucleotide metabolism, may contribute to the unique status of anophelines as P. falciparum vectors. While the fitness effects of nucleotide co-option by P. falciparum parasites on their anopheline hosts are not yet known, our results suggest that anopheline genome content may be responding to selection pressure from P. falciparum. Whether this response is defensive, in an attempt to redress improper nucleotide balance resulting from P. falciparum infection, or perhaps symbiotic, resulting from an as-yet-unknown mutualism between anophelines and P. falciparum, is an open question that deserves further study.

Conclusions

Clearly, there is a wealth of functional information to be gained from detailed manual genome annotation, yet the rapid increase in the number of available sequences means that most researchers will not have the time or resources to manually annotate all the sequence data they generate. We believe that efforts to maximize the amount of information obtained from automated annotation can help address the functional annotation deficit that most evolutionary biologists now face, and here demonstrate the value of such an approach.

Electronic supplementary material

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

Gancidin W, a potential low-toxicity antimalarial agent isolated from an endophytic Streptomyces SUK10

Endophytic Streptomyces strains are potential sources for novel bioactive molecules. In this study, the diketopiperazine gancidin W (GW) was isolated from the endophytic actinobacterial genus Streptomyces, SUK10, obtained from the bark of Shorea ovalis tree, and it was tested in vivo against Plasmodium berghei PZZ1/100. GW exhibited an inhibition rate of nearly 80% at 6.25 and 3.125 μg kg-1 body weight on day four using the 4-day suppression test method on male ICR strain mice. Comparing GW at both concentrations with quinine hydrochloride and normal saline as positive and negative controls, respectively, 50% of the mice treated with 3.125 μg kg-1 body weight managed to survive for more than 11 months after infection, which almost reached the life span of normal mice. Biochemical tests of selected enzymes and proteins in blood samples of mice treated with GW were also within normal levels; in addition, no abnormalities or injuries were found on internal vital organs. These findings indicated that this isolated bioactive compound from Streptomyces SUK10 exhibits very low toxicity and is a good candidate for potential use as an antimalarial agent in an animal model.

Population Parameters Underlying an Ongoing Soft Sweep in Southeast Asian Malaria Parasites

Multiple kelch13 alleles conferring artemisinin resistance (ART-R) are currently spreading through Southeast Asian malaria parasite populations, providing a unique opportunity to observe an ongoing soft selective sweep, investigate why resistance alleles have evolved multiple times and determine fundamental population genetic parameters for Plasmodium. We sequenced kelch13 (n = 1,876), genotyped 75 flanking SNPs, and measured clearance rate (n = 3,552) in parasite infections from Western Thailand (2001–2014). We describe 32 independent coding mutations including common mutations outside the kelch13 propeller associated with significant reductions in clearance rate. Mutations were first observed in 2003 and rose to 90% by 2014, consistent with a selection coefficient of ∼0.079. ART-R allele diversity rose until 2012 and then dropped as one allele (C580Y) spread to high frequency. The frequency with which adaptive alleles arise is determined by the rate of mutation and the population size. Two factors drive this soft sweep: (1) multiple kelch13 amino-acid mutations confer resistance providing a large mutational target—we estimate the target is 87–163 bp. (2) The population mutation parameter (Θ = 2N_eμ) can be estimated from the frequency distribution of ART-R alleles and is ∼5.69, suggesting that short term effective population size is 88 thousand to 1.2 million. This is 52–705 times greater than N_e estimated from fluctuation in allele frequencies, suggesting that we have previously underestimated the capacity for adaptive evolution in Plasmodium. Our central conclusions are that retrospective studies may underestimate the complexity of selective events and the N_e relevant for adaptation for malaria is considerably higher than previously estimated.

Mosquito Vectors and the Globalization of Plasmodium falciparum Malaria

Plasmodium falciparum malaria remains a devastating public health problem. Recent discoveries have shed light on the origin and evolution of Plasmodium parasites and their interactions with their vertebrate and mosquito hosts. P. falciparum malaria originated in Africa from a single horizontal transfer between an infected gorilla and a human, and became global as the result of human migration. Today, P. falciparum malaria is transmitted worldwide by more than 70 different anopheline mosquito species. Recent studies indicate that the mosquito immune system can be a barrier to malaria transmission and that the P. falciparum Pfs47 gene allows the parasite to evade mosquito immune detection. Here, we review the origin and globalization of P. falciparum and integrate this history with analysis of the biology, evolution, and dispersal of the main mosquito vectors. This new perspective broadens our understanding of P. falciparum population structure and the dispersal of important parasite genetic traits.

Mitochondrial DNA from the eradicated European Plasmodium vivax and P. falciparum from 70-year-old slides from the Ebro Delta in Spain

Phylogenetic analysis of Plasmodium parasites has indicated that their modern-day distribution is a result of a series of human-mediated dispersals involving transport between Africa, Europe, America, and Asia. A major outstanding question is the phylogenetic affinity of the malaria causing parasites Plasmodium vivax and falciparum in historic southern Europe-where it was endemic until the mid-20th century, after which it was eradicated across the region. Resolving the identity of these parasites will be critical for answering several hypotheses on the malaria dispersal. Recently, a set of slides with blood stains of malaria-affected people from the Ebro Delta (Spain), dated between 1942 and 1944, have been found in a local medical collection. We extracted DNA from three slides, two of them stained with Giemsa (on which Plasmodium parasites could still be seen under the microscope) and another one consisting of dried blood spots. We generated the data using Illumina sequencing after using several strategies aimed at increasing the Plasmodium DNA yield: depletion of the human genomic (g)DNA content through hybridization with human gDNA baits, and capture-enrichment using gDNA derived from P. falciparum Plasmodium mitochondrial genome sequences were subsequently reconstructed from the resulting data. Phylogenetic analysis of the eradicated European P. vivax mtDNA genome indicates that the European isolate is closely related to the most common present-day American haplotype and likely entered the American continent post-Columbian contact. Furthermore, the European P. falciparum mtDNA indicates a link with current Indian strains that is in agreement with historical accounts.

A comparative molecular survey of malaria prevalence among Eastern chimpanzee populations in Issa Valley (Tanzania) and Kalinzu (Uganda)

BACKGROUND

Habitat types can affect vector and pathogen distribution and transmission dynamics. The prevalence and genetic diversity of Plasmodium spp. in two eastern chimpanzee populations-Kalinzu Forest Reserve, Uganda and Issa Valley, Tanzania-inhabiting different habitat types was investigated. As a follow up study the effect of host sex and age on infections patterns in Kalinzu Forest Reserve chimpanzees was determined.

METHODS

Molecular methods were employed to detect Plasmodium DNA from faecal samples collected from savanna-woodland (Issa Valley) and forest (Kalinzu Forest Reserve) chimpanzee populations.

RESULTS

Based on a Cytochrome-b PCR assay, 32 out of 160 Kalinzu chimpanzee faecal samples were positive for Plasmodium DNA, whilst no positive sample was detected in 171 Issa Valley chimpanzee faecal samples. Sequence analysis revealed that previously known Laverania species (Plasmodium reichenowi, Plasmodium billbrayi and Plasmodium billcollinsi) are circulating in the Kalinzu chimpanzees. A significantly higher proportion of young individuals were tested positive for infections, and switching of Plasmodium spp. was reported in one individual. Amongst the positive individuals sampled more than once, the success of amplification of Plasmodium DNA from faeces varied over sampling time.

CONCLUSION

The study showed marked differences in the prevalence of malaria parasites among free ranging chimpanzee populations living in different habitats. In addition, a clear pattern of Plasmodium infections with respect to host age was found. The results presented in this study contribute to understanding the ecological aspects underlying the malaria infections in the wild. Nevertheless, integrative long-term studies on vector abundance, Plasmodium diversity during different seasons between sites would provide more insight on the occurrence, distribution and ecology of these pathogens.

Phylogeographic Evidence for 2 Genetically Distinct Zoonotic Plasmodium knowlesi Parasites, Malaysia

Infections of humans with the zoonotic simian malaria parasite Plasmodium knowlesi occur throughout Southeast Asia, although most cases have occurred in Malaysia, where P. knowlesi is now the dominant malaria species. This apparently skewed distribution prompted an investigation of the phylogeography of this parasite in 2 geographically separated regions of Malaysia, Peninsular Malaysia and Malaysian Borneo. We investigated samples collected from humans and macaques in these regions. Haplotype network analyses of sequences from 2 P. knowlesi genes, type A small subunit ribosomal 18S RNA and cytochrome c oxidase subunit I, showed 2 genetically distinct divergent clusters, 1 from each of the 2 regions of Malaysia. We propose that these parasites represent 2 distinct P. knowlesi types that independently became zoonotic. These types would have evolved after the sea-level rise at the end of the last ice age, which separated Malaysian Borneo from Peninsular Malaysia.

Genomic epidemiology of artemisinin resistant malaria

The current epidemic of artemisinin resistant Plasmodium falciparum in Southeast Asia is the result of a soft selective sweep involving at least 20 independent kelch13 mutations. In a large global survey, we find that kelch13 mutations which cause resistance in Southeast Asia are present at low frequency in Africa. We show that African kelch13 mutations have originated locally, and that kelch13 shows a normal variation pattern relative to other genes in Africa, whereas in Southeast Asia there is a great excess of non-synonymous mutations, many of which cause radical amino-acid changes. Thus, kelch13 is not currently undergoing strong selection in Africa, despite a deep reservoir of variations that could potentially allow resistance to emerge rapidly. The practical implications are that public health surveillance for artemisinin resistance should not rely on kelch13 data alone, and interventions to prevent resistance must account for local evolutionary conditions, shown by genomic epidemiology to differ greatly between geographical regions.

DOI:http://dx.doi.org/10.7554/eLife.08714.001

Malaria is an infectious disease caused by a microscopic parasite called Plasmodium, which is transferred between humans by mosquitos. One species of malaria parasite called Plasmodium falciparum can cause particularly severe and life-threatening forms of the disease. Currently, the most widely used treatment for P. falciparum infections is artemisinin combination therapy, a treatment that combines the drug artemisinin (or a closely related molecule) with another antimalarial drug. However, resistance to artemisinin has started to spread throughout Southeast Asia.

Artemisinin resistance is caused by mutations in a parasite gene called kelch13, and researchers have identified over 20 different mutations in P. falciparum that confer artemisinin resistance. The diversity of mutations involved, and the fact that the same mutation can arise independently in different locations, make it difficult to track the spread of resistance using conventional molecular marker approaches.

Here, Amato, Miotto et al. sequenced the entire genomes of more than 3,000 clinical samples of P. falciparum from Southeast Asia and Africa, collected as part of a global network of research groups called the MalariaGEN Plasmodium falciparum Community Project. Amato, Miotto et al. found that African parasites had independently acquired many of the same kelch13 mutations that are known to cause resistance to artemisinin in Southeast Asia. However the kelch13 mutations seen in Africa remained at low levels in the parasite population, and appeared to be under much less pressure for evolutionary selection than those found in Southeast Asia.

These findings demonstrate that the emergence and spread of resistance to antimalarial drugs does not depend solely on the mutational process, but also on other factors that influence whether the mutations will spread in the population. Understanding how this is affected by different patterns of drug treatments and other environmental conditions will be important in developing more effective strategies for combating malaria.

DOI:http://dx.doi.org/10.7554/eLife.08714.002

Polymorphisms of the artemisinin resistant marker (K13) in Plasmodium falciparum parasite populations of Grande Comore Island 10 years after artemisinin combination therapy

Background

Plasmodium falciparum malaria is a significant public health problem in Comoros, and artemisinin combination therapy (ACT) remains the first choice for treating acute uncomplicated P. falciparum. The emergence and spread of artemisinin-resistant P. falciparum in Southeast Asia, associated with mutations in K13-propeller gene, poses a potential threat to ACT efficacy. Detection of mutations in the P. falciparum K13-propeller gene may provide the first-hand information on changes in parasite susceptibility to artemisinin. The objective of this study is to determinate the prevalence of mutant K13-propeller gene among the P. falciparum isolates collected from Grande Comore Island, Union of Comoros, where ACT has been in use since 2004.

Methods

A total of 207 P. falciparum clinical isolates were collected from the island during March 2006 and October 2007 (n = 118) and March 2013 and December 2014 (n = 89). All isolates were analysed for single nucleotide polymorphisms (SNPs) and haplotypes in the K13-propeller gene using nested PCR and DNA sequencing.

Results

Only three 2006–2007 samples carried SNPs in the K13-propeller gene, one having a synonymous (G538G) and the other having two non-synonymous (S477Y and D584E) substitutions leading to two mutated haplotypes (2.2 %, 2/95). Three synonymous mutations (R471R, Y500Y, and G538G) (5.9 %, 5/85) and 7 non-synonymous substitutions (21.2 %, 18/85) with nine mutated haplotypes (18.8 %, 16/85) were found in isolates from 2013 to 2014. However, none of the polymorphisms associated with artemisinin-resistance in Southeast Asia was detected from any of the parasites examined.

Conclusion

This study showed increased K13-propeller gene diversity among P. falciparum populations on the Island over the course of 8 years (2006–2014). Nevertheless, none of the polymorphisms known to be associated with artemisinin resistance in Asia was detected in the parasite populations examined. Our data suggest that P. falciparum populations in Grande Comore are still effectively susceptible to artemisinin. Our results provide insights into P. falciparum populations regarding mutations in the gene associated with artemisinin resistance and will be useful for developing and updating anti-malarial guidance in Comoros.

Differential Spleen Remodeling Associated with Different Levels of Parasite Virulence Controls Disease Outcome in Malaria Parasite Infections

Infections by malaria parasites can lead to very different clinical outcomes, ranging from mild symptoms to death. Differences in the ability of the spleen to deal with the infected red blood cells (iRBCs) are linked to differences in virulence. Using virulent and avirulent strains of the rodent malaria parasite Plasmodium yoelii, we investigated how parasite virulence modulates overall spleen function. Following parasite invasion, a difference in parasite virulence was observed in association with different levels of spleen morphology and iRBC rigidity, both of which contributed to enhanced parasite clearance. Moreover, iRBC rigidity as modulated by the spleen was demonstrated to correlate with disease outcome and thus can be used as a robust indicator of virulence. The data indicate that alterations in the biomechanical properties of iRBCs are the result of the complex interaction between host and parasite. Furthermore, we confirmed that early spleen responses are a key factor in directing the clinical outcome of an infection. IMPORTANCE The spleen and its response to parasite infection are important in eliminating parasites in malaria. By comparing P. yoelii parasite lines with different disease outcomes in mice that had either intact spleens or had had their spleens removed, we showed that upon parasite infection, the spleen exhibits dramatic changes that can affect parasite clearance. The spleen itself directly impacts RBC deformability independently of parasite genetics. The data indicated that the changes in the biomechanical properties of malaria parasite-infected RBCs are the result of the complex interaction between host and parasite, and RBC deformability itself can serve as a novel predictor of clinical outcome. The results also suggest that early responses in the spleen are a key factor directing the clinical outcome of an infection.

In vivo antimalarial activity of the endophytic actinobacteria, Streptomyces SUK 10

Endophytic bacteria, such as Streptomyces, have the potential to act as a source for novel bioactive molecules with medicinal properties. The present study was aimed at assessing the antimalarial activity of crude extract isolated from various strains of actinobacteria living endophytically in some Malaysian medicinal plants. Using the four day suppression test method on male ICR strain mice, compounds produced from three strains of Streptomyces (SUK8, SUK10, and SUK27) were tested in vivo against Plasmodium berghei PZZ1/100 in an antimalarial screen using crude extracts at four different concentrations. One of these extracts, isolated from Streptomyces SUK10 obtained from the bark of Shorea ovalis tree, showed inhibition of the test organism and was further tested against P. berghei-infected mice for antimalarial activity at different concentrations. There was a positive relationship between the survival of the infected mouse group treated with 50 µg/kg body weight (bw) of ethyl acetate-SUK10 crude extract and the ability to inhibit the parasites growth. The parasite inhibition percentage for this group showed that 50% of the mice survived for more than 90 days after infection with the parasite. The nucleotide sequence and phylogenetic tree suggested that Streptomyces SUK10 may constitute a new species within the Streptomyces genus. As part of the drug discovery process, these promising finding may contribute to the medicinal and pharmaceutical field for malarial treatment.

Population Structure Shapes Copy Number Variation in Malaria Parasites

If copy number variants (CNVs) are predominantly deleterious, we would expect them to be more efficiently purged from populations with a large effective population size (Ne) than from populations with a small Ne. Malaria parasites (Plasmodium falciparum) provide an excellent organism to examine this prediction, because this protozoan shows a broad spectrum of population structures within a single species, with large, stable, outbred populations in Africa, small unstable inbred populations in South America and with intermediate population characteristics in South East Asia. We characterized 122 single-clone parasites, without prior laboratory culture, from malaria-infected patients in seven countries in Africa, South East Asia and South America using a high-density single-nucleotide polymorphism/CNV microarray. We scored 134 high-confidence CNVs across the parasite exome, including 33 deletions and 102 amplifications, which ranged in size from <500 bp to 59 kb, as well as 10,107 flanking, biallelic single-nucleotide polymorphisms. Overall, CNVs were rare, small, and skewed toward low frequency variants, consistent with the deleterious model. Relative to African and South East Asian populations, CNVs were significantly more common in South America, showed significantly less skew in allele frequencies, and were significantly larger. On this background of low frequency CNV, we also identified several high-frequency CNVs under putative positive selection using an FST outlier analysis. These included known adaptive CNVs containing rh2b and pfmdr1, and several other CNVs (e.g., DNA helicase and three conserved proteins) that require further investigation. Our data are consistent with a significant impact of genetic structure on CNV burden in an important human pathogen.

Potential antimalarials from African natural products: A reviw

Malaria remains an overwhelming infectious disease with significant health challenges in African and other endemic countries globally. Resistance to antimalarial drugs has become one of the most momentous challenges to human health, and thus has necessitated the hunt for new and effective drugs. Consequently, few decades have witnessed a surfeit of research geared to validate the effectiveness of commonly used traditionally medicines against malaria fever. The present review work focuses on documenting natural products from African whose activity has been reported in vivo or in vitro against malaria parasite. Literature was collected using electronic search of published articles (Google Scholar, PubMed, Medline, Sciencedirect, and Science domain) that report on antiplasmodial activity of natural products from differernts Africa region. A total of 652 plant taxa from 146 families, 134 isolated antimalarial compounds from 39 plants species, 2 herbal formulations and 4 insect/products were found to be reported in literature from 1996 to 2015. Plants species from family Asteraceae (11.04%), Fababceae (8.128%), Euphorbiaceae (5.52%), Rubiaceas (5.52%), and Apocyanaceae (5.214%), have received more scientific validation than others. African natural products possess remarkable healing properties as revealed in the various citations as promising antimalarial agents. Some of these natural products from Africa demonstrate high, promising or low activities against Plasmodium parasite. This study also shows that natural products from Africa have a huge amount of novel antimalarial compounds that could serve as a leads for the development of new and effective antiplasmodial drugs. However, in a view of bridging the gap in knowledge, clinical validation of these natural products are of paramount importance.

Genetic variability and population structure of Plasmodium falciparum parasite populations from different malaria ecological regions of Kenya

Transmission intensity, movement of human and vector hosts, biogeographical features, and malaria control measures are some of the important factors that determine Plasmodium falciparum parasite genetic variability and population structure. Kenya has different malaria ecologies which might require different disease intervention methods. Refined parasite population genetic studies are critical for informing malaria control and elimination strategies. This study describes the genetic diversity and population structure of P. falciparum parasites from the different malaria ecological zones in Kenya. Twelve multi-locus microsatellite (MS) loci previously described were genotyped in 225 P. falciparum isolates collected between 2012 and 2013 from five sites; three in lowland endemic regions (Kisumu, Kombewa, and Malindi) and two in highland, epidemic regions (Kisii and Kericho). Parasites from the lowland endemic and highland epidemic regions of western Kenya had high genetic diversity compared to coastal lowland endemic region of Kenya [Malindi]. The Kenyan parasites had a mean genetic differentiation index (FST) of 0.072 (p=0.011). The multi-locus genetic analysis of the 12 MS revealed all the parasites had unique haplotypes. Significant linkage disequilibrium (LD) was observed in all the five parasite populations. Kisumu had the most significant index of association values (0.16; p<0.0001) whereas Kisii had the least significant index of association values (0.03; p<0.0001). Our data suggest high genetic diversity in Kenyan parasite population with the exception of parasite from Malindi where malaria has been on the decline. The presence of significant LD suggests that there is occurrence of inbreeding in the parasite population. Parasite populations from Kisii showed the strongest evidence for epidemic population structure whereas the rest of the regions showed panmixia. Defining the genetic diversity of the parasites in different ecological regions of Kenya after introduction of the artemether-lumefantrine is important in refining the spread of drug resistant strains and malaria transmission for more effective control and eventual elimination of malaria in Kenya.

A novel locus of resistance to severe malaria in a region of ancient balancing selection

The high prevalence of sickle haemoglobin in Africa shows that malaria has been a major force for human evolutionary selection, but surprisingly few other polymorphisms have been proven to confer resistance to malaria in large epidemiological studies. To address this problem, we conducted a multi-centre genome-wide association study (GWAS) of life-threatening Plasmodium falciparum infection (severe malaria) in over 11,000 African children, with replication data in a further 14,000 individuals. Here we report a novel malaria resistance locus close to a cluster of genes encoding glycophorins that are receptors for erythrocyte invasion by P. falciparum. We identify a haplotype at this locus that provides 33% protection against severe malaria (odds ratio = 0.67, 95% confidence interval = 0.60-0.76, P value = 9.5 × 10(-11)) and is linked to polymorphisms that have previously been shown to have features of ancient balancing selection, on the basis of haplotype sharing between humans and chimpanzees. Taken together with previous observations on the malaria-protective role of blood group O, these data reveal that two of the strongest GWAS signals for severe malaria lie in or close to genes encoding the glycosylated surface coat of the erythrocyte cell membrane, both within regions of the genome where it appears that evolution has maintained diversity for millions of years. These findings provide new insights into the host-parasite interactions that are critical in determining the outcome of malaria infection.

The remarkable journey of adaptation of the Plasmodium falciparum malaria parasite to New World anopheline mosquitoes

Plasmodium falciparum originated in Africa, dispersed around the world as a result of human migration and had to adapt to several different indigenous anopheline mosquitoes. Anophelines from the New World are evolutionary distant form African ones and this probably resulted in a more stringent selection of Plasmodium as it adapted to these vectors. It is thought that Plasmodium has been genetically selected by some anopheline species through unknown mechanisms. The mosquito immune system can greatly limit infection and P. falciparum evolved a strategy to evade these responses, at least in part mediated by Pfs47, a highly polymorphic gene. We propose that adaptation of P. falciparum to new vectors may require evasion of their immune system. Parasites with a Pfs47 haplotype compatible with the indigenous mosquito vector would be able to survive and be transmitted. The mosquito antiplasmodial response could be an important determinant of P. falciparum population structure and could affect malaria transmission in the Americas.

Mitochondrial genome sequence diversity of Indian Plasmodium falciparum isolates

We have analysed the whole mitochondrial (mt) genome sequences (each ~6 kilo nucleotide base pairs in length) of four field isolates of the malaria parasite Plasmodium falciparum collected from different locations in India. Comparative genomic analyses of mt genome sequences revealed three novel India-specific single nucleotide polymorphisms. In general, high mt genome diversity was found in Indian P. falciparum, at a level comparable to African isolates. A population phylogenetic tree placed the presently sequenced Indian P. falciparum with the global isolates, while a previously sequenced Indian isolate was an outlier. Although this preliminary study is limited to a few numbers of isolates, the data have provided fundamental evidence of the mt genome diversity and evolutionary relationships of Indian P. falciparum with that of global isolates.

Mitochondrial population genomic analyses reveal population structure and demography of Indian Plasmodium falciparum

Inference on the genetic diversity of Plasmodium falciparum populations could help in better management of malaria. A very recent study with mitochondrial (mt) genomes in global P. falciparum had revealed interesting evolutionary genetic patterns of Indian isolates in comparison to global ones. However, no population genetic study using the whole mt genome sequences of P. falciparum isolates collected in the entire distribution range in India has yet been performed. We herewith have analyzed 85 whole mt genomes (48 already published and 37 entirely new) sampled from eight differentially endemic Indian locations to estimate genetic diversity and infer population structure and historical demography of Indian P. falciparum. We found 19 novel Indian-specific Single Nucleotide Polymorphisms (SNPs) and 22 novel haplotypes segregating in Indian P. falciparum. Accordingly, high haplotype and nucleotide diversities were detected in Indian P. falciparum in comparison to many other global isolates. Indian P. falciparum populations were found to be moderately sub-structured with four different genetic clusters. Interestingly, group of local populations aggregate to form each cluster; while samples from Jharkhand and Odisha formed a single cluster, P. falciparum isolates from Asom formed an independent one. Similarly, Surat, Bilaspur and Betul formed a single cluster and Goa and Mangalore formed another. Interestingly, P. falciparum isolates from the two later populations were significantly genetically differentiated from isolates collected in other six Indian locations. Signature of historical population expansion was evident in five population samples, and the onset of expansion event was found to be very similar to African P. falciparum. In agreement with the previous finding, the estimated Time to Most Recent Common Ancestor (TMRCA) and the effective population size were high in Indian P. falciparum. All these genetic features of Indian P. falciparum with high mt genome diversity are somehow similar to Africa, but quite different from other Asian population samples.

The evolutionary origins of Southeast Asian Ovalocytosis

Southeast Asian Ovalocytosis (SAO) is a common red blood cell disorder that is maintained as a balanced polymorphism in human populations. In individuals heterozygous for the SAO-causing mutation there are minimal detrimental effects and well-documented protection from severe malaria caused by Plasmodium vivax and Plasmodium falciparum; however, the SAO-causing mutation is fully lethal in utero when homozygous. The present-day high frequency of SAO in Island Southeast Asia indicates the trait is maintained by strong heterozygote advantage. Our study elucidates the evolutionary origin of SAO by characterizing DNA sequence variation in a 9.5 kilobase region surrounding the causal mutation in the SLC4A1 gene. We find substantial haplotype diversity among SAO chromosomes and estimate the age of the trait to be approximately 10,005 years (95% CI: 4930-23,200 years). This date is far older than any other human malaria-resistance trait examined previously in Southeast Asia, and considerably pre-dates the widespread adoption of agriculture associated with the spread of speakers of Austronesian languages some 4000 years ago. Using a genealogy-based method we find no evidence of historical positive selection acting on SAO (s=0.0, 95% CI: 0.0-0.03), in sharp contrast to the strong present-day selection coefficient (e.g., 0.09) estimated from the frequency of this recessively lethal trait. This discrepancy may be due to a recent increase in malaria-driven selection pressure following the spread of agriculture, with SAO targeted as a standing variant by positive selection in malarial populations.

Plasmodium vivax populations are more genetically diverse and less structured than sympatric Plasmodium falciparum populations

Introduction

The human malaria parasite, Plasmodium vivax, is proving more difficult to control and eliminate than Plasmodium falciparum in areas of co-transmission. Comparisons of the genetic structure of sympatric parasite populations may provide insight into the mechanisms underlying the resilience of P. vivax and can help guide malaria control programs.

Methodology/Principle findings

P. vivax isolates representing the parasite populations of four areas on the north coast of Papua New Guinea (PNG) were genotyped using microsatellite markers and compared with previously published microsatellite data from sympatric P. falciparum isolates. The genetic diversity of P. vivax (H_e = 0.83–0.85) was higher than that of P. falciparum (H_e = 0.64–0.77) in all four populations. Moderate levels of genetic differentiation were found between P. falciparum populations, even over relatively short distances (less than 50 km), with 21–28% private alleles and clear geospatial genetic clustering. Conversely, very low population differentiation was found between P. vivax catchments, with less than 5% private alleles and no genetic clustering observed. In addition, the effective population size of P. vivax (30353; 13043–69142) was larger than that of P. falciparum (18871; 8109–42986).

Conclusions/Significance

Despite comparably high prevalence, P. vivax had higher diversity and a panmictic population structure compared to sympatric P. falciparum populations, which were fragmented into subpopulations. The results suggest that in comparison to P. falciparum, P. vivax has had a long-term large effective population size, consistent with more intense and stable transmission, and limited impact of past control and elimination efforts. This underlines suggestions that more intensive and sustained interventions will be needed to control and eventually eliminate P. vivax. This research clearly demonstrates how population genetic analyses can reveal deeper insight into transmission patterns than traditional surveillance methods.

The neglected human malaria parasite Plasmodium vivax is responsible for a large proportion of the global malaria burden. Efforts to control malaria have revealed that P. vivax is more resilient than the other major human malaria parasite, Plasmodium falciparum. This study utilised population genetics to compare patterns of P. vivax and P. falciparum transmission in Papua New Guinea, a region where infection rates of the two species are similar. The results demonstrated that P. vivax populations are more genetically diverse than those of P. falciparum suggestive of a parasite population that is more resilient to environmental challenges, undergoing higher levels of interbreeding locally and between distant parasite populations. Unique characteristics of P. vivax such as relapse, which allows different strains from past infections to produce subsequent infections, may provide more opportunities for the exchange and dissemination of genetic material. The contrasting patterns observed for the two species may be the result of a differential impact of past elimination attempts and indicate that more rigorous interventions will be needed in efforts to control and eventually eliminate P. vivax.

Imputation-based population genetics analysis of Plasmodium falciparum malaria parasites

Whole-genome sequencing technologies are being increasingly applied to Plasmodium falciparum clinical isolates to identify genetic determinants of malaria pathogenesis. However, genome-wide discovery methods, such as haplotype scans for signatures of natural selection, are hindered by missing genotypes in sequence data. Poor correlation between single nucleotide polymorphisms (SNPs) in the P. falciparum genome complicates efforts to apply established missing-genotype imputation methods that leverage off patterns of linkage disequilibrium (LD). The accuracy of state-of-the-art, LD-based imputation methods (IMPUTE, Beagle) was assessed by measuring allelic r2 for 459 P. falciparum samples from malaria patients in 4 countries: Thailand, Cambodia, Gambia, and Malawi. In restricting our analysis to 86 k high-quality SNPs across the populations, we found that the complete-case analysis was restricted to 21k SNPs (24.5%), despite no single SNP having more than 10% missing genotypes. The accuracy of Beagle in filling in missing genotypes was consistently high across all populations (allelic r2, 0.87-0.96), but the performance of IMPUTE was mixed (allelic r2, 0.34-0.99) depending on reference haplotypes and population. Positive selection analysis using Beagle-imputed haplotypes identified loci involved in resistance to chloroquine (crt) in Thailand, Cambodia, and Gambia, sulfadoxine-pyrimethamine (dhfr, dhps) in Cambodia, and artemisinin (kelch13) in Cambodia. Tajima's D-based analysis identified genes under balancing selection that encode well-characterized vaccine candidates: apical merozoite antigen 1 (ama1) and merozoite surface protein 1 (msp1). In contrast, the complete-case analysis failed to identify any well-validated drug resistance or candidate vaccine loci, except kelch13. In a setting of low LD and modest levels of missing genotypes, using Beagle to impute P. falciparum genotypes is a viable strategy for conducting accurate large-scale population genetics and association analyses, and supporting global surveillance for drug resistance markers and candidate vaccine antigens.

Triosephosphate isomerase gene promoter variation: -5G/A and -8G/A polymorphisms in clinical malaria groups in two African populations

TPI1 promoter polymorphisms occur in high prevalence in individuals from African origin. Malaria-patients from Angola and Mozambique were screened for the TPI1 gene promoter variants rs1800200A>G, (-5G>A), rs1800201G>A, (-8G>A), rs1800202T>G, (-24T>G), and for the intron 5 polymorphism rs2071069G>A, (2262G>A). -5G>A and -8G>A variants occur in 47% and 53% in Angola and Mozambique, respectively while -24T>G was monomorphic for the wild-type T allele. Six haplotypes were identified and -8A occurred in 45% of the individuals, especially associated with the GAG haplotype and more frequent in non-severe malaria groups, although not significantly. The arising and dispersion of -5G>A and -8G>A polymorphisms is controversial. Their age was estimated by analyses of two microsatellite loci, CD4 and ATN1, adjacent to TPI1 gene. The -5G>A is older than -8G>A, with an average estimate of approximately 35,000 years. The -8A variant arose in two different backgrounds, suggesting independent mutational events. The first, on the -5G background, may have occurred in East Africa around 20,800 years ago; the second, on the -5A background, may have occurred in West Africa some 7500 years ago. These estimates are within the period of spread of agriculture and the malaria mosquito vector in Africa, which could has been a possible reason for the selection of -8A polymorphism in malaria endemic countries.

Hypothesis: dynamics of classical malaria epidemics show Plasmodium falciparum's survival strategy

Areas of marginal transmission can generate enormous lethal falciparum malaria epidemics when factors favoring the parasite shift only slightly. Although usually described in terms of vectorial capacity, medical scientists working in India in the early 20th century came to the conclusion that "an epidemic of relapses" was the key triggering event of malaria epidemics. This explanation has been largely discarded, because the biology of Plasmodium falciparum recrudescence has since been differentiated from P. vivax relapse. Using data from the Punjab in 1908 and Ceylon in 1934-1935, the genesis of malaria epidemics has been re-examined to inform current control efforts. The epidemics were focused geographically depending on recent rainfall or drought. Epidemics arose very suddenly and simultaneously in several places. Malaria spleen surveys indicated very little recent malaria transmission, and blood smears showed very few gametocytes just before the epidemic. Population stress as indicated by high grain prices because of a poor harvest caused by drought the previous year was a risk factor for malaria epidemics. Although increased female Anopheline survival because of increased humidity played an important part in the magnification of the epidemic, it does not explain its genesis. Human population stress triggering a shift toward gametocytogenesis is hypothesized as the key initiation factor for malaria epidemics. Its evolutionary significance may be that it allows the parasite to match the tropical agricultural cycle.

A new method for estimating species age supports the coexistence of malaria parasites and their Mammalian hosts

Species in the genus Plasmodium cause malaria in humans and infect a variety of mammals and other vertebrates. Currently, estimated ages for several mammalian Plasmodium parasites differ by as much as one order of magnitude, an inaccuracy that frustrates reliable estimation of evolutionary rates of disease-related traits. We developed a novel statistical approach to dating the relative age of evolutionary lineages, based on Total Least Squares regression. We validated this lineage dating approach by applying it to the genus Drosophila. Using data from the Drosophila 12 Genomes project, our approach accurately reconstructs the age of well-established Drosophila clades, including the speciation event that led to the subgenera Drosophila and Sophophora, and age of the melanogaster species subgroup. We applied this approach to hundreds of loci from seven mammalian Plasmodium species. We demonstrate the existence of a molecular clock specific to individual Plasmodium proteins, and estimate the relative age of mammalian-infecting Plasmodium. These analyses indicate that: 1) the split between the human parasite Plasmodium vivax and P. knowlesi, from Old World monkeys, occurred 6.1 times earlier than that between P. falciparum and P. reichenowi, parasites of humans and chimpanzees, respectively; and 2) mammalian Plasmodium parasites originated 22 times earlier than the split between P. falciparum and P. reichenowi. Calibrating the absolute divergence times for Plasmodium with eukaryotic substitution rates, we show that the split between P. falciparum and P. reichenowi occurred 3.0-5.5 Ma, and that mammalian Plasmodium parasites originated over 64 Ma. Our results indicate that mammalian-infecting Plasmodium evolved contemporaneously with their hosts, with little evidence for parasite host-switching on an evolutionary scale, and provide a solid timeframe within which to place the evolution of new Plasmodium species.

Arab gene geography: From population diversities to personalized medical genomics

Genetic disorders are not equally distributed over the geography of the Arab region. While a number of disorders have a wide geographical presence encompassing 10 or more Arab countries, almost half of these disorders occur in a single Arab country or population. Nearly, one-third of the genetic disorders in Arabs result from congenital malformations and chromosomal abnormalities, which are also responsible for a significant proportion of neonatal and perinatal deaths in Arab populations. Strikingly, about two-thirds of these diseases in Arab patients follow an autosomal recessive mode of inheritance. High fertility rates together with increased consanguineous marriages, generally noticed in Arab populations, tend to increase the rates of genetic and congenital abnormalities. Many of the nearly 500 genes studied in Arab people revealed striking spectra of heterogeneity with many novel and rare mutations causing large arrays of clinical outcomes. In this review we provided an overview of Arab gene geography, and various genetic abnormalities in Arab populations, including disorders of blood, metabolic, circulatory and neoplasm, and also discussed their associated molecules or genes responsible for the cause of these disorders. Although studying Arab-specific genetic disorders resulted in a high value knowledge base, approximately 35% of genetic diseases in Arabs do not have a defined molecular etiology. This is a clear indication that comprehensive research is required in this area to understand the molecular pathologies causing diseases in Arab populations.

Pooled sequencing and rare variant association tests for identifying the determinants of emerging drug resistance in malaria parasites

We explored the potential of pooled sequencing to swiftly and economically identify selective sweeps due to emerging artemisinin (ART) resistance in a South-East Asian malaria parasite population. ART resistance is defined by slow parasite clearance from the blood of ART-treated patients and mutations in the kelch gene (chr. 13) have been strongly implicated to play a role. We constructed triplicate pools of 70 slow-clearing (resistant) and 70 fast-clearing (sensitive) infections collected from the Thai–Myanmar border and sequenced these to high (∼150-fold) read depth. Allele frequency estimates from pools showed almost perfect correlation (Lin’s concordance = 0.98) with allele frequencies at 93 single nucleotide polymorphisms measured directly from individual infections, giving us confidence in the accuracy of this approach. By mapping genome-wide divergence (F_ST) between pools of drug-resistant and drug-sensitive parasites, we identified two large (>150 kb) regions (on chrs. 13 and 14) and 17 smaller candidate genome regions. To identify individual genes within these genome regions, we resequenced an additional 38 parasite genomes (16 slow and 22 fast-clearing) and performed rare variant association tests. These confirmed kelch as a major molecular marker for ART resistance (P = 6.03 × 10⁻⁶). This two-tier approach is powerful because pooled sequencing rapidly narrows down genome regions of interest, while targeted rare variant association testing within these regions can pinpoint the genetic basis of resistance. We show that our approach is robust to recurrent mutation and the generation of soft selective sweeps, which are predicted to be common in pathogen populations with large effective population sizes, and may confound more traditional gene mapping approaches.

Microsatellite markers reveal low levels of population sub-structuring of Plasmodium falciparum in southwestern Nigeria

BACKGROUND

Genetic diversity studies provide evidence of Plasmodium falciparum differentiation that could affect fitness and adaptation to drugs and target antigens for vaccine development. This study describes the genetic structure of P. falciparum populations in urban and rural sites from southwestern Nigeria.

METHODOLOGY

Ten neutral microsatellite loci were genotyped in 196 P. falciparum infections from three localities: Aramoko-Ekiti, a rural community; Lekki, an urban location and Badagry, a peri-urban border settlement. Analysis was performed on the genetic diversity, linkage disequilibrium, population structure and inter-population differentiation.

RESULTS

Allelic diversity values were similar across all populations, with mean expected heterozygosity (HE) values between 0.65 and 0.79. No matching multilocus haplotypes were found and analysis of multilocus LD showed no significant index of association. Genetic differentiation between populations was low (ΦPT = 0.017).

CONCLUSION

The absence of detectable population structure of P. falciparum in southwestern Nigeria is evident in the lack of significant differentiation between populations separated by about 200 km. This implies that a fairly uniform malaria control strategy may be effective over a wide geographic range in this highly endemic region. However, more wide-scale survey across the country will be required to inform malaria control in this large and densely populated endemic region.

Five-year tracking of Plasmodium falciparum allele frequencies in a holoendemic area with indistinct seasonal transitions

BACKGROUND

The renewed malaria eradication efforts require an understanding of the seasonal patterns of frequency of polymorphic variants in order to focus limited funds productively. Although cross-sectional studies in holoendemic areas spanning a single year could be useful in describing parasite genotype status at a given point, such information is inadequate in describing temporal trends in genotype polymorphisms. For Plasmodium falciparum isolates from Kisumu District Hospital, Plasmodium falciparum chloroquine-resistance transporter gene (Pfcrt-K76T) and P. falciparum multidrug resistance gene 1 (PfMDR1-N86Y), were analyzed for polymorphisms and parasitemia changes in the 53 months from March 2008 to August 2012. Observations were compared with prevailing climatic factors, including humidity, rainfall, and temperature.

METHODS

Parasitemia (the percentage of infected red blood cells per total red blood cells) was established by microscopy for P. falciparum malaria-positive samples. P. falciparum DNA was extracted from whole blood using a Qiagen DNA Blood Mini Kit. Single nucleotide polymorphism identification at positions Pfcrt-K76T and PfMDR1-N86Y was performed using real-time polymerase chain reaction and/or sequencing. Data on climatic variables were obtained from http://www.tutiempo.net/en/.

RESULTS

A total of 895 field isolates from 2008 (n=169), 2009 (n=161), 2010 (n=216), 2011 (n=223), and 2012 (n=126) showed large variations in monthly frequency of PfMDR1-N86Y and Pfcrt-K76T as the mutant genotypes decreased from 68.4%±15% and 38.1%±13% to 29.8%±18% and 13.3%±9%, respectively. The mean percentage of parasitemia was 2.61%±1.01% (coefficient of variation 115.86%; n=895). There was no correlation between genotype or parasitemia and climatic factors.

CONCLUSION

This study shows variability in the frequency of Pfcrt-K76T and PfMDR1-N86Y polymorphisms during the study period, bringing into focus the role of cross-sectional studies in describing temporal genotype trends. The lack of correlation between genotypes and climatic changes, especially precipitation, emphasizes the cost of investment in genotype change.

Origin and expansion of the Yunnan Shoot Borer, Tomicus yunnanensis (coleoptera: scolytinae): a mixture of historical natural expansion and contemporary human-mediated relocation

The Yunnan shoot borer, Tomicus yunnanensis, is a recently-discovered, aggressive pest of the Yunnan pine stands in southwestern China. Despite many bionomics studies and massive controlling efforts, research on its population genetics is extremely limited. The present study, aimed at investigating the origin and dispersal of this important forestry pest, analyzed the population genetic structure and demographic history using a mitochondrial cox1 gene fragment. Our results showed that T. yunnanensis most likely originated from the Central-Yunnan Altiplano, and the divergence time analysis placed the origin approximately 0.72 million-years ago. Host separation and specialization might have caused the speciation of T. yunnanensis. Genetic structure analyses identified two population groups, with six populations near the origin area forming one group and the remaining six populations from western and eastern Yunnan and southwestern Sichuan comprising the other. Divergence time analysis placed the split of the two groups at approximately 0.60 million-years ago, and haplotype phylogenetic tree, network, as well as migration rate suggested that populations of the latter group were established via a small number of individuals from the former one. Migration analysis also showed a certain degree of recent expansion from southwestern Sichuan to eastern Yunnan. Our findings implied that T. yunnanensis underwent both historical expansion and recent dispersal. The historical expansion may relate to the oscillation of regional climate due to glacial and interglacial periods in the Pleistocene, while human-mediated transportation of pine-wood material might have assisted the relocation and establishment of this pest in novel habitats.

A barcode of organellar genome polymorphisms identifies the geographic origin of Plasmodium falciparum strains

Malaria is a major public health problem that is actively being addressed in a global eradication campaign. Increased population mobility through international air travel has elevated the risk of re-introducing parasites to elimination areas and dispersing drug-resistant parasites to new regions. A simple genetic marker that quickly and accurately identifies the geographic origin of infections would be a valuable public health tool for locating the source of imported outbreaks. Here we analyse the mitochondrion and apicoplast genomes of 711 Plasmodium falciparum isolates from 14 countries, and find evidence that they are non-recombining and co-inherited. The high degree of linkage produces a panel of relatively few single-nucleotide polymorphisms (SNPs) that is geographically informative. We design a 23-SNP barcode that is highly predictive (~92%) and easily adapted to aid case management in the field and survey parasite migration worldwide.

Tracing the source of malarial infections is an important step towards monitoring and controlling the disease. Here, Preston et al. analyse sequence data from 711 isolates and design a genetic barcode based on combined mitochondrial and apicoplast genomes that is able to distinguish between malaria parasites isolated from different geographical regions.

New insights into the evolutionary history of Plasmodium falciparum from mitochondrial genome sequence analyses of Indian isolates

Estimating genetic diversity and inferring the evolutionary history of Plasmodium falciparum could be helpful in understanding origin and spread of virulent and drug-resistant forms of the malaria pathogen and therefore contribute to malaria control programme. Genetic diversity of the whole mitochondrial (mt) genome of P. falciparum sampled across the major distribution ranges had been reported, but no Indian P. falciparum isolate had been analysed so far, even though India is highly endemic to P. falciparum malaria. We have sequenced the whole mt genome of 44 Indian field isolates and utilized published data set of 96 genome sequences to present global genetic diversity and to revisit the evolutionary history of P. falciparum. Indian P. falciparum presents high genetic diversity with several characteristics of ancestral populations and shares many of the genetic features with African and to some extent Papua New Guinean (PNG) isolates. Similar to African isolates, Indian P. falciparum populations have maintained high effective population size and undergone rapid expansion in the past with oldest time to the most recent common ancestor (TMRCA). Interestingly, one of the four single nucleotide polymorphisms (SNPs) that differentiates P. falciparum from P. falciparum-like isolates (infecting non-human primates in Africa) was found to be segregating in five Indian P. falciparum isolates. This SNP was in tight linkage with other two novel SNPs that were found exclusively in these five Indian isolates. The results on the mt genome sequence analyses of Indian isolates on the whole add to the current understanding on the evolutionary history of P. falciparum.

Intracerebral bleed, right haemiparesis and seizures: an atypical presentation of vivax malaria

Falciparum malaria is notoriously known to produce life-threatening complications. Despite growing reports of chloroquine resistance and severe disease, vivax malaria continues to be viewed as a benign disease. We report a rare case of a 47-year-old healthy man from a malaria-endemic region, presenting with intracerebral bleed, right haemiparesis, aphasia and seizures following vivax malaria. This was successfully managed conservatively, without any neurosurgical intervention, with combination therapy of intravenous artesunate, oral hydroxychloquine and primaquine. In a country where Plasmodium vivax is responsible for majority of cases of malaria, it is high time the national malaria control programmes focus on the elimination of P. vivax in addition to its more dangerous counterpart, P. falciparum.

PlasmoView: a web-based resource to visualise global Plasmodium falciparum genomic variation

Malaria is a global public health challenge, with drug resistance a major barrier to disease control and elimination. To meet the urgent need for better treatments and vaccines, a deeper knowledge of Plasmodium biology and malaria epidemiology is required. An improved understanding of the genomic variation of malaria parasites, especially the most virulent Plasmodium falciparum (Pf) species, has the potential to yield new insights in these areas. High-throughput sequencing and genotyping is generating large amounts of genomic data across multiple parasite populations. The resulting ability to identify informative variants, particularly single-nucleotide polymorphisms (SNPs), will lead to the discovery of intra- and inter-population differences and thus enable the development of genetic barcodes for diagnostic assays and clinical studies. Knowledge of genetic variability underlying drug resistance and other differential phenotypes will also facilitate the identification of novel mutations and contribute to surveillance and stratified medicine applications. The PlasmoView interactive web-browsing tool enables the research community to visualise genomic variation and annotation (eg, biological function) in a geographic setting. The first release contains over 600,000 high-quality SNPs in 631 Pf isolates from laboratory strains and four malaria-endemic regions (West Africa, East Africa, Southeast Asia and Oceania).

Reprint of "Pharmacokinetic modelling of the anti-malarial drug artesunate and its active metabolite dihydroartemisinin"

A four compartment mechanistic mathematical model is developed for the pharmacokinetics of the commonly used anti-malarial drug artesunate and its principle metabolite dihydroartemisinin following oral administration of artesunate. The model is structurally unidentifiable unless additional constraints are imposed. Combinations of mechanistically derived constraints are considered to assess their effects on structural identifiability and on model fits. Certain combinations of the constraints give rise to locally or globally identifiable model structures. Initial validation of the model under various combinations of the constraints leading to identifiable model structures was performed against a dataset of artesunate and dihydroartemisinin concentration-time profiles of 19 malaria patients. When all the discussed constraints were imposed on the model, the resulting globally identifiable model structure was found to fit reasonably well to those patients with normal drug absorption profiles. However, there is wide variability in the fitted parameters and further investigation is warranted.

Genome-wide analysis of selection on the malaria parasite Plasmodium falciparum in West African populations of differing infection endemicity

Locally varying selection on pathogens may be due to differences in drug pressure, host immunity, transmission opportunities between hosts, or the intensity of between-genotype competition within hosts. Highly recombining populations of the human malaria parasite Plasmodium falciparum throughout West Africa are closely related, as gene flow is relatively unrestricted in this endemic region, but markedly varying ecology and transmission intensity should cause distinct local selective pressures. Genome-wide analysis of sequence variation was undertaken on a sample of 100 P. falciparum clinical isolates from a highly endemic region of the Republic of Guinea where transmission occurs for most of each year and compared with data from 52 clinical isolates from a previously sampled population from The Gambia, where there is relatively limited seasonal malaria transmission. Paired-end short-read sequences were mapped against the 3D7 P. falciparum reference genome sequence, and data on 136,144 single nucleotide polymorphisms (SNPs) were obtained. Within-population analyses identifying loci showing evidence of recent positive directional selection and balancing selection confirm that antimalarial drugs and host immunity have been major selective agents. Many of the signatures of recent directional selection reflected by standardized integrated haplotype scores were population specific, including differences at drug resistance loci due to historically different antimalarial use between the countries. In contrast, both populations showed a similar set of loci likely to be under balancing selection as indicated by very high Tajima’s D values, including a significant overrepresentation of genes expressed at the merozoite stage that invades erythrocytes and several previously validated targets of acquired immunity. Between-population F_ST analysis identified exceptional differentiation of allele frequencies at a small number of loci, most markedly for five SNPs covering a 15-kb region within and flanking the gdv1 gene that regulates the early stages of gametocyte development, which is likely related to the extreme differences in mosquito vector abundance and seasonality that determine the transmission opportunities for the sexual stage of the parasite.

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

Although the geographic origin of malaria cases imported into the United States can often be inferred from travel histories, these histories may be lacking or incomplete. We hypothesized that mitochondrial haplotypes could provide region-specific molecular barcodes for tracing the origin of imported Plasmodium vivax infections. An analysis of 348 mitochondrial genomes from worldwide parasites and new sequences from 69 imported malaria cases diagnosed across the United States allowed for a geographic assignment of most infections originating from the Americas, southeast Asia, east Asia, and Melanesia. However, mitochondrial lineages from Africa, south Asia, central Asia, and the Middle East, which altogether contribute the vast majority of imported malaria cases in the United States, were closely related to each other and could not be reliably assigned to their geographic origins. More mitochondrial genomes are required to characterize molecular barcodes of P. vivax from these regions.

Recurrent bottlenecks in the malaria life cycle obscure signals of positive selection

Detecting signals of selection in the genome of malaria parasites is a key to identify targets for drug and vaccine development. Malaria parasites have a unique life cycle alternating between vector and host organism with a population bottleneck at each transition. These recurrent bottlenecks could influence the patterns of genetic diversity and the power of existing population genetic tools to identify sites under positive selection. We therefore simulated the site-frequency spectrum of a beneficial mutant allele through time under the malaria life cycle. We investigated the power of current population genetic methods to detect positive selection based on the site-frequency spectrum as well as temporal changes in allele frequency. We found that a within-host selective advantage is difficult to detect using these methods. Although a between-host transmission advantage could be detected, the power is decreased when compared with the classical Wright–Fisher (WF) population model. Using an adjusted null site-frequency spectrum that takes the malaria life cycle into account, the power of tests based on the site-frequency spectrum to detect positive selection is greatly improved. Our study demonstrates the importance of considering the life cycle in genetic analysis, especially in parasites with complex life cycles.

Whole mitochondrial genome sequence of an Indian Plasmodium falciparum field isolate

Mitochondrial genome sequence of malaria parasites has served as a potential marker for inferring evolutionary history of the Plasmodium genus. In Plasmodium falciparum, the mitochondrial genome sequences from around the globe have provided important evolutionary understanding, but no Indian sequence has yet been utilized. We have sequenced the whole mitochondrial genome of a single P. falciparum field isolate from India using novel primers and compared with the 3D7 reference sequence and 1 previously reported Indian sequence. While the 2 Indian sequences were highly divergent from each other, the presently sequenced isolate was highly similar to the reference 3D7 strain.

Genetic structure of Plasmodium vivax using the merozoite surface protein 1 icb5-6 fragment reveals new hybrid haplotypes in southern Mexico

BACKGROUND

Plasmodium vivax is a protozoan parasite with an extensive worldwide distribution, being highly prevalent in Asia as well as in Mesoamerica and South America. In southern Mexico, P. vivax transmission has been endemic and recent studies suggest that these parasites have unique biological and genetic features. The msp1 gene has shown high rate of nucleotide substitutions, deletions, insertions, and its mosaic structure reveals frequent events of recombination, maybe between highly divergent parasite isolates.

METHODS

The nucleotide sequence variation in the polymorphic icb5-6 fragment of the msp1 gene of Mexican and worldwide isolates was analysed. To understand how genotype diversity arises, disperses and persists in Mexico, the genetic structure and genealogical relationships of local isolates were examined. To identify new sequence hybrids and their evolutionary relationships with other P. vivax isolates circulating worldwide two haplotype networks were constructed questioning that two portions of the icb5-6 have different evolutionary history.

RESULTS

Twelve new msp1 icb5-6 haplotypes of P. vivax from Mexico were identified. These nucleotide sequences show mosaic structure comprising three partially conserved and two variable subfragments and resulted into five different sequence types. The variable subfragment sV1 has undergone recombination events and resulted in hybrid sequences and the haplotype network allocated the Mexican haplotypes to three lineages, corresponding to the Sal I and Belem types, and other more divergent group. In contrast, the network from icb5-6 fragment but not sV1 revealed that the Mexican haplotypes belong to two separate lineages, none of which are closely related to Sal I or Belem sequences.

CONCLUSIONS

These results suggest that the new hybrid haplotypes from southern Mexico were the result of at least three different recombination events. These rearrangements likely resulted from the recombination between haplotypes of highly divergent lineages that are frequently distributed in South America and Asia and diversified rapidly.

Pharmacokinetic modelling of the anti-malarial drug artesunate and its active metabolite dihydroartemisinin

A four compartment mechanistic mathematical model is developed for the pharmacokinetics of the commonly used anti-malarial drug artesunate and its principle metabolite dihydroartemisinin following oral administration of artesunate. The model is structurally unidentifiable unless additional constraints are imposed. Combinations of mechanistically derived constraints are considered to assess their effects on structural identifiability and on model fits. Certain combinations of the constraints give rise to locally or globally identifiable model structures. Initial validation of the model under various combinations of the constraints leading to identifiable model structures was performed against a dataset of artesunate and dihydroartemisinin concentration-time profiles of 19 malaria patients. When all the discussed constraints were imposed on the model, the resulting globally identifiable model structure was found to fit reasonably well to those patients with normal drug absorption profiles. However, there is wide variability in the fitted parameters and further investigation is warranted.

Plasmodium falciparum: genetic diversity and complexity of infections in an isolated village in western Thailand

Genetic diversity of Plasmodium falciparum is intimately associated with morbidity, mortality and malaria control strategies. It is therefore imperative to study genetic makeup and population structure of this parasite in endemic areas. In Kong Mong Tha, an isolated village in western Thailand, the majority of P. falciparum infections are asymptomatic. In this study we investigated complexity of infections and single nucleotide polymorphisms (SNPs) in the P. falciparum population of Kong Mong Tha, and compared results with those previously obtained from Mae Sod, in northwestern Thailand, where the majority of infections were symptomatic. Using PCR-based determination of the 5' merozoite surface protein 1 gene (msp1) recombinant types, we found that 39% of 59 P. falciparum isolates from Kong Mong Tha had multiple 5' recombinant types with a mean number of 1.54. These values were much lower than those obtained from Mae Sod: 96% for multiple infections and with a mean number of 3.61. Analysis of full-length sequences of two housekeeping genes, the P-type Ca(2+)-transporting ATPase gene (n=33) plus adenylosuccinate lyase gene (n=33), and three vaccine candidate antigen genes, msp1 (n=26), the circumsporozoite protein gene, csp (n=30) and the apical membrane antigen 1 gene, ama 1 (n=32), revealed that in all of these genes within-population SNP diversity was at similar levels between Kong Mong Tha and Mae Sod, suggesting that the extent of MOI and clinical manifestations of malaria are not strongly associated with genetic diversity. Additionally, we did not detect significant genetic differentiation between the two parasite populations, as estimated by the Wright's fixation index of inter-population variance in allele frequencies, suggesting that gene flow prevented the formation of population structuring. Thus, this study highlights unique features of P. falciparum populations in Thailand. The implications of these finding are discussed.

Malarial parasite diversity in chimpanzees: the value of comparative approaches to ascertain the evolution of Plasmodium falciparum antigens

BACKGROUND

Plasmodium falciparum shares its most recent common ancestor with parasites found in African apes; these species constitute the so-called Laverania clade. In this investigation, the evolutionary history of Plasmodium lineages found in chimpanzees (Pan troglodytes) was explored.

METHODS

Here, the remainders of 74 blood samples collected as part of the chimpanzees' routine health examinations were studied. For all positive samples with parasite lineages belonging to the Laverania clade, the complete mitochondrial genome (mtDNA), the gene encoding dihydrofolate reductase-thymidylate synthase (dhfr-ts), the chloroquine resistance transporter (Pfcrt), the circumsporozoite protein (csp), merozoite surface protein 2 (msp2), and the DBL-1 domain from var2CSA were amplified, cloned, and sequenced. Other Plasmodium species were included in the mtDNA, dhfr-ts, and csp analyses. Phylogenetic and evolutionary genetic analyses were performed, including molecular clock analyses on the mtDNA.

RESULTS/CONCLUSIONS

Nine chimpanzees were malaria positive (12.2%); four of those infections were identified as P. falciparum, two as a Plasmodium reichenowi-like parasite or Plasmodium sp., one as Plasmodium gaboni, and two as Plasmodium malariae. All P. falciparum isolates were resistant to chloroquine indicating that the chimpanzees acquired such infections from humans in recent times. Such findings, however, are not sufficient for implicating chimpanzees as an animal reservoir for P. falciparum.Timing estimates support that the Laverania clade has co-existed with hominids for a long-period of time. The proposed species P. gaboni, Plasmodium billbrayi, and Plasmodium billcollinsi are monophyletic groups supporting that they are indeed different species.An expanded CSP phylogeny is presented, including all the Laverania species and other malarial parasites. Contrasting with other Plasmodium, the Laverania csp exhibits great conservation at the central tandem repeat region. Msp2 and var2CSA, however, show extended recent polymorphism in P. falciparum that likely originated after the P. reichenowi-P. falciparum split. The accumulation of such diversity may indicate adaptation to the human host. These examples support the notion that comparative approaches among P. falciparum and its related species will be of great value in understanding the evolution of proteins that are important in parasite invasion of the human red blood cell, as well as those involved in malaria pathogenesis.