Transcriptional variation in malaria parasites: why and how

Expression of the MSPDBL2 antigen in a discrete subset of Plasmodium falciparum schizonts is regulated by GDV1 but may not be linked to sexual commitment

The Plasmodium falciparum merozoite surface protein MSPDBL2 is a polymorphic antigen targeted by acquired immune responses, and normally expressed in only a minority of mature schizonts. The potential relationship of MSPDBL2 to sexual commitment is examined, as variable mspdbl2 transcript levels and proportions of MSPDBL2-positive mature schizonts in clinical isolates have previously correlated with levels of many sexual stage parasite gene transcripts, although not with the master regulator ap2-g. It is demonstrated that conditional overexpression of the gametocyte development protein GDV1, which promotes sexual commitment, also substantially increases the proportion of MSPDBL2-positive schizonts in culture. Conversely, truncation of the gdv1 gene is shown to prevent any expression of MSPDBL2. However, across diverse P. falciparum cultured lines, the variable proportions of MSPDBL2 positivity in schizonts do not correlate significantly with variable gametocyte conversion rates, indicating it is not involved in sexual commitment. Confirming this, examining a line with endogenous hemagglutinin-tagged AP2-G showed that the individual schizonts expressing MSPDBL2 are mostly different from those expressing AP2-G. Using a selection-linked integration system, modified P. falciparum lines were engineered to express an intact or disrupted version of MSPDBL2, showing the protein is not required for sexual commitment or early gametocyte development. Asexual parasite multiplication rates were also not affected by expression of either intact or disrupted MSPDBL2 in a majority of schizonts. Occurring alongside sexual commitment, the role of the discrete MSPDBL2-positive schizont subpopulation requires further investigation in natural infections where it is under immune selection.

IMPORTANCE

Malaria parasites in the blood are remarkably variable, able to switch antigenic targets so they may survive within humans who have already developed specific immune responses. This is one of the challenges in developing vaccines against malaria. MSPDBL2 is a target of naturally acquired immunity expressed in minority proportions of schizonts, the end stages of each 2-day replication cycle in red blood cells which contain merozoites prepared to invade new red blood cells. Results show that the proportion of schizonts expressing MSPDBL2 is positively controlled by the expression of the regulatory gametocyte development protein GDV1. It was previously known that expression of GDV1 leads to increased expression of AP2-G which causes parasites to switch to sexual development, so a surprising finding here is that MSPDBL2-positive parasites are mostly distinct from those that express AP2-G. This discrete antigenic subpopulation of mostly asexual parasites is regulated alongside sexually committed parasites, potentially enabling survival under stress conditions.

Comparative analysis of peripheral whole blood transcriptome from asymptomatic carriers reveals upregulation of subsets of surface proteins implicated in Plasmodium falciparum phenotypic plasticity

The molecular mechanism underlying Plasmodium falciparum's persistence in the asymptomatic phase of infection remains largely unknown. However, large-scale shifts in the parasites' gene expression during asymptomatic infections may enhance phenotypic plasticity, maximizing their fitness and leading to the persistence of the asymptomatic infections. To uncover these mechanisms, we aimed to identify parasite genetic factors implicated in asymptomatic infections through whole transcriptome analysis. We analyzed publicly available transcriptome datasets containing asymptomatic malaria (ASM), uncomplicated malaria (SM), and malaria-naïve (NSM) samples from 35 subjects for differentially expressed genes (DEGs) and long noncoding RNAs. Our analysis identified 755 and 1773 DEGs in ASM vs SM and NSM, respectively. These DEGs revealed sets of genes coding for proteins of unknown functions (PUFs) upregulated in ASM vs SM and ASM, suggesting their role in underlying fundamental molecular mechanisms during asymptomatic infections. Upregulated genes in ASM vs SM revealed a subset of 24 clonal variant genes (CVGs) involved in host-parasite and symbiotic interactions and modulation of the symbiont of host erythrocyte aggregation pathways. Moreover, we identified 237 differentially expressed noncoding RNAs in ASM vs SM, of which 11 were found to interact with CVGs, suggesting their possible role in regulating the expression of CVGs. Our results suggest that P. falciparum utilizes phenotypic plasticity as an adaptive mechanism during asymptomatic infections by upregulating clonal variant genes, with long noncoding RNAs possibly playing a crucial role in their regulation. Thus, our study provides insights into the parasites' genetic factors that confer a fitness advantage during asymptomatic infections.

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.

