Falciparum Genome Research Articles

A transcription factor helps Plasmodium falciparum gametocytogenesis take shape

The recent study by Campelo Morillo et al. has shown that one of the small number of non-ApiAP2 DNA-binding proteins in the Plasmodium falciparum genome acts as a transcription factor in the gametocytogenesis cascade and is responsible for the gametocyte's distinctive morphology.

Distinct Histone Post-translational Modifications during Plasmodium falciparum Gametocyte Development.

Histones are the building units of nucleosomes, which constitute chromatin. Histone post-translational modifications (PTMs) play an essential role in epigenetic gene regulation. The Plasmodium falciparum genome encodes canonical and variant histones and a collection of conserved enzymes for histone PTMs and chromatin remodeling. Herein, we profiled the P. falciparum histone PTMs during the development of gametocytes, the obligatory stage for parasite transmission. Mass spectrometric analysis of histones extracted from the early, middle, and late stages of gametocytes identified 457 unique histone peptides with 90 PTMs, of which 50% were novel. The gametocyte histone PTMs display distinct patterns from asexual stages, with many new methylation sites in histones H3 and H3.3 (e.g., K14, K18, and K37). Quantitative analyses revealed a high abundance of acetylation in H3 and H4, mono-methylation of H3/H3.3 K37, and ubiquitination of H3BK112, suggesting that these PTMs play critical roles in gametocytes. Gametocyte histones also showed extensive and unique combinations of PTMs. These data indicate that the parasite harbors distinct transcription regulation mechanisms during gametocyte development and lay the foundation for further characterization of epigenetic regulation in the life cycle of the malaria parasite.

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.

Chromatin structure can introduce systematic biases in genome-wide analyses of Plasmodium falciparum

Background: The maintenance, regulation, and dynamics of heterochromatin in the human malaria parasite, Plasmodium falciparum, has drawn increasing attention due to its regulatory role in mutually exclusive virulence gene expression and the silencing of key developmental regulators. The advent of genome-wide analyses such as chromatin-immunoprecipitation followed by sequencing (ChIP-seq) has been instrumental in understanding chromatin composition; however, even in model organisms, ChIP-seq experiments are susceptible to intrinsic experimental biases arising from underlying chromatin structure. Methods: We performed a control ChIP-seq experiment, re-analyzed previously published ChIP-seq datasets and compared different analysis approaches to characterize biases of genome-wide analyses in P. falciparum. Results: We found that heterochromatic regions in input control samples used for ChIP-seq normalization are systematically underrepresented in regard to sequencing coverage across the P. falciparum genome. This underrepresentation, in combination with a non-specific or inefficient immunoprecipitation, can lead to the identification of false enrichment and peaks across these regions. We observed that such biases can also be seen at background levels in specific and efficient ChIP-seq experiments. We further report on how different read mapping approaches can also skew sequencing coverage within highly similar subtelomeric regions and virulence gene families. To ameliorate these issues, we discuss orthogonal methods that can be used to characterize bona fide chromatin-associated proteins. Conclusions: Our results highlight the impact of chromatin structure on genome-wide analyses in the parasite and the need for caution when characterizing chromatin-associated proteins and features.

Structural and Biophysical Characterization of Plasmodium falciparum Bromodomain Protein 1

