Search

Background: Plasmodium vivax causes the majority of malaria episodes outside Africa, but remains a relatively understudied pathogen. The pathology of P. vivax infection depends critically on the parasite's ability to recognize and invade human erythrocytes. This invasion process involves an interaction between P. vivax Duffy Binding Protein (PvDBP) in merozoites and the Duffy antigen receptor for chemokines (DARC) on the erythrocyte surface. Whole-genome sequencing of clinical isolates recently established that some P. vivax genomes contain two copies of the PvDBP gene. The frequency of this duplication is particularly high in Madagascar, where there is also evidence for P. vivax infection in DARC-negative individuals. The functional significance and global prevalence of this duplication, and whether there are other copy number variations at the PvDBP locus, is unknown., Methodology/principal Findings: Using whole-genome sequencing and PCR to study the PvDBP locus in P. vivax clinical isolates, we found that PvDBP duplication is widespread in Cambodia. The boundaries of the Cambodian PvDBP duplication differ from those previously identified in Madagascar, meaning that current molecular assays were unable to detect it. The Cambodian PvDBP duplication did not associate with parasite density or DARC genotype, and ranged in prevalence from 20% to 38% over four annual transmission seasons in Cambodia. This duplication was also present in P. vivax isolates from Brazil and Ethiopia, but not India., Conclusions/significance: PvDBP duplications are much more widespread and complex than previously thought, and at least two distinct duplications are circulating globally. The same duplication boundaries were identified in parasites from three continents, and were found at high prevalence in human populations where DARC-negativity is essentially absent. It is therefore unlikely that PvDBP duplication is associated with infection of DARC-negative individuals, but functional tests will be required to confirm this hypothesis., Competing Interests: The authors have declared that no competing interests exist.

In regions of coendemicity for Plasmodium falciparum and Plasmodium vivax where mefloquine is used to treat P. falciparum infection, drug pressure mediated by increased copy numbers of the multidrug resistance 1 gene (pvmdr1) may select for mefloquine-resistant P. vivax Surveillance is not undertaken routinely owing in part to methodological challenges in detection of gene amplification. Using genomic data on 88 P. vivax samples from western Thailand, we identified pvmdr1 amplification in 17 isolates, all exhibiting tandem copies of a 37.6-kilobase pair region with identical breakpoints. A novel breakpoint-specific polymerase chain reaction assay was designed to detect the amplification. The assay demonstrated high sensitivity, identifying amplifications in 13 additional, polyclonal infections. Application to 132 further samples identified the common breakpoint in all years tested (2003-2015), with a decline in prevalence after 2012 corresponding to local discontinuation of mefloquine regimens. Assessment of the structure of pvmdr1 amplification in other geographic regions will yield information about the population-specificity of the breakpoints and underlying amplification mechanisms., (© The Author 2016. Published by Oxford University Press for the Infectious Diseases Society of America.)

Background: Severe anemia is a major cause of sickness and death in African children, yet the causes of anemia in this population have been inadequately studied., Methods: We conducted a case-control study of 381 preschool children with severe anemia (hemoglobin concentration, <5.0 g per deciliter) and 757 preschool children without severe anemia in urban and rural settings in Malawi. Causal factors previously associated with severe anemia were studied. The data were examined by multivariate analysis and structural equation modeling., Results: Bacteremia (adjusted odds ratio, 5.3; 95% confidence interval [CI], 2.6 to 10.9), malaria (adjusted odds ratio, 2.3; 95% CI, 1.6 to 3.3), hookworm (adjusted odds ratio, 4.8; 95% CI, 2.0 to 11.8), human immunodeficiency virus infection (adjusted odds ratio, 2.0; 95% CI, 1.0 to 3.8), the G6PD -202/-376 genetic disorder (adjusted odds ratio, 2.4; 95% CI, 1.3 to 4.4), vitamin A deficiency (adjusted odds ratio, 2.8; 95% CI, 1.3 to 5.8), and vitamin B 12 deficiency (adjusted odds ratio, 2.2; 95% CI, 1.4 to 3.6) were associated with severe anemia. Folate deficiency, sickle cell disease, and laboratory signs of an abnormal inflammatory response were uncommon. Iron deficiency was not prevalent in case patients (adjusted odds ratio, 0.37; 95% CI, 0.22 to 0.60) and was negatively associated with bacteremia. Malaria was associated with severe anemia in the urban site (with seasonal transmission) but not in the rural site (where malaria was holoendemic). Seventy-six percent of hookworm infections were found in children under 2 years of age., Conclusions: There are multiple causes of severe anemia in Malawian preschool children, but folate and iron deficiencies are not prominent among them. Even in the presence of malaria parasites, additional or alternative causes of severe anemia should be considered.

The widespread distribution and relapsing nature of Plasmodium vivax infection present major challenges for the elimination of malaria. To characterize the genetic diversity of this parasite in individual infections and across the population, we performed deep genome sequencing of >200 clinical samples collected across the Asia-Pacific region and analyzed data on >300,000 SNPs and nine regions of the genome with large copy number variations. Individual infections showed complex patterns of genetic structure, with variation not only in the number of dominant clones but also in their level of relatedness and inbreeding. At the population level, we observed strong signals of recent evolutionary selection both in known drug resistance genes and at new loci, and these varied markedly between geographical locations. These findings demonstrate a dynamic landscape of local evolutionary adaptation in the parasite population and provide a foundation for genomic surveillance to guide effective strategies for control and elimination of P. vivax., Competing Interests: The authors declare no competing financial interests.

Diversity at pathogen genetic loci can be driven by host adaptive immune selection pressure and may reveal proteins important for parasite biology. Population-based genome sequencing of Plasmodium falciparum, the parasite responsible for the most severe form of malaria, has highlighted two related polymorphic genes called dblmsp and dblmsp2, which encode Duffy binding-like (DBL) domain-containing proteins located on the merozoite surface but whose function remains unknown. Using recombinant proteins and transgenic parasites, we show that DBLMSP and DBLMSP2 directly and avidly bind human IgM via their DBL domains. We used whole genome sequence data from over 400 African and Asian P. falciparum isolates to show that dblmsp and dblmsp2 exhibit extreme protein polymorphism in their DBL domain, with multiple variants of two major allelic classes present in every population tested. Despite this variability, the IgM binding function was retained across diverse sequence representatives. Although this interaction did not seem to have an effect on the ability of the parasite to invade red blood cells, binding of DBLMSP and DBLMSP2 to IgM inhibited the overall immunoreactivity of these proteins to IgG from patients who had been exposed to the parasite. This suggests that IgM binding might mask these proteins from the host humoral immune system., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

Similarity between two individuals in the combination of genetic markers along their chromosomes indicates shared ancestry and can be used to identify historical connections between different population groups due to admixture. We use a genome-wide, haplotype-based, analysis to characterise the structure of genetic diversity and gene-flow in a collection of 48 sub-Saharan African groups. We show that coastal populations experienced an influx of Eurasian haplotypes over the last 7000 years, and that Eastern and Southern Niger-Congo speaking groups share ancestry with Central West Africans as a result of recent population expansions. In fact, most sub-Saharan populations share ancestry with groups from outside of their current geographic region as a result of gene-flow within the last 4000 years. Our in-depth analysis provides insight into haplotype sharing across different ethno-linguistic groups and the recent movement of alleles into new environments, both of which are relevant to studies of genetic epidemiology.

