Malaria is a major public health problem that is actively being addressed in a global eradication campaign. Increased population mobility through international air travel has elevated the risk of re-introducing parasites to elimination areas and dispersing drug-resistant parasites to new regions. A simple genetic marker that quickly and accurately identifies the geographic origin of infections would be a valuable public health tool for locating the source of imported outbreaks. Here we analyse the mitochondrion and apicoplast genomes of 711 Plasmodium falciparum isolates from 14 countries, and find evidence that they are non-recombining and co-inherited. The high degree of linkage produces a panel of relatively few single-nucleotide polymorphisms (SNPs) that is geographically informative. We design a 23-SNP barcode that is highly predictive (~92%) and easily adapted to aid case management in the field and survey parasite migration worldwide.

Previous genome-wide analyses of single nucleotide variation in Plasmodium falciparum identified evidence of an extended haplotype region on chromosome 6 in West Africa, suggesting recent positive selection. Such a pattern is not seen in samples from East Africa or South East Asia, so it could be marking a selective process specific to West Africa. Analyses of the haplotype structure in samples taken at different times could give clues to possible causes of selection.

Many studies report associations between human genetic factors and immunity to malaria but few have been reliably replicated. These studies are usually country-specific, use small sample sizes and are not directly comparable due to differences in methodologies. This study brings together samples and data collected from multiple sites across Africa and Asia to use standardized methods to look for consistent genetic effects on anti-malarial antibody levels.

It is important to understand intrinsic variation in asexual blood stage multiplication rates of the most virulent human malaria parasite, Plasmodium falciparum. Here, multiplication rates of long-term laboratory adapted parasite clones and new clinical isolates were measured, using a newly standardised assay of growth from low starting density in replicate parallel cultures with erythrocytes from multiple different donors, across multiple cycles. Multiplication rates of long-term established clones were between 7.6 and 10.5 fold per 48 hours, with clone Dd2 having a higher rate than others (clones 3D7, HB3 and D10). Parasite clone-specific growth was then analysed in co-culture assays with all possible heterologous pairwise combinations. This showed that co-culture of different parasites did not affect their replication rates, indicating that there were no suppressive interactions operating between parasites. Multiplication rates of eleven new clinical isolates were measured after a few weeks of culture, and showed a spectrum of replication rates between 2.3 and 6.0 fold per 48 hours, the entire range being lower than for the long-term laboratory adapted clones. Multiplication rate estimates remained stable over time for several isolates tested repeatedly up to three months after culture initiation, indicating considerable persistence of this important trait variation.

Asymptomatic malaria parasites are significant sources of infections for onward malaria transmission. Conventional tools for malaria diagnosis such as microscopy and rapid diagnostic test kits (RDT) have relatively low sensitivity, hence the need for alternative tools for active screening of such low-density infections.

Selection pressure from continued exposure to insecticides drives development of insecticide resistance and changes in resting behaviour of malaria vectors. There is need to understand how resistance drives changes in resting behaviour within vector species. The association between insecticide resistance and resting behaviour of Anopheles gambiae sensu lato (s.l.) in Northern Ghana was examined.

The spread of drug resistance to antimalarial treatments poses a serious public health risk globally. To combat this risk, molecular surveillance of drug resistance is imperative. We report the prevalence of mutations in the Plasmodium falciparum kelch 13 propeller domain associated with partial artemisinin resistance, which we determined by using Sanger sequencing samples from patients enrolled in therapeutic efficacy studies from 9 sub-Saharan countries during 2014-2018. Of the 2,865 samples successfully sequenced before treatment (day of enrollment) and on the day of treatment failure, 29 (1.0%) samples contained 11 unique nonsynonymous mutations and 83 (2.9%) samples contained 27 unique synonymous mutations. Two samples from Kenya contained the S522C mutation, which has been associated with delayed parasite clearance; however, no samples contained validated or candidate artemisinin-resistance mutations.

