Early origin and recent expansion of Plasmodium falciparum

Plasmodium chabaudi chabaudi malaria parasites can develop stable resistance to atovaquone with a mutation in the cytochrome b gene

BACKGROUND

Plasmodium falciparum, has developed resistance to many of the drugs in use. The recommended treatment policy is now to use drug combinations. The atovaquone-proguanil (AP) drug combination, is one of the treatment and prophylaxis options. Atovaquone (ATQ) exerts its action by inhibiting plasmodial mitochondria electron transport at the level of the cytochrome bc1 complex. Plasmodium falciparum in vitro resistance to ATQ has been associated with specific point mutations in the region spanning codons 271-284 of the cytochrome b gene. ATQ -resistant Plasmodium yoelii and Plasmodium berghei lines have been obtained and resistant lines have amino acid mutations in their CYT b protein sequences. Plasmodium chabaudi model for studying drug-responses and drug-resistance selection is a very useful rodent malaria model but no ATQ resistant parasites have been reported so far. The aim of this study was to determine the ATQ sensitivity of the P. chabaudi clones, to select a resistant parasite line and to perform genotypic characterization of the cytb gene of these clones.

METHODS

To select for ATQ resistance, Plasmodium. chabaudi chabaudi clones were exposed to gradually increasing concentrations of ATQ during several consecutive passages in mice. Plasmodium chabaudi cytb gene was amplified and sequenced.

RESULTS

ATQ resistance was selected from the clone AS-3CQ. In order to confirm whether an heritable genetic mutation underlies the response of AS-ATQ to ATQ, the stability of the drug resistance phenotype in this clone was evaluated by measuring drug responses after (i) multiple blood passages in the absence of the drug, (ii) freeze/thawing of parasites in liquid nitrogen and (iii) transmission through a mosquito host, Anopheles stephensi. ATQ resistance phenotype of the drug-selected parasite clone kept unaltered. Therefore, ATQ resistance in clone AS-ATQ is genetically encoded. The Minimum Curative Dose of AS-ATQ showed a six-fold increase in MCD to ATQ relative to AS-3CQ.

CONCLUSIONS

A mutation was found on the P. chabaudi cytb gene from the AS-ATQ sample a substitution at the residue Tyr268 for an Asn, this mutation is homologous to the one found in P. falciparum isolates resistant to ATQ.

Detecting signatures of balancing selection to identify targets of anti-parasite immunity

Parasite antigen genes might evolve under frequency-dependent immune selection. The distinctive patterns of polymorphism that result can be detected using population genetic methods that test for signatures of balancing selection, allowing genes encoding important targets of immunity to be identified. Analyses can be complicated by population structures, histories and features of a parasite's genome. However, new sequencing technologies facilitate scans of polymorphism throughout parasite genomes to identify the most exceptional gene specific signatures. We focus on malaria parasites to illustrate challenges and opportunities for detecting targets of frequency-dependent immune selection to discover new potential vaccine candidates.

Evidence for negative selection on the gene encoding rhoptry-associated protein 1 (RAP-1) in Plasmodium spp

Assessing how natural selection, negative or positive, operates on genes with low polymorphism is challenging. We investigated the genetic diversity of orthologous genes encoding the rhoptry-associated protein 1 (RAP-1), a low polymorphic protein of malarial parasites that is involved in erythrocyte invasion. We applied evolutionary genetic methods to study the polymorphism in RAP-1 from Plasmodium falciparum (n=32) and Plasmodium vivax (n=6), the two parasites responsible for most human malaria morbidity and mortality, as well as RAP-1 orthologous in closely related malarial species found in non-human primates (NHPs). Overall, genes encoding RAP-1 are highly conserved in all Plasmodium spp. included in this investigation. We found no evidence for natural selection, positive or negative, acting on the gene encoding RAP-1 in P. falciparum or P. vivax. However, we found evidence that the orthologous genes in non-human primate parasites (Plasmodium cynomolgi, Plasmodium inui, and Plasmodium knowlesi) are under purifying (negative) selection. We discuss the importance of considering negative selection while studying genes encoding proteins with low polymorphism and how selective pressures may differ among orthologous genes in closely related malarial parasites species.

Malaria infection and the anthropological evolution

During the evolution of the genus Homo, with regard to species habilis, erectus and sapiens, malaria infection played a key biological role, influencing the anthropological development too. Plasmodia causing malaria developed two kinds of evolution, according to a biological and philogenetical point of view. In particular, Plasmodium vivax, Plasmodium malariae, and Plasmodium ovale, would have either coevolved with human mankind (coevolution), or reached human species during the most ancient phases of genus Homo evolution. On the other hand, Plasmodium falciparum has been transmitted to humans by monkeys in a more recent period, probably between the end of Mesolithic and the beginning of Neolithic age. The authors show both direct and indirect biomolecular evidences of malaria infection, detected in buried subjects, dating to the Ancient World, and brought to light in the course of archeological excavations in some relevant Mediterranean sites. In this literature review the Authors organize present scientific evidences: these confirm the malarial role in affecting the evolution of populations in Mediterranean countries. The people living in several different regions on the Mediterranean Sea sides, the cradle of western civilization, have been progressively influenced by malaria, in the course of the spread of this endemic disease during the last millennia. In addition, populations affected by endemic malaria developed cultural, dietary and behaviour adaptations, contributing to decrease the risk of disease. These habits were not probably fully conscious. Nevertheless it may be thought that both these customs and biological modifications, caused by malarial plasmodia, favoured the emergence of groups of people with a greater resistance against malaria. All these considered factors decreased demographical impact, influencing in a favourable way the general development and growth of civilization.

Geographic structuring of the Plasmodium falciparum sarco(endo)plasmic reticulum Ca2+ ATPase (PfSERCA) gene diversity

Artemisinin, a thapsigargin-like sesquiterpene has been shown to inhibit the Plasmodium falciparum sarco/endoplasmic reticulum calcium-ATPase PfSERCA. To collect baseline pfserca sequence information before field deployment of Artemisinin-based Combination therapies that may select mutant parasites, we conducted a sequence analysis of 100 isolates from multiple sites in Africa, Asia and South America. Coding sequence diversity was large, with 29 mutated codons, including 32 SNPs (average of one SNP/115 bp), of which 19 were novel mutations. Most SNP detected in this study were clustered within a region in the cytosolic head of the protein. The PfSERCA functional domains were very well conserved, with non synonymous mutations located outside the functional domains, except for the S769N mutation associated in French Guiana with elevated IC₅₀ for artemether. The S769N mutation is located close to the hinge of the headpiece, which in other species modulates calcium affinity and in consequence efficacy of inhibitors, possibly linking calcium homeostasis to drug resistance. Genetic diversity was highest in Senegal, Brazil and French Guiana, and few mutations were identified in Asia. Population genetic analysis was conducted for a partial fragment of the gene encompassing nucleotide coordinates 87-2862 (unambiguous sequence available for 96 isolates). This supported a geographic clustering, with a separation between Old and New World samples and one dominant ancestral haplotype. Genetic drift alone cannot explain the observed polymorphism, suggesting that other evolutionary mechanisms are operating. One possible contributor could be the frequency of haemoglobinopathies that are associated with calcium dysregulation in the erythrocyte.

Lineage-specific positive selection at the merozoite surface protein 1 (msp1) locus of Plasmodium vivax and related simian malaria parasites

Background

