Label-free imaging and classification of live P. falciparum enables high performance parasitemia quantification without fixation or staining.

Manual microscopic inspection of fixed and stained blood smears has remained the gold standard for Plasmodium parasitemia analysis for over a century. Unfortunately, smear preparation consumes time and reagents, while manual microscopy is skill-dependent and labor-intensive. Here, we demonstrate that deep learning enables both life stage classification and accurate parasitemia quantification of ordinary brightfield microscopy images of live, unstained red blood cells. We tested our method using both a standard light microscope equipped with visible and near-ultraviolet (UV) illumination, and a custom-built microscope employing deep-UV illumination. While using deep-UV light achieved an overall four-category classification of Plasmodium falciparum blood stages of greater than 99% and a recall of 89.8% for ring-stage parasites, imaging with near-UV light on a standard microscope resulted in 96.8% overall accuracy and over 90% recall for ring-stage parasites. Both imaging systems were tested extrinsically by parasitemia titration, revealing superior performance over manually-scored Giemsa-stained smears, and a limit of detection below 0.1%. Our results establish that label-free parasitemia analysis of live cells is possible in a biomedical laboratory setting without the need for complex optical instrumentation. We anticipate future extensions of this work could enable label-free clinical diagnostic measurements, one day eliminating the need for conventional blood smear analysis.

Nonparametric Binary Classification to Distinguish Closely Related versus Unrelated Plasmodium falciparum Parasites.

Assessing genetic relatedness of Plasmodium falciparum genotypes is a key component of antimalarial efficacy trials. Previous methods have focused on determining a priori definitions of the level of genetic similarity sufficient to classify two infections as sharing the same strain. However, factors such as mixed-strain infections, allelic suppression, imprecise typing methods, and heterozygosity complicate comparisons of apicomplexan genotypes. Here, we introduce a novel method for nonparametric statistical testing of relatedness for P. falciparum parasites. First, the background distribution of genetic distance between unrelated strains is computed. Second, a threshold genetic distance is computed from this empiric distribution of distances to demarcate genetic distances that are unlikely to have arisen by chance. Third, the genetic distance between paired samples is computed, and paired samples with genetic distances below the threshold are classified as related. The method is designed to work with any arbitrary genetic distance definition. We validated this procedure using two independent approaches to calculating genetic distance. We assessed the concordance of the novel nonparametric classification with a gold-standard Bayesian approach for 175 pairs of recurrent P. falciparum episodes from previously published malaria efficacy trials with microsatellite data from five studies in Guinea and Angola. The novel nonparametric approach was 98% sensitive and 84-89% specific in correctly identifying related genotypes compared with a gold-standard Bayesian algorithm. The approach provides a unified and systematic method to statistically assess relatedness of P. falciparum parasites using arbitrary genetic distance methodologies.

Genotyping cognate Plasmodium falciparum in humans and mosquitoes to estimate onward transmission of asymptomatic infections.

Malaria control may be enhanced by targeting reservoirs of Plasmodium falciparum transmission. One putative reservoir is asymptomatic malaria infections and the scale of their contribution to transmission in natural settings is not known. We assess the contribution of asymptomatic malaria to onward transmission using a 14-month longitudinal cohort of 239 participants in a high transmission site in Western Kenya. We identify P. falciparum in asymptomatically- and symptomatically-infected participants and naturally-fed mosquitoes from their households, genotype all parasites using deep sequencing of the parasite genes pfama1 and pfcsp, and use haplotypes to infer participant-to-mosquito transmission through a probabilistic model. In 1,242 infections (1,039 in people and 203 in mosquitoes), we observe 229 (pfcsp) and 348 (pfama1) unique parasite haplotypes. Using these to link human and mosquito infections, compared with symptomatic infections, asymptomatic infections more than double the odds of transmission to a mosquito among people with both infection types (Odds Ratio: 2.56; 95% Confidence Interval (CI): 1.36-4.81) and among all participants (OR 2.66; 95% CI: 2.05-3.47). Overall, 94.6% (95% CI: 93.1-95.8%) of mosquito infections likely resulted from asymptomatic infections. In high transmission areas, asymptomatic infections are the major contributor to mosquito infections and may be targeted as a component of transmission reduction.

Development of a new barcode-based, multiplex-PCR, next-generation-sequencing assay and data processing and analytical pipeline for multiplicity of infection detection of Plasmodium falciparum.