The malaria parasite Plasmodium falciparum utilizes multiple alternative receptor-ligand interactions for the invasion of human erythrocytes. While some P. falciparum clones make use of sialic acid (SA) residues on the surface of the human glycophorin receptors to invade the erythrocyte, others use alternative receptors independent of sialic acid residues. We hypothesized that over the years, intensified malaria control interventions and declining prevalence in The Gambia have resulted in a selection of parasites with a dominant invasion pathways and ligand expression profiles.

Carriers of persistent asymptomatic Plasmodium falciparum infections constitute an infectious reservoir that maintains malaria transmission. Understanding the extent of carriage and characteristics of carriers specific to endemic areas could guide use of interventions to reduce infectious reservoir.

Genetic diversity studies provide evidence of Plasmodium falciparum differentiation that could affect fitness and adaptation to drugs and target antigens for vaccine development. This study describes the genetic structure of P. falciparum populations in urban and rural sites from southwestern Nigeria.

Sequence diversity in pathogen antigens is an obstacle to the development of interventions against many infectious diseases. In malaria caused by Plasmodium falciparum, the PfEMP1 family of variant surface antigens encoded by var genes are adhesion molecules that play a pivotal role in malaria pathogenesis and clinical disease. PfEMP1 is a major target of protective immunity, however, development of drugs or vaccines based on PfEMP1 is problematic due to extensive sequence diversity within the PfEMP1 family. Here we identified the PfEMP1 variants transcribed by P. falciparum strains selected for a virulence-associated adhesion phenotype (IgM-positive rosetting). The parasites transcribed a subset of Group A PfEMP1 variants characterised by an unusual PfEMP1 architecture and a distinct N-terminal domain (either DBLα1.5 or DBLα1.8 type). Antibodies raised in rabbits against the N-terminal domains showed functional activity (surface reactivity with live infected erythrocytes (IEs), rosette inhibition and induction of phagocytosis of IEs) down to low concentrations (<10 µg/ml of total IgG) against homologous parasites. Furthermore, the antibodies showed broad cross-reactivity against heterologous parasite strains with the same rosetting phenotype, including clinical isolates from four sub-Saharan African countries that showed surface reactivity with either DBLα1.5 antibodies (variant HB3var6) or DBLα1.8 antibodies (variant TM284var1). These data show that parasites with a virulence-associated adhesion phenotype share IE surface epitopes that can be targeted by strain-transcending antibodies to PfEMP1. The existence of shared surface epitopes amongst functionally similar disease-associated P. falciparum parasite isolates suggests that development of therapeutic interventions to prevent severe malaria is a realistic goal.

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

Super-resistant Plasmodium falciparum threatens the effectiveness of sulfadoxine-pyrimethamine in intermittent preventive treatment for malaria during pregnancy. It is characterized by the A581G Pfdhps mutation on a background of the double-mutant Pfdhps and the triple-mutant Pfdhfr. Using samples collected during 2004-2008, we investigated the evolutionary origin of the A581G mutation by characterizing microsatellite diversity flanking Pfdhps triple-mutant (437G+540E+581G) alleles from 3 locations in eastern Africa and comparing it with double-mutant (437G+540E) alleles from the same area. In Ethiopia, both alleles derived from 1 lineage that was distinct from those in Uganda and Tanzania. Uganda and Tanzania triple mutants derived from the previously characterized southeastern Africa double-mutant lineage. The A581G mutation has occurred multiple times on local Pfdhps double-mutant backgrounds; however, a novel microsatellite allele incorporated into the Tanzania lineage since 2004 illustrates the local expansion of emergent triple-mutant lineages.

Mutations in the Plasmodium falciparum K13-propeller domain have recently been shown to be important determinants of artemisinin resistance in Southeast Asia. This study investigated the prevalence of K13-propeller polymorphisms across sub-Saharan Africa. A total of 1212 P. falciparum samples collected from 12 countries were sequenced. None of the K13-propeller mutations previously reported in Southeast Asia were found, but 22 unique mutations were detected, of which 7 were nonsynonymous. Allele frequencies ranged between 1% and 3%. Three mutations were observed in >1 country, and the A578S was present in parasites from 5 countries. This study provides the baseline prevalence of K13-propeller mutations in sub-Saharan Africa.