Epigenetic regulatory mechanisms are important for the growth and survival of eukaryotic organisms. Post‐translational modifications found on histone proteins in the nucleosome function as a dynamic cell signaling mechanism that controls many essential functions such as cell growth, DNA replication, and gene expression. Bromodomains are conserved protein interaction modules that specifically recognize acetylated lysine residues on histones. Among five Plasmodium species, P. falciparum is the major cause of malaria cases worldwide. According to a WHO report in 2019 there were 409,000 deaths reported worldwide and the most affected age group was of children under 5 years of age. Interestingly, the P. falciparum genome encodes eight bromodomain‐containing proteins, and bromodomain protein 1 (PfBDP1) has been shown to play an important role in malaria pathogenesis. The domain architecture of PfBDP1 includes a single C‐terminal bromodomain and several ankyrin repeats. PfBDP1 is thought to act by tethering a transcriptional activator complex to acetylated histone H3 to control genes required for parasite invasion. Additionally, PfBDP1 forms a complex with the bromodomain protein 2 (PfBDP2) to carry out its cellular functions. We hypothesized that the bromodomain is crucial for the function of PfBDP1 via recognition of acetylated lysine modifications. We carried out several histone binding assays to identify novel histone interactions of PfBDP1 bromodomain. We also utilized a combination of structural biology and biophysical approaches to identify amino acid residues critical for ligand coordination. The rise of drug‐resistance towards P. falciparum has created an urgent need to outline the molecular mechanisms driving pathogenesis in hopes of uncovering new potential drug targets. Importantly, the results from this study may improve the development of new therapeutics to treat malaria.

Prediction and comparative analysis of CTCF binding sites based on a first principle approach

We calculated the patterns for the CCCTC transcription factor (CTCF) binding sites across many genomes on a first principle approach. The validation of the first principle method was done on the human as well as on the mouse genome. The predicted human CTCF binding sites are consistent with the consensus sequence, ChIP-seq data for the K562 cell, nucleosome positions for IMR90 cell as well as the CTCF binding sites in the mouse HOXA gene. The analysis of Homo sapiens, Mus musculus, Sus scrofa, Capra hircus and Drosophila melanogaster whole genomes shows: binding sites are organized in cluster-like groups, where two consecutive sites obey a power-law with coefficient ranging from 0.3292 ± 0.0068 to 0.5409 ± 0.0064; the distance between these groups varies from 18.08 ± 0.52 kbp to 42.1 ± 2.0 kbp. The genome of Aedes aegypti does not show a power law, but 19.9% of binding sites are 144 ± 4 and 287 ± 5 bp distant of each other. We run negative tests, confirming the under-representation of CTCF binding sites in Caenorhabditis elegans, Plasmodium falciparum and Arabidopsis thaliana complete genomes.

Prefoldin subunit 6 of Plasmodium falciparum binds merozoite surface protein-1.

Malaria is a human disease caused by eukaryotic protozoan parasites of the Plasmodium genus. Plasmodium falciparum (Pf) causes the most lethal form of human malaria and is responsible for widespread mortality worldwide. Prefoldin is a heterohexameric molecular complex that binds and delivers unfolded proteins to chaperonin for correct folding. The prefoldin PFD6 is predicted to interact with merozoite surface protein-1 (MSP-1), a protein well known to play a pivotal role in erythrocyte binding and invasion by Plasmodium merozoites. We previously found that the P. falciparum (Pf) genome contains six prefoldin genes and a prefoldin-like gene whose molecular functions are unidentified. Here, we analyzed the expression of PfPFD-6 during the asexual blood stages of the parasite and investigated its interacting partners. PfPFD-6 was found to be significantly expressed at the trophozoite and schizont stages. Pull-down assays suggest PfPFD-6 interacts with MSP-1. In silico analysis suggested critical residues involved in the PfPFD-6-MSP-1 interaction. Our data suggest PfPFD-6 may play a role in stabilizing or trafficking MSP-1.

GeneTargeter: Automated In Silico Design for Genome Editing in the Malaria Parasite, Plasmodium falciparum.

Functional characterization of the multitude of poorly described proteins in the human malarial pathogen, Plasmodium falciparum, requires tools to enable genome-scale perturbation studies. Here, we present GeneTargeter (genetargeter.mit.edu), a software tool for automating the design of homology-directed repair donor vectors to achieve gene knockouts, conditional knockdowns, and epitope tagging of P. falciparum genes. We demonstrate GeneTargeter-facilitated genome-scale design of six different types of knockout and conditional knockdown constructs for the P. falciparum genome and validate the computational design process experimentally with successful donor vector assembly and transfection. The software's modular nature accommodates arbitrary destination vectors and allows customizable designs that extend the genome manipulation outcomes attainable in Plasmodium and other organisms.