The 200 kDa merozoite surface protein 1 (MSP-1) of malaria parasites, a strong vaccine candidate, plays a key role during erythrocyte invasion and is a target of host protective immune response. Plasmodium vivax, the most widespread human malaria parasite, is closely related to parasites that infect Asian Old World monkeys, and has been considered to have become a parasite of man by host switch from a macaque malaria parasite. Several Asian monkey parasites have a range of natural hosts. The same parasite species shows different disease manifestations among host species. This suggests that host immune responses to P. vivax-related malaria parasites greatly differ among host species (albeit other factors). It is thus tempting to invoke that a major immune target parasite protein such as MSP-1 underwent unique evolution, depending on parasite species that exhibit difference in host range and host specificity.

Results

We performed comparative phylogenetic and population genetic analyses of the gene encoding MSP-1 (msp1) from P. vivax and nine P. vivax-related simian malaria parasites. The inferred phylogenetic tree of msp1 significantly differed from that of the mitochondrial genome, with a striking displacement of P. vivax from a position close to P. cynomolgi in the mitochondrial genome tree to an outlier of Asian monkey parasites. Importantly, positive selection was inferred for two ancestral branches, one leading to P. inui and P. hylobati and the other leading to P. vivax, P. fieldi and P. cynomolgi. This ancestral positive selection was estimated to have occurred three to six million years ago, coinciding with the period of radiation of Asian macaques. Comparisons of msp1 polymorphisms between P. vivax, P. inui and P. cynomolgi revealed that while some positively selected amino acid sites or regions are shared by these parasites, amino acid changes greatly differ, suggesting that diversifying selection is acting species-specifically on msp1.

Conclusions

The present results indicate that the msp1 locus of P. vivax and related parasite species has lineage-specific unique evolutionary history with positive selection. P. vivax and related simian malaria parasites offer an interesting system toward understanding host species-dependent adaptive evolution of immune-target surface antigen genes such as msp1.

On the diversity of malaria parasites in African apes and the origin of Plasmodium falciparum from Bonobos

The origin of Plasmodium falciparum, the etiological agent of the most dangerous forms of human malaria, remains controversial. Although investigations of homologous parasites in African Apes are crucial to resolve this issue, studies have been restricted to a chimpanzee parasite related to P. falciparum, P. reichenowi, for which a single isolate was available until very recently. Using PCR amplification, we detected Plasmodium parasites in blood samples from 18 of 91 individuals of the genus Pan, including six chimpanzees (three Pan troglodytes troglodytes, three Pan t. schweinfurthii) and twelve bonobos (Pan paniscus). We obtained sequences of the parasites' mitochondrial genomes and/or from two nuclear genes from 14 samples. In addition to P. reichenowi, three other hitherto unknown lineages were found in the chimpanzees. One is related to P. vivax and two to P. falciparum that are likely to belong to distinct species. In the bonobos we found P. falciparum parasites whose mitochondrial genomes indicated that they were distinct from those present in humans, and another parasite lineage related to P. malariae. Phylogenetic analyses based on this diverse set of Plasmodium parasites in African Apes shed new light on the evolutionary history of P. falciparum. The data suggested that P. falciparum did not originate from P. reichenowi of chimpanzees (Pan troglodytes), but rather evolved in bonobos (Pan paniscus), from which it subsequently colonized humans by a host-switch. Finally, our data and that of others indicated that chimpanzees and bonobos maintain malaria parasites, to which humans are susceptible, a factor of some relevance to the renewed efforts to eradicate malaria.

The changing disease-scape in the third epidemiological transition

The epidemiological transition model describes the changing relationship between humans and their diseases. The first transition occurred with the shift to agriculture about 10,000 YBP, resulting in a pattern of infectious and nutritional diseases still evident today. In the last two centuries, some populations have undergone a second transition, characterized by a decline in infectious disease and rise in degenerative disease. We are now in the throes of a third epidemiological transition, in which a resurgence of familiar infections is accompanied by an array of novel diseases, all of which have the potential to spread rapidly due to globalization.

Positively selected G6PD-Mahidol mutation reduces Plasmodium vivax density in Southeast Asians

Glucose-6-phosphate dehydrogenase (G6PD) deficiency--the most common known enzymopathy--is associated with neonatal jaundice and hemolytic anemia usually after exposure to certain infections, foods, or medications. Although G6PD-deficient alleles appear to confer a protective effect against malaria, the link with clinical protection from Plasmodium infection remains unclear. We investigated the effect of a common G6PD deficiency variant in Southeast Asia--the G6PD-Mahidol(487A) variant--on human survival related to vivax and falciparum malaria. Our results show that strong and recent positive selection has targeted the Mahidol variant over the past 1500 years. We found that the G6PD-Mahidol(487A) variant reduces vivax, but not falciparum, parasite density in humans, which indicates that Plasmodium vivax has been a driving force behind the strong selective advantage conferred by this mutation.

Contrasting population structures of the genes encoding ten leading vaccine-candidate antigens of the human malaria parasite, Plasmodium falciparum

The extensive diversity of Plasmodium falciparum antigens is a major obstacle to a broadly effective malaria vaccine but population genetics has rarely been used to guide vaccine design. We have completed a meta-population genetic analysis of the genes encoding ten leading P. falciparum vaccine antigens, including the pre-erythrocytic antigens csp, trap, lsa1 and glurp; the merozoite antigens eba175, ama1, msp's 1, 3 and 4, and the gametocyte antigen pfs48/45. A total of 4553 antigen sequences were assembled from published data and we estimated the range and distribution of diversity worldwide using traditional population genetics, Bayesian clustering and network analysis. Although a large number of distinct haplotypes were identified for each antigen, they were organized into a limited number of discrete subgroups. While the non-merozoite antigens showed geographically variable levels of diversity and geographic restriction of specific subgroups, the merozoite antigens had high levels of diversity globally, and a worldwide distribution of each subgroup. This shows that the diversity of the non-merozoite antigens is organized by physical or other location-specific barriers to gene flow and that of merozoite antigens by features intrinsic to all populations, one important possibility being the immune response of the human host. We also show that current malaria vaccine formulations are based upon low prevalence haplotypes from a single subgroup and thus may represent only a small proportion of the global parasite population. This study demonstrates significant contrasts in the population structure of P. falciparum vaccine candidates that are consistent with the merozoite antigens being under stronger balancing selection than non-merozoite antigens and suggesting that unique approaches to vaccine design will be required. The results of this study also provide a realistic framework for the diversity of these antigens to be incorporated into the design of next-generation malaria vaccines.

From evolutionary genetics to human immunology: how selection shapes host defence genes

Pathogens have always been a major cause of human mortality, so they impose strong selective pressure on the human genome. Data from population genetic studies, including genome-wide scans for selection, are providing important insights into how natural selection has shaped immunity and host defence genes in specific human populations and in the human species as a whole. These findings are helping to delineate genes that are important for host defence and to increase our understanding of how past selection has had an impact on disease susceptibility in modern populations. A tighter integration between population genetic studies and immunological phenotype studies is now necessary to reveal the mechanisms that have been crucial for our past and present survival against infection.

The dynamics of mutations associated with anti-malarial drug resistance in Plasmodium falciparum

The evolution of resistance in Plasmodium falciparum against safe and affordable drugs such as chloroquine (CQ) and sulfadoxine-pyrimethamine (SP) is a major global health threat. Investigating the dynamics of resistance against these antimalarial drugs will lead to approaches for addressing the problem of resistance in malarial parasites that are solidly based in evolutionary genetics and population biology. In this article, we discuss current developments in population biology modeling and evolutionary genetics. Despite great advancements achieved in the past decade, understanding the complex dynamics of mutations conferring drug resistance in P. falciparum requires approaches that consider the parasite population structure among other demographic processes.