Vaccines are the most reliable alternative to elicit sterile immunity against malaria but their development has been hindered by polymorphisms and strain-specificity in previously studied antigens. New vaccine candidates are therefore urgently needed. Highly conserved Plasmodium falciparum reticulocyte-binding protein homologue-5 (PfRH5) has been identified as a potential candidate for anti-disease vaccine development. PfRH5 is essential for erythrocyte invasion by merozoites and crucial for parasite survival. However, there is paucity of data on the extent of genetic variations on PfRH5 in field isolates of Plasmodium falciparum. This study described genetic polymorphisms at the high affinity binding polypeptides (HABPs) 36718, 36727, 36728 of PfRH5 in Nigerian isolates of P. falciparum. This study tested the hypothesis that only specific conserved B and T cell epitopes on PfRH5 HABPs are crucial for vaccine development.

Background: Although thousands of clinical isolates of Plasmodium falciparum are being sequenced and analysed by short read technology, the data do not resolve the highly variable subtelomeric regions of the genomes that contain polymorphic gene families involved in immune evasion and pathogenesis. There is also no current standard definition of the boundaries of these variable subtelomeric regions. Methods: Using long-read sequence data (Pacific Biosciences SMRT technology), we assembled and annotated the genomes of 15 P. falciparum isolates, ten of which are newly cultured clinical isolates. We performed comparative analysis of the entire genome with particular emphasis on the subtelomeric regions and the internal var genes clusters.   Results: The nearly complete sequence of these 15 isolates has enabled us to define a highly conserved core genome, to delineate the boundaries of the subtelomeric regions, and to compare these across isolates. We found highly structured variable regions in the genome. Some exported gene families purportedly involved in release of merozoites show copy number variation. As an example of ongoing genome evolution, we found a novel CLAG gene in six isolates.  We also found a novel gene that was relatively enriched in the South East Asian isolates compared to those from Africa. Conclusions: These 15 manually curated new reference genome sequences with their nearly complete subtelomeric regions and fully assembled genes are an important new resource for the malaria research community. We report the overall conserved structure and pattern of important gene families and the more clearly defined subtelomeric regions.

Monitoring of Plasmodium falciparum sensitivity to antimalarial drugs in Africa is vital for malaria elimination. However, the commonly used ex vivo/in vitro 50% inhibitory concentration (IC50) test gives inconsistent results for several antimalarials, while the alternative ring-stage survival assay (RSA) for artemisinin derivatives has not been widely adopted. Here, we applied an alternative two-color flow cytometry-based parasite survival rate assay (PSRA) to detect ex vivo antimalarial tolerance in P. falciparum isolates from The Gambia. The PSRA infers parasite viability by quantifying reinvasion of uninfected cells following 3 consecutive days of drug exposure (10-fold the IC50 of drug for field isolates). The drug survival rate is obtained for each isolate from the slope of the growth/death curve. We obtained parasite survival rates of 41 isolates for dihydroartemisinin (DHA) and lumefantrine (LUM) out of 51 infections tested by ring-stage survival assay (RSA) against DHA. We also determined the genotypes for known drug resistance genetic loci in the P. falciparum genes Pfdhfr, Pfdhps, Pfmdr, Pfcrt, and Pfk13 The PSRA results determined for 41 Gambian isolates showed faster killing and lower variance after treatment with DHA than after treatment with LUM, despite a strong correlation between the two drugs. Four and three isolates were tolerant to DHA and LUM, respectively, with continuous growth during drug exposure. Isolates with the PfMDR1-Y184F mutant variant showed increased LUM survival, though the results were not statistically significant. Sulfadoxine/pyrimethamine (SP) resistance markers were fixed, while all other antimalarial variants were prevalent in more than 50% of the population. The PSRA detected ex vivo antimalarial tolerance in Gambian P. falciparum This calls for its wider application and for increased vigilance against resistance to artemisinin combination therapies (ACTs) in this population.

