A genome-wide map of diversity in Plasmodium falciparum

Molecular approaches for a better understanding of the epidemiology and population genetics ofLeishmania

Molecular approaches are being used increasingly for epidemiological studies of visceral and cutaneous leishmaniases. Several molecular markers resolving genetic differences betweenLeishmaniaparasites at species and strain levels have been developed to address key epidemiological and population genetic questions. The current gold standard, multilocus enzyme typing (MLEE), needs cultured parasites and lacks discriminatory power. PCR assays identifying species directly with clinical samples have proven useful in numerous field studies. Multilocus sequence typing (MLST) is potentially the most powerful phylogenetic approach and will, most probably, replace MLEE in the future. Multilocus microsatellite typing (MLMT) is able to discriminate below the zymodeme level and seems to be the best candidate for becoming the gold standard for distinction of strains. Population genetic studies by MLMT revealed geographical and hierarchic population structure inL. tropica, L. majorand theL. donovanicomplex. The existence of hybrids and gene flow betweenLeishmaniapopulations suggests that sexual recombination is more frequent than previously thought. However, typing and analytical tools need to be further improved. Accessible databases should be created and sustained for integrating data obtained by different researchers. This would allow for global analyses and help to avoid biases in analyses due to small sample sizes.

Whole-genome sequencing and microarray analysis of ex vivo Plasmodium vivax reveal selective pressure on putative drug resistance genes

Plasmodium vivax causes 25-40% of malaria cases worldwide, yet research on this human malaria parasite has been neglected. Nevertheless, the recent publication of the P. vivax reference genome now allows genomics and systems biology approaches to be applied to this pathogen. We show here that whole-genome analysis of the parasite can be achieved directly from ex vivo-isolated parasites, without the need for in vitro propagation. A single isolate of P. vivax obtained from a febrile patient with clinical malaria from Peru was subjected to whole-genome sequencing (30× coverage). This analysis revealed over 18,261 single-nucleotide polymorphisms (SNPs), 6,257 of which were further validated using a tiling microarray. Within core chromosomal genes we find that one SNP per every 985 bases of coding sequence distinguishes this recent Peruvian isolate, designated IQ07, from the reference Salvador I strain obtained in 1972. This full-genome sequence of an uncultured P. vivax isolate shows that the same regions with low numbers of aligned sequencing reads are also highly variable by genomic microarray analysis. Finally, we show that the genes containing the largest ratio of nonsynonymous-to-synonymous SNPs include two AP2 transcription factors and the P. vivax multidrug resistance-associated protein (PvMRP1), an ABC transporter shown to be associated with quinoline and antifolate tolerance in Plasmodium falciparum. This analysis provides a data set for comparative analysis with important potential for identifying markers for global parasite diversity and drug resistance mapping studies.

A novel multifunctional oligonucleotide microarray for Toxoplasma gondii

BACKGROUND

Microarrays are invaluable tools for genome interrogation, SNP detection, and expression analysis, among other applications. Such broad capabilities would be of value to many pathogen research communities, although the development and use of genome-scale microarrays is often a costly undertaking. Therefore, effective methods for reducing unnecessary probes while maintaining or expanding functionality would be relevant to many investigators.

RESULTS

Taking advantage of available genome sequences and annotation for Toxoplasma gondii (a pathogenic parasite responsible for illness in immunocompromised individuals) and Plasmodium falciparum (a related parasite responsible for severe human malaria), we designed a single oligonucleotide microarray capable of supporting a wide range of applications at relatively low cost, including genome-wide expression profiling for Toxoplasma, and single-nucleotide polymorphism (SNP)-based genotyping of both T. gondii and P. falciparum. Expression profiling of the three clonotypic lineages dominating T. gondii populations in North America and Europe provides a first comprehensive view of the parasite transcriptome, revealing that ~49% of all annotated genes are expressed in parasite tachyzoites (the acutely lytic stage responsible for pathogenesis) and 26% of genes are differentially expressed among strains. A novel design utilizing few probes provided high confidence genotyping, used here to resolve recombination points in the clonal progeny of sexual crosses. Recent sequencing of additional T. gondii isolates identifies >620 K new SNPs, including ~11 K that intersect with expression profiling probes, yielding additional markers for genotyping studies, and further validating the utility of a combined expression profiling/genotyping array design. Additional applications facilitating SNP and transcript discovery, alternative statistical methods for quantifying gene expression, etc. are also pursued at pilot scale to inform future array designs.

CONCLUSIONS

In addition to providing an initial global view of the T. gondii transcriptome across major lineages and permitting detailed resolution of recombination points in a historical sexual cross, the multifunctional nature of this array also allowed opportunities to exploit probes for purposes beyond their intended use, enhancing analyses. This array is in widespread use by the T. gondii research community, and several aspects of the design strategy are likely to be useful for other pathogens.

Allele frequency-based and polymorphism-versus-divergence indices of balancing selection in a new filtered set of polymorphic genes in Plasmodium falciparum

