The power and promise of genetic mapping from Plasmodium falciparum crosses utilizing human liver-chimeric mice

MalKinID: A classification model for identifying malaria parasite genealogical relationships using identity-by-descent

Pathogen genomics is a powerful tool for tracking infectious disease transmission. In malaria, identity-by-descent is used to assess the genetic relatedness between parasites and has been used to study transmission and importation. In theory, identity-by-descent can be used to distinguish genealogical relationships to reconstruct transmission history or identify parasites for QTL experiments. MalKinID (Malaria Kinship Identifier) is a new classification model designed to identify genealogical relationships among malaria parasites based on genome-wide identity-by-descent proportions and identity-by-descent segment distributions. MalKinID was calibrated to the genomic data from 3 laboratory-based genetic crosses (yielding 440 parent-child and 9060 full-sibling comparisons). MalKinID identified lab-generated F1 progeny with >80% sensitivity and showed that 0.39 (95% CI 0.28, 0.49) of the second-generation progeny of a NF54 and NHP4026 cross were F1s and 0.56 (0.45, 0.67) were backcrosses of an F1 with the parental NF54 strain. In simulated outcrossed importations, MalKinID reconstructs genealogy history with high precision and sensitivity, with F1-scores exceeding 0.84. However, when importation involves inbreeding, such as during serial co-transmission, the precision and sensitivity of MalKinID declined, with F1-scores (the harmonic mean of precision and sensitivity) of 0.76 (0.56, 0.92) and 0.23 (0.0, 0.4) for parent-child and full-sibling and <0.05 for second-degree and third-degree relatives. Disentangling inbred relationships required adapting MalKinID to perform multisample comparisons. Genealogical inference is most powered when (1) outcrossing is the norm or (2) multisample comparisons based on a predefined pedigree are used. MalKinID lays the foundations for using identity-by-descent to track parasite transmission history and for separating progeny for quantitative-trait-locus experiments.

Systematic in vitro evolution in Plasmodium falciparum reveals key determinants of drug resistance

Surveillance of drug resistance and the discovery of novel targets-key objectives in the fight against malaria-rely on identifying resistance-conferring mutations in Plasmodium parasites. Current approaches, while successful, require laborious experimentation or large sample sizes. To elucidate shared determinants of antimalarial resistance that can empower in silico inference, we examined the genomes of 724 Plasmodium falciparum clones, each selected in vitro for resistance to one of 118 compounds. We identified 1448 variants in 128 recurrently mutated genes, including drivers of antimalarial multidrug resistance. In contrast to naturally occurring variants, those selected in vitro are more likely to be missense or frameshift, involve bulky substitutions, and occur in conserved, ordered protein domains. Collectively, our dataset reveals mutation features that predict drug resistance in eukaryotic pathogens.

Identification of the drug/metabolite transporter 1 as a marker of quinine resistance in a NF54×Cam3.II P. falciparum genetic cross

The genetic basis of Plasmodium falciparum resistance to quinine (QN), a drug used to treat severe malaria, has long been enigmatic. To gain further insight, we used FRG-NOD human liver-chimeric mice to conduct a P. falciparum genetic cross between QN-sensitive and QN-resistant parasites, which also differ in their susceptibility to chloroquine (CQ). By applying different selective conditions to progeny pools prior to cloning, we recovered 120 unique recombinant progeny. These progeny were subjected to drug profiling and QTL analyses with QN, CQ, and monodesethyl-CQ (md-CQ, the active metabolite of CQ), which revealed predominant peaks on chromosomes 7 and 12, consistent with a multifactorial mechanism of resistance. A shared chromosome 12 region mapped to resistance to all three antimalarials and was preferentially co-inherited with pfcrt. We identified an ATP-dependent zinc metalloprotease (FtsH1) as one of the top candidates and observed using CRISPR/Cas9 SNP-edited lines that ftsh1 is a potential mediator of QN resistance and a modulator of md-CQ resistance. As expected, CQ and md-CQ resistance mapped to a chromosome 7 region harboring pfcrt. However, for QN, high-grade resistance mapped to a chromosome 7 peak centered 295kb downstream of pfcrt. We identified the drug/metabolite transporter 1 (DMT1) as the top candidate due to its structural similarity to PfCRT and proximity to the peak. Deleting DMT1 in QN-resistant Cam3.II parasites significantly sensitized the parasite to QN but not to the other drugs tested, suggesting that DMT1 mediates QN response specifically. We localized DMT1 to structures associated with vesicular trafficking, as well as the parasitophorous vacuolar membrane, lipid bodies, and the digestive vacuole. We also observed that mutant DMT1 transports more QN than the wild-type isoform in vitro. Our study demonstrates that DMT1 is a novel marker of QN resistance and a new chromosome 12 locus associates with CQ and QN response, with ftsh1 is a potential candidate, suggesting these genes should be genotyped in surveillance and clinical settings.

