Quantification of multiple infections of Plasmodium falciparum in vitro

Cultivation of Asexual Intraerythrocytic Stages of Plasmodium falciparum

Successfully developed in 1976, the continuous in vitro culture of Plasmodium falciparum has many applications in the field of malaria research. It has become an important experimental model that directly uses a human pathogen responsible for a high prevalence of morbidity and mortality in many parts of the world and is a major source of biological material for immunological, biochemical, molecular, and pharmacological studies. Until present, the basic techniques described by Trager and Jensen and Haynes et al. remain unchanged in many malaria research laboratories. Nonetheless, different factors, including culture media, buffers, serum substitutes and supplements, sources of erythrocytes, and conditions of incubation (especially oxygen concentration), have been modified by different investigators to adapt the original technique in their laboratories or enhance the in vitro growth of the parasites. The possible effects and benefits of these modifications for the continuous cultivation of asexual intraerythrocytic stages of P. falciparum, as well as future challenges in developing a serum-free cultivation system and axenic cultures, are discussed.

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.

Deconstructing the parasite multiplication rate of Plasmodium falciparum

Epidemiological indicators describing population-level malaria transmission dynamics are widely used to guide policy recommendations. However, the determinants of malaria outcomes within individuals are still poorly understood. This conceptual gap partly reflects the fact that there are few indicators that robustly predict the trajectory of individual infections or clinical outcomes. The parasite multiplication rate (PMR) is a widely used indicator for the Plasmodium intraerythrocytic development cycle (IDC), for example, but its relationship to clinical outcomes is complex. Here, we review its calculation and use in P. falciparum malaria research, as well as the parasite and host factors that impact it. We also provide examples of metrics that can help to link within-host dynamics to malaria clinical outcomes when used alongside the PMR.

Pairwise growth competitions identify relative fitness relationships among artemisinin resistant Plasmodium falciparum field isolates

Background

Competitive outcomes between co-infecting malaria parasite lines can reveal fitness disparities in blood stage growth. Blood stage fitness costs often accompany the evolution of drug resistance, with the expectation that relatively fitter parasites will be more likely to spread in populations. With the recent emergence of artemisinin resistance, it is important to understand the relative competitive fitness of the metabolically active asexual blood stage parasites. Genetically distinct drug resistant parasite clones with independently evolved sets of mutations are likely to vary in asexual proliferation rate, contributing to their chance of transmission to the mosquito vector.

Methods

An optimized in vitro 96-well plate-based protocol was used to quantitatively measure-head-to-head competitive fitness during blood stage development between seven genetically distinct field isolates from a hotspot of emerging artemisinin resistance and the laboratory strain, NF54. These field isolates were isolated from patients in Southeast Asia carrying different alleles of kelch13 and included both artemisinin-sensitive and artemisinin-resistant isolates. Fluorescent labeled microsatellite markers were used to track the relative densities of each parasite throughout the co-growth period of 14–60 days. All-on-all competitions were conducted for the panel of eight parasite lines (28 pairwise competitions) to determine their quantitative competitive fitness relationships.

Results

Twenty-eight pairwise competitive growth outcomes allowed for an unambiguous ranking among a set of seven genetically distinct parasite lines isolated from patients in Southeast Asia displaying a range of both kelch13 alleles and clinical clearance times and a laboratory strain, NF54. This comprehensive series of assays established the growth relationships among the eight parasite lines. Interestingly, a clinically artemisinin resistant parasite line that carries the wild-type form of kelch13 outcompeted all other parasites in this study. Furthermore, a kelch13 mutant line (E252Q) was competitively more fit without drug than lines with other resistance-associated kelch13 alleles, including the C580Y allele that has expanded to high frequencies under drug pressure in Southeast Asian resistant populations.

Conclusions

This optimized competitive growth assay can be employed for assessment of relative growth as an index of fitness during the asexual blood stage growth between natural lines carrying different genetic variants associated with artemisinin resistance. Improved understanding of the fitness costs of different parasites proliferating in human blood and the role different resistance mutations play in the context of specific genetic backgrounds will contribute to an understanding of the potential for specific mutations to spread in populations, with the potential to inform targeted strategies for malaria therapy.

Large and finite sample properties of a maximum-likelihood estimator for multiplicity of infection

Reliable measures of transmission intensities can be incorporated into metrics for monitoring disease-control interventions. Genetic (molecular) measures like multiplicity of infection (MOI) have several advantages compared with traditional measures, e.g., R0. Here, we investigate the properties of a maximum-likelihood approach to estimate MOI and pathogen-lineage frequencies. By verifying regulatory conditions, we prove asymptotical unbiasedness, consistency and efficiency of the estimator. Finite sample properties concerning bias and variance are evaluated over a comprehensive parameter range by a systematic simulation study. Moreover, the estimator's sensitivity to model violations is studied. The estimator performs well for realistic sample sizes and parameter ranges. In particular, the lineage-frequency estimates are almost unbiased independently of sample size. The MOI estimate's bias vanishes with increasing sample size, but might be substantial if sample size is too small. The estimator's variance matrix agrees well with the Cramér-Rao lower bound, even for small sample size. The numerical and analytical results of this study can be used for study design. This is exemplified by a malaria data set from Venezuela. It is shown how the results can be used to determine the necessary sample size to achieve certain performance goals. An implementation of the likelihood method and a simulation algorithm for study design, implemented as an R script, is available as S1 File alongside a documentation (S2 File) and example data (S3 File).

Multiplication rate variation in the human malaria parasite Plasmodium falciparum

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.

Efficacy of OZ439 (artefenomel) against early Plasmodium falciparum blood-stage malaria infection in healthy volunteers

