Transcriptional heterogeneity and tightly regulated changes in gene expression during Plasmodium berghei sporozoite development.

Despite the critical role of Plasmodium sporozoites in malaria transmission, we still know little about the mechanisms underlying their development in mosquitoes. Here, we use single-cell RNA sequencing to characterize the gene expression profiles of 16,038 Plasmodium berghei sporozoites isolated throughout their development from midgut oocysts to salivary glands, and from forced salivation experiments. Our results reveal a succession of tightly regulated changes in gene expression occurring during the maturation of sporozoites and highlight candidate genes that could play important roles in oocyst egress, sporozoite motility, and the mechanisms underlying the invasion of mosquito salivary glands and mammalian hepatocytes. In addition, the single-cell data reveal extensive transcriptional heterogeneity among parasites isolated from the same anatomical site, suggesting that Plasmodium development in mosquitoes is asynchronous and regulated by intrinsic as well as environmental factors. Finally, our analyses show a decrease in transcriptional activity preceding the translational repression observed in mature sporozoites and associated with their quiescent state in salivary glands, followed by a rapid reactivation of the transcriptional machinery immediately upon salivation.

A Plasmodium falciparum RING Finger E3 Ubiquitin Ligase Modifies the Roles of PfMDR1 and PfCRT in Parasite Drug Responses.

Protein ubiquitination is an important posttranslational regulation mechanism that mediates Plasmodium development and modifies parasite responses to antimalarial drugs. Although mutations in several parasite ubiquitination enzymes have been linked to increased drug tolerance, the molecular mechanisms by which ubiquitination pathways mediate these parasite responses remain largely unknown. Here, we investigate the roles of a Plasmodium falciparum ring finger ubiquitin ligase (PfRFUL) in parasite development and in responses to antimalarial drugs. We engineered a transgenic parasite having the Pfrful gene tagged with an HA-2A-NeoR-glmS sequence to knockdown (KD) Pfrful expression using glucosamine (GlcN). A Western blot analysis of the proteins from GlcN-treated pSLI-HA-NeoR-glmS-tagged (PfRFULg) parasites, relative to their wild-type (Dd2) controls, showed changes in the ubiquitination of numerous proteins. PfRFUL KD rendered the parasites more sensitive to multiple antimalarial drugs, including mefloquine, piperaquine, amodiaquine, and dihydroartemisinin. PfRFUL KD also decreased the protein level of the P. falciparum multiple drug resistance 1 protein (PfMDR1) and altered the ratio of two bands of the P. falciparum chloroquine resistance transporter (PfCRT), suggesting contributions to the changed drug responses by the altered ubiquitination of these two molecules. The inhibition of proteasomal protein degradation by epoxomicin increased the PfRFUL level, suggesting the degradation of PfRFUL by the proteasome pathways, whereas the inhibition of E3 ubiquitin ligase activities by JNJ26854165 reduced the PfRFUL level. This study reveals the potential mechanisms of PfRFUL in modifying the expression of drug transporters and their roles in parasite drug responses. PfRFUL could be a potential target for antimalarial drug development.

No Association between the Plasmodium vivax crt-o MS334 or In9pvcrt Polymorphisms and Chloroquine Failure in a Pre-Elimination Clinical Cohort from Malaysia with a Large Clonal Expansion.