Is the AT-rich DNA of malaria parasites a drug target?

A recent antimalarial screen by Alder and colleagues has uncovered a natural product, PDE-I2, with DNA-binding and schizonticidal activity against Plasmodium falciparum. Parasite specificity is likely conferred by the extremely high A+T content of the P. falciparum genome. We discuss here this character as a potential target for future drugs.

Analysis of var Gene Transcript Patterns by Quantitative Real-Time PCR.

Quantitative real-time PCR (qPCR) is a simple and sensitive method for determining the amount of a specific target DNA sequence present in a sample. Compared to RNA-seq, reverse transcription qPCR (RT-qPCR) is fast, requires only low input material and is easy to analyze. Therefore, qPCR is widely used to analyze gene expressionin P. falciparum, including analyses of the multicopy gene families encoding variant surface antigens (VSAs), whose expression is clonally variant and prone to changes over time. In the recent years, several P. falciparum genomes of culture-adapted strains have been sequenced, providing the knowledge to design variable gene family-specific qPCR primers for each P. falciparum genetic background. Here, we describe the required materials, methods and key factors to perform RT-qPCR experiments to determine VSA transcript abundances in the P. falciparum clones 3D7/NF54, IT4, HB3, and 7G8.

Zinc finger proteins of Plasmodium falciparum.

Zinc finger proteins (ZFPs) are a large diverse family of proteins with one or more zinc finger domains in which zinc is important in stabilising the domain. ZFPs can interact with DNA, RNA, lipids or even other proteins and therefore contribute to diverse cellular processes including transcriptional regulation, ubiquitin-mediated protein degradation, mRNA decay and stability. In this review, we provide the first comprehensive classification of ZFPs of the malaria parasite Plasmodium falciparum and provide a state of knowledge on the main ZFPs in the parasite, which include the C2H2, CCCH, RING finger and the PHD finger proteins. TAKE AWAYS: The Plasmodium falciparum genome encodes 170 putative Zinc finger proteins (ZFPs). The C2H2, CCCH, RING finger and PHD finger subfamilies of ZFPs are most represented. Known ZFP functions include the regulation of mRNA metabolism and proteostasis.

Identifying transcript 5'capped ends in Plasmodium falciparum.

BackgroundThe genome of the human malaria parasite Plasmodium falciparum is poorly annotated, in particular, the 5′ capped ends of its mRNA transcripts. New approaches are needed to fully catalog P. falciparum transcripts for understanding gene function and regulation in this organism.MethodsWe developed a transcriptomic method based on next-generation sequencing of complementary DNA (cDNA) enriched for full-length fragments using eIF4E, a 5′ cap-binding protein, and an unenriched control. DNA sequencing adapter was added after enrichment of full-length cDNA using two different ligation protocols. From the mapped sequence reads, enrichment scores were calculated for all transcribed nucleotides and used to calculate P-values of 5′ capped nucleotide enrichment. Sensitivity and accuracy were increased by combining P-values from replicate experiments. Data were obtained for P. falciparum ring, trophozoite and schizont stages of intra-erythrocytic development.Results5′ capped nucleotide signals were mapped to 17,961 non-overlapping P. falciparum genomic intervals. Analysis of the dominant 5′ capped nucleotide in these genomic intervals revealed the presence of two groups with distinctive epigenetic features and sequence patterns. A total of 4,512 transcripts were annotated as 5′ capped based on the correspondence of 5′ end with 5′ capped nucleotide annotated from full-length cDNA data.DiscussionThe presence of two groups of 5′ capped nucleotides suggests that alternative mechanisms may exist for producing 5′ capped transcript ends in P. falciparum. The 5′ capped transcripts that are antisense, outside of, or partially overlapping coding regions may be important regulators of gene function in P. falciparum.

Prediction of gene expression regulation by human microRNAs in Plasmodium falciparum