Background: In regions where malaria is endemic, individuals are often infected with multiple distinct parasite genotypes, a situation that may impact on evolution of parasite virulence and drug resistance. Most approaches to studying genotypic diversity have involved analysis of a modest number of polymorphic loci, although whole genome sequencing enables a broader characterisation of samples., Methods: PCR-based microsatellite typing of a panel of ten loci was performed on Plasmodium falciparum in 95 clinical isolates from a highly endemic area in the Republic of Guinea, to characterize within-isolate genetic diversity. Separately, single nucleotide polymorphism (SNP) data from genome-wide short-read sequences of the same samples were used to derive within-isolate fixation indices (F ws), an inverse measure of diversity within each isolate compared to overall local genetic diversity. The latter indices were compared with the microsatellite results, and also with indices derived by randomly sampling modest numbers of SNPs., Results: As expected, the number of microsatellite loci with more than one allele in each isolate was highly significantly inversely correlated with the genome-wide F ws fixation index (r = -0.88, P < 0.001). However, the microsatellite analysis revealed that most isolates contained mixed genotypes, even those that had no detectable genome sequence heterogeneity. Random sampling of different numbers of SNPs showed that an F ws index derived from ten or more SNPs with minor allele frequencies of >10 % had high correlation (r > 0.90) with the index derived using all SNPs., Conclusions: Different types of data give highly correlated indices of within-infection diversity, although PCR-based analysis detects low-level minority genotypes not apparent in bulk sequence analysis. When whole-genome data are not obtainable, quantitative assay of ten or more SNPs can yield a reasonably accurate estimate of the within-infection fixation index (F ws).

Genome-wide searches for loci involved in human resistance to malaria are currently being conducted on a large scale in Africa using case-control studies. Here, we explore the utility of an alternative approach-"environmental correlation analysis, ECA," which tests for clines in allele frequencies across a gradient of an environmental selection pressure-to identify genes that have historically protected against death from malaria. We collected genotype data from 12,425 newborns on 57 candidate malaria resistance loci and 9,756 single nucleotide polymorphisms (SNPs) selected at random from across the genome, and examined their allele frequencies for geographic correlations with long-term malaria prevalence data based on 84,042 individuals living under different historical selection pressures from malaria in coastal Kenya. None of the 57 candidate SNPs showed significant (P < 0.05) correlations in allele frequency with local malaria transmission intensity after adjusting for population structure and multiple testing. In contrast, two of the random SNPs that had highly significant correlations (P < 0.01) were in genes previously linked to malaria resistance, namely, CDH13, encoding cadherin 13, and HS3ST3B1, encoding heparan sulfate 3-O-sulfotransferase 3B1. Both proteins play a role in glycoprotein-mediated cell-cell adhesion which has been widely implicated in cerebral malaria, the most life-threatening form of this disease. Other top genes, including CTNND2 which encodes δ-catenin, a molecular partner to cadherin, were significantly enriched in cadherin-mediated pathways affecting inflammation of the brain vascular endothelium. These results demonstrate the utility of ECA in the discovery of novel genes and pathways affecting infectious disease., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Transmission of the malaria parasite Plasmodium falciparum from humans to the mosquito vector requires differentiation of a sub-population of asexual forms replicating within red blood cells into non-dividing male and female gametocytes. The nature of the molecular mechanism underlying this key differentiation event required for malaria transmission is not fully understood., Methods: Whole genome sequencing was used to examine the genomic diversity of the gametocyte non-producing 3D7-derived lines F12 and A4. These lines were used in the recent detection of the PF3D7&#95;1222600 locus (encoding PfAP2-G), which acts as a genetic master switch that triggers gametocyte development., Results: The evolutionary changes from the 3D7 parental strain through its derivatives F12 (culture-passage derived cloned line) and A4 (transgenic cloned line) were identified. The genetic differences including the formation of chimeric var genes are presented., Conclusion: A genomics resource is provided for the further study of gametocytogenesis or other phenotypes using these parasite lines.

Michael Parker and Dominic Kwiatkowski discuss important ethical considerations for sustainable genomics research in Africa.

Background: Glucose-6-phosphate dehydrogenase (G6PD) deficiency exhibits considerable allelic heterogeneity which manifests with variable biochemical and clinical penetrance. It has long been thought that G6PD deficiency confers partial protection against severe malaria, however prior genetic association studies have disagreed with regard to the strength and specificity of a protective effect, which might reflect differences in the host genetic background, environmental influences, or in the specific clinical phenotypes considered., Methods: A case-control association study of severe malaria was conducted in The Gambia, a region in West Africa where there is considerable allelic heterogeneity underlying expression of G6PD deficiency trait, evaluating the three major nonsynonymous polymorphisms known to be associated with enzyme deficiency (A968G, T542A, and C202T) in a cohort of 3836 controls and 2379 severe malaria cases., Results: Each deficiency allele exhibited a similar trend toward protection against severe malaria overall (15-26% reduced risk); however, in stratifying severe malaria to two of its constituent clinical subphenotypes, severe malarial anaemia (SMA) and cerebral malaria (CM), the three deficiency alleles exhibited trends of opposing effect, with risk conferred to SMA and protection with respect to CM. To assess the overall effect of G6PD deficiency trait, deficiency alleles found across all three loci were pooled. G6PD deficiency trait was found to be significantly associated with protection from severe malaria overall (OR 0.83 [0.75-0.92], P = 0.0006), but this was limited to CM (OR 0.73 [0.61-0.87], P = 0.0005), with a trend toward increased risk for SMA, especially in fully-deficient individuals (OR 1.43 [0.99-2.08], P = 0.056). Sex-stratified testing largely comported with these results, with evidence suggesting that protection by G6PD deficiency trait is conferred to both males and females, though susceptibility to SMA may be restricted to fully-deficient male hemizygotes., Conclusions: In a part of Africa where multiple alleles contribute to expression of G6PD deficiency trait, these findings clarify and extend previous work done in populations where a single variant predominates, and taken together suggest a causal role for G6PD deficiency trait itself with respect to severe malaria, with opposing effects seen on two major clinical subphenotypes.

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

Background: The global prevalence of X-linked glucose-6-phosphate dehydrogenase (G6PD) deficiency is thought to be a result of selection by malaria, but epidemiological studies have yielded confusing results. We investigated the relationships between G6PD deficiency and both malaria and non-malarial illnesses among children in Kenya., Methods: We did this study in Kilifi County, Kenya, where the G6PD c.202T allele is the only significant cause of G6PD deficiency. We tested the associations between G6PD deficiency and severe and complicated Plasmodium falciparum malaria through a case-control study of 2220 case and 3940 control children. Cases were children aged younger than 14 years, who visited the high dependency ward of Kilifi County Hospital with severe malaria between March 1, 1998, and Feb 28, 2010. Controls were children aged between 3-12 months who were born within the same study area between August 2006, and September 2010. We assessed the association between G6PD deficiency and both uncomplicated malaria and other common diseases of childhood in a cohort study of 752 children aged younger than 10 years. Participants of this study were recruited from a representative sample of households within the Ngerenya and Chonyi areas of Kilifi County between Aug 1, 1998, and July 31, 2001. The primary outcome measure for the case-control study was the odds ratio for hospital admission with severe malaria (computed by logistic regression) while for the cohort study it was the incidence rate ratio for uncomplicated malaria and non-malaria illnesses (computed by Poisson regression), by G6PD deficiency category., Findings: 2863 (73%) children in the control group versus 1643 (74%) in the case group had the G6PD normal genotype, 639 (16%) versus 306 (14%) were girls heterozygous for G6PD c.202T, and 438 (11%) versus 271 (12%) children were either homozygous girls or hemizygous boys. Compared with boys and girls without G6PD deficiency, we found significant protection from severe malaria (odds ratio [OR] 0·82, 95% CI 0·70-0·97; p=0·020) among G6PD c.202T heterozygous girls but no evidence for protection among G6PD c.202T hemizygous boys and homozygous girls (OR 1·18, 0·99-1·40; p=0·056). Median follow-up for the mild disease cohort study was 2·24 years (IQR 2·22-2·85). G6PD c.202T had no effect on other common diseases of childhood in heterozygous girls (incidence rate ratio 0·98, 95% CI 0·86-1·11; p=0·82) or homozygous girls or hemizygous boys (0·93, 0·82-1·04; p=0·25), with the sole exception of a marginally significant increase in the incidence of helminth infections among heterozygous girls., Interpretation: Heterozygous girls might be the driving force for the positive selection of G6PD deficiency alleles. Further studies are needed to definitively establish the mechanisms by which G6PD deficiency confers an advantage against malaria in heterozygous individuals. Such studies could lead to the development of new treatments., Funding: Wellcome Trust, UK Medical Research Council, European Union, and Foundation for the National Institutes of Health (as part of the Bill & Melinda Gates Grand Challenges in Global Health Initiative)., (Copyright © 2015 Uyoga et al. Open Access article distributed under the terms of CC BY. Published by Elsevier Ltd.. All rights reserved.)