A Plasmodium falciparum genetic cross reveals the contributions of pfcrt and plasmepsin II/III to piperaquine drug resistance

Piperaquine (PPQ) is widely used in combination with dihydroartemisinin as a first-line treatment against malaria. Multiple genetic drivers of PPQ resistance have been reported, including mutations in the Plasmodium falciparum chloroquine resistance transporter (pfcrt) and increased copies of plasmepsin II/III (pm2/3). We generated a cross between a Cambodia-derived multidrug-resistant KEL1/PLA1 lineage isolate (KH004) and a drug-susceptible Malawian parasite (Mal31). Mal31 harbors a wild-type (3D7-like) pfcrt allele and a single copy of pm2/3, while KH004 has a chloroquine-resistant (Dd2-like) pfcrt allele with an additional G367C substitution and multiple copies of pm2/3. We recovered 104 unique recombinant parasites and examined a targeted set of progeny representing all possible combinations of variants at pfcrt and pm2/3. We performed a detailed analysis of competitive fitness and a range of PPQ susceptibility phenotypes with these progenies, including PPQ survival assay, area under the dose response curve, and a limited point IC50. We find that inheritance of the KH004 pfcrt allele is required for reduced PPQ sensitivity, whereas copy number variation in pm2/3 further decreases susceptibility but does not confer resistance in the absence of additional mutations in pfcrt. A deep investigation of genotype-phenotype relationships demonstrates that progeny clones from experimental crosses can be used to understand the relative contributions of pfcrt, pm2/3, and parasite genetic background to a range of PPQ-related traits. Additionally, we find that the resistance phenotype associated with parasites inheriting the G367C substitution in pfcrt is consistent with previously validated PPQ resistance mutations in this transporter.IMPORTANCEResistance to piperaquine, used in combination with dihydroartemisinin, has emerged in Cambodia and threatens to spread to other malaria-endemic regions. Understanding the causal mutations of drug resistance and their impact on parasite fitness is critical for surveillance and intervention and can also reveal new avenues to limiting the evolution and spread of drug resistance. An experimental genetic cross is a powerful tool for pinpointing the genetic determinants of key drug resistance and fitness phenotypes and has the distinct advantage of quantifying the effects of naturally evolved genetic variation. Our study was strengthened since the full range of copies of KH004 pm2/3 was inherited among the progeny clones, allowing us to directly test the role of the pm2/3 copy number on resistance-related phenotypes in the context of a unique pfcrt allele. Our multigene model suggests an important role for both loci in the evolution of this multidrug-resistant parasite lineage.

MalKinID: A Likelihood-Based Model for Identifying Malaria Parasite Genealogical Relationships Using Identity-by-Descent

Pathogen genomics is a powerful tool for tracking infectious disease transmission. In malaria, identity-by-descent (IBD) is used to assess the genetic relatedness between parasites and has been used to study transmission and importation. In theory, IBD can be used to distinguish genealogical relationships to reconstruct transmission history or identify parasites for genotype-to-phenotype quantitative-trait-locus experiments. MalKinID (Malaria Kinship Identifier) is a new likelihood-based classification model designed to identify genealogical relationships among malaria parasites based on genome-wide IBD proportions and IBD segment distributions. MalKinID was calibrated to the genomic data from three laboratory-based genetic crosses (yielding 440 parent-child and 9060 full-sibling comparisons). MalKinID identified lab generated F1 progeny with >80% sensitivity and showed that 0.39 (95% CI 0.28, 0.49) of the second-generation progeny of a NF54 and NHP4026 cross were F1s and 0.56 (0.45, 0.67) were backcrosses of an F1 with the parental NF54 strain. In simulated outcrossed importations, MalKinID accurately reconstructs genealogy history with high precision and sensitivity, with F1-scores exceeding 0.84. However, when importation involves inbreeding, such as during serial co-transmission, the precision and sensitivity of MalKinID declined, with F1-scores of 0.76 (0.56, 0.92) and 0.23 (0.0, 0.4) for PC and FS and <0.05 for second-degree and third-degree relatives. Genealogical inference is most powered 1) when outcrossing is the norm or 2) when multi-sample comparisons based on a predefined pedigree are used. MalKinID lays the foundations for using IBD to track parasite transmission history and for separating progeny for quantitative-trait-locus experiments.

Mapping the genomic landscape of multidrug resistance in Plasmodium falciparum and its impact on parasite fitness