BackgroundMalaria is a disease annually causing over 400,000 deaths. Deep understanding of molecular and genetic processes underlying its life cycle and pathogenicity is required to efficiently resist it. RNA interference is a mechanism of the gene expression regulation typical for a wide variety of species. Even though the existence of this phenomenon in Plasmodium falciparum has long been rejected, several recent works pose hypotheses and provide direct and indirect evidence of the existence of mechanisms similar to RNA interference in this organism. In particular, the possibility of regulation of P. falciparum gene expression through human microRNAs is of great importance both for fundamental biology and for medicine. In the present work we address the problem of possibility of the existence in the P. falciparum genome of the nucleotide sequences such that mRNAs transcribed from genes containing these sequences could form duplexes with human microRNAs. Using bioinformatics methods we have analysed genomes of 15 P. falciparum isolates for sequences homological to these microRNAs.ResultsThe analysis has demonstrated the existence of a vast number of genes that could potentially be regulated by the human microRNAs in the plasmodial genome.ConclusionsDespite the fact that the numbers of homological intervals vary significantly between isolates, the hsa-miR-451a and hsa-miR-223-3p microRNAs are expected to make the most notable contribution to the pathogenesis of P. falciparum malaria. The majority of homological intervals occur in genes encoding cell adhesion proteins.

Complement C1s cleaves PfEMP1 at interdomain conserved sites inhibiting Plasmodium falciparum cytoadherence

Cytoadhesion of Plasmodium falciparum-infected erythrocytes (IEs) to the endothelial lining of blood vessels protects parasites from splenic destruction, but also leads to detrimental inflammation and vessel occlusion. Surface display of the P. falciparum erythrocyte membrane protein 1 (PfEMP1) adhesion ligands exposes them to host antibodies and serum proteins. PfEMP1 are important targets of acquired immunity to malaria, and through evolution, the protein family has expanded and diversified to bind a select set of host receptors through antigenically diversified receptor-binding domains. Here, we show that complement component 1s (C1s) in serum cleaves PfEMP1 at semiconserved arginine motifs located at interdomain regions between the receptor-binding domains, rendering the IE incapable of binding the two main PfEMP1 receptors, CD36 and endothelial protein C receptor (EPCR). Bioinformatic analyses of PfEMP1 protein sequences from 15 P. falciparum genomes found the C1s motif was present in most PfEMP1 variants. Prediction of C1s cleavage and loss of binding to endothelial receptors was further corroborated by testing of several different parasite lines. These observations suggest that the parasites have maintained susceptibility for cleavage by the serine protease, C1s, and provides evidence for a complex relationship between the complement system and the P. falciparum cytoadhesion virulence determinant.

CAMP-Dependent Signaling Pathways as Potential Targets for Inhibition of Plasmodium falciparum Blood Stages.

We review the role of signaling pathways in regulation of the key processes of merozoite egress and red blood cell invasion by Plasmodium falciparum and, in particular, the importance of the second messengers, cAMP and Ca2+, and cyclic nucleotide dependent kinases. cAMP-dependent protein kinase (PKA) is comprised of cAMP-binding regulatory, and catalytic subunits. The less well conserved cAMP-binding pockets should make cAMP analogs attractive drug leads, but this approach is compromised by the poor membrane permeability of cyclic nucleotides. We discuss how the conserved nature of ATP-binding pockets makes ATP analogs inherently prone to off-target effects and how ATP analogs and genetic manipulation can be useful research tools to examine this. We suggest that targeting PKA interaction partners as well as substrates, or developing inhibitors based on PKA interaction sites or phosphorylation sites in PKA substrates, may provide viable alternative approaches for the development of anti-malarial drugs. Proximity of PKA to a substrate is necessary for substrate phosphorylation, but the P. falciparum genome encodes few recognizable A-kinase anchor proteins (AKAPs), suggesting the importance of PKA-regulatory subunit myristylation and membrane association in determining substrate preference. We also discuss how Pf14-3-3 assembles a phosphorylation-dependent signaling complex that includes PKA and calcium dependent protein kinase 1 (CDPK1) and how this complex may be critical for merozoite invasion, and a target to block parasite growth. We compare altered phosphorylation levels in intracellular and egressed merozoites to identify potential PKA substrates. Finally, as host PKA may have a critical role in supporting intracellular parasite development, we discuss its role at other stages of the life cycle, as well as in other apicomplexan infections. Throughout our review we propose possible new directions for the therapeutic exploitation of cAMP-PKA-signaling in malaria and other diseases caused by apicomplexan parasites.