The origin of malignant malaria

Plasmodium falciparum, the causative agent of malignant malaria, is among the most severe human infectious diseases. The closest known relative of P. falciparum is a chimpanzee parasite, Plasmodium reichenowi, of which one single isolate was previously known. The co-speciation hypothesis suggests that both parasites evolved separately from a common ancestor over the last 5-7 million years, in parallel with the divergence of their hosts, the hominin and chimpanzee lineages. Genetic analysis of eight new isolates of P. reichenowi, from wild and wild-born captive chimpanzees in Cameroon and Côte d'Ivoire, shows that P. reichenowi is a geographically widespread and genetically diverse chimpanzee parasite. The genetic lineage comprising the totality of global P. falciparum is fully included within the much broader genetic diversity of P. reichenowi. This finding is inconsistent with the co-speciation hypothesis. Phylogenetic analysis indicates that all extant P. falciparum populations originated from P. reichenowi, likely by a single host transfer, which may have occurred as early as 2-3 million years ago, or as recently as 10,000 years ago. The evolutionary history of this relationship may be explained by two critical genetic mutations. First, inactivation of the CMAH gene in the human lineage rendered human ancestors unable to generate the sialic acid Neu5Gc from its precursor Neu5Ac, and likely made humans resistant to P. reichenowi. More recently, mutations in the dominant invasion receptor EBA 175 in the P. falciparum lineage provided the parasite with preference for the overabundant Neu5Ac precursor, accounting for its extreme human pathogenicity.

Evolutionary analysis of the dynamics of viral infectious disease

Many organisms that cause infectious diseases, particularly RNA viruses, mutate so rapidly that their evolutionary and ecological behaviours are inextricably linked. Consequently, aspects of the transmission and epidemiology of these pathogens are imprinted on the genetic diversity of their genomes. Large-scale empirical analyses of the evolutionary dynamics of important pathogens are now feasible owing to the increasing availability of pathogen sequence data and the development of new computational and statistical methods of analysis. In this Review, we outline the questions that can be answered using viral evolutionary analysis across a wide range of biological scales.

Paraphyly of organelle DNAs in Cycas Sect. Asiorientales due to ancient ancestral polymorphisms

Background

This study addresses the apportionment of genetic diversity between Cycas revoluta and C. taitungensis, species that constitute the section Asiorientales and represent a unique, basal lineage of the Laurasian genus Cycas. Fossil evidence indicates divergence of the section from the rest of Cycas at least 30 million years ago. Geographically, C. taitungensis is limited to Taiwan whereas C. revoluta is found in the Ryukyu Archipelago and on mainland China.

Results

The phylogenies of ribosomal ITS region of mtDNA and the intergenic spacer between atpB and rbcL genes of cpDNA were reconstructed. Phylogenetic analyses revealed paraphyly of both loci in the two species and also in the section Asiorientales. The lack of reciprocal monophyly between these long isolated sections is likely due to persistent shared ancestral polymorphisms. Molecular dating estimated that mt- and cp DNA lineages coalesced to the most recent common ancestors (TMRCA) about 327 (mt) and 204 MYA (cp), corresponding with the divergence of cycad sections in the Mesozoic.

Conclusion

Fates of newly derived mutations of cycads follow Klopfstein et al.'s surfing model where the majority of new mutations do not spread geographically and remain at low frequencies or are eventually lost by genetic drift. Only successful 'surfing mutations' reach very high frequencies and occupy a large portion of a species range. These mutations exist as dominant cytotypes across populations and species. Geographical subdivision is lacking in both species, even though recurrent gene flow by both pollen and seed is severely limited. In total, the contrasting levels between historical and ongoing gene flow, large population sizes, a long lifespan, and slow mutation rates in both organelle DNAs have all likely contributed to the unusually long duration of paraphyly in cycads.

A new malaria agent in African hominids

Plasmodium falciparum is the major human malaria agent responsible for 200 to 300 million infections and one to three million deaths annually, mainly among African infants. The origin and evolution of this pathogen within the human lineage is still unresolved. A single species, P. reichenowi, which infects chimpanzees, is known to be a close sister lineage of P. falciparum. Here we report the discovery of a new Plasmodium species infecting Hominids. This new species has been isolated in two chimpanzees (Pan troglodytes) kept as pets by villagers in Gabon (Africa). Analysis of its complete mitochondrial genome (5529 nucleotides including Cyt b, Cox I and Cox III genes) reveals an older divergence of this lineage from the clade that includes P. falciparum and P. reichenowi (∼21±9 Myrs ago using Bayesian methods and considering that the divergence between P. falciparum and P. reichenowi occurred 4 to 7 million years ago as generally considered in the literature). This time frame would be congruent with the radiation of hominoids, suggesting that this Plasmodium lineage might have been present in early hominoids and that they may both have experienced a simultaneous diversification. Investigation of the nuclear genome of this new species will further the understanding of the genetic adaptations of P. falciparum to humans. The risk of transfer and emergence of this new species in humans must be now seriously considered given that it was found in two chimpanzees living in contact with humans and its close relatedness to the most virulent agent of malaria.

In 2002, the publication of the genome of Plasmodium falciparum, the most malignant agent of malaria, generated hopes in the fight against this deadly disease by the opportunities it offered to discover new drug targets. Since then results have not lived up to the expectations. The development of comparative genomics to further understanding of P. falciparum has indeed been hindered by a lack of knowledge of closely related species' genomes. Only one species, P. reichenowi, infecting chimpanzees, was hitherto known as a sister lineage of P. falciparum. Here we describe a new Plasmodium species infecting chimpanzees in Africa. Based on its whole mitochondrial genome, we demonstrate that this species is a relative of P. falciparum and P. reichenowi. The analysis of its genome should thus offer the opportunity to explore P. falciparum specific adaptations to humans. Our results bring new elements to the debate surrounding the origin of this lineage. They suggest that it may have been present in early hominoids and may have experienced a radiation congruent with that of its hosts. Our discovery highlights the paucity of our knowledge on the richness of Plasmodium species infecting primates and calls for more research in this direction.

A variant in the gene FUT9 is associated with susceptibility to placental malaria infection

Malaria in pregnancy forms a substantial part of the worldwide burden of malaria, with an estimated annual death toll of up to 200 000 infants, as well as increased maternal morbidity and mortality. Studies of genetic susceptibility to malaria have so far focused on infant malaria, with only a few studies investigating the genetic basis of placental malaria, focusing only on a limited number of candidate genes. The aim of this study therefore was to identify novel host genetic factors involved in placental malaria infection. To this end we carried out a nested case-control study on 180 Mozambican pregnant women with placental malaria infection, and 180 controls within an intervention trial of malaria prevention. We genotyped 880 SNPs in a set of 64 functionally related genes involved in glycosylation and innate immunity. A single nucleotide polymorphism (SNP) located in the gene FUT9, rs3811070, was significantly associated with placental malaria infection (odds ratio = 2.31, permutation P-value=0.028). Haplotypic analysis revealed a similarly strong association of a common haplotype of four SNPs including rs3811070. FUT9 codes for a fucosyl-transferase that is catalyzing the last step in the biosynthesis of the Lewis-x antigen, which forms part of the Lewis blood group-related antigens. These results therefore suggest an involvement of this antigen in the pathogenesis of placental malaria infection.

Prospective identification of malaria parasite genes under balancing selection

Background

Endemic human pathogens are subject to strong immune selection, and interrogation of pathogen genome variation for signatures of balancing selection can identify important target antigens. Several major antigen genes in the malaria parasite Plasmodium falciparum have shown such signatures in polymorphism-versus-divergence indices (comparing with the chimpanzee parasite P. reichenowi), and in allele frequency based indices.

