Silence, activate, poise and switch! Mechanisms of antigenic variation in Plasmodium falciparum

Deciphering the Plasmodium falciparum perinuclear var gene expression site

The protozoan parasite Plasmodium falciparum, responsible for the deadliest form of human malaria, employs antigenic variation via monoallelic expression as a key survival strategy. The selective activation of one out of the 60-member var gene family is key to understanding the parasite's ability to cause severe disease and evade the host immune response. var gene activation is initiated by its relocation to a specialized expression site. While the perinuclear expression site (PES) plays a crucial role in enabling the expression of a single allele, the characteristics of this PES remain largely obscure. Recent breakthroughs in genome editing tools and the discovery of regulatory noncoding RNAs have shed light on this intriguing biological feature, offering significant insights into the mechanisms of pathogen virulence.

RNA polymerase III is involved in regulating Plasmodium falciparum virulence

While often undetected and untreated, persistent seasonal asymptomatic malaria infections remain a global public health problem. Despite the presence of parasites in the peripheral blood, no symptoms develop. Disease severity is correlated with the levels of infected red blood cells (iRBCs) adhering within blood vessels. Changes in iRBC adhesion capacity have been linked to seasonal asymptomatic malaria infections, however how this is occurring is still unknown. Here, we present evidence that RNA polymerase III (RNA Pol III) transcription in Plasmodium falciparum is downregulated in field isolates obtained from asymptomatic individuals during the dry season. Through experiments with in vitro cultured parasites, we have uncovered an RNA Pol III-dependent mechanism that controls pathogen proliferation and expression of a major virulence factor in response to external stimuli. Our findings establish a connection between P. falciparum cytoadhesion and a non-coding RNA family transcribed by Pol III. Additionally, we have identified P. falciparum Maf1 as a pivotal regulator of Pol III transcription, both for maintaining cellular homeostasis and for responding adaptively to external signals. These results introduce a novel perspective that contributes to our understanding of P. falciparum virulence. Furthermore, they establish a connection between this regulatory process and the occurrence of seasonal asymptomatic malaria infections.

A Study of Cross-Protection between Eimeria maxima Immunovariants

For reasons unknown, Eimeria maxima is unique among Eimeria species infecting chickens in the immunovariability it displays among isolates from different geographical areas. Eimeria maxima oocysts (named EmaxAPU3) were isolated late in grow-out (6 weeks) from litter in a commercial broiler operation that was using Eimeria vaccination as the coccidiosis control program. Cross-protection studies (n = 4) were conducted in immunologically naïve chickens between EmaxAPU3 and two E. maxima lab strains (EmaxAPU1, EmaxAPU2) by immunizing with one E. maxima strain and challenging with either the homologous or heterologous E. maxima. As measured by oocyst output, immunization with EmaxAPU1 protected against homologous challenge (EmaxAPU1) and against heterologous challenge with EmaxAPU3, but not against EmaxAPU2. Similarly, immunization with EmaxAPU3 protected against homologous challenge (EmaxAPU3) and against heterologous challenge with EmaxAPU1, but not against EmaxAPU2. Immunization of chickens with EmaxAPU2 elicited a protective response against homologous challenge (EmaxAPU2), but not against EmaxAPU1 nor EmaxAPU3. The most plausible explanation for the appearance of this immunovariant late in grow-out is that E. maxima APU3 escaped immunity directed to E. maxima antigenic types in the commercial vaccine.

Nonsense-mediated decay machinery in Plasmodium falciparum is inefficient and non-essential

Nonsense-mediated decay (NMD) is a conserved mRNA quality control process that eliminates transcripts bearing a premature termination codon. In addition to its role in removing erroneous transcripts, NMD is involved in post-transcriptional regulation of gene expression via programmed intron retention in metazoans. The apicomplexan parasite Plasmodium falciparum shows relatively high levels of intron retention, but it is unclear whether these variant transcripts are functional targets of NMD. In this study, we use CRISPR-Cas9 to disrupt and epitope-tag the P. falciparum orthologs of two core NMD components: PfUPF1 (PF3D7_1005500) and PfUPF2 (PF3D7_0925800). We localize both PfUPF1 and PfUPF2 to puncta within the parasite cytoplasm and show that these proteins interact with each other and other mRNA-binding proteins. Using RNA-seq, we find that although these core NMD orthologs are expressed and interact in P. falciparum, they are not required for degradation of nonsense transcripts. Furthermore, our work suggests that the majority of intron retention in P. falciparum has no functional role and that NMD is not required for parasite growth ex vivo. IMPORTANCE In many organisms, the process of destroying nonsense transcripts is dependent on a small set of highly conserved proteins. We show that in the malaria parasite, these proteins do not impact the abundance of nonsense transcripts. Furthermore, we demonstrate efficient CRISPR-Cas9 editing of the malaria parasite using commercial Cas9 nuclease and synthetic guide RNA, streamlining genomic modifications in this genetically intractable organism.

Knobs, Adhesion, and Severe Falciparum Malaria

Plasmodium falciparum can cause a severe disease with high mortality. A major factor contributing to the increased virulence of P. falciparum, as compared to other human malarial parasites, is the sequestration of infected erythrocytes in the capillary beds of organs and tissues. This sequestration is due to the cytoadherence of infected erythrocytes to endothelial cells. Cytoadherence is primarily mediated by a parasite protein expressed on the surface of the infected erythrocyte called P. falciparum erythrocyte membrane protein-1 (PfEMP1). PfEMP1 is embedded in electron-dense protuberances on the surface of the infected erythrocytes called knobs. These knobs are assembled on the erythrocyte membrane via exported parasite proteins, and the knobs function as focal points for the cytoadherence of infected erythrocytes to endothelial cells. PfEMP1 is a member of the var gene family, and there are approximately 60 antigenically distinct PfEMP1 alleles per parasite genome. Var gene expression exhibits allelic exclusion, with only a single allele being expressed by an individual parasite. This results in sequential waves of antigenically distinct infected erythrocytes and this antigenic variation allows the parasite to establish long-term chronic infections. A wide range of endothelial cell receptors can bind to the various PfEMP1 alleles, and thus, antigenic variation also results in a change in the cytoadherence phenotype. The cytoadherence phenotype may result in infected erythrocytes sequestering in different tissues and this difference in sequestration may explain the wide range of possible clinical manifestations associated with severe falciparum malaria.

Antibodies to variable surface antigens induce antigenic variation in the intestinal parasite Giardia lamblia

The genomes of most protozoa encode families of variant surface antigens. In some parasitic microorganisms, it has been demonstrated that mutually exclusive changes in the expression of these antigens allow parasites to evade the host's immune response. It is widely assumed that antigenic variation in protozoan parasites is accomplished by the spontaneous appearance within the population of cells expressing antigenic variants that escape antibody-mediated cytotoxicity. Here we show, both in vitro and in animal infections, that antibodies to Variant-specific Surface Proteins (VSPs) of the intestinal parasite Giardia lamblia are not cytotoxic, inducing instead VSP clustering into liquid-ordered phase membrane microdomains that trigger a massive release of microvesicles carrying the original VSP and switch in expression to different VSPs by a calcium-dependent mechanism. This novel mechanism of surface antigen clearance throughout its release into microvesicles coupled to the stochastic induction of new phenotypic variants not only changes current paradigms of antigenic switching but also provides a new framework for understanding the course of protozoan infections as a host/parasite adaptive process.

