Regulatory motifs uncovered among gene expression clusters in 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.

Chromatin Accessibility-Based Characterization of the Gene Regulatory Network Underlying Plasmodium falciparum Blood-Stage Development

Underlying the development of malaria parasites within erythrocytes and the resulting pathogenicity is a hardwired program that secures proper timing of gene transcription and production of functionally relevant proteins. How stage-specific gene expression is orchestrated in vivo remains unclear. Here, using the assay for transposase accessible chromatin sequencing (ATAC-seq), we identified ∼4,000 regulatory regions in P. falciparum intraerythrocytic stages. The vast majority of these sites are located within 2 kb upstream of transcribed genes and their chromatin accessibility pattern correlates positively with abundance of the respective mRNA transcript. Importantly, these regions are sufficient to drive stage-specific reporter gene expression and DNA motifs enriched in stage-specific sets of regulatory regions interact with members of the P. falciparum AP2 transcription factor family. Collectively, this study provides initial insights into the in vivo gene regulatory network of P. falciparum intraerythrocytic stages and should serve as a valuable resource for future studies.

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.

Expression of exogenous human hepatic nuclear factor-1α by a lentiviral vector and its interactions with Plasmodium falciparum subtilisin-like protease 2

The onset, severity, and ultimate outcome of malaria infection are influenced by parasite-expressed virulence factors as well as by individual host responses to these determinants. In both humans and mice, liver injury follows parasite entry, persisting to the erythrocytic stage in the case of infection with the fatal strain of Plasmodium falciparum. Hepatic nuclear factor (HNF)-1α is a master regulator of not only the liver damage and adaptive responses but also diverse metabolic functions. In this study, we analyzed the expression of host HNF-1α in relation to malaria infection and evaluated its interaction with the 5'-untranslated region of subtilisin-like protease 2 (subtilase, Sub2). Recombinant human HNF-1α expressed by a lentiviral vector (LV HNF-1α) was introduced into mice. Interestingly, differences in the activity of the 5'-untranslated region of the Pf-Sub2 promoter were detected in 293T cells, and LV HNF-1α was observed to influence promoter activity, suggesting that host HNF-1α interacts with the Sub2 gene.

Chromatin-driven de novo discovery of DNA binding motifs in the human malaria parasite

Background

Despite extensive efforts to discover transcription factors and their binding sites in the human malaria parasite Plasmodium falciparum, only a few transcription factor binding motifs have been experimentally validated to date. As a consequence, gene regulation in P. falciparum is still poorly understood. There is now evidence that the chromatin architecture plays an important role in transcriptional control in malaria.

Results

We propose a methodology for discovering cis-regulatory elements that uses for the first time exclusively dynamic chromatin remodeling data. Our method employs nucleosome positioning data collected at seven time points during the erythrocytic cycle of P. falciparum to discover putative DNA binding motifs and their transcription factor binding sites along with their associated clusters of target genes. Our approach results in 129 putative binding motifs within the promoter region of known genes. About 75% of those are novel, the remaining being highly similar to experimentally validated binding motifs. About half of the binding motifs reported show statistically significant enrichment in functional gene sets and strong positional bias in the promoter region.

Conclusion

Experimental results establish the principle that dynamic chromatin remodeling data can be used in lieu of gene expression data to discover binding motifs and their transcription factor binding sites. Our approach can be applied using only dynamic nucleosome positioning data, independent from any knowledge of gene function or expression.

MicroRNA genes preferentially expressed in dendritic cells contain sites for conserved transcription factor binding motifs in their promoters

Background

MicroRNAs (miRNAs) play a fundamental role in the regulation of gene expression by translational repression or target mRNA degradation. Regulatory elements in miRNA promoters are less well studied, but may reveal a link between their expression and a specific cell type.

Results

To explore this link in myeloid cells, miRNA expression profiles were generated from monocytes and dendritic cells (DCs). Differences in miRNA expression among monocytes, DCs and their stimulated progeny were observed. Furthermore, putative promoter regions of miRNAs that are significantly up-regulated in DCs were screened for Transcription Factor Binding Sites (TFBSs) based on TFBS motif matching score, the degree to which those TFBSs are over-represented in the promoters of the up-regulated miRNAs, and the extent of conservation of the TFBSs in mammals.

Conclusions