BACKGROUND: Simultaneous infection with multiple malaria parasite strains is common in high transmission areas. Quantifying the number of strains per host, or the multiplicity of infection (MOI), provides additional parasite indices for assessing transmission levels but it is challenging to measure accurately with current tools. This paper presents new laboratory and analytical methods for estimating the MOI of Plasmodium falciparum. METHODS: Based on 24 single nucleotide polymorphisms (SNPs) previously identified as stable, unlinked targets across 12 of the 14 chromosomes within P. falciparum genome, three multiplex PCRs of short target regions and subsequent next generation sequencing (NGS) of the amplicons were developed. A bioinformatics pipeline including B4Screening pathway removed spurious amplicons to ensure consistent frequency calls at each SNP location, compiled amplicons by SNP site diversity, and performed algorithmic haplotype and strain reconstruction. The pipeline was validated by 108 samples generated from cultured-laboratory strain mixtures in different proportions and concentrations, with and without pre-amplification, and using whole blood and dried blood spots (DBS). The pipeline was applied to 273 smear-positive samples from surveys conducted in western Kenya, then providing results into StrainRecon Thresholding for Infection Multiplicity (STIM), a novel MOI estimator. RESULTS: The 24 barcode SNPs were successfully identified uniformly across the 12 chromosomes of P. falciparum in a sample using the pipeline. Pre-amplification and parasite concentration, while non-linearly associated with SNP read depth, did not influence the SNP frequency calls. Based on consistent SNP frequency calls at targeted locations, the algorithmic strain reconstruction for each laboratory-mixed sample had 98.5% accuracy in dominant strains. STIM detected up to 5 strains in field samples from western Kenya and showed declining MOI over time (q < 0.02), from 4.32 strains per infected person in 1996 to 4.01, 3.56 and 3.35 in 2001, 2007 and 2012, and a reduction in the proportion of samples with 5 strains from 57% in 1996 to 18% in 2012. CONCLUSION: The combined approach of new multiplex PCRs and NGS, the unique bioinformatics pipeline and STIM could identify 24 barcode SNPs of P. falciparum correctly and consistently. The methodology could be applied to field samples to reliably measure temporal changes in MOI.

Case Report: The First Case of Genotypically Confirmed Plasmodium falciparum Kelch 13 Propeller Mutation in Sri Lanka and Its Implications on the Elimination Status of Malaria.

This case report discusses recrudescence of imported Plasmodium falciparum malaria, in the presence of P. falciparum Kelch13 (PfK13) propeller mutation, in a patient diagnosed and fully treated with artemether-lumefantrine under direct observation in Sri Lanka. This patient presented with a history of 5 days of fever following his arrival from the Democratic Republic of Congo (DRC). He had visited Rwanda 1 week before arrival to Sri Lanka. Treatment was commenced with artemisinin-based combination therapy, artemether-lumefantrine, which is the first-line drug recommended for uncomplicated falciparum malaria. Blood smears were negative for parasites by the third day of treatment. Approximately 2 weeks later, he developed fever again and was diagnosed as having a recrudescence of falciparum malaria. He was treated and responded to the second-line antimalarial dihydroartemisinin-piperaquine. Molecular testing of blood taken from the first infection revealed the presence of amino acid substitutions K189T and R561H within the PfK13 gene. R561H mutation is associated with delayed parasite clearance in Southeast Asia. Although seldom reported from DRC, an emergence and clonal expansion of parasites harboring R561H allele has been reported from Rwanda recently; thus, it is likely that this patient may have got the infection from Rwanda. Sri Lanka eliminated malaria in 2016. However, in the backdrop of continuing imported malaria cases, early diagnosis and prompt treatment is essential to prevent the re-establishment of the disease.

Genetic Markers of Adaptation of Plasmodium falciparum to Transmission by American Vectors Identified in the Genomes of Parasites from Haiti and South America.