Naturally Acquired Kelch13 Mutations in Plasmodium falciparum Strains Modulate In Vitro Ring-Stage Artemisinin-Based Drug Tolerance and Parasite Survival in Response to Hyperoxia

The ring-stage survival assay was utilized to assess the impact of physiological hyperoxic stress on dihydroartemisinin (DHA) tolerance for a panel of Plasmodium falciparum strains with and without Kelch13 mutations. Strains without naturally acquired Kelch13 mutations or the postulated genetic background associated with delayed parasite clearance time demonstrated reduced proliferation under hyperoxic conditions in the subsequent proliferation cycle. Dihydroartemisinin tolerance in three isolates with naturally acquired Kelch13 mutations but not two genetically manipulated laboratory strains was modulated by in vitro hyperoxic stress exposure of early-ring-stage parasites in the cycle before drug exposure. Reduced parasite tolerance to additional derivatives, including artemisinin, artesunate, and OZ277, was observed within the second proliferation cycle. OZ439 and epoxomicin completely prevented parasite survival under both hyperoxia and normoxic in vitro culture conditions, highlighting the unique relationship between DHA tolerance and Kelch13 mutation-associated genetic background. IMPORTANCE Artemisinin-based combination therapy (ACT) for treating malaria is under intense scrutiny following treatment failures in the Greater Mekong subregion of Asia. This is further compounded by the potential for extensive loss of life if treatment failures extend to the African continent. Although Plasmodium falciparum has become resistant to all antimalarial drugs, artemisinin "resistance" does not present in the same way as resistance to other antimalarial drugs. Instead, a partial resistance or tolerance is demonstrated, associated with the parasite's genetic profile and linked to a molecular marker referred to as K13. It is suggested that parasites may have adapted to drug treatment, as well as the presence of underlying population health issues such as hemoglobinopathies, and/or environmental pressures, resulting in parasite tolerance to ACT. Understanding parasite evolution and control of artemisinin tolerance will provide innovative approaches to mitigate the development of artemisinin tolerance and thereby artemisinin-based drug treatment failure and loss of life globally to malaria infections.

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.

Kinetic Tracking of Plasmodium falciparum Antigens on Infected Erythrocytes with a Novel Reporter of Protein Insertion and Surface Exposure

Intracellular malaria parasites export many proteins into their host cell, inserting several into the erythrocyte plasma membrane to enable interactions with their external environment. While static techniques have identified some surface-exposed proteins, other candidates have eluded definitive localization and membrane topology determination. Moreover, both export kinetics and the mechanisms of membrane insertion remain largely unexplored. We introduce Reporter of Insertion and Surface Exposure (RISE), a method for continuous nondestructive tracking of antigen exposure on infected cells. RISE utilizes a small 11-amino acid (aa) HiBit fragment of NanoLuc inserted into a target protein and detects surface exposure through high-affinity complementation to produce luminescence. We tracked the export and surface exposure of CLAG3, a parasite protein linked to nutrient uptake, throughout the Plasmodiumfalciparum cycle in human erythrocytes. Our approach revealed key determinants of trafficking and surface exposure. Removal of a C-terminal transmembrane domain aborted export. Unexpectedly, certain increases in the exposed reporter size improved the luminescence signal, but other changes abolished the surface signal, revealing that both size and charge of the extracellular epitope influence membrane insertion. Marked cell-to-cell variation with larger inserts containing multiple HiBit epitopes suggests complex regulation of CLAG3 insertion at the host membrane. Quantitative, continuous tracking of CLAG3 surface exposure thus reveals multiple factors that determine this protein’s trafficking and insertion at the host erythrocyte membrane. The RISE assay will enable study of surface antigens from divergent intracellular pathogens.

