Gene regulation in the intraerythrocytic cycle of Plasmodium falciparum

What functional genomics has taught us about transcriptional regulation in malaria parasites

Malaria parasites are characterized by a complex life cycle that is accompanied by dynamic gene expression patterns. The factors and mechanisms that regulate gene expression in these parasites have been searched for even before the advent of next generation sequencing technologies. Functional genomics approaches have substantially boosted this area of research and have yielded significant insights into the interplay between epigenetic, transcriptional and post-transcriptional mechanisms. Recently, considerable progress has been made in identifying sequence-specific transcription factors and DNA-encoded regulatory elements. Here, we review the insights obtained from these efforts including the characterization of core promoters, the involvement of sequence-specific transcription factors in life cycle progression and the mapping of gene regulatory elements. Furthermore, we discuss recent developments in the field of functional genomics and how they might contribute to further characterization of this complex gene regulatory network.

Parasite-Related Genetic and Epigenetic Aspects and Host Factors Influencing Plasmodium falciparum Invasion of Erythrocytes

Malaria, a disease caused by Plasmodium parasites, is widespread throughout tropical and sub-tropical regions worldwide; it mostly affects children and pregnant woman. Eradication has stalled despite effective prevention measures and medication being available for this disease; this has mainly been due to the parasite's resistance to medical treatment and the mosquito vector's resistance to insecticides. Tackling such resistance involves using renewed approaches and techniques for accruing a deep understanding of the parasite's biology, and developing new drugs and vaccines. Studying the parasite's invasion of erythrocytes should shed light on its ability to switch between invasion phenotypes related to the expression of gene sets encoding proteins acting as ligands during target cell invasion, thereby conferring mechanisms for evading a particular host's immune response and adapting to changes in target cell surface receptors. This review considers some factors influencing the expression of such phenotypes, such as Plasmodium's genetic, transcriptional and epigenetic characteristics, and explores some host-related aspects which could affect parasite phenotypes, aiming at integrating knowledge regarding this topic and the possible relationship between the parasite's biology and host factors playing a role in erythrocyte invasion.

Comparative transcriptomics of female and male gametocytes in Plasmodium berghei and the evolution of sex in alveolates

BACKGROUND

The clinical symptoms of malaria are caused by the asexual replication of Plasmodium parasites in the blood of the vertebrate host. To spread to new hosts, however, the malaria parasite must differentiate into sexual forms, termed gametocytes, which are ingested by a mosquito vector. Sexual differentiation produces either female or male gametocytes, and involves significant morphological and biochemical changes. These transformations prepare gametocytes for the rapid progression to gamete formation and fertilisation, which occur within 20 min of ingestion. Here we present the transcriptomes of asexual, female, and male gametocytes in P. berghei, and a comprehensive statistically-based differential-expression analysis of the transcriptional changes that underpin this sexual differentiation.

RESULTS

RNA-seq analysis revealed numerous differences in the transcriptomes of female and male gametocytes compared to asexual stages. Overall, there is net downregulation of transcripts in gametocytes compared to asexual stages, with this trend more marked in female gametocytes. Our analysis identified transcriptional changes in previously-characterised gametocyte-specific pathways, which validated our approach. We also detected many previously-unreported female- and male-specific pathways and genes. Transcriptional biases in stage and gender were then used to investigate sex-specificity and sexual dimorphism of Plasmodium in an evolutionary context. Sex-related gene expression is well conserved between Plasmodium species, but relatively poorly conserved in related organisms outside this genus. This pattern of conservation is most evident in genes necessary for both male and female gametocyte formation. However, this trend is less pronounced for male-specific genes, which are more highly conserved outside the genus than genes specific to female development.

CONCLUSIONS

We characterised the transcriptional changes that are integral to the development of the female and male sexual forms of Plasmodium. These differential-expression patterns provide a vital insight into understanding the gender-specific characteristics of this essential stage that is the primary target for treatments that block parasite transmission. Our results also offer insight into the evolution of sex genes through Alveolata, and suggest that many Plasmodium sex genes evolved within the genus. We further hypothesise that male gametocytes co-opted pre-existing cellular machinery in their evolutionary history, whereas female gametocytes evolved more through the development of novel, parasite-specific pathways.

Cluster analysis of Plasmodium RNA-seq time-course data identifies stage-specific co-regulated biological processes and regulatory elements

In this study, we interpreted RNA-seq time-course data of three developmental stages of Plasmodium species by clustering genes based on similarities in their expression profile without prior knowledge of the gene function. Functional enrichment of clusters of upregulated genes at specific time-points reveals potential targetable biological processes with information on their timings. We identified common consensus sequences that these clusters shared as potential points of coordinated transcriptional control. Five cluster groups showed upregulated profile patterns of biological interest. This included two clusters from the Intraerythrocytic Developmental Cycle (cluster 4 = 16 genes, and cluster 9 = 32 genes), one from the sexual development stage (cluster 2 = 851 genes), and two from the gamete-fertilization stage in the mosquito host (cluster 4 = 153 genes, and cluster 9 = 258 genes). The IDC expressed the least numbers of genes with only 1448 genes showing any significant activity of the 5020 genes (~29%) in the experiment. Gene ontology (GO) enrichment analysis of these clusters revealed a total of 671 uncharacterized genes implicated in 14 biological processes and components associated with these stages, some of which are currently being investigated as drug targets in on-going research. Five putative transcription regulatory binding motifs shared by members of each cluster were also identified, one of which was also identified in a previous study by separate researchers. Our study shows stage-specific genes and biological processes that may be important in antimalarial drug research efforts. In addition, timed-coordinated control of separate processes may explain the paucity of factors in parasites.