The malaria parasite, Plasmodium falciparum, was introduced into Hispaniola and other regions of the Americas through the slave trade spanning the 16th through the 19th centuries. During this period, more than 12 million Africans were brought across the Atlantic to the Caribbean and other regions of the Americas. Since malaria is holoendemic in West Africa, a substantial percentage of these individuals carried the parasite. St. Domingue on Hispaniola, now modern-day Haiti, was a major port of disembarkation, and malaria is still actively transmitted there. We undertook a detailed study of the phylogenetics of the Haitian parasites and those from Colombia and Peru utilizing whole-genome sequencing. Principal-component and phylogenetic analyses, based upon single nucleotide polymorphisms (SNPs) in protein coding regions, indicate that, despite the potential for millions of introductions from Africa, the Haitian parasites share an ancestral relationship within a well-supported monophyletic clade with parasites from South America, while belonging to a distinct lineage. This result, in stark contrast to the historical record of parasite introductions, is best explained by a severe population bottleneck experienced by the parasites introduced into the Americas. Here, evidence is presented for targeted selection of rare African alleles in genes which are expressed in the mosquito stages of the parasite's life cycle. These genetic markers support the hypothesis that the severe population bottleneck was caused by the required adaptation of the parasite to transmission by new definitive hosts among the Anopheles (Nyssorhynchus) spp. found in the Caribbean and South America.IMPORTANCE Historical data suggest that millions of P. falciparum parasite lineages were introduced into the Americas during the trans-Atlantic slave trade, which would suggest a paraphyletic origin of the extant isolates in the Western Hemisphere. Our analyses of whole-genome variants show that the American parasites belong to a well-supported monophyletic clade. We hypothesize that the required adaptation to American vectors created a severe bottleneck, reducing the effective introduction to a few lineages. In support of this hypothesis, we discovered genes expressed in the mosquito stages of the life cycle that have alleles with multiple, high-frequency or fixed, nonsynonymous mutations in the American populations which are rarely found in African isolates. These alleles appear to be in gene products critical for transmission through the anopheline vector. Thus, these results may inform efforts to develop novel transmission-blocking vaccines by identifying parasite proteins functionally interacting with the vector that are important for successful transmission. Further, to the best of our knowledge, these are the first whole-genome data available from Haitian P. falciparum isolates. Defining the genome of these parasites provides genetic markers useful for mapping parasite populations and monitoring parasite movements/introductions.

Misdiagnosis of imported falciparum malaria from African areas due to an increased prevalence of pfhrp2/pfhrp3 gene deletion: the Djibouti case.

Following the diagnosis of a falciparum malaria case imported from Djibouti and not detected by a pfHRP2-based rapid diagnostic test (RDT), we investigated the prevalence of the pfhrp2/pfhrp3-deleted parasites in Djibouti using 378 blood samples collected between January and May 2019, from Djiboutian patients with suspected malaria. Malaria diagnosis by quantitative PCR confirmed the presence of Plasmodium falciparum for 20.9% (79/378) samples while RDTs did not detect HRP2 antigen in 83.5% (66/79) of these samples. Quantitative PCRs targeting the pfhrp2/pfhrp3 genes confirmed the absence of both genes for 86.5% of P. falciparum strains. The very large number (86.5%) of falciparum parasites lacking the pfhrp2/pfhrp3 genes observed in this study, now justifies the use of non-HRP2 alternative RDTs in Djibouti. In this area and in most countries where HRP2-based RDTs constitute the main arsenal for falciparum malaria diagnosis, it is important to implement a systematic surveillance and to inform biologists and clinicians about the risk of malaria misdiagnosis. Further investigations are needed to better understand the mechanism of selection and diffusion of the pfhrp2/pfhrp3-deleted parasites.

A modified two-color flow cytometry assay to quantify in-vitro reinvasion and determine invasion phenotypes at low Plasmodium falciparum parasitemia.

Invasion of human red blood cells (RBCs) by Plasmodium parasites is a crucial yet poorly characterised phenotype. Two-color flow cytometry (2cFCM) promises to be a very sensitive and high throughput method for phenotyping parasite invasion. However, current protocols require high (~1.0%) parasitemia for assay set-up and need to be adapted for low parasitemia samples, which are becoming increasingly common in low transmission settings. Background fluorescence from nuclei-containing uninfected RBCs and high autologous reinvasion rates (merozoite invasion of donor uninfected RBCs present at 50% assay volume) are some of the limitations to the method's sensitivity to enumerate low parasitemia (<0.5%) with nucleic acid-based stains. Here, we describe modifications for plating unlabeled donor to labeled target RBCs per assay well and for gating parasitemia, that produces accurate quantifications of low reinvasion parasitemia. Plasmodium falciparum 3D7, Dd2 and field isolates at various low and high parasitemia (0.05%-2.0%) were used to set-up SyBr Green 1-based 2cFCM invasion assays. Target RBCs were labeled with CTFR proliferation dye. We show that this dye combination allowed for efficient parasite invasion into target RBCs and that a 1:3 ratio of unlabeled to labeled RBCs per assay greatly skewed autologous reinvasion (p < 0.001). Accuracy of quantifying reinvasion was limited to an assay parasitemia of 0.02% with minimal background interference. Invasion inhibition by enzymatic treatments increased averagely by 10% (p<0.05) across the entire parasitemia range. The effect was greater for samples with <0.5% parasitemia. Overall, a more sensitive method for phenotyping invasion of low P. falciparum parasitemia is described.