Drug-resistant Plasmodium falciparum parasites have swept across Southeast Asia and now threaten Africa. By implementing a P. falciparum genetic cross using humanized mice, we report the identification of key determinants of resistance to artemisinin (ART) and piperaquine (PPQ) in the dominant Asian KEL1/PLA1 lineage. We mapped k13 as the central mediator of ART resistance in vitro and identified secondary markers. Applying bulk segregant analysis, quantitative trait loci mapping using 34 recombinant haplotypes, and gene editing, our data reveal an epistatic interaction between mutant PfCRT and multicopy plasmepsins 2/3 in mediating high-grade PPQ resistance. Susceptibility and parasite fitness assays implicate PPQ as a driver of selection for KEL1/PLA1 parasites. Mutant PfCRT enhanced susceptibility to lumefantrine, the first-line partner drug in Africa, highlighting a potential benefit of opposing selective pressures with this drug and PPQ. We also identified that the ABCI3 transporter can operate in concert with PfCRT and plasmepsins 2/3 in mediating multigenic resistance to antimalarial agents.

A Plasmodium falciparum genetic cross reveals the contributions of pfcrt and plasmepsin II/III to piperaquine drug resistance

Piperaquine (PPQ) is widely used in combination with dihydroartemisinin (DHA) as a first-line treatment against malaria parasites. Multiple genetic drivers of PPQ resistance have been reported, including mutations in the Plasmodium falciparum chloroquine resistance transporter (pfcrt) and increased copies of plasmepsin II/III (pm2/3). We generated a cross between a Cambodia-derived multi-drug resistant KEL1/PLA1 lineage isolate (KH004) and a drug susceptible parasite isolated in Malawi (Mal31). Mal31 harbors a wild-type (3D7-like) pfcrt allele and a single copy of pm2/3, while KH004 has a chloroquine-resistant (Dd2-like) pfcrt allele with an additional G367C substitution and four copies of pm2/3. We recovered 104 unique recombinant progeny and examined a targeted set of progeny representing all possible combinations of variants at pfcrt and pm2/3 for detailed analysis of competitive fitness and a range of PPQ susceptibility phenotypes, including PPQ survival assay (PSA), area under the dose-response curve (AUC), and a limited point IC50 (LP-IC50). We find that inheritance of the KH004 pfcrt allele is required for PPQ resistance, whereas copy number variation in pm2/3 further enhances resistance but does not confer resistance in the absence of PPQ-R-associated mutations in pfcrt. Deeper investigation of genotype-phenotype relationships demonstrates that progeny clones from experimental crosses can be used to understand the relative contributions of pfcrt, pm2/3, and parasite genetic background, to a range of PPQ-related traits and confirm the critical role of the PfCRT G367C substitution in PPQ resistance.

Chloroquine resistance evolution in Plasmodium falciparum is mediated by the putative amino acid transporter AAT1

Malaria parasites break down host haemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter, pfcrt, located in the digestive vacuole membrane, confer CQ resistance in Plasmodium falciparum, and typically also affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on a putative amino acid transporter (pfaat1) variant S258L, which increased from 0% to 97% in frequency between 1984 and 2014 in parallel with the pfcrt1 K76T variant. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containing pfaat1 that is selected by CQ treatment. Gene editing demonstrated that pfaat1 S258L potentiates CQ resistance but at a cost of reduced fitness, while pfaat1 F313S, a common southeast Asian polymorphism, reduces CQ resistance while restoring fitness. Our analyses reveal hidden complexity in CQ resistance evolution, suggesting that pfaat1 may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport.

Mapping the genomic landscape of multidrug resistance in Plasmodium falciparum and its impact on parasite fitness

Drug-resistant Plasmodium falciparum parasites have swept across Southeast Asia and now threaten Africa. By implementing a P. falciparum genetic cross using humanized mice, we report the identification of key determinants of resistance to artemisinin (ART) and piperaquine (PPQ) in the dominant Asian KEL1/PLA1 lineage. We mapped k13 as the central mediator of ART resistance and identified secondary markers. Applying bulk segregant analysis, quantitative trait loci mapping and gene editing, our data reveal an epistatic interaction between mutant PfCRT and multicopy plasmepsins 2/3 in mediating high-grade PPQ resistance. Susceptibility and parasite fitness assays implicate PPQ as a driver of selection for KEL1/PLA1 parasites. Mutant PfCRT enhanced susceptibility to lumefantrine, the first-line partner drug in Africa, highlighting a potential benefit of opposing selective pressures with this drug and PPQ. We also identified that the ABCI3 transporter can operate in concert with PfCRT and plasmepsins 2/3 in mediating multigenic resistance to antimalarial agents.

Bulk segregant linkage mapping for rodent and human malaria parasites