Babesia duncani multi-omics identifies virulence factors and drug targets

Babesiosis is a malaria-like disease in humans and animals that is caused by Babesia species, which are tick-transmitted apicomplexan pathogens. Babesia duncani causes severe to lethal infection in humans, but despite the risk that this parasite poses as an emerging pathogen, little is known about its biology, metabolic requirements or pathogenesis. Unlike other apicomplexan parasites that infect red blood cells, B. duncani can be continuously cultured in vitro in human erythrocytes and can infect mice resulting in fulminant babesiosis and death. We report comprehensive, detailed molecular, genomic, transcriptomic and epigenetic analyses to gain insights into the biology of B. duncani. We completed the assembly, 3D structure and annotation of its nuclear genome, and analysed its transcriptomic and epigenetics profiles during its asexual life cycle stages in human erythrocytes. We used RNA-seq data to produce an atlas of parasite metabolism during its intraerythrocytic life cycle. Characterization of the B. duncani genome, epigenome and transcriptome identified classes of candidate virulence factors, antigens for diagnosis of active infection and several attractive drug targets. Furthermore, metabolic reconstitutions from genome annotation and in vitro efficacy studies identified antifolates, pyrimethamine and WR-99210 as potent inhibitors of B. duncani to establish a pipeline of small molecules that could be developed as effective therapies for the treatment of human babesiosis.

The RRM-mediated RNA binding activity in T. brucei RAP1 is essential for VSG monoallelic expression

Trypanosoma brucei is a protozoan parasite that causes human African trypanosomiasis. Its major surface antigen VSG is expressed from subtelomeric loci in a strictly monoallelic manner. We previously showed that the telomere protein TbRAP1 binds dsDNA through its 737RKRRR741 patch to silence VSGs globally. How TbRAP1 permits expression of the single active VSG is unknown. Through NMR structural analysis, we unexpectedly identify an RNA Recognition Motif (RRM) in TbRAP1, which is unprecedented for RAP1 homologs. Assisted by the 737RKRRR741 patch, TbRAP1 RRM recognizes consensus sequences of VSG 3'UTRs in vitro and binds the active VSG RNA in vivo. Mutating conserved RRM residues abolishes the RNA binding activity, significantly decreases the active VSG RNA level, and derepresses silent VSGs. The competition between TbRAP1's RNA and dsDNA binding activities suggests a VSG monoallelic expression mechanism in which the active VSG's abundant RNA antagonizes TbRAP1's silencing effect, thereby sustaining its full-level expression.

Patterns of Heterochromatin Transitions Linked to Changes in the Expression of Plasmodium falciparum Clonally Variant Genes

The survival of malaria parasites in the changing human blood environment largely depends on their ability to alter gene expression by epigenetic mechanisms. The active state of Plasmodium falciparum clonally variant genes (CVGs) is associated with euchromatin characterized by the histone mark H3K9ac, whereas the silenced state is characterized by H3K9me3-based heterochromatin. Expression switches are linked to euchromatin-heterochromatin transitions, but these transitions have not been characterized for the majority of CVGs. To define the heterochromatin distribution patterns associated with the alternative transcriptional states of CVGs, we compared H3K9me3 occupancy at a genome-wide level among several parasite subclones of the same genetic background that differed in the transcriptional state of many CVGs. We found that de novo heterochromatin formation or the complete disruption of a heterochromatin domain is a relatively rare event, and for the majority of CVGs, expression switches can be explained by the expansion or retraction of heterochromatin domains. We identified different modalities of heterochromatin changes linked to transcriptional differences, but despite this complexity, heterochromatin distribution patterns generally enable the prediction of the transcriptional state of specific CVGs. We also found that in some subclones, several var genes were simultaneously in an active state. Furthermore, the heterochromatin levels in the putative regulatory region of the gdv1 antisense noncoding RNA, a regulator of sexual commitment, varied between parasite lines with different sexual conversion rates. IMPORTANCE The malaria parasite P. falciparum is responsible for more than half a million deaths every year. P. falciparum clonally variant genes (CVGs) mediate fundamental host-parasite interactions and play a key role in parasite adaptation to fluctuations in the conditions of the human host. The expression of CVGs is regulated at the epigenetic level by changes in the distribution of a type of chromatin called heterochromatin. Here, we describe at a genome-wide level the changes in the heterochromatin distribution associated with the different transcriptional states of CVGs. Our results also reveal a likely role for heterochromatin at a particular locus in determining the parasite investment in transmission to mosquitoes. Additionally, this data set will enable the prediction of the transcriptional state of CVGs from epigenomic data, which is important for the study of parasite adaptation to the conditions of the host in natural malaria infections.

Global gene expression of human malaria parasite liver stages throughout intrahepatocytic development

Plasmodium falciparum (Pf) is causing the greatest malaria burden, yet the liver stages (LS) of this most important parasite species have remained poorly studied. Here, we used a human liver-chimeric mouse model in combination with a novel fluorescent PfNF54 parasite line (PfNF54cspGFP) to isolate PfLS-infected hepatocytes and generate transcriptomes that cover the major LS developmental phases in human hepatocytes. RNA-seq analysis of early Pf LS trophozoites two days after infection, revealed a central role of translational regulation in the transformation of the extracellular invasive sporozoite into intracellular LS. The developmental time course gene expression analysis indicated that fatty acid biosynthesis, isoprenoid biosynthesis and iron metabolism are sustaining LS development along with amino acid metabolism and biosynthesis. Countering oxidative stress appears to play an important role during intrahepatic LS development. Furthermore, we observed expression of the variant PfEMP1 antigen-encoding var genes, and we confirmed expression of PfEMP1 protein during LS development. Transcriptome comparison of the late Pf liver stage schizonts with P. vivax (Pv) late liver stages revealed highly conserved gene expression profiles among orthologous genes. A notable difference however was the expression of genes regulating sexual stage commitment. While Pv schizonts expressed markers of sexual commitment, the Pf LS parasites were not sexually committed and showed expression of gametocytogenesis repression factors. Our results provide the first comprehensive gene expression profile of the human malaria parasite Pf LS isolated during in vivo intrahepatocytic development. This data will inform biological studies and the search for effective intervention strategies that can prevent infection.

Detection of positive selection acting on protein surfaces at the whole-genome scale in the human malaria parasite Plasmodium falciparum

The host-parasite evolutionary arms race is a fundamental process with medical implications. During this process, the host develops parasite resistance, and the parasite develops host immune evasion strategies. Thus, this process accelerates relevant protein evolution. This study test hypothesizes that proteins subject to sequence evolution structural constraints play a crucial role and that these constraints hinder the modification of such proteins in this process. These hypotheses were tested using Plasmodium falciparum model and evaluated protein structures predicted for the entire proteome by the AlphaFold method. Based on dN/dS test results and P. falciparum and P. reichenowi comparisons, the presented approach identified proteins subject to purifying selection acting on the whole sequence and buried residues (dN < dS) and positive selection on nonburied residues. Of the 26 proteins, some known antigens (ring-exported protein 3, RAP protein, erythrocyte binding antigen-140, and protein P47) targeted by the host immune system are promising vaccine candidates. The set also contained 11 enzymes, including FIKK kinase, which modifies host proteins. This set was compared with genes for which the dN/dS test suggested that positive selection acts on the whole gene (i.e., dN > dS). The present study found that such genes encode enzymes and antigenic vaccine candidates less frequently than genes for which evolution is not subject to selection constraints and positive selection acts on only exposed residues. The analysis was repeated comparing P. falciparum with P. alderi, which is more distantly related. The study discusses the potential implications of the presented methodology for rational vaccine design and the parasitology and evolutionary biology fields.