Refining the transcriptome of the human malaria parasite Plasmodium falciparum using amplification-free RNA-seq.

BACKGROUND: Plasmodium parasites undergo several major developmental transitions during their complex lifecycle, which are enabled by precisely ordered gene expression programs. Transcriptomes from the 48-h blood stages of the major human malaria parasite Plasmodium falciparum have been described using cDNA microarrays and RNA-seq, but these assays have not always performed well within non-coding regions, where the AT-content is often 90-95%. RESULTS: We developed a directional, amplification-free RNA-seq protocol (DAFT-seq) to reduce bias against AT-rich cDNA, which we have applied to three strains of P. falciparum (3D7, HB3 and IT). While strain-specific differences were detected, overall there is strong conservation between the transcriptional profiles. For the 3D7 reference strain, transcription was detected from 89% of the genome, with over 78% of the genome transcribed into mRNAs. We also find that transcription from bidirectional promoters frequently results in non-coding, antisense transcripts. These datasets allowed us to refine the 5' and 3' untranslated regions (UTRs), which can be variable, long (> 1000 nt), and often overlap those of adjacent transcripts. CONCLUSIONS: The approaches applied in this study allow a refined description of the transcriptional landscape of P. falciparum and demonstrate that very little of the densely packed P. falciparum genome is inactive or redundant. By capturing the 5' and 3' ends of mRNAs, we reveal both constant and dynamic use of transcriptional start sites across the intraerythrocytic developmental cycle that will be useful in guiding the definition of regulatory regions for use in future experimental gene expression studies.

Efficacy and safety of artemisinin-based combination therapy and the implications of Pfkelch13 and Pfcoronin molecular markers in treatment failure in Senegal.

In 2006, Senegal adopted artemisinin-based combination therapy (ACT) as first-line treatment in the management of uncomplicated malaria. This study aimed to update the status of antimalarial efficacy more than ten years after their first introduction. This was a randomized, three-arm, open-label study to evaluate the efficacy and safety of artemether-lumefantrine (AL), artesunate-amodiaquine (ASAQ) and dihydroartemisinin-piperaquine (DP) in Senegal. Malaria suspected patients were screened, enrolled, treated, and followed for 28 days for AL and ASAQ arms or 42 days for DP arm. Clinical and parasitological responses were assessed following antimalarial treatment. Genotyping (msp1, msp2 and 24 SNP-based barcode) were done to differentiate recrudescence from re-infection; in case of PCR-confirmed treatment failure, Pfk13 propeller and Pfcoronin genes were sequenced. Data was entered and analyzed using the WHO Excel-based application. A total of 496 patients were enrolled. In Diourbel, PCR non-corrected/corrected adequate clinical and parasitological responses (ACPR) was 100.0% in both the AL and ASAQ arms. In Kedougou, PCR corrected ACPR values were 98.8%, 100% and 97.6% in AL, ASAQ and DP arms respectively. No Pfk13 or Pfcoronin mutations associated with artemisinin resistance were found. This study showed that AL, ASAQ and DP remain efficacious and well-tolerated in the treatment of uncomplicated P. falciparum malaria in Senegal.

Molecular surveillance over 14 years confirms reduction of Plasmodium vivax and falciparum transmission after implementation of Artemisinin-based combination therapy in Papua, Indonesia.