Malaria parasites do respond to heat

The capacity of malaria parasites to respond to changes in their environment at the transcriptional level has been the subject of debate, but recent evidence has unambiguously demonstrated that Plasmodium spp. can produce adaptive transcriptional responses when exposed to some specific types of stress. These include metabolic conditions and febrile temperature. The Plasmodium falciparum protective response to thermal stress is similar to the response in other organisms, but it is regulated by a transcription factor evolutionarily unrelated to the conserved transcription factor that drives the heat shock (HS) response in most eukaryotes. Of the many genes that change expression during HS, only a subset constitutes an authentic response that contributes to parasite survival.

Satellite-Like W-Elements: Repetitive, Transcribed, and Putative Mobile Genetic Factors with Potential Roles for Biology and Evolution of Schistosoma mansoni

A large portion of animal and plant genomes consists of noncoding DNA. This part includes tandemly repeated sequences and gained attention because it offers exciting insights into genome biology. We investigated satellite-DNA elements of the platyhelminth Schistosoma mansoni, a parasite with remarkable biological features. Schistosoma mansoni lives in the vasculature of humans causing schistosomiasis, a disease of worldwide importance. Schistosomes are the only trematodes that have evolved separate sexes, and the sexual maturation of the female depends on constant pairing with the male. The schistosome karyotype comprises eight chromosome pairs, males are homogametic (ZZ) and females are heterogametic (ZW). Part of the repetitive DNA of S. mansoni are W-elements (WEs), originally discovered as female-specific satellite DNAs in the heterochromatic block of the W-chromosome. Based on new genome and transcriptome data, we performed a reanalysis of the W-element families (WEFs). Besides a new classification of 19 WEFs, we provide first evidence for stage-, sex-, pairing-, gonad-, and strain-specific/preferential transcription of WEs as well as their mobile nature, deduced from autosomal copies of full-length and partial WEs. Structural analyses suggested roles as sources of noncoding RNA-like hammerhead ribozymes, for which we obtained functional evidence. Finally, the variable WEF occurrence in different schistosome species revealed remarkable divergence. From these results, we propose that WEs potentially exert enduring influence on the biology of S. mansoni. Their variable occurrence in different strains, isolates, and species suggests that schistosome WEs may represent genetic factors taking effect on variability and evolution of the family Schistosomatidae.

A heat-shock response regulated by the PfAP2-HS transcription factor protects human malaria parasites from febrile temperatures

Periodic fever is a characteristic clinical feature of human malaria, but how parasites survive febrile episodes is not known. Although the genomes of Plasmodium species encode a full set of chaperones, they lack the conserved eukaryotic transcription factor HSF1, which activates the expression of chaperones following heat shock. Here, we show that PfAP2-HS, a transcription factor in the ApiAP2 family, regulates the protective heat-shock response in Plasmodium falciparum. PfAP2-HS activates the transcription of hsp70-1 and hsp90 at elevated temperatures. The main binding site of PfAP2-HS in the entire genome coincides with a tandem G-box DNA motif in the hsp70-1 promoter. Engineered parasites lacking PfAP2-HS have reduced heat-shock survival and severe growth defects at 37 °C but not at 35 °C. Parasites lacking PfAP2-HS also have increased sensitivity to imbalances in protein homeostasis (proteostasis) produced by artemisinin, the frontline antimalarial drug, or the proteasome inhibitor epoxomicin. We propose that PfAP2-HS contributes to the maintenance of proteostasis under basal conditions and upregulates specific chaperone-encoding genes at febrile temperatures to protect the parasite against protein damage.

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.

The ApiAP2 factor PfAP2-HC is an integral component of heterochromatin in the malaria parasite Plasmodium falciparum