Repurposing of Plasmodium falciparum var genes beyond the blood stage

A commonly observed survival strategy in protozoan parasites is the sequential expression of clonally variant-surface antigens to avoid elimination by the host's immune response. In malaria-causing P. falciparum, the immunovariant erythrocyte-membrane protein-1 (PfEMP1) adhesin family, encoded by var genes, is responsible for both antigenic variation and cytoadherence of infected erythrocytes to the microvasculature. Until recently, the biological function of these variant genes was believed to be restricted to intraerythrocytic developmental stages. With the advent of new technologies, var gene expression has been confirmed in transmission and pre-erythrocytic stages. Here, we discuss how repurposing of var gene expression beyond chronic blood-stage infection may be critical for successful transmission.

Discovery of RUF6 ncRNA-interacting proteins involved in P. falciparum immune evasion

Non-coding RNAs (ncRNAs) are emerging regulators of immune evasion and transmission of Plasmodium falciparum RUF6 is an ncRNA gene family that is transcribed by RNA polymerase III but actively regulates the Pol II-transcribed var virulence gene family. Understanding how RUF6 ncRNA connects to downstream effectors is lacking. We developed an RNA-directed proteomic discovery (ChIRP-MS) protocol to identify in vivo RUF6 ncRNA-protein interactions. The RUF6 ncRNA interactome was purified with biotinylated antisense oligonucleotides. Quantitative label-free mass spectrometry identified several unique proteins linked to gene transcription including RNA Pol II subunits, nucleosome assembly proteins, and a homologue of DEAD box helicase 5 (DDX5). Affinity purification of Pf-DDX5 identified proteins originally found by our RUF6-ChIRP protocol, validating the technique's robustness for identifying ncRNA interactomes in P. falciparum Inducible displacement of nuclear Pf-DDX5 resulted in significant down-regulation of the active var gene. Our work identifies a RUF6 ncRNA-protein complex that interacts with RNA Pol II to sustain the var gene expression, including a helicase that may resolve G-quadruplex secondary structures in var genes to facilitate transcriptional activation and progression.

Epigenetic and Epitranscriptomic Gene Regulation in Plasmodium falciparum and How We Can Use It against Malaria

Malaria, caused by Plasmodium parasites, is still one of the biggest global health challenges. P. falciparum is the deadliest species to humans. In this review, we discuss how this parasite develops and adapts to the complex and heterogenous environments of its two hosts thanks to varied chromatin-associated and epigenetic mechanisms. First, one small family of transcription factors, the ApiAP2 proteins, functions as master regulators of spatio-temporal patterns of gene expression through the parasite life cycle. In addition, chromatin plasticity determines variable parasite cell phenotypes that link to parasite growth, virulence and transmission, enabling parasite adaptation within host conditions. In recent years, epitranscriptomics is emerging as a new regulatory layer of gene expression. We present evidence of the variety of tRNA and mRNA modifications that are being characterized in Plasmodium spp., and the dynamic changes in their abundance during parasite development and cell fate. We end up outlining that new biological systems, like the mosquito model, to decipher the unknowns about epigenetic mechanisms in vivo; and novel methodologies, to study the function of RNA modifications; are needed to discover the Achilles heel of the parasite. With this new knowledge, future strategies manipulating the epigenetics and epitranscriptomic machinery of the parasite have the potential of providing new weapons against malaria.

Emerging biology of noncoding RNAs in malaria parasites

In eukaryotic organisms, noncoding RNAs (ncRNAs) have been implicated as important regulators of multifaceted biological processes, including transcriptional, posttranscriptional, and epigenetic regulation of gene expression. In recent years, it is becoming clear that protozoan parasites encode diverse ncRNA transcripts; however, little is known about their cellular functions. Recent advances in high-throughput “omic” studies identified many novel long ncRNAs (lncRNAs) in apicomplexan parasites, some of which undergo splicing, polyadenylation, and encode small proteins. To date, only a few of them are characterized, leaving a big gap in our understanding regarding their origin, mode of action, and functions in parasite biology. In this review, we focus on lncRNAs of the human malaria parasite Plasmodium falciparum and highlight their cellular functions and possible mechanisms of action.

Mechanisms for Host Immune Evasion Mediated by Plasmodium falciparum-Infected Erythrocyte Surface Antigens

Plasmodium falciparum infection causes the most severe form of malaria. It has been hypothesized that P. falciparum directly suppresses host immune responses because sufficient acquired immunity is often not induced even by repeated P. falciparum infections in malaria-endemic areas. It is known that many kinds of P. falciparum-derived proteins are expressed on the surface of P. falciparum-infected erythrocytes (IEs), and these proteins have long been thought to be a key to the elucidation of the host immune evasion mechanisms. Our recent studies have revealed that the P. falciparum-derived erythrocyte surface antigen, RIFIN, the largest multiple gene family protein in the P. falciparum genome, suppresses host immune cell activation through direct interaction with human inhibitory immune receptors. In this review, we will discuss the molecular mechanisms for host immune evasion by P. falciparum-infected erythrocyte surface antigens. In addition, we will discuss the recently identified host immune response to P. falciparum using specialized antibodies that target host-P. falciparum-derived molecule interactions.

P. falciparum PfRUVBL proteins binds at TARE region and var gene promoter located in subtelomeric region

In order to survive and establish infection, Plasmodium parasite employ various strategies to evade host immune response. Var genes family, a repertoire of 60 genes, express parasite-specific protein PfEMP1, a variable surface antigen, on the membrane of infected erythrocyte, and by continuously switching the variants of PfEMP1, helps the parasite to avoid detection and destruction by host immune system during intra-erythrocytic developmental cycle. Although chromatin modifications are recognized to be a prominent phenomenon in regulation of mono-allelic expression of these var genes, however the precise histone codes and molecular players & mechanisms guiding these modifications are yet to be unravelled in depth. In this study, we have functionally characterized RUVBL proteins of P. falciparum and shown that PfMYST (an essential lysine acetyl transferase) and PfRUVBL protein complex occupy the TARE region and var gene promoter in ring stage of the parasite. Further we have demonstrated that PfMYST/PfRUVBL complex interact with core histone, H3 & H4. Overall the findings of this study adds a layer by identifying the potential role of epigenetic regulators, PfMYST & PfRUVBL in regulation of monoallelic expression of var genes in malaria parasite.

Plasmodium falciparum Acetyl-CoA Synthetase Is Essential for Parasite Intraerythrocytic Development and Chromatin Modification