Methodology/Principal Findings

To compare methods for prospective identification of genes under balancing selection, 26 additional genes known or predicted to encode surface-exposed proteins of the invasive blood stage merozoite were first sequenced from a panel of 14 independent P. falciparum cultured lines and P. reichenowi. Six genes at the positive extremes of one or both of the Hudson-Kreitman-Aguade (HKA) and McDonald-Kreitman (MK) indices were identified. Allele frequency based analysis was then performed on a Gambian P. falciparum population sample for these six genes and three others as controls. Tajima's D (TjD) index was most highly positive for the msp3/6-like PF10_0348 (TjD = 1.96) as well as the positive control ama1 antigen gene (TjD = 1.22). Across the genes there was a strong correlation between population TjD values and the relative HKA indices (whether derived from the population or the panel of cultured laboratory isolates), but no correlation with the MK indices.

Conclusions/Significance

Although few individual parasite genes show significant evidence of balancing selection, analysis of population genomic and comparative sequence data with the HKA and TjD indices should discriminate those that do, and thereby identify likely targets of immunity.

Molecular population genetics of SLC4A1 and Southeast Asian ovalocytosis

Southeast Asian ovalocytosis (SAO) is an erythrocyte abnormality that protects affected individuals from cerebral malaria. This trait is caused by a 27-bp deletion in the SLC4A1 gene, which is lethal when homozygous. We reseqeunced approximately 5 kb of SLC4A1 in an Indonesian population where SAO is prevalent to better understand the evolution of this clinically important trait. The four SAO chromosomes we resequenced share a single haplotype that differs from a sampled non-SAO haplotype only by the 27-bp deletion. Comparison of Indonesian sequence data to that from two other Asian populations (aboriginal Taiwanese and Japanese) shows Indonesian SLC4A1 to be strongly differentiated from the Taiwanese, but not the Japanese. Indeed, the Taiwanese sample contains only chromosomes that are highly divergent from all sampled SAO chromosomes. Because earlier studies have found an association between Austronesian-speakers (who most likely originated in Taiwan) and SAO, our failure to find SAO-like chromosomes in Taiwan is unexpected. Finally, our data find a strong excess of high-frequency derived alleles in all three populations. These alleles include the non-synonymous 'Memphis' variant, which is known to affect anion transport across the erythrocyte membrane. Our data suggest a role for recent natural selection acting on Memphis or a linked variant.

Pleistocene glacial refugia across the Appalachian Mountains and coastal plain in the millipede genus Narceus: evidence from population genetic, phylogeographic, and paleoclimatic data

Background

Species that are widespread throughout historically glaciated and currently non-glaciated areas provide excellent opportunities to investigate the role of Pleistocene climatic change on the distribution of North American biodiversity. Many studies indicate that northern animal populations exhibit low levels of genetic diversity over geographically widespread areas whereas southern populations exhibit relatively high levels. Recently, paleoclimatic data have been combined with niche-based distribution modeling to locate possible refugia during the Last Glacial Maximum. Using phylogeographic, population, and paleoclimatic data, we show that the distribution and mitochondrial data for the millipede genus Narceus are consistent with classical examples of Pleistocene refugia and subsequent post-glacial population expansion seen in other organismal groups.

Results

The phylogeographic structure of Narceus reveals a complex evolutionary history with signatures of multiple refugia in southeastern North America followed by two major northern expansions. Evidence for refugial populations were found in the southern Appalachian Mountains and in the coastal plain. The northern expansions appear to have radiated from two separate refugia, one from the Gulf Coastal Plain area and the other from the mid-Atlantic coastal region. Distributional models of Narceus during the Last Glacial Maximum show a dramatic reduction from the current distribution, with suitable ecological zones concentrated along the Gulf and Atlantic coastal plain. We found a strong correlation between these zones of ecological suitability inferred from our paleo-model with levels of genetic diversity derived from phylogenetic and population estimates of genetic structuring.

Conclusion

The signature of climatic change, during and after the Pleistocene, on the distribution of the millipede genus Narceus is evident in the genetic data presented. Niche-based historical distribution modeling strengthens the conclusions drawn from the genetic data and proves useful in identifying probable refugia. Such interdisciplinary biogeographic studies provide a comprehensive approach to understanding these processes that generate and maintain biodiversity as well as the framework necessary to explore questions regarding evolutionary diversification of taxa.

Population structure of Toxoplasma gondii: clonal expansion driven by infrequent recombination and selective sweeps

Toxoplasma gondii is among the most successful parasites. It is capable of infecting all warm-blooded animals and causing opportunistic disease in humans. T. gondii has a striking clonal population structure consisting of three predominant lineages in North America and Europe. Clonality is associated with the recent emergence of a monomorphic version of Chr1a, which drove a selective genetic sweep within the past 10,000 years. Strains from South America diverged from those in North America some 1-2 mya; recently, however, the monomorphic Chr1a has extended into regions of South America, where it is also associated with clonality. The recent spread of a few dominant lineages has dramatically shaped the population structure of T. gondii and has resulted in most lineages sharing a highly pathogenic nature. Understanding the factors that have shaped the population structure of T. gondii has implications for the emergence and transmission of human pathogens.

African genetic diversity: implications for human demographic history, modern human origins, and complex disease mapping

Comparative studies of ethnically diverse human populations, particularly in Africa, are important for reconstructing human evolutionary history and for understanding the genetic basis of phenotypic adaptation and complex disease. African populations are characterized by greater levels of genetic diversity, extensive population substructure, and less linkage disequilibrium (LD) among loci compared to non-African populations. Africans also possess a number of genetic adaptations that have evolved in response to diverse climates and diets, as well as exposure to infectious disease. This review summarizes patterns and the evolutionary origins of genetic diversity present in African populations, as well as their implications for the mapping of complex traits, including disease susceptibility.

Comparative evolutionary genomics of human malaria parasites

The parasites Plasmodium falciparum and Plasmodium vivax are responsible for the majority of human malaria cases worldwide. Despite many similarities in their biology, they frequently are studied in isolation. With the completion of the P. vivax genome and the generation of an initial P. falciparum genetic diversity map, attempts are being made to infer inter- and intra-species genome evolution. Here, we briefly review our current knowledge of comparative evolutionary genomics of the two species in the light of several presentations at the Molecular Approaches to Malaria 2008 meeting in Lorne, Australia and ask the question: can evolutionary genomics of one species inform the other?

Malaria research in the post-genomic era

For many pathogens the availability of genome sequence, permitting genome-dependent methods of research, can partially substitute for powerful forward genetic methods (genome-independent) that have advanced model organism research for decades. In 2002 the genome sequence of Plasmodium falciparum, the parasite causing the most severe type of human malaria, was completed, eliminating many of the barriers to performing state-of-the-art molecular biological research on malaria parasites. Although new, licensed therapies may not yet have resulted from genome-dependent experiments, they have produced a wealth of new observations about the basic biology of malaria parasites, and it is likely that these will eventually lead to new therapeutic approaches. This review will focus on the basic research discoveries that have depended, in part, on the availability of the Plasmodium genome sequences.

Big bang in the evolution of extant malaria parasites