Background: Many studies report associations between human genetic factors and immunity to malaria but few have been reliably replicated. These studies are usually country-specific, use small sample sizes and are not directly comparable due to differences in methodologies. This study brings together samples and data collected from multiple sites across Africa and Asia to use standardized methods to look for consistent genetic effects on anti-malarial antibody levels., Methods: Sera, DNA samples and clinical data were collected from 13,299 individuals from ten sites in Senegal, Mali, Burkina Faso, Sudan, Kenya, Tanzania, and Sri Lanka using standardized methods. DNA was extracted and typed for 202 Single Nucleotide Polymorphisms with known associations to malaria or antibody production, and antibody levels to four clinical grade malarial antigens [AMA1, MSP1, MSP2, and (NANP)4] plus total IgE were measured by ELISA techniques. Regression models were used to investigate the associations of clinical and genetic factors with antibody levels., Results: Malaria infection increased levels of antibodies to malaria antigens and, as expected, stable predictors of anti-malarial antibody levels included age, seasonality, location, and ethnicity. Correlations between antibodies to blood-stage antigens AMA1, MSP1 and MSP2 were higher between themselves than with antibodies to the (NANP)4 epitope of the pre-erythrocytic circumsporozoite protein, while there was little or no correlation with total IgE levels. Individuals with sickle cell trait had significantly lower antibody levels to all blood-stage antigens, and recessive homozygotes for CD36 (rs321198) had significantly lower anti-malarial antibody levels to MSP2., Conclusion: Although the most significant finding with a consistent effect across sites was for sickle cell trait, its effect is likely to be via reducing a microscopically positive parasitaemia rather than directly on antibody levels. However, this study does demonstrate a framework for the feasibility of combining data from sites with heterogeneous malaria transmission levels across Africa and Asia with which to explore genetic effects on anti-malarial immunity.

Background: Genome wide sequence analyses of malaria parasites from widely separated areas of the world have identified contrasting population structures and signatures of selection. To compare relatively closely situated but ecologically contrasting regions within an endemic African country, population samples of Plasmodium falciparum clinical isolates were collected in Ghana from Kintampo in the central forest-savannah area, and Navrongo in a drier savannah area ~350 km to the north with more seasonally-restricted transmission. Parasite DNA was sequenced and paired-end reads mapped to the P. falciparum reference genome., Results: High coverage genome wide sequence data for 85 different clinical isolates enabled analysis of 121,712 single nucleotide polymorphisms (SNPs). The local populations had similar proportions of mixed genotype infections, similar SNP allele frequency distributions, and eleven chromosomal regions had elevated integrated haplotype scores (|iHS|) in both. A between-population Rsb metric comparing extended haplotype homozygosity indicated a stronger signal within Kintampo for one of these regions (on chromosome 14) and in Navrongo for two of these regions (on chromosomes 10 and 13). At least one gene in each of these identified regions is a potential target of locally varying selection. The candidates include genes involved in parasite development in mosquitoes, members of variant-expressed multigene families, and a leading vaccine-candidate target of immunity., Conclusions: Against a background of very similar population structure and selection signatures in the P. falciparum populations of Ghana, three narrow genomic regions showed evidence indicating local differences in historical timing or intensity of selection. Sampling of closely situated populations across heterogeneous environments has potential to refine the mapping of important loci under temporally or spatially varying selection.

