Multilocus genotyping of microbial pathogens has revealed a range of population structures, with some bacteria showing extensive recombination and others showing almost complete clonality. The population structure of the protozoan parasite Plasmodium falciparum has been harder to evaluate, since most studies have used a limited number of antigen-encoding loci that are known to be under strong selection. We describe length variation at 12 microsatellite loci in 465 infections collected from 9 locations worldwide. These data reveal dramatic differences in parasite population structure in different locations. Strong linkage disequilibrium (LD) was observed in six of nine populations. Significant LD occurred in all locations with prevalence <1% and in only two of five of the populations from regions with higher transmission intensities. Where present, LD results largely from the presence of identical multilocus genotypes within populations, suggesting high levels of self-fertilization in populations with low levels of transmission. We also observed dramatic variation in diversity and geographical differentiation in different regions. Mean heterozygosities in South American countries (0.3-0.4) were less than half those observed in African locations (0. 76-0.8), with intermediate heterozygosities in the Southeast Asia/Pacific samples (0.51-0.65). Furthermore, variation was distributed among locations in South America (F:(ST) = 0.364) and within locations in Africa (F:(ST) = 0.007). The intraspecific patterns of diversity and genetic differentiation observed in P. falciparum are strikingly similar to those seen in interspecific comparisons of plants and animals with differing levels of outcrossing, suggesting that similar processes may be involved. The differences observed may also reflect the recent colonization of non-African populations from an African source, and the relative influences of epidemiology and population history are difficult to disentangle. These data reveal a range of population structures within a single pathogen species and suggest intimate links between patterns of epidemiology and genetic structure in this organism.

Multiple, selectively neutral genetic markers are the most appropriate tools for analysis of parasite population structure and epidemiology, but yet existing methods for characterization of malaria field samples utilize a limited number of antigen encoding genes, which appear to be under strong selection. We describe protocols for characterization of 12 microsatellite markers from finger-prick blood samples infected with Plasmodium falciparum. A two-step, heminested strategy was used to amplify all loci, and products were visualized by fluorescent end-labelling of internal primers. This procedure allows amplification from low levels of template, while eliminating the problem of spurious products due to primer carry over from the primary round of PCR. The loci can be conveniently multiplexed, while accurate sizing and quantification of PCR products can be automated using the GENOTYPER software. The primers do not amplify co-infecting malaria species such as P. vivax and P. malariae. To demonstrate the utility of these markers, we characterized 57 infected finger-prick blood samples from the village of Mebat in Papua New Guinea for all 12 loci, and all samples were genotyped a second time to measure reproducibility. Numbers of alleles per locus range from 4 to 10 in this population, while heterozygosities range from 0.21 to 0.87. Reproducibility (measured as concordance between predominant alleles detected in replicate samples) ranged from 92 to 98% for the 12 loci. The composition of PCR products from infections containing multiple malaria clones could also be defined using strict criteria and scored in a highly repeatable manner.

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.