Malaria parasites (genus Plasmodium) infect all classes of terrestrial vertebrates and display host specificity in their infections. It is therefore assumed that malaria parasites coevolved intimately with their hosts. Here, we propose a novel scenario of malaria parasite-host coevolution. A phylogenetic tree constructed using the malaria parasite mitochondrial genome reveals that the extant primate, rodent, bird, and reptile parasite lineages rapidly diverged from a common ancestor during an evolutionary short time period. This rapid diversification occurred long after the establishment of the primate, rodent, bird, and reptile host lineages, which implies that host-switch events contributed to the rapid diversification of extant malaria parasite lineages. Interestingly, the rapid diversification coincides with the radiation of the mammalian genera, suggesting that adaptive radiation to new mammalian hosts triggered the rapid diversification of extant malaria parasite lineages.

Host associations and evolutionary relationships of avian blood parasites from West Africa

The host specificity of blood parasites recovered from a survey of 527 birds in Cameroon and Gabon was examined at several levels within an evolutionary framework. Unique mitochondrial lineages of Haemoproteus were recovered from an average of 1.3 host species (maximum=3) and 1.2 host families (maximum=3) while lineages of Plasmodium were recovered from an average of 2.5 species (maximum=27) and 1.6 families (maximum=9). Averaged within genera, lineages of both Plasmodium and Haemoproteus were constrained in their host distribution relative to random expectations. However, while several individual lineages within both genera exhibited significant host constraint, host breadth varied widely among related lineages, particularly within the genus Plasmodium. Several lineages of Plasmodium exhibited extreme generalist host-parasitism strategies while other lineages appeared to have been constrained to certain host families over recent evolutionary history. Sequence data from two nuclear genes recovered from a limited sample of Plasmodium parasites indicated that, at the resolution of this study, inferences regarding host breadth were unlikely to be grossly affected by the use of parasite mitochondrial lineages as a proxy for biological species. The use of divergent host-parasitism strategies among closely related parasite lineages suggests that host range is a relatively labile character. Since host specificity may also influence parasite virulence, these results argue for considering the impact of haematozoa on avian hosts on a lineage-specific basis.

Complex population genetic and demographic history of the Salangid, Neosalanx taihuensis, based on cytochrome b sequences

BACKGROUND

The Salangid icefish Neosalanx taihuensis (Salangidae) is an economically important fish, which is endemic to China, restricted to large freshwater systems (e.g. lakes, large rivers and estuaries) and typically exhibit low vagility. The continuous distribution ranges from the temperate region of the Huai and Yellow River basins to the subtropical region of the Pearl River basin. This wide ranging distribution makes the species an ideal model for the study of palaeoclimatic effects on population genetic structure and phylogeography. Here, we aim to analyze population genetic differentiation within and between river basins and demographic history in order to understand how this species responded to severe climatic oscillations, decline of the sea levels during the Pleistocene ice ages and tectonic activity.

RESULTS

We obtained the complete mtDNA cytochrome b sequences (1141 bp) of 354 individuals from 13 populations in the Pearl River, the Yangze River and the Huai River basin. Thirty-six haplotypes were detected. Haplotype frequency distributions were strongly skewed, with most haplotypes (n = 24) represented only in single samples each and thus restricted to a single population. The most common haplotype (H36) was found in 49.15% of all individuals. Analysis of molecular variance (AMOVA) revealed a random pattern in the distribution of genetic diversity, which is inconsistent with contemporary hydrological structure. Significant levels of genetic subdivision were detected among populations within basins rather than between the three basins. Demographic analysis revealed that the population size in the Pearl River basin has remained relatively constant whereas the populations in the Yangze River and the Huai River basins expanded about 221 and 190 kyr ago, respectively, with the majority of mutations occurring after the last glacial maximum (LGM).

CONCLUSION

The observed complex genetic pattern of N. taihuensis is coherent with a scenario of multiple unrelated founding events by long-distance colonization and dispersal combined with contiguous population expansion and locally restricted gene flow. We also found that this species was likely severely impacted by past glaciations. More favourable climate and the formation of large suitable habitations together facilitated population expansion after the late Quaternary (especially the LGM). We proposed that all populations should be managed and conserved separately, especially for habitat protection.

Evidence of gene conversion in the evolutionary process of the codon 41/42 (-CTTT) mutation causing beta-thalassemia in southern China

The 4-bp deletion (-CTTT) at codon 41/42 (CD41/42) of the human beta-globin gene represents one of the most common beta-thalassemia mutations in East Asia and Southeast Asia, which is historically afflicted with endemic malaria, thus hypothetically evolving under natural selection by malaria infection. To understand the evolutionary process of generating the beta(CD41/42) allele and its maintenance, including the effect of natural selection on the pattern of linkage disequilibrium (LD), we sequenced a 15.933-kb region spanning 20.693 kb of the beta-globin cluster surrounding the 4-bp deletion using a sample from a Chinese population consisting of 24 normal individuals and 16 heterozygotes for the deletion. Forty-nine polymorphic sites were found. Analysis of the data, using a variety of methods including formal population genetics analysis and visual approaches, suggests that the spread of the CD41/42 (-CTTT) deletion is most likely mediated by interallelic gene conversion, although independent deletions in different lineages are also possible. The neutrality test resulted in a significant positive Tajima's D for the beta-globin locus, which is consistent with the existence of balancing selection. This suggests that the 4-bp deletion that occurred at this locus may be an event that is subject to natural selection, due to malaria, which leads to the heterozygote advantage, spread widely with further help by conversion and migration. The evolutionary process of this mutant through gene conversion that could conceivably take place between the 4-bp deletion and the normal sequence in the respective region is discussed in detail.

On the origin and spread of the Scab disease of apple: out of central Asia

BACKGROUND

Venturia inaequalis is an ascomycete fungus responsible for apple scab, a disease that has invaded almost all apple growing regions worldwide, with the corresponding adverse effects on apple production. Monitoring and predicting the effectiveness of intervention strategies require knowledge of the origin, introduction pathways, and population biology of pathogen populations. Analysis of the variation of genetic markers using the inferential framework of population genetics offers the potential to retrieve this information.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we present a population genetic analysis of microsatellite variation in 1,273 strains of V. inaequalis representing 28 orchard samples from seven regions in five continents. Analysis of molecular variance revealed that most of the variation (88%) was distributed within localities, which is consistent with extensive historical migrations of the fungus among and within regions. Despite this shallow population structure, clustering analyses partitioned the data set into separate groups corresponding roughly to geography, indicating that each region hosts a distinct population of the fungus. Comparison of the levels of variability among populations, along with coalescent analyses of migration models and estimates of genetic distances, was consistent with a scenario in which the fungus emerged in Central Asia, where apple was domesticated, before its introduction into Europe and, more recently, into other continents with the expansion of apple growing. Across the novel range, levels of variability pointed to multiple introductions and all populations displayed signatures of significant post-introduction increases in population size. Most populations exhibited high genotypic diversity and random association of alleles across loci, indicating recombination both in native and introduced areas.

CONCLUSIONS/SIGNIFICANCE

Venturia inaequalis is a model of invasive phytopathogenic fungus that has now reached the ultimate stage of the invasion process with a broad geographic distribution and well-established populations displaying high genetic variability, regular sexual reproduction, and demographic expansion.

Sequence analysis of Plasmodium falciparum cytochrome b in multiple geographic sites

Background

The antimalarial drug atovaquone specifically targets Plasmodium falciparum cytochrome b (Pfcytb), a mitochondrial gene with uniparental inheritance. Cases of resistance to atovaquone associated with mutant Pfcytb have been reported, justifying efforts to better document the natural polymorphism of this gene. To this end, a large molecular survey was conducted in several malaria endemic areas where atovaquone was not yet in regular use.

Methods