Analysis of evolutionarily conserved TFBSs in DC promoters revealed preferential clustering of sites within 500 bp upstream of the precursor miRNAs and that many mRNAs of cognate TFs of the conserved TFBSs were indeed expressed in the DCs. Taken together, our data provide evidence that selected miRNAs expressed in DCs have evolutionarily conserved TFBSs relevant to DC biology in their promoters.

Transcriptional activity of Plasmodium subtilisin-like protease 2 (Pf-Sub2) 5'untranslated regions and its interaction with hepatocyte growth factor

The onset, severity, and ultimate outcome of malaria infection are influenced by parasite-expressed virulence factors and individual host responses to these determinants. In both humans and mice, liver injury is involved after parasite entry, which persists until the erythrocyte stage after infection with the fatal strain Plasmodium falciparum (Pf). Hepatocyte growth factor (HGF) has strong anti-apoptotic effects in various kinds of cells, and also has diverse metabolic functions. In this work, Pf-subtilisin-like protease 2 (Pf-Sub2) 5'untranslated region (UTR) was analyzed and its transcriptional activity was estimated by luciferase expression. Fourteen TATA boxes were observed but only one Oct-1 and c-Myb were done. In addition, host HGF interaction with Pf-Sub2 was evaluated by co-transfection of HGF- and Pf-Sub2-cloned vector. Interestingly, -1,422/+12 UTR exhibited the strongest luciferase activity but -329 to +12 UTR did not exhibit luciferase activity. Moreover, as compared with the control of unexpressed HGF, the HGF protein suppressed luciferase expression driven by the 5'untranslated region of the Pf-Sub2 promoter. Taken together, it is suggested that HGF controls and interacts with the promoter region of the Pf-Sub2 gene.

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.

In silico and biological survey of transcription-associated proteins implicated in the transcriptional machinery during the erythrocytic development of Plasmodium falciparum

Background

Malaria is the most important parasitic disease in the world with approximately two million people dying every year, mostly due to Plasmodium falciparum infection. During its complex life cycle in the Anopheles vector and human host, the parasite requires the coordinated and modulated expression of diverse sets of genes involved in epigenetic, transcriptional and post-transcriptional regulation. However, despite the availability of the complete sequence of the Plasmodium falciparum genome, we are still quite ignorant about Plasmodium mechanisms of transcriptional gene regulation. This is due to the poor prediction of nuclear proteins, cognate DNA motifs and structures involved in transcription.

Results

A comprehensive directory of proteins reported to be potentially involved in Plasmodium transcriptional machinery was built from all in silico reports and databanks. The transcription-associated proteins were clustered in three main sets of factors: general transcription factors, chromatin-related proteins (structuring, remodelling and histone modifying enzymes), and specific transcription factors. Only a few of these factors have been molecularly analysed. Furthermore, from transcriptome and proteome data we modelled expression patterns of transcripts and corresponding proteins during the intra-erythrocytic cycle. Finally, an interactome of these proteins based either on in silico or on 2-yeast-hybrid experimental approaches is discussed.

Conclusion

This is the first attempt to build a comprehensive directory of potential transcription-associated proteins in Plasmodium. In addition, all complete transcriptome, proteome and interactome raw data were re-analysed, compared and discussed for a better comprehension of the complex biological processes of Plasmodium falciparum transcriptional regulation during the erythrocytic development.

Gene regulation in the intraerythrocytic cycle of Plasmodium falciparum

MOTIVATION

To date, there is little knowledge about one of the processes fundamental to the biology of Plasmodium falciparum, gene regulation including transcriptional control. We use noisy threshold models to identify regulatory sequence elements explaining membership to a gene expression cluster where each cluster consists of genes active during the part of the developmental cycle inside a red blood cell. Our approach is both able to capture the combinatorial nature of gene regulation and to incorporate uncertainty about the functionality of putative regulatory sequence elements.

RESULTS

We find a characteristic pattern where the most common motifs tend to be absent upstream of genes active in the first half of the cycle and present upstream of genes active in the second half. We find no evidence that motif's score, orientation, location and multiplicity improves prediction of gene expression. Through comparative genome analysis, we find a list of potential transcription factors and their associated motifs.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Patterns of gene-specific and total transcriptional activity during the Plasmodium falciparum intraerythrocytic developmental cycle