Increasing reports of resistance to a frontline malaria blood-stage treatment, chloroquine (CQ), raises concerns for the elimination of Plasmodium vivax. The absence of an effective molecular marker of CQ resistance in P. vivax greatly constrains surveillance of this emerging threat. A recent genetic cross between CQ sensitive (CQS) and CQ resistant (CQR) NIH-1993 strains of P. vivax linked a moderate CQR phenotype with two candidate markers in P. vivax CQ resistance transporter gene (pvcrt-o): MS334 and In9pvcrt. Longer TGAAGH motif lengths at MS334 were associated with CQ resistance, as were shorter motifs at the In9pvcrt locus. In this study, high-grade CQR clinical isolates of P. vivax from a low endemic setting in Malaysia were used to investigate the association between the MS334 and In9pvcrt variants and treatment efficacy. Among a total of 49 independent monoclonal P. vivax isolates assessed, high-quality MS334 and In9pvcrt sequences could be derived from 30 (61%) and 23 (47%), respectively. Five MS334 and six In9pvcrt alleles were observed, with allele frequencies ranging from 2 to 76% and 3 to 71%, respectively. None of the clinical isolates had the same variant as the NIH-1993 CQR strain, and none of the variants were associated with CQ treatment failure (all P > 0.05). Multi-locus genotypes (MLGs) at 9 neutral microsatellites revealed a predominant P. vivax strain (MLG6) accounting for 52% of Day 0 infections. The MLG6 strain comprised equal proportions of CQS and CQR infections. Our study reveals complexity in the genetic basis of CQ resistance in the Malaysian P. vivax pre-elimination setting and suggests that the proposed pvcrt-o MS334 and In9pvcrt markers are not reliable markers of CQ treatment efficacy in this setting. Further studies are needed in other endemic settings, applying hypothesis-free genome-wide approaches, and functional approaches to understand the biological impact of the TGAAGH repeats linked to CQ response in a cross are warranted to comprehend and track CQR P. vivax.

Detection of simple and complex de novo mutations with multiple reference sequences.

The characterization of de novo mutations in regions of high sequence and structural diversity from whole-genome sequencing data remains highly challenging. Complex structural variants tend to arise in regions of high repetitiveness and low complexity, challenging both de novo assembly, in which short reads do not capture the long-range context required for resolution, and mapping approaches, in which improper alignment of reads to a reference genome that is highly diverged from that of the sample can lead to false or partial calls. Long-read technologies can potentially solve such problems but are currently unfeasible to use at scale. Here we present Corticall, a graph-based method that combines the advantages of multiple technologies and prior data sources to detect arbitrary classes of genetic variant. We construct multisample, colored de Bruijn graphs from short-read data for all samples, align long-read-derived haplotypes and multiple reference data sources to restore graph connectivity information, and call variants using graph path-finding algorithms and a model for simultaneous alignment and recombination. We validate and evaluate the approach using extensive simulations and use it to characterize the rate and spectrum of de novo mutation events in 119 progeny from four Plasmodium falciparum experimental crosses, using long-read data on the parents to inform reconstructions of the progeny and to detect several known and novel nonallelic homologous recombination events.

Single-cell transcription analysis of Plasmodium vivax blood-stage parasites identifies stage- and species-specific profiles of expression.

Plasmodium vivax and P. falciparum, the parasites responsible for most human malaria worldwide, exhibit striking biological differences, which have important clinical consequences. Unfortunately, P. vivax, unlike P. falciparum, cannot be cultivated continuously in vitro, which limits our understanding of its biology and, consequently, our ability to effectively control vivax malaria. Here, we describe single-cell gene expression profiles of 9,215 P. vivax parasites from bloodstream infections of Aotus and Saimiri monkeys. Our results show that transcription of most P. vivax genes occurs during short periods of the intraerythrocytic cycle and that this pattern of gene expression is conserved in other Plasmodium species. However, we also identify a strikingly high proportion of species-specific transcripts in late schizonts, possibly associated with the specificity of erythrocyte invasion. Our findings provide new and robust markers of blood-stage parasites, including some that are specific to the elusive P. vivax male gametocytes, and will be useful for analyzing gene expression data from laboratory and field samples.

Transcriptome profiling of Plasmodium vivax in Saimiri monkeys identifies potential ligands for invasion.

