Tim J C Anderson

We present an optimized probe design for copy number variation (CNV) and SNP genotyping in the Plasmodium falciparum genome. We demonstrate that variable length and isothermal probes are superior to static length probes. We show that sample preparation and hybridization conditions mitigate the effects of host DNA contamination in field samples. The microarray and workflow presented can be used to identify CNVs and SNPs with 95% accuracy in a single hybridization, in field samples containing up to 92% human DNA contamination.

This article outlines genome-scale approaches that can be used to identify mutations in malaria (Plasmodium) parasites that underlie drug resistance and contribute to treatment failure. These approaches include genetic mapping by linkage or genome-wide association studies, drug selection and characterization of resistant mutants, and the identification of genome regions under strong recent selection. While these genomic approaches can identify candidate resistance loci, genetic manipulation is needed to demonstrate causality. We therefore also describe the growing arsenal of available transfection approaches for direct incrimination of mutations suspected to play a role in resistance. Our intention is both to review past progress and highlight promising approaches for future investigations.

In western Cambodia, malaria parasites clear slowly from the blood after treatment with artemisinin derivatives, but it is unclear whether this resul

Malaria parasites vary in phenotypic traits of biomedical or biological interest such as growth rate, virulence, sex ratio and drug resistance, and there is considerable interest in identifying the genes that underlie this variation. An important first step is to determine trait heritability ( H ² ). We evaluate two approaches to measuring H ² in natural parasite populations using relatedness inferred from genetic marker data. We collected single-clone Plasmodium falciparum infections from 185 patients from the Thailand–Burma border, monitored parasite clearance following treatment with artemisinin combination therapy (ACT), measured resistance to six antimalarial drugs and genotyped parasites using 335 microsatellites. We found strong relatedness structure. There were 27 groups of two to eight clonally identical (CI) parasites, and 74 per cent of parasites showed significant relatedness to one or more other parasites. Initially, we used matrices of allele sharing and variance components (VC) methods to estimate H ² . Inhibitory concentrations (IC ₅₀ ) for six drugs showed significant H ² (0.24 to 0.79, p = 0.06 to 2.85 × 10 ⁻⁹ ), demonstrating that this study design has adequate power. However, a phenotype of current interest—parasite clearance following ACT—showed no detectable heritability ( H ² = 0–0.09, ns) in this population. The existence of CI parasites allows the use of a simple ANOVA approach for quantifying H ² , analogous to that used in human twin studies. This gave similar results to the VC method and requires considerably less genotyping information. We conclude (i) that H ² can be effectively measured in malaria parasite populations using minimal genotype data, allowing rational design of genome-wide association studies; and (ii) while drug response (IC ₅₀ ) shows significant H ² , parasite clearance following ACT was not heritable in the population studied.

Background: Schistosoma mansoni is a blood fluke that infects approximately 90 million people. The complete life cycle of this parasite can be maintained in the laboratory, making this one of the few experimentally tractable human helminth infections, and a rich literature reveals heritable variation in important biomedical traits such as virulence, host-specificity, transmission and drug resistance. However, there is a current lack of tools needed to study S. mansoni's molecular, quantitative, and population genetics. Our goal was to construct a genetic linkage map for S. mansoni, and thus provide a new resource that will help stimulate research on this neglected pathogen. Results: We genotyped grandparents, parents and 88 progeny to construct a 5.6 cM linkage map containing 243 microsatellites positioned on 203 of the largest scaffolds in the genome sequence. The map allows 70% of the estimated 300 Mb genome to be ordered on chromosomes, and highlights where scaffolds have been incorrectly assembled. The markers fall into eight main linkage groups, consistent with seven pairs of autosomes and one pair of sex chromosomes, and we were able to anchor linkage groups to chromosomes using fluorescent in situ hybridization. The genome measures 1,228.6 cM. Marker segregation reveals higher female recombination, confirms ZW inheritance patterns, and identifies recombination hotspots and regions of segregation distortion. Conclusions: The genetic linkage map presented here is the first for S. mansoni and the first for a species in the phylum Platyhelminthes. The map provides the critical tool necessary for quantitative genetic analysis, aids genome assembly, and furnishes a framework for comparative flatworm genomics and field-based molecular epidemiological studies.

Tim Anderson discusses a new study of molecular variation in alleles at the dihydropteroate synthase locus, which underlies resistance to sulfadoxine, in over 5,000 parasites from 50 locations.

Multiple displacement amplification (MDA) using Phi29 has proved to be an efficient, high-fidelity method for whole genome amplification in many organisms. This project was designed to evaluate this approach for use with the malaria parasite Plasmodium falciparum. In particular, we were concerned that the AT richness and presence of contaminating human DNA could limit efficiency of MDA in this system. We amplified 60 DNA samples using phi29 and scored 14 microsatellites, 9 single-nucleotide polymorphisms (SNPs), and gene copy number at GTP-cyclohydrolase I both before and after MDA. We observed 100% concordance in 829 microsatellite genotypes and in 499 SNP genotypes. Furthermore, copy number estimates for the GTP-cyclohydrolase I gene were correlated (r² = 0.67) in pre- and postamplification samples. These data confirm that MDA permits scoring of a range of different types of polymorphisms in P. falciparum malaria and can be used to extend the life of valuable DNA stocks.

Copy number polymorphism (CNP) is ubiquitous in eukaryotic genomes, but the degree to which this reflects the action of positive selection is poorly understood. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive CNP. We provide compelling evidence that gch1 CNP is an adaptive consequence of selection by antifolate drugs, which target enzymes downstream in this pathway. (1) We compared gch1 CNP in parasites from Thailand (strong historical antifolate selection) with those from neighboring Laos (weak antifolate selection). Two percent of chromosomes had amplified copy number in Laos, while 72% carried multiple (2–11) copies in Thailand, and differentiation exceeded that observed at 73 synonymous SNPs. (2) We found five amplicon types containing one to greater than six genes and spanning 1 to >11 kb, consistent with parallel evolution and strong selection for this gene amplification. gch1 was the only gene occurring in all amplicons suggesting that this locus is the target of selection. (3) We observed reduced microsatellite variation and increased linkage disequilibrium (LD) in a 900-kb region flanking gch1 in parasites from Thailand, consistent with rapid recent spread of chromosomes carrying multiple copies of gch1. (4) We found that parasites bearing dhfr-164L, which causes high-level resistance to antifolate drugs, carry significantly (p = 0.00003) higher copy numbers of gch1 than parasites bearing 164I, indicating functional association between genes located on different chromosomes but linked in the same biochemical pathway. These results demonstrate that CNP at gch1 is adaptive and the associations with dhfr-164L strongly suggest a compensatory function. More generally, these data demonstrate how selection affects multiple enzymes in a single biochemical pathway, and suggest that investigation of structural variation may provide a fast-track to locating genes underlying adaptation. Recent comparative genomic hybridization studies have revealed extensive copy number variation in eukaryotic genomes. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive copy number polymorphism (CNP). We provide compelling evidence that gch1 CNP is adaptive and most likely results from selection by antifolate drugs, which target enzymes downstream in this pathway. Gch1 CNP shows extreme geographical differentiation; hitchhiking reduces diversity and increases LD in flanking sequence, indicating recent rapid spread within Thailand, while amplicon structure reveals multiple origins and parallel evolution. Furthermore, strong association between elevated copy number and a critical mutation dhfr-I164L that underlies high-level antifolate resistance indicates functional linkage and fitness epistasis between genes on different chromosomes. These data reveal hidden complexity in the evolutionary response to antifolate treatment and demonstrate that analysis of structural variation can provide a fast-track to locating genes that underlie adaptation.