Genetic epidemiology can provide important insights into parasite transmission that can inform public health interventions. The current study compared long-term changes in the genetic diversity and structure of co-endemic Plasmodium falciparum and P. vivax populations. The study was conducted in Papua Indonesia, where high-grade chloroquine resistance in P. falciparum and P. vivax led to a universal policy of Artemisinin-based Combination Therapy (ACT) in 2006. Microsatellite typing and population genetic analyses were undertaken on available isolates collected between 2004 and 2017 from patients with uncomplicated malaria (n = 666 P. falciparum and n = 615 P. vivax). The proportion of polyclonal P. falciparum infections fell from 28% (38/135) before policy change (2004-2006) to 18% (22/125) at the end of the study (2015-2017); p<0.001. Over the same period, polyclonal P. vivax infections fell from 67% (80/119) to 35% (33/93); p<0.001. P. falciparum strains persisted for up to 9 years compared to 3 months for P. vivax, reflecting higher rates of outbreeding in the latter. Sub-structure was observed in the P. falciparum population, but not in P. vivax, confirming different patterns of outbreeding. The P. falciparum population exhibited 4 subpopulations that changed in frequency over time. Notably, a sharp rise was observed in the frequency of a minor subpopulation (K2) in the late post-ACT period, accounting for 100% of infections in late 2016-2017. The results confirm epidemiological evidence of reduced P. falciparum and P. vivax transmission over time. The smaller change in P. vivax population structure is consistent with greater outbreeding associated with relapsing infections and highlights the need for radical cure to reduce recurrent infections. The study emphasizes the challenge in disrupting P. vivax transmission and demonstrates the potential of molecular data to inform on the impact of public health interventions.

Genetic diversity of Plasmodium vivax and Plasmodium falciparum lactate dehydrogenases in Myanmar isolates.

BACKGROUND: Plasmodium lactate dehydrogenase (pLDH) is a major target in diagnosing the erythrocytic stage of malaria parasites because it is highly expressed during blood-stage parasites and is distinguished from human LDH. Rapid diagnostic tests (RDTs) for malaria use pLDH as a target antigen; however, genetic variations in pLDH within the natural population threaten the efficacy of pLDH-based RDTs. METHODS: Genetic polymorphisms of Plasmodium vivax LDH (PvLDH) and Plasmodium falciparum LDH (PfLDH) in Myanmar isolates were analysed by nucleotide sequencing analysis. Genetic polymorphisms and the natural selection of PvLDH and PfLDH were analysed using DNASTAR, MEGA6, and DnaSP ver. 5.10.00 programs. The genetic diversity and natural selection of global PvLDH and PfLDH were also analysed. The haplotype network of global PvLDH and PfLDH was constructed using NETWORK ver. 5.0.0.3. Three-dimensional structures of PvLDH and PfLDH were built with YASARA Structure ver. 18.4.24 and the impact of mutations on structural change and stability was evaluated with SDM ver. 2, CUPSAT and MAESTROweb. RESULTS: Forty-nine PvLDH and 52 PfLDH sequences were obtained from Myanmar P. vivax and P. falciparum isolates. Non-synonymous nucleotide substitutions resulting in amino acid changes were identified in both Myanmar PvLDH and PfLDH. Amino acid changes were also identified in the global PvLDH and PfLDH populations, but they did not produce structural alterations in either protein. Low genetic diversity was observed in global PvLDH and PfLDH, which may be maintained by a strong purifying selection. CONCLUSION: This study extends knowledge for genetic diversity and natural selection of global PvLDH and PfLDH. Although amino acid changes were observed in global PvLDH and PfLDH, they did not alter the conformational structures of the proteins. These suggest that PvLDH and PfLDH are genetically well-conserved in global populations, which indicates that they are suitable antigens for diagnostic purpose and attractive targets for drug development.

Falciparum malaria from coastal Tanzania and Zanzibar remains highly connected despite effective control efforts on the archipelago.

BACKGROUND: Tanzania's Zanzibar archipelago has made significant gains in malaria control over the last decade and is a target for malaria elimination. Despite consistent implementation of effective tools since 2002, elimination has not been achieved. Importation of parasites from outside of the archipelago is thought to be an important cause of malaria's persistence, but this paradigm has not been studied using modern genetic tools. METHODS: Whole-genome sequencing (WGS) was used to investigate the impact of importation, employing population genetic analyses of Plasmodium falciparum isolates from both the archipelago and mainland Tanzania. Ancestry, levels of genetic diversity and differentiation, patterns of relatedness, and patterns of selection between these two populations were assessed by leveraging recent advances in deconvolution of genomes from polyclonal malaria infections. RESULTS: Significant decreases in the effective population sizes were inferred in both populations that coincide with a period of decreasing malaria transmission in Tanzania. Identity by descent analysis showed that parasites in the two populations shared long segments of their genomes, on the order of 5 cM, suggesting shared ancestry within the last 10 generations. Even with limited sampling, two of isolates between the mainland and Zanzibar were identified that are related at the expected level of half-siblings, consistent with recent importation. CONCLUSIONS: These findings suggest that importation plays an important role for malaria incidence on Zanzibar and demonstrate the value of genomic approaches for identifying corridors of parasite movement to the island.