The malaria parasite Plasmodium falciparum possesses a unique Acetyl-CoA Synthetase (PfACS), which provides acetyl moieties for different metabolic and regulatory cellular pathways. We characterized PfACS and studied its role focusing on epigenetic modifications using the var gene family as reporter genes. For this, mutant lines to modulate plasmodial ACS expression by degron-mediated protein degradation and ribozyme-induced transcript decay were created. Additionally, an inhibitor of the human Acetyl-CoA Synthetase 2 was tested for its effectiveness in interfering with PfACS. The knockdown of PfACS or its inhibition resulted in impaired parasite growth. Decreased levels of PfACS also led to differential histone acetylation patterns, altered variant gene expression, and concomitantly decreased cytoadherence of infected red blood cells containing knocked-down parasites. Further, ChIP analysis revealed the presence of PfACS in many loci in ring stage parasites, underscoring its involvement in the regulation of chromatin. Due to its central function in the plasmodial metabolism and significant differences to human ACS, PfACS is an interesting target for drug development.

Expression Patterns of Plasmodium falciparum Clonally Variant Genes at the Onset of a Blood Infection in Malaria-Naive Humans

Clonally variant genes (CVGs) play fundamental roles in the adaptation of Plasmodium falciparum to fluctuating conditions of the human host. However, their expression patterns under the natural conditions of the blood circulation have been characterized in detail for only a few specific gene families. Here, we provide a detailed characterization of the complete P. falciparum transcriptome across the full intraerythrocytic development cycle (IDC) at the onset of a blood infection in malaria-naive human volunteers. We found that the vast majority of transcriptional differences between parasites obtained from the volunteers and the parental parasite line maintained in culture occurred in CVGs. In particular, we observed a major increase in the transcript levels of most genes of the pfmc-2tm and gbp families and of specific genes of other families, such as phist, hyp10, rif, or stevor, in addition to previously reported changes in var and clag3 gene expression. Increased transcript levels of individual pfmc-2tm, rif, and stevor genes involved activation in small subsets of parasites. Large transcriptional differences correlated with changes in the distribution of heterochromatin, confirming their epigenetic nature. Furthermore, the similar expression of several CVGs between parasites collected at different time points along the blood infection suggests that the epigenetic memory for multiple CVG families is lost during transmission stages, resulting in a reset of their transcriptional state. Finally, the CVG expression patterns observed in a volunteer likely infected by a single sporozoite suggest that new epigenetic patterns are established during liver stages.

Structural basis of LAIR1 targeting by polymorphic Plasmodium RIFINs

RIFIN, a large family of Plasmodium variant surface antigens, plays a crucial role in malaria pathogenesis by mediating immune suppression through activation of inhibitory receptors such as LAIR1, and antibodies with LAIR1 inserts have been identified that bind infected erythrocytes through RIFIN. However, details of RIFIN-mediated LAIR1 recognition and receptor activation have been unclear. Here, we use negative-stain EM to define the architecture of LAIR1-inserted antibodies and determine crystal structures of RIFIN-variable 2 (V2) domain in complex with a LAIR1 domain. These structures reveal the LAIR1-binding region of RIFIN to be hydrophobic and membrane-distal, to exhibit extensive structural diversity, and to interact with RIFIN-V2 in a one-to-one fashion. Through structural and sequence analysis of various LAIR1 constructs, we identify essential elements of RIFIN-binding on LAIR1. Furthermore, a structure-derived LAIR1-binding sequence signature ascertained >20 LAIR1-binding RIFINs, including some from P. falciparum field strains and Plasmodium species infecting gorillas and chimpanzees.

A single-cell atlas of Plasmodium falciparum transmission through the mosquito

Malaria parasites have a complex life cycle featuring diverse developmental strategies, each uniquely adapted to navigate specific host environments. Here we use single-cell transcriptomics to illuminate gene usage across the transmission cycle of the most virulent agent of human malaria - Plasmodium falciparum. We reveal developmental trajectories associated with the colonization of the mosquito midgut and salivary glands and elucidate the transcriptional signatures of each transmissible stage. Additionally, we identify both conserved and non-conserved gene usage between human and rodent parasites, which point to both essential mechanisms in malaria transmission and species-specific adaptations potentially linked to host tropism. Together, the data presented here, which are made freely available via an interactive website, provide a fine-grained atlas that enables intensive investigation of the P. falciparum transcriptional journey. As well as providing insights into gene function across the transmission cycle, the atlas opens the door for identification of drug and vaccine targets to stop malaria transmission and thereby prevent disease.

Here the authors use single-cell RNA-seq to profile the transmission stages of the human malaria parasite Plasmodium falciparum as it progresses through the Anopheles mosquito. They highlight unique patterns of gene usage throughout this development and identify potential pleiotropic genes that function at multiple life cycle stages.

Unlocking G-Quadruplexes as Antiviral Targets

Guanine-rich DNA and RNA sequences can fold into noncanonical nucleic acid structures called G-quadruplexes (G4s). Since the discovery that these structures may act as scaffolds for the binding of specific ligands, G4s aroused the attention of a growing number of scientists. The versatile roles of G4 structures in viral replication, transcription, and translation suggest direct applications in therapy or diagnostics. G4-interacting molecules (proteins or small molecules) may also affect the balance between latent and lytic phases, and increasing evidence reveals that G4s are implicated in generally suppressing viral processes, such as replication, transcription, translation, or reverse transcription. In this review, we focus on the discovery of G4s in viruses and the role of G4 ligands in the antiviral drug discovery process. After assessing the role of viral G4s, we argue that host G4s participate in immune modulation, viral tumorigenesis, cellular pathways involved in virus maturation, and DNA integration of viral genomes, which can be potentially employed for antiviral therapeutics. Furthermore, we scrutinize the impediments and shortcomings in the process of studying G4 ligands and drug discovery. Finally, some unanswered questions regarding viral G4s are highlighted for prospective future projects. SIGNIFICANCE STATEMENT: G-quadruplexes (G4s) are noncanonical nucleic acid structures that have gained increasing recognition during the last few decades. First identified as relevant targets in oncology, their importance in virology is now increasingly clear. A number of G-quadruplex ligands are known: viral transcription and replication are the main targets of these ligands. Both viral and cellular G4s may be targeted; this review embraces the different aspects of G-quadruplexes in both host and viral contexts.

Trypanosoma brucei ATR Links DNA Damage Signaling during Antigenic Variation with Regulation of RNA Polymerase I-Transcribed Surface Antigens

Trypanosoma brucei evades mammalian immunity by using recombination to switch its surface-expressed variant surface glycoprotein (VSG), while ensuring that only one of many subtelomeric multigene VSG expression sites are transcribed at a time. DNA repair activities have been implicated in the catalysis of VSG switching by recombination, not transcriptional control. How VSG switching is signaled to guide the appropriate reaction or to integrate switching into parasite growth is unknown. Here, we show that the loss of ATR, a DNA damage-signaling protein kinase, is lethal, causing nuclear genome instability and increased VSG switching through VSG-localized damage. Furthermore, ATR loss leads to the increased transcription of silent VSG expression sites and expression of mixed VSGs on the cell surface, effects that are associated with the altered localization of RNA polymerase I and VEX1. This work shows that ATR acts in antigenic variation both through DNA damage signaling and surface antigen expression control.

•

Loss of the repair protein kinase ATR in Trypanosoma brucei is lethal

•

Loss of T. brucei ATR alters VSG coat expression needed for immune evasion

•

Monoallelic RNA polymerase I VSG expression is undermined by ATR loss

•

ATR loss leads to expression of subtelomeric VSGs, indicative of recombination

Mapping immune variation and var gene switching in naive hosts infected with Plasmodium falciparum