Triaging informative cis-regulatory elements for the combinatorial control of temporal gene expression during Plasmodium falciparum intraerythrocytic development

Background

Over 2700 genes are subject to stage-specific regulation during the intraerythrocytic development of the human malaria parasite Plasmodium falciparum. Bioinformatic analyses have identified a large number of over-represented motifs in the 5′ flanking regions of these genes that may act as cis-acting factors in the promoter-based control of temporal expression. Triaging these lists to provide candidates most likely to play a role in regulating temporal expression is challenging, but important if we are to effectively design in vitro studies to validate this role.

Methods

We report here the application of a repeated search of variations of 5′ flanking sequences from P. falciparum using the Finding Informative Regulatory Elements (FIRE) algorithm.

Results

Our approach repeatedly found a short-list of high scoring DNA motifs, for which cognate specific transcription factors were available, that appear to be typically associated with upregulation of mRNA accumulation during the first half of intraerythrocytic development.

Conclusions

We propose these cis-trans interactions may provide a combinatorial promoter-based control of gene expression to complement more global mechanisms of gene regulation that can account for temporal control during the second half of intraerythrocytic development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0701-0) contains supplementary material, which is available to authorized users.

Analysis of the spatial and temporal arrangement of transcripts over intergenic regions in the human malarial parasite Plasmodium falciparum

Background

The ability of the human malarial parasite Plasmodium falciparum to invade, colonise and multiply within diverse host environments, as well as to manifest its virulence within the human host, are activities tightly linked to the temporal and spatial control of gene expression. Yet, despite the wealth of high throughput transcriptomic data available for this organism there is very little information regarding the location of key transcriptional landmarks or their associated cis-acting regulatory elements. Here we provide a systematic exploration of the size and organisation of transcripts within intergenic regions to yield surrogate information regarding transcriptional landmarks, and to also explore the spatial and temporal organisation of transcripts over these poorly characterised genomic regions.

Results

Utilising the transcript data for a cohort of 105 genes we demonstrate that the untranscribed regions of mRNA are large and apportioned predominantly to the 5′ end of the open reading frame. Given the relatively compact size of the P. falciparum genome, we suggest that whilst transcriptional units are likely to spatially overlap, temporal co-transcription of adjacent transcriptional units is actually limited. Critically, the size of intergenic regions is directly dependent on the orientation of the two transcriptional units arrayed over them, an observation we extend to an analysis of the complete sequences of twelve additional organisms that share moderately compact genomes.

Conclusions

Our study provides a theoretical framework that extends our current understanding of the transcriptional landscape across the P. falciparum genome. Demonstration of a consensus gene-spacing rule that is shared between P. falciparum and ten other moderately compact genomes of apicomplexan parasites reveals the potential for our findings to have a wider impact across a phylum that contains many organisms important to human and veterinary health.

Dynamic clustering of gene expression

It is well accepted that genes are simultaneously involved in multiple biological processes and that genes are coordinated over the duration of such events. Unfortunately, clustering methodologies that group genes for the purpose of novel gene discovery fail to acknowledge the dynamic nature of biological processes and provide static clusters, even when the expression of genes is assessed across time or developmental stages. By taking advantage of techniques and theories from time frequency analysis, periodic gene expression profiles are dynamically clustered based on the assumption that different spectral frequencies characterize different biological processes. A two-step cluster validation approach is proposed to statistically estimate both the optimal number of clusters and to distinguish significant clusters from noise. The resulting clusters reveal coordinated coexpressed genes. This novel dynamic clustering approach has broad applicability to a vast range of sequential data scenarios where the order of the series is of interest.

Analysis of the molecular mechanisms governing the stage-specific expression of a prototypical housekeeping gene during intraerythrocytic development of P. falciparum

Gene expression during the intraerythrocytic development cycle of the human malarial parasite Plasmodium falciparum is subject to tight temporal control, resulting in a cascade of gene expression to meet the physiological demands of growth, replication, and reinvasion. The roles of the different molecular mechanisms that drive this temporal program of gene expression are poorly understood. Here we report the use of the bxb1 integrase system to reconstitute all aspects of the absolute and temporal control of the prototypical housekeeping gene encoding the proliferating cell nuclear antigen (Pfpcna) around an integrated luciferase reporter cassette. A quantitative analysis of the effect of the serial deletion of 5′ and 3′ genetic elements and sublethal doses of histone deacetylase inhibitors demonstrates that while the absolute control of gene expression could be perturbed, no effect on the temporal control of gene expression was observed. These data provide support for a novel model for the temporal control of potentially hundreds of genes during the intraerythrocytic development of this important human pathogen.

