Mapping a 50-spin-qubit network through correlated sensing

Spins associated to optically accessible solid-state defects have emerged as a versatile platform for exploring quantum simulation, quantum sensing and quantum communication. Pioneering experiments have shown the sensing, imaging, and control of multiple nuclear spins surrounding a single electron spin defect. However, the accessible size of these spin networks has been constrained by the spectral resolution of current methods. Here, we map a network of 50 coupled spins through high-resolution correlated sensing schemes, using a single nitrogen-vacancy center in diamond. We develop concatenated double-resonance sequences that identify spin-chains through the network. These chains reveal the characteristic spin frequencies and their interconnections with high spectral resolution, and can be fused together to map out the network. Our results provide new opportunities for quantum simulations by increasing the number of available spin qubits. Additionally, our methods might find applications in nano-scale imaging of complex spin systems external to the host crystal.

Drug resistance and vaccine target surveillance of Plasmodium falciparum using nanopore sequencing in Ghana

Malaria results in over 600,000 deaths annually, with the highest burden of deaths in young children living in sub-Saharan Africa. Molecular surveillance can provide important information for malaria control policies, including detection of antimalarial drug resistance. However, genome sequencing capacity in malaria-endemic countries is limited. We designed and implemented an end-to-end workflow to detect Plasmodium falciparum antimalarial resistance markers and diversity in the vaccine target circumsporozoite protein (csp) using nanopore sequencing in Ghana. We analysed 196 clinical samples and showed that our method is rapid, robust, accurate and straightforward to implement. Importantly, our method could be applied to dried blood spot samples, which are readily collected in endemic settings. We report that P. falciparum parasites in Ghana are mostly susceptible to chloroquine, with persistent sulfadoxine-pyrimethamine resistance and no evidence of artemisinin resistance. Multiple single nucleotide polymorphisms were identified in csp, but their significance is uncertain. Our study demonstrates the feasibility of nanopore sequencing for malaria genomic surveillance in endemic countries.

Genomic analysis of Plasmodium vivax describes patterns of connectivity and putative drivers of adaptation in Ethiopia

Ethiopia has the greatest burden of Plasmodium vivax in Africa, but little is known about the epidemiological landscape of parasites across the country. We analysed the genomic diversity of 137 P. vivax isolates collected nine Ethiopian districts from 2012 to 2016. Signatures of selection were detected by cross-country comparisons with isolates from Thailand (n = 104) and Indonesia (n = 111), representing regions with low and high chloroquine resistance respectively. 26% (35/137) of Ethiopian infections were polyclonal, and 48.5% (17/35) of these comprised highly related clones (within-host identity-by-descent > 25%), indicating frequent co-transmission and superinfection. Parasite gene flow between districts could not be explained entirely by geographic distance, with economic and cultural factors hypothesised to have an impact on connectivity. Amplification of the duffy binding protein gene (pvdbp1) was prevalent across all districts (16-75%). Cross-population haplotype homozygosity revealed positive selection in a region proximal to the putative chloroquine resistance transporter gene (pvcrt-o). An S25P variant in amino acid transporter 1 (pvaat1), whose homologue has recently been implicated in P. falciparum chloroquine resistance evolution, was prevalent in Ethiopia (96%) but not Thailand or Indonesia (35-53%). The genomic architecture in Ethiopia highlights circulating variants of potential public health concern in an endemic setting with evidence of stable transmission.

Genomics of Plasmodium vivax in Colombia reveals evidence of local bottle-necking and inter-country connectivity in the Americas

Colombia aims to eliminate malaria by 2030 but remains one of the highest burden countries in the Americas. Plasmodium vivax contributes half of all malaria cases, with its control challenged by relapsing parasitaemia, drug resistance and cross-border spread. Using 64 Colombian P. vivax genomes collected between 2013 and 2017, we explored diversity and selection in two major foci of transmission: Chocó and Córdoba. Open-access data from other countries were used for comparative assessment of drug resistance candidates and to assess cross-border spread. Across Colombia, polyclonal infections were infrequent (12%), and infection connectivity was relatively high (median IBD = 5%), consistent with low endemicity. Chocó exhibited a higher frequency of polyclonal infections (23%) than Córdoba (7%), although the difference was not significant (P = 0.300). Most Colombian infections carried double pvdhfr (95%) and single pvdhps (71%) mutants, but other drug resistance mutations were less prevalent (< 10%). There was no evidence of selection at the pvaat1 gene, whose P. falciparum orthologue has recently been implicated in chloroquine resistance. Global population comparisons identified other putative adaptations. Within the Americas, low-level connectivity was observed between Colombia and Peru, highlighting potential for cross-border spread. Our findings demonstrate the potential of molecular data to inform on infection spread and adaptation.

Non-O ABO blood group genotypes differ in their associations with Plasmodium falciparum rosetting and severe malaria

Blood group O is associated with protection against severe malaria and reduced size and stability of P. falciparum-host red blood cell (RBC) rosettes compared to non-O blood groups. Whether the non-O blood groups encoded by the specific ABO genotypes AO, BO, AA, BB and AB differ in their associations with severe malaria and rosetting is unknown. The A and B antigens are host RBC receptors for rosetting, hence we hypothesized that the higher levels of A and/or B antigen on RBCs from AA, BB and AB genotypes compared to AO/BO genotypes could lead to larger rosettes, increased microvascular obstruction and higher risk of malaria pathology. We used a case-control study of Kenyan children and in vitro adhesion assays to test the hypothesis that "double dose" non-O genotypes (AA, BB, AB) are associated with increased risk of severe malaria and larger rosettes than "single dose" heterozygotes (AO, BO). In the case-control study, compared to OO, the double dose genotypes consistently had higher odds ratios (OR) for severe malaria than single dose genotypes, with AB (OR 1.93) and AO (OR 1.27) showing most marked difference (p = 0.02, Wald test). In vitro experiments with blood group A-preferring P. falciparum parasites showed that significantly larger rosettes were formed with AA and AB host RBCs compared to OO, whereas AO and BO genotypes rosettes were indistinguishable from OO. Overall, the data show that ABO genotype influences P. falciparum rosetting and support the hypothesis that double dose non-O genotypes confer a greater risk of severe malaria than AO/BO heterozygosity.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae and Anopheles coluzzii

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of malaria control, but the genetics of resistance are only partially understood. We performed a large scale multi-country genome-wide association study of resistance to two widely used insecticides: deltamethrin and pirimiphos-methyl, using sequencing data from An. gambiae and An. coluzzii from ten locations in West Africa. Resistance was highly multi-genic, multi-allelic and variable between populations. While the strongest and most consistent association with deltamethrin resistance came from Cyp6aa1, this was based on several independent copy number variants (CNVs) in An. coluzzii, and on a non-CNV haplotype in An. gambiae. For pirimiphos-methyl, signals included Ace1, cytochrome P450s, glutathione S-transferases and the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes showed evidence of cross-resistance to both insecticides. These locally-varying, multi-allelic patterns highlight the challenges involved in genomic monitoring of resistance, and may form the basis for improved surveillance methods.

Malaria outbreak in Laos driven by a selective sweep for Plasmodium falciparum kelch13 R539T mutants: a genetic epidemiology analysis

BACKGROUND

Malaria outbreaks are important public health concerns that can cause resurgence in endemic regions approaching elimination. We investigated a Plasmodium falciparum outbreak in Attapeu Province, Laos, during the 2020-21 malaria season, using genomic epidemiology methods to elucidate parasite population dynamics and identify its causes.

METHODS

In this genetic analysis, 2164 P falciparum dried blood spot samples were collected from southern Laos between Jan 1, 2017, and April 1, 2021, which included 249 collected during the Attapeu outbreak between April 1, 2020, and April 1, 2021, by routine surveillance. Genetic barcodes obtained from these samples were used to investigate epidemiological changes underpinning the outbreak, estimate population diversity, and analyse population structure. Whole-genome sequencing data from additional historical samples were used to reconstruct the ancestry of outbreak strains using identity-by-descent analyses.

FINDINGS

The outbreak parasite populations were characterised by unprecedented loss of genetic diversity, primarily caused by rapid clonal expansion of a multidrug-resistant strain (LAA1) carrying the kelch13 Arg539Thr (R539T) mutation. LAA1 replaced kelch13 Cys580Tyr (C580Y) mutants resistant to dihydroartemisinin-piperaquine (KEL1/PLA1) as the dominant strain. LAA1 inherited 58·8% of its genome from a strain circulating in Cambodia in 2008. A secondary outbreak strain (LAA2) carried the kelch13 C580Y allele, and a genome that is essentially identical to a Cambodian parasite from 2009. A third, low-frequency strain (LAA7) was a recombinant of KEL1/PLA1 with a kelch13 R539T mutant.

INTERPRETATION