Unlike the case in Asia and Latin America, Plasmodium vivax infections are rare in sub-Saharan Africa due to the absence of the Duffy blood group antigen (Duffy antigen), the only known erythrocyte receptor for the P. vivax merozoite invasion ligand, Duffy binding protein 1 (DBP1). However, P. vivax infections have been documented in Duffy-negative individuals throughout Africa, suggesting that P. vivax may use ligands other than DBP1 to invade Duffy-negative erythrocytes through other receptors. To identify potential P. vivax ligands, we compared parasite gene expression in Saimiri and Aotus monkey erythrocytes infected with P. vivax Salvador I (Sal I). DBP1 binds Aotus but does not bind to Saimiri erythrocytes; thus, P. vivax Sal I must invade Saimiri erythrocytes independent of DBP1. Comparing RNA sequencing (RNAseq) data for late-stage infections in Saimiri and Aotus erythrocytes when invasion ligands are expressed, we identified genes that belong to tryptophan-rich antigen and merozoite surface protein 3 (MSP3) families that were more abundantly expressed in Saimiri infections compared with Aotus infections. These genes may encode potential ligands responsible for P. vivax infections of Duffy-negative Africans.

Activity of Plasmodium vivax promoter elements in Plasmodium knowlesi, and a centromere-containing plasmid that expresses NanoLuc throughout the parasite life cycle.

BACKGROUND: Plasmodium knowlesi is now the major cause of human malaria in Malaysia, complicating malaria control efforts that must attend to the elimination of multiple Plasmodium species. Recent advances in the cultivation of P. knowlesi erythrocytic-stage parasites in vitro, transformation with exogenous DNA, and infection of mosquitoes with gametocytes from culture have opened up studies of this pathogen without the need for resource-intensive and costly non-human primate (NHP) models. For further understanding and development of methods for parasite transformation in malaria research, this study examined the activity of various trans-species transcriptional control sequences and the influence of Plasmodium vivax centromeric (pvcen) repeats in plasmid-transfected P. knowlesi parasites. METHODS: In vitro cultivated P. knowlesi parasites were transfected with plasmid constructs that incorporated Plasmodium vivax or Plasmodium falciparum 5' UTRs driving the expression of bioluminescence markers (firefly luciferase or Nanoluc). Promoter activities were assessed by bioluminescence, and parasites transformed with human resistant allele dihydrofolate reductase-expressing plasmids were selected using antifolates. The stability of transformants carrying pvcen-stabilized episomes was assessed by bioluminescence over a complete parasite life cycle through a rhesus macaque monkey, mosquitoes, and a second rhesus monkey. RESULTS: Luciferase expression assessments show that certain P. vivax promoter regions, not functional in the more evolutionarily-distant P. falciparum, can drive transgene expression in P. knowlesi. Further, pvcen repeats may improve the stability of episomal plasmids in P. knowlesi and support detection of NanoLuc-expressing elements over the full parasite life cycle from rhesus macaque monkeys to Anopheles dirus mosquitoes and back again to monkeys. In assays of drug responses to chloroquine, G418 and WR9910, anti-malarial half-inhibitory concentration (IC50) values of blood stages measured by NanoLuc activity proved comparable to IC50 values measured by the standard SYBR Green method. CONCLUSION: All three P. vivax promoters tested in this study functioned in P. knowlesi, whereas two of the three were inactive in P. falciparum. NanoLuc-expressing, centromere-stabilized plasmids may support high-throughput screenings of P. knowlesi for new anti-malarial agents, including compounds that can block the development of mosquito- and/or liver-stage parasites.

Artemisinin resistance phenotypes and K13 inheritance in a Plasmodium falciparum cross and Aotus model.