In 2005 Richard Carter's group surprised the malaria genetics community with an elegant approach to rapidly mapping the genetic basis of phenotypic traits in rodent malaria parasites. This approach, which he termed "linkage group selection", utilized bulk pools of progeny, rather than individual clones, and exploited simple selection schemes to identify genome regions underlying resistance to drug treatment (or other phenotypes). This work was the first application of "bulk segregant" methodologies for genetic mapping in microbes: this approach is now widely used in yeast, and across multiple recombining pathogens ranging from Aspergillus fungi to Schistosome parasites. Genetic crosses of human malaria parasites (for which Richard Carter was also a pioneer) can now be conducted in humanized mice, providing new opportunities for exploiting bulk segregant approaches for a wide variety of malaria parasite traits. We review the application of bulk segregant approaches to mapping malaria parasite traits and suggest additional developments that may further expand the utility of this powerful approach.

Plasmodium falciparum resistance to artemisinin-based combination therapies

Multidrug-resistant Plasmodium falciparum parasites are a major threat to public health in intertropical regions. Understanding the mechanistic basis, origins, and spread of resistance can inform strategies to mitigate its impact and reduce the global burden of malaria. The recent emergence in Africa of partial resistance to artemisinins, the core component of first-line combination therapies, is particularly concerning. Here, we review recent advances in elucidating the mechanistic basis of artemisinin resistance, driven primarily by point mutations in P. falciparum Kelch13, a key regulator of hemoglobin endocytosis and parasite response to artemisinin-induced stress. We also review resistance to partner drugs, including piperaquine and mefloquine, highlighting a key role for plasmepsins 2/3 and the drug and solute transporters P. falciparum chloroquine-resistance transporter and P. falciparum multidrug-resistance protein-1.

Keystone Malaria Symposium 2022: a vibrant discussion of progress made and challenges ahead from drug discovery to treatment

In recent years, the field of malaria research has made substantial progress in the areas of antimalarial drug resistance and discovery. These efforts are essential to combatting the devastating impact of malaria, which, in 2020, resulted in an estimated 241 million cases and 627 000 deaths. Recent advances in this area were presented at a Keystone Symposium entitled ‘Malaria: Confronting Challenges from Drug Discovery to Treatment’, held in person in Breckenridge, Colorado, in April 2022. Herein, we present a summary of the proceedings of this vibrant scientific exchange, which brought together a superb group of faculty, postdocs, and students from around the globe.

A Malaria Parasite Cross Reveals Genetic Determinants of Plasmodium falciparum Growth in Different Culture Media

What genes determine in vitro growth and nutrient utilization in asexual blood-stage malaria parasites? Competition experiments between NF54, clone 3D7, a lab-adapted African parasite, and a recently isolated Asian parasite (NHP4026) reveal contrasting outcomes in different media: 3D7 outcompetes NHP4026 in media containing human serum, while NHP4026 outcompetes 3D7 in media containing AlbuMAX, a commercial lipid-rich bovine serum formulation. To determine the basis for this polymorphism, we conducted parasite genetic crosses using humanized mice and compared genome-wide allele frequency changes in three independent progeny populations cultured in media containing human serum or AlbuMAX. This bulk segregant analysis detected three quantitative trait loci (QTL) regions [on chromosome (chr) 2 containing aspartate transaminase AST; chr 13 containing EBA-140; and chr 14 containing cysteine protease ATG4] linked with differential growth in serum or AlbuMAX in each of the three independent progeny pools. Selection driving differential growth was strong (s = 0.10 – 0.23 per 48-hour lifecycle). We conducted validation experiments for the strongest QTL on chr 13: competition experiments between ΔEBA-140 and 3D7 wildtype parasites showed fitness reversals in the two medium types as seen in the parental parasites, validating this locus as the causative gene. These results (i) demonstrate the effectiveness of bulk segregant analysis for dissecting fitness traits in P. falciparum genetic crosses, and (ii) reveal intimate links between red blood cell invasion and nutrient composition of growth media. Use of parasite crosses combined with bulk segregant analysis will allow systematic dissection of key nutrient acquisition/metabolism and red blood cell invasion pathways in P. falciparum.

Optimizing bulk segregant analysis of drug resistance using Plasmodium falciparum genetic crosses conducted in humanized mice

Classical malaria parasite genetic crosses involve isolation, genotyping, and phenotyping of progeny parasites, which is time consuming and laborious. We tested a rapid alternative approach-bulk segregant analysis (BSA)-that utilizes sequencing of bulk progeny populations with and without drug selection for rapid identification of drug resistance loci. We used dihydroartemisinin (DHA) selection in two genetic crosses and investigated how synchronization, cryopreservation, and the drug selection regimen impacted BSA success. We detected a robust quantitative trait locus (QTL) at kelch13 in both crosses but did not detect QTLs at four other candidate loci. QTLs were detected using synchronized, but not unsynchronized progeny pools, consistent with the stage-specific action of DHA. We also successfully applied BSA to cryopreserved progeny pools, expanding the utility of this approach. We conclude that BSA provides a powerful approach for investigating the genetic architecture of drug resistance in Plasmodium falciparum.