Signatures of balancing selection operating on specific gene loci in endemic pathogens can identify candidate targets of naturally acquired immunity. In malaria parasites, several leading vaccine candidates convincingly show such signatures when subjected to several tests of neutrality, but the discovery of new targets affected by selection to a similar extent has been slow. A small minority of all genes are under such selection, as indicated by a recent study of 26 Plasmodium falciparum merozoite-stage genes that were not previously prioritized as vaccine candidates, of which only one (locus PF10_0348) showed a strong signature. Therefore, to focus discovery efforts on genes that are polymorphic, we scanned all available shotgun genome sequence data from laboratory lines of P. falciparum and chose six loci with more than five single nucleotide polymorphisms per kilobase (including PF10_0348) for in-depth frequency-based analyses in a Kenyan population (allele sample sizes >50 for each locus) and comparison of Hudson-Kreitman-Aguade (HKA) ratios of population diversity (π) to interspecific divergence (K) from the chimpanzee parasite Plasmodium reichenowi. Three of these (the msp3/6-like genes PF10_0348 and PF10_0355 and the surf(4.1) gene PFD1160w) showed exceptionally high positive values of Tajima's D and Fu and Li's F indices and have the highest HKA ratios, indicating that they are under balancing selection and should be prioritized for studies of their protein products as candidate targets of immunity. Combined with earlier results, there is now strong evidence that high HKA ratio (as well as the frequency-independent ratio of Watterson's /K) is predictive of high values of Tajima's D. Thus, the former offers value for use in genome-wide screening when numbers of genome sequences within a species are low or in combination with Tajima's D as a 2D test on large population genomic samples.

Plasmodium immunomics

The Plasmodium parasite, the causative agent of malaria, is an excellent model for immunomic-based approaches to vaccine development. The Plasmodium parasite has a complex life cycle with multiple stages and stage-specific expression of ∼5300 putative proteins. No malaria vaccine has yet been licensed. Many believe that an effective vaccine will need to target several antigens and multiple stages, and will require the generation of both antibody and cellular immune responses. Vaccine efforts to date have been stage-specific and based on only a very limited number of proteins representing <0.5% of the genome. The recent availability of comprehensive genomic, proteomic and transcriptomic datasets from human and selected non-human primate and rodent malarias provide a foundation to exploit for vaccine development. This information can be mined to identify promising vaccine candidate antigens, by proteome-wide screening of antibody and T cell reactivity using specimens from individuals exposed to malaria and technology platforms such as protein arrays, high throughput protein production and epitope prediction algorithms. Such antigens could be incorporated into a rational vaccine development process that targets specific stages of the Plasmodium parasite life cycle with immune responses implicated in parasite elimination and control. Immunomic approaches which enable the selection of the best possible targets by prioritising antigens according to clinically relevant criteria may overcome the problem of poorly immunogenic, poorly protective vaccines that has plagued malaria vaccine developers for the past 25 years. Herein, current progress and perspectives regarding Plasmodium immunomics are reviewed.

Genome scanning of Amazonian Plasmodium falciparum shows subtelomeric instability and clindamycin-resistant parasites

Here, we fully characterize the genomes of 14 Plasmodium falciparum patient isolates taken recently from the Iquitos region using genome scanning, a microarray-based technique that delineates the majority of single-base changes, indels, and copy number variants distinguishing the coding regions of two clones. We show that the parasite population in the Peruvian Amazon bears a limited number of genotypes and low recombination frequencies. Despite the essentially clonal nature of some isolates, we see high frequencies of mutations in subtelomeric highly variable genes and internal var genes, indicating mutations arising during self-mating or mitotic replication. The data also reveal that one or two meioses separate different isolates, showing that P. falciparum clones isolated from different individuals in defined geographical regions could be useful in linkage analyses or quantitative trait locus studies. Through pairwise comparisons of different isolates we discovered point mutations in the apicoplast genome that are close to known mutations that confer clindamycin resistance in other species, but which were hitherto unknown in malaria parasites. Subsequent drug sensitivity testing revealed over 100-fold increase of clindamycin EC(50) in strains harboring one of these mutations. This evidence of clindamycin-resistant parasites in the Amazon suggests that a shift should be made in health policy away from quinine + clindamycin therapy for malaria in pregnant women and infants, and that the development of new lincosamide antibiotics for malaria should be reconsidered.

Spiroindolones, a potent compound class for the treatment of malaria

Recent reports of increased tolerance to artemisinin derivatives--the most recently adopted class of antimalarials--have prompted a need for new treatments. The spirotetrahydro-beta-carbolines, or spiroindolones, are potent drugs that kill the blood stages of Plasmodium falciparum and Plasmodium vivax clinical isolates at low nanomolar concentration. Spiroindolones rapidly inhibit protein synthesis in P. falciparum, an effect that is ablated in parasites bearing nonsynonymous mutations in the gene encoding the P-type cation-transporter ATPase4 (PfATP4). The optimized spiroindolone NITD609 shows pharmacokinetic properties compatible with once-daily oral dosing and has single-dose efficacy in a rodent malaria model.

Low-complexity regions in Plasmodium falciparum: missing links in the evolution of an extreme genome

Over the past decade, attempts to explain the unusual size and prevalence of low-complexity regions (LCRs) in the proteins of the human malaria parasite Plasmodium falciparum have used both neutral and adaptive models. This past research has offered conflicting explanations for LCR characteristics and their role in, and influence on, the evolution of genome structure. Here we show that P. falciparum LCRs (PfLCRs) are not a single phenomenon, but rather consist of at least three distinct types of sequence, and this heterogeneity is the source of the conflict in the literature. Using molecular and population genetics, we show that these families of PfLCRs are evolving by different mechanisms. One of these families, named here the HighGC family, is of particular interest because these LCRs act as recombination hotspots, both in genes under positive selection for high levels of diversity which can be created by recombination (antigens) and those likely to be evolving neutrally or under negative selection (metabolic enzymes). We discuss how the discovery of these distinct species of PfLCRs helps to resolve previous contradictory studies on LCRs in malaria and contributes to our understanding of the evolution of the of the parasite's unusual genome.

Variable numbers of tandem repeats in Plasmodium falciparum genes

Genome variation studies in Plasmodium falciparum have focused on SNPs and, more recently, large-scale copy number polymorphisms and ectopic rearrangements. Here, we examine another source of variation: variable number tandem repeats (VNTRs). Interspersed low complexity features, including the well-studied P. falciparum microsatellite sequences, are commonly classified as VNTRs; however, this study is focused on longer coding VNTR polymorphisms, a small class of copy number variations. Selection against frameshift mutation is a main constraint on tandem repeats (TRs) in coding regions, while limited propagation of TRs longer than 975 nt total length is a minor restriction in coding regions. Comparative analysis of three P. falciparum genomes reveals that more than 9% of all P. falciparum ORFs harbor VNTRs, much more than has been reported for any other species. Moreover, genotyping of VNTR loci in a drug-selected line, progeny of a genetic cross, and 334 field isolates demonstrates broad variability in these sequences. Functional enrichment analysis of ORFs harboring VNTRs identifies stress and DNA damage responses along with chromatin modification activities, suggesting an influence on genome mutability and functional variation. Analysis of the repeat units and their flanking regions in both P. falciparum and Plasmodium reichenowi sequences implicates a replication slippage mechanism in the generation of TRs from an initially unrepeated sequence. VNTRs can contribute to rapid adaptation by localized sequence duplication. They also can confound SNP-typing microarrays or mapping short-sequence reads and therefore must be accounted for in such analyses.

Functional diversification between two related Plasmodium falciparum merozoite invasion ligands is determined by changes in the cytoplasmic domain

The pathogenesis of Plasmodium falciparum depends on efficient invasion into host erythrocytes. Parasite ligands encoded by multi-gene families interact with erythrocyte receptors. P. falciparum reticulocyte binding protein homologues (PfRhs) are expressed at the apical surface of invasive merozoites and have divergent ectodomains that are postulated to bind different erythrocyte receptors. Variant expression of these paralogues results in the use of alternative invasion pathways. Two PfRh proteins, PfRh2a and PfRh2b, are identical for 2700 N-terminal amino acids and differ only in a C-terminal 500 amino acid region, which includes a unique ectodomain, transmembrane domain and cytoplasmic domain. Despite their similarity, PfRh2b is required for a well-defined invasion pathway while PfRh2a is not required or sufficient for this pathway. Mapping the genomic region encoding these proteins revealed a recombinogenic locus with PfRh2a and PfRh2b in a head-to-head orientation. We have generated viable PfRh2a/2b chimeric parasites to identify the regions required for alternative invasion pathway utilization. We find that the differential ability to use these pathways is conferred by the cytoplasmic domains of PfRh2a and PfRh2b, not the ectodomain or transmembrane regions. Our results highlight the importance of the cytoplasmic domain for functional diversification of a major adhesive ligand family in malaria parasites.

Uncertainty in mapping malaria epidemiology: implications for control

Malaria is a location-specific, dynamic infectious disease transmitted by mosquitoes to humans and is influenced by environmental, vector, parasite, and host factors. The principal purposes of malarial epidemiology are 1) to describe the malarial distribution in space and time along with the physical, biologic, and social etiologic factors and 2) to guide control objectives for either modeling impact or measuring progress of control tactics. Mapping malaria and many of its causative factors has been achieved on many different levels from global distribution to biologic quantitative trait localization in humans, parasites, and mosquitoes. Despite these important achievements, a large degree of uncertainty still exists on the annual burden of malarial cases. Accurate, sensitive detection and treatment of asymptomatic reservoirs important to infectious transmission are additional components necessary for future control measures. Presently spurred by the leadership and funding of Bill and Melinda Gates, the malarial community is developing and implementing plans for elimination of malaria. The challenge for malariologists is to digitally integrate and map epidemiologic factors and intervention measures in space and time to target effective, sustainable control alongside research efforts.