Concerns about malaria parasite resistance to treatment with artemisinin drugs (ARTs) have grown with findings of prolonged parasite clearance t1/2s (>5 h) and their association with mutations in Plasmodium falciparum Kelch-propeller protein K13. Here, we describe a P. falciparum laboratory cross of K13 C580Y mutant with C580 wild-type parasites to investigate ART response phenotypes in vitro and in vivo. After genotyping >400 isolated progeny, we evaluated 20 recombinants in vitro: IC50 measurements of dihydroartemisinin were at similar low nanomolar levels for C580Y- and C580-type progeny (mean ratio, 1.00; 95% CI, 0.62-1.61), whereas, in a ring-stage survival assay, the C580Y-type progeny had 19.6-fold (95% CI, 9.76-39.2) higher average counts. In splenectomized Aotus monkeys treated with three daily doses of i.v. artesunate, t1/2 calculations by three different methods yielded mean differences of 0.01 h (95% CI, -3.66 to 3.67), 0.80 h (95% CI, -0.92 to 2.53), and 2.07 h (95% CI, 0.77-3.36) between C580Y and C580 infections. Incidences of recrudescence were 57% in C580Y (4 of 7) versus 70% in C580 (7 of 10) infections (-13% difference; 95% CI, -58% to 35%). Allelic substitution of C580 in a C580Y-containing progeny clone (76H10) yielded a transformant (76H10C580Rev) that, in an infected monkey, recrudesced regularly 13 times over 500 d. Frequent recrudescences of ART-treated P. falciparum infections occur with or without K13 mutations and emphasize the need for improved partner drugs to effectively eliminate the parasites that persist through the ART component of combination therapy.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution.

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Kelch Mutations in Plasmodium falciparum Protein K13 Do Not Modulate Dormancy after Artemisinin Exposure and Sorbitol Selection In Vitro.

Some Kelch mutations of the Plasmodium falciparum K13 protein confer increased survival to dihydroartemisinin (DHA)-treated ring-stage parasites. Here, we asked if K13 mutations affect a dormancy phenotype allowing parasites to survive DHA exposure and then sorbitol selection. Although recrudescence from dormancy differed between two distinct parasite lines, it was similar for isogenic lines carrying single-site substitutions in K13. Therefore, K13 mutations do not alter the DHA-sorbitol combined dormancy phenotype; rather, traits from other loci likely determine this phenotype.

A set of microsatellite markers to differentiate Plasmodium falciparum progeny of four genetic crosses.

BACKGROUND: Four Plasmodium falciparum genetic crosses (HB3×3D7, HB3×Dd2, 7G8×GB4, and 803×GB4) have produced sets of recombinant progeny that are widely used for malaria research, including investigations of anti-malarial drug resistance. It is critical to maintain the progeny free from cross-contamination. Microsatellite polymorphisms can be used to validate parasite identity. RESULTS: A set of 12 markers was developed that differentiates the parents of the four P. falciparum crosses. This typing set identified distinguishing patterns of inheritance (fingerprints) in segregant collections of 15 progeny clones from HB3×3D7, 32 from HB3×Dd2, 33 from 7G8×GB4, and 81 from 803×GB4. Stronger amplification was observed with shorter relative to longer alleles of individual microsatellites. In experiments with mixed parental DNAs, electropherograms showed that signals of cross-contamination can be missed when minor peaks less than 1/4 or 1/3 the height of the major peak are disregarded by threshold settings commonly used for population studies. CONCLUSIONS: Microsatellite typing is an effective method to check the identity of P. falciparum lines and detect parasite cross-contamination in cultures; however, care must be taken not to ignore minor peaks that can be overlooked. The 12 microsatellite markers presented here provide a rapid and efficient means to distinguish the segregants of laboratory crosses. Fingerprint patterns from these markers are useful to maintain the integrity of diverse parasite lines in and between research laboratories.

A single nucleotide polymorphism in the Plasmodium falciparum atg18 gene associates with artemisinin resistance and confers enhanced parasite survival under nutrient deprivation.