Objectives

OZ439, or artefenomel, is an investigational synthetic ozonide antimalarial with similar potency, but a significantly improved pharmacokinetic profile, compared with artemisinins. We wished to measure key pharmacokinetic and pharmacodynamic parameters and the pharmacokinetic/pharmacodynamic relationship of artefenomel in humans to guide the drug's further development as combination therapy in patients.

Patients and methods

We tested artefenomel in the human induced blood-stage malaria (IBSM) model. Plasmodium infection was monitored by quantitative PCR (qPCR) and upon reaching 1000 parasites/mL single doses of 100, 200 and 500 mg of artefenomel were administered orally with evaluation of drug exposure and parasitaemia until rescue treatment after 16 days or earlier, if required.

Results

A single 100 mg dose had only a transient effect, while the 200 mg dose resulted in a significant reduction in parasitaemia before early recrudescence. At the highest (500 mg) dose, initial clearance of parasites below the limit of detection of qPCR was observed, with a 48 h parasite reduction ratio (PRR₄₈) >10 000 and a parasite clearance half-life of 3.6 h (95% CI 3.4–3.8 h). However, at this dose, recrudescence was seen in four of eight subjects 6–10 days after treatment. Pharmacokinetic/pharmacodynamic modelling predicted an MIC of 4.1 ng/mL.

Conclusions

These results confirm the antimalarial potential of artefenomel for use in a single-exposure combination therapy. The observations from this study support and will assist further clinical development of artefenomel.

The potential impact of coinfection on antimicrobial chemotherapy and drug resistance

Across a range of pathogens, resistance to chemotherapy is a growing problem in both public health and animal health. Despite the ubiquity of coinfection, and its potential effects on within-host biology, the role played by coinfecting pathogens on the evolution of resistance and efficacy of antimicrobial chemotherapy is rarely considered. In this review, we provide an overview of the mechanisms of interaction of coinfecting pathogens, ranging from immune modulation and resource modulation, to drug interactions. We discuss their potential implications for the evolution of resistance, providing evidence in the rare cases where it is available. Overall, our review indicates that the impact of coinfection has the potential to be considerable, suggesting that this should be taken into account when designing antimicrobial drug treatments.

A likelihood approach to estimate the number of co-infections

The number of co-infections of a pathogen (multiplicity of infection or MOI) is a relevant parameter in epidemiology as it relates to transmission intensity. Notably, such quantities can be built into a metric in the context of disease control and prevention. Having applications to malaria in mind, we develop here a maximum-likelihood (ML) framework to estimate the quantities of interest at low computational and no additional costs to study designs or data collection. We show how the ML estimate for the quantities of interest and corresponding confidence-regions are obtained from multiple genetic loci. Assuming specifically that infections are rare and independent events, the number of infections per host follows a conditional Poisson distribution. Under this assumption, we show that a unique ML estimate for the parameter () describing MOI exists which is found by a simple recursion. Moreover, we provide explicit formulas for asymptotic confidence intervals, and show that profile-likelihood-based confidence intervals exist, which are found by a simple two-dimensional recursion. Based on the confidence intervals we provide alternative statistical tests for the MOI parameter. Finally, we illustrate the methods on three malaria data sets. The statistical framework however is not limited to malaria.

Multi-strain infections and 'relapse' of Leucocytozoon sabrazesi gametocytes in domestic chickens in southern China

Leucocytozoon parasites infect many species of avian hosts, including domestic chicken, and can inflict heavy economic loss to the poultry industry. Although the prevalence and distribution of two Leucocytozoon species (L. sabrazesi and L. caulleryi) have been reported in China previously, there are many questions related to the parasite infection that remain unanswered, including population diversity and transmission dynamics in domestic chickens. Here we surveyed chicken blood samples from seven sites in four provinces of China to identify Leucocytozoon infection, characterized parasite diversity within individual infected hosts and between sampling sites, and investigated the dynamics of gametocytemia in chickens over time. We found high infection rates in three of the seven sites. Clustering parasite sequences of the mitochondrial cytochrome oxidase III (coxIII) and cytochrome b (cytb) genes showed lack of grouping according to geographic origins and individual hosts carrying large numbers of L. sabrazesi strains. Monitoring gametocytemia in blood samples from infected chickens over time showed ‘relapse’ or persistence of low-level gametocytemia for 4–5 months, which could be explored as an in vivo model for testing drugs against liver stages of Apicomplexan parasites. This study provides important information on population diversity and transmission dynamics of L. sabrazesi and for disease control.

The role of asymptomatic P. falciparum parasitaemia in the evolution of antimalarial drug resistance in areas of seasonal transmission

In areas with seasonal transmission, proper management of acute malaria cases that arise in the transmission season can markedly reduce the disease burden. However, asymptomatic carriage of Plasmodium falciparum sustains a long-lasting reservoir in the transmission-free dry season that seeds cyclical malaria outbreaks. Clinical trials targeting asymptomatic parasitaemia in the dry season failed to interrupt the malaria epidemics that follow annual rains. These asymptomatic infections tend to carry multiple-clones, capable of producing gametocytes and infecting Anopheles mosquitoes. Different clones within an infection fluctuate consistently, indicative of interaction between clones during the long course of asymptomatic carriage. However, the therapy-free environment that prevails in the dry season dis-advantages the drug resistant lineages and favors the wild-type parasites. This review highlights some biological and epidemiological characteristics of asymptomatic parasitaemia and calls for consideration of policies to diminish parasite exposure to drugs "therapy-free" and allow natural selection to curb drug resistance in the above setting.