The Apicomplexan AP2 family: integral factors regulating Plasmodium development

Malaria is caused by protozoan parasites of the genus Plasmodium and involves infection of multiple hosts and cell types during the course of an infection. To complete its complex life cycle the parasite requires strict control of gene regulation for survival and successful propagation. Thus far, the Apicomplexan AP2 (ApiAP2) family of DNA-binding proteins is the sole family of proteins to have surfaced as candidate transcription factors in all apicomplexan species. Work from several laboratories is beginning to shed light on how the ApiAP2 proteins from Plasmodium spp. contribute to the regulation of gene expression at various stages of parasite development. Here we highlight recent progress toward understanding the role of Plasmodium ApiAP2 proteins in DNA related regulatory processes including transcriptional regulation and gene silencing.

Identification and genome-wide prediction of DNA binding specificities for the ApiAP2 family of regulators from the malaria parasite

The molecular mechanisms underlying transcriptional regulation in apicomplexan parasites remain poorly understood. Recently, the Apicomplexan AP2 (ApiAP2) family of DNA binding proteins was identified as a major class of transcriptional regulators that are found across all Apicomplexa. To gain insight into the regulatory role of these proteins in the malaria parasite, we have comprehensively surveyed the DNA-binding specificities of all 27 members of the ApiAP2 protein family from Plasmodium falciparum revealing unique binding preferences for the majority of these DNA binding proteins. In addition to high affinity primary motif interactions, we also observe interactions with secondary motifs. The ability of a number of ApiAP2 proteins to bind multiple, distinct motifs significantly increases the potential complexity of the transcriptional regulatory networks governed by the ApiAP2 family. Using these newly identified sequence motifs, we infer the trans-factors associated with previously reported plasmodial cis-elements and provide evidence that ApiAP2 proteins modulate key regulatory decisions at all stages of parasite development. Our results offer a detailed view of ApiAP2 DNA binding specificity and take the first step toward inferring comprehensive gene regulatory networks for P. falciparum.

Long- and short-term selective forces on malaria parasite genomes

Plasmodium parasites, the causal agents of malaria, result in more than 1 million deaths annually. Plasmodium are unicellular eukaryotes with small ∼23 Mb genomes encoding ∼5200 protein-coding genes. The protein-coding genes comprise about half of these genomes. Although evolutionary processes have a significant impact on malaria control, the selective pressures within Plasmodium genomes are poorly understood, particularly in the non-protein-coding portion of the genome. We use evolutionary methods to describe selective processes in both the coding and non-coding regions of these genomes. Based on genome alignments of seven Plasmodium species, we show that protein-coding, intergenic and intronic regions are all subject to purifying selection and we identify 670 conserved non-genic elements. We then use genome-wide polymorphism data from P. falciparum to describe short-term selective processes in this species and identify some candidate genes for balancing (diversifying) selection. Our analyses suggest that there are many functional elements in the non-genic regions of these genomes and that adaptive evolution has occurred more frequently in the protein-coding regions of the genome.

Revisiting the Plasmodium falciparum RIFIN family: from comparative genomics to 3D-model prediction

Background

Subtelomeric RIFIN genes constitute the most abundant multigene family in Plasmodium falciparum. RIFIN products are targets for the human immune response and contribute to the antigenic variability of the parasite. They are transmembrane proteins grouped into two sub-families (RIF_A and RIF_B). Although recent data show that RIF_A and RIF_B have different sub-cellular localisations and possibly different functions, the same structural organisation has been proposed for members of the two sub-families. Despite recent advances, our knowledge of the regulation of RIFIN gene expression is still poor and the biological role of the protein products remain obscure.

Results

Comparative studies on RIFINs in three clones of P. falciparum (3D7, HB3 and Dd2) by Multidimensional scaling (MDS) showed that gene sequences evolve differently in the 5'upstream, coding, and 3'downstream regions, and suggested a possible role of highly conserved 3' downstream sequences. Despite the expected polymorphism, we found that the overall structure of RIFIN repertoires is conserved among clones suggesting a balance between genetic drift and homogenisation mechanisms which guarantees emergence of novel variants but preserves the functionality of genes. Protein sequences from a bona fide set of 3D7 RIFINs were submitted to predictors of secondary structure elements. In contrast with the previously proposed structural organisation, no signal peptide and only one transmembrane helix were predicted for the majority of RIF_As. Finally, we developed a strategy to obtain a reliable 3D-model for RIF_As. We generated 265 possible structures from 53 non-redundant sequences, from which clustering and quality assessments selected two models as the most representative for putative RIFIN protein structures.

Conclusion

First, comparative analyses of RIFIN repertoires in different clones of P. falciparum provide insights on evolutionary mechanisms shaping the multigene family. Secondly, we found that members of the two sub-families RIF_As and RIF_Bs have different structural organization in accordance with recent experimental results. Finally, representative models for RIF_As have an "Armadillo-like" fold which is known to promote protein-protein interactions in diverse contexts.