Antimalarial interventions have yielded a significant decline in malaria prevalence in The Gambia, where artemether-lumefantrine (AL) has been used as a first-line antimalarial for a decade. Clinical Plasmodium falciparum isolates collected from 2012 to 2015 were analyzed ex vivo for antimalarial susceptibility and genotyped for drug resistance markers (pfcrt K76T, pfmdr1 codons 86, 184, and 1246, and pfk13) and microsatellite variation. Additionally, allele frequencies of single nucleotide polymorphisms (SNPs) from other drug resistance-associated genes were compared from genomic sequence data sets from 2008 (n = 79) and 2014 (n = 168). No artemisinin resistance-associated pfk13 mutation was found, and only 4% of the isolates tested in 2015 showed significant growth after exposure to dihydroartemisinin. Conversely, the 50% inhibitory concentrations (IC50s) of amodiaquine and lumefantrine increased within this period. pfcrt 76T and pfmdr1 184F mutants remained at a prevalence above 80%. pfcrt 76T was positively associated with higher IC50s to chloroquine. pfmdr1 NYD increased in frequency between 2012 and 2015 due to lumefantrine selection. The TNYD (pfcrt 76T and pfmdr1 NYD wild-type haplotype) also increased in frequency following AL implementation in 2008. These results suggest selection for pfcrt and pfmdr1 genotypes that enable tolerance to lumefantrine. Increased tolerance to lumefantrine calls for sustained chemotherapeutic monitoring in The Gambia to minimize complete artemisinin combination therapy (ACT) failure in the future.

Development of resistance to deployed antimalarial drugs is inevitable and needs prompt and continuous discovery of novel candidate drugs. Therefore, the antimalarial activity of 125 compounds from the Medicine for Malaria Ventures (MMV) pathogen box was determined. Combining standard IC50 and normalised growth rate inhibition (GR50) analyses, we found 16 and 22 compounds had higher potencies than CQ respectively. Seven compounds with relatively high potencies (low GR50 and IC50) against P. falciparum 3D7 were further analysed. Three of these were tested on 10 natural P. falciparum isolates from The Gambia using our newly developed parasite survival rate assay (PSRA). According to the IC50, GR50 and PSRA analyses, compound MMV667494 was most potent and highly cytotoxic to parasites. MMV010576 was slow acting but more potent than dihydroartemisinin (DHA) 72 h after exposure. MMV634140 was potent against the laboratory-adapted 3D7 isolate, but 4 out of 10 natural Gambian isolates survived and replicated slowly despite 72 h of exposure to the compound, suggesting potential drug tolerance and risk of resistance development. These results emphasise the usefulness of in vitro testing as a starting point for drug discovery. Improved approaches to data analyses and the use of natural isolates will facilitate the prioritisation of compounds for further clinical development.

Inferences from genetic association studies rely largely on the definition and description of the underlying populations that highlight their genetic similarities and differences. The clustering of human populations into subgroups (population structure) can significantly confound disease associations. This study investigated the fine-scale genetic structure within Cameroon that may underlie disparities observed with Cameroonian ethnicities in malaria genome-wide association studies in sub-Saharan Africa. Genotype data of 1073 individuals from three regions and three ethnic groups in Cameroon were analyzed using measures of genetic proximity to ascertain fine-scale genetic structure. Model-based clustering revealed distinct ancestral proportions among the Bantu, Semi-Bantu and Foulbe ethnic groups, while haplotype-based coancestry estimation revealed possible longstanding and ongoing sympatric differentiation among individuals of the Foulbe ethnic group, and their Bantu and Semi-Bantu counterparts. A genome scan found strong selection signatures in the HLA gene region, confirming longstanding knowledge of natural selection on this genomic region in African populations following immense disease pressure. Signatures of selection were also observed in the HBB gene cluster, a genomic region known to be under strong balancing selection in sub-Saharan Africa due to its co-evolution with malaria. This study further supports the role of evolution in shaping genomes of Cameroonian populations and reveals fine-scale hierarchical structure among and within Cameroonian ethnicities that may impact genetic association studies in the country.