Falciparum malaria is clinically heterogeneous and the relative contribution of parasite and host in shaping disease severity remains unclear. We explored the interaction between inflammation and parasite variant surface antigen (VSA) expression, asking whether this relationship underpins the variation observed in controlled human malaria infection (CHMI). We uncovered marked heterogeneity in the host response to blood challenge; some volunteers remained quiescent, others triggered interferon-stimulated inflammation and some showed transcriptional evidence of myeloid cell suppression. Significantly, only inflammatory volunteers experienced hallmark symptoms of malaria. When we tracked temporal changes in parasite VSA expression to ask whether variants associated with severe disease rapidly expand in naive hosts, we found no transcriptional evidence to support this hypothesis. These data indicate that parasite variants that dominate severe malaria do not have an intrinsic growth or survival advantage; instead, they presumably rely upon infection-induced changes in their within-host environment for selection.

Investigation of Heterochromatin Protein 1 Function in the Malaria Parasite Plasmodium falciparum Using a Conditional Domain Deletion and Swapping Approach

The human malaria parasite Plasmodium falciparum encodes a single ortholog of heterochromatin protein 1 (PfHP1) that plays a crucial role in the epigenetic regulation of various survival-related processes. PfHP1 is essential for parasite proliferation and the heritable silencing of genes linked to antigenic variation, host cell invasion, and sexual conversion. Here, we employed CRISPR/Cas9-mediated genome editing combined with the DiCre/loxP system to investigate how the PfHP1 chromodomain (CD), hinge domain, and chromoshadow domain (CSD) contribute to overall PfHP1 function. We show that the 76 C-terminal residues are responsible for targeting PfHP1 to the nucleus. Furthermore, we reveal that each of the three functional domains of PfHP1 are required for heterochromatin formation, gene silencing, and mitotic parasite proliferation. Finally, we discovered that the hinge domain and CSD of HP1 are functionally conserved between P. falciparum and P. berghei, a related malaria parasite infecting rodents. In summary, our study provides new insights into PfHP1 function and offers a tool for further studies on epigenetic regulation and life cycle decision in malaria parasites.IMPORTANCE Malaria is caused by unicellular Plasmodium species parasites that repeatedly invade and replicate inside red blood cells. Some blood-stage parasites exit the cell cycle and differentiate into gametocytes that are essential for malaria transmission via the mosquito vector. Epigenetic control mechanisms allow the parasites to alter the expression of surface antigens and to balance the switch between parasite multiplication and gametocyte production. These processes are crucial to establish chronic infection and optimize parasite transmission. Here, we performed a mutational analysis of heterochromatin protein 1 (HP1) in P. falciparum We demonstrate that all three domains of this protein are indispensable for the proper function of HP1 in parasite multiplication, heterochromatin formation, and gene silencing. Moreover, expression of chimeric proteins revealed the functional conservation of HP1 proteins between different Plasmodium species. These results provide new insight into the function and evolution of HP1 as an essential epigenetic regulator of parasite survival.

Histone acetyltransferase PfGCN5 regulates stress responsive and artemisinin resistance related genes in Plasmodium falciparum

Plasmodium falciparum has evolved resistance to almost all front-line drugs including artemisinin, which threatens malaria control and elimination strategies. Oxidative stress and protein damage responses have emerged as key players in the generation of artemisinin resistance. In this study, we show that PfGCN5, a histone acetyltransferase, binds to the stress-responsive genes in a poised state and regulates their expression under stress conditions. Furthermore, we show that upon artemisinin exposure, genome-wide binding sites for PfGCN5 are increased and it is directly associated with the genes implicated in artemisinin resistance generation like BiP and TRiC chaperone. Interestingly, expression of genes bound by PfGCN5 was found to be upregulated during stress conditions. Moreover, inhibition of PfGCN5 in artemisinin-resistant parasites increases the sensitivity of the parasites to artemisinin treatment indicating its role in drug resistance generation. Together, these findings elucidate the role of PfGCN5 as a global chromatin regulator of stress-responses with a potential role in modulating artemisinin drug resistance and identify PfGCN5 as an important target against artemisinin-resistant parasites.

Exploring the virulence gene interactome with CRISPR/dCas9 in the human malaria parasite

Mutually exclusive expression of the var multigene family is key to immune evasion and pathogenesis in Plasmodium falciparum, but few factors have been shown to play a direct role. We adapted a CRISPR-based proteomics approach to identify novel factors associated with var genes in their natural chromatin context. Catalytically inactive Cas9 ("dCas9") was targeted to var gene regulatory elements, immunoprecipitated, and analyzed with mass spectrometry. Known and novel factors were enriched including structural proteins, DNA helicases, and chromatin remodelers. Functional characterization of PfISWI, an evolutionarily divergent putative chromatin remodeler enriched at the var gene promoter, revealed a role in transcriptional activation. Proteomics of PfISWI identified several proteins enriched at the var gene promoter such as acetyl-CoA synthetase, a putative MORC protein, and an ApiAP2 transcription factor. These findings validate the CRISPR/dCas9 proteomics method and define a new var gene-associated chromatin complex. This study establishes a tool for targeted chromatin purification of unaltered genomic loci and identifies novel chromatin-associated factors potentially involved in transcriptional control and/or chromatin organization of virulence genes in the human malaria parasite.

The second life of Plasmodium in the mosquito host: gene regulation on the move

Malaria parasites face dynamically changing environments and strong selective constraints within human and mosquito hosts. To survive such hostile and shifting conditions, Plasmodium switches transcriptional programs during development and has evolved mechanisms to adjust its phenotype through heterogeneous patterns of gene expression. In vitro studies on culture-adapted isolates have served to set the link between chromatin structure and functional gene expression. Yet, experimental evidence is limited to certain stages of the parasite in the vertebrate, i.e. blood, while the precise mechanisms underlying the dynamic regulatory landscapes during development and in the adaptation to within-host conditions remain poorly understood. In this review, we discuss available data on transcriptional and epigenetic regulation in Plasmodium mosquito stages in the context of sporogonic development and phenotypic variation, including both bet-hedging and environmentally triggered direct transcriptional responses. With this, we advocate the mosquito offers an in vivo biological model to investigate the regulatory networks, transcription factors and chromatin-modifying enzymes and their modes of interaction with regulatory sequences, which might be responsible for the plasticity of the Plasmodium genome that dictates stage- and cell type-specific blueprints of gene expression.

Transcriptional variation in malaria parasites: why and how

Transcriptional differences enable the generation of alternative phenotypes from the same genome. In malaria parasites, transcriptional plasticity plays a major role in the process of adaptation to fluctuations in the environment. Multiple studies with culture-adapted parasites and field isolates are starting to unravel the different transcriptional alternatives available to Plasmodium falciparum and the underlying molecular mechanisms. Here we discuss how epigenetic variation, directed transcriptional responses and also genetic changes that affect transcript levels can all contribute to transcriptional variation and, ultimately, parasite survival. Some transcriptional changes are driven by stochastic events. These changes can occur spontaneously, resulting in heterogeneity within parasite populations that provides the grounds for adaptation by dynamic natural selection. However, transcriptional changes can also occur in response to external cues. A better understanding of the mechanisms that the parasite has evolved to alter its transcriptome may ultimately contribute to the design of strategies to combat malaria to which the parasite cannot adapt.

