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

The F204S mutation in adrenodoxin oxidoreductase drives salinomycin resistance in Eimeria tenella

Salinomycin is a polyether ionophore widely used for the treatment of coccidiosis in poultry. However, the emergence of coccidia strains resistant to salinomycin presents challenges for control efforts, and the mechanisms underlying this resistance in Eimeria remain inadequately understood. In this study, 78 stable salinomycin-resistant strains were generated through experimental evolution approaches. Whole-genome sequencing of salinomycin-resistant Eimeria tenella isolates revealed single nucleotide polymorphisms (SNPs), with 12 candidate genes harboring nonsynonymous mutations identified. To confirm the candidate gene responsible for conferring salinomycin resistance, we leveraged reverse genetic strategies and identified a key amino acid substitution (F204S) in adrenodoxin oxidoreductase (EtADR), which markedly reduced susceptibility to salinomycin. Our results elucidate the complex interactions among salinomycin resistance, parasite fitness, point mutations, and the structure of EtADR, laying the foundation for future studies on drug resistance in Eimeria and contributing to the development of targeted control strategies.

Measuring Growth, Resistance, and Recovery after Artemisinin Treatment of Plasmodium falciparum in a single semi-high-throughput Assay

Background

Artemisinin partial resistance (ART-R) has spread throughout Southeast Asia and mutations in pfKelch13, the molecular marker of resistance, are widely reported in East Africa. Effective in vitro assays and robust phenotypes are crucial for monitoring populations for the emergence and spread of resistance. The recently developed extended Recovery Ring-stage Survival Assay used a qPCR-based readout to reduce the labor intensiveness for in vitro phenotyping of ART-R and improved correlation with the clinical phenotype of ART-R. Here, we extend and refine this assay to include measurements of parasite growth and recovery after drug exposure. Clinical isolates and progeny from two genetic crosses were used to optimize and validate the reliability of a straight-from-blood, SYBR Green-based qPCR protocol in a 96-well plate format to accurately measure phenotypes for Growth, Resistance, and Recovery.

Results

The assay determined growth between 6 h and 96 h, resistance at 120 h, and recovery from 120 h and 192 h. Growth can be accurately captured by qPCR and is shown by reproduction of previous growth phenotypes from HB3 × Dd2. Resistance measured at 120 h continually shows the most consistent phenotype for ring stage susceptibility. Recovery identifies an additional response to drug than parasites that are determined sensitive by Fold Change at 120 h. Comparison of progeny phenotypes for Growth vs Resistance showed a minor but significant correlation, whereas Growth vs Recovery and Resistance vs Recovery showed no significant correlation. Additionally, dried blood spot (DBS) samples matched Fold Change measured from liquid samples demonstrating Resistance can be easily quantified using either storage method.

Conclusions

The qPCR-based methodology provides the throughput needed to quickly measure large numbers of parasites for multiple relevant phenotypes. Growth can reveal fitness defects and illuminate relationships between proliferation rates and drug response. Recovery serves as a complementary phenotype to resistance that quantifies the ability of sensitive parasites to tolerate drug exposure. All three phenotypes offer a comprehensive assessment of parasite-drug interaction each with independent genetic determinants of main effect and overlapping secondary effects that should be further. By adapting our method to include DBS, readouts can be easily extended to ex vivo surveillance applications.

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.

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.

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.

Barcoding Genetically Distinct Plasmodium falciparum Strains for Comparative Assessment of Fitness and Antimalarial Drug Resistance

The repeated emergence of antimalarial drug resistance in Plasmodium falciparum, including to the current frontline antimalarial artemisinin, is a perennial problem for malaria control. Next-generation sequencing has greatly accelerated the identification of polymorphisms in resistance-associated genes but has also highlighted the need for more sensitive and accurate laboratory tools to profile current and future antimalarials and to quantify the impact of drug resistance acquisition on parasite fitness. The interplay of fitness and drug response is of fundamental importance in understanding why particular genetic backgrounds are better at driving the evolution of drug resistance in natural populations, but the impact of parasite fitness landscapes on the epidemiology of drug resistance has typically been laborious to accurately quantify in the lab, with assays being limited in accuracy and throughput. Here we present a scalable method to profile fitness and drug response of genetically distinct P. falciparum strains with well-described sensitivities to several antimalarials. We leverage CRISPR/Cas9 genome-editing and barcode sequencing to track unique barcodes integrated into a nonessential gene (pfrh3). We validate this approach in multiplex competitive growth assays of three strains with distinct geographical origins. Furthermore, we demonstrate that this method can be a powerful approach for tracking artemisinin response as it can identify an artemisinin resistant strain within a mix of multiple parasite lines, suggesting an approach for scaling the laborious ring-stage survival assay across libraries of barcoded parasite lines. Overall, we present a novel high-throughput method for multiplexed competitive growth assays to evaluate parasite fitness and drug response. IMPORTANCE The complex interplay between antimalarial resistance and parasite fitness has important implications for understanding the development and spread of drug resistance alleles and the impact of genetic background on transmission. One limitation with current methodologies to measure parasite fitness is the ability to scale this beyond simple head-to-head competition experiments between a wildtype control line and test line, with a need for a scalable approach that allows tracking of parasite growth in complex mixtures. In our study, we have used CRISPR editing to insert unique DNA barcodes into a safe-harbor genomic locus to tag multiple parasite strains and use next-generation sequencing to read out strain dynamics. We observe inherent fitness differences between the strains, as well as sensitive modulation of responses to challenge with clinically relevant antimalarials, including artemisinin.