Plasmodium falciparum accompanied the human expansion out of Africa

Plasmodium falciparum is distributed throughout the tropics and is responsible for an estimated 230 million cases of malaria every year, with a further 1.4 billion people at risk of infection. Little is known about the genetic makeup of P. falciparum populations, despite variation in genetic diversity being a key factor in morbidity, mortality, and the success of malaria control initiatives. Here we analyze a worldwide sample of 519 P. falciparum isolates sequenced for two housekeeping genes (63 single nucleotide polymorphisms from around 5000 nucleotides per isolate). We observe a strong negative correlation between within-population genetic diversity and geographic distance from sub-Saharan Africa (R(2) = 0.95) over Africa, Asia, and Oceania. In contrast, regional variation in transmission intensity seems to have had a negligible impact on the distribution of genetic diversity. The striking geographic patterns of isolation by distance observed in P. falciparum mirror the ones previously documented in humans and point to a joint sub-Saharan African origin between the parasite and its host. Age estimates for the expansion of P. falciparum further support that anatomically modern humans were infected prior to their exit out of Africa and carried the parasite along during their colonization of the world.

Malaria parasite sequences from chimpanzee support the co-speciation hypothesis for the origin of virulent human malaria (Plasmodium falciparum)

Phylogenetic analyses of the mitochondrial cytochrome b (cytb), apicoplast caseinolytic protease C (clpC), and 18S rRNA sequences of Plasmodium isolates from chimpanzees along with those of the virulent human malaria parasite P. falciparum showed that the common chimpanzee (Pan troglodytes) malaria parasites, assigned by Rich et al. (2009) to P. reichenowi, constitute a paraphyletic assemblage. The assumption that P. falciparum diverged from P. reichenowi as recently as 5000-50,000 years ago would require a rate of synonymous substitution/site/year in cytb and clpC on the order of 10(-5)-10(-6), several orders of magnitude higher than any known from eukaryotic organelle genomes, and would imply an unrealistically recent timing of the most recent common ancestor of P. falciparum mitochondrial genomes. The available data are thus most consistent with the hypothesis that P. reichenowi (in the strict sense) and P. falciparum co-speciated with their hosts about 5-7 million years ago.

The selection landscape of malaria parasites

Malaria parasites have to survive and transmit within a highly selective and ever-changing host environment. Because immunity to malaria is nonsterilizing and builds up slowly through repeated infections, commonly the parasite invades a host that is immunologically and physiologically different from its previous host. During the course of infection, the parasite must also keep pace with changes in host immune responses and red-blood-cell physiology. Here, we describe the "selection landscape" of the most virulent of the human malaria parasites, Plasmodium falciparum, and the adaptive mechanisms it uses to navigate through that landscape. Taking a cost-benefit view of parasite fitness, we consider the evolutionary outcomes of the most important forces of selection operating on the parasite, namely immunity, host death, drugs, mosquito availability, and coinfection. Given the huge potential for malaria parasite evolution in the context of the recently renewed effort to eradicate malaria, a deeper understanding of P. falciparum adaptation is essential.

Malaria parasite genome scan: insights into antimalarial resistance

Genome scan and genotype-phenotype association study offer excellent opportunities to unearth drug/vaccine targets in human pathogens including malaria parasites. A recently conducted such study in worldwide isolates in the most devastating malaria parasite Plasmodium falciparum has reported important genomic information on genetic basis of antimalarial resistance. Several unknown genes were also found to be under strong influence of natural selection. The findings provide important insights into the malaria parasite genome evolution in general and to use this information to develop more focused malaria control strategies, in particular.

Identification of a mutant PfCRT-mediated chloroquine tolerance phenotype in Plasmodium falciparum