We report a large multicenter genome-wide association study of Plasmodium falciparum resistance to artemisinin, the frontline antimalarial drug. Across 15 locations in Southeast Asia, we identified at least 20 mutations in kelch13 (PF3D7&#95;1343700) affecting the encoded propeller and BTB/POZ domains, which were associated with a slow parasite clearance rate after treatment with artemisinin derivatives. Nonsynonymous polymorphisms in fd (ferredoxin), arps10 (apicoplast ribosomal protein S10), mdr2 (multidrug resistance protein 2) and crt (chloroquine resistance transporter) also showed strong associations with artemisinin resistance. Analysis of the fine structure of the parasite population showed that the fd, arps10, mdr2 and crt polymorphisms are markers of a genetic background on which kelch13 mutations are particularly likely to arise and that they correlate with the contemporary geographical boundaries and population frequencies of artemisinin resistance. These findings indicate that the risk of new resistance-causing mutations emerging is determined by specific predisposing genetic factors in the underlying parasite population.

Background: The emergence of artemisinin-resistant Plasmodium falciparum in Southeast Asia threatens malaria treatment efficacy. Mutations in a kelch protein encoded on P. falciparum chromosome 13 (K13) have been associated with resistance in vitro and in field samples from Cambodia., Methods: P. falciparum infections from artesunate efficacy trials in Bangladesh, Cambodia, Laos, Myanmar, and Vietnam were genotyped at 33 716 genome-wide single-nucleotide polymorphisms (SNPs). Linear mixed models were used to test associations between parasite genotypes and parasite clearance half-lives following artesunate treatment. K13 mutations were tested for association with artemisinin resistance, and extended haplotypes on chromosome 13 were examined to determine whether mutations arose focally and spread or whether they emerged independently., Results: The presence of nonreference K13 alleles was associated with prolonged parasite clearance half-life (P = 1.97 × 10(-12)). Parasites with a mutation in any of the K13 kelch domains displayed longer parasite clearance half-lives than parasites with wild-type alleles. Haplotype analysis revealed both population-specific emergence of mutations and independent emergence of the same mutation in different geographic areas., Conclusions: K13 appears to be a major determinant of artemisinin resistance throughout Southeast Asia. While we found some evidence of spreading resistance, there was no evidence of resistance moving westward from Cambodia into Myanmar., (© The Author 2014. Published by Oxford University Press on behalf of the Infectious Diseases Society of America. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Background: Glucose-6-phosphate dehydrogenase (G6PD) is an enzyme that plays an important role in many cellular functions. Deficiency of this enzyme results from point mutations in the coding region of the G6PD gene. G6PD-deficiency is important in malaria, as certain anti-malarial drugs could induce haemolysis in such patients and mutations in this gene may influence the susceptibility or resistance to the disease. Detailed information on genetic variations in the G6PD gene for Sri Lankan populations is yet to be revealed. This study describes a set of G6PD mutations present in a Sri Lankan population and their association with uncomplicated malaria., Methods: DNA was extracted from 1,051 individuals. Sixty-eight SNPs in the region of the G6PD gene were genotyped. A database created during the 1992-1993 malaria epidemic for the same individuals was used to assess the associations between the G6PD SNPs and parasite density or disease severity of uncomplicated malaria infections. Linkage disequilibrium for SNPs and haplotype structures were identified., Results: Seventeen genetic variants were polymorphic in this population. The mutant allele was the major allele in 9 SNPs. Common G6PD variants already described in Asians or South-Asians seemed to be absent or rare in this population. Both the severity of disease in uncomplicated malaria infections and parasitaemia were significantly lower in males infected with Plasmodium falciparum carrying the ancestral allele of rs915942 compared to those carrying the mutant allele. The parasite density of males infected with P. falciparum was significantly lower also in those who possessed the mutant alleles of rs5986877, rs7879049 and rs7053878. Two haplotype blocks were identified, where the recombination rates were higher in males with no history of malaria when compared to those who have experienced the disease in the past., Conclusions: This is the most detailed survey of G6PD SNPs in a Sri Lankan population undertaken so far that enabled novel description of single nucleotide polymorphisms within the G6PD gene. A few of these genetic variations identified, demonstrated a tendency to be associated with either disease severity or parasite density in uncomplicated disease in males. Known G6PD gene polymorphisms already described from elsewhere were either absent or rare in the local study population.

Background: Selection by host immunity and antimalarial drugs has driven extensive adaptive evolution in Plasmodium falciparum and continues to produce ever-changing landscapes of genetic variation., Methods: We performed whole-genome sequencing of 69 P. falciparum isolates from Malawi and used population genetics approaches to investigate genetic diversity and population structure and identify loci under selection., Results: High genetic diversity (π = 2.4 × 10(-4)), moderately high multiplicity of infection (2.7), and low linkage disequilibrium (500-bp) were observed in Chikhwawa District, Malawi, an area of high malaria transmission. Allele frequency-based tests provided evidence of recent population growth in Malawi and detected potential targets of host immunity and candidate vaccine antigens. Comparison of the sequence variation between isolates from Malawi and those from 5 geographically dispersed countries (Kenya, Burkina Faso, Mali, Cambodia, and Thailand) detected population genetic differences between Africa and Asia, within Southeast Asia, and within Africa. Haplotype-based tests of selection to sequence data from all 6 populations identified signals of directional selection at known drug-resistance loci, including pfcrt, pfdhps, pfmdr1, and pfgch1., Conclusions: The sequence variations observed at drug-resistance loci reflect differences in each country's historical use of antimalarial drugs and may be useful in formulating local malaria treatment guidelines., (© The Author 2014. Published by Oxford University Press on behalf of the Infectious Diseases Society of America.)

Pathogen genome sequencing directly from clinical samples is quickly gaining importance in genetic and medical research studies. However, low DNA yield from blood-borne pathogens is often a limiting factor. The problem worsens in extremely base-biased genomes such as the AT-rich Plasmodium falciparum. We present a strategy for whole-genome amplification (WGA) of low-yield samples from P. falciparum prior to short-read sequencing. We have developed WGA conditions that incorporate tetramethylammonium chloride for improved amplification and coverage of AT-rich regions of the genome. We show that this method reduces amplification bias and chimera formation. Our data show that this method is suitable for as low as 10 pg input DNA, and offers the possibility of sequencing the parasite genome from small blood samples., (© The Author 2014. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.)

The African continent continues to bear the greatest burden of malaria and the greatest diversity of parasites, mosquito vectors, and human victims. The evolutionary plasticity of malaria parasites and their vectors is a major obstacle to eliminating the disease. Of current concern is the recently reported emergence of resistance to the front-line drug, artemisinin, in South-East Asia in Plasmodium falciparum, which calls for preemptive surveillance of the African parasite population for genetic markers of emerging drug resistance. Here we describe the Plasmodium Diversity Network Africa (PDNA), which has been established across 11 countries in sub-Saharan Africa to ensure that African scientists are enabled to work together and to play a key role in the global effort for tracking and responding to this public health threat., (Copyright © 2014, American Association for the Advancement of Science.)

Background: The relationship between glucose-6-phosphate dehydrogenase (G6PD) deficiency and clinical phenomena such as primaquine-sensitivity and protection from severe malaria remains poorly defined, with past association studies yielding inconsistent and conflicting results. One possibility is that examination of a single genetic variant might underestimate the presence of true effects in the presence of unrecognized functional allelic diversity., Methods: We systematically examined this possibility in Kenya, conducting a fine-mapping association study of erythrocyte G6PD activity in 1828 Kenyan children across 30 polymorphisms at or around the G6PD locus., Results: We demonstrate a strong functional role for c.202G>A (rs1050828), which accounts for the majority of variance in enzyme activity observed (P=1.5×10⁻²⁰⁰, additive model). Additionally, we identify other common variants that exert smaller, intercorrelated effects independent of c.202G>A, and haplotype analyses suggest that each variant tags one of two haplotype motifs that are opposite in sequence identity and effect direction. We posit that these effects are of biological and possible clinical significance, specifically noting that c.376A>G (rs1050829) augments 202AG heterozygote risk for deficiency trait by two-fold (OR = 2.11 [1.12 - 3.84], P=0.014)., Conclusions: Our results suggest that c.202G>A is responsible for the majority of the observed prevalence of G6PD deficiency trait in Kenya, but also identify a novel role for c.376A>G as a genetic modifier which marks a common haplotype that augments the risk conferred to 202AG heterozygotes, suggesting that variation at both loci merits consideration in genetic association studies probing G6PD deficiency-associated clinical phenotypes.

Background: The practice of making datasets publicly available for use by the wider scientific community has become firmly integrated in genomic science. One significant gap in literature around data sharing concerns how it impacts on scientists' ability to preserve values and ethical standards that form an essential component of scientific collaborations. We conducted a qualitative sociological study examining the potential for harm to ethnic groups, and implications of such ethical concerns for data sharing. We focused our empirical work on the MalariaGEN Consortium, one of the first international collaborative genomics research projects in Africa., Methods: We conducted a study in three MalariaGEN project sites in Kenya, the Gambia, and the United Kingdom. The study entailed analysis of project documents and 49 semi-structured interviews with fieldworkers, researchers and ethics committee members., Results: Concerns about how best to address the potential for harm to ethnic groups in MalariaGEN crystallised in discussions about the development of a data sharing policy. Particularly concerning for researchers was how best to manage the sharing of genomic data outside of the original collaboration. Within MalariaGEN, genomic data is accompanied by information about the locations of sample collection, the limitations of consent and ethics approval, and the values and relations that accompanied sample collection. For interviewees, this information and context were of important ethical value in safeguarding against harmful uses of data, but is not customarily shared with secondary data users. This challenged the ability of primary researchers to protect against harmful uses of 'their' data., Conclusion: We identified three protective mechanisms--trust, the existence of a shared morality, and detailed contextual understanding--which together might play an important role in preventing the use of genomic data in ways that could harm the ethnic groups included in the study. We suggest that the current practice of sharing of datasets as isolated objects rather than as embedded within a particular scientific culture, without regard for the normative context within which samples were collected, may cause ethical tensions to emerge that could have been prevented or addressed had the 'ethical metadata' that accompanies genomic data also been shared.

Background: Artemisinin resistance in Plasmodium falciparum has emerged in Southeast Asia and now poses a threat to the control and elimination of malaria. Mapping the geographic extent of resistance is essential for planning containment and elimination strategies., Methods: Between May 2011 and April 2013, we enrolled 1241 adults and children with acute, uncomplicated falciparum malaria in an open-label trial at 15 sites in 10 countries (7 in Asia and 3 in Africa). Patients received artesunate, administered orally at a daily dose of either 2 mg per kilogram of body weight per day or 4 mg per kilogram, for 3 days, followed by a standard 3-day course of artemisinin-based combination therapy. Parasite counts in peripheral-blood samples were measured every 6 hours, and the parasite clearance half-lives were determined., Results: The median parasite clearance half-lives ranged from 1.9 hours in the Democratic Republic of Congo to 7.0 hours at the Thailand-Cambodia border. Slowly clearing infections (parasite clearance half-life >5 hours), strongly associated with single point mutations in the "propeller" region of the P. falciparum kelch protein gene on chromosome 13 (kelch13), were detected throughout mainland Southeast Asia from southern Vietnam to central Myanmar. The incidence of pretreatment and post-treatment gametocytemia was higher among patients with slow parasite clearance, suggesting greater potential for transmission. In western Cambodia, where artemisinin-based combination therapies are failing, the 6-day course of antimalarial therapy was associated with a cure rate of 97.7% (95% confidence interval, 90.9 to 99.4) at 42 days., Conclusions: Artemisinin resistance to P. falciparum, which is now prevalent across mainland Southeast Asia, is associated with mutations in kelch13. Prolonged courses of artemisinin-based combination therapies are currently efficacious in areas where standard 3-day treatments are failing. (Funded by the U.K. Department of International Development and others; ClinicalTrials.gov number, NCT01350856.).

Background: Plasmodium falciparum malaria is one of the most widespread and deadliest infectious diseases in children under five years in endemic areas. The disease has been a strong force for evolutionary selection in the human genome, and uncovering the critical host genetic factors that confer resistance to the disease would provide clues to the molecular basis of protective immunity and improve vaccine development initiatives., Methods: The effect of single nucleotide polymorphisms (SNPs) and plasma transforming growth factor beta (TGF-β) and interleukin 10 (IL-10) levels on malaria pathology was investigated in a case-control study of 1862 individuals from two major ethnic groups in three regions with intense perennial P. falciparum transmission in Cameroon. Thirty-four malaria candidate polymorphisms, including the sickle cell trait (HbS), were assayed on the Sequenom iPLEX platform while plasma TGF-β and IL-10 levels were measured by sandwich ELISA., Results: The study confirms the known protective effect of HbS against severe malaria and also reveals a protective effect of SNPs in the nitrogen oxide synthase 2 (NOS2) gene against malaria infection, anaemia and uncomplicated malaria. Furthermore, ADCY9 rs10775349 (additive G) and ABO rs8176746 AC individuals were associated with protection from hyperpyrexia and hyperparasitaemia, respectively. Meanwhile, individuals with the EMR1 rs373533 GT, EMR1 rs461645 CT and RTN3 rs542998 (additive C) genotypes were more susceptible to hyperpyrexia while both females and males with the rs1050828 and rs1050829 SNPs of G6PD, respectively, were more vulnerable to anaemia. Plasma TGF-β levels were strongly correlated with heterozygosity for the ADCY9 rs2230739 and HBB rs334 SNPs while individuals with the ABO rs8176746 AC genotype had lower IL-10 levels., Conclusion: Taken together, this study suggests that some rare polymorphisms in candidate genes may have important implications for the susceptibility of Cameroonians to severe malaria. Moreover using the uncomplicated malaria phenotype may permit the identification of novel pathways in the early development of the disease.

The ancient biological 'arms race' between microbial pathogens and humans has shaped genetic variation in modern populations, and this has important implications for the growing field of medical genomics. As humans migrated throughout the world, populations encountered distinct pathogens, and natural selection increased the prevalence of alleles that are advantageous in the new ecosystems in both host and pathogens. This ancient history now influences human infectious disease susceptibility and microbiome homeostasis, and contributes to common diseases that show geographical disparities, such as autoimmune and metabolic disorders. Using new high-throughput technologies, analytical methods and expanding public data resources, the investigation of natural selection is leading to new insights into the function and dysfunction of human biology.

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

Background:  Analysis of genome-wide polymorphism in many organisms has potential to identify genes under recent selection. However, data on historical allele frequency changes are rarely available for direct confirmation., Methods:  We genotyped single nucleotide polymorphisms (SNPs) in 4 Plasmodium falciparum drug resistance genes in 668 archived parasite-positive blood samples of a Gambian population between 1984 and 2008. This covered a period before antimalarial resistance was detected locally, through subsequent failure of multiple drugs until introduction of artemisinin combination therapy. We separately performed genome-wide sequence analysis of 52 clinical isolates from 2008 to prospect for loci under recent directional selection., Results:  Resistance alleles increased from very low frequencies, peaking in 2000 for chloroquine resistance-associated crt and mdr1 genes and at the end of the survey period for dhfr and dhps genes respectively associated with pyrimethamine and sulfadoxine resistance. Temporal changes fit a model incorporating likely selection coefficients over the period. Three of the drug resistance loci were in the top 4 regions under strong selection implicated by the genome-wide analysis., Conclusions:  Genome-wide polymorphism analysis of an endemic population sample robustly identifies loci with detailed documentation of recent selection, demonstrating power to prospectively detect emerging drug resistance genes.

Haptoglobin (Hp) scavenges free hemoglobin following malaria-induced hemolysis. Few studies have investigated the relationship between the common Hp variants and the risk of severe malaria, and their results are inconclusive. We conducted a case-control study of 996 children with severe Plasmodium falciparum malaria and 1220 community controls and genotyped for Hp, hemoglobin (Hb) S heterozygotes, and α(+)thalassemia. Hb S heterozygotes and α(+)thalassemia homozygotes were protected from severe malaria (odds ratio [OR], 0.12; 95% confidence interval [CI], 0.07-0.18 and OR, 0.69; 95% CI, 0.53-0.91, respectively). The risk of severe malaria also varied by Hp genotype: Hp2-1 was associated with the greatest protection against severe malaria and Hp2-2 with the greatest risk. Meta-analysis of the current and published studies suggests that Hp2-2 is associated with increased risk of severe malaria compared with Hp2-1. We found a significant interaction between Hp genotype and α(+)thalassemia in predicting risk of severe malaria: Hp2-1 in combination with heterozygous or homozygous α(+)thalassemia was associated with protection from severe malaria (OR, 0.73; 95% CI, 0.54-0.99 and OR, 0.48; 95% CI, 0.32-0.73, respectively), but α(+)thalassemia in combination with Hp2-2 was not protective. This epistatic interaction together with varying frequencies of α(+)thalassemia across Africa may explain the inconsistent relationship between Hp genotype and malaria reported in previous studies.

The life cycles of many parasites involve transitions between disparate host species, requiring these parasites to go through multiple developmental stages adapted to each of these specialized niches. Transmission of malaria parasites (Plasmodium spp.) from humans to the mosquito vector requires differentiation from asexual stages replicating within red blood cells into non-dividing male and female gametocytes. Although gametocytes were first described in 1880, our understanding of the molecular mechanisms involved in commitment to gametocyte formation is extremely limited, and disrupting this critical developmental transition remains a long-standing goal. Here we show that expression levels of the DNA-binding protein PfAP2-G correlate strongly with levels of gametocyte formation. Using independent forward and reverse genetics approaches, we demonstrate that PfAP2-G function is essential for parasite sexual differentiation. By combining genome-wide PfAP2-G cognate motif occurrence with global transcriptional changes resulting from PfAP2-G ablation, we identify early gametocyte genes as probable targets of PfAP2-G and show that their regulation by PfAP2-G is critical for their wild-type level expression. In the asexual blood-stage parasites pfap2-g appears to be among a set of epigenetically silenced loci prone to spontaneous activation. Stochastic activation presents a simple mechanism for a low baseline of gametocyte production. Overall, these findings identify PfAP2-G as a master regulator of sexual-stage development in malaria parasites and mark the first discovery of a transcriptional switch controlling a differentiation decision in protozoan parasites.

P. falciparum malaria is one of the most widespread and deadliest infectious diseases in children under five years in endemic areas. The disease has been a strong force for evolutionary selection in the human genome, and uncovering the critical human genetic factors that confer resistance to the disease would provide clues to the molecular basis of protective immunity that would be invaluable for vaccine development. We investigated the effect of single nucleotide polymorphisms (SNPs) on malaria pathology in a case- control study of 1862 individuals from two major ethnic groups in three regions with intense perennial P. falciparum transmission in Cameroon. Twenty nine polymorphisms in cytokine and toll-like receptor (TLR) genes as well as the sickle cell trait (HbS) were assayed on the Sequenom iPLEX platform. Our results confirm the known protective effect of HbS against severe malaria and also reveal a protective effect of SNPs in interleukin-10 (IL10) cerebral malaria and hyperpyrexia. Furthermore, IL17RE rs708567 GA and hHbS rs334 AT individuals were associated with protection from uncomplicated malaria and anaemia respectively in this study. Meanwhile, individuals with the hHbS rs334 TT, IL10 rs3024500 AA, and IL17RD rs6780995 GA genotypes were more susceptible to severe malarial anaemia, cerebral malaria, and hyperpyrexia respectively. Taken together, our results suggest that polymorphisms in some immune response genes may have important implications for the susceptibility to severe malaria in Cameroonians. Moreover using uncomplicated malaria may allow us to identify novel pathways in the early development of the disease.

Early identification of causal genetic variants underlying antimalarial drug resistance could provide robust epidemiological tools for timely public health interventions. Using a novel natural genetics strategy for mapping novel candidate genes we analyzed >75,000 high quality single nucleotide polymorphisms selected from high-resolution whole-genome sequencing data in 27 isolates of Plasmodium falciparum. We identified genetic variants associated with susceptibility to dihydroartemisinin that implicate one region on chromosome 13, a candidate gene on chromosome 1 (PFA0220w, a UBP1 ortholog) and others (PFB0560w, PFB0630c, PFF0445w) with putative roles in protein homeostasis and stress response. There was a strong signal for positive selection on PFA0220w, but not the other candidate loci. Our results demonstrate the power of full-genome sequencing-based association studies for uncovering candidate genes that determine parasite sensitivity to artemisinins. Our study provides a unique reference for the interpretation of results from resistant infections.

Natural selection is a significant force that shapes the architecture of the human genome and introduces diversity across global populations. The question of whether advantageous mutations have arisen in the human genome as a result of single or multiple mutation events remains unanswered except for the fact that there exist a handful of genes such as those that confer lactase persistence, affect skin pigmentation, or cause sickle cell anemia. We have developed a long-range-haplotype method for identifying genomic signatures of positive selection to complement existing methods, such as the integrated haplotype score (iHS) or cross-population extended haplotype homozygosity (XP-EHH), for locating signals across the entire allele frequency spectrum. Our method also locates the founder haplotypes that carry the advantageous variants and infers their corresponding population frequencies. This presents an opportunity to systematically interrogate the whole human genome whether a selection signal shared across different populations is the consequence of a single mutation process followed subsequently by gene flow between populations or of convergent evolution due to the occurrence of multiple independent mutation events either at the same variant or within the same gene. The application of our method to data from 14 populations across the world revealed that positive-selection events tend to cluster in populations of the same ancestry. Comparing the founder haplotypes for events that are present across different populations revealed that convergent evolution is a rare occurrence and that the majority of shared signals stem from the same evolutionary event., (Copyright © 2013 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

We describe an analysis of genome variation in 825 P. falciparum samples from Asia and Africa that identifies an unusual pattern of parasite population structure at the epicenter of artemisinin resistance in western Cambodia. Within this relatively small geographic area, we have discovered several distinct but apparently sympatric parasite subpopulations with extremely high levels of genetic differentiation. Of particular interest are three subpopulations, all associated with clinical resistance to artemisinin, which have skewed allele frequency spectra and high levels of haplotype homozygosity, indicative of founder effects and recent population expansion. We provide a catalog of SNPs that show high levels of differentiation in the artemisinin-resistant subpopulations, including codon variants in transporter proteins and DNA mismatch repair proteins. These data provide a population-level genetic framework for investigating the biological origins of artemisinin resistance and for defining molecular markers to assist in its elimination.

Combining data from genome-wide association studies (GWAS) conducted at different locations, using genotype imputation and fixed-effects meta-analysis, has been a powerful approach for dissecting complex disease genetics in populations of European ancestry. Here we investigate the feasibility of applying the same approach in Africa, where genetic diversity, both within and between populations, is far more extensive. We analyse genome-wide data from approximately 5,000 individuals with severe malaria and 7,000 population controls from three different locations in Africa. Our results show that the standard approach is well powered to detect known malaria susceptibility loci when sample sizes are large, and that modern methods for association analysis can control the potential confounding effects of population structure. We show that pattern of association around the haemoglobin S allele differs substantially across populations due to differences in haplotype structure. Motivated by these observations we consider new approaches to association analysis that might prove valuable for multicentre GWAS in Africa: we relax the assumptions of SNP-based fixed effect analysis; we apply Bayesian approaches to allow for heterogeneity in the effect of an allele on risk across studies; and we introduce a region-based test to allow for heterogeneity in the location of causal alleles., Competing Interests: The authors have declared that no competing interests exist.

Wild-living chimpanzees and gorillas harbor a multitude of Plasmodium species, including six of the subgenus Laverania, one of which served as the progenitor of Plasmodium falciparum. Despite the magnitude of this reservoir, it is unknown whether apes represent a source of human infections. Here, we used Plasmodium species-specific PCR, single-genome amplification, and 454 sequencing to screen humans from remote areas of southern Cameroon for ape Laverania infections. Among 1,402 blood samples, we found 1,000 to be Plasmodium mitochondrial DNA (mtDNA) positive, all of which contained human parasites as determined by sequencing and/or restriction enzyme digestion. To exclude low-abundance infections, we subjected 514 of these samples to 454 sequencing, targeting a region of the mtDNA genome that distinguishes ape from human Laverania species. Using algorithms specifically developed to differentiate rare Plasmodium variants from 454-sequencing error, we identified single and mixed-species infections with P. falciparum, Plasmodium malariae, and/or Plasmodium ovale. However, none of the human samples contained ape Laverania parasites, including the gorilla precursor of P. falciparum. To characterize further the diversity of P. falciparum in Cameroon, we used single-genome amplification to amplify 3.4-kb mtDNA fragments from 229 infected humans. Phylogenetic analysis identified 62 new variants, all of which clustered with extant P. falciparum, providing further evidence that P. falciparum emerged following a single gorilla-to-human transmission. Thus, unlike Plasmodium knowlesi-infected macaques in southeast Asia, African apes harboring Laverania parasites do not seem to serve as a recurrent source of human malaria, a finding of import to ongoing control and eradication measures.

The cost of whole-genome sequencing (WGS) is decreasing rapidly as next-generation sequencing technology continues to advance, and the prospect of making WGS available for public health applications is becoming a reality. So far, a number of studies have demonstrated the use of WGS as an epidemiological tool for typing and controlling outbreaks of microbial pathogens. Success of these applications is hugely dependent on efficient generation of clean genetic material that is free from host DNA contamination for rapid preparation of sequencing libraries. The presence of large amounts of host DNA severely affects the efficiency of characterizing pathogens using WGS and is therefore a serious impediment to clinical and epidemiological sequencing for health care and public health applications. We have developed a simple enzymatic treatment method that takes advantage of the methylation of human DNA to selectively deplete host contamination from clinical samples prior to sequencing. Using malaria clinical samples with over 80% human host DNA contamination, we show that the enzymatic treatment enriches Plasmodium falciparum DNA up to ∼9-fold and generates high-quality, nonbiased sequence reads covering >98% of 86,158 catalogued typeable single-nucleotide polymorphism loci.

The recent emergence of artemisinin-resistant Plasmodium falciparum malaria in western Cambodia could threaten prospects for malaria elimination. Identification of the genetic basis of resistance would provide tools for molecular surveillance, aiding efforts to contain resistance. Clinical trials of artesunate efficacy were conducted in Bangladesh, in northwestern Thailand near the Myanmar border, and at two sites in western Cambodia. Parasites collected from trial participants were genotyped at 8,079 single nucleotide polymorphisms (SNPs) using a P. falciparum-specific SNP array. Parasite genotypes were examined for signatures of recent positive selection and association with parasite clearance phenotypes to identify regions of the genome associated with artemisinin resistance. Four SNPs on chromosomes 10 (one), 13 (two), and 14 (one) were significantly associated with delayed parasite clearance. The two SNPs on chromosome 13 are in a region of the genome that appears to be under strong recent positive selection in Cambodia. The SNPs on chromosomes 10 and 13 lie in or near genes involved in postreplication repair, a DNA damage-tolerance pathway. Replication and validation studies are needed to refine the location of loci responsible for artemisinin resistance and to understand the mechanism behind it; however, two SNPs on chromosomes 10 and 13 may be useful markers of delayed parasite clearance in surveillance for artemisinin resistance in Southeast Asia.

The lack of an effective licensed vaccine remains one of the most significant gaps in the portfolio of tools being developed to eliminate Plasmodium falciparum malaria. Vaccines targeting erythrocyte invasion - an essential step for both parasite development and malaria pathogenesis - have faced the particular challenge of genetic diversity. Immunity-driven balancing selection pressure on parasite invasion proteins often results in the presence of multiple, antigenically distinct, variants within a population, leading to variant-specific immune responses. Such variation makes it difficult to design a vaccine that covers the full range of diversity, and could potentially facilitate the evolution of vaccine-resistant parasite strains. In this study, we investigate the effect of genetic diversity on invasion inhibition by antibodies to a high priority P. falciparum invasion candidate antigen, P. falciparum Reticulocyte Binding Protein Homologue 5 (PfRH5). Previous work has shown that virally delivered PfRH5 can induce antibodies that protect against a wide range of genetic variants. Here, we show that a full-length recombinant PfRH5 protein expressed in mammalian cells is biochemically active, as judged by saturable binding to its receptor, basigin, and is able to induce antibodies that strongly inhibit P. falciparum growth and invasion. Whole genome sequencing of 290 clinical P. falciparum isolates from across the world identifies only five non-synonymous PfRH5 SNPs that are present at frequencies of 10% or more in at least one geographical region. Antibodies raised against the 3D7 variant of PfRH5 were able to inhibit nine different P. falciparum strains, which between them included all of the five most common PfRH5 SNPs in this dataset, with no evidence for strain-specific immunity. We conclude that protein-based PfRH5 vaccines are an urgent priority for human efficacy trials., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Whole genome sequencing (WGS) of Plasmodium vivax is problematic due to the reliance on clinical isolates which are generally low in parasitaemia and sample volume. Furthermore, clinical isolates contain a significant contaminating background of host DNA which confounds efforts to map short read sequence of the target P. vivax DNA. Here, we discuss a methodology to significantly improve the success of P. vivax WGS on natural (non-adapted) patient isolates. Using 37 patient isolates from Indonesia, Thailand, and travellers, we assessed the application of CF11-based white blood cell filtration alone and in combination with short term ex vivo schizont maturation. Although CF11 filtration reduced human DNA contamination in 8 Indonesian isolates tested, additional short-term culture increased the P. vivax DNA yield from a median of 0.15 to 6.2 ng µl(-1) packed red blood cells (pRBCs) (p = 0.001) and reduced the human DNA percentage from a median of 33.9% to 6.22% (p = 0.008). Furthermore, post-CF11 and culture samples from Thailand gave a median P. vivax DNA yield of 2.34 ng µl(-1) pRBCs, and 2.65% human DNA. In 22 P. vivax patient isolates prepared with the 2-step method, we demonstrate high depth (median 654X coverage) and breadth (≥89%) of coverage on the Illumina GAII and HiSeq platforms. In contrast to the A+T-rich P. falciparum genome, negligible bias was observed in coverage depth between coding and non-coding regions of the P. vivax genome. This uniform coverage will greatly facilitate the detection of SNPs and copy number variants across the genome, enabling unbiased exploration of the natural diversity in P. vivax populations.

The genetic basis for susceptibility to malaria has been studied widely in African populations but less is known of the contribution of specific genetic variants in Asian populations. We genotyped 67 single-nucleotide polymorphisms (SNPs) in 1030 severe malaria cases and 2840 controls from Vietnam. After data quality control, genotyping data of 956 cases and 2350 controls were analysed for 65 SNPs (3 gender confirmation, 62 positioned in/near 42 malarial candidate genes). A total of 14 SNPs were monomorphic and 2 (rs8078340 and rs33950507) were not in Hardy-Weinberg equilibrium in controls (P<0.01). In all, 7/46 SNPs in 6 genes (ICAM1, IL1A, IL17RC, IL13, LTA and TNF) were associated with severe malaria, with 3/7 SNPs in the TNF/LTA region. Genotype-phenotype correlations between SNPs and clinical parameters revealed that genotypes of rs708567 (IL17RC) correlate with parasitemia (P=0.028, r(2)=0.0086), with GG homozygotes having the lowest parasite burden. Additionally, rs708567 GG homozygotes had a decreased risk of severe malaria (P=0.007, OR=0.78 (95% CI; 0.65-0.93)) and death (P=0.028, OR=0.58 (95% CI; 0.37-0.93)) than those with AA and AG genotypes. In summary, variants in six genes encoding adhesion and proinflammatory molecules are associated with severe malaria in the Vietnamese. Further replicative studies in independent populations will be necessary to confirm these findings.

Malaria elimination strategies require surveillance of the parasite population for genetic changes that demand a public health response, such as new forms of drug resistance. Here we describe methods for the large-scale analysis of genetic variation in Plasmodium falciparum by deep sequencing of parasite DNA obtained from the blood of patients with malaria, either directly or after short-term culture. Analysis of 86,158 exonic single nucleotide polymorphisms that passed genotyping quality control in 227 samples from Africa, Asia and Oceania provides genome-wide estimates of allele frequency distribution, population structure and linkage disequilibrium. By comparing the genetic diversity of individual infections with that of the local parasite population, we derive a metric of within-host diversity that is related to the level of inbreeding in the population. An open-access web application has been established for the exploration of regional differences in allele frequency and of highly differentiated loci in the P. falciparum genome.

Background: Genome and transcriptome studies of Plasmodium nucleic acids obtained from parasitized whole blood are greatly improved by depletion of human DNA or enrichment of parasite DNA prior to next-generation sequencing and microarray hybridization. The most effective method currently used is a two-step procedure to deplete leukocytes: centrifugation using density gradient media followed by filtration through expensive, commercially available columns. This method is not easily implemented in field studies that collect hundreds of samples and simultaneously process samples for multiple laboratory analyses. Inexpensive syringes, hand-packed with CF11 cellulose powder, were recently shown to improve ex vivo cultivation of Plasmodium vivax obtained from parasitized whole blood. This study was undertaken to determine whether CF11 columns could be adapted to isolate Plasmodium falciparum DNA from parasitized whole blood and achieve current quantity and purity requirements for Illumina sequencing., Methods: The CF11 procedure was compared with the current two-step standard of leukocyte depletion using parasitized red blood cells cultured in vitro and parasitized blood obtained ex vivo from Cambodian patients with malaria. Procedural variations in centrifugation and column size were tested, along with a range of blood volumes and parasite densities., Results: CF11 filtration reliably produces 500 nanograms of DNA with less than 50% human DNA contamination, which is comparable to that obtained by the two-step method and falls within the current quality control requirements for Illumina sequencing. In addition, a centrifuge-free version of the CF11 filtration method to isolate P. falciparum DNA at remote and minimally equipped field sites in malaria-endemic areas was validated., Conclusions: CF11 filtration is a cost-effective, scalable, one-step approach to remove human DNA from P. falciparum-infected whole blood samples., (© 2012 Venkatesan et al; BioMed Central Ltd.)

Background: Massively parallel sequencing technology is revolutionizing approaches to genomic and genetic research. Since its advent, the scale and efficiency of Next-Generation Sequencing (NGS) has rapidly improved. In spite of this success, sequencing genomes or genomic regions with extremely biased base composition is still a great challenge to the currently available NGS platforms. The genomes of some important pathogenic organisms like Plasmodium falciparum (high AT content) and Mycobacterium tuberculosis (high GC content) display extremes of base composition. The standard library preparation procedures that employ PCR amplification have been shown to cause uneven read coverage particularly across AT and GC rich regions, leading to problems in genome assembly and variation analyses. Alternative library-preparation approaches that omit PCR amplification require large quantities of starting material and hence are not suitable for small amounts of DNA/RNA such as those from clinical isolates. We have developed and optimized library-preparation procedures suitable for low quantity starting material and tolerant to extremely high AT content sequences., Results: We have used our optimized conditions in parallel with standard methods to prepare Illumina sequencing libraries from a non-clinical and a clinical isolate (containing ~53% host contamination). By analyzing and comparing the quality of sequence data generated, we show that our optimized conditions that involve a PCR additive (TMAC), produces amplified libraries with improved coverage of extremely AT-rich regions and reduced bias toward GC neutral templates., Conclusion: We have developed a robust and optimized Next-Generation Sequencing library amplification method suitable for extremely AT-rich genomes. The new amplification conditions significantly reduce bias and retain the complexity of either extremes of base composition. This development will greatly benefit sequencing clinical samples that often require amplification due to low mass of DNA starting material.

No vaccine has yet proven effective against the blood-stages of Plasmodium falciparum, which cause the symptoms and severe manifestations of malaria. We recently found that PfRH5, a P. falciparum-specific protein expressed in merozoites, is efficiently targeted by broadly-neutralizing, vaccine-induced antibodies. Here we show that antibodies against PfRH5 efficiently inhibit the in vitro growth of short-term-adapted parasite isolates from Cambodia, and that the EC(50) values of antigen-specific antibodies against PfRH5 are lower than those against PfAMA1. Since antibody responses elicited by multiple antigens are speculated to improve the efficacy of blood-stage vaccines, we conducted detailed assessments of parasite growth inhibition by antibodies against PfRH5 in combination with antibodies against seven other merozoite antigens. We found that antibodies against PfRH5 act synergistically with antibodies against certain other merozoite antigens, most notably with antibodies against other erythrocyte-binding antigens such as PfRH4, to inhibit the growth of a homologous P. falciparum clone. A combination of antibodies against PfRH4 and basigin, the erythrocyte receptor for PfRH5, also potently inhibited parasite growth. This methodology provides the first quantitative evidence that polyclonal vaccine-induced antibodies can act synergistically against P. falciparum antigens and should help to guide the rational development of future multi-antigen vaccines.

Acquired immunity in vertebrates maintains polymorphisms in endemic pathogens, leading to identifiable signatures of balancing selection. To comprehensively survey for genes under such selection in the human malaria parasite Plasmodium falciparum, we generated paired-end short-read sequences of parasites in clinical isolates from an endemic Gambian population, which were mapped to the 3D7 strain reference genome to yield high-quality genome-wide coding sequence data for 65 isolates. A minority of genes did not map reliably, including the hypervariable var, rifin, and stevor families, but 5,056 genes (90.9% of all in the genome) had >70% sequence coverage with minimum read depth of 5 for at least 50 isolates, of which 2,853 genes contained 3 or more single nucleotide polymorphisms (SNPs) for analysis of polymorphic site frequency spectra. Against an overall background of negatively skewed frequencies, as expected from historical population expansion combined with purifying selection, the outlying minority of genes with signatures indicating exceptionally intermediate frequencies were identified. Comparing genes with different stage-specificity, such signatures were most common in those with peak expression at the merozoite stage that invades erythrocytes. Members of clag, PfMC-2TM, surfin, and msp3-like gene families were highly represented, the strongest signature being in the msp3-like gene PF10&#95;0355. Analysis of msp3-like transcripts in 45 clinical and 11 laboratory adapted isolates grown to merozoite-containing schizont stages revealed surprisingly low expression of PF10&#95;0355. In diverse clonal parasite lines the protein product was expressed in a minority of mature schizonts (<1% in most lines and ∼10% in clone HB3), and eight sub-clones of HB3 cultured separately had an intermediate spectrum of positive frequencies (0.9 to 7.5%), indicating phase variable expression of this polymorphic antigen. This and other identified targets of balancing selection are now prioritized for functional study., Competing Interests: The authors have declared that no competing interests exist.

Studies on DNA from pathogenic organisms, within clinical samples, are often complicated by the presence of large amounts of host, e.g., human DNA. Isolation of pathogen DNA from these samples would improve the efficiency of next-generation sequencing (NGS) and pathogen identification. Here we describe a solution-based hybridisation method for isolation of pathogen DNA from a mixed population. This straightforward and inexpensive technique uses probes made from whole-genome DNA and off-the-shelf reagents. In this study, Escherichia coli DNA was successfully enriched from a mixture of E.coli and human DNA. After enrichment, genome coverage following NGS was significantly higher and the evenness of coverage and GC content were unaffected. This technique was also applied to samples containing a mixture of human and Plasmodium falciparum DNA. The P.falciparum genome is particularly difficult to sequence due to its high AT content (80.6%) and repetitive nature. Post enrichment, a bias in the recovered DNA was observed, with a poorer representation of the AT-rich non-coding regions. This uneven coverage was also observed in pre-enrichment samples, but to a lesser degree. Despite the coverage bias in enriched samples, SNP (single-nucleotide polymorphism) calling in coding regions was unaffected and the majority of samples had over 90% of their coding region covered at 5× depth. This technique shows significant promise as an effective method to enrich pathogen DNA from samples with heavy human contamination, particularly when applied to GC-neutral genomes.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is an X-linked enzymopathy that affects hundreds of millions of people worldwide, conferring increased risk of neonatal jaundice and oxidant-induced hemolytic anemia. Screening and diagnosis of G6PD deficiency is currently performed using genetic or biochemical assays, the former being cost ineffective in populations with significant allelic heterogeneity, and the latter being limited in ability to detect female heterozygotes. Cytochemical assays can obviate these shortcomings, but at the expense of added technical complexity and labor. We describe here a simple, novel cytofluorometric method that extends the classic methemoglobin reduction test, assessing G6PD deficiency at the level of an individual erythrocyte. In preliminary testing in Malian children, there was strong concordance between our method and established genetic and biochemical techniques. The assay is robust and economical, and could serve as a screening method as well as a research tool, especially for high-throughput applications such as flow cytometry.

Our understanding of the composition of multi-clonal malarial infections and the epidemiological factors which shape their diversity remain poorly understood. Traditionally within-host diversity has been defined in terms of the multiplicity of infection (MOI) derived by PCR-based genotyping. Massively parallel, single molecule sequencing technologies now enable individual read counts to be derived on genome-wide datasets facilitating the development of new statistical approaches to describe within-host diversity. In this class of measures the F(WS) metric characterizes within-host diversity and its relationship to population level diversity. Utilizing P. falciparum field isolates from patients in West Africa we here explore the relationship between the traditional MOI and F(WS) approaches. F(WS) statistics were derived from read count data at 86,158 SNPs in 64 samples sequenced on the Illumina GA platform. MOI estimates were derived by PCR at the msp-1 and -2 loci. Significant correlations were observed between the two measures, particularly with the msp-1 locus (P = 5.92×10(-5)). The F(WS) metric should be more robust than the PCR-based approach owing to reduced sensitivity to potential locus-specific artifacts. Furthermore the F(WS) metric captures information on a range of parameters which influence out-crossing risk including the number of clones (MOI), their relative proportions and genetic divergence. This approach should provide novel insights into the factors which correlate with, and shape within-host diversity.

The Brazilian Amazon is a hypo-endemic malaria region with nearly 300,000 cases each year. A variety of genetic polymorphisms, particularly in erythrocyte receptors and immune response related genes, have been described to be associated with susceptibility and resistance to malaria. In order to identify polymorphisms that might be associated with malaria clinical outcomes in a Brazilian Amazonian population, sixty-four human single nucleotide polymorphisms in 37 genes were analyzed using a Sequenom massARRAY iPLEX platform. A total of 648 individuals from two malaria endemic areas were studied, including 535 malaria cases (113 individuals with clinical mild malaria, 122 individuals with asymptomatic infection and 300 individuals with history of previous mild malaria) and 113 health controls with no history of malaria. The data revealed significant associations (p<0.003) between one SNP in the IL10 gene (rs1800896) and one SNP in the TLR4 gene (rs4986790) with reduced risk for clinical malaria, one SNP in the IRF1 gene (rs2706384) with increased risk for clinical malaria, one SNP in the LTA gene (rs909253) with protection from clinical malaria and one SNP in the TNF gene (RS1800750) associated with susceptibility to clinical malaria. Also, a new association was found between a SNP in the CTL4 gene (rs2242665), located at the major histocompatibility complex III region, and reduced risk for clinical malaria. This study represents the first association study from an Amazonian population involving a large number of host genetic polymorphisms with susceptibility or resistance to Plasmodium infection and malaria outcomes. Further studies should include a larger number of individuals, refined parameters and a fine-scale map obtained through DNA sequencing to increase the knowledge of the Amazonian population genetic diversity.