Malaria parasites undergo a complex life cycle in the human host and the mosquito vector. The ApiAP2 family of DNA-binding proteins plays a dominant role in parasite development and life cycle progression. Most ApiAP2 factors studied to date act as transcription factors regulating stage-specific gene expression. Here, we characterized an ApiAP2 factor in Plasmodium falciparum that we termed PfAP2-HC. We demonstrate that PfAP2-HC specifically binds to heterochromatin throughout the genome. Intriguingly, PfAP2-HC does not bind DNA in vivo and recruitment of PfAP2-HC to heterochromatin is independent of its DNA-binding domain but strictly dependent on heterochromatin protein 1. Furthermore, our results suggest that PfAP2-HC functions neither in the regulation of gene expression nor in heterochromatin formation or maintenance. In summary, our findings reveal PfAP2-HC as a core component of heterochromatin in malaria parasites and identify unexpected properties and substantial functional divergence among the members of the ApiAP2 family of regulatory proteins.

•

The ApiAP2 factor AP2-HC is a core component of heterochromatin in malaria parasites

•

Binding of AP2-HC to heterochromatin strictly depends on heterochromatin protein 1

•

The AP2 DNA-binding domain of AP2-HC is dispensable for heterochromatin association

Direct Nanopore Sequencing of mRNA Reveals Landscape of Transcript Isoforms in Apicomplexan Parasites

Alternative splicing is a widespread phenomenon in metazoans by which single genes are able to produce multiple isoforms of the gene product. However, this has been poorly characterized in apicomplexans, a major phylum of some of the most important global parasites. Efforts have been hampered by atypical transcriptomic features, such as the high AU content of Plasmodium RNA, but also the limitations of short-read sequencing in deciphering complex splicing events. In this study, we utilized the long read direct RNA sequencing platform developed by Oxford Nanopore Technologies to survey the alternative splicing landscape of Toxoplasma gondii and Plasmodium falciparum. We find that while native RNA sequencing has a reduced throughput, it allows us to obtain full-length or nearly full-length transcripts with comparable quantification to Illumina sequencing. By comparing these data with available gene models, we find widespread alternative splicing, particularly intron retention, in these parasites. Most of these transcripts contain premature stop codons, suggesting that in these parasites, alternative splicing represents a pathway to transcriptomic diversity, rather than expanding proteomic diversity. Moreover, alternative splicing rates are comparable between parasites, suggesting a shared splicing machinery, despite notable transcriptomic differences between the parasites. This study highlights a strategy in using long-read sequencing to understand splicing events at the whole-transcript level and has implications in the future interpretation of transcriptome sequencing studies.

IMPORTANCE We have used a novel nanopore sequencing technology to directly analyze parasite transcriptomes. The very long reads of this technology reveal the full-length genes of the parasites that cause malaria and toxoplasmosis. Gene transcripts must be processed in a process called splicing before they can be translated to protein. Our analysis reveals that these parasites very frequently only partially process their gene products, in a manner that departs dramatically from their human hosts.

Using scRNA-seq to Identify Transcriptional Variation in the Malaria Parasite Ookinete Stage

The crossing of the mosquito midgut epithelium by the malaria parasite motile ookinete form represents the most extreme population bottleneck in the parasite life cycle and is a prime target for transmission blocking strategies. However, we have little understanding of the clonal variation that exists in a population of ookinetes in the vector, partially because the parasites are difficult to access and are found in low numbers. Within a vector, variation may result as a response to specific environmental cues or may exist independent of those cues as a potential bet-hedging strategy. Here we use single-cell RNA-seq to profile transcriptional variation in Plasmodium berghei ookinetes across different vector species, and between and within individual midguts. We then compare our results to low-input transcriptomes from individual Anopheles coluzzii midguts infected with the human malaria parasite Plasmodium falciparum. Although the vast majority of transcriptional changes in ookinetes are driven by development, we have identified candidate genes that may be responding to environmental cues or are clonally variant within a population. Our results illustrate the value of single-cell and low-input technologies in understanding clonal variation of parasite populations.

Transcriptional Analysis of Tightly Synchronized Plasmodium falciparum Intraerythrocytic Stages by RT-qPCR