Mutant forms of the Plasmodium falciparum transporter PfCRT constitute the key determinant of parasite resistance to chloroquine (CQ), the former first-line antimalarial, and are ubiquitous to infections that fail CQ treatment. However, treatment can often be successful in individuals harboring mutant pfcrt alleles, raising questions about the role of host immunity or pharmacokinetics vs. the parasite genetic background in contributing to treatment outcomes. To examine whether the parasite genetic background dictates the degree of mutant pfcrt-mediated CQ resistance, we replaced the wild type pfcrt allele in three CQ-sensitive strains with mutant pfcrt of the 7G8 allelic type prevalent in South America, the Oceanic region and India. Recombinant clones exhibited strain-dependent CQ responses that ranged from high-level resistance to an incremental shift that did not meet CQ resistance criteria. Nonetheless, even in the most susceptible clones, 7G8 mutant pfcrt enabled parasites to tolerate CQ pressure and recrudesce in vitro after treatment with high concentrations of CQ. 7G8 mutant pfcrt was found to significantly impact parasite responses to other antimalarials used in artemisinin-based combination therapies, in a strain-dependent manner. We also report clinical isolates from French Guiana that harbor mutant pfcrt, identical or related to the 7G8 haplotype, and manifest a CQ tolerance phenotype. One isolate, H209, harbored a novel PfCRT C350R mutation and demonstrated reduced quinine and artemisinin susceptibility. Our data: 1) suggest that high-level CQR is a complex biological process dependent on the presence of mutant pfcrt; 2) implicate a role for variant pfcrt alleles in modulating parasite susceptibility to other clinically important antimalarials; and 3) uncover the existence of a phenotype of CQ tolerance in some strains harboring mutant pfcrt.

Plasmodium falciparum resistance to the antimalarial drug chloroquine has been found to result primarily from point mutations in PfCRT, which provide a highly sensitive marker of in vivo treatment failure and in vitro resistance. Debate has nonetheless continued about the singular role of mutant PfCRT and the contribution of the parasite genetic background. To address this, we have generated recombinant P. falciparum lines expressing a mutant pfcrt allele, or the reference wild type allele, in three distinct chloroquine-sensitive strains. Their analysis reveals a spectrum of responses ranging from high-level resistance to a previously unrecognized tolerance phenotype. The latter is characterized by virtually unchanged chloroquine IC₅₀ values, significantly elevated IC₉₀ values, and the ability to recrudesce after exposure to drug concentrations that are lethal to chloroquine-sensitive parasites. This tolerance phenotype was also observed in an isolate from French Guiana, confirming its presence in malaria-endemic regions. Mutant PfCRT significantly affected parasite responses to other antimalarials, including ones used in artemisinin-based combination therapies, in a strain-dependent manner. Our data suggest that successful CQ treatment of drug-resistant parasites is dependent on both host immunity and the strain-dependent extent to which mutant pfcrt imparts resistance.

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.

Restriction Site Tiling Analysis: accurate discovery and quantitative genotyping of genome-wide polymorphisms using nucleotide arrays

A method for the simultaneous identification of polymorphic loci and the quantitative genotyping of thousands of loci in individuals is presented.

High-throughput genotype data can be used to identify genes important for local adaptation in wild populations, phenotypes in lab stocks, or disease-related traits in human medicine. Here we advance microarray-based genotyping for population genomics with Restriction Site Tiling Analysis. The approach simultaneously discovers polymorphisms and provides quantitative genotype data at 10,000s of loci. It is highly accurate and free from ascertainment bias. We apply the approach to uncover genomic differentiation in the purple sea urchin.

Three members of the 6-cys protein family of Plasmodium play a role in gamete fertility

The process of fertilization is critically dependent on the mutual recognition of gametes and in Plasmodium, the male gamete surface protein P48/45 is vital to this process. This protein belongs to a family of 10 structurally related proteins, the so called 6-cys family. To identify the role of additional members of this family in Plasmodium fertilisation, we performed genetic and functional analysis on the five members of the 6-cys family that are transcribed during the gametocyte stage of P. berghei. This analysis revealed that in addition to P48/45, two members (P230 and P47) also play an essential role in the process of parasite fertilization. Mating studies between parasites lacking P230, P48/45 or P47 demonstrate that P230, like P48/45, is a male fertility factor, consistent with the previous demonstration of a protein complex containing both P48/45 and P230. In contrast, disruption of P47 results in a strong reduction of female fertility, while males remain unaffected. Further analysis revealed that gametes of mutants lacking expression of p48/45 or p230 or p47 are unable to either recognise or attach to each other. Disruption of the paralog of p230, p230p, also specifically expressed in gametocytes, had no observable effect on fertilization. These results indicate that the P. berghei 6-cys family contains a number of proteins that are either male or female specific ligands that play an important role in gamete recognition and/or attachment. The implications of low levels of fertilisation that exist even in the absence of these proteins, indicating alternative pathways of fertilisation, as well as positive selection acting on these proteins, are discussed in the context of targeting these proteins as transmission blocking vaccine candidates.

Sexual reproduction for malaria parasites is an essential process and is necessary for parasite transmission between hosts. Fertilisation between female and male gametes occurs in the midgut of the mosquito and proteins on the surface of gametes are principle targets in transmission blocking strategies. Despite their importance, relatively little is known about the malaria proteins involved in fertilisation. In this study we show that two gamete proteins, one expressed on the surface of males, the other on the surface of females, have important roles in the mutual recognition and attachment of gametes. Mutant parasites that lack the presence of these surface proteins show a strong reduction in fertility. Comparison of these gamete surface proteins in different malaria parasites showed that these proteins are evolving rapidly either across their length or at discreet regions/domains. We found, that despite the drastic reduction in zygote formation, low levels of fertilisation can still occur in the absence of these surface proteins, indicating that gametes can use alternative proteins to recognize each other. Both genetic variation of gamete surface proteins and the presence of different fertilisation pathways have important implications for transmission blocking vaccines targeting gamete surface proteins.

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.

Origins of the recent emergence of Plasmodium falciparum pyrimethamine resistance alleles in Madagascar

The combination of sulfadoxine-pyrimethamine is recommended for use as intermittent preventive treatment of malaria during pregnancy and is deployed in Africa. The emergence and the spread of resistant parasites are major threats to such an intervention. We have characterized the Plasmodium falciparum dhfr (pfdhfr) haplotypes and flanking microsatellites in 322 P. falciparum isolates collected from the Comoros Islands and Madagascar. One hundred fifty-six (48.4%) carried the wild-type pfdhfr allele, 19 (5.9%) carried the S108N single-mutation allele, 30 (9.3%) carried the I164L single-mutation allele, 114 (35.4%) carried the N51I/C59R/S108N triple-mutation allele, and 3 (1.0%) carried the N51I/C59R/S108N/I164L quadruple-mutation allele. Microsatellite analysis showed the introduction from the Comoros Islands of the ancestral pfdhfr triple mutant allele of Asian origin and its spread in Madagascar. Evidence for the emergence on multiple occasions of the I164L single-mutation pfdhfr allele in Madagascar was also obtained. Thus, the conditions required to generate mutants with quadruple mutations are met in Madagascar, representing a serious threat to current drug policy.

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.

