- Open Access
Triaging informative cis-regulatory elements for the combinatorial control of temporal gene expression during Plasmodium falciparum intraerythrocytic development
Parasites & Vectors volume 8, Article number: 81 (2015)
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.
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.
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.
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.
The human malarial parasite Plasmodium falciparum adopts numerous morphologically distinct forms as it completes its complex life cycle in the human host and mosquito vector. As the parasite invades, colonises and multiplies within these diverse host environments a complex programme of developmentally-linked gene expression, utilising a diverse range of molecular mechanisms to exert control, has been described; for reviews see [1-3]. These are perhaps best exemplified during asexual intraerythrocytic development, where morphological transition from the newly invaded ring form progresses over a 48 hour period, through trophozoites and schizonts, to produce merozoites ready to reinitiate invasion in a new host erythrocyte. Over this 48 hr period, a well-defined cascade of peak mRNA steady-state accumulation has been described for some 50% of the parasite’s genome, with temporally- and functionally-linked clusters of genes being expressed in time to meet their biological demand [4-7].
With little apparent inter-strain variation in mRNA profiles during intraerythrocytic development, and minimal changes resulting from drug perturbations, this transcriptional cascade has been described as “hard-wired” [7-11]. Analyses of the molecular mechanisms that govern this developmentally-linked gene expression suggest that this “hard-wiring” is likely the result of globally-acting regulatory mechanisms, specifically; stage-specific variations in nucleosome positioning, processivity of the RNA polymerase II complex and stage-specific variations in the stability of the mRNA transcript [12-19]. Hypotheses that considered regulation of stage-specific gene expression exerted at the level of individual promoters, through specific transcription factor biding to cis-regulatory DNA motifs, fell out of favour in the early 2000′s due to the apparent absence of transcription factors in the P. falciparum genome [1,20,21]. In 2008, however, a restricted number of specific transcription factors, sharing the apetela 2 (AP2) DNA binding motif, were found in P. falciparum, with homologues quickly identified throughout all apicomplexans, leading to their designation as ApiAP2 transcription factors [22-25]. ApiAP2 have subsequently been shown to be critical regulators in the regulation of gene expression throughout the Plasmodium spp. life cycle as well as potentially playing a role in the monoallelic expression of the PfEMP1 virulence protein family through modulation of the local chromatin environment [26-31]. In 2010, using protein binding arrays, the cognate cis-acting DNA motif for 24 of the 27 P. falciparum ApiAP2 were determined . Interestingly, these DNA motifs are widely distributed within intergenic regions, with many intergenic regions sharing multiple ApiAP2 binding sites. Whilst this multiplicity of ApiAP2 binding sites may represent the means for a model of multifactorial control (a point that will be picked up later), whether all predicted DNA binding sites actually act as cis-regulatory sites remains to be addressed. In the absence of well-defined transcription start sites for P. falciparum, our inability to relate the position of a predicted ApiAP2 to this key transcriptional landmark hampers our efforts to design functional studies to explore their role in the control of transcription initiation.
In silico approaches have also been used to identify DNA motifs enriched within the flanking sequence of genes that share temporal peak mRNA profiles, function (utilising Gene Ontology terms) or share homologues in other Plasmodium spp [33-38]. Unfortunately, the catalogue of motifs predicted by each approach poorly overlap. Moreover, searches typically take an arbitrary length of flanking sequence for analysis. Our recent work exploring the size of flanking sequences in P. falciparum, highlight the challenge with such an arbitrary approach as we showed that the size of intergenic regions flanking a gene varies according to the nature of the transcriptional activity that takes place over this region . In this same study we also predict that transcription start sites lie further upstream of the start of the open reading frame than has previously been suggested and thus, key information may have been missed in these studies.
Recognising the challenge in defining DNA motifs that are most likely playing a role in the promoter-based control of transcription initiation in the absence of transcription start site data, we established a programme of work to; i) identify high-scoring DNA motifs that are repeatedly linked with genes that share the same temporal profile of peak mRNA accumulation and ii) undertake a search for potential new DNA motifs that lie further upstream from regions of intergenic sequences explored to date. To carry out this study we utilised the Finding Informative Regulatory Elements (FIRE) algorithm to explore correlations between DNA motifs located in intergenic sequences upstream of genes that share the same temporal profile of steady-state mRNA levels .
The source code and P. falciparum accessory files for the FIRE algorithm were obtained from the authors of the original FIRE study  and utilised on a PC operating a UNIX environment using the default sensitivity and stringency settings. These files are currently hosted, and freely available, online at https://tavazoielab.c2b2.columbia.edu/FIRE/. 5′ gene flanking sequences were obtained from a bespoke PERL script (intergenic.dist.2FASTA.pl available from https://sites.google.com/site/emesbioinformatics/group-software) using the P. falciparum General Feature Format (GFF) and genome sequence file downloaded from PlasmoDB5.5 (https://www.plasmoDB.org/plasmo). The intergenic.dist.2FASTA.pl program allows the user to specify the windows of 5′ flanking sequence to occur (-1000 to 0 and-1500 to-500 bp upstream of the start codon) and whether to capture sequences up to adjacent flanking genes if they fall within this window, or only when a full 1000 bp intergenic sequence can be captured. The FIRE output files for each search secured in separate folders. The FIRE motif heat maps and FIRE interaction heat maps resulting from the search of Groups A to D are attached in the Additional file 1. Analysis of the distribution of mutual information score(s) for the same motif discovered in one search (singleton) or multiple searches were performed using a Kruskall-Wallace one way analysis of variance with Dunn’s post-test (GraphPad Prism v5.1). WebLogos of DNA-binding specificities of all 27 members of the ApiAP2 protein family from P. falciparum along with their mRNA abundance profiles during intraerythrocytic development were sourced from the protein binding array study of Campbell et al. .
Results and discussion
Repeated discovery of DNA motifs associated with the temporal cascade of transcription during intraerythrocytic development
The Finding Informative Regulatory Elements (FIRE) algorithm discovers DNA motifs whose presence or absence in gene flanking sequences provides the most information about the expression profile of the associated flanking gene. For P. falciparum, peak mRNA accumulation data for some 2700 genes is available from a published temporal microarray study undertaken in 2 hr increments over the entire 48 hr of intraerythrocytic development . These data provide a continuous expression profile that can be used to discover overrepresented DNA motifs in the flanking intergenic sequences of genes that share the same temporal profile of peak mRNA accumulation. Using such an approach, FIRE has previously been used to discover 21 DNA motifs in a search of 1000 bp of 5′ flanking sequence in P. falciparum . We adopt here an approach that searches different permutations of a more recent annotation of P. falciparum gene 5′ flanking sequences to identify DNA motifs that are repeatedly discovered – thus offering an insight into their likelihood as cis-regulatory elements. Using our own recently published observations relating to the likely placement of transcription start sites between 600–1350 bp upstream of P. falciparum open reading frames , we also use an additional, but same sized, window to search further upstream than the original FIRE study to explore whether any potential new informative regulatory sites can be determined.
We elected to use search windows of 1000 bp. Not only did this allow a comparison to the original FIRE study, but we have also recently shown that this distance represents approximately half of the median size of an intergenic space that contains two promoter regions in P. falciparum . In total, four groups of 5′ flanking sequences (groups A to D) were secured for our analysis (Figure 1A). Group A most closely represents the sequences secured in the original FIRE report, i.e. 1000 bp of the most immediate flanking sequence. Based on our prediction that transcription start sites likely lie further upstream than considered in the original FIRE study, group C sequences were secured from a 1000 bp window located between 500 and 1500 bp upstream of each open reading frame. For both groups A and C, if the 1000 bp window overlapped with an adjacent open reading frame, the sequence captured was truncated to ensure only intergenic sequences were selected. Thus, two sets of sequences of up to 1000 bp for each gene were secured. Given our interest in repeatedly searching for the enrichment of the same DNA motif, two additional sets of upstream flanking sequences were secured. Whereas groups A and C captured up to 1000 bp of sequence, groups B and D secure the corresponding windows of 1000 bp sequence, respectively, but only when the entire 1000 bp sequence could be obtained. We hypothesised that those DNA motifs more likely associated with the control of stage-specific expression would be repeatedly identified in each of the groups, albeit with slightly different scores based on the different amount of sequences secured. Here, groups A to D consisted of 5579, 4300, 5297 and 3099 upstream flanking sequences, respectively.
FIRE analysis was performed on groups A to D, with 8–17 DNA motifs reported from each search. The algorithm produces a FIRE motif heat-map (see Figure 1B for example, see Additional file 1 for all files) for each search that provides a range of information for each DNA motif discovered. A colour-map is used to describe the correlation between either the over-representation (yellow) or under-representation (blue) of the DNA motif in genes that sharing the same peak mRNA accumulation profile (with the morphological stage representative of these timepoints indicated in Figure 1B). Correlations where this data is significantly over or under-represented (p < 0.05 after a Bonefori correction) are highlighted by bold red or blue surrounding lines, respectively. To the right of the heat map, the sequences of the seed motif for the search and the final optimized motif (as a WebLogo image) are shown alongside qualitative and quantitative evaluations of this DNA motif. Qualitatively, the location of the DNA motif in the 5′ flanking sequence is indicated for all motifs discovered here as well as any evidence of a positional or orientation bias following randomization trails. In the absence of well mapped transcription start sites in P. falciparum with which to correlate with these data, and the relatively few biases observed, no further analysis of these qualitative outcomes was performed here. The quantitative data reported includes; (i) mutual information, which indicates the extent of the association of the DNA motif with genes that share the same temporal profile of peak mRNA accumulation, (ii) the statistical significance (Z-score) of this association when compared to 10,000 randomizations of the input sequences, (iii) robustness of the association, i.e. how often the same motif is found in 10 separate jack-knife trials that remove one third of input sequences and (iv) the conservation index, indicates the shared presence of this motif in the flanking regions of orthologous genes in the murine malaria parasite P. yoelii (with indices of >0.95 considered significant).
Inspection of the lists of motifs identified in these searches reveal a total of 28 distinct DNA motifs (see Additional file 1). Of these, 14 had been previously described in the original FIRE study. The remaining 14 novel DNA motifs all share the common feature of each being discovered only once across groups A to D. A similar representation of singleton motif discovery in the original FIRE report can now be drawn by comparison to the searches performed here. Here, seven of the 21 motifs were not rediscovered in our analysis. Comparison of the mutual information scores between motifs discovered in two or more of the five groups (A to D and the original study) and those only discovered in a single search revealed a significantly lower score (one way analysis of variance with Dunns post-test, p < 0.05) in the singleton group. A second aspect of the search addressed whether searches for motifs in sequences located between 500 and 1500 bp upstream of the open reading frame would identify new motifs. Only four motifs were uniquely discovered in this region; all as singletons with low mutual information scores (0.027 to 0.031). Whilst it was hoped that this approach may have discovered additional motifs, it was recognised that the efficiency of the search algorithm in discovering motifs is dependent on the total sequences available for analysis. The use of windows located further upstream of the open reading frame will increase the likelihood over overlap with an adjacent open reading frame, thus limiting the total amount of sequences captured for such an analysis.
This outcome supports the approach adopted here in using repeated rediscovery of DNA motifs; those DNA motifs that are repeatedly discovered have a higher mutual information score. Taking the distribution of mutual information scores in the singleton DNA motifs (0.030 ± 0.004), a cut-off of 0.04 was established for the mutual information score for DNA motifs to be taken forward here. Thus, a short-list of 11 FIRE motifs (Fm1–11) was created, each motif being found in at least two of groups A to D as well as in the original FIRE analysis (Figure 2).
Fm1-11: a network of cis-acting motifs regulating ring-stage expression in P. falciparum?
To explore whether Fm1–11 represent likely cis-acting regulatory motifs, they were compared to consensus high affinity DNA binding motifs determined for the P. falciparum AP2 specific transcription factors [24,32]. Comparison of these AP2 DNA binding motifs against those of Fm1–11 revealed that six of these (Fm1, 2, 4, 8, 9 and 11) could be unambiguously attributed to a specific AP2 protein. Two further Fm (Fm5 and 6), sharing a degenerate CACA sequence, could not be attributed to a single AP2 transcription factor; instead a cluster of three AP2 transcription factors sharing affinity for these motifs were identified. Thus, of the 11 Fm identified here, eight appear to have a cognate specific transcription factor(s) available to bind them. Intriguingly, for the two Fm (Fm7 and 10) with the highest mutual information scores we could not identify a cognate AP trans-acting factor. These two motifs, therefore, may represent either cis-acting sites for non-AP2 transcription factors or other factors within the RNA polymerase II complex. Of note is that no other DNA motif identified from the analysis of groups A to D, or from the original FIRE study, had a reliably identifiable cognate AP2 binding partner.
Ranking Fm1–11 by the time during intraerythrocytic development their over-representation correlates with the peak of mRNA accumulation identifies an interesting common temporal property. Fm1–11 are overrepresented in the 5′ flanking sequence of genes that share a peak of mRNA accumulation within the first 24 hours of intraerythrocytic development – correlating with the ring and early trophozoite morphological stages. This contrasts with nuclear transcription run-on data that indicates that overall transcriptional activity during intraerythrocytic development is low during the first third of the cycle (ring stages) . This then increases gradually as the parasites mature in trophozoites and peaks in mid-schizont stages – some 12–16 hours past the latest timepoint linked with Fm1–11. Two additional global processes that also contribute to nucleic acid metabolism also appear to be at play later during intraerythrocytic development. First, mRNA half-life increase as intraerythrocytic development progresses; ranging from a mean of 9.5 minutes in ring stage parasites to 65 minutes in mature schizonts . Second, nucleosome occupancy over intergenic regions is most compact in ring and late-schizont stage parasites [12,13,15,16,19]. The lowest level of nucleosome occupancy is in the mature trophozoite stages, and presumably reflects an increased accessibility of the genomic DNA for nucleic acid metabolism, i.e. transcription and replication. These global mechanisms would appear to provide a reasonable explanation for the characteristic temporal transcriptional upregulation of hundreds of genes during the later stages of intraerythrocytic development. They do not, however, provide a clear insight into the temporal control of gene expression during stages. Our data, to this point, would lead us to suggest that Fm1–11 are cis-acting factors that play a role in the promoter-based control of genes expressed during the first 24 hours of intraerythocytic development.
Our last observation relevant to this evolving model of cis-trans promoter based control of ring-stage expression comes from a second output file from the FIRE algorithm - a motif interaction heat map. This colour-map illustrates any co-localization of the identified motifs within the same 5′ flanking region. Taking Fm1–11, we determined whether we could repeatedly find the same co-localization of these motifs in each of the searches we performed. Thus, by excluding co-localization of Fm1–11 that occur in only one search of groups A to D, we developed a qualitative network of interactions that were repeatedly discovered between Fm1–11 in two or more searches. This is illustrated in Figure 3 where the thickness of the line emphasizes how often the co-localization was discovered. The strongest link in the network was between Fm4 and Fm8, which was found in all four groups analysed. Interestingly, the AP2 transcription factors associated with Fm4 and Fm8 both share the same ring-stage profile of expression as do the flanking regions of genes in which these motifs are over-represented. This suggests that binding of AP2 to Fm4 and Fm8 may function as positive regulators in the upregulation of transcription of these genes within the ring-stage parasite. As a contrasting observation, expression of the cognate AP2 partner for motifs Fm1, 2, 5 and 6 (as determined from transcriptional and proteomic profiles) is actually upregulated in mature trophozoite stages. This observation could be rationalised if we consider AP2 binding to these Fm motifs acts as a negative regulator of gene expression. That is, a corollary of ring-stage specific expression is that these genes are not subject to global mechanisms that upregulate gene expression in the mature trophozoite stages – thus, AP2 binding to these Fm DNA motifs may act as an isolating negative regulator in maintaining the developmentally-linked expression pattern for these genes. Of note, however, is that whilst the potential for negative regulation though cis-trans promoter interactions has been suggested from promotor deletion studies, no direct demonstration for such a role for AP2 has been demonstrated thus far [40,41].
A second interesting feature of this interaction network is the triad of Fm7, 8 and 10. Found in searches of three of the four groups, these Fm represent the highest scoring motifs by mutual information score. Disappointingly, without a cognate AP2 binding partner for Fm 7 and 10, any further discussion relating to a role in directing stage-specific patterns of expression could not be made. However, as a whole, the evidence presented here for a network of colocalised cis-acting motifs, with cognate partners that may act in a positive and negative regulatory role, resonates with a model of combinatorial gene control originally proposed by van Noort and Huynen . They hypothesised that in the apparent absence of a large number of well-defined specific transcription factors in P. falciparum, that the necessary complexity necessary to drive the observed cascade of temporally-linked mRNA accumulation could be provided through the combination of a smaller number of transcription factors. The description of such a small number of specific transcription factors, the AP2 family, occurred subsequent to their report in 2006. Importantly, key elements of their model are indicated here, specifically; evidence for multiple cis-trans interactions within the same 5′ flanking region which may positively or negatively regulate promoter function. A refinement we suggest here is that this molecular mechanism would appear to be particularly important during the first 24 hours of intraerythrocytic development. Subsequent work that indicates multiple binding affinities for AP2 domains, the presence of multiple AP2 domains within a single protein and the potential for AP2 heterodimers suggest that there are additional layers of complexity to explore in these cis-trans interactions [23,32,43].
The occurrence and position of eight of the Fm motifs (Fm1, 2, 4–6, 8, 9 and 11) in all P. falciparum 5′ flanking regions has been previously mapped and reported (Additional file 1 for ). The frequency of Fm motif incidence per gene varies (0.2 to 196.7), but provides some one million occurrences in total over these intergenic regions. As such, it would appear that context, both in terms of position relative to a transcription start site and availability (accounted in part through nucleosome occupancy) to interact with a cognate AP2 partner is important in determining whether a mapped motif actually acts as a cis-acting regulatory site. Although transcription start sites have been bioinformatically predicted in P. falciparum, few sites have been confirmed experimentally . Recent work, using an improved directional, amplification free, RNAseq approach should shortly provide these key transcriptional landmarks (Chappell, Rayner and Berriman, pers comm.). Access of trans-acting factors to the DNA motifs is perhaps less on an issue, with P. falciparum intergenic regions being relatively depleted of nucleosomes compared to open reading frames [16,19]. Initial analysis of nucleosome binding over the predicted AP2 binding motifs suggests that some 65–97% of the total of all motifs are nucleosome free at some point during intraerythocytic development . With improved resolution stage-specific nucleosome occupancy maps now available [12,13], determining the temporal availability of Fm motifs, specifically those spatially organised around well mapped transcription start sites, offers an opportunity to test the hypothesis that specific transcription factor interactions with these motifs direct stage-specific transcription early during intraerythrocytic development.
Here we report a repeated bioinformatics search for over-represented DNA motifs within 5′ flanking intergenic regions of P. falciparum that we consider most likely to play a role in the stage-specific regulation of genes during intraerythrocytic development. Our search repeatedly identified 11 high scoring DNA motifs, and, significantly, we could identify a likely cognate AP2 trans-acting partner for 8 of these. Evidence of preference for regulation of mRNA accumulation during ring-stage development as well as an apparent interaction network between several of these motifs has led us to propose a nuanced modification to the combinatorial gene control model originally proposed by van Noort and Huynen . We propose that that cis-trans control of promoter function appears to offer a model for stage-specific expression during ring-stage development, complementing more global mechanisms regulating gene expression during the latter stages of intraerythrocytic development.
Apetela-2 domain transcription factor
Finding Informative Regulatory Elements
Deitsch K, Duraisingh M, Dzikowski R, Gunasekera A, Khan S, Le Roch K, et al. Mechanisms of gene regulation in Plasmodium. Am J Trop Med Hyg. 2007;77:201–8.
Horrocks P, Wong E, Russell K, Emes RD. Control of gene expression in Plasmodium falciparum - ten years on. Mol Biochem Parasitol. 2009;164:9–25.
Llinás M, Deitsch KW, Voss TS. Plasmodium gene regulation: far more to factor in. Trends Parasitol. 2008;24:551–6.
Bozdech Z, Llinás M, Pulliam BL, Wong ED, Zhu J, DeRisi JL. The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum. PLoS Biol. 2003;1:e5.
Le Roch KG, Johnson JR, Florens L, Zhou Y, Santrosyan A, Grainger M, et al. Global analysis of transcript and protein levels across the Plasmodium falciparum life cycle. Genome Res. 2004;14:2308–18.
Le Roch KG, Zhou YY, Blair PL, Grainger M, Moch JK, Haynes JD, et al. Discovery of gene function by expression profiling of the malaria parasite life cycle. Science. 2003;301:1503–8.
Llinas M, Bozdech Z, Wong ED, Adai AT, DeRisi JL. Comparative whole genome transcriptome analysis of three Plasmodium falciparum strains. Nuc Acids Res. 2006;34:1166–73.
Ganesan K, Ponmee N, Jiang L, Fowble JW, White J, Kamchonwongpaisan S, et al. A genetically hard-wired metabolic transcriptome in Plasmodium falciparum fails to mount protective responses to lethal antifolates. PLoS Pathog. 2008;4:e1000214.
Gunasekera AM, Myrick A, Le Roch K, Winzeler E, Wirth DF. Plasmodium falciparum: genome wide perturbations in transcript profiles among mixed stage cultures after chloroquine treatment. Exp Parasitol. 2007;117:87–92.
Hu G, Cabrera A, Kono M, Mok S, Chaal BK, Haase S, et al. Transcriptional profiling of growth perturbations of the human malaria parasite Plasmodium falciparum. Nat Biotechnol. 2010;28:91–8.
Natalang O, Bischoff E, Deplaine G, Proux C, Dillies MA, Sismeiro O, et al. Dynamic RNA profiling in Plasmodium falciparum synchronized blood stages exposed to lethal doses of artesunate. BMC Genomics. 2008;9:388.
Ay F, Bunnik EM, Varoquaux N, Bol SM, Prudhomme J, Vert JP, et al. Three-dimensional modeling of the P. falciparum genome during the erythrocytic cycle reveals a strong connection between genome architecture and gene expression. Genome Res. 2014;24:974–88.
Bunnik EM, Polishko A, Prudhomme J, Ponts N, Gill SS, Lonardi S, et al. DNA-encoded nucleosome occupancy is associated with transcription levels in the human malaria parasite Plasmodium falciparum. BMC Genomics. 2014;15:347–59.
Gopalakrishnan AM, Nyindodo LA, Ross Fergus M, Lopez-Estrano C. Plasmodium falciparum: preinitiation complex occupancy of active and inactive promoters during erythrocytic stage. Exp Parasitol. 2009;121:46–54.
Ponts N, Harris EY, Lonardi S, Le Roch KG. Nucleosome occupancy at transcription start sites in the human malaria parasite: a hard-wired evolution of virulence? Infect Genet Evol. 2011;11:716–24.
Ponts N, Harris EY, Prudhomme J, Wick I, Eckhardt-Ludka C, Hicks GR, et al. Nucleosome landscape and control of transcription in the human malaria parasite. Genome Res. 2010;20:228–38.
Shock JL, Fischer KF, DeRisi JL. Whole-genome analysis of mRNA decay in Plasmodium falciparum reveals a global lengthening of mRNA half-life during the intra-erythrocytic development cycle. Genome Biol. 2007;8:R134.
Sims JS, Militello KT, Sims PA, Patel VP, Kasper JM, Wirth DF. Patterns of gene-specific and total transcriptional activity during the Plasmodium falciparum intraerythrocytic developmental cycle. Eukaryot Cell. 2009;8:327–38.
Westenberger SJ, Cui L, Dharia N, Winzeler E, Cui L. Genome-wide nucleosome mapping of Plasmodium falciparum reveals histone-rich coding and histone-poor intergenic regions and chromatin remodeling of core and subtelomeric genes. BMC Genomics. 2009;10:610–21.
Aravind L, Iyer LM, Wellems TE, Miller LH. Plasmodium biology: genomic gleanings. Cell. 2003;115:771–85.
Coulson RMR, Hall N, Ouzounis CA. Comparative genomics of transcriptional control in the human malaria parasite Plasmodium falciparum. Genome Res. 2004;14:1548–54.
De Silva EK, Gehrke AR, Olszewski K, Leon I, Chahal JS, Bulyk ML, et al. Specific DNA-binding by Apicomplexan AP2 transcription factors. Proc Natl Acad Sci U S A. 2008;105:8393–8.
Lindner SE, De Silva EK, Keck JL, Llinas M. Structural determinants of DNA binding by a P. falciparum ApiAP2 transcriptional regulator. J Mol Biol. 2010;395:558–67.
Painter HJ, Campbell TL, Llinas M. The Apicomplexan AP2 family: integral factors regulating Plasmodium development. Mol Biochem Parasitol. 2011;176:1–7.
Balaji S, Babu MM, Iyer LM, Aravind L. Discovery of the principal specific transcription factors of apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains. Nucl Acids Res. 2005;33:3994–4006.
Flueck C, Bartfai R, Niederwieser I, Witmer K, Alako BT, Moes S, et al. A major role for the Plasmodium falciparum ApiAP2 protein PfSIP2 in chromosome end biology. PLoS Pathog. 2010;6:e1000784.
Iwanaga S, Kaneko I, Kato T, Yuda M. Identification of an AP2-family protein that is critical for malaria liver stage development. PLoS One. 2012;7:e47557.
Kafsack BF, Rovira-Graells N, Clark TG, Bancells C, Crowley VM, Campino SG, et al. A transcriptional switch underlies commitment to sexual development in malaria parasites. Nature. 2014;507:248–52.
Sinha A, Hughes KR, Modrzynska KK, Otto TD, Pfander C, Dickens NJ, et al. A cascade of DNA-binding proteins for sexual commitment and development in Plasmodium. Nature. 2014;507:253–7.
Yuda M, Iwanaga S, Shigenobu S, Kato T, Kaneko I. Transcription factor AP2-Sp and its target genes in malarial sporozoites. Mol Microbiol. 2010;75:854–63.
Yuda M, Iwanaga S, Shigenobu S, Mair GR, Janse CJ, Waters AP, et al. Identification of a transcription factor in the mosquito-invasive stage of malaria parasites. Mol Microbiol. 2009;71:1402–14.
Campbell TL, De Silva EK, Olszewski KL, Elemento O, Llinas M. Identification and genome-wide prediction of DNA binding specificities for the ApiAP2 family of regulators from the malaria parasite. PLoS Pathog. 2010;6:e1001165.
Elemento O, Slonim N, Tavazoie S. A universal framework for regulatory element discovery across all genomes and data types. Mol Cell. 2007;28:337–50.
Gunasekera AM, Myrick A, Militello KT, Sims JS, Dong CK, Gierahn T, et al. Regulatory motifs uncovered among gene expression clusters in Plasmodium falciparum. Mol Biochem Parasitol. 2007;153:19–30.
Jurgelenaite R, Dijkstra TM, Kocken CH, Heskes T. Gene regulation in the intraerythrocytic cycle of Plasmodium falciparum. Bioinformatics. 2009;25:1484–91.
Wu J, Sieglaff DH, Gervin J, Xie XS. Discovering regulatory motifs in the Plasmodium genome using comparative genomics. Bioinformatics. 2008;24:1843–9.
Young J, Johnson J, Benner C, Yan SF, Chen K, Le Roch K, et al. In silico discovery of transcription regulatory elements in Plasmodium falciparum. BMC Genomics. 2008;9:70.
Young JA, Fivelman QL, Blair PL, de la Vega P, Le Roch KG, Zhou YY, et al. The Plasmodium falciparum sexual development transcriptome: a microarray analysis using ontology-based pattern identification. Mol BiochemParasitol. 2005;143:67–79.
Russell K, Hasenkamp S, Emes R, Horrocks P. Analysis of the spatial and temporal arrangement of transcripts over intergenic regions in the human malarial parasite Plasmodium falciparum. BMC Genomics. 2013;14:267–77.
Horrocks P, Lanzer M. Mutational analysis identifies a five base pair cis-acting sequence essential for GBP130 promoter activity in Plasmodium falciparum. Mol Biochem Parasitol. 1999;99:77–87.
Porter ME. Positive and negative effects of deletions and mutations within the 5′ flanking sequences of Plasmodium falciparum DNA polymerase delta. Mol Biochem Parasitol. 2002;122:9–19.
van Noort V, Huynen MA. Combinatorial gene regulation in Plasmodium falciparum. Trends Genet. 2006;22:73–8.
Bougdor A, Braun L, Cannella D, Hakimi MA. Chromatin modifications: implications in the regulation of gene expression in Toxoplasma gondii. Cell Microbiol. 2010;12:413–23.
Brick K, Watanabe J, Pizzi E. Core promoters are predicted by their distinct physicochemical properties in the genome of Plasmodium falciparum. Genome Biol. 2008;9:R178.
This work was supported by a Biotechnology & Biological Sciences Research Council (BBSRC, BB/H002405/1) New Investigator Award to PH and BBSRC PhD award to KR.
The authors declare that they have no competing interests.
KR carried out the bioinformatics, interpreted the data described and drafted the manuscript. RE helped design the study and wrote the PERL scripts to secure P. falciparum sequence information. PH conceived and coordinated the study, participated in the interpretation of the data and drafted the manuscript. All authors read and approved the final version of the manuscript.
Additional file 1:
FIRE motif heat maps and FIRE interaction heat maps for groups A to D of P. falciparum 5′ gene flanking regions. Full lists of FIRE Weblogo motifs for groups A to D, as well as the original Elemento et al.  study are shown.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Russell, K., Emes, R. & Horrocks, P. Triaging informative cis-regulatory elements for the combinatorial control of temporal gene expression during Plasmodium falciparum intraerythrocytic development. Parasites Vectors 8, 81 (2015). https://doi.org/10.1186/s13071-015-0701-0
- AP2 transcription factor
- Cis-acting DNA motifs
- Combinatorial control
- Finding informative regulatory elements
- Stage-specific expression