Complex nutrient channel phenotypes despite Mendelian inheritance in a Plasmodium falciparum genetic cross

Malaria parasites activate a broad-selectivity ion channel on their host erythrocyte membrane to obtain essential nutrients from the bloodstream. This conserved channel, known as the plasmodial surface anion channel (PSAC), has been linked to parasite clag3 genes in P. falciparum, but epigenetic switching between the two copies of this gene hinders clear understanding of how the encoded protein determines PSAC activity. Here, we used linkage analysis in a P. falciparum cross where one parent carries a single clag3 gene to overcome the effects of switching and confirm a primary role of the clag3 product with high confidence. Despite Mendelian inheritance, CLAG3 conditional knockdown revealed remarkably preserved nutrient and solute uptake. Even more surprisingly, transport remained sensitive to a CLAG3 isoform-specific inhibitor despite quantitative knockdown, indicating that low doses of the CLAG3 transgene are sufficient to confer block. We then produced a complete CLAG3 knockout line and found it exhibits an incomplete loss of transport activity, in contrast to rhoph2 and rhoph3, two PSAC-associated genes that cannot be disrupted because nutrient uptake is abolished in their absence. Although the CLAG3 knockout did not incur a fitness cost under standard nutrient-rich culture conditions, this parasite could not be propagated in a modified medium that more closely resembles human plasma. These studies implicate oligomerization of CLAG paralogs encoded by various chromosomes in channel formation. They also reveal that CLAG3 is dispensable under standard in vitro conditions but required for propagation under physiological conditions.

CRISPR Interference of a Clonally Variant GC-Rich Noncoding RNA Family Leads to General Repression of var Genes in Plasmodium falciparum

Plasmodium falciparum is the deadliest malaria parasite species, accounting for the vast majority of disease cases and deaths. The virulence of this parasite is reliant upon the mutually exclusive expression of cytoadherence proteins encoded by the 60-member var gene family. Antigenic variation of this multigene family serves as an immune evasion mechanism, ultimately leading to chronic infection and pathogenesis. Understanding the regulation mechanism of antigenic variation is key to developing new therapeutic and control strategies. Our study uncovers a novel layer in the epigenetic regulation of transcription of this family of virulence genes by means of a multigene-targeting CRISPR interference approach.

The human malaria parasite Plasmodium falciparum uses mutually exclusive expression of the PfEMP1-encoding var gene family to evade the host immune system. Despite progress in the molecular understanding of the default silencing mechanism, the activation mechanism of the uniquely expressed var member remains elusive. A GC-rich noncoding RNA (ncRNA) gene family has coevolved with Plasmodium species that express var genes. Here, we show that this ncRNA family is transcribed in a clonally variant manner, with predominant transcription of a single member occurring when the ncRNA is located adjacent to and upstream of an active var gene. We developed a specific CRISPR interference (CRISPRi) strategy that allowed for the transcriptional repression of all GC-rich members. A lack of GC-rich ncRNA transcription led to the downregulation of the entire var gene family in ring-stage parasites. Strikingly, in mature blood-stage parasites, the GC-rich ncRNA CRISPRi affected the transcription patterns of other clonally variant gene families, including the downregulation of all Pfmc-2TM members. We provide evidence for the key role of GC-rich ncRNA transcription in var gene activation and discovered a molecular link between the transcriptional control of various clonally variant multigene families involved in parasite virulence. This work opens new avenues for elucidating the molecular processes that control immune evasion and pathogenesis in P. falciparum.

Plasmodium falciparum-CD36 Structure-Function Relationships Defined by Ortholog Scanning Mutagenesis

BACKGROUND

The interaction of Plasmodium falciparum-infected erythrocytes (IEs) with the host receptor CD36 is among the most studied host-parasite interfaces. CD36 is a scavenger receptor that binds numerous ligands including the cysteine-rich interdomain region (CIDR)α domains of the erythrocyte membrane protein 1 family (PfEMP1) expressed on the surface of IEs. CD36 is conserved across species, but orthologs display differential binding of IEs.

METHODS

In this study, we exploited these differences, combined with the recent crystal structure and 3-dimensional modeling of CD36, to investigate malaria-CD36 structure-function relationships and further define IE-CD36 binding interactions.

RESULTS

We show that a charged surface in the membrane-distal region of CD36 is necessary for IE binding. Moreover, IE interaction with this binding surface is influenced by additional CD36 domains, both proximal to and at a distance from this site.

CONCLUSIONS

Our data indicate that subtle sequence and spatial differences in these domains modify receptor conformation and regulate the ability of CD36 to selectively interact with its diverse ligands.

Mapping and functional analysis of heterochromatin protein 1 phosphorylation in the malaria parasite Plasmodium falciparum

Previous studies in model eukaryotes have demonstrated that phosphorylation of heterochromatin protein 1 (HP1) is important for dynamically regulating its various functions. However, in the malaria parasite Plasmodium falciparum both the function of HP1 phosphorylation and the identity of the protein kinases targeting HP1 are still elusive. In order to functionally analyze phosphorylation of P. falciparum HP1 (PfHP1), we first mapped PfHP1 phosphorylation sites by liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis of native PfHP1, which identified motifs from which potential kinases could be predicted; in particular, several phosphorylated residues were embedded in motifs rich in acidic residues, reminiscent of targets for P. falciparum casein kinase 2 (PfCK2). Secondly, we tested recombinant PfCK2 and a number of additional protein kinases for their ability to phosphorylate PfHP1 in in vitro kinase assays. These experiments validated our prediction that PfHP1 acts as a substrate for PfCK2. Furthermore, LC-MS/MS analysis showed that PfCK2 phosphorylates three clustered serine residues in an acidic motif within the central hinge region of PfHP1. To study the role of PfHP1 phosphorylation in live parasites we used CRISPR/Cas9-mediated genome editing to generate a number of conditional PfHP1 phosphomutants based on the DiCre/LoxP system. Our studies revealed that neither PfCK2-dependent phosphorylation of PfHP1, nor phosphorylation of the hinge domain in general, affect PfHP1's ability to localize to heterochromatin, and that PfHP1 phosphorylation in this region is dispensable for the proliferation of P. falciparum blood stage parasites.

Structural insights into diverse modes of ICAM-1 binding by Plasmodium falciparum-infected erythrocytes

A major determinant of pathogenicity in malaria caused by Plasmodium falciparum is the adhesion of parasite-infected erythrocytes to the vasculature or tissues of infected individuals. This occludes blood flow, leads to inflammation, and increases parasitemia by reducing spleen-mediated clearance of the parasite. This adhesion is mediated by PfEMP1, a multivariant family of around 60 proteins per parasite genome which interact with specific host receptors. One of the most common of these receptors is intracellular adhesion molecule-1 (ICAM-1), which is bound by 2 distinct groups of PfEMP1, A-type and B or C (BC)-type. Here, we present the structure of a domain from a B-type PfEMP1 bound to ICAM-1, revealing a complex binding site. Comparison with the existing structure of an A-type PfEMP1 bound to ICAM-1 shows that the 2 complexes share a globally similar architecture. However, while the A-type PfEMP1 bind ICAM-1 through a highly conserved binding surface, the BC-type PfEMP1 use a binding site that is more diverse in sequence, similar to how PfEMP1 interact with other human receptors. We also show that A- and BC-type PfEMP1 present ICAM-1 at different angles, perhaps influencing the ability of neighboring PfEMP1 domains to bind additional receptors. This illustrates the deep diversity of the PfEMP1 and demonstrates how variations in a single domain architecture can modulate binding to a specific ligand to control function and facilitate immune evasion.