Population genetic analysis of large sequence polymorphisms in Plasmodium falciparum blood-stage antigens

Plasmodium falciparum, the causative agent of human malaria, invades host erythrocytes using several proteins on the surface of the invasive merozoite, which have been proposed as potential vaccine candidates. Members of the multi-gene PfRh family are surface antigens that have been shown to play a central role in directing merozoites to alternative erythrocyte receptors for invasion. Recently, we identified a large structural polymorphism, a 0.58Kb deletion, in the C-terminal region of the PfRh2b gene, present at a high frequency in parasite populations from Senegal. We hypothesize that this region is a target of humoral immunity. Here, by analyzing 371 P. falciparum isolates we show that this major allele is present at varying frequencies in different populations within Senegal, Africa, and throughout the world. For allelic dimorphisms in the asexual stage antigens, Msp-2 and EBA-175, we find minimal geographic differentiation among parasite populations from Senegal and other African localities, suggesting extensive gene flow among these populations and/or immune-mediated frequency-dependent balancing selection. In contrast, we observe a higher level of inter-population divergence (as measured by F(st)) for the PfRh2b deletion, similar to that observed for SNPs from the sexual stage Pfs45/48 loci, which is postulated to be under directional selection. We confirm that the region containing the PfRh2b polymorphism is a target of humoral immune responses by demonstrating antibody reactivity of endemic sera. Our analysis of inter-population divergence suggests that in contrast to the large allelic dimorphisms in EBA-175 and Msp-2, the presence or absence of the large PfRh2b deletion may not elicit frequency-dependent immune selection, but may be under positive immune selection, having important implications for the development of these proteins as vaccine candidates.

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.

Simultaneous transcription of duplicated var2csa gene copies in individual Plasmodium falciparum parasites

BACKGROUND

Single nucleotide polymorphisms are common in duplicated genes, causing functional preservation, alteration or silencing. The Plasmodium falciparum genes var2csa and Pf332 are duplicated in the haploid genome of the HB3 parasite line. Whereas the molecular function of Pf332 remains to be elucidated, VAR2CSA is known to be the main adhesin in placental parasite sequestration. Sequence variations introduced upon duplication of these genes provide discriminative possibilities to analyze allele-specific transcription with a bearing towards understanding gene dosage impact on parasite biology.

RESULTS

We demonstrate an approach combining real-time PCR allelic discrimination and discriminative RNA-FISH to distinguish between highly similar gene copies in P. falciparum parasites. The duplicated var2csa variants are simultaneously transcribed, both on a population level and intriguingly also in individual cells, with nuclear co-localization of the active genes and corresponding transcripts. This indicates transcriptional functionality of duplicated genes, challenges the dogma of mutually exclusive var gene transcription and suggests mechanisms behind antigenic variation, at least in respect to the duplicated and highly similar var2csa genes.

CONCLUSIONS

Allelic discrimination assays have traditionally been applied to study zygosity in diploid genomes. The assays presented here are instead successfully applied to the identification and evaluation of transcriptional activity of duplicated genes in the haploid genome of the P. falciparum parasite. Allelic discrimination and gene or transcript localization by FISH not only provide insights into transcriptional regulation of genes such as the virulence associated var genes, but also suggest that this sensitive and precise approach could be used for further investigation of genome dynamics and gene regulation.

Discovery, mechanisms of action and combination therapy of artemisinin

Despite great international efforts, malaria still inflicts an enormous toll on human lives, especially in Africa. Throughout history, antimalarial medicines have been one of the most powerful tools in malaria control. However, the acquisition and spread of parasite strains that are resistant to multiple antimalarial drugs have become one of the greatest challenges to malaria treatment, and are associated with the increase in morbidity and mortality in many malaria-endemic countries. To deal with this grave situation, artemisinin-based combinatory therapies (ACTs) have been introduced and widely deployed in malarious regions. Artemisinin is a new class of antimalarial compounds discovered by Chinese scientists from the sweet wormwood Artemisia annua. The potential development of resistance to artemisinins by Plasmodium falciparum threatens the usable lifespan of ACTs, and therefore is a subject of close surveillance and extensive research. Studies at the Thai-Cambodian border, a historical epicenter of multidrug resistance, have detected reduced susceptibility to artemisinins as manifested by prolonged parasite-clearance times, raising considerable concerns on resistance development. Despite this significance, there is still controversy on the mode of action of artemisinins. Although a number of potential cellular targets of artemisinins have been proposed, they remain to be verified experimentally. Here, we review the history of artemisinin discovery, discuss the mode of action and potential drug targets, and present strategies to elucidate resistance mechanisms.

Genetic and genomic analyses of host-pathogen interactions in malaria

The Plasmodium parasite successfully infects and replicates in both human and insect vectors. Population studies in humans have long detected the enormous selective pressure placed by the parasite on its human host, revealing the footprints of co-evolution. Available complete genomic sequences for the human and insect hosts, and additional sequences from multiple field isolates of Plasmodiumfalciparum have identified a wide array of protein and gene families that play a crucial role at the interface of host-parasite interaction. Selected examples of such interactions will be reviewed herein.

Assisted assembly: how to improve a de novo genome assembly by using related species

A method is described for improving low sequence coverage genome assemblies

We describe a new assembly algorithm, where a genome assembly with low sequence coverage, either throughout the genome or locally, due to cloning bias, is considerably improved through an assisting process via a related genome. We show that the information provided by aligning the whole-genome shotgun reads of the target against a reference genome can be used to substantially improve the quality of the resulting assembly.

From parasite genomes to one healthy world: Are we having fun yet?

In 1990, the Human Genome Sequencing Project was established. This laid the ground work for an explosion of sequence data that has since followed. As a result of this effort, the first complete genome of an animal, Caenorhabditis elegans was published in 1998. The sequence of Drosophila melanogaster was made available in March, 2000 and in the following year, working drafts of the human genome were generated with the completed sequence (92%) being released in 2003. Recent advancements and next-generation technologies have made sequencing common place and have infiltrated every aspect of biological research, including parasitology. To date, sequencing of 32 apicomplexa and 24 nematode genomes are either in progress or near completion, and over 600k nematode EST and 200k apicomplexa EST submissions fill the databases. However, the winds have shifted and efforts are now refocusing on how best to store, mine and apply these data to problem solving. Herein we tend not to summarize existing X-omics datasets or present new technological advances that promise future benefits. Rather, the information to follow condenses up-to-date-applications of existing technologies to problem solving as it relates to parasite research. Advancements in non-parasite systems are also presented with the proviso that applications to parasite research are in the making.