The polymorphism of the Pfcytb was analysed by direct sequencing of PCR products corresponding to the full length coding region. Sequence was generated for 671 isolates originating from three continents: Africa (Senegal, Ivory Coast, Central African Republic and Madagascar), Asia (Cambodia) and South America (French Guiana).

Results

Overall, 11 polymorphic sites were observed, of which eight were novel mutations. There was a large disparity in the geographic distribution of the mutants. All isolates from Senegal, Central African Republic and Madagascar displayed a Camp/3D7 wild type Pfcytb sequence, as did most samples originating from Cambodia and Ivory Coast. One synonymous (t759a at codon V253V) and two non-synonymous (t553g and a581g at codons F185V and H194R, respectively) singletons were detected in Ivory Coast. Likewise, two synonymous (a126t and c793t at codons -T42T and L265L, respectively) singletons were observed in Cambodia. In contrast, seven mutated sites, affecting seven codons and defining four mutant haplotypes were observed in French Guiana. The wild type allele was observed in only 14% of the French Guiana isolates. The synonymous c688t mutation at position L230L was highly prevalent; the most frequent allele was the c688t single mutant, observed in 84% of the isolates. The other alleles were singletons (a126t/a165c, a4g/a20t/a1024c and a20t/t341c/c688t corresponding to T42T/S55S, N2D/N71I/I342L, N71I/L114S/L230L, respectively" please replace with ' corresponding to T42T/S55S, N2D/N71I/I342L and N71I/L114S/L230L, respectively). The codon 268 polymorphisms associated with atovaquone resistance were not observed in the panel the isolates studied. Overall, the wild type PfCYTb protein isoform was highly predominant in all study areas, including French Guiana, suggesting stringent functional constraints.

Conclusion

These data along with previously identified Pfcytb field polymorphisms indicate a clustering of molecular signatures, suggesting different ancestral types in South America and other continents. The absence of mutations associated with most atovaquone-proguanil clinical failures indicates that the atovaquone-proguanil association is an interesting treatment option in the study areas.

A quantitative synthesis of the medicinal ethnobotany of the Malinké of Mali and the Asháninka of Peru, with a new theoretical framework

BACKGROUND

Although ethnomedically and taxonomically guided searches for new medicinal plants can improve the percentage of plants found containing active compounds when compared to random sampling, ethnobotany has fulfilled little of its promise in the last few decades to deliver a bounty of new, laboratory-proven medicinal plants and compounds. It is quite difficult to test, isolate, and elucidate the structure and mechanism of compounds from the plethora of new medicinal plant uses described each year with limited laboratory time and resources and the high cost of clinical trials of new drug candidates.

METHODS

A new quantitative theoretical framework of mathematical formulas called "relational efficacy" is proposed that should narrow down this search for new plant-derived medicines based on the hypothesis that closely related plants used to treat closely related diseases in distantly related cultures have a higher probability of being effective because they are more likely to be independent discoveries of similar plant compounds and disease mechanisms. A prerequisite to this hypothesis, the idea that empirical testing in traditional medicine will lead to choosing similar medicinal plants and therefore the medicinal flora of two distant cultures will prove to be more similar than their general flora, is tested using resampling statistics on cross-cultural field data of the plants used by the Malinké of Mali and the Asháninka of Peru to treat the diseases malaria, African sleeping sickness, Chagas' disease, leishmaniasis, diabetes, eczema, asthma, and uterine fibroids.

RESULTS

In this case, the similarity of the medicinal floras is found to be significantly greater than the similarity of the general floras, but only when the diseases in question are grouped into the categories of parasitic and autoimmune diseases.

CONCLUSION

If the central theoretical framework of this hypothesis is shown to be true, it will allow the synthesis of medicinal plant information from around the world to pinpoint the species with the highest potential efficacy to take into the laboratory and analyze further, ultimately saving much field and laboratory time and resources.

Evidence of non-neutral polymorphism in Plasmodium falciparum gamete surface protein genes Pfs47 and Pfs48/45

Targeted disruption of particular members of the Plasmodium 6-cys protein gene family, including Ps47, Ps48/45 and Ps230, is known to dramatically affect parasite fertility. Because loci critical to fertility in many eukaryote species have been shown to be under strong positive selection, we examined sequence variation in four members of the 6-cys protein gene family in Plasmodium falciparum (Pfs36, Pfs38, Pfs47 and Pfs48/45) to determine whether genetic variation in these loci may be of functional significance. Sequence polymorphism among 11 laboratory isolates of P. falciparum was compared with divergence from the respective orthologues in the closely related species P. reichenowi, showing an almost significant skew towards within-species non-synonymous polymorphism in Pfs47 and Pfs48/45 (by the McDonald-Kreitman test) but clearly non-significant results for Pfs36 and Pfs38. A preliminary analysis of Pfs47 sequence polymorphism in field isolates of P. falciparum showed exceptionally high fixation indices (F(ST)) among geographically distinct populations, similar to results seen previously for Pfs48/45. Therefore, both Pfs47 and Pfs48/45 were further analysed by sequencing polymorphic parts of the genes from a Tanzanian population sample of oocysts (a means of analysing diploid genotypes). Both genes displayed higher inbreeding coefficients (F(IS)) compared with the average of 11 unlinked microsatellite loci. These results suggest that allelic variation in these two genes may be functionally significant in influencing mating interactions, a hypothesis that could be tested by fertilization experiments with targeted allelic replacement.

Genetic linkage and association analyses for trait mapping in Plasmodium falciparum

Genetic studies of Plasmodium falciparum laboratory crosses and field isolates have produced valuable insights into determinants of drug responses, antigenic variation, disease virulence, cellular development and population structures of these virulent human malaria parasites. Full-genome sequences and high-resolution haplotype maps of SNPs and microsatellites are now available for all 14 parasite chromosomes. Rapidly increasing genetic and genomic information on Plasmodium parasites, mosquitoes and humans will combine as a rich resource for new advances in our understanding of malaria, its transmission and its manifestations of disease.

Global phylogeographic limits of Hawaii's avian malaria

The introduction of avian malaria (Plasmodium relictum) to Hawaii has provided a model system for studying the influence of exotic disease on naive host populations. Little is known, however, about the origin or the genetic variation of Hawaii's malaria and traditional classification methods have confounded attempts to place the parasite within a global ecological and evolutionary context. Using fragments of the parasite mitochondrial gene cytochrome b and the nuclear gene dihydrofolate reductase-thymidylate synthase obtained from a global survey of greater than 13000 avian samples, we show that Hawaii's avian malaria, which can cause high mortality and is a major limiting factor for many species of native passerines, represents just one of the numerous lineages composing the morphological parasite species. The single parasite lineage detected in Hawaii exhibits a broad host distribution worldwide and is dominant on several other remote oceanic islands, including Bermuda and Moorea, French Polynesia. The rarity of this lineage in the continental New World and the restriction of closely related lineages to the Old World suggest limitations to the transmission of reproductively isolated parasite groups within the morphological species.

Population genomics of the immune evasion (var) genes of Plasmodium falciparum

Var genes encode the major surface antigen (PfEMP1) of the blood stages of the human malaria parasite Plasmodium falciparum. Differential expression of up to 60 diverse var genes in each parasite genome underlies immune evasion. We compared the diversity of the DBLα domain of var genes sampled from 30 parasite isolates from a malaria endemic area of Papua New Guinea (PNG) and 59 from widespread geographic origins (global). Overall, we obtained over 8,000 quality-controlled DBLα sequences. Within our sampling frame, the global population had a total of 895 distinct DBLα “types” and negligible overlap among repertoires. This indicated that var gene diversity on a global scale is so immense that many genomes would need to be sequenced to capture its true extent. In contrast, we found a much lower diversity in PNG of 185 DBLα types, with an average of approximately 7% overlap among repertoires. While we identify marked geographic structuring, nearly 40% of types identified in PNG were also found in samples from different countries showing a cosmopolitan distribution for much of the diversity. We also present evidence to suggest that recombination plays a key role in maintaining the unprecedented levels of polymorphism found in these immune evasion genes. This population genomic framework provides a cost effective molecular epidemiological tool to rapidly explore the geographic diversity of var genes.

Malaria parasites live in red blood cells of the human host for part of the life cycle, during which a family of diverse antigens known as PfEMP1 are placed on the surface. PfEMP1 variants switch by sequential expression of up to 60 var genes. This allows the parasite to evade immune detection within an individual host, enhancing its chances to be transmitted to the mosquito vector in situations where mosquitoes are seasonally available. Methods to rapidly assess var gene diversity in parasite populations are needed to measure antigenic diversity and define relationships with malaria transmission. Using a specialized framework, we completed the first systematic sampling of var genes from parasite genomes obtained from the same (Papua New Guinea [PNG]) and different (global) populations. Globally, there was no limit to the number of var genes because parasites rarely shared var genes. In PNG, var gene numbers were restricted due to high levels of sharing, and most were only found in that population. Recombination was important to the evolution of var genes in PNG. The data suggest there are distinct var genes in different populations, which may have consequences for the spread of malarial disease from one geographic area to another.

Bayesian analysis of new and old malaria parasite DNA sequence data demonstrates the need for more phylogenetic signal to clarify the descent of Plasmodium falciparum

Molecular systematic studies published during the last 15 years to clarify the phylogenetic relationships among the malaria parasites have led to two major hypotheses on the descent of Plasmodium falciparum: One supports an avian origin as a result of a relatively recent host switch, and another one favours the evolutionary development of P. falciparum together with its human host from primate ancestors. In this paper, we present phylogenetic analyses of three different Plasmodium genes, the nuclear 18 small sub-unit (SSU) ribosomal ribonucleic acid (rRNA), the mitochondrial cytochrome b (cyt b) and the plastid caseinolytic protease C (ClpC) gene, using numerous haemosporidian parasite DNA sequences obtained from the GenBank as well as several new sequences for major malaria parasites including the avian one Plasmodium cathemerium, which has never been considered in molecular phylogenetic analyses before. Most modern and sophisticated DNA substitution models based on Bayesian inference analysis were applied to estimate the cyt b and ClpC phylogenetic trees, whereas the 18 SSU rRNA gene was examined with regards to its secondary structure using PHASE software. Our results indicate that the data presently available are generally neither sufficient in number nor in information to solve the problem of the phylogenetic origin of P. falciparum.

Molecular epidemiology of malaria

Malaria persists as an undiminished global problem, but the resources available to address it have increased. Many tools for understanding its biology and epidemiology are well developed, with a particular richness of comparative genome sequences. Targeted genetic manipulation is now effectively combined with in vitro culture assays on the most important human parasite, Plasmodium falciparum, and with in vivo analysis of rodent and monkey malaria parasites in their laboratory hosts. Studies of the epidemiology, prevention, and treatment of human malaria have already been influenced by the availability of molecular methods, and analyses of parasite polymorphisms have long had useful and highly informative applications. However, the molecular epidemiology of malaria is currently undergoing its most substantial revolution as a result of the genomic information and technologies that are available in well-resourced centers. It is a challenge for research agendas to face the real needs presented by a disease that largely exists in extremely resource-poor settings, but it is one that there appears to be an increased willingness to undertake. To this end, developments in the molecular epidemiology of malaria are reviewed here, emphasizing aspects that may be current and future priorities.

Blood parasites of chickens in Uganda and Cameroon with molecular descriptions of Leucocytozoon schoutedeni and Trypanosoma gallinarum

Using microscopy and PCR, we determined the prevalence of blood parasites in village chickens in Uganda and Cameroon. Of 148 individuals tested, 18.3% were infected with Leucocytozoon schoutedeni (Haemosporida, Leucocytozoidae) and 4.1% were infected with Trypanosoma gallinarum (Kinetoplastida, Trypanosomatidae). No other blood parasites were detected. Subsequent phylogenetic analysis of the cytochrome b gene of L. schoutedeni identified 2 distinct lineages that were found at all 3 sampling locations in Uganda. The sequence divergence between these 2 lineages is 1.5%. One of these lineages was also found in chickens in Cameroon, nearly 2,000 km distant. There are no morphological differences between blood stages of the parasites represented by the 2 different lineages, suggesting that cytochrome b gene sequence divergence can be as high as 1.5% within a single well-defined morphospecies of Leucocytozoon. We sequenced a portion of the small subunit ribosomal RNA gene (SSU rRNA) of T. gallinarum, and redescribe T. gallinarum for the first time since its discovery in 1911. These are the first assignments of DNA sequence data to these morphospecies of Leucocytozoon and Trypanosoma and may represent an example of intraspecific sequence divergence.

A genome-wide map of diversity in Plasmodium falciparum

Genetic variation allows the malaria parasite Plasmodium falciparum to overcome chemotherapeutic agents, vaccines and vector control strategies and remain a leading cause of global morbidity and mortality. Here we describe an initial survey of genetic variation across the P. falciparum genome. We performed extensive sequencing of 16 geographically diverse parasites and identified 46,937 SNPs, demonstrating rich diversity among P. falciparum parasites (pi = 1.16 x 10(-3)) and strong correlation with gene function. We identified multiple regions with signatures of selective sweeps in drug-resistant parasites, including a previously unidentified 160-kb region with extremely low polymorphism in pyrimethamine-resistant parasites. We further characterized 54 worldwide isolates by genotyping SNPs across 20 genomic regions. These data begin to define population structure among African, Asian and American groups and illustrate the degree of linkage disequilibrium, which extends over relatively short distances in African parasites but over longer distances in Asian parasites. We provide an initial map of genetic diversity in P. falciparum and demonstrate its potential utility in identifying genes subject to recent natural selection and in understanding the population genetics of this parasite.

The origin and age of Plasmodium vivax

The evolutionary history of Plasmodium vivax has recently been addressed in terms of its origin as a parasite of humans and the age of extant populations. The consensus is that P. vivax originated as a result of a host switch from a non-human primate to hominids and that the extant populations did not originate as recently as previously proposed. Here, we show that, in a comparison of parasite isolates from across the world, Asian populations of P. vivax are the oldest. We discuss how this result, together with the phylogenetic evidence that P. vivax derived from Plasmodium found in Southeast Asian macaques, is most simply explained by assuming an Asian origin of this parasite. Nevertheless, the available data show only the tip of the iceberg. We discuss how sampling might affect time estimates to the most recent common ancestor for P. vivax populations and suggest that spatially explicit estimates are needed to understand the demographic history of this parasite better.

AA. A comparative study of the genetic diversity of the 42kDa fragment of the merozoite surface protein 1 in Plasmodium falciparum and P. vivax

We investigated the genetic diversity of the 42kDa fragment of the merozoite surface protein 1 (MSP-1) antigen in Plasmodium falciparum and P. vivax, as well as in non-human primate malarial parasites. This fragment undergoes a proteolytic cleavage generating two fragments of 19kDa (MSP-1(19)) and 33kDa (MSP-1(33)) that are critical in erythrocyte invasion. We found that overall the MSP-1(33) fragment exhibits greater genetic diversity than the MSP-1(19) regardless of the species. We have found evidence for positive natural selection only in the human malaria parasites by comparing the rate of non-synonymous versus synonymous substitutions. In addition, we found clear differences between the two major human malaria parasites. In the case of P. falciparum, positive natural selection is acting on the MSP-1(19) region while the MSP-1(33) is neutral or under purifying selection. The opposite pattern was observed in P. vivax. Our results suggest different roles of this antigen in the host-parasite immune interaction in each of the major human malarial parasites.

Differential evidence of natural selection on two leading sporozoite stage malaria vaccine candidate antigens

Experimental malaria vaccines based on two sporozoite stage candidate antigens of Plasmodium falciparum, the circumsporozoite protein (CSP) and thrombospondin-related adhesive protein (TRAP), have undergone clinical trials of efficacy. The relevance of naturally existing polymorphism in these molecules remains unknown. Sequence polymorphism in the genes encoding these antigens was studied in a Gambian population (sample of 48 trap and 44 csp gene sequences) to test for signatures of selection that would result from naturally acquired immunity. Allele frequency distributions were analyzed and compared with data from another population (in Thailand). Patterns of non-synonymous and synonymous polymorphism in P. falciparum and in Plasmodium vivax were compared with divergence from related species. Results indicate that polymorphism in TRAP is under strong selection for amino acid sequence diversity and that allele frequencies are under balancing selection within the Gambian P. falciparum population. There was no such evidence for CSP, calling into question the idea that most polymorphisms in this gene are under immune selection. There was a weak trend for regions known to encode T cell epitopes to have slightly higher indices suggesting balancing selection. Overall, the results predict more allele-specific immunity to TRAP than to CSP and should be considered in design and efficacy testing of vaccine candidates based on these antigens.

Haplotype trees and modern human origins

A haplotype is a multisite haploid genotype at two or more polymorphic sites on the same chromosome in a defined DNA region. An evolutionary tree of the haplotypes can be estimated if the DNA region had little to no recombination. Haplotype trees can be used to reconstruct past human gene-flow patterns and historical events, but any single tree captures only a small portion of evolutionary history, and is subject to error. A fuller view of human evolution requires multiple DNA regions, and errors can be minimized by cross-validating inferences across loci. An analysis of 25 DNA regions reveals an out-of-Africa expansion event at 1.9 million years ago. Gene flow with isolation by distance was established between African and Eurasian populations by about 1.5 million years ago, with no detectable interruptions since. A second out-of-Africa expansion occurred about 700,000 years ago, and involved interbreeding with at least some Eurasian populations. A third out-of-Africa event occurred around 100,000 years ago, and was also characterized by interbreeding, with the hypothesis of a total Eurasian replacement strongly rejected (P < 10(-17)). This does not preclude the possibility that some Eurasian populations could have been replaced, and the status of Neanderthals is indecisive. Demographic inferences from haplotype trees have been inconsistent, so few definitive conclusions can be made at this time. Haplotype trees from human parasites offer additional insights into human evolution and raise the possibility of an Asian isolate of humanity, but once again not in a definitive fashion. Haplotype trees can also indicate which genes were subject to positive selection in the lineage leading to modern humans. Genetics provides many insights into human evolution, but those insights need to be integrated with fossil and archaeological data to yield a fuller picture of the origin of modern humans.

Variable SNP density in aspartyl-protease genes of the malaria parasite Plasmodium falciparum

An analysis of the diversity of the aspartyl proteases of Plasmodium falciparum, known as plasmepsins (PMs), was completed in view of their possible role as drug targets. DNA sequence polymorphisms were identified in nine pm genes including their non-coding (introns and 5' flanking) sequences. All genes contained at least one single nucleotide polymorphism (SNP). Extensive microsatellite diversity was observed predominantly in non-coding sequences. All but one non-synonymous polymorphism (a conservative substitution) were mapped to the surface of the predicted protein, contradicting a possible role in enzymatic activity. The distribution of SNPs was found to be non-random among pm genes, with pm6 and pm10 having significantly higher SNP densities, suggesting they were under selection. For pm6 the majority of the SNPs were in introns and some of these may contribute to splice site variation. SNPs were found at a high density in both the coding and non-coding sequences of pm10. Recombination was important in generating additional diversity at this locus. Although direct selection for pm10 mutations could not be ruled out, the presence of balancing selection and a high density of SNPs in non-coding sequence led us to propose that another gene under selection may be influencing the diversity in the region. By sequencing short DNA tags in a 200 kb region flanking pm10 we show that a cluster of antigen genes, known to be under diversifying selection, may contribute to the observed diversity. We discuss the importance of diversity and local selection effects when choosing drug targets for intervention strategies.

Contrasting signatures of selection on the Plasmodium falciparum erythrocyte binding antigen gene family

Erythrocyte binding antigens of Plasmodium falciparum are involved in erythrocyte invasion, and may be targets of acquired immunity. Of the five eba genes, protein products have been detected for eba-175, eba-181 and eba-140, but not for psieba-165 or ebl-1, providing opportunity for comparative analysis of genetic variation to identify selection. Region II of each of these genes was sequenced from a cross-sectional sample of parasites in an endemic Kenyan population, and the frequency distributions of polymorphisms analysed. A positive value of Tajima's D was observed for eba-175 (D=1.13) indicating an excess of intermediate frequency polymorphisms, while all other genes had negative values, the most negative being ebl-1 (D=-2.35) followed by psieba-165 (D=-1.79). The eba-175 and ebl-1 genes were then studied in a sample of parasites from Thailand, for which a positive Tajima's D value was again observed for eba-175 (D=1.79), and a negative value for ebl-1 (D=-1.85). This indicates that eba-175 is under balancing selection in each population, in strong contrast to the other members of the gene family, particularly ebl-1 and psieba-165 that may have been under recent directional selection. Population expansion simulations were performed under a neutral model, further supporting the departures from neutrality of these genes.

Rangewide population genetic structure of the African malaria vector Anopheles funestus

Anopheles funestus is a primary vector of malaria in Africa south of the Sahara. We assessed its rangewide population genetic structure based on samples from 11 countries, using 10 physically mapped microsatellite loci, two per autosome arm and the X (N = 548), and 834 bp of the mitochondrial ND5 gene (N = 470). On the basis of microsatellite allele frequencies, we found three subdivisions: eastern (coastal Tanzania, Malawi, Mozambique and Madagascar), western (Burkina Faso, Mali, Nigeria and western Kenya), and central (Gabon, coastal Angola). A. funestus from the southwest of Uganda had affinities to all three subdivisions. Mitochondrial DNA (mtDNA) corroborated this structure, although mtDNA gene trees showed less resolution. The eastern subdivision had significantly lower diversity, similar to the pattern found in the codistributed malaria vector Anopheles gambiae. This suggests that both species have responded to common geographic and/or climatic constraints. The western division showed signatures of population expansion encompassing Kenya west of the Rift Valley through Burkina Faso and Mali. This pattern also bears similarity to A. gambiae, and may reflect a common response to expanding human populations following the development of agriculture. Due to the presumed recent population expansion, the correlation between genetic and geographic distance was weak. Mitochondrial DNA revealed further cryptic subdivision in A. funestus, not detected in the nuclear genome. Mozambique and Madagascar samples contained two mtDNA lineages, designated clade I and clade II, that were separated by two fixed differences and an average of 2% divergence, which implies that they have evolved independently for approximately 1 million years. Clade I was found in all 11 locations, whereas clade II was sampled only on Madagascar and Mozambique. We suggest that the latter clade may represent mtDNA capture by A. funestus, resulting from historical gene flow either among previously isolated and divergent populations or with a related species.