PCR correction strategies for malaria drug trials: updates and clarifications.

Malaria drug trials conducted in endemic areas face a major challenge in their analysis because it is difficult to establish whether parasitaemia in blood samples collected after treatment indicate drug failure or a new infection acquired after treatment. It is therefore vital to reliably distinguish drug failures from new infections in order to obtain accurate estimates of drug failure rates. This distinction can be achieved for Plasmodium falciparum by comparing parasite genotypes obtained at the time of treatment (the baseline) and on the day of recurring parasitaemia. Such PCR correction is required to obtain accurate failure rates, even for new effective drugs. Despite the routine use of PCR correction in surveillance of drug resistance and in clinical drug trials, limitations inherent to the molecular genotyping methods have led some researchers to question the validity of current PCR correction strategies. Here we describe and discuss recent developments in these genotyping approaches, with a particular focus on method validation and limitations of the genotyping strategies. Our aim is to update scientists from public and private bodies who are working on the development, deployment, and surveillance of new malaria drugs. We aim to promote discussion around these issues and argue for the adoption of improved standardised PCR correction methodologies.

Resurrection of the ancestral RH5 invasion ligand provides a molecular explanation for the origin of P. falciparum malaria in humans.

Many important infectious diseases are the result of zoonoses, in which pathogens that normally infect animals acquire mutations that enable the breaching of species barriers to permit the infection of humans. Our understanding of the molecular events that enable host switching are often limited, and yet this is a fundamentally important question. Plasmodium falciparum, the etiological agent of severe human malaria, evolved following a zoonotic transfer of parasites from gorillas. One gene-rh5-which encodes an essential ligand for the invasion of host erythrocytes, is suspected to have played a critical role in this host switch. Genome comparisons revealed an introgressed sequence in the ancestor of P. falciparum containing rh5, which likely allowed the ancestral parasites to infect both gorilla and human erythrocytes. To test this hypothesis, we resurrected the ancestral introgressed reticulocyte-binding protein homologue 5 (RH5) sequence and used quantitative protein interaction assays to demonstrate that this ancestral protein could bind the basigin receptor from both humans and gorillas. We also showed that this promiscuous receptor binding phenotype of RH5 was shared with the parasite clade that transferred its genome segment to the ancestor of P. falciparum, while the other lineages exhibit host-specific receptor binding, confirming the central importance of this introgression event for Plasmodium host switching. Finally, since its transfer to humans, P. falciparum, and also the RH5 ligand, have evolved a strong human specificity. We show that this subsequent restriction to humans can be attributed to a single amino acid mutation in the RH5 sequence. Our findings reveal a molecular pathway for the origin and evolution of human P. falciparum malaria and may inform molecular surveillance to predict future zoonoses.

Prevalence of Plasmodium falciparum Pfcrt and Pfmdr1 alleles in settings with different levels of Plasmodium vivax co-endemicity in Ethiopia.

Plasmodium falciparum and P. vivax co-exist at different endemicity levels across Ethiopia. For over two decades Artemether-Lumefantrine (AL) is the first line treatment for uncomplicated P. falciparum, while chloroquine (CQ) is still used to treat P. vivax. It is currently unclear whether a shift from CQ to AL for P. falciparum treatment has implications for AL efficacy and results in a reversal of mutations in genes associated to CQ resistance, given the high co-endemicity of the two species and the continued availability of CQ for the treatment of P. vivax. This study thus assessed the prevalence of Pfcrt-K76T and Pfmdr1-N86Y point mutations in P. falciparum. 18S RNA gene based nested PCR confirmed P. falciparum samples (N = 183) collected through community and health facility targeted cross-sectional surveys from settings with varying P. vivax and P. falciparum endemicity were used. The proportion of Plasmodium infections that were P. vivax was 62.2% in Adama, 41.4% in Babile, 30.0% in Benishangul-Gumuz to 6.9% in Gambella. The Pfcrt-76T mutant haplotype was observed more from samples with higher endemicity of P. vivax as being 98.4% (61/62), 100% (31/31), 65.2% (15/23) and 41.5% (22/53) in samples from Adama, Babile, Benishangul-Gumuz and Gambella, respectively. However, a relatively higher proportion of Pfmdr1-N86 allele (77.3-100%) were maintained in all sites. The observed high level of the mutant Pfcrt-76T allele in P. vivax co-endemic sites might require that utilization of CQ needs to be re-evaluated in settings co-endemic for the two species. A country-wide assessment is recommended to clarify the implication of the observed level of variation in drug resistance markers on the efficacy of AL-based treatment against uncomplicated P. falciparum malaria.