The relationships among gene regulatory mechanisms in the malaria parasite Plasmodium falciparum throughout its asexual intraerythrocytic developmental cycle (IDC) remain poorly understood. To investigate the level and nature of transcriptional activity and its role in controlling gene expression during the IDC, we performed nuclear run-on on whole-transcriptome samples from time points throughout the IDC and found a peak in RNA polymerase II-dependent transcriptional activity related to both the number of nuclei per parasite and variable transcriptional activity per nucleus over time. These differential total transcriptional activity levels allowed the calculation of the absolute transcriptional activities of individual genes from gene-specific nuclear run-on hybridization data. For half of the genes analyzed, sense-strand transcriptional activity peaked at the same time point as total activity. The antisense strands of several genes were substantially transcribed. Comparison of the transcriptional activity of the sense strand of each gene to its steady-state RNA abundance across the time points assayed revealed both correlations and discrepancies, implying transcriptional and posttranscriptional regulation, respectively. Our results demonstrate that such comparisons can effectively indicate gene regulatory mechanisms in P. falciparum and suggest that genes with diverse transcriptional activity levels and patterns combine to produce total transcriptional activity levels tied to parasite development during the IDC.

Identification of putative regulatory motifs in the upstream regions of co-expressed functional groups of genes in Plasmodium falciparum

BACKGROUND

Regulation of gene expression in Plasmodium falciparum (Pf) remains poorly understood. While over half the genes are estimated to be regulated at the transcriptional level, few regulatory motifs and transcription regulators have been found.

RESULTS

The study seeks to identify putative regulatory motifs in the upstream regions of 13 functional groups of genes expressed in the intraerythrocytic developmental cycle of Pf. Three motif-discovery programs were used for the purpose, and motifs were searched for only on the gene coding strand. Four motifs -- the 'G-rich', the 'C-rich', the 'TGTG' and the 'CACA' motifs -- were identified, and zero to all four of these occur in the 13 sets of upstream regions. The 'CACA motif' was absent in functional groups expressed during the ring to early trophozoite transition. For functional groups expressed in each transition, the motifs tended to be similar. Upstream motifs in some functional groups showed 'positional conservation' by occurring at similar positions relative to the translational start site (TLS); this increases their significance as regulatory motifs. In the ribonucleotide synthesis, mitochondrial, proteasome and organellar translation machinery genes, G-rich, C-rich, CACA and TGTG motifs, respectively, occur with striking positional conservation. In the organellar translation machinery group, G-rich motifs occur close to the TLS. The same motifs were sometimes identified for multiple functional groups; differences in location and abundance of the motifs appear to ensure different modes of action.

CONCLUSION

The identification of positionally conserved over-represented upstream motifs throws light on putative regulatory elements for transcription in Pf.

Control of gene expression in Plasmodium falciparum - ten years on

Ten years ago this journal published a review with an almost identical title detailing how the then recent introduction of transfection technology had advanced our understanding of the molecular control of transcriptional processes in Plasmodium falciparum, particularly in terms of promoter structure and function. In the succeeding years, sequencing of several Plasmodium spp. genomes and application of high throughput global postgenomic technologies have proven as significant, if not more, as has the ability to genetically manipulate these parasites in dissecting the molecular control of gene expression. Here we aim to review our current understanding of the control of gene expression in P. falciparum, including evidence available from other Plasmodium spp. and apicomplexan parasites. Specifically, however, we will address the current polarised debate regarding the level at which control is mediated, and attempt to identify some of the challenges this field faces in the next 10 years.

Regulatory hotspots in the malaria parasite genome dictate transcriptional variation

The determinants of transcriptional regulation in malaria parasites remain elusive. The presence of a well-characterized gene expression cascade shared by different Plasmodium falciparum strains could imply that transcriptional regulation and its natural variation do not contribute significantly to the evolution of parasite drug resistance. To clarify the role of transcriptional variation as a source of stain-specific diversity in the most deadly malaria species and to find genetic loci that dictate variations in gene expression, we examined genome-wide expression level polymorphisms (ELPs) in a genetic cross between phenotypically distinct parasite clones. Significant variation in gene expression is observed through direct co-hybridizations of RNA from different P. falciparum clones. Nearly 18% of genes were regulated by a significant expression quantitative trait locus. The genetic determinants of most of these ELPs resided in hotspots that are physically distant from their targets. The most prominent regulatory locus, influencing 269 transcripts, coincided with a Chromosome 5 amplification event carrying the drug resistance gene, pfmdr1, and 13 other genes. Drug selection pressure in the Dd2 parental clone lineage led not only to a copy number change in the pfmdr1 gene but also to an increased copy number of putative neighboring regulatory factors that, in turn, broadly influence the transcriptional network. Previously unrecognized transcriptional variation, controlled by polymorphic regulatory genes and possibly master regulators within large copy number variants, contributes to sweeping phenotypic evolution in drug-resistant malaria parasites.