The molecular machinery of translational control in malaria parasites

Translational control regulates the levels of protein synthesized from its transcript and is key for the rapid adjustment of gene expression in response to environmental stimuli. The regulation of translation is of special importance for malaria parasites, which pass through a complex life cycle that includes various replication phases in the different organs of the human and mosquito hosts and a sexual reproduction phase in the mosquito midgut. In particular, the quiescent transmission stages rely on translational control to rapidly adapt to the new environment, once they switch over from the human to the mosquito and vice versa. Three control mechanisms are currently proposed in Plasmodium, (1) global regulation that acts on the translation initiation complex; (2) mRNA-specific regulation, involving cis control elements, mRNA-binding proteins and translational repressors; and (3) induced mRNA decay by the Ccr4-Not and the RNA exosome complex. The main molecules controlling translation are highly conserved in malaria parasites and an increasing number of studies shed light on the interwoven pathways leading to the up or downregulation of protein synthesis in the diverse plasmodial stages. We here highlight recent findings on translational control during life cycle progression of Plasmodium and discuss the molecules involved in regulating protein synthesis.

A Specific PfEMP1 Is Expressed in P. falciparum Sporozoites and Plays a Role in Hepatocyte Infection

Heterochromatin plays a central role in the process of immune evasion, pathogenesis, and transmission of the malaria parasite Plasmodium falciparum during blood stage infection. Here, we use ChIP sequencing to demonstrate that sporozoites from mosquito salivary glands expand heterochromatin at subtelomeric regions to silence blood-stage-specific genes. Our data also revealed that heterochromatin enrichment is predictive of the transcription status of clonally variant genes members that mediate cytoadhesion in blood stage parasites. A specific member (here called NF54var^sporo) of the var gene family remains euchromatic, and the resultant PfEMP1 (NF54_SpzPfEMP1) is expressed at the sporozoite surface. NF54_SpzPfEMP1-specific antibodies efficiently block hepatocyte infection in a strain-specific manner. Furthermore, human volunteers immunized with infective sporozoites developed antibodies against NF54_SpzPfEMP1. Overall, we show that the epigenetic signature of var genes is reset in mosquito stages. Moreover, the identification of a strain-specific sporozoite PfEMP1 is highly relevant for vaccine design based on sporozoites.

•

Sporozoites expand subtelomeric heterochromatin to silence blood-stage-specific genes

•

A strain-specific PfEMP1 is expressed on the surface of sporozoites

•

NF54_SpzPfEMP1 is immunogenic in sporozoite-infected human volunteers

•

Antibodies against NF54_SpzPfEMP1 block sporozoite infection of hepatocytes

The Malaria Cell Atlas: Single parasite transcriptomes across the complete Plasmodium life cycle

Malaria parasites adopt a remarkable variety of morphological life stages as they transition through multiple mammalian host and mosquito vector environments. We profiled the single-cell transcriptomes of thousands of individual parasites, deriving the first high-resolution transcriptional atlas of the entire Plasmodium berghei life cycle. We then used our atlas to precisely define developmental stages of single cells from three different human malaria parasite species, including parasites isolated directly from infected individuals. The Malaria Cell Atlas provides both a comprehensive view of gene usage in a eukaryotic parasite and an open-access reference dataset for the study of malaria parasites.

Monoallelic expression and epigenetic inheritance sustained by a Trypanosoma brucei variant surface glycoprotein exclusion complex

The largest gene families in eukaryotes are subject to allelic exclusion, but mechanisms underpinning single allele selection and inheritance remain unclear. Here, we describe a protein complex sustaining variant surface glycoprotein (VSG) allelic exclusion and antigenic variation in Trypanosoma brucei parasites. The VSG-exclusion-1 (VEX1) protein binds both telomeric VSG-associated chromatin and VEX2, an ortholog of nonsense-mediated-decay helicase, UPF1. VEX1 and VEX2 assemble in an RNA polymerase-I transcription-dependent manner and sustain the active, subtelomeric VSG-associated transcription compartment. VSG transcripts and VSG coats become highly heterogeneous when VEX proteins are depleted. Further, the DNA replication-associated chromatin assembly factor, CAF-1, binds to and specifically maintains VEX1 compartmentalisation following DNA replication. Thus, the VEX-complex controls VSG-exclusion, while CAF-1 sustains VEX-complex inheritance in association with the active-VSG. Notably, the VEX2-orthologue and CAF-1 in mammals are also implicated in exclusion and inheritance functions. In trypanosomes, these factors sustain a highly effective and paradigmatic immune evasion strategy.

Phosphorylation of the VAR2CSA extracellular region is associated with enhanced adhesive properties to the placental receptor CSA

Plasmodium falciparum is the main cause of disease and death from malaria. P. falciparum virulence resides in the ability of infected erythrocytes (IEs) to sequester in various tissues through the interaction between members of the polymorphic P. falciparum erythrocyte membrane protein 1 (PfEMP1) adhesin family to various host receptors. Here, we investigated the effect of phosphorylation of variant surface antigen 2-CSA (VAR2CSA), a member of the PfEMP1 family associated to placental sequestration, on its capacity to adhere to chondroitin sulfate A (CSA) present on the placental syncytium. We showed that phosphatase treatment of IEs impairs cytoadhesion to CSA. MS analysis of recombinant VAR2CSA phosphosites prior to and after phosphatase treatment, as well as of native VAR2CSA expressed on IEs, identified critical phosphoresidues associated with CSA binding. Site-directed mutagenesis on recombinant VAR2CSA of 3 phosphoresidues localised within the CSA-binding region confirmed in vitro their functional importance. Furthermore, using clustered regularly interspaced short palindromic repeats/CRISPR-associated protein-9 nuclease (CRISPR/Cas9), we generated a parasite line in which the phosphoresidue T934 is changed to alanine and showed that this mutation strongly impairs IEs cytoadhesion to CSA. Taken together, these results demonstrate that phosphorylation of the extracellular region of VAR2CSA plays a major role in IEs cytoadhesion to CSA and provide new molecular insights for strategies aiming to reduce the morbidity and mortality of PM.

Host-parasite interactions and ecology of the malaria parasite-a bioinformatics approach

Malaria remains one of the highest mortality infectious diseases. Malaria is caused by parasites from the genus Plasmodium. Most deaths are caused by infections involving Plasmodium falciparum, which has a complex life cycle. Malaria parasites are extremely well adapted for interactions with their host and their host's immune system and are able to suppress the human immune system, erase immunological memory and rapidly alter exposed antigens. Owing to this rapid evolution, parasites develop drug resistance and express novel forms of antigenic proteins that are not recognized by the host immune system. There is an emerging need for novel interventions, including novel drugs and vaccines. Designing novel therapies requires knowledge about host-parasite interactions, which is still limited. However, significant progress has recently been achieved in this field through the application of bioinformatics analysis of parasite genome sequences. In this review, we describe the main achievements in 'malarial' bioinformatics and provide examples of successful applications of protein sequence analysis. These examples include the prediction of protein functions based on homology and the prediction of protein surface localization via domain and motif analysis. Additionally, we describe PlasmoDB, a database that stores accumulated experimental data. This tool allows data mining of the stored information and will play an important role in the development of malaria science. Finally, we illustrate the application of bioinformatics in the development of population genetics research on malaria parasites, an approach referred to as reverse ecology.