Gene copy number variation throughout the Plasmodium falciparum genome

Background

Gene copy number variation (CNV) is responsible for several important phenotypes of the malaria parasite Plasmodium falciparum, including drug resistance, loss of infected erythrocyte cytoadherence and alteration of receptor usage for erythrocyte invasion. Despite the known effects of CNV, little is known about its extent throughout the genome.

Results

We performed a whole-genome survey of CNV genes in P. falciparum using comparative genome hybridisation of a diverse set of 16 laboratory culture-adapted isolates to a custom designed high density Affymetrix GeneChip array. Overall, 186 genes showed hybridisation signals consistent with deletion or amplification in one or more isolate. There is a strong association of CNV with gene length, genomic location, and low orthology to genes in other Plasmodium species. Sub-telomeric regions of all chromosomes are strongly associated with CNV genes independent from members of previously described multigene families. However, ~40% of CNV genes were located in more central regions of the chromosomes. Among the previously undescribed CNV genes, several that are of potential phenotypic relevance are identified.

Conclusion

CNV represents a major form of genetic variation within the P. falciparum genome; the distribution of gene features indicates the involvement of highly non-random mutational and selective processes. Additional studies should be directed at examining CNV in natural parasite populations to extend conclusions to clinical settings.

Stepwise acquisition of pyrimethamine resistance in the malaria parasite

The spread of high-level pyrimethamine resistance in Africa threatens to curtail the therapeutic lifetime of antifolate antimalarials. We studied the possible evolutionary pathways in the evolution of pyrimethamine resistance using an approach in which all possible mutational intermediates were created by site-directed mutagenesis and assayed for their level of drug resistance. The coding sequence for dihydrofolate reductase (DHFR) from the malaria parasite Plasmodium falciparum was mutagenized, and tests were carried out in Escherichia coli under conditions in which the endogenous bacterial enzyme was selectively inhibited. We studied 4 key amino acid replacements implicated in pyrimethamine resistance: N51I, C59R, S108N, and I164L. Using empirical estimates of the mutational spectrum in P. falciparum and probabilities of fixation based on the relative levels of resistance, we found that the predicted favored pathways of drug resistance are consistent with those reported in previous kinetic studies, as well as DHFR polymorphisms observed in natural populations. We found that 3 pathways account for nearly 90% of the simulated realizations of the evolution of pyrimethamine resistance. The most frequent pathway (S108N and then C59R, N51I, and I164L) accounts for more than half of the simulated realizations. Our results also suggest an explanation for why I164L is detected in Southeast Asia and South America, but not at significant frequencies in Africa.

A comparative proteomic analysis of the simple amino acid repeat distributions in Plasmodia reveals lineage specific amino acid selection

Background

Microsatellites have been used extensively in the field of comparative genomics. By studying microsatellites in coding regions we have a simple model of how genotypic changes undergo selection as they are directly expressed in the phenotype as altered proteins. The simplest of these tandem repeats in coding regions are the tri-nucleotide repeats which produce a repeat of a single amino acid when translated into proteins. Tri-nucleotide repeats are often disease associated, and are also known to be unstable to both expansion and contraction. This makes them sensitive markers for studying proteome evolution, in closely related species.

Results

The evolutionary history of the family of malarial causing parasites Plasmodia is complex because of the life-cycle of the organism, where it interacts with a number of different hosts and goes through a series of tissue specific stages. This study shows that the divergence between the primate and rodent malarial parasites has resulted in a lineage specific change in the simple amino acid repeat distribution that is correlated to A–T content. The paper also shows that this altered use of amino acids in SAARs is consistent with the repeat distributions being under selective pressure.

Conclusions

The study shows that simple amino acid repeat distributions can be used to group related species and to examine their phylogenetic relationships. This study also shows that an outgroup species with a similar A–T content can be distinguished based only on the amino acid usage in repeats, and suggest that this might be a useful feature for proteome clustering. The lineage specific use of amino acids in repeat regions suggests that comparative studies of SAAR distributions between proteomes gives an insight into the mechanisms of expansion and the selective pressures acting on the organism.

Characterization of the repertoire diversity of the Plasmodium falciparum stevor multigene family in laboratory and field isolates

Background

The evasion of host immune response by the human malaria parasite Plasmodium falciparum has been linked to expression of a range of variable antigens on the infected erythrocyte surface. Several genes are potentially involved in this process with the var, rif and stevor multigene families being the most likely candidates and coding for rapidly evolving proteins. The high sequence diversity of proteins encoded by these gene families may have evolved as an immune evasion strategy that enables the parasite to establish long lasting chronic infections. Previous findings have shown that the hypervariable region (HVR) of STEVOR has significant sequence diversity both within as well as across different P. falciparum lines. However, these studies did not address whether or not there are ancestral stevor that can be found in different parasites.

Methods

DNA and RNA sequences analysis as well as phylogenetic approaches were used to analyse the stevor sequence repertoire and diversity in laboratory lines and Kilifi (Kenya) fresh isolates.

Results

Conserved stevor genes were identified in different P. falciparum isolates from different global locations. Consistent with previous studies, the HVR of the stevor gene family was found to be highly divergent both within and between isolates. Importantly phylogenetic analysis shows some clustering of stevor sequences both within a single parasite clone as well as across different parasite isolates.

Conclusion

This indicates that the ancestral P. falciparum parasite genome already contained multiple stevor genes that have subsequently diversified further within the different P. falciparum populations. It also confirms that STEVOR is under strong selection pressure.

Gene copy number and malaria biology

Alteration in gene copy number provides a simple way to change expression levels and alter phenotype. This was fully appreciated by bacteriologists more than 25 years ago, but the extent and implications of copy number polymorphism (CNP) have only recently become apparent in other organisms. New methods demonstrate the ubiquity of CNPs in eukaryotes and their medical importance in humans. CNP is also widespread in the Plasmodium falciparum genome and has an important and underappreciated role in determining phenotype. In this review, we summarize the distribution of CNP, its evolutionary dynamics within populations, its functional importance and its mode of evolution.