BACKGROUND: Artemisinin-resistant Plasmodium falciparum has been reported throughout the Greater Mekong subregion and threatens to disrupt current malaria control efforts worldwide. Polymorphisms in kelch13 have been associated with clinical and in vitro resistance phenotypes; however, several studies suggest that the genetic determinants of resistance may involve multiple genes. Current proposed mechanisms of resistance conferred by polymorphisms in kelch13 hint at a connection to an autophagy-like pathway in P. falciparum. RESULTS: A SNP in autophagy-related gene 18 (atg18) was associated with long parasite clearance half-life in patients following artemisinin-based combination therapy. This gene encodes PfAtg18, which is shown to be similar to the mammalian/yeast homologue WIPI/Atg18 in terms of structure, binding abilities, and ability to form puncta in response to stress. To investigate the contribution of this polymorphism, the atg18 gene was edited using CRISPR/Cas9 to introduce a T38I mutation into a k13-edited Dd2 parasite. The presence of this SNP confers a fitness advantage by enabling parasites to grow faster in nutrient-limited settings. The mutant and parent parasites were screened against drug libraries of 6349 unique compounds. While the SNP did not modulate the parasite's susceptibility to any of the anti-malarial compounds using a 72-h drug pulse, it did alter the parasite's susceptibility to 227 other compounds. CONCLUSIONS: These results suggest that the atg18 T38I polymorphism may provide additional resistance against artemisinin derivatives, but not partner drugs, even in the absence of kelch13 mutations, and may also be important in parasite survival during nutrient deprivation.

Alternative methods for the Plasmodium falciparum artemisinin ring-stage survival assay with increased simplicity and parasite stage-specificity.

BACKGROUND: Artemisinin-based combination therapy is recommended to treat Plasmodium falciparum worldwide, but observations of longer artemisinin (ART) parasite clearance times (PCTs) in Southeast Asia are widely interpreted as a sign of potential ART resistance. In search of an in vitro correlate of in vivo PCT after ART treatment, a ring-stage survival assay (RSA) of 0-3 h parasites was developed and linked to polymorphisms in the Kelch propeller protein (K13). However, RSA remains a laborious process, involving heparin, Percoll gradient, and sorbitol treatments to obtain rings in the 0-3 h window. Here two alternative RSA protocols are presented and compared to the standard Percoll-based method, one highly stage-specific and one streamlined for laboratory application. METHODS: For all protocols, P. falciparum cultures were synchronized with 5 % sorbitol treatment twice over two intra-erythrocytic cycles. For a filtration-based RSA, late-stage schizonts were passed through a 1.2 µm filter to isolate merozoites, which were incubated with uninfected erythrocytes for 45 min. The erythrocytes were then washed to remove lysis products and further incubated until 3 h post-filtration. Parasites were pulsed with either 0.1 % dimethyl sulfoxide (DMSO) or 700 nM dihydroartemisinin in 0.1 % DMSO for 6 h, washed twice in drug-free media, and incubated for 66-90 h, when survival was assessed by microscopy. For a sorbitol-only RSA, synchronized young (0-3 h) rings were treated with 5 % sorbitol once more prior to the assay and adjusted to 1 % parasitaemia. The drug pulse, incubation, and survival assessment were as described above. RESULTS: Ring-stage survival of P. falciparum parasites containing either the K13 C580 or C580Y polymorphism (associated with low and high RSA survival, respectively) were assessed by the described filtration and sorbitol-only methods and produced comparable results to the reported Percoll gradient RSA. Advantages of both new methods include: fewer reagents, decreased time investment, and fewer procedural steps, with enhanced stage-specificity conferred by the filtration method. CONCLUSIONS: Assessing P. falciparum ART sensitivity in vitro via RSA can be streamlined and accurately evaluated in the laboratory by filtration or sorbitol synchronization methods, thus increasing the accessibility of the assay to research groups.

Editing the Plasmodium vivax genome, using zinc-finger nucleases.

Plasmodium vivax is a major cause of malaria morbidity worldwide yet has remained genetically intractable. To stably modify this organism, we used zinc-finger nucleases (ZFNs), which take advantage of homology-directed DNA repair mechanisms at the site of nuclease action. Using ZFNs specific to the gene encoding P. vivax dihydrofolate reductase (pvdhfr), we transfected blood specimens from Saimiri boliviensis monkeys infected with the pyrimethamine (Pyr)-susceptible Chesson strain with a ZFN plasmid carrying a Pyr-resistant mutant pvdhfr sequence. We obtained Pyr-resistant parasites in vivo that carried mutant pvdhfr and additional silent mutations designed to confirm editing. These results herald the era of stable P. vivax genetic modifications.