Host-parasite interaction explains variation in the prevalence of avian haemosporidians at the community level

Parasites are a selective force that shape host community structure and dynamics, but host communities can also influence parasitism. Understanding the dual nature from host-parasite interactions can be facilitated by quantifying the variation in parasite prevalence among host species and then comparing that variation to other ecological factors that are known to also shape host communities. Avian haemosporidian parasites (e.g. Plasmodium and Haemoproteus) are abundant and widespread representing an excellent model for the study of host-parasite interactions. Several geographic and environmental factors have been suggested to determine prevalence of avian haemosporidians in bird communities. However, it remains unknown whether host and parasite traits, represented by phylogenetic distances among species and degree of specialization in host-parasite relationships, can influence infection status. The aims of this study were to analyze factors affecting infection status in a bird community and to test whether the degree of parasite specialization on their hosts is determined by host traits. Our statistical analyses suggest that infection status is mainly determined by the interaction between host species and parasite lineages where tolerance and/or susceptibility to parasites plays an essential role. Additionally, we found that although some of the parasite lineages infected a low number of bird individuals, the species they infected were distantly related and therefore the parasites themselves should not be considered typical host specialists. Infection status was higher for generalist than for specialist parasites in some, but not all, host species. These results suggest that detected prevalence in a species mainly results from the interaction between host immune defences and parasite exploitation strategies wherein the result of an association between particular parasite lineages and particular host species is idiosyncratic.

Exploring parasite heterogeneity using single-cell RNA-seq reveals a gene signature among sexual stage Plasmodium falciparum parasites

The malaria parasite has a complex lifecycle, including several events of differentiation and stage progression, while actively evading immunity in both its mosquito and human hosts. Important parasite gene expression and regulation during these events remain hidden in rare populations of cells. Here, we combine a capillary-based platform for cell isolation with single-cell RNA-sequencing to transcriptionally profile 165 single infected red blood cells (iRBCs) during the intra-erythrocytic developmental cycle (IDC). Unbiased analyses of single-cell data grouped the cells into eight transcriptional states during IDC. Interestingly, we uncovered a gene signature from the single iRBC analyses that can successfully discriminate between developing asexual and sexual stage parasites at cellular resolution, and we verify five, previously undefined, gametocyte stage specific genes. Moreover, we show the capacity of detecting expressed genes from the variable gene families in single parasites, despite the sparse nature of data. In total, the single parasite transcriptomics holds promise for molecular dissection of rare parasite phenotypes throughout the malaria lifecycle.

Independent regulation of Plasmodium falciparum rif gene promoters

All Plasmodium species express variant antigens which may mediate immune escape in the vertebrate host. In Plasmodium falciparum, the rif gene family encodes variant antigens which are partly exposed on the infected red blood cell surface and may function as virulence factors. Not all rif genes are expressed at the same time and it is unclear what controls rif gene expression. In this work, we addressed global rif transcription using plasmid vectors with two drug resistance markers, one controlled by a rif 5′ upstream region and the second by a constitutively active promoter. After spontaneous integration into the genome of one construct, we observed that the resistance marker controlled by the rif 5′ upstream region was expressed dependent on the applied drug pressure. Then, the global transcription of rif genes in these transfectants was compared in the presence or absence of drugs. The relative transcript quantities of all rif loci did not change profoundly between strains grown with or without drug. We conclude that either there is no crosstalk between rif loci or that the elusive system of allelic exclusion of rif gene transcription is not controlled by their 5′ upstream region alone.

Single-cell RNA-seq reveals hidden transcriptional variation in malaria parasites

Single-cell RNA-sequencing is revolutionising our understanding of seemingly homogeneous cell populations but has not yet been widely applied to single-celled organisms. Transcriptional variation in unicellular malaria parasites from the Plasmodium genus is associated with critical phenotypes including red blood cell invasion and immune evasion, yet transcriptional variation at an individual parasite level has not been examined in depth. Here, we describe the adaptation of a single-cell RNA-sequencing (scRNA-seq) protocol to deconvolute transcriptional variation for more than 500 individual parasites of both rodent and human malaria comprising asexual and sexual life-cycle stages. We uncover previously hidden discrete transcriptional signatures during the pathogenic part of the life cycle, suggesting that expression over development is not as continuous as commonly thought. In transmission stages, we find novel, sex-specific roles for differential expression of contingency gene families that are usually associated with immune evasion and pathogenesis.

Heterologous Expression of a Novel Drug Transporter from the Malaria Parasite Alters Resistance to Quinoline Antimalarials

Antimalarial drug resistance hampers effective malaria treatment. Critical SNPs in a particular, putative amino acid transporter were recently linked to chloroquine (CQ) resistance in malaria parasites. Here, we show that this conserved protein (PF3D7_0629500 in Plasmodium falciparum; AAT1 in P. chabaudi) is a structural homologue of the yeast amino acid transporter Tat2p, which is known to mediate quinine uptake and toxicity. Heterologous expression of PF3D7_0629500 in yeast produced CQ hypersensitivity, coincident with increased CQ uptake. PF3D7_0629500-expressing cultures were also sensitized to related antimalarials; amodiaquine, mefloquine and particularly quinine. Drug sensitivity was reversed by introducing a SNP linked to CQ resistance in the parasite. Like Tat2p, PF3D7_0629500-dependent quinine hypersensitivity was suppressible with tryptophan, consistent with a common transport mechanism. A four-fold increase in quinine uptake by PF3D7_0629500 expressing cells was abolished by the resistance SNP. The parasite protein localised primarily to the yeast plasma membrane. Its expression varied between cells and this heterogeneity was used to show that high-expressing cell subpopulations were the most drug sensitive. The results reveal that the PF3D7_0629500 protein can determine the level of sensitivity to several major quinine-related antimalarials through an amino acid-inhibitable drug transport function. The potential clinical relevance is discussed.

Plasmodium knowlesi: a superb in vivo nonhuman primate model of antigenic variation in malaria

Antigenic variation in malaria was discovered in Plasmodium knowlesi studies involving longitudinal infections of rhesus macaques (M. mulatta). The variant proteins, known as the P. knowlesi Schizont Infected Cell Agglutination (SICA) antigens and the P. falciparum Erythrocyte Membrane Protein 1 (PfEMP1) antigens, expressed by the SICAvar and var multigene families, respectively, have been studied for over 30 years. Expression of the SICA antigens in P. knowlesi requires a splenic component, and specific antibodies are necessary for variant antigen switch events in vivo. Outstanding questions revolve around the role of the spleen and the mechanisms by which the expression of these variant antigen families are regulated. Importantly, the longitudinal dynamics and molecular mechanisms that govern variant antigen expression can be studied with P. knowlesi infection of its mammalian and vector hosts. Synchronous infections can be initiated with established clones and studied at multi-omic levels, with the benefit of computational tools from systems biology that permit the integration of datasets and the design of explanatory, predictive mathematical models. Here we provide an historical account of this topic, while highlighting the potential for maximizing the use of P. knowlesi - macaque model systems and summarizing exciting new progress in this area of research.