The role of geography in human adaptation

Various observations argue for a role of adaptation in recent human evolution, including results from genome-wide studies and analyses of selection signals at candidate genes. Here, we use genome-wide SNP data from the HapMap and CEPH-Human Genome Diversity Panel samples to study the geographic distributions of putatively selected alleles at a range of geographic scales. We find that the average allele frequency divergence is highly predictive of the most extreme F_ST values across the whole genome. On a broad scale, the geographic distribution of putatively selected alleles almost invariably conforms to population clusters identified using randomly chosen genetic markers. Given this structure, there are surprisingly few fixed or nearly fixed differences between human populations. Among the nearly fixed differences that do exist, nearly all are due to fixation events that occurred outside of Africa, and most appear in East Asia. These patterns suggest that selection is often weak enough that neutral processes—especially population history, migration, and drift—exert powerful influences over the fate and geographic distribution of selected alleles.

Since the beginning of the study of evolution, people have been fascinated by recent human evolution and adaptation. Despite great progress in our understanding of human history, we still know relatively little about the selection pressures and historical factors that have been important over the past 100,000 years. In that time human populations have spread around the world and adapted in a wide variety of ways to the new environments they have encountered. Here, we investigate the genomic signal of these adaptations using a large set of geographically diverse human populations typed at thousands of genetic markers across the genome. We find that patterns at selected loci are predictable from the patterns found at all markers genome-wide. On the basis of this, we argue that selection has been strongly constrained by the historical relationships and gene flow between populations.

Large-scale genotyping and genetic mapping in Plasmodium parasites

The completion of many malaria parasite genomes provides great opportunities for genomewide characterization of gene expression and high-throughput genotyping. Substantial progress in malaria genomics and genotyping has been made recently, particularly the development of various microarray platforms for large-scale characterization of the Plasmodium falciparum genome. Microarray has been used for gene expression analysis, detection of single nucleotide polymorphism (SNP) and copy number variation (CNV), characterization of chromatin modifications, and other applications. Here we discuss some recent advances in genetic mapping and genomic studies of malaria parasites, focusing on the use of high-throughput arrays for the detection of SNP and CNV in the P. falciparum genome. Strategies for genetic mapping of malaria traits are also discussed.

Sequence conservation in Plasmodium falciparum alpha-helical coiled coil domains proposed for vaccine development

Background

The availability of the P. falciparum genome has led to novel ways to identify potential vaccine candidates. A new approach for antigen discovery based on the bioinformatic selection of heptad repeat motifs corresponding to α-helical coiled coil structures yielded promising results. To elucidate the question about the relationship between the coiled coil motifs and their sequence conservation, we have assessed the extent of polymorphism in putative α-helical coiled coil domains in culture strains, in natural populations and in the single nucleotide polymorphism data available at PlasmoDB.

Methodology/Principal Findings

14 α-helical coiled coil domains were selected based on preclinical experimental evaluation. They were tested by PCR amplification and sequencing of different P. falciparum culture strains and field isolates. We found that only 3 out of 14 α-helical coiled coils showed point mutations and/or length polymorphisms. Based on promising immunological results 5 of these peptides were selected for further analysis. Direct sequencing of field samples from Papua New Guinea and Tanzania showed that 3 out of these 5 peptides were completely conserved. An in silico analysis of polymorphism was performed for all 166 putative α-helical coiled coil domains originally identified in the P. falciparum genome. We found that 82% (137/166) of these peptides were conserved, and for one peptide only the detected SNPs decreased substantially the probability score for α-helical coiled coil formation. More SNPs were found in arrays of almost perfect tandem repeats. In summary, the coiled coil structure prediction was rarely modified by SNPs. The analysis revealed a number of peptides with strictly conserved α-helical coiled coil motifs.

Conclusion/Significance

We conclude that the selection of α-helical coiled coil structural motifs is a valuable approach to identify potential vaccine targets showing a high degree of conservation.

Sequence diversity and evolutionary dynamics of the dimorphic antigen merozoite surface protein-6 and other Msp genes of Plasmodium falciparum

Immune evasion by Plasmodium falciparum is favored by extensive allelic diversity of surface antigens. Some of them, most notably the vaccine-candidate merozoite surface protein (MSP)-1, exhibit a poorly understood pattern of allelic dimorphism, in which all observed alleles group into two highly diverged allelic families with few or no inter-family recombinants. Here we describe contrasting levels and patterns of sequence diversity in genes encoding three MSP-1-associated surface antigens of P. falciparum, ranging from an ancient allelic dimorphism in the Msp-6 gene to a near lack of allelic divergence in Msp-9 to a more classical multi-allele polymorphism in Msp-7. Other members of the Msp-7 gene family exhibit very little polymorphism in non-repetitive regions. A comparison of P. falciparum Msp-6 sequences to an orthologous sequence from P. reichenowi provided evidence for distinct evolutionary histories of the 5' and 3' segments of the dimorphic region in PfMsp-6, consistent with one dimorphic lineage having arisen from recombination between now-extinct ancestral alleles. In addition, we uncovered two surprising patterns of evolution in repetitive sequence. First, in Msp-6, large deletions are associated with (nearly) identical sequence motifs at their borders. Second, a comparison of PfMsp-9 with the P. reichenowi ortholog indicated retention of a significant inter-unit diversity within an 18-base pair repeat within the coding region of P. falciparum, but homogenization in P. reichenowi.

An unusual recent expansion of the C-terminal domain of RNA polymerase II in primate malaria parasites features a motif otherwise found only in mammalian polymerases

The tail of the enzyme RNA polymerase II is responsible for integrating the diverse events of gene expression in eukaryotes and is indispensable for life in yeast, fruit flies, and mice. The tail features a C-terminal domain (CTD), which is comprised of tandemly repeated Y(1)-S(2)-P(3)-T(4)-S(5)-P(6)-S(7) amino acid heptads that are highly conserved across evolutionary lineages, with all mammalian polymerases featuring 52 identical heptad repeats. However, the composition and function of protozoan CTDs remain less well understood. We find that malaria parasites (genus Plasmodium) display an unprecedented plasticity within the length and composition of their CTDs. The CTD in malaria parasites which infect human and nonhuman primates has expanded compared to closely related species that infect rodents or birds. In addition, this variability extends to different isolates within a single species, such as isolates of the human malaria parasite, Plasmodium falciparum. Our results indicate that expanded CTD heptads in malaria parasites correlates with parasitism of primates and provide the first demonstration of polymorphism of the RNA polymerase II CTD within a single species. The expanded set of CTD heptads feature lysine in the seventh position (Y(1)-S(2)-P(3)-T(4)-S(5)-P(6)-K(7)), a sequence only seen otherwise in the distal portion of mammalian polymerases. These observations raise new questions for the radiation of malaria parasites into diverse hosts and for the molecular evolution of RNA polymerase II.

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.