Prevalence of Plasmodium falciparum anti-malarial resistance-associated polymorphisms in pfcrt, pfmdr1 and pfnhe1 in Muheza, Tanzania, prior to introduction of artemisinin combination therapy.

BACKGROUND: A report of the chloroquine and amodiaquine resistance pfcrt-SVMNT haplotype in Tanzania raises concern about high-level resistance to the artesunate-amodiaquine combination treatment widely employed in Africa. Mutations in the pfmdr1 multi-drug resistance gene may also be associated with resistance, and a highly polymorphic microsatellite (ms-4760) of the pfnhe1 gene involved in quinine susceptibility has not been surveyed in Tanzania. METHODS: A total of 234 samples collected between 2003 - 2006 from an observational birth cohort of young children in Muheza, Tanzania were analysed. In these children, 141 cases of P. falciparum infections were treated with AQ and 93 episodes were treated with QN. Haplotypes of pfcrt and pfmdr1 were determined by a Taqman assay, and ms-4760 repeats in pfnhe1 were assessed by nested PCR amplification and direct sequencing. Parasite population diversity was evaluated using microsatellite markers on five different chromosomes. RESULTS: The pfcrt-CVIET haplotype was present alone in 93.6% (219/234) of the samples over the study period; the wild-type chloroquine- and amodiaquine-sensitive haplotype pfcrt-CVMNK was present in 4.3% (10/234) of the samples; and both haplotypes were present in 2.1% (5/234) of the samples. No significant change in wild-type pfcrt-CVMNK prevalence was evident over the 4-year period of the study. The pfcrt-SVMNT haplotype associated with high-level amodiaquine resistance was not detected in this study. The pfmdr1 locus was genotyped in 178 of these samples. The pfmdr1-YYNY haplotype predominated in 67.4% (120/178) of infections and was significantly associated with the pfcrt-CVIET haplotype. All samples carried the wild-type pfmdr1-N1042 codon. The ms-4760 repeat on pfnhe1 locus displayed 12 distinct haplotypes with ms-4760-1 predominating in the population. Analysis of these haplotypes showed no association of a particular haplotype with quinine treatment outcome. CONCLUSION: The pfcrt-CVIET chloroquine resistance haplotype dominated in the collection of P. falciparum samples from Muheza. The pfcrt-SVMNT haplotype, which threatens the efficacy of amodiaquine and was reported in the same time period from Korogwe, Tanzania, 40 Km from Muheza, was not detected. Relative low prevalence of pfcrt-SVMNT in Africa may result from genetic or other factors rendering P. falciparum less supportive of this haplotype than in South America or other regions. TRIAL REGISTRATION: Trial Protocol Number: 08-I-N064.

Plasmodium falciparum clearance rates in response to artesunate in Malian children with malaria: effect of acquired immunity.

BACKGROUND: Artemisinin resistance, a long parasite clearance half-life in response to artemisinin, has been described in patients with Plasmodium falciparum malaria in southeast Asia. Few baseline half-lives have been reported from Africa, where artemisinins were recently introduced. METHODS: We treated P. falciparum malaria in 215 Malian children aged 0.5-15 years with artesunate (0, 24, 48 hours) and amodiaquine (72, 96, 120 hours). We estimated half-life by measuring parasite density every 6 hours until undetectable and evaluated the effects of age, sex, ethnicity, and red blood cell (RBC) polymorphisms on half-life. We quantified the proportion of parasitized RBCs recognized by autologous immunoglobulin G (IgG). RESULTS: The geometric mean half-life was 1.9 hours (95% confidence interval, 1.8-2.0) and did not correlate with parasite ex vivo susceptibility to artemisinins. In a linear model accounting for host factors, half-life decreased by 4.1 minutes for every 1-year increase in age. The proportion of parasitized RBCs recognized by IgG correlated inversely with half-life (r = -0.475; P = .0006). CONCLUSIONS: Parasite clearance in response to artesunate is faster in Mali than in southeast Asia. IgG responses to parasitized RBCs shorten half-life and may influence this parameter in areas where age is not an adequate surrogate of immunity and correlates of parasite-clearing immunity have not been identified. CLINICAL TRIALS REGISTRATION: NCT00669084.

