Single-genome sequencing reveals within-host evolution of human malaria parasites

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.

MalariaSED: a deep learning framework to decipher the regulatory contributions of noncoding variants in malaria parasites

Malaria remains one of the deadliest infectious diseases. Transcriptional regulation effects of noncoding variants in this unusual genome of malaria parasites remain elusive. We developed a sequence-based, ab initio deep learning framework, MalariaSED, for predicting chromatin profiles in malaria parasites. The MalariaSED performance was validated by published ChIP-qPCR and TF motifs results. Applying MalariaSED to ~ 1.3 million variants shows that geographically differentiated noncoding variants are associated with parasite invasion and drug resistance. Further analysis reveals chromatin accessibility changes at Plasmodium falciparum rings are partly associated with artemisinin resistance. MalariaSED illuminates the potential functional roles of noncoding variants in malaria parasites.

Intra-host Trypanosoma cruzi strain dynamics shape disease progression: the missing link in Chagas disease pathogenesis

Chronic Chagasic cardiomyopathy develops years after infection in 20-40% of patients, but disease progression is poorly understood. Here, we assessed Trypanosoma cruzi parasite dynamics and pathogenesis over a 2.5-year period in naturally infected rhesus macaques. Individuals with better control of parasitemia were infected with a greater diversity of parasite strains compared to those with increasing parasitemia over time. Also, the in vivo parasite multiplication rate decreased with increasing parasite diversity, suggesting competition among strains or a stronger immune response in multiple infections. Significant differences in electrocardiographic (ECG) profiles were observed in Chagasic macaques compared to uninfected controls, suggesting early conduction defects, and changes in ECG patterns over time were observed only in macaques with increasing parasitemia and lower parasite diversity. Disease progression was also associated with plasma fibronectin degradation, which may serve as a biomarker. These data provide a novel framework for the understanding of Chagas disease pathogenesis, with parasite diversity shaping disease progression.IMPORTANCEChagas disease progression remains poorly understood, and patients at increased risk of developing severe cardiac disease cannot be distinguished from those who may remain asymptomatic. Monitoring of Trypanosoma cruzi strain dynamics and pathogenesis over 2-3 years in naturally infected macaques shows that increasing parasite diversity in hosts is detrimental to parasite multiplication and Chagasic cardiomyopathy disease progression. This provides a novel framework for the understanding of Chagas disease pathogenesis.

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.

Epigenetically regulated RNA-binding proteins signify malaria hypnozoite dormancy

Dormancy enables relapsing malaria parasites, such as Plasmodium vivax and cynomolgi, to survive unfavorable conditions. It is enabled by hypnozoites, parasites remaining quiescent inside hepatocytes before reactivating and establishing blood-stage infection. We integrate omics approaches to explore gene-regulatory mechanisms underlying hypnozoite dormancy. Genome-wide profiling of activating and repressing histone marks identifies a few genes that get silenced by heterochromatin during hepatic infection of relapsing parasites. By combining single-cell transcriptomics, chromatin accessibility profiling, and fluorescent in situ RNA hybridization, we show that these genes are expressed in hypnozoites and that their silencing precedes parasite development. Intriguingly, these hypnozoite-specific genes mainly encode proteins with RNA-binding domains. We hence hypothesize that these likely repressive RNA-binding proteins keep hypnozoites in a developmentally competent but dormant state and that heterochromatin-mediated silencing of the corresponding genes aids reactivation. Exploring the regulation and exact function of these proteins hence could provide clues for targeted reactivation and killing of these latent pathogens.

Stochastic expression of invasion genes in Plasmodium falciparum schizonts

Genetically identical cells are known to exhibit differential phenotypes in the same environmental conditions. These phenotypic variants are linked to transcriptional stochasticity and have been shown to contribute towards adaptive flexibility of a wide range of unicellular organisms. Here, we investigate transcriptional heterogeneity and stochastic gene expression in Plasmodium falciparum by performing the quasilinear multiple annealing and looping based amplification cycles (MALBAC) based amplification and single cell RNA sequencing of blood stage schizonts. Our data reveals significant transcriptional variations in the schizont stage with a distinct group of highly variable invasion gene transcripts being identified. Moreover, the data reflects several diversification processes including putative developmental “checkpoint”; transcriptomically distinct parasite sub-populations and transcriptional switches in variable gene families (var, rifin, phist). Most of these features of transcriptional variability are preserved in isogenic parasite cell populations (albeit with a lesser amplitude) suggesting a role of epigenetic factors in cell-to-cell transcriptional variations in human malaria parasites. Lastly, we apply quantitative RT-PCR and RNA-FISH approach and confirm stochastic expression of key invasion genes, such as, msp1, msp3, msp7, eba181 and ama1 which represent prime candidates for invasion-blocking vaccines.

Genetically identical cells can be phenotypically diverse to allow adaptive flexibility in a given environment. This phenotypic diversity is driven by epigenetic and transcriptional variability. Here, Tripathi et al. perform scRNA-seq of isogenic and non-isogenic Plasmodium falciparum schizont populations to explore transcriptional heterogeneity and stochastic gene expression during the course of development.

An open dataset of Plasmodium vivax genome variation in 1,895 worldwide samples

This report describes the MalariaGEN Pv4 dataset, a new release of curated genome variation data on 1,895 samples of Plasmodium vivax collected at 88 worldwide locations between 2001 and 2017. It includes 1,370 new samples contributed by MalariaGEN and VivaxGEN partner studies in addition to previously published samples from these and other sources. We provide genotype calls at over 4.5 million variable positions including over 3 million single nucleotide polymorphisms (SNPs), as well as short indels and tandem duplications. This enlarged dataset highlights major compartments of parasite population structure, with clear differentiation between Africa, Latin America, Oceania, Western Asia and different parts of Southeast Asia. Each sample has been classified for drug resistance to sulfadoxine, pyrimethamine and mefloquine based on known markers at the dhfr, dhps and mdr1 loci. The prevalence of all of these resistance markers was much higher in Southeast Asia and Oceania than elsewhere. This open resource of analysis-ready genome variation data from the MalariaGEN and VivaxGEN networks is driven by our collective goal to advance research into the complex biology of P. vivax and to accelerate genomic surveillance for malaria control and elimination.

Decrypting the complexity of the human malaria parasite biology through systems biology approaches

The human malaria parasite, Plasmodium falciparum, is a unicellular protozoan responsible for over half a million deaths annually. With a complex life cycle alternating between human and invertebrate hosts, this apicomplexan is notoriously adept at evading host immune responses and developing resistance to all clinically administered treatments. Advances in omics-based technologies, increased sensitivity of sequencing platforms and enhanced CRISPR based gene editing tools, have given researchers access to more in-depth and untapped information about this enigmatic micro-organism, a feat thought to be infeasible in the past decade. Here we discuss some of the most important scientific achievements made over the past few years with a focus on novel technologies and platforms that set the stage for subsequent discoveries. We also describe some of the systems-based methods applied to uncover gaps of knowledge left through single-omics applications with the hope that we will soon be able to overcome the spread of this life-threatening disease.

Complex infections in vivax malaria: the more you look, the more you find

The human malaria parasite Plasmodium vivax commonly causes complex multiclonal infections. Recently, Dia et al. have developed innovative methods for single-cell sequencing (SCS) of P. vivax infections by adapting an approach used previously for Plasmodium falciparum. Their studies provide fascinating new insights into P. vivax intrahost diversity and evolution.