Case-control approach to identify Plasmodium falciparum polymorphisms associated with severe malaria

Background

Studies to identify phenotypically-associated polymorphisms in the Plasmodium falciparum 23 Mb genome will require a dense array of marker loci. It was considered promising to undertake initial allelic association studies to prospect for virulence polymorphisms in Thailand, as the low endemicity would allow higher levels of linkage disequilibrium (LD) than would exist in more highly endemic areas.

Methodology/Principal Findings

Assessment of LD was first made with 11 microsatellite loci widely dispersed in the parasite genome, and 16 microsatellite loci covering a ∼140 kb region of chromosome 2 (an arbitrarily representative non-telomeric part of the genome), in a sample of 100 P. falciparum isolates. The dispersed loci showed minimal LD (Index of Association, I^S_A = 0.013, P = 0.10), while those on chromosome 2 showed significant LD values mostly between loci <5 kb apart. A disease association study was then performed comparing parasites in 113 severe malaria cases and 245 mild malaria controls. Genotyping was performed on almost all polymorphisms in the binding domains of three erythrocyte binding antigens (eba175, eba140 and eba181), and repeat sequence polymorphisms ∼2 kb apart in each of three reticulocyte binding homologues (Rh1, Rh2a/b, and Rh4). Differences between cases and controls were seen for (i) codons 388-90 in eba175, and (ii) a repeat sequence centred on Rh1 codon 667.

Conclusions/Significance

Allelic association studies on P. falciparum require dense genotypic markers, even in a population of only moderate endemicity that has more extensive LD than highly endemic populations. Disease-associated polymorphisms in the eba175 and Rh1 genes encode differences in the middle of previously characterised erythrocyte binding domains, marking these for further investigation.

Hundreds of microsatellites for genotyping Plasmodium yoelii parasites

Genetic crosses have been employed to study various traits of rodent malaria parasites and to locate loci that contribute to drug resistance, immune protection, and disease virulence. Compared with human malaria parasites, genetic crossing of rodent malaria parasites is more easily performed; however, genotyping methods using microsatellites (MSs) or large-scale single nucleotide polymorphisms (SNPs) that have been widely used in typing Plasmodium falciparum are not available for rodent malaria species. Here we report a genome-wide search of the Plasmodium yoelii yoelii (P. yoelii) genome for simple sequence repeats (SSRs) and the identification of nearly 600 polymorphic MS markers for typing the genomes of P. yoelii and Plasmodium berghei. The MS markers are randomly distributed across the 14 physical chromosomes assembled from genome sequences of three rodent malaria species, although some variations in the numbers of MS expected according to chromosome size exist. The majority of the MS markers are AT-rich repeats, similar to those found in the P. falciparum genome. The MS markers provide an important resource for genotyping, lay a foundation for developing linkage maps, and will greatly facilitate genetic studies of P. yoelii.

Optimizing comparative genomic hybridization probes for genotyping and SNP detection in Plasmodium falciparum

Microarray-based comparative genomic hybridizations (CGH) interrogate genomic DNA to identify structural differences such as amplifications and deletions that are easily detected as large signal aberrations. Subtle signal deviations caused by single nucleotide polymorphisms (SNPs) can also be detected but is challenged by a high AT content (81%) in P. falciparum. We compared genome-wide CGH signal to sequence polymorphisms between parasite strains 3D7, HB3, and Dd2 using NimbleGen microarrays. From 23,191 SNPs (excluding var/rif/stevor genes), our CGH probe set detected SNPs with >99.9% specificity but low (<10%) sensitivity. Probe length, melting temperature, GC content, SNP location in the probe, mutation type, and hairpin structures affected SNP sensitivity. Previously unrecognized variable number tandem repeats (VNTRs) also were detected by this method. These findings will guide the redesign of a probe set to optimize an openly available CGH microarray platform for high-resolution genotyping suitable for population genomics studies.

Plasmodium genomics: latest milestone

Our knowledge on comparative genomics of the malaria parasites has advanced a step forward with the publication of the genomes of two primate-infecting malaria parasites: Plasmodium vivax and Plasmodium knowlesi. Even though the genomes of these organisms are the fifth and sixth Plasmodium genomes to be sequenced, respectively, both have revealed previously unknown features, which are discussed in this month's Genome Watch.

Use of high-density tiling microarrays to identify mutations globally and elucidate mechanisms of drug resistance in Plasmodium falciparum

Using tiling microarrays, a mechanism by which Plasmodium falciparum parasites acquire resistance to the antimalarial fosmidomycin has been elucidated

Background

The identification of genetic changes that confer drug resistance or other phenotypic changes in pathogens can help optimize treatment strategies, support the development of new therapeutic agents, and provide information about the likely function of genes. Elucidating mechanisms of phenotypic drug resistance can also assist in identifying the mode of action of uncharacterized but potent antimalarial compounds identified in high-throughput chemical screening campaigns against Plasmodium falciparum.

Results

Here we show that tiling microarrays can detect de novo a large proportion of the genetic changes that differentiate one genome from another. We show that we detect most single nucleotide polymorphisms or small insertion deletion events and all known copy number variations that distinguish three laboratory isolates using readily accessible methods. We used the approach to discover mutations that occur during the selection process after transfection. We also elucidated a mechanism by which parasites acquire resistance to the antimalarial fosmidomycin, which targets the parasite isoprenoid synthesis pathway. Our microarray-based approach allowed us to attribute in vitro derived fosmidomycin resistance to a copy number variation event in the pfdxr gene, which enables the parasite to overcome fosmidomycin-mediated inhibition of isoprenoid biosynthesis.

Conclusions

We show that newly emerged single nucleotide polymorphisms can readily be detected and that malaria parasites can rapidly acquire gene amplifications in response to in vitro drug pressure. The ability to define comprehensively genetic variability in P. falciparum with a single overnight hybridization creates new opportunities to study parasite evolution and improve the treatment and control of malaria.