Development of a new barcode-based, multiplex-PCR, next-generation-sequencing assay and data processing and analytical pipeline for multiplicity of infection detection of Plasmodium falciparum

BackgroundSimultaneous infection with multiple malaria parasite strains is common in high transmission areas. Quantifying the number of strains per host, or the multiplicity of infection (MOI), provides additional parasite indices for assessing transmission levels but it is challenging to measure accurately with current tools. This paper presents new laboratory and analytical methods for estimating the MOI of Plasmodium falciparum.MethodsBased on 24 single nucleotide polymorphisms (SNPs) previously identified as stable, unlinked targets across 12 of the 14 chromosomes within P. falciparum genome, three multiplex PCRs of short target regions and subsequent next generation sequencing (NGS) of the amplicons were developed. A bioinformatics pipeline including B4Screening pathway removed spurious amplicons to ensure consistent frequency calls at each SNP location, compiled amplicons by SNP site diversity, and performed algorithmic haplotype and strain reconstruction. The pipeline was validated by 108 samples generated from cultured-laboratory strain mixtures in different proportions and concentrations, with and without pre-amplification, and using whole blood and dried blood spots (DBS). The pipeline was applied to 273 smear-positive samples from surveys conducted in western Kenya, then providing results into StrainRecon Thresholding for Infection Multiplicity (STIM), a novel MOI estimator.ResultsThe 24 barcode SNPs were successfully identified uniformly across the 12 chromosomes of P. falciparum in a sample using the pipeline. Pre-amplification and parasite concentration, while non-linearly associated with SNP read depth, did not influence the SNP frequency calls. Based on consistent SNP frequency calls at targeted locations, the algorithmic strain reconstruction for each laboratory-mixed sample had 98.5% accuracy in dominant strains. STIM detected up to 5 strains in field samples from western Kenya and showed declining MOI over time (q < 0.02), from 4.32 strains per infected person in 1996 to 4.01, 3.56 and 3.35 in 2001, 2007 and 2012, and a reduction in the proportion of samples with 5 strains from 57% in 1996 to 18% in 2012.ConclusionThe combined approach of new multiplex PCRs and NGS, the unique bioinformatics pipeline and STIM could identify 24 barcode SNPs of P. falciparum correctly and consistently. The methodology could be applied to field samples to reliably measure temporal changes in MOI.

An integrated platform for genome engineering and gene expression perturbation in Plasmodium falciparum

Establishing robust genome engineering methods in the malarial parasite, Plasmodium falciparum, has the potential to substantially improve the efficiency with which we gain understanding of this pathogen’s biology to propel treatment and elimination efforts. Methods for manipulating gene expression and engineering the P. falciparum genome have been validated. However, a significant barrier to fully leveraging these advances is the difficulty associated with assembling the extremely high AT content DNA constructs required for modifying the P. falciparum genome. These are frequently unstable in commonly-used circular plasmids. We address this bottleneck by devising a DNA assembly framework leveraging the improved reliability with which large AT-rich regions can be efficiently manipulated in linear plasmids. This framework integrates several key functional genetics outcomes via CRISPR/Cas9 and other methods from a common, validated framework. Overall, this molecular toolkit enables P. falciparum genetics broadly and facilitates deeper interrogation of parasite genes involved in diverse biological processes.

Elucidation of an essential function of the unique charged domain of Plasmodium topoisomerase III.