Development of the malaria parasite, Plasmodium falciparum, in the blood is driven by a number of different genes expressed at different times and at different levels. Exactly what influences such transcriptional changes remains elusive, particularly in regard to important phenotypes like drug resistance. Using cDNA microarray hybridizations from the progeny of a Plasmodium genetic cross, we mapped gene expression quantitative trait loci (eQTLs) in an experimental population of malaria parasites. Each gene's transcript level was used as a segregating phenotype to identify regions of the Plasmodium genome dictating transcriptional variation. Several regulatory hotspots controlled the majority of gene expression variation, mostly via trans-acting mechanisms. One, influencing the largest number of transcripts, coincided with an amplified region of the genome traditionally associated with multiple drug resistance (MDR). Overall, integration of two functional genomic tools (gene mapping and transcript quantitation) has revealed a system-wide rewiring of the parasite transcription network: pleiotropic phenotypic variation, driven by drug selection on genome structure that may be attributed in large part to adaptive copy number polymorphisms in the parasite. These structural variants alter the expression of genes within the amplicon as well as many genes scattered across the Plasmodium genome.

Heritable levels of transcriptional variation, predominantly regulated by genomic copy number variants via trans mechanisms, are surprisingly abundant in drug-resistant malaria parasites.

Transcriptional control and gene silencing in Plasmodium falciparum

Infection with the apicomplexan parasite Plasmodium falciparum is associated with a high burden of morbidity and mortality across the developing world, yet the mechanisms of transcriptional control in this organism are poorly understood. While P. falciparum possesses many of the characteristics common to eukaryotic transcription, including much of the canonical machinery, it also demonstrates unique patterns of gene expression and possesses unusually AT-rich intergenic sequences. Importantly, several biological processes that are critical to parasite virulence involve highly regulated patterns of gene expression and silencing. The relative scarcity of transcription-associated proteins and specific cis-regulatory motifs recognized in the P. falciparum genome have been thought to reflect a reduced role for transcription factors in transcriptional control in these parasites. New approaches and technologies, however, have led to the discovery of many more of these elements, including an expanded family of DNA-binding proteins, and a re-assessment of this hypothesis is required. We review the current understanding of transcriptional control in P. falciparum, specifically highlighting promoter-driven and epigenetic mechanisms involved in the control of transcription initiation.

Discovering regulatory motifs in the Plasmodium genome using comparative genomics

MOTIVATION

Understanding gene regulation in Plasmodium, the causative agent of malaria, is an important step in deciphering its complex life cycle as well as leading to possible new targets for therapeutic applications. Very little is known about gene regulation in Plasmodium, and in particular, few regulatory elements have been identified. Such discovery has been significantly hampered by the high A-T content of some of the genomes of Plasmodium species, as well as the challenge in associating discovered regulatory elements to gene regulatory cascades due to Plasmodium's complex life cycle.

RESULTS

We report a new method of using comparative genomics to systematically discover motifs in Plasmodium without requiring any functional data. Different from previous methods, our method does not depend on sequence alignments, and thus is particularly suitable for highly divergent genomes. We applied our method to discovering regulatory motifs between the human parasite, P.falciparum, and its rodent-infectious relative, P.yoelii. We also tested our procedure against comparisons between P.falciparum and the primate-infectious, P.knowlesi. Our computational effort leads to an initial catalog of 38 distinct motifs, corresponding to over 16 200 sites in the Plasmodium genome. The functionality of these motifs was further supported by their defined distribution within the genome as well as a correlation with gene expression patterns. This initial map provides a systematic view of gene regulation in Plasmodium, which can be refined as additional genomes become available.

AVAILABILITY

The new algorithm, named motif discovery using orthologous sequences (MDOS), is available at http://www.ics.uci.edu/ approximately xhx/project/mdos/.

A computational model of gene expression reveals early transcriptional events at the subtelomeric regions of the malaria parasite, Plasmodium falciparum

A mathematical model of the intraerythrocytic developmental cycle identifies a delay between subtelomeric and central chromosomal gene activities in the malaria parasite, Plasmodium falciparum.