The origins and relatedness structure of mixed infections vary with local prevalence of P. falciparum malaria.

Individual malaria infections can carry multiple strains of Plasmodium falciparum with varying levels of relatedness. Yet, how local epidemiology affects the properties of such mixed infections remains unclear. Here, we develop an enhanced method for strain deconvolution from genome sequencing data, which estimates the number of strains, their proportions, identity-by-descent (IBD) profiles and individual haplotypes. Applying it to the Pf3k data set, we find that the rate of mixed infection varies from 29% to 63% across countries and that 51% of mixed infections involve more than two strains. Furthermore, we estimate that 47% of symptomatic dual infections contain sibling strains likely to have been co-transmitted from a single mosquito, and find evidence of mixed infections propagated over successive infection cycles. Finally, leveraging data from the Malaria Atlas Project, we find that prevalence correlates within Africa, but not Asia, with both the rate of mixed infection and the level of IBD.

Estimating Relatedness Between Malaria Parasites.

Understanding the relatedness of individuals within or between populations is a common goal in biology. Increasingly, relatedness features in genetic epidemiology studies of pathogens. These studies are relatively new compared to those in humans and other organisms, but are important for designing interventions and understanding pathogen transmission. Only recently have researchers begun to routinely apply relatedness to apicomplexan eukaryotic malaria parasites, and to date have used a range of different approaches on an ad hoc basis. Therefore, it remains unclear how to compare different studies and which measures to use. Here, we systematically compare measures based on identity-by-state (IBS) and identity-by-descent (IBD) using a globally diverse data set of malaria parasites, Plasmodium falciparum and P. vivax, and provide marker requirements for estimates based on IBD. We formally show that the informativeness of polyallelic markers for relatedness inference is maximized when alleles are equifrequent. Estimates based on IBS are sensitive to allele frequencies, which vary across populations and by experimental design. For portability across studies, we thus recommend estimates based on IBD. To generate estimates with errors below an arbitrary threshold of 0.1, we recommend ∼100 polyallelic or 200 biallelic markers. Marker requirements are immediately applicable to haploid malaria parasites and other haploid eukaryotes. C.I.s facilitate comparison when different marker sets are used. This is the first attempt to provide rigorous analysis of the reliability of, and requirements for, relatedness inference in malaria genetic epidemiology. We hope it will provide a basis for statistically informed prospective study design and surveillance strategies.

Impact of Multiplicity of Plasmodium falciparum Infection on Clinical Disease in Malawi.

Multiplicity of infection (MOI), the number of unique Plasmodium falciparum parasite genotypes found in one infected individual, may contribute to the development of clinical malaria disease. However, the independent contribution of MOI and parasite density to clinical disease has not been well characterized. We conducted a two-year longitudinal cohort study of adults and children in a high-transmission setting in Malawi to test the hypothesis that increased MOI was independently associated with clinical disease, after accounting for parasite density. Of 1,062 episodes of infection, 477 (44.9%) were associated with symptoms. After controlling for repeated measures within an individual, key demographic factors, and parasite density, there was no association between MOI and clinical disease (OR = 1.02, 95% CI: 0.70-1.51). Although the limited ability to discern MOI in low-density asymptomatic infections may have impacted our results, we conclude that MOI is not an independent risk factor for clinical disease.

Evolutionary sex allocation theory explains sex ratios in natural Plasmodium falciparum infections.

Malaria transmission is achieved by sexual stages, called gametocytes, and the proportion of gametocytes that are male versus female (sex ratio) influences transmission success. In malaria model systems, variation in gametocyte sex ratios can be explained by the predictions of evolutionary sex allocation theory. We test these predictions using natural Plasmodium falciparum infections. The predicted negative correlation between sex ratio and gametocyte density holds: the sex ratio increases when gametocyte densities decrease, and this is most apparent in single genotype infections and in the dry season. We do not observe higher gametocyte sex ratios in mixed compared with single genotype infections.