Topoisomerase III (TopoIII) along with RecQ helicases are required for the resolution of abnormal DNA structures that result from the stalling of replication forks. Sequence analyses have identified a putative TopoIII in the Plasmodium falciparum genome (PfTopoIII). PfTopoIII shows dual nuclear and mitochondrial localization. The expression and association of PfTopoIII with mtDNA is tightly linked to the asexual replication of the parasite. In this study, we observed that PfTopoIII physically interacts with PfBlm and PfWrn. Sequence alignment and domain analyses have revealed that it contains a unique positively charged region, spanning 85 amino acids, within domain II. A molecular dynamics simulation study revealed that this unstructured domain communicates with DNA and attains a thermodynamically stable state upon DNA binding. Here, we found that the association between PfTopoIII and the mitochondrial genome is negatively affected by the absence of the charged domain. Our study shows that PfTOPOIII can completely rescue the slow growth phenotype of the ΔtopoIII strain in Saccharomyces cerevisiae, but neither PfY421FtopoIII (catalytic-active site mutant) nor Pf(Δ259-337)topoIII (charged region deletion mutant) can functionally complement ScTOPOIII. Hydroxyurea (HU) led to stalling of the replication fork during the S phase, caused moderate toxicity to the growth of P. falciparum, and was associated with concomitant transcriptional upregulation of PfTOPOIII. In addition, ectopic expression of PfTOPOIII reversed HU-induced toxicity. Interestingly, the expression of Pf(Δ259-337)topoIII failed to reverse HU-mediated toxicity. Taken together, our results establish the importance of TopoIII during Plasmodium replication and emphasize the essential requirement of the charged domain in PfTopoIII function.

Purification and initial characterization of Plasmodium falciparum K+ channels, PfKch1 and PfKch2 produced in Saccharomyces cerevisiae

Resistance towards known antimalarial drugs poses a significant problem, urging for novel drugs that target vital proteins in the malaria parasite Plasmodium falciparum. However, recombinant production of malaria proteins is notoriously difficult. To address this, we have investigated two putative K+ channels, PfKch1 and PfKch2, identified in the P. falciparum genome. We show that PfKch1 and PfKch2 and a C-terminally truncated version of PfKch1 (PfKch11−1094) could indeed be functionally expressed in vivo, since a K+-uptake deficient Saccharomyces cerevisiae strain was complemented by the P. falciparum cDNAs. PfKch11−1094-GFP and GFP-PfKch2 fusion proteins were overexpressed in yeast, purified and reconstituted in lipid bilayers to determine their electrophysiological activity. Single channel conductance amounted to 16 ± 1 pS for PfKch11−1094-GFP and 28 ± 2 pS for GFP-PfKch2. We predicted regulator of K+-conductance (RCK) domains in the C-terminals of both channels, and we accordingly measured channel activity in the presence of Ca2+.

The human malaria parasite genome is configured into thousands of coexpressed linear regulatory units

The human malaria parasite Plasmodium falciparum thrives in radically different host environments in mosquitoes and humans, with only a limited set of transcription factors. The nature of regulatory elements or their target genes in the P. falciparum genome remains elusive. Here, we found that this eukaryotic parasite uses an efficient way to maximally use genetic and epigenetic regulation to form regulatory units (RUs) during blood infections. Genes located in the same RU tend to have the same pattern of expression over time and are associated with open chromatin along regulatory elements. To precisely define and quantify these RUs, a novel hidden Markov model was developed to capture the regulatory structure in a genome-wide fashion by integrating expression and epigenetic evidence. We successfully identified thousands of RUs and cross-validated with previous findings. We found more genes involved in red blood cell (RBC) invasion located in the same RU as the PfAP2-I (AP2-I) transcription factor, demonstrating that AP2-I is responsible for regulating RBC invasion. Our study has provided a regulatory mechanism for a compact eukaryotic genome and offers new insights into the in vivo transcriptional regulation of the P. falciparum intraerythrocytic stage.