These results strongly suggest that the outbreak was driven by a selective sweep, possibly associated with multidrug-resistant phenotypes of the outbreak strains. Established resistant populations can circulate at low frequencies for years before suddenly overwhelming dominant strains when the conditions for selection become favourable-eg, when front-line therapies change. Genetic surveillance can support elimination by characterising key properties of outbreaks such as population diversity, drug resistance marker prevalence, and the origins of outbreak strains.

FUNDING

Bill & Melinda Gates Foundation; The Global Fund to Fight AIDS, Tuberculosis and Malaria; Wellcome Trust.

TRANSLATION

For the Lao translation of the abstract see Supplementary Materials section.

Genomic variation during culture adaptation of genetically complex Plasmodium falciparum clinical isolates

Experimental studies on the biology of malaria parasites have mostly been based on laboratory-adapted lines, but there is limited understanding of how these may differ from parasites in natural infections. Loss-of-function mutants have previously been shown to emerge during culture of some Plasmodium falciparum clinical isolates in analyses focusing on single-genotype infections. The present study included a broader array of isolates, mostly representing multiple-genotype infections, which are more typical in areas where malaria is highly endemic. Genome sequence data from multiple time points over several months of culture adaptation of 28 West African isolates were analysed, including previously available sequences along with new genome sequences from additional isolates and time points. Some genetically complex isolates eventually became fixed over time to single surviving genotypes in culture, whereas others retained diversity, although proportions of genotypes varied over time. Drug resistance allele frequencies did not show overall directional changes, suggesting that resistance-associated costs are not the main causes of fitness differences among parasites in culture. Loss-of-function mutants emerged during culture in several of the multiple-genotype isolates, affecting genes (including AP2-HS, EPAC and SRPK1) for which loss-of-function mutants were previously seen to emerge in single-genotype isolates. Parasite clones were derived by limiting dilution from six of the isolates, and sequencing identified de novo variants not detected in the bulk isolate sequences. Interestingly, several of these were nonsense mutants and frameshifts disrupting the coding sequence of EPAC, the gene with the largest number of independent nonsense mutants previously identified in laboratory-adapted lines. Analysis of genomic identity by descent to explore relatedness among clones revealed co-occurring non-identical sibling parasites, illustrative of the natural genetic structure within endemic populations.

Lineage-informative microhaplotypes for spatio-temporal surveillance of Plasmodium vivax malaria parasites