Background

The malaria parasite, Plasmodium falciparum, replicates asexually in a well-defined infection cycle within human erythrocytes (red blood cells). The intra-erythrocytic developmental cycle (IDC) proceeds with a 48 hour periodicity.

Results

Based on available malaria microarray data, which monitored gene expression over one complete IDC in one-hour time intervals, we built a mathematical model of the IDC using a circular variant of non-linear principal component analysis. This model enables us to identify rates of expression change within the data and reveals early transcriptional events at the subtelomeres of the parasite's nuclear chromosomes.

Conclusion

A delay between subtelomeric and central gene activities suggests that key events of the IDC are initiated at the subtelomeric regions of the P. falciparum nuclear chromosomes.

Distinct and stage specific nuclear factors regulate the expression of falcipains, Plasmodium falciparum cysteine proteases

Background

Plasmodium falciparum cysteine proteases (falcipains) play indispensable roles in parasite infection and development, especially in the process of host erythrocyte rupture/invasion and hemoglobin degradation. No detailed molecular analysis of transcriptional regulation of parasite proteases especially cysteine proteases has yet been reported. In this study, using a combination of transient transfection assays and electrophoretic mobility shift assays (EMSA), we demonstrate the presence of stage specific nuclear factors that bind to unique sequence elements in the 5'upstream regions of the falcipains and probably modulate the expression of cysteine proteases.

Results

Falcipains differ in their timing of expression and exhibit ability to compensate each other's functions at asexual blood stages of the parasite. Present study was undertaken to study the transcriptional regulation of falcipains. Transient transfection assay employing firefly luciferase as a reporter revealed that a ~1 kb sequence upstream of translational start site is sufficient for the functional transcriptional activity of falcipain-1 gene, while falcipain-2, -2' and -3 genes that exist within 12 kb stretch on chromosome 11 require ~2 kb upstream sequences for the expression of reporter luciferase activity. EMSA analysis elucidated binding of distinct nuclear factors to specific sequences within the 5'upstream regions of falcipain genes. Analysis of falcipains' 5'upstream regulatory regions did not reveal the presence of sequences known to bind general eukaryotic factors. However, we did find parasite specific sequence elements such as poly(dA) poly(dT) tracts, CCAAT boxes and a single 7 bp-G rich sequence, (A/G)NGGGG(C/A) in the 5' upstream regulatory regions of these genes, thereby suggesting the role(s) of Plasmodium specific transcriptional factors in the regulation of falcipain genes.

Conclusion

Taken together, these results suggest that expression of Plasmodium cysteine proteases is regulated at the transcriptional level and parasite specific factors regulate the expression of falcipain genes. These findings open new venues for further studies in identification of parasite specific transcription factors.

Transcriptional responses of Plasmodium falciparum to alpha-difluoromethylornithine-induced polyamine depletion

Abstract Polyamines are essential polycationic molecules involved in multiple cellular events, including cell differentiation, division and death. Inhibition of polyamine biosynthesis has been considered in diverse therapeutic strategies ranging from tumour suppressors to anti-parasitic agents. In the human malaria parasite, Plasmodium falciparum, inhibition of ornithine decarboxylase (ODC) results in the arrest of schizogony due to polyamine depletion. However, the exact physiological role of the polyamines in the parasite is unknown. Here, we present results of the depletion of polyamines in the malaria parasite by alpha-difluoromethylornithine inhibition of ODC, as observed with differential transcriptome profiling. Upon depletion of their endogenous polyamines, the up- and downregulated parasite transcripts were selected with suppression subtractive hybridisation and differences were detected using blots or DNA microarrays. A direct linkage between polyamine depletion and the differential expression of two distinct transcripts was observed, indicating the existence of a transcriptional feedback response in the P. falciparum transcriptome upon drug challenge. The data presented provide input into the role of the polyamines in the cellular biology of P. falciparum and contribute towards the validation of polyamine biosynthesis as an antimalarial target.

Exploring functional genomics for drug target and therapeutics discovery in Plasmodia