Progressive heterogeneity of enlarged and irregularly shaped apicoplasts in P. falciparum persister blood stages after drug treatment.

Morphological modifications and shifts in organelle relationships are hallmarks of dormancy in eukaryotic cells. Communications between altered mitochondria and nuclei are associated with metabolic quiescence of cancer cells that can survive chemotherapy. In plants, changes in the pathways between nuclei, mitochondria, and chloroplasts are associated with cold stress and bud dormancy. Plasmodium falciparum parasites, the deadliest agent of malaria in humans, contain a chloroplast-like organelle (apicoplast) derived from an ancient photosynthetic symbiont. Antimalarial treatments can fail because a small fraction of the blood stage parasites enter dormancy and recrudesce after drug exposure. Altered mitochondrial-nuclear interactions in these persisters have been described for P. falciparum, but interactions of the apicoplast remained to be characterized. In the present study, we examined the apicoplasts of dormant persisters obtained after exposure to dihydroartemisinin (a first-line antimalarial drug) followed by sorbitol treatment, or after exposure to sorbitol treatment alone. As previously observed, the mitochondrion of persisters was consistently enlarged and in close association with the nucleus. In contrast, the apicoplast varied from compact and oblate, like those of active ring stage parasites, to enlarged and irregularly shaped. Enlarged apicoplasts became more prevalent later in dormancy, but regular size apicoplasts subsequently predominated when actively replicating parasites recrudesced. All three organelles, nucleus, mitochondrion, and apicoplast, became closer during dormancy. Understanding their relationships in erythrocytic-stage persisters may lead to new strategies to prevent recrudescences and protect the future of malaria chemotherapy.

Single-cell analyses of polyclonal Plasmodium vivax infections and their consequences on parasite transmission.

Most Plasmodium vivax infections contain genetically distinct parasites, but the consequences of this polyclonality on the development of asexual parasites, their sexual differentiation, and their transmission remain unknown. We describe infections of Saimiri monkeys with two strains of P. vivax and the analyses of 117,350 parasites characterized by single cell RNA sequencing and individually genotyped. In our model, consecutive inoculations fail to establish polyclonal infections. By contrast, simultaneous inoculations of two strains lead to sustained polyclonal infections, although without detectable differences in parasite regulation or sexual commitment. Analyses of sporozoites dissected from mosquitoes fed on coinfected monkeys show that all genotypes are successfully transmitted to mosquitoes. However, after sporozoite inoculation, not all genotypes contribute to the subsequent blood infections, highlighting an important bottleneck during pre-erythrocytic development. Overall, these studies provide new insights on the mechanisms regulating the establishment of polyclonal P. vivax infections and their consequences for disease transmission.

Higher microsatellite diversity in Plasmodium vivax than in sympatric Plasmodium falciparum populations in Pursat, Western Cambodia.

Previous microsatellite analyses of sympatric populations of Plasmodium vivax and Plasmodium falciparum in Brazil revealed higher diversity in the former species. However, it remains unclear whether regional species-specific differences in prevalence and transmission levels might account for these findings. Here, we examine sympatric populations of P. vivax (n=87) and P. falciparum (n=164) parasites from Pursat province, Western Cambodia, where both species are similarly prevalent. Using 10 genome-wide microsatellites for P. falciparum and 13 for P. vivax, we found that the P. vivax population was more diverse than the sympatric P. falciparum population (average virtual heterozygosity [HE], 0.87 vs. 0.66, P=0.003), with more multiple-clone infections (89.6% vs. 47.6%) and larger mean number of alleles per marker (16.2 vs. 11.1, P=0.07). Both populations showed significant multi-locus linkage disequilibrium suggestive of a predominantly clonal mode of parasite reproduction. The higher microsatellite diversity found in P. vivax isolates, compared to sympatric P. falciparum isolates, does not necessarily result from local differences in transmission level and may reflect differences in population history between species or increased mutation rates in P. vivax.