Challenges in understanding the origin of recurrent Plasmodium vivax infections constrains the surveillance of antimalarial efficacy and transmission of this neglected parasite. Recurrent infections within an individual may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or new inoculations (reinfection). Molecular inference of familial relatedness (identity-by-descent or IBD) based on whole genome sequence data, together with analysis of the intervals between parasitaemic episodes ("time-to-event" analysis), can help resolve the probable origin of recurrences. Whole genome sequencing of predominantly low-density P. vivax infections is challenging, so an accurate and scalable genotyping method to determine the origins of recurrent parasitaemia would be of significant benefit. We have developed a P. vivax genome-wide informatics pipeline to select specific microhaplotype panels that can capture IBD within small, amplifiable segments of the genome. Using a global set of 615 P. vivax genomes, we derived a panel of 100 microhaplotypes, each comprising 3-10 high frequency SNPs within <200 bp sequence windows. This panel exhibits high diversity in regions of the Asia-Pacific, Latin America and the horn of Africa (median HE = 0.70-0.81) and it captured 89% (273/307) of the polyclonal infections detected with genome-wide datasets. Using data simulations, we demonstrate lower error in estimating pairwise IBD using microhaplotypes, relative to traditional biallelic SNP barcodes. Our panel exhibited high accuracy in predicting the country of origin (median Matthew's correlation coefficient >0.9 in 90% countries tested) and it also captured local infection outbreak and bottlenecking events. The informatics pipeline is available open-source and yields microhaplotypes that can be readily transferred to high-throughput amplicon sequencing assays for surveillance in malaria-endemic regions.

Pf7: an open dataset of Plasmodium falciparum genome variation in 20,000 worldwide samples

We describe the MalariaGEN Pf7 data resource, the seventh release of Plasmodium falciparum genome variation data from the MalariaGEN network.  It comprises over 20,000 samples from 82 partner studies in 33 countries, including several malaria endemic regions that were previously underrepresented.  For the first time we include dried blood spot samples that were sequenced after selective whole genome amplification, necessitating new methods to genotype copy number variations.  We identify a large number of newly emerging crt mutations in parts of Southeast Asia, and show examples of heterogeneities in patterns of drug resistance within Africa and within the Indian subcontinent.  We describe the profile of variations in the C-terminal of the csp gene and relate this to the sequence used in the RTS,S and R21 malaria vaccines.  Pf7 provides high-quality data on genotype calls for 6 million SNPs and short indels, analysis of large deletions that cause failure of rapid diagnostic tests, and systematic characterisation of six major drug resistance loci, all of which can be freely downloaded from the MalariaGEN website.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae s.l

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of the most widespread tools currently used to control malaria. The genetic underpinnings of resistance are still only partially understood, with much of the variance in resistance phenotype left unexplained. We performed a multi-country large scale genome-wide association study of resistance to two insecticides widely used in malaria control: deltamethrin and pirimiphos-methyl. Using a bioassay methodology designed to maximise the phenotypic difference between resistant and susceptible samples, we sequenced 969 phenotyped female An. gambiae and An. coluzzii from ten locations across four countries in West Africa (Benin, Côte d'Ivoire, Ghana and Togo), identifying single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) segregating in the populations. The patterns of resistance association were highly multiallelic and variable between populations, with different genomic regions contributing to resistance, as well as different mutations within a given region. While the strongest and most consistent association with deltamethrin resistance came from the region around Cyp6aa1 , this resistance was based on a combination of several independent CNVs in An. coluzzii , and on a non-CNV bearing haplotype in An. gambiae . Further signals involved a range of cytochrome P450, mitochondrial, and immunity genes. Similarly, for pirimiphos-methyl, while the strongest signal came from the region of Ace1 , more widespread signals included cytochrome P450s, glutathione S-transferases, and a subunit of the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes were associated with resistance to both insecticide classes, suggesting possible cross-resistance mechanisms. These locally-varying, multigenic and multiallelic patterns highlight the challenges involved in genomic monitoring and surveillance of resistance, and form the basis for improvement of methods used to detect and predict resistance. Based on simulations of resistance variants, we recommend that yet larger scale studies, exceeding 500 phenotyped samples per population, are required to better identify associated genomic regions.

A molecular barcode and web-based data analysis tool to identify imported Plasmodium vivax malaria

Traditionally, patient travel history has been used to distinguish imported from autochthonous malaria cases, but the dormant liver stages of Plasmodium vivax confound this approach. Molecular tools offer an alternative method to identify, and map imported cases. Using machine learning approaches incorporating hierarchical fixation index and decision tree analyses applied to 799 P. vivax genomes from 21 countries, we identified 33-SNP, 50-SNP and 55-SNP barcodes (GEO33, GEO50 and GEO55), with high capacity to predict the infection's country of origin. The Matthews correlation coefficient (MCC) for an existing, commonly applied 38-SNP barcode (BR38) exceeded 0.80 in 62% countries. The GEO panels outperformed BR38, with median MCCs > 0.80 in 90% countries at GEO33, and 95% at GEO50 and GEO55. An online, open-access, likelihood-based classifier framework was established to support data analysis (vivaxGEN-geo). The SNP selection and classifier methods can be readily amended for other use cases to support malaria control programs.

An open dataset of Plasmodium vivax genome variation in 1,895 worldwide samples

This report describes the MalariaGEN Pv4 dataset, a new release of curated genome variation data on 1,895 samples of Plasmodium vivax collected at 88 worldwide locations between 2001 and 2017. It includes 1,370 new samples contributed by MalariaGEN and VivaxGEN partner studies in addition to previously published samples from these and other sources. We provide genotype calls at over 4.5 million variable positions including over 3 million single nucleotide polymorphisms (SNPs), as well as short indels and tandem duplications. This enlarged dataset highlights major compartments of parasite population structure, with clear differentiation between Africa, Latin America, Oceania, Western Asia and different parts of Southeast Asia. Each sample has been classified for drug resistance to sulfadoxine, pyrimethamine and mefloquine based on known markers at the dhfr, dhps and mdr1 loci. The prevalence of all of these resistance markers was much higher in Southeast Asia and Oceania than elsewhere. This open resource of analysis-ready genome variation data from the MalariaGEN and VivaxGEN networks is driven by our collective goal to advance research into the complex biology of P. vivax and to accelerate genomic surveillance for malaria control and elimination.

BIRC6 modifies risk of invasive bacterial infection in Kenyan children

Invasive bacterial disease is a major cause of morbidity and mortality in African children. Despite being caused by diverse pathogens, children with sepsis are clinically indistinguishable from one another. In spite of this, most genetic susceptibility loci for invasive infection that have been discovered to date are pathogen specific and are not therefore suggestive of a shared genetic architecture of bacterial sepsis. Here, we utilise probabilistic diagnostic models to identify children with a high probability of invasive bacterial disease among critically unwell Kenyan children with Plasmodium falciparum parasitaemia. We construct a joint dataset including 1445 bacteraemia cases and 1143 severe malaria cases, and population controls, among critically unwell Kenyan children that have previously been genotyped for human genetic variation. Using these data, we perform a cross-trait genome-wide association study of invasive bacterial infection, weighting cases according to their probability of bacterial disease. In doing so, we identify and validate a novel risk locus for invasive infection secondary to multiple bacterial pathogens, that has no apparent effect on malaria risk. The locus identified modifies splicing of BIRC6 in stimulated monocytes, implicating regulation of apoptosis and autophagy in the pathogenesis of sepsis in Kenyan children.

Malaria protection due to sickle haemoglobin depends on parasite genotype

Host genetic factors can confer resistance against malaria¹, raising the question of whether this has led to evolutionary adaptation of parasite populations. Here we searched for association between candidate host and parasite genetic variants in 3,346 Gambian and Kenyan children with severe malaria caused by Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and three regions of the parasite genome, which is not explained by population structure or other covariates, and which is replicated in additional samples. The HbS-associated alleles include nonsynonymous variants in the gene for the acyl-CoA synthetase family member^2–4PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The alleles are in strong linkage disequilibrium and have frequencies that covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

A strong association has been found between three regions of the Plasmodium falciparum genome and sickle haemoglobin in children with severe malaria, suggesting parasites have adapted to overcome natural host immunity.

Genetic surveillance in the Greater Mekong subregion and South Asia to support malaria control and elimination

Background

National Malaria Control Programmes (NMCPs) currently make limited use of parasite genetic data. We have developed GenRe-Mekong, a platform for genetic surveillance of malaria in the Greater Mekong Subregion (GMS) that enables NMCPs to implement large-scale surveillance projects by integrating simple sample collection procedures in routine public health procedures.

Methods

Samples from symptomatic patients are processed by SpotMalaria, a high-throughput system that produces a comprehensive set of genotypes comprising several drug resistance markers, species markers and a genomic barcode. GenRe-Mekong delivers Genetic Report Cards, a compendium of genotypes and phenotype predictions used to map prevalence of resistance to multiple drugs.

Results

GenRe-Mekong has worked with NMCPs and research projects in eight countries, processing 9623 samples from clinical cases. Monitoring resistance markers has been valuable for tracking the rapid spread of parasites resistant to the dihydroartemisinin-piperaquine combination therapy. In Vietnam and Laos, GenRe-Mekong data have provided novel knowledge about the spread of these resistant strains into previously unaffected provinces, informing decision-making by NMCPs.

Conclusions

GenRe-Mekong provides detailed knowledge about drug resistance at a local level, and facilitates data sharing at a regional level, enabling cross-border resistance monitoring and providing the public health community with valuable insights. The project provides a rich open data resource to benefit the entire malaria community.

Funding

The GenRe-Mekong project is funded by the Bill and Melinda Gates Foundation (OPP11188166, OPP1204268). Genotyping and sequencing were funded by the Wellcome Trust (098051, 206194, 203141, 090770, 204911, 106698/B/14/Z) and Medical Research Council (G0600718). A proportion of samples were collected with the support of the UK Department for International Development (201900, M006212), and Intramural Research Program of the National Institute of Allergy and Infectious Diseases.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Malaria is a cause of iron deficiency in African children

Malaria and iron deficiency (ID) are common and interrelated public health problems in African children. Observational data suggest that interrupting malaria transmission reduces the prevalence of ID1. To test the hypothesis that malaria might cause ID, we used sickle cell trait (HbAS, rs334 ), a genetic variant that confers specific protection against malaria2, as an instrumental variable in Mendelian randomization analyses. HbAS was associated with a 30% reduction in ID among children living in malaria-endemic countries in Africa (n = 7,453), but not among individuals living in malaria-free areas (n = 3,818). Genetically predicted malaria risk was associated with an odds ratio of 2.65 for ID per unit increase in the log incidence rate of malaria. This suggests that an intervention that halves the risk of malaria episodes would reduce the prevalence of ID in African children by 49%.

High-throughput genotyping assays for identification of glycophorin B deletion variants in population studies

Glycophorins are the most abundant sialoglycoproteins on the surface of human erythrocyte membranes. Genetic variation in glycophorin region of human chromosome 4 (containing GYPA, GYPB, and GYPE genes) is of interest because the gene products serve as receptors for pathogens of major public health interest, including Plasmodiumsp., Babesiasp., Influenza virus, Vibrio cholerae El Tor Hemolysin, and Escherichia coli. A large structural rearrangement and hybrid glycophorin variant, known as Dantu, which was identified in East African populations, has been linked with a 40% reduction in risk for severe malaria. Apart from Dantu, other large structural variants exist, with the most common being deletion of the whole GYPB gene and its surrounding region, resulting in multiple different deletion forms. In West Africa particularly, these deletions are estimated to account for between 5 and 15% of the variation in different populations, mostly attributed to the forms known as DEL1 and DEL2. Due to the lack of specific variant assays, little is known of the distribution of these variants. Here, we report a modification of a previous GYPB DEL1 assay and the development of a novel GYPB DEL2 assay as high-throughput PCR-RFLP assays, as well as the identification of the crossover/breakpoint for GYPB DEL2. Using 393 samples from three study sites in Ghana as well as samples from HapMap and 1000 G projects for validation, we show that our assays are sensitive and reliable for genotyping GYPB DEL1 and DEL2. To the best of our knowledge, this is the first report of such high-throughput genotyping assays by PCR-RFLP for identifying specific GYPB deletion types in populations. These assays will enable better identification of GYPB deletions for large genetic association studies and functional experiments to understand the role of this gene cluster region in susceptibility to malaria and other diseases.

Emergence of artemisinin-resistant Plasmodium falciparum with kelch13 C580Y mutations on the island of New Guinea

The rapid and aggressive spread of artemisinin-resistant Plasmodium falciparum carrying the C580Y mutation in the kelch13 gene is a growing threat to malaria elimination in Southeast Asia, but there is no evidence of their spread to other regions. We conducted cross-sectional surveys in 2016 and 2017 at two clinics in Wewak, Papua New Guinea (PNG) where we identified three infections caused by C580Y mutants among 239 genotyped clinical samples. One of these mutants exhibited the highest survival rate (6.8%) among all parasites surveyed in ring-stage survival assays (RSA) for artemisinin. Analyses of kelch13 flanking regions, and comparisons of deep sequencing data from 389 clinical samples from PNG, Indonesian Papua and Western Cambodia, suggested an independent origin of the Wewak C580Y mutation, showing that the mutants possess several distinctive genetic features. Identity by descent (IBD) showed that multiple portions of the mutants' genomes share a common origin with parasites found in Indonesian Papua, comprising several mutations within genes previously associated with drug resistance, such as mdr1, ferredoxin, atg18 and pnp. These findings suggest that a P. falciparum lineage circulating on the island of New Guinea has gradually acquired a complex ensemble of variants, including kelch13 C580Y, which have affected the parasites' drug sensitivity. This worrying development reinforces the need for increased surveillance of the evolving parasite populations on the island, to contain the spread of resistance.

Multi-locus genotyping reveals established endemicity of a geographically distinct Plasmodium vivax population in Mauritania, West Africa

BACKGROUND

Plasmodium vivax has been recently discovered as a significant cause of malaria in Mauritania, although very rare elsewhere in West Africa. It has not been known if this is a recently introduced or locally remnant parasite population, nor whether the genetic structure reflects epidemic or endemic transmission.

METHODOLOGY/PRINCIPAL FINDINGS

To investigate the P. vivax population genetic structure in Mauritania and compare with populations previously analysed elsewhere, multi-locus genotyping was undertaken on 100 clinical isolates, using a genome-wide panel of 38 single nucleotide polymorphisms (SNPs), plus seven SNPs in drug resistance genes. The Mauritanian P. vivax population is shown to be genetically diverse and divergent from populations elsewhere, indicated consistently by genetic distance matrix analysis, principal components analyses, and fixation indices. Only one isolate had a genotype clearly indicating recent importation, from a southeast Asian source. There was no linkage disequilibrium in the local parasite population, and only a small number of infections appeared to be closely genetically related, indicating that there is ongoing genetic recombination consistent with endemic transmission. The P. vivax diversity in a remote mining town was similar to that in the capital Nouakchott, with no indication of local substructure or of epidemic population structure. Drug resistance alleles were virtually absent in Mauritania, in contrast with P. vivax in other areas of the world.

CONCLUSIONS/SIGNIFICANCE

The molecular epidemiology indicates that there is long-standing endemic transmission that will be very challenging to eliminate. The virtual absence of drug resistance alleles suggests that most infections have been untreated, and that this endemic infection has been more neglected in comparison to P. vivax elsewhere.

Molecular epidemiology of resistance to antimalarial drugs in the Greater Mekong subregion: an observational study

BACKGROUND

The Greater Mekong subregion is a recurrent source of antimalarial drug resistance in Plasmodium falciparum malaria. This study aimed to characterise the extent and spread of resistance across this entire region between 2007 and 2018.

METHODS

P falciparum isolates from Myanmar, Thailand, Laos, and Cambodia were obtained from clinical trials and epidemiological studies done between Jan 1, 2007, and Dec 31, 2018, and were genotyped for molecular markers (pfkelch, pfcrt, pfplasmepsin2, and pfmdr1) of antimalarial drug resistance. Genetic relatedness was assessed using microsatellite and single nucleotide polymorphism typing of flanking sequences around target genes.

FINDINGS

10 632 isolates were genotyped. A single long pfkelch Cys580Tyr haplotype (from -50 kb to +31·5 kb) conferring artemisinin resistance (PfPailin) now dominates across the eastern Greater Mekong subregion. Piperaquine resistance associated with pfplasmepsin2 gene amplification and mutations in pfcrt downstream of the Lys76Thr chloroquine resistance locus has also developed. On the Thailand-Myanmar border a different pfkelch Cys580Tyr lineage rose to high frequencies before it was eliminated. Elsewhere in Myanmar the Cys580Tyr allele remains widespread at low allele frequencies. Meanwhile a single artemisinin-resistant pfkelch Phe446Ile haplotype has spread across Myanmar. Despite intense use of dihydroartemisinin-piperaquine in Kayin state, eastern Myanmar, both in treatment and mass drug administrations, no selection of piperaquine resistance markers was observed. pfmdr1 amplification, a marker of resistance to mefloquine, remains at low prevalence across the entire region.

INTERPRETATION

Artemisinin resistance in P falciparum is now prevalent across the Greater Mekong subregion. In the eastern Greater Mekong subregion a multidrug resistant P falciparum lineage (PfPailin) dominates. In Myanmar a long pfkelch Phe446Ile haplotype has spread widely but, by contrast with the eastern Greater Mekong subregion, there is no indication of artemisinin combination therapy (ACT) partner drug resistance from genotyping known markers, and no evidence of spread of ACT resistant P falciparum from the east to the west. There is still a window of opportunity to prevent global spread of ACT resistance.

FUNDING

Thailand Science Research and Innovation, Initiative 5%, Expertise France, Wellcome Trust.

Interferon-gamma polymorphisms and risk of iron deficiency and anaemia in Gambian children

Background: Anaemia is a major public health concern especially in African children living in malaria-endemic regions. Interferon-gamma (IFN-γ) is elevated during malaria infection and is thought to influence erythropoiesis and iron status. Genetic variants in the IFN-γ gene (IFNG) are associated with increased IFN-γ production. We investigated putative functional single nucleotide polymorphisms (SNPs) and haplotypes of IFNG in relation to nutritional iron status and anaemia in Gambian children over a malaria season. Methods: We used previously available data from Gambian family trios to determine informative SNPs and then used the Agena Bioscience MassArray platform to type five SNPs from the IFNG gene in a cohort of 780 Gambian children aged 2-6 years. We also measured haemoglobin and biomarkers of iron status and inflammation at the start and end of a malaria season. Results: We identified five IFNG haplotype-tagging SNPs ( IFNG-1616 [rs2069705], IFNG+874 [rs2430561], IFNG+2200 [rs1861493], IFNG+3234 [rs2069718] and IFNG+5612 [rs2069728]). The IFNG+2200C [rs1861493] allele was associated with reduced haemoglobin concentrations (adjusted β -0.44 [95% CI -0.75, -0.12]; Bonferroni adjusted P = 0.03) and a trend towards iron deficiency compared to wild-type at the end of the malaria season in multivariable models adjusted for potential confounders. A haplotype uniquely identified by IFNG+2200C was similarly associated with reduced haemoglobin levels and trends towards iron deficiency, anaemia and iron deficiency anaemia at the end of the malaria season in models adjusted for age, sex, village, inflammation and malaria parasitaemia. Conclusion: We found limited statistical evidence linking IFNG polymorphisms with a risk of developing iron deficiency and anaemia in Gambian children. More definitive studies are needed to investigate the effects of genetically influenced IFN-γ levels on the risk of iron deficiency and anaemia in children living in malaria-endemic areas.

Development of copy number assays for detection and surveillance of piperaquine resistance associated plasmepsin 2/3 copy number variation in Plasmodium falciparum

BACKGROUND

Long regarded as an epicenter of drug-resistant malaria, Southeast Asia continues to provide new challenges to the control of Plasmodium falciparum malaria. Recently, resistance to the artemisinin combination therapy partner drug piperaquine has been observed in multiple locations across Southeast Asia. Genetic studies have identified single nucleotide polymorphisms as well as copy number variations in the plasmepsin 2 and plasmepsin 3 genes, which encode haemoglobin-degrading proteases that associate with clinical and in vitro piperaquine resistance.

RESULTS

To accurately and quickly determine the presence of copy number variations in the plasmepsin 2/3 genes in field isolates, this study developed a quantitative PCR assay using TaqMan probes. Copy number estimates were validated using a separate SYBR green-based quantitative PCR assay as well as a novel PCR-based breakpoint assay to detect the hybrid gene product. Field samples from 2012 to 2015 across three sites in Cambodia were tested using DNA extracted from dried blood spots and whole blood to monitor the extent of plasmepsin 2/3 gene amplifications, as well as amplifications in the multidrug resistance transporter 1 gene (pfmdr1), a marker of mefloquine resistance. This study found high concordance across all methods of copy number detection. For samples derived from dried blood spots, a success rate greater than 80% was found in each assay, with more recent samples performing better. Evidence of extensive plasmepsin 2/3 copy number amplifications was observed in Pursat (94%, 2015) (Western Cambodia) and Preah Vihear (87%, 2014) (Northern Cambodia), and lower levels in Ratanakiri (16%, 2014) (Eastern Cambodia). A shift was observed from two copies of plasmepsin 2 in Pursat in 2013 to three copies in 2014-2015 (25% to 64%). Pfmdr1 amplifications were absent in all samples from Preah Vihear and Ratanakiri in 2014 and absent in Pursat in 2015.

CONCLUSIONS

The multiplex TaqMan assay is a robust tool for monitoring both plasmepsin 2/3 and pfmdr1 copy number variations in field isolates, and the SYBR-green and breakpoint assays are useful for monitoring plasmepsin 2/3 amplifications. This study shows increasing levels of plasmepsin 2 copy numbers across Cambodia from 2012 to 2015 and a complete reversion of multicopy pfmdr1 parasites to single copy parasites in all study locations.

Triple artemisinin-based combination therapies versus artemisinin-based combination therapies for uncomplicated Plasmodium falciparum malaria: a multicentre, open-label, randomised clinical trial

BACKGROUND

Artemisinin and partner-drug resistance in Plasmodium falciparum are major threats to malaria control and elimination. Triple artemisinin-based combination therapies (TACTs), which combine existing co-formulated ACTs with a second partner drug that is slowly eliminated, might provide effective treatment and delay emergence of antimalarial drug resistance.

METHODS

In this multicentre, open-label, randomised trial, we recruited patients with uncomplicated P falciparum malaria at 18 hospitals and health clinics in eight countries. Eligible patients were aged 2-65 years, with acute, uncomplicated P falciparum malaria alone or mixed with non-falciparum species, and a temperature of 37·5°C or higher, or a history of fever in the past 24 h. Patients were randomly assigned (1:1) to one of two treatments using block randomisation, depending on their location: in Thailand, Cambodia, Vietnam, and Myanmar patients were assigned to either dihydroartemisinin-piperaquine or dihydroartemisinin-piperaquine plus mefloquine; at three sites in Cambodia they were assigned to either artesunate-mefloquine or dihydroartemisinin-piperaquine plus mefloquine; and in Laos, Myanmar, Bangladesh, India, and the Democratic Republic of the Congo they were assigned to either artemether-lumefantrine or artemether-lumefantrine plus amodiaquine. All drugs were administered orally and doses varied by drug combination and site. Patients were followed-up weekly for 42 days. The primary endpoint was efficacy, defined by 42-day PCR-corrected adequate clinical and parasitological response. Primary analysis was by intention to treat. A detailed assessment of safety and tolerability of the study drugs was done in all patients randomly assigned to treatment. This study is registered at ClinicalTrials.gov, NCT02453308, and is complete.

FINDINGS

Between Aug 7, 2015, and Feb 8, 2018, 1100 patients were given either dihydroartemisinin-piperaquine (183 [17%]), dihydroartemisinin-piperaquine plus mefloquine (269 [24%]), artesunate-mefloquine (73 [7%]), artemether-lumefantrine (289 [26%]), or artemether-lumefantrine plus amodiaquine (286 [26%]). The median age was 23 years (IQR 13 to 34) and 854 (78%) of 1100 patients were male. In Cambodia, Thailand, and Vietnam the 42-day PCR-corrected efficacy after dihydroartemisinin-piperaquine plus mefloquine was 98% (149 of 152; 95% CI 94 to 100) and after dihydroartemisinin-piperaquine was 48% (67 of 141; 95% CI 39 to 56; risk difference 51%, 95% CI 42 to 59; p<0·0001). Efficacy of dihydroartemisinin-piperaquine plus mefloquine in the three sites in Myanmar was 91% (42 of 46; 95% CI 79 to 98) versus 100% (42 of 42; 95% CI 92 to 100) after dihydroartemisinin-piperaquine (risk difference 9%, 95% CI 1 to 17; p=0·12). The 42-day PCR corrected efficacy of dihydroartemisinin-piperaquine plus mefloquine (96% [68 of 71; 95% CI 88 to 99]) was non-inferior to that of artesunate-mefloquine (95% [69 of 73; 95% CI 87 to 99]) in three sites in Cambodia (risk difference 1%; 95% CI -6 to 8; p=1·00). The overall 42-day PCR-corrected efficacy of artemether-lumefantrine plus amodiaquine (98% [281 of 286; 95% CI 97 to 99]) was similar to that of artemether-lumefantrine (97% [279 of 289; 95% CI 94 to 98]; risk difference 2%, 95% CI -1 to 4; p=0·30). Both TACTs were well tolerated, although early vomiting (within 1 h) was more frequent after dihydroartemisinin-piperaquine plus mefloquine (30 [3·8%] of 794) than after dihydroartemisinin-piperaquine (eight [1·5%] of 543; p=0·012). Vomiting after artemether-lumefantrine plus amodiaquine (22 [1·3%] of 1703) and artemether-lumefantrine (11 [0·6%] of 1721) was infrequent. Adding amodiaquine to artemether-lumefantrine extended the electrocardiogram corrected QT interval (mean increase at 52 h compared with baseline of 8·8 ms [SD 18·6] vs 0·9 ms [16·1]; p<0·01) but adding mefloquine to dihydroartemisinin-piperaquine did not (mean increase of 22·1 ms [SD 19·2] for dihydroartemisinin-piperaquine vs 20·8 ms [SD 17·8] for dihydroartemisinin-piperaquine plus mefloquine; p=0·50).

INTERPRETATION

Dihydroartemisinin-piperaquine plus mefloquine and artemether-lumefantrine plus amodiaquine TACTs are efficacious, well tolerated, and safe treatments of uncomplicated P falciparum malaria, including in areas with artemisinin and ACT partner-drug resistance.

FUNDING

UK Department for International Development, Wellcome Trust, Bill & Melinda Gates Foundation, UK Medical Research Council, and US National Institutes of Health.

Interferon-gamma polymorphisms and risk of iron deficiency and anaemia in Gambian children

Background: Anaemia is a major public health concern especially in African children living in malaria-endemic regions. Interferon-gamma (IFN-γ) is elevated during malaria infection and is thought to influence erythropoiesis and iron status. Genetic variants in the IFN-γ gene (IFNG) are associated with increased IFN-γ production. We investigated putative functional single nucleotide polymorphisms (SNPs) and haplotypes of IFNG in relation to nutritional iron status and anaemia in Gambian children over a malaria season. Methods: We used previously available data from Gambian family trios to determine informative SNPs and then used the Agena Bioscience MassArray platform to type five SNPs from the IFNG gene in a cohort of 780 Gambian children. We also measured haemoglobin and biomarkers of iron status and inflammation at the start and end of a malaria season. Results: We identified five IFNG haplotype-tagging SNPs ( IFNG-1616 [rs2069705], IFNG+874 [rs2430561], IFNG+2200 [rs1861493], IFNG+3234 [rs2069718] and IFNG+5612 [rs2069728]). The IFNG+2200C [rs1861493] allele was associated with reduced haemoglobin concentrations (adjusted β -0.44 [95% CI -0.75, -0.12]; Bonferroni adjusted P = 0.03) and a trend towards iron deficiency compared to wild-type at the end of the malaria season in multivariable models adjusted for potential confounders. A haplotype uniquely identified by IFNG+2200C was similarly associated with reduced haemoglobin levels and trends towards iron deficiency, anaemia and iron deficiency anaemia at the end of the malaria season in models adjusted for age, sex, village, inflammation and malaria parasitaemia. Conclusion: We found limited statistical evidence linking IFNG polymorphisms with a risk of developing iron deficiency and anaemia in Gambian children. More definitive studies are needed to investigate the effects of genetically influenced IFN-γ levels on the risk of iron deficiency and anaemia in children living in malaria-endemic areas.

Genomic Analysis of Plasmodium vivax in Southern Ethiopia Reveals Selective Pressures in Multiple Parasite Mechanisms

The Horn of Africa harbors the largest reservoir of Plasmodium vivax in the continent. Most of sub-Saharan Africa has remained relatively vivax-free due to a high prevalence of the human Duffy-negative trait, but the emergence of strains able to invade Duffy-negative reticulocytes poses a major public health threat. We undertook the first population genomic investigation of P. vivax from the region, comparing the genomes of 24 Ethiopian isolates against data from Southeast Asia to identify important local adaptions. The prevalence of the Duffy binding protein amplification in Ethiopia was 79%, potentially reflecting adaptation to Duffy negativity. There was also evidence of selection in a region upstream of the chloroquine resistance transporter, a putative chloroquine-resistance determinant. Strong signals of selection were observed in genes involved in immune evasion and regulation of gene expression, highlighting the need for a multifaceted intervention approach to combat P. vivax in the region.

A high throughput multi-locus insecticide resistance marker panel for tracking resistance emergence and spread in Anopheles gambiae

The spread of resistance to insecticides in disease-carrying mosquitoes poses a threat to the effectiveness of control programmes, which rely largely on insecticide-based interventions. Monitoring mosquito populations is essential, but obtaining phenotypic measurements of resistance is laborious and error-prone. High-throughput genotyping offers the prospect of quick and repeatable estimates of resistance, while also allowing resistance markers to be tracked and studied. To demonstrate the potential of highly-mulitplexed genotypic screening for measuring resistance-association of mutations and tracking their spread, we developed a panel of 28 known or putative resistance markers in the major malaria vector Anopheles gambiae, which we used to screen mosquitoes from a wide swathe of Sub-Saharan Africa (Burkina Faso, Ghana, Democratic Republic of Congo (DRC) and Kenya). We found resistance association in four markers, including a novel mutation in the detoxification gene Gste2 (Gste2-119V). We also identified a duplication in Gste2 combining a resistance-associated mutation with its wild-type counterpart, potentially alleviating the costs of resistance. Finally, we describe the distribution of the multiple origins of kdr resistance, finding unprecedented diversity in the DRC. This panel represents the first step towards a quantitative genotypic model of insecticide resistance that can be used to predict resistance status in An. gambiae.

The ferroportin Q248H mutation protects from anemia, but not malaria or bacteremia

Iron acquisition is critical for life. Ferroportin (FPN) exports iron from mature erythrocytes, and deletion of the Fpn gene results in hemolytic anemia and increased fatality in malaria-infected mice. The FPN Q248H mutation (glutamine to histidine at position 248) renders FPN partially resistant to hepcidin-induced degradation and was associated with protection from malaria in human studies of limited size. Using data from cohorts including over 18,000 African children, we show that the Q248H mutation is associated with modest protection against anemia, hemolysis, and iron deficiency, but we found little evidence of protection against severe malaria or bacteremia. We additionally observed no excess Plasmodium growth in Q248H erythrocytes ex vivo, nor evidence of selection driven by malaria exposure, suggesting that the Q248H mutation does not protect from malaria and is unlikely to deprive malaria parasites of iron essential for their growth.

Evolution and expansion of multidrug-resistant malaria in southeast Asia: a genomic epidemiology study

BACKGROUND

A multidrug-resistant co-lineage of Plasmodium falciparum malaria, named KEL1/PLA1, spread across Cambodia in 2008-13, causing high rates of treatment failure with the frontline combination therapy dihydroartemisinin-piperaquine. Here, we report on the evolution and spread of KEL1/PLA1 in subsequent years.

METHODS

For this genomic epidemiology study, we analysed whole genome sequencing data from P falciparum clinical samples collected from patients with malaria between 2007 and 2018 from Cambodia, Laos, northeastern Thailand, and Vietnam, through the MalariaGEN P falciparum Community Project. Previously unpublished samples were provided by two large-scale multisite projects: the Tracking Artemisinin Resistance Collaboration II (TRAC2) and the Genetic Reconnaissance in the Greater Mekong Subregion (GenRe-Mekong) project. By investigating genome-wide relatedness between parasites, we inferred patterns of shared ancestry in the KEL1/PLA1 population.

FINDINGS

We analysed 1673 whole genome sequences that passed quality filters, and determined KEL1/PLA1 status in 1615. Before 2009, KEL1/PLA1 was only found in western Cambodia; by 2016-17 its prevalence had risen to higher than 50% in all of the surveyed countries except for Laos. In northeastern Thailand and Vietnam, KEL1/PLA1 exceeded 80% of the most recent P falciparum parasites. KEL1/PLA1 parasites maintained high genetic relatedness and low diversity, reflecting a recent common origin. Several subgroups of highly related parasites have recently emerged within this co-lineage, with diverse geographical distributions. The three largest of these subgroups (n=84, n=79, and n=47) mostly emerged since 2016 and were all present in Cambodia, Laos, and Vietnam. These expanding subgroups carried new mutations in the crt gene, which arose on a specific genetic background comprising multiple genomic regions. Four newly emerging crt mutations were rare in the early period and became more prevalent by 2016-17 (Thr93Ser, rising to 19·8%; His97Tyr to 11·2%; Phe145Ile to 5·5%; and Ile218Phe to 11·1%).

INTERPRETATION

After emerging and circulating for several years within Cambodia, the P falciparum KEL1/PLA1 co-lineage diversified into multiple subgroups and acquired new genetic features, including novel crt mutations. These subgroups have rapidly spread into neighbouring countries, suggesting enhanced fitness. These findings highlight the urgent need for elimination of this increasingly drug-resistant parasite co-lineage, and the importance of genetic surveillance in accelerating malaria elimination efforts.

FUNDING

Wellcome Trust, Bill & Melinda Gates Foundation, UK Medical Research Council, and UK Department for International Development.

Determinants of dihydroartemisinin-piperaquine treatment failure in Plasmodium falciparum malaria in Cambodia, Thailand, and Vietnam: a prospective clinical, pharmacological, and genetic study

BACKGROUND

The emergence and spread of resistance in Plasmodium falciparum malaria to artemisinin combination therapies in the Greater Mekong subregion poses a major threat to malaria control and elimination. The current study is part of a multi-country, open-label, randomised clinical trial (TRACII, 2015-18) evaluating the efficacy, safety, and tolerability of triple artemisinin combination therapies. A very high rate of treatment failure after treatment with dihydroartemisinin-piperaquine was observed in Thailand, Cambodia, and Vietnam. The immediate public health importance of our findings prompted us to report the efficacy data on dihydroartemisinin-piperaquine and its determinants ahead of the results of the overall trial, which will be published later this year.

METHODS

Patients aged between 2 and 65 years presenting with uncomplicated P falciparum or mixed species malaria at seven sites in Thailand, Cambodia, and Vietnam were randomly assigned to receive dihydroartemisinin-piperaquine with or without mefloquine, as part of the TRACII trial. The primary outcome was the PCR-corrected efficacy at day 42. Next-generation sequencing was used to assess the prevalence of molecular markers associated with artemisinin resistance (kelch13 mutations, in particular Cys580Tyr) and piperaquine resistance (plasmepsin-2 and plasmepsin-3 amplifications and crt mutations). This study is registered with ClinicalTrials.gov, number NCT02453308.

FINDINGS

Between Sept 28, 2015, and Jan 18, 2018, 539 patients with acute P falciparum malaria were screened for eligibility, 292 were enrolled, and 140 received dihydroartemisinin-piperaquine. The overall Kaplan-Meier estimate of PCR-corrected efficacy of dihydroartemisinin-piperaquine at day 42 was 50·0% (95% CI 41·1-58·3). PCR-corrected efficacies for individual sites were 12·7% (2·2-33·0) in northeastern Thailand, 38·2% (15·9-60·5) in western Cambodia, 73·4% (57·0-84·3) in Ratanakiri (northeastern Cambodia), and 47·1% (33·5-59·6) in Binh Phuoc (southwestern Vietnam). Treatment failure was associated independently with plasmepsin2/3 amplification status and four mutations in the crt gene (Thr93Ser, His97Tyr, Phe145Ile, and Ile218Phe). Compared with the results of our previous TRACI trial in 2011-13, the prevalence of molecular markers of artemisinin resistance (kelch13 Cys580Tyr mutations) and piperaquine resistance (plasmepsin2/3 amplifications and crt mutations) has increased substantially in the Greater Mekong subregion in the past decade.

INTERPRETATION

Dihydroartemisinin-piperaquine is not treating malaria effectively across the eastern Greater Mekong subregion. A highly drug-resistant P falciparum co-lineage is evolving, acquiring new resistance mechanisms, and spreading. Accelerated elimination of P falciparum malaria in this region is needed urgently, to prevent further spread and avoid a potential global health emergency.

FUNDING

UK Department for International Development, Wellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, and National Institutes of Health.

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae, caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so.

The origins of malaria artemisinin resistance defined by a genetic and transcriptomic background

The predisposition of parasites acquiring artemisinin resistance still remains unclear beyond the mutations in Pfk13 gene and modulation of the unfolded protein response pathway. To explore the chain of casualty underlying artemisinin resistance, we reanalyze 773 P. falciparum isolates from TRACI-study integrating TWAS, GWAS, and eQTL analyses. We find the majority of P. falciparum parasites are transcriptomically converged within each geographic site with two broader physiological profiles across the Greater Mekong Subregion (GMS). We report 8720 SNP-expression linkages in the eastern GMS parasites and 4537 in the western. The minimal overlap between them suggests differential gene regulatory networks facilitating parasite adaptations to their unique host environments. Finally, we identify two genetic and physiological backgrounds associating with artemisinin resistance in the GMS, together with a farnesyltransferase protein and a thioredoxin-like protein which may act as vital intermediators linking the Pfk13 C580Y mutation to the prolonged parasite clearance time.

A forward genetic screen reveals a primary role for Plasmodium falciparum Reticulocyte Binding Protein Homologue 2a and 2b in determining alternative erythrocyte invasion pathways

Invasion of human erythrocytes is essential for Plasmodium falciparum parasite survival and pathogenesis, and is also a complex phenotype. While some later steps in invasion appear to be invariant and essential, the earlier steps of recognition are controlled by a series of redundant, and only partially understood, receptor-ligand interactions. Reverse genetic analysis of laboratory adapted strains has identified multiple genes that when deleted can alter invasion, but how the relative contributions of each gene translate to the phenotypes of clinical isolates is far from clear. We used a forward genetic approach to identify genes responsible for variable erythrocyte invasion by phenotyping the parents and progeny of previously generated experimental genetic crosses. Linkage analysis using whole genome sequencing data revealed a single major locus was responsible for the majority of phenotypic variation in two invasion pathways. This locus contained the PfRh2a and PfRh2b genes, members of one of the major invasion ligand gene families, but not widely thought to play such a prominent role in specifying invasion phenotypes. Variation in invasion pathways was linked to significant differences in PfRh2a and PfRh2b expression between parasite lines, and their role in specifying alternative invasion was confirmed by CRISPR-Cas9-mediated genome editing. Expansion of the analysis to a large set of clinical P. falciparum isolates revealed common deletions, suggesting that variation at this locus is a major cause of invasion phenotypic variation in the endemic setting. This work has implications for blood-stage vaccine development and will help inform the design and location of future large-scale studies of invasion in clinical isolates.

Plasmodium parasites cause more than 200 million cases of malaria each year. All the symptoms of malaria are caused after Plasmodium parasites invade human red blood cells. Once inside, they grow, multiply and break open the red blood cells to release new parasites. This cycle is repeated every 48 hours, rapidly amplifying the number of parasites and causing severe anemia and other complications. Plasmodium falciparum, the parasite species responsible for almost all malaria deaths, can use multiple different pathways to invade human red blood cells, but the relative importance of each is not well understood. We tested the invasion pathways used by a collection of closely related parasites and compared their genome sequences to identify the genes responsible. This analysis revealed that expression differences in two neighboring genes of the Reticulocyte Binding Homologue family are responsible for most of the variation in two invasion pathways. P. falciparum may use variation in these genes to avoid the immune system or adapt to specific blood groups, which has important implications for vaccine development against malaria.

Multi-population genomic analysis of malaria parasites indicates local selection and differentiation at the gdv1 locus regulating sexual development

Parasites infect hosts in widely varying environments, encountering diverse challenges for adaptation. To identify malaria parasite genes under locally divergent selection across a large endemic region with a wide spectrum of transmission intensity, genome sequences were obtained from 284 clinical Plasmodium falciparum infections from four newly sampled locations in Senegal, The Gambia, Mali and Guinea. Combining these with previous data from seven other sites in West Africa enabled a multi-population analysis to identify discrete loci under varying local selection. A genome-wide scan showed the most exceptional geographical divergence to be at the early gametocyte gene locus gdv1 which is essential for parasite sexual development and transmission. We identified a major structural dimorphism with alternative 1.5 kb and 1.0 kb sequence deletions at different positions of the 3'-intergenic region, in tight linkage disequilibrium with the most highly differentiated single nucleotide polymorphism, one of the alleles being very frequent in Senegal and The Gambia but rare in the other locations. Long non-coding RNA transcripts were previously shown to include the entire antisense of the gdv1 coding sequence and the portion of the intergenic region with allelic deletions, suggesting adaptive regulation of parasite sexual development and transmission in response to local conditions.

Human candidate gene polymorphisms and risk of severe malaria in children in Kilifi, Kenya: a case-control association study

BACKGROUND

Human genetic factors are important determinants of malaria risk. We investigated associations between multiple candidate polymorphisms-many related to the structure or function of red blood cells-and risk for severe Plasmodium falciparum malaria and its specific phenotypes, including cerebral malaria, severe malaria anaemia, and respiratory distress.

METHODS

We did a case-control study in Kilifi County, Kenya. We recruited as cases children presenting with severe malaria to the high-dependency ward of Kilifi County Hospital. We included as controls infants born in the local community between Aug 1, 2006, and Sept 30, 2010, who were part of a genetics study. We tested for associations between a range of candidate malaria-protective genes and risk for severe malaria and its specific phenotypes. We used a permutation approach to account for multiple comparisons between polymorphisms and severe malaria. We judged p values less than 0·005 significant for the primary analysis of the association between candidate genes and severe malaria.

FINDINGS

Between June 11, 1995, and June 12, 2008, 2244 children with severe malaria were recruited to the study, and 3949 infants were included as controls. Overall, 263 (12%) of 2244 children with severe malaria died in hospital, including 196 (16%) of 1233 with cerebral malaria. We investigated 121 polymorphisms in 70 candidate severe malaria-associated genes. We found significant associations between risk for severe malaria overall and polymorphisms in 15 genes or locations, of which most were related to red blood cells: ABO, ATP2B4, ARL14, CD40LG, FREM3, INPP4B, G6PD, HBA (both HBA1 and HBA2), HBB, IL10, LPHN2 (also known as ADGRL2), LOC727982, RPS6KL1, CAND1, and GNAS. Combined, these genetic associations accounted for 5·2% of the variance in risk for developing severe malaria among individuals in the general population. We confirmed established associations between severe malaria and sickle-cell trait (odds ratio [OR] 0·15, 95% CI 0·11-0·20; p=2·61 × 10-58), blood group O (0·74, 0·66-0·82; p=6·26 × 10-8), and -α3·7-thalassaemia (0·83, 0·76-0·90; p=2·06 × 10-6). We also found strong associations between overall risk of severe malaria and polymorphisms in both ATP2B4 (OR 0·76, 95% CI 0·63-0·92; p=0·001) and FREM3 (0·64, 0·53-0·79; p=3·18 × 10-14). The association with FREM3 could be accounted for by linkage disequilibrium with a complex structural mutation within the glycophorin gene region (comprising GYPA, GYPB, and GYPE) that encodes for the rare Dantu blood group antigen. Heterozygosity for Dantu was associated with risk for severe malaria (OR 0·57, 95% CI 0·49-0·68; p=3·22 × 10-11), as was homozygosity (0·26, 0·11-0·62; p=0·002).

INTERPRETATION

Both ATP2B4 and the Dantu blood group antigen are associated with the structure and function of red blood cells. ATP2B4 codes for plasma membrane calcium-transporting ATPase 4 (the major calcium pump on red blood cells) and the glycophorins are ligands for parasites to invade red blood cells. Future work should aim at uncovering the mechanisms by which these polymorphisms can result in severe malaria protection and investigate the implications of these associations for wider health.

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.

Genomic analysis of a pre-elimination Malaysian Plasmodium vivax population reveals selective pressures and changing transmission dynamics

The incidence of Plasmodium vivax infection has declined markedly in Malaysia over the past decade despite evidence of high-grade chloroquine resistance. Here we investigate the genetic changes in a P. vivax population approaching elimination in 51 isolates from Sabah, Malaysia and compare these with data from 104 isolates from Thailand and 104 isolates from Indonesia. Sabah displays extensive population structure, mirroring that previously seen with the emergence of artemisinin-resistant P. falciparum founder populations in Cambodia. Fifty-four percent of the Sabah isolates have identical genomes, consistent with a rapid clonal expansion. Across Sabah, there is a high prevalence of loci known to be associated with antimalarial drug resistance. Measures of differentiation between the three countries reveal several gene regions under putative selection in Sabah. Our findings highlight important factors pertinent to parasite resurgence and molecular cues that can be used to monitor low-endemic populations at the end stages of P. vivax elimination.

Two complement receptor one alleles have opposing associations with cerebral malaria and interact with α+thalassaemia

Malaria has been a major driving force in the evolution of the human genome. In sub-Saharan African populations, two neighbouring polymorphisms in the Complement Receptor One (CR1) gene, named Sl2 and McCb, occur at high frequencies, consistent with selection by malaria. Previous studies have been inconclusive. Using a large case-control study of severe malaria in Kenyan children and statistical models adjusted for confounders, we estimate the relationship between Sl2 and McCb and malaria phenotypes, and find they have opposing associations. The Sl2 polymorphism is associated with markedly reduced odds of cerebral malaria and death, while the McCb polymorphism is associated with increased odds of cerebral malaria. We also identify an apparent interaction between Sl2 and α+thalassaemia, with the protective association of Sl2 greatest in children with normal α-globin. The complex relationship between these three mutations may explain previous conflicting findings, highlighting the importance of considering genetic interactions in disease-association studies.

Origins of the current outbreak of multidrug-resistant malaria in southeast Asia: a retrospective genetic study

BACKGROUND

Antimalarial resistance is rapidly spreading across parts of southeast Asia where dihydroartemisinin-piperaquine is used as first-line treatment for Plasmodium falciparum malaria. The first published reports about resistance to antimalarial drugs came from western Cambodia in 2013. Here, we analyse genetic changes in the P falciparum population of western Cambodia in the 6 years before those reports.

METHODS

We analysed genome sequence data on 1492 P falciparum samples from 11 locations across southeast Asia, including 464 samples collected in western Cambodia between 2007 and 2013. Different epidemiological origins of resistance were identified by haplotypic analysis of the kelch13 artemisinin resistance locus and the plasmepsin 2-3 piperaquine resistance locus.

FINDINGS

We identified more than 30 independent origins of artemisinin resistance, of which the KEL1 lineage accounted for 140 (91%) of 154 parasites resistant to dihydroartemisinin-piperaquine. In 2008, KEL1 combined with PLA1, the major lineage associated with piperaquine resistance. By 2013, the KEL1/PLA1 co-lineage had reached a frequency of 63% (24/38) in western Cambodia and had spread to northern Cambodia.

INTERPRETATION

The KEL1/PLA1 co-lineage emerged in the same year that dihydroartemisinin-piperaquine became the first-line antimalarial drug in western Cambodia and spread rapidly thereafter, displacing other artemisinin-resistant parasite lineages. These findings have important implications for management of the global health risk associated with the current outbreak of multidrug-resistant malaria in southeast Asia.

FUNDING

Wellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, UK Department for International Development, and the Intramural Research Program of the National Institute of Allergy and Infectious Diseases.

Resistance to malaria through structural variation of red blood cell invasion receptors

The malaria parasite Plasmodium falciparum invades human red blood cells by a series of interactions between host and parasite surface proteins. By analyzing genome sequence data from human populations, including 1269 individuals from sub-Saharan Africa, we identify a diverse array of large copy-number variants affecting the host invasion receptor genes GYPA and GYPB We find that a nearby association with severe malaria is explained by a complex structural rearrangement involving the loss of GYPB and gain of two GYPB-A hybrid genes, which encode a serologically distinct blood group antigen known as Dantu. This variant reduces the risk of severe malaria by 40% and has recently increased in frequency in parts of Kenya, yet it appears to be absent from west Africa. These findings link structural variation of red blood cell invasion receptors with natural resistance to severe malaria.

Massive introgression drives species radiation at the range limit of Anopheles gambiae

Impacts of introgressive hybridisation may range from genomic erosion and species collapse to rapid adaptation and speciation but opportunities to study these dynamics are rare. We investigated the extent, causes and consequences of a hybrid zone between Anopheles coluzzii and Anopheles gambiae in Guinea-Bissau, where high hybridisation rates appear to be stable at least since the 1990s. Anopheles gambiae was genetically partitioned into inland and coastal subpopulations, separated by a central region dominated by A. coluzzii. Surprisingly, whole genome sequencing revealed that the coastal region harbours a hybrid form characterised by an A. gambiae-like sex chromosome and massive introgression of A. coluzzii autosomal alleles. Local selection on chromosomal inversions may play a role in this process, suggesting potential for spatiotemporal stability of the coastal hybrid form and providing resilience against introgression of medically-important loci and traits, found to be more prevalent in inland A. gambiae.

Population genetic structure and adaptation of malaria parasites on the edge of endemic distribution

To determine whether the major human malaria parasite Plasmodium falciparum exhibits fragmented population structure or local adaptation at the northern limit of its African distribution where the dry Sahel zone meets the Sahara, samples were collected from diverse locations within Mauritania over a range of ~1000 km. Microsatellite genotypes were obtained for 203 clinical infection samples from eight locations, and Illumina paired‐end sequences were obtained to yield high coverage genomewide single nucleotide polymorphism (SNP) data for 65 clinical infection samples from four locations. Most infections contained single parasite genotypes, reflecting low rates of transmission and superinfection locally, in contrast to the situation seen in population samples from countries further south. A minority of infections shared related or identical genotypes locally, indicating some repeated transmission of parasite clones without recombination. This caused some multilocus linkage disequilibrium and local divergence, but aside from the effect of repeated genotypes there was minimal differentiation between locations. Several chromosomal regions had elevated integrated haplotype scores (|iHS|) indicating recent selection, including those containing drug resistance genes. A genomewide F_ST scan comparison with previous sequence data from an area in West Africa with higher infection endemicity indicates that regional gene flow prevents genetic isolation, but revealed allele frequency differentiation at three drug resistance loci and an erythrocyte invasion ligand gene. Contrast of extended haplotype signatures revealed none to be unique to Mauritania. Discrete foci of infection on the edge of the Sahara are genetically highly connected to the wider continental parasite population, and local elimination would be difficult to achieve without very substantial reduction in malaria throughout the region.

Plasmodium malariae and P. ovale genomes provide insights into malaria parasite evolution

Elucidation of the evolutionary history and interrelatedness of Plasmodium species that infect humans has been hampered by a lack of genetic information for three human-infective species: P. malariae and two P. ovale species (P. o. curtisi and P. o. wallikeri). These species are prevalent across most regions in which malaria is endemic and are often undetectable by light microscopy, rendering their study in human populations difficult. The exact evolutionary relationship of these species to the other human-infective species has been contested. Using a new reference genome for P. malariae and a manually curated draft P. o. curtisi genome, we are now able to accurately place these species within the Plasmodium phylogeny. Sequencing of a P. malariae relative that infects chimpanzees reveals similar signatures of selection in the P. malariae lineage to another Plasmodium lineage shown to be capable of colonization of both human and chimpanzee hosts. Molecular dating suggests that these host adaptations occurred over similar evolutionary timescales. In addition to the core genome that is conserved between species, differences in gene content can be linked to their specific biology. The genome suggests that P. malariae expresses a family of heterodimeric proteins on its surface that have structural similarities to a protein crucial for invasion of red blood cells. The data presented here provide insight into the evolution of the Plasmodium genus as a whole.

Genetic markers associated with dihydroartemisinin-piperaquine failure in Plasmodium falciparum malaria in Cambodia: a genotype-phenotype association study

BACKGROUND

As the prevalence of artemisinin-resistant Plasmodium falciparum malaria increases in the Greater Mekong subregion, emerging resistance to partner drugs in artemisinin combination therapies seriously threatens global efforts to treat and eliminate this disease. Molecular markers that predict failure of artemisinin combination therapy are urgently needed to monitor the spread of partner drug resistance, and to recommend alternative treatments in southeast Asia and beyond.

METHODS

We did a genome-wide association study of 297 P falciparum isolates from Cambodia to investigate the relationship of 11 630 exonic single-nucleotide polymorphisms (SNPs) and 43 copy number variations (CNVs) with in-vitro piperaquine 50% inhibitory concentrations (IC50s), and tested whether these genetic variants are markers of treatment failure with dihydroartemisinin-piperaquine. We then did a survival analysis of 133 patients to determine whether candidate molecular markers predicted parasite recrudescence following dihydroartemisinin-piperaquine treatment.

FINDINGS

Piperaquine IC50s increased significantly from 2011 to 2013 in three Cambodian provinces (2011 vs 2013 median IC50s: 20·0 nmol/L [IQR 13·7-29·0] vs 39·2 nmol/L [32·8-48·1] for Ratanakiri, 19·3 nmol/L [15·1-26·2] vs 66·2 nmol/L [49·9-83·0] for Preah Vihear, and 19·6 nmol/L [11·9-33·9] vs 81·1 nmol/L [61·3-113·1] for Pursat; all p≤10-3; Kruskal-Wallis test). Genome-wide analysis of SNPs identified a chromosome 13 region that associates with raised piperaquine IC50s. A non-synonymous SNP (encoding a Glu415Gly substitution) in this region, within a gene encoding an exonuclease, associates with parasite recrudescence following dihydroartemisinin-piperaquine treatment. Genome-wide analysis of CNVs revealed that a single copy of the mdr1 gene on chromosome 5 and a novel amplification of the plasmepsin 2 and plasmepsin 3 genes on chromosome 14 also associate with raised piperaquine IC50s. After adjusting for covariates, both exo-E415G and plasmepsin 2-3 markers significantly associate (p=3·0 × 10-8 and p=1·7 × 10-7, respectively) with decreased treatment efficacy (survival rates 0·38 [95% CI 0·25-0·51] and 0·41 [0·28-0·53], respectively).

INTERPRETATION

The exo-E415G SNP and plasmepsin 2-3 amplification are markers of piperaquine resistance and dihydroartemisinin-piperaquine failures in Cambodia, and can help monitor the spread of these phenotypes into other countries of the Greater Mekong subregion, and elucidate the mechanism of piperaquine resistance. Since plasmepsins are involved in the parasite's haemoglobin-to-haemozoin conversion pathway, targeted by related antimalarials, plasmepsin 2-3 amplification probably mediates piperaquine resistance.

FUNDING

Intramural Research Program of the US National Institute of Allergy and Infectious Diseases, National Institutes of Health, Wellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, and UK Department for International Development.

Culture adaptation of malaria parasites selects for convergent loss-of-function mutants

Cultured human pathogens may differ significantly from source populations. To investigate the genetic basis of laboratory adaptation in malaria parasites, clinical Plasmodium falciparum isolates were sampled from patients and cultured in vitro for up to three months. Genome sequence analysis was performed on multiple culture time point samples from six monoclonal isolates, and single nucleotide polymorphism (SNP) variants emerging over time were detected. Out of a total of five positively selected SNPs, four represented nonsense mutations resulting in stop codons, three of these in a single ApiAP2 transcription factor gene, and one in SRPK1. To survey further for nonsense mutants associated with culture, genome sequences of eleven long-term laboratory-adapted parasite strains were examined, revealing four independently acquired nonsense mutations in two other ApiAP2 genes, and five in Epac. No mutants of these genes exist in a large database of parasite sequences from uncultured clinical samples. This implicates putative master regulator genes in which multiple independent stop codon mutations have convergently led to culture adaptation, affecting most laboratory lines of P. falciparum. Understanding the adaptive processes should guide development of experimental models, which could include targeted gene disruption to adapt fastidious malaria parasite species to culture.

Characterisation of the opposing effects of G6PD deficiency on cerebral malaria and severe malarial anaemia

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is believed to confer protection against Plasmodium falciparum malaria, but the precise nature of the protective effecthas proved difficult to define as G6PD deficiency has multiple allelic variants with different effects in males and females, and it has heterogeneous effects on the clinical outcome of P. falciparum infection. Here we report an analysis of multiple allelic forms of G6PD deficiency in a large multi-centre case-control study of severe malaria, using the WHO classification of G6PD mutations to estimate each individual's level of enzyme activity from their genotype. Aggregated across all genotypes, we find that increasing levels of G6PD deficiency are associated with decreasing risk of cerebral malaria, but with increased risk of severe malarial anaemia. Models of balancing selection based on these findings indicate that an evolutionary trade-off between different clinical outcomes of P. falciparum infection could have been a major cause of the high levels of G6PD polymorphism seen in human populations.

THE REAL McCOIL: A method for the concurrent estimation of the complexity of infection and SNP allele frequency for malaria parasites

As many malaria-endemic countries move towards elimination of Plasmodium falciparum, the most virulent human malaria parasite, effective tools for monitoring malaria epidemiology are urgent priorities. P. falciparum population genetic approaches offer promising tools for understanding transmission and spread of the disease, but a high prevalence of multi-clone or polygenomic infections can render estimation of even the most basic parameters, such as allele frequencies, challenging. A previous method, COIL, was developed to estimate complexity of infection (COI) from single nucleotide polymorphism (SNP) data, but relies on monogenomic infections to estimate allele frequencies or requires external allele frequency data which may not available. Estimates limited to monogenomic infections may not be representative, however, and when the average COI is high, they can be difficult or impossible to obtain. Therefore, we developed THE REAL McCOIL, Turning HEterozygous SNP data into Robust Estimates of ALelle frequency, via Markov chain Monte Carlo, and Complexity Of Infection using Likelihood, to incorporate polygenomic samples and simultaneously estimate allele frequency and COI. This approach was tested via simulations then applied to SNP data from cross-sectional surveys performed in three Ugandan sites with varying malaria transmission. We show that THE REAL McCOIL consistently outperforms COIL on simulated data, particularly when most infections are polygenomic. Using field data we show that, unlike with COIL, we can distinguish epidemiologically relevant differences in COI between and within these sites. Surprisingly, for example, we estimated high average COI in a peri-urban subregion with lower transmission intensity, suggesting that many of these cases were imported from surrounding regions with higher transmission intensity. THE REAL McCOIL therefore provides a robust tool for understanding the molecular epidemiology of malaria across transmission settings.