Functional genomics approaches are indispensable tools in the drug discovery arena and have recently attained increased attention in antibacterial drug discovery research. However, the application of functional genomics to post-genomics research of Plasmodia is still in comparatively early stages. Nonetheless, with this genus having the most species sequenced of any eukaryotic organism so far, the Plasmodia could provide unique opportunities for the study of intracellular eukaryotic pathogens. This review presents the status quo of functional genomics of the malaria parasite including descriptions of the transcriptome, proteome and interactome. We provide examples for the in silico mining of the X-ome data sets and illustrate how X-omic data from drug challenged parasites might be used in elucidating amongst others, the mode-of-action of inhibitory compounds, validate potential targets and discover novel targets/therapeutics.

Structural RNAs of known and unknown function identified in malaria parasites by comparative genomics and RNA analysis

As the genomes of more eukaryotic pathogens are sequenced, understanding how molecular differences between parasite and host might be exploited to provide new therapies has become a major focus. Central to cell function are RNA-containing complexes involved in gene expression, such as the ribosome, the spliceosome, snoRNAs, RNase P, and telomerase, among others. In this article we identify by comparative genomics and validate by RNA analysis numerous previously unknown structural RNAs encoded by the Plasmodium falciparum genome, including the telomerase RNA, U3, 31 snoRNAs, as well as previously predicted spliceosomal snRNAs, SRP RNA, MRP RNA, and RNAse P RNA. Furthermore, we identify six new RNA coding genes of unknown function. To investigate the relationships of the RNA coding genes to other genomic features in related parasites, we developed a genome browser for P. falciparum (http://areslab.ucsc.edu/cgi-bin/hgGateway). Additional experiments provide evidence supporting the prediction that snoRNAs guide methylation of a specific position on U4 snRNA, as well as predicting an snRNA promoter element particular to Plasmodium sp. These findings should allow detailed structural comparisons between the RNA components of the gene expression machinery of the parasite and its vertebrate hosts.

Sequences conserved by selection across mouse and human malaria species

Background

Little is known, either experimentally or computationally, about the genomic sequence features that regulate malaria genes. A sequence conservation analysis of the malaria species P. falciparum, P. berghei, P. yoelii, and P. chabaudi could significantly advance knowledge of malaria gene regulation.

Results

We computationally identify intergenic sequences conserved beyond neutral expectations, using a conservation algorithm that accounts for the strong compositional biases in malaria genomes. We first quantify the composition-specific divergence at silent positions in coding sequence. Using this as a background, we examine gene 5' regions, identifying 610 blocks conserved far beyond neutral expectations across the three mouse malariae, and 81 blocks conserved as strongly across all four species (p < 10^-6). Detailed analysis of these blocks indicates that only a minor fraction are likely to be previously unknown coding sequences. Analogous noncoding conserved blocks have been shown to regulate adjacent genes in other phylogenies, making the predicted blocks excellent candidates for novel regulatory functions. We also find three potential transcription factor binding motifs which exhibit strong conservation and overrepresentation among the rodent malariae.

Conclusion

A broader finding of our analysis is that less malaria intergenic sequence has been conserved by selection than in yeast or vertebrate genomes. This supports the hypothesis that transcriptional regulation is simpler in malaria than other eukaryotic species. We have built a public database containing all sequence alignments and functional predictions, and we expect this to be a valuable resource to the malaria research community.

An enhancer-like region regulates hrp3 promoter stage-specific gene expression in the human malaria parasite Plasmodium falciparum

The asexual blood stage of Plasmodium falciparum is comprised of morphologically distinct ring, trophozoite and schizont stages. Each of these developmental stages possesses a distinct pattern of gene expression. Regulation of P. falciparum gene expression is thought to occur, at least in part, at the promoter level. Previously, we have found that although the hrp3 mRNA is only seen in ring-stage parasites, deletion of a specific sequence in the 5' end of the promoter region decreased ring-stage expression of hrp3 and enabled detection of its transcripts in trophozoite-stage parasites. In order to investigate this stage specific regulation of gene expression, we employed a series of nested deletions of the 1.7-kb hrp3 promoter. Firefly luciferase gene was used as a reporter to evaluate the role of promoter sequences in gene regulation. Using this approach, we identified a ring-stage specific regulatory region on the hrp3 promoter located between -1.7 kb and -1.1 kb from the ATG initiation codon. Small 100-150 bp truncations on this enhancer-like region failed to uncover discrete regulatory sequences, suggesting the multipartite nature of this element. The data presented in this study demonstrate that stage specific promoter activity of the hrp3 gene in P. falciparum blood stage parasites is supported, at least in-part, by a small promoter region that can function in the absence of a larger chromosomal context.