In Plasmodium falciparum, the parasite responsible for the most severe forms of human malaria, many fundamental processes are controlled at the transcriptional level. Studies on diverse aspects of basic parasite biology as well as molecular epidemiology studies often rely on the ability to accurately measure transcript levels, but this is complicated by the cyclic expression patterns of the majority of malaria parasite genes. Here, we provide a complete workflow to measure transcript levels in P. falciparum intraerythrocytic blood stages, overcoming the confounding factors that are commonly encountered. The method described covers all the steps from synchronization of parasite cultures to reverse transcriptase quantitative PCR (RT-qPCR) analysis.

Intrinsic multiplication rate variation and plasticity of human blood stage malaria parasites

Pathogen multiplication rate is theoretically an important determinant of virulence, although often poorly understood and difficult to measure accurately. We show intrinsic asexual blood stage multiplication rate variation of the major human malaria parasite Plasmodium falciparum to be associated with blood-stage infection intensity in patients. A panel of clinical isolates from a highly endemic West African population was analysed repeatedly during five months of continuous laboratory culture, showing a range of exponential multiplication rates at all timepoints tested, mean rates increasing over time. All isolates had different genome sequences, many containing within-isolate diversity that decreased over time in culture, but increases in multiplication rates were not primarily attributable to genomic selection. New mutants, including premature stop codons emerging in a few isolates, did not attain sufficiently high frequencies to substantially affect overall multiplication rates. Significantly, multiplication rate variation among the isolates at each of the assayed culture timepoints robustly correlated with parasite levels seen in patients at clinical presentation, indicating innate parasite control of multiplication rate that contributes to virulence.

Lindsay Stewart et al. analyze clinical isolates of the human malaria parasite Plasmodium falciparum from a highly endemic West African population and show that intrinsic multiplication rate variation is associated with blood-stage infection intensity. Their results indicate that parasite control of multiplication contributes to virulence.

Assays for quantification of male and female gametocytes in human blood by qRT-PCR in the absence of pure sex-specific gametocyte standards

Background

Malaria transmission from humans to Anopheles mosquitoes requires the presence of gametocytes in human peripheral circulation, and the dynamics of transmission are determined largely by the density and sex ratio of the gametocytes. Molecular methods are thus employed to measure gametocyte densities, particularly when assessing transmission epidemiology and the efficacy of transmission-blocking interventions. However, accurate quantification of male and female gametocytes with molecular methods requires pure male and female gametocytes as reference standards, which are not widely available.

Methods

qRT-PCR assays were used to quantify levels of sex-specific mRNA transcripts in Plasmodium falciparum female and male gametocytes (pfs25 and pfMGET, respectively) using synthetic complimentary RNA standards and in vitro cultured gametocytes. Assays were validated and assay performance was investigated in blood samples of clinical trial participants using these standards and compared to absolute quantification by droplet digital PCR (ddPCR).

Results

The number of transcript copies per gametocyte were determined to be 279.3 (95% CI 253.5–307.6) for the female-specific transcript pfs25, and 12.5 (95% CI 10.6–14.9) for the male-specific transcript pfMGET. These numbers can be used to convert from transcript copies/mL to gametocyte/mL. The reportable range was determined to be 5.71 × 10⁶ to 5.71 female gametocytes/mL for pfs25, and 1.73 × 10⁷ to 1.73 × 10¹ male gametocytes/mL for pfMGET. The limit of detection was 3.9 (95% CI 2.5–8.2) female gametocytes/mL for pfs25, and 26.9 (95% CI 19.3–51.7) male gametocytes/mL for PfMGET. Both assays showed minimal intra-assay and inter-assay variability with coefficient of variation < 3%. No cross-reactivity was observed in both assays in uninfected human blood samples. Comparison of results from ddPCR to qRT-PCR assays on clinical blood samples indicated a high-level agreement (ICC = 0.998 for pfs25 and 0.995 for pfMGET).

Conclusions

This study reports the validation of qRT-PCR assays that are able to accurately quantify female and male P. falciparum gametocytes at sub-microscopic densities. The assays showed excellent reproducibility, sensitivity, precision, specificity, and accuracy. The methodology will enable the estimation of gametocyte density in the absence of pure female and male gametocyte standards, and will facilitate clinical trials and epidemiological studies.