Prediction of the general transcription factors associated with RNA polymerase II in Plasmodium falciparum: conserved features and differences relative to other eukaryotes

DNA sequence and chromatin differentiate sequence-specific transcription factor binding in the human malaria parasite Plasmodium falciparum

Development of the malaria parasite, Plasmodium falciparum, is regulated by a limited number of sequence-specific transcription factors (TFs). However, the mechanisms by which these TFs recognize genome-wide binding sites is largely unknown. To address TF specificity, we investigated the binding of two TF subsets that either bind CACACA or GTGCAC DNA sequence motifs and further characterized two additional ApiAP2 TFs, PfAP2-G and PfAP2-EXP, which bind unique DNA motifs (GTAC and TGCATGCA). We also interrogated the impact of DNA sequence and chromatin context on P. falciparum TF binding by integrating high-throughput in vitro and in vivo binding assays, DNA shape predictions, epigenetic post-translational modifications, and chromatin accessibility. We found that DNA sequence context minimally impacts binding site selection for paralogous CACACA-binding TFs, while chromatin accessibility, epigenetic patterns, co-factor recruitment, and dimerization correlate with differential binding. In contrast, GTGCAC-binding TFs prefer different DNA sequence context in addition to chromatin dynamics. Finally, we determined that TFs that preferentially bind divergent DNA motifs may bind overlapping genomic regions due to low-affinity binding to other sequence motifs. Our results demonstrate that TF binding site selection relies on a combination of DNA sequence and chromatin features, thereby contributing to the complexity of P. falciparum gene regulatory mechanisms.

Epigenetic Regulation and Chromatin Remodeling in Malaria Parasites

Plasmodium falciparum, the human malaria parasite, infects two hosts and various cell types, inducing distinct morphological and physiological changes in the parasite in response to different environmental conditions. These variations required the parasite to adapt and develop elaborate molecular mechanisms to ensure its spread and transmission. Recent findings have significantly improved our understanding of the regulation of gene expression in P. falciparum. Here, we provide an up-to-date overview of technologies used to highlight the transcriptomic adjustments occurring in the parasite throughout its life cycle. We also emphasize the complementary and complex epigenetic mechanisms regulating gene expression in malaria parasites. This review concludes with an outlook on the chromatin architecture, the remodeling systems, and how this 3D genome organization is critical in various biological processes.

Plasmodium falciparum Development from Gametocyte to Oocyst: Insight from Functional Studies

Malaria elimination may never succeed without the implementation of transmission-blocking strategies. The transmission of Plasmodium spp. parasites from the human host to the mosquito vector depends on circulating gametocytes in the peripheral blood of the vertebrate host. Once ingested by the mosquito during blood meals, these sexual forms undergo a series of radical morphological and metabolic changes to survive and progress from the gut to the salivary glands, where they will be waiting to be injected into the vertebrate host. The design of effective transmission-blocking strategies requires a thorough understanding of all the mechanisms that drive the development of gametocytes, gametes, sexual reproduction, and subsequent differentiation within the mosquito. The drastic changes in Plasmodium falciparum shape and function throughout its life cycle rely on the tight regulation of stage-specific gene expression. This review outlines the mechanisms involved in Plasmodium falciparum sexual stage development in both the human and mosquito vector, and zygote to oocyst differentiation. Functional studies unravel mechanisms employed by P. falciparum to orchestrate the expression of stage-specific functional products required to succeed in its complex life cycle, thus providing us with potential targets for developing new therapeutics. These mechanisms are based on studies conducted with various Plasmodium species, including predominantly P. falciparum and the rodent malaria parasites P. berghei. However, the great potential of epigenetics, genomics, transcriptomics, proteomics, and functional genetic studies to improve the understanding of malaria as a disease remains partly untapped because of limitations in studies using human malaria parasites and field isolates.

The role of long noncoding RNAs in malaria parasites

The human malaria parasites, including Plasmodium falciparum, persist as a major cause of global morbidity and mortality. The recent stalling of progress toward malaria elimination substantiates a need for novel interventions. Controlled gene expression is central to the parasite's numerous life cycle transformations and adaptation. With few specific transcription factors (TFs) identified, crucial roles for chromatin states and epigenetics in parasite transcription have become evident. Although many chromatin-modifying enzymes are known, less is known about which factors mediate their impacts on transcriptional variation. Like those of higher eukaryotes, long noncoding RNAs (lncRNAs) have recently been shown to have integral roles in parasite gene regulation. This review aims to summarize recent developments and key findings on the role of lncRNAs in P. falciparum.

A sequence-based foldability score combined with AlphaFold2 predictions to disentangle the protein order/disorder continuum

Order and disorder govern protein functions, but there is a great diversity in disorder, from regions that are-and stay-fully disordered to conditional order. This diversity is still difficult to decipher even though it is encoded in the amino acid sequences. Here, we developed an analytic Python package, named pyHCA, to estimate the foldability of a protein segment from the only information of its amino acid sequence and based on a measure of its density in regular secondary structures associated with hydrophobic clusters, as defined by the hydrophobic cluster analysis (HCA) approach. The tool was designed by optimizing the separation between foldable segments from databases of disorder (DisProt) and order (SCOPe [soluble domains] and OPM [transmembrane domains]). It allows to specify the ratio between order, embodied by regular secondary structures (either participating in the hydrophobic core of well-folded 3D structures or conditionally formed in intrinsically disordered regions) and disorder. We illustrated the relevance of pyHCA with several examples and applied it to the sequences of the proteomes of 21 species ranging from prokaryotes and archaea to unicellular and multicellular eukaryotes, for which structure models are provided in the AlphaFold protein structure database. Cases of low-confidence scores related to disorder were distinguished from those of sequences that we identified as foldable but are still excluded from accurate modeling by AlphaFold2 due to a lack of sequence homologs or to compositional biases. Overall, our approach is complementary to AlphaFold2, providing guides to map structural innovations through evolutionary processes, at proteome and gene scales.

Identification of Co-Existing Mutations and Gene Expression Trends Associated With K13-Mediated Artemisinin Resistance in Plasmodium falciparum

Plasmodium falciparum infects millions and kills thousands of people annually the world over. With the emergence of artemisinin and/or multidrug resistant strains of the pathogen, it has become even more challenging to control and eliminate the disease. Multiomics studies of the parasite have started to provide a glimpse into the confounding genetics and mechanisms of artemisinin resistance and identified mutations in Kelch13 (K13) as a molecular marker of resistance. Over the years, thousands of genomes and transcriptomes of artemisinin-resistant/sensitive isolates have been documented, supplementing the search for new genes/pathways to target artemisinin-resistant isolates. This meta-analysis seeks to recap the genetic landscape and the transcriptional deregulation that demarcate artemisinin resistance in the field. To explore the genetic territory of artemisinin resistance, we use genomic single-nucleotide polymorphism (SNP) datasets from 2,517 isolates from 15 countries from the MalariaGEN Network (The Pf3K project, pilot data release 4, 2015) to dissect the prevalence, geographical distribution, and co-existing patterns of genetic markers associated with/enabling artemisinin resistance. We have identified several mutations which co-exist with the established markers of artemisinin resistance. Interestingly, K13-resistant parasites harbor α-ß hydrolase and putative HECT domain-containing protein genes with the maximum number of SNPs. We have also explored the multiple, publicly available transcriptomic datasets to identify genes from key biological pathways whose consistent deregulation may be contributing to the biology of resistant parasites. Surprisingly, glycolytic and pentose phosphate pathways were consistently downregulated in artemisinin-resistant parasites. Thus, this meta-analysis highlights the genetic and transcriptomic features of resistant parasites to propel further exploratory studies in the community to tackle artemisinin resistance.

Sources of transcription variation in Plasmodium falciparum

Variation in transcript abundance can contribute to both short-term environmental response and long-term evolutionary adaptation. Most studies are designed to assess differences in mean transcription levels and do not consider other potentially important and confounding sources of transcriptional variation. Detailed quantification of variation sources will improve our ability to detect and identify the mechanisms that contribute to genome-wide transcription changes that underpin adaptive responses. To quantify innate levels of expression variation, we measured mRNA levels for more than 5000 genes in the malaria parasite, Plasmodium falciparum, among clones derived from two parasite strains across biologically and experimentally replicated batches. Using a mixed effects model, we partitioned the total variation among four sources - between strain, within strain, environmental batch effects, and stochastic noise. We found 646 genes with significant variation attributable to at least one of these sources. These genes were categorized by their predominant variation source and further examined using gene ontology enrichment analysis to associate function with each source of variation. Genes with environmental batch effect and within strain transcript variation may contribute to phenotypic plasticity, while genes with between strain variation may contribute to adaptive responses and processes that lead to parasite strain-specific survival under varied conditions.

A unique GCN5 histone acetyltransferase complex controls erythrocyte invasion and virulence in the malaria parasite Plasmodium falciparum

The histone acetyltransferase GCN5-associated SAGA complex is evolutionarily conserved from yeast to human and functions as a general transcription co-activator in global gene regulation. In this study, we identified a divergent GCN5 complex in Plasmodium falciparum, which contains two plant homeodomain (PHD) proteins (PfPHD1 and PfPHD2) and a plant apetela2 (AP2)-domain transcription factor (PfAP2-LT). To dissect the functions of the PfGCN5 complex, we generated parasite lines with either the bromodomain in PfGCN5 or the PHD domain in PfPHD1 deleted. The two deletion mutants closely phenocopied each other, exhibiting significantly reduced merozoite invasion of erythrocytes and elevated sexual conversion. These domain deletions caused dramatic decreases not only in histone H3K9 acetylation but also in H3K4 trimethylation, indicating synergistic crosstalk between the two euchromatin marks. Domain deletion in either PfGCN5 or PfPHD1 profoundly disturbed the global transcription pattern, causing altered expression of more than 60% of the genes. At the schizont stage, these domain deletions were linked to specific down-regulation of merozoite genes involved in erythrocyte invasion, many of which contain the AP2-LT binding motif and are also regulated by AP2-I and BDP1, suggesting targeted recruitment of the PfGCN5 complex to the invasion genes by these specific factors. Conversely, at the ring stage, PfGCN5 or PfPHD1 domain deletions disrupted the mutually exclusive expression pattern of the entire var gene family, which encodes the virulent factor PfEMP1. Correlation analysis between the chromatin state and alteration of gene expression demonstrated that up- and down-regulated genes in these mutants are highly correlated with the silent and active chromatin states in the wild-type parasite, respectively. Collectively, the PfGCN5 complex represents a novel HAT complex with a unique subunit composition including an AP2 transcription factor, which signifies a new paradigm for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

Epigenetic regulation of gene expression plays essential roles in orchestrating the general and parasite-specific cellular pathways in the malaria parasite Plasmodium falciparum. To better understand the epigenetic mechanisms in this parasite, we characterized the histone acetyltransferase GCN5-mediated transcription regulation during intraerythrocytic development of the parasite. Using tandem affinity purification and proteomic characterization, we identified that the PfGCN5-associated complex contains nine core components, including two PHD domain proteins (PfPHD1 and PfPHD2) and an AP2-domain transcription factor, which is divergent from the canonical GCN5 complexes evolutionarily conserved from yeast to human. To understand the functions of the PfGCN5 complex, we performed domain deletions in two subunits of this complex, PfGCN5 and PfPHD1. We found that the two deletion mutants displayed very similar growth phenotypes, including significantly reduced merozoite invasion rates and elevated sexual conversion. These two mutants were associated with dramatic decreases in histone H3K9 acetylation and H3K4 trimethylation, which led to global changes in chromatin states and gene expression. Consistent with the phenotypes, genes significantly affected by the PfGCN5 and PfPHD1 gene disruption include those participating in parasite-specific pathways such as invasion, virulence, and sexual development. In conclusion, this study presents a new model of the PfGCN5 complex for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

From Genes to Transcripts, a Tightly Regulated Journey in Plasmodium

Over the past decade, we have witnessed significant progresses in understanding gene regulation in Apicomplexa including the human malaria parasite, Plasmodium falciparum. This parasite possesses the ability to convert in multiple stages in various hosts, cell types, and environments. Recent findings indicate that P. falciparum is talented at using efficient and complementary molecular mechanisms to ensure a tight control of gene expression at each stage of its life cycle. Here, we review the current understanding on the contribution of the epigenome, atypical transcription factors, and chromatin organization to regulate stage conversion in P. falciparum. The adjustment of these regulatory mechanisms occurring during the progression of the life cycle will be extensively discussed.

Dynamic Chromatin Structure and Epigenetics Control the Fate of Malaria Parasites

Multiple hosts and various life cycle stages prompt the human malaria parasite, Plasmodium falciparum, to acquire sophisticated molecular mechanisms to ensure its survival, spread, and transmission to its next host. To face these environmental challenges, increasing evidence suggests that the parasite has developed complex and complementary layers of regulatory mechanisms controlling gene expression. Here, we discuss the recent developments in the discovery of molecular components that contribute to cell replication and differentiation and highlight the major contributions of epigenetics, transcription factors, and nuclear architecture in controlling gene regulation and life cycle progression in Plasmodium spp.

The highly diverse TATA box-binding proteins among protists: A review

Transcription is the first step of gene expression regulation and is a fundamental mechanism for establishing the viability and development of a cell. The TATA box-binding protein (TBP) interaction with a TATA box in a promoter is one of the best studied mechanisms in transcription initiation. TBP is a transcription factor that is highly conserved from archaea to humans and is essential for the transcription initiated by each of the three RNA polymerases. In addition, the discovery of TBP-related factor 1 (TRF1) and other factors related to TBP shed light on the variability among transcription initiation complexes, thus demonstrating that the compositions of these complexes are, in fact, more complicated than originally believed. Despite these facts, the majority of studies on transcription have been performed on animal, plant and fungal cells, which serve as canonical models, and information regarding protist cells is relatively scarce. The aim of this work is to review the diversity of the TBPs that have been documented in protists and describe some of the specific features that differentiate them from their counterparts in higher eukaryotes.

Epigenetic reader complexes of the human malaria parasite, Plasmodium falciparum

Epigenetic regulatory mechanisms are central to the development and survival of all eukaryotic organisms. These mechanisms critically depend on the marking of chromatin domains with distinctive histone tail modifications (PTMs) and their recognition by effector protein complexes. Here we used quantitative proteomic approaches to unveil interactions between PTMs and associated reader protein complexes of Plasmodium falciparum, a unicellular parasite causing malaria. Histone peptide pull-downs with the most prominent and/or parasite-specific PTMs revealed the binding preference for 14 putative and novel reader proteins. Amongst others, they highlighted the acetylation-level-dependent recruitment of the BDP1/BDP2 complex and identified an PhD-finger protein (PHD 1, PF3D7_1008100) that could mediate a cross-talk between H3K4me2/3 and H3K9ac marks. Tagging and interaction proteomics of 12 identified proteins unveiled the composition of 5 major epigenetic complexes, including the elusive TBP-associated-factor complex as well as two distinct GCN5/ADA2 complexes. Furthermore, it has highlighted a remarkable degree of interaction between these five (sub)complexes. Collectively, this study provides an extensive inventory of PTM-reader interactions and composition of epigenetic complexes. It will not only fuel further explorations of gene regulation amongst ancient eukaryotes, but also provides a stepping stone for exploration of PTM-reader interactions for antimalarial drug development.

Developments in drug design strategies for bromodomain protein inhibitors to target Plasmodium falciparum parasites

Introduction: Bromodomains (BRDs) bind to acetylated lysine residues, often on histones. The BRD proteins can contribute to gene regulation either directly through enzymatic activity or indirectly through recruitment of chromatin-modifying complexes or transcription factors. There is no evidence of direct orthologues of the Plasmodium falciparum BRD proteins (PfBDPs) outside the apicomplexans. PfBDPs are expressed during the parasite's life cycle in both the human host's blood and in the mosquito. PfBDPs could also prove to be promising targets for novel antimalarials, which are urgently required to address increasing drug resistance.Areas covered: This review discusses recent studies of the biology of PfBDPs, current target-based strategies for PfBDP inhibitor discovery, and different approaches to the important step of validating the specificity of hit compounds for PfBDPs.Expert opinion: The novelty of Plasmodium BRDs suggests that they could be targeted by selective compounds. Chemical series that showed promise in screens against human BRDs could be leveraged to create targeted compound libraries, as could hits from P. falciparum phenotypic screens. These targeted libraries and hits could be screened in target-based strategies aimed at discovery and optimization of novel inhibitors of PfBDPs. A key task for the field is to generate parasite assays to validate the hit compounds' specificity for PfBDPs.

The role of epigenetics and chromatin structure in transcriptional regulation in malaria parasites

Due to the unique selective pressures and extreme changes faced by the human malaria parasite Plasmodium falciparum throughout its life cycle, the parasite has evolved distinct features to alter its gene expression patterns. Along with classical gene regulation by transcription factors (TFs), of which only one family, the AP2 TFs, has been described in the parasite genome, a large body of evidence points toward chromatin structure and epigenetic factors mediating the changes in gene expression associated with parasite life cycle stages. These attributes may be critically important for immune evasion, host cell invasion and development of the parasite in its two hosts, the human and the Anopheles vector. Thus, the factors involved in the maintenance and regulation of chromatin and epigenetic features represent potential targets for antimalarial drugs. In this review, we discuss the mechanisms in P. falciparum that regulate chromatin structure, nucleosome landscape, the 3-dimensional structure of the genome and additional distinctive features created by parasite-specific genes and gene families. We review conserved traits of chromatin in eukaryotes in order to highlight what is unique in the parasite.

Transcriptional regulatory module analysis reveals that bridge proteins reconcile multiple signals in extracellular electron transfer pathways

Shewanella oneidensis MR-1 shows remarkable respiratory versatility with a large variety of extracellular electron acceptors (termed extracellular electron transfer, EET). To utilize the various electron acceptors, the bacterium must employ complex regulatory mechanisms to elicit the relevant EET pathways. To investigate the relevant mechanisms, we integrated EET genes and related transcriptional factors (TFs) into transcriptional regulatory modules (TRMs) and showed that many bridge proteins in these modules were signal proteins, which generally contained one or more signal processing domains (eg, GGDEF, EAL, PAS, etc.). Since Shewanella has to respond to diverse environmental conditions despite encoding few EET-relevant TFs, the overabundant signal proteins involved in the TRMs can help decipher the mechanism by which these microbes elicit a wide range of condition-specific responses. By combining proteomic data and protein bioinformatic analysis, we demonstrated that diverse signal proteins reconciled the different EET pathways, and we discussed the functional roles of signal proteins involved in the well-known MtrCAB pathway. Additionally, we showed that the signal proteins SO_2145 and SO_1417 played central roles in triggering EET pathways in anaerobic environments. Taken together, our results suggest that signal proteins have a profound impact on the transcriptional regulation of EET genes and thus have potential applications in microbial fuel cells.

Order in Disorder as Observed by the "Hydrophobic Cluster Analysis" of Protein Sequences

Hydrophobic cluster analysis (HCA) is an original approach for protein sequence analysis, which provides access to the foldable repertoire of the protein universe, including yet unannotated protein segments ("dark proteome"). Foldable segments correspond to ordered regions, as well as to intrinsically disordered regions (IDRs) undergoing disorder to order transitions. In this review, how HCA can be used to give insight into this last category of foldable segments is illustrated, with examples matching known 3D structures. After reviewing the HCA principles, examples of short foldable segments are given, which often contain short linear motifs, typically matching hydrophobic clusters. These segments become ordered upon contact with partners, with secondary structure preferences generally corresponding to those observed in the 3D structures within the complexes. Such small foldable segments are sometimes larger than the segments of known 3D structures, including flanking hydrophobic clusters that may be critical for interaction specificity or regulation, as well as intervening sequences allowing fuzziness. Cases of larger conditionally disordered domains are also presented, with lower density in hydrophobic clusters than well-folded globular domains or with exposed hydrophobic patches, which are stabilized by interaction with partners.

Unraveling the 3D genome of human malaria parasites

The chromosomes within the eukaryotic cell nucleus are highly dynamic and adopt complex hierarchical structures. Understanding how this three-dimensional (3D) nuclear architectureaffects gene regulation, cell cycle progression and disease pathogenesis are important biological questions in development and disease. Recently, many genome-wide technologies including chromosome conformation capture (3C) and 3C-based methodologies (4C, 5C, and Hi-C) have been developed to investigate 3D chromatin structure. In this review, we introduce 3D genome methodologies, with a focus on their application for understanding the nuclear architecture of the human malaria parasite, Plasmodium falciparum. An increasing amount of evidence now suggests that gene regulation in the parasite is largely regulated by epigenetic mechanisms and nuclear reorganization. Here, we explore the 3D genome architecture of P. falciparum, including local and global chromatin structure. In addition, molecular components important for maintaining 3D chromatin organization including architectural proteins and long non-coding RNAs are discussed. Collectively, these studies contribute to our understanding of how the plasticity of 3D genome architecture regulates gene expression and cell cycle progression in this deadly parasite.

Network-Based Methods for Identifying Key Active Proteins in the Extracellular Electron Transfer Process in Shewanella oneidensis MR-1

Shewanella oneidensis MR-1 can transfer electrons from the intracellular environment to the extracellular space of the cells to reduce the extracellular insoluble electron acceptors (Extracellular Electron Transfer, EET). Benefiting from this EET capability, Shewanella has been widely used in different areas, such as energy production, wastewater treatment, and bioremediation. Genome-wide proteomics data was used to determine the active proteins involved in activating the EET process. We identified 1012 proteins with decreased expression and 811 proteins with increased expression when the EET process changed from inactivation to activation. We then networked these proteins to construct the active protein networks, and identified the top 20 key active proteins by network centralization analysis, including metabolism- and energy-related proteins, signal and transcriptional regulatory proteins, translation-related proteins, and the EET-related proteins. We also constructed the integrated protein interaction and transcriptional regulatory networks for the active proteins, then found three exclusive active network motifs involved in activating the EET process-Bi-feedforward Loop, Regulatory Cascade with a Feedback, and Feedback with a Protein-Protein Interaction (PPI)-and identified the active proteins involved in these motifs. Both enrichment analysis and comparative analysis to the whole-genome data implicated the multiheme c-type cytochromes and multiple signal processing proteins involved in the process. Furthermore, the interactions of these motif-guided active proteins and the involved functional modules were discussed. Collectively, by using network-based methods, this work reported a proteome-wide search for the key active proteins that potentially activate the EET process.

Improving pairwise comparison of protein sequences with domain co-occurrence

Comparing and aligning protein sequences is an essential task in bioinformatics. More specifically, local alignment tools like BLAST are widely used for identifying conserved protein sub-sequences, which likely correspond to protein domains or functional motifs. However, to limit the number of false positives, these tools are used with stringent sequence-similarity thresholds and hence can miss several hits, especially for species that are phylogenetically distant from reference organisms. A solution to this problem is then to integrate additional contextual information to the procedure. Here, we propose to use domain co-occurrence to increase the sensitivity of pairwise sequence comparisons. Domain co-occurrence is a strong feature of proteins, since most protein domains tend to appear with a limited number of other domains on the same protein. We propose a method to take this information into account in a typical BLAST analysis and to construct new domain families on the basis of these results. We used Plasmodium falciparum as a case study to evaluate our method. The experimental findings showed an increase of 14% of the number of significant BLAST hits and an increase of 25% of the proteome area that can be covered with a domain. Our method identified 2240 new domains for which, in most cases, no model of the Pfam database could be linked. Moreover, our study of the quality of the new domains in terms of alignment and physicochemical properties show that they are close to that of standard Pfam domains. Source code of the proposed approach and supplementary data are available at: https://gite.lirmm.fr/menichelli/pairwise-comparison-with-cooccurrence

Deciphering the functions of the different proteins of an organism constitutes a first step toward the understanding of its biology. Because they provide strong clues regarding protein functions, domains occupy a key position among the relevant annotations that can be assigned to a protein. Protein domains are sequential motifs that are conserved along evolution and are found in different proteins and in different combinations. One common approach for identifying the domains of a protein is to run sequence-sequence comparisons with local alignment tools as BLAST. However these approaches sometimes miss several hits, especially for species that are phylogenetically distant from reference organisms. We propose here an approach to increase the sensitivity of pairwise sequence comparisons. This approach makes use of the fact that protein domains tend to appear with a limited number of other domains on the same protein (the domain co-occurrence property). On P. falciparum, our approach allows identifying 2240 new domains for which, in most cases, no domain of the Pfam database could be linked.

Genome-wide analysis of the human malaria parasite Plasmodium falciparum transcription factor PfNF-YB shows interaction with a CCAAT motif

Little is known about transcription factor regulation during the Plasmodium falciparum intraerythrocytic cycle. In order to elucidate the role of the P. falciparum (Pf)NF-YB transcription factor we searched for target genes in the entire genome. PfNF-YB mRNA is highly expressed in late trophozoite and schizont stages relative to the ring stage. In order to determine the candidate genes bound by PfNF-YB a ChIP-on-chip assay was carried out and 297 genes were identified. Ninety nine percent of PfNF-YB binding was to putative promoter regions of protein coding genes of which only 16% comprise proteins of known function. Interestingly, our data reveal that PfNF-YB binding is not exclusively to a canonical CCAAT box motif. PfNF-YB binds to genes coding for proteins implicated in a range of different biological functions, such as replication protein A large subunit (DNA replication), hypoxanthine phosphoribosyltransferase (nucleic acid metabolism) and multidrug resistance protein 2 (intracellular transport).

Genomics of apicomplexan parasites

The increasing prevalence of infections involving intracellular apicomplexan parasites such as Plasmodium, Toxoplasma, and Cryptosporidium (the causative agents of malaria, toxoplasmosis, and cryptosporidiosis, respectively) represent a significant global healthcare burden. Despite their significance, few treatments are available; a situation that is likely to deteriorate with the emergence of new resistant strains of parasites. To lay the foundation for programs of drug discovery and vaccine development, genome sequences for many of these organisms have been generated, together with large-scale expression and proteomic datasets. Comparative analyses of these datasets are beginning to identify the molecular innovations supporting both conserved processes mediating fundamental roles in parasite survival and persistence, as well as lineage-specific adaptations associated with divergent life-cycle strategies. The challenge is how best to exploit these data to derive insights into parasite virulence and identify those genes representing the most amenable targets. In this review, we outline genomic datasets currently available for apicomplexans and discuss biological insights that have emerged as a consequence of their analysis. Of particular interest are systems-based resources, focusing on areas of metabolism and host invasion that are opening up opportunities for discovering new therapeutic targets.

The Role of Chromatin Structure in Gene Regulation of the Human Malaria Parasite

The human malaria parasite, Plasmodium falciparum, depends on a coordinated regulation of gene expression for development and propagation within the human host. Recent developments suggest that gene regulation in the parasite is largely controlled by epigenetic mechanisms. Here, we discuss recent advancements contributing to our understanding of the mechanisms controlling gene regulation in the parasite, including nucleosome landscape, histone modifications, and nuclear architecture. In addition, various processes involved in regulation of parasite-specific genes and gene families are examined. Finally, we address the use of epigenetic processes as targets for novel antimalarial therapies. Collectively, these topics highlight the unique biology of P. falciparum, and contribute to our understanding of mechanisms regulating gene expression in this deadly parasite.

Bromodomains in Protozoan Parasites: Evolution, Function, and Opportunities for Drug Development

Parasitic infections remain one of the most pressing global health concerns of our day, affecting billions of people and producing unsustainable economic burdens. The rise of drug-resistant parasites has created an urgent need to study their biology in hopes of uncovering new potential drug targets. It has been established that disrupting gene expression by interfering with lysine acetylation is detrimental to survival of apicomplexan (Toxoplasma gondii and Plasmodium spp.) and kinetoplastid (Leishmania spp. and Trypanosoma spp.) parasites. As "readers" of lysine acetylation, bromodomain proteins have emerged as key gene expression regulators and a promising new class of drug target. Here we review recent studies that demonstrate the essential roles played by bromodomain-containing proteins in parasite viability, invasion, and stage switching and present work showing the efficacy of bromodomain inhibitors as novel antiparasitic agents. In addition, we performed a phylogenetic analysis of bromodomain proteins in representative pathogens, some of which possess unique features that may be specific to parasite processes and useful in future drug development.

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.

Improvement in Protein Domain Identification Is Reached by Breaking Consensus, with the Agreement of Many Profiles and Domain Co-occurrence

Traditional protein annotation methods describe known domains with probabilistic models representing consensus among homologous domain sequences. However, when relevant signals become too weak to be identified by a global consensus, attempts for annotation fail. Here we address the fundamental question of domain identification for highly divergent proteins. By using high performance computing, we demonstrate that the limits of state-of-the-art annotation methods can be bypassed. We design a new strategy based on the observation that many structural and functional protein constraints are not globally conserved through all species but might be locally conserved in separate clades. We propose a novel exploitation of the large amount of data available: 1. for each known protein domain, several probabilistic clade-centered models are constructed from a large and differentiated panel of homologous sequences, 2. a decision-making protocol combines outcomes obtained from multiple models, 3. a multi-criteria optimization algorithm finds the most likely protein architecture. The method is evaluated for domain and architecture prediction over several datasets and statistical testing hypotheses. Its performance is compared against HMMScan and HHblits, two widely used search methods based on sequence-profile and profile-profile comparison. Due to their closeness to actual protein sequences, clade-centered models are shown to be more specific and functionally predictive than the broadly used consensus models. Based on them, we improved annotation of Plasmodium falciparum protein sequences on a scale not previously possible. We successfully predict at least one domain for 72% of P. falciparum proteins against 63% achieved previously, corresponding to 30% of improvement over the total number of Pfam domain predictions on the whole genome. The method is applicable to any genome and opens new avenues to tackle evolutionary questions such as the reconstruction of ancient domain duplications, the reconstruction of the history of protein architectures, and the estimation of protein domain age. Website and software: http://www.lcqb.upmc.fr/CLADE.

Current sequence databases contain hundreds of billions of nucleotides coding for genes and a classification of these sequences is a primary problem in genomics. A reasonable way to organize these sequences is through their predicted domains, but the identification of domains in very divergent sequences, spanning the entire phylogenetic tree of species, is a difficult problem. By generating multiple probabilistic models for a domain, describing the spread of evolutionary patterns in different phylogenetic clades, we can effectively explore domains that are likely to be coded in gene sequences. Through a machine learning approach and optimization techniques, coding for expected evolutionary constraints, we filter the many possibilities of domain identification found for a gene and propose the most likely domain architecture associated to it. The application of this novel approach to the full genome of Plasmodium falciparum, to a dataset of sequences from three SCOP datasets highlights the interest of exploring multiple pathways of domain evolution in the aim of extracting biological information from genomic sequences. Our new computational approach was developed with the hope of providing a novel tier of accurate and precise tools that complement existing tools such as HMMer, HHblits and PSI-BLAST, by exploring in a novel way the large amount of sequence data available. The existence of powerful databases for sequences, domains and architectures help make this hope a reality.

Landscape and Dynamics of Transcription Initiation in the Malaria Parasite Plasmodium falciparum

A comprehensive map of transcription start sites (TSSs) across the highly AT-rich genome of P. falciparum would aid progress toward deciphering the molecular mechanisms that underlie the timely regulation of gene expression in this malaria parasite. Using high-throughput sequencing technologies, we generated a comprehensive atlas of transcription initiation events at single-nucleotide resolution during the parasite intra-erythrocytic developmental cycle. This detailed analysis of TSS usage enabled us to define architectural features of plasmodial promoters. We demonstrate that TSS selection and strength are constrained by local nucleotide composition. Furthermore, we provide evidence for coordinate and stage-specific TSS usage from distinct sites within the same transcription unit, thereby producing transcript isoforms, a subset of which are developmentally regulated. This work offers a framework for further investigations into the interactions between genomic sequences and regulatory factors governing the complex transcriptional program of this major human pathogen.

Plasmodium falciparum XPD translocates in 5' to 3' direction, is expressed throughout the blood stages, and interacts with p44

XPD helicase, a TFIIH subunit, is essential for several processes including transcription, NER, cell cycle regulation, and apoptosis in eukaryotes. Another component of TFIIH, namely p44, is among the well-known interacting partners of XPD and is vital in regulating the helicase activities of latter. However, none of the above mentioned proteins have been functionally characterized in Plasmodium falciparum. Consequently, in this study, we performed detailed studies on XPD and its interacting partner, p44, from P. falciparum 3D7 strain. Accordingly, we expressed and purified recombinant PfXPD and its fragments and Pfp44 proteins and characterized the enzymatic activities of PfXPD and its fragments. The in vivo stage-specific expression and subcellular localizations of PfXPD and Pfp44 proteins were studied using the specific antibodies in the intraerythrocytic developmental stages of P. falciparum 3D7 strain. Our results suggest that PfXPD displays the characteristic ssDNA-dependent ATPase and 5'-3' DNA helicase activities. We also report the existence of two high molecular weight forms of p44 in P. falciparum 3D7 strain. Both PfXPD and Pfp44 colocalize in the nucleus and interact with each other, which suggest that they are most likely components of the same complex apparently, TFIIH. Furthermore, during trophozoite and schizont stages, both proteins exhibit a distinct cytoplasmic distribution pattern which implies that PfXPD and Pfp44 might also be involved in other functions. These studies will aid in understanding the basic biology of malaria parasite.

Comparative insight into nucleotide excision repair components of Plasmodium falciparum

Nucleotide excision repair (NER) is one of the DNA repair pathways crucial for maintenance of genome integrity and deals with repair of DNA damages arising due to exogenous and endogenous factors. The multi-protein transcription initiation factor TFIIH plays a critical role in NER and transcription and is highly conserved throughout evolution. The malaria parasite Plasmodium falciparum has been a challenge for the researchers for a long time because of emergence of drug resistance. The availability of its genome sequence has opened new avenues for research. Antimalarial drugs like chloroquine and mefloquine have been reported to inhibit NER pathway mediated repair reactions and thus promote mutagenesis. Previous studies have validated existence and implied possible association of defective or altered DNA repair pathways with development of drug resistant phenotype in certain P. falciparum strains. We conjecture that a compromised NER pathway in combination with other DNA repair pathways might be conducive for the emergence and sustenance of drug resistance in P. falciparum. Therefore we decided to unravel the components of NER pathway in P. falciparum and using bioinformatics based approaches here we report a genome wide in silico analysis of NER components from P. falciparum and their comparison with the human host. Our results reveal that P. falciparum genome contains almost all the components of NER but we were unable to find clear homologue for p62 and XPC in its genome. The structure modeling of all the components further suggests that their structures are significantly conserved. Furthermore this study lays a foundation to perform similar comparative studies between drug resistant and drug sensitive strains of parasite in order to understand DNA repair-related mechanisms of drug resistance.

Development of diaminoquinazoline histone lysine methyltransferase inhibitors as potent blood-stage antimalarial compounds

Modulating epigenetic mechanisms in malarial parasites is an emerging avenue for the discovery of novel antimalarial drugs. Previously we demonstrated the potent in vitro and in vivo antimalarial activity of (1-benzyl-4-piperidyl)[6,7-dimethoxy-2-(4-methyl-1,4-diazepin-1-yl)-4-quinazolinyl]amine (BIX01294; 1), a known human G9a inhibitor, together with its dose-dependent effects on histone methylation in the malarial parasite. This work describes our initial medicinal chemistry efforts to optimise the diaminoquinazoline chemotype for antimalarial activity. A variety of analogues were designed by substituting the 2 and 4 positions of the quinazoline core, and these molecules were tested against Plasmodium falciparum (3D7 strain). Several analogues with IC50 values as low as 18.5 nM and with low mammalian cell toxicity (HepG2) were identified. Certain pharmacophoric features required for antimalarial activity were found to be analogous to the previously published SAR of these analogues for G9a inhibition, thereby suggesting potential similarities between the malarial and human HKMT targets of this chemotype. Physiochemical, in vitro activity, and in vitro metabolism studies were also performed for a select set of potent analogues to evaluate their potential as antimalarial leads.

Nuclear glycolytic enzyme enolase of Toxoplasma gondii functions as a transcriptional regulator

Apicomplexan parasites including Toxoplasma gondii have complex life cycles within different hosts and their infectivity relies on their capacity to regulate gene expression. However, little is known about the nuclear factors that regulate gene expression in these pathogens. Here, we report that T. gondii enolase TgENO2 is targeted to the nucleus of actively replicating parasites, where it specifically binds to nuclear chromatin in vivo. Using a ChIP-Seq technique, we provide evidence for TgENO2 enrichment at the 5′ untranslated gene regions containing the putative promoters of 241 nuclear genes. Ectopic expression of HA-tagged TgENO1 or TgENO2 led to changes in transcript levels of numerous gene targets. Targeted disruption of TgENO1 gene results in a decrease in brain cyst burden of chronically infected mice and in changes in transcript levels of several nuclear genes. Complementation of this knockout mutant with ectopic TgENO1-HA fully restored normal transcript levels. Our findings reveal that enolase functions extend beyond glycolytic activity and include a direct role in coordinating gene regulation in T. gondii.

The Babesia bovis gene and promoter model: an update from full-length EST analysis

BACKGROUND

Babesia bovis is an apicomplexan parasite that causes babesiosis in infected cattle. Genomes of pathogens contain promising information that can facilitate the development of methods for controlling infections. Although the genome of B. bovis is publically available, annotated gene models are not highly reliable prior to experimental validation. Therefore, we validated a preproposed gene model of B. bovis and extended the associated annotations on the basis of experimentally obtained full-length expressed sequence tags (ESTs).

RESULTS

From in vitro cultured merozoites, 12,286 clones harboring full-length cDNAs were sequenced from both ends using the Sanger method, and 6,787 full-length cDNAs were assembled. These were then clustered, and a nonredundant referential data set of 2,115 full-length cDNA sequences was constructed. The comparison of the preproposed gene model with our data set identified 310 identical genes, 342 almost identical genes, 1,054 genes with potential structural inconsistencies, and 409 novel genes. The median length of 5' untranslated regions (UTRs) was 152 nt. Subsequently, we identified 4,086 transcription start sites (TSSs) and 2,023 transcriptionally active regions (TARs) by examining 5' ESTs. We identified ATGGGG and CCCCAT sites as consensus motifs in TARs that were distributed around -50 bp from TSSs. In addition, we found ACACA, TGTGT, and TATAT sites, which were distributed periodically around TSSs in cycles of approximately 150 bp. Moreover, related periodical distributions were not observed in mammalian promoter regions.

CONCLUSIONS

The observations in this study indicate the utility of integrated bioinformatics and experimental data for improving genome annotations. In particular, full-length cDNAs with one-base resolution for TSSs enabled the identification of consensus motifs in promoter sequences and demonstrated clear distributions of identified motifs. These observations allowed the illustration of a model promoter composition, which supports the differences in transcriptional regulation frameworks between apicomplexan parasites and mammals.

Comprehensive repertoire of foldable regions within whole genomes

In order to get a comprehensive repertoire of foldable domains within whole proteomes, including orphan domains, we developed a novel procedure, called SEG-HCA. From only the information of a single amino acid sequence, SEG-HCA automatically delineates segments possessing high densities in hydrophobic clusters, as defined by Hydrophobic Cluster Analysis (HCA). These hydrophobic clusters mainly correspond to regular secondary structures, which together form structured or foldable regions. Genome-wide analyses revealed that SEG-HCA is opposite of disorder predictors, both addressing distinct structural states. Interestingly, there is however an overlap between the two predictions, including small segments of disordered sequences, which undergo coupled folding and binding. SEG-HCA thus gives access to these specific domains, which are generally poorly represented in domain databases. Comparison of the whole set of SEG-HCA predictions with the Conserved Domain Database (CDD) also highlighted a wide proportion of predicted large (length >50 amino acids) segments, which are CDD orphan. These orphan sequences may either correspond to highly divergent members of already known families or belong to new families of domains. Their comprehensive description thus opens new avenues to investigate new functional and/or structural features, which remained so far uncovered. Altogether, the data described here provide new insights into the protein architecture and organization throughout the three kingdoms of life.

Spontaneous or induced folding into a specific 3D structure is a key property of proteins to perform their biological functions. Folded 3D structures of proteins perform specific functions, including interactions with other proteins. Intrinsically disordered regions also mediate interaction, gaining structure only when bound to a target protein. In both cases, hydrophobicity generally plays a major role in the protein segment “foldability”. Here, we developed an original procedure to identify foldable segments from only the information of a single amino acid sequence and to explore protein structures at a proteomic scale. Our approach goes beyond the simple consideration of mean hydrophobicity, by including the secondary structure information through the use of a two-dimensional transposition of the sequence. The developed procedure, combined with disorder predictors, may facilitate the specific identification of small segments that undergo coupled folding and binding. Combined with the analysis of specific domain databases, it also highlights orphan foldable segments, which remain yet uncharacterized.

Identification of hidden relationships from the coupling of hydrophobic cluster analysis and domain architecture information

MOTIVATION

Describing domain architecture is a critical step in the functional characterization of proteins. However, some orphan domains do not match any profile stored in dedicated domain databases and are thereby difficult to analyze.

RESULTS

We present here an original novel approach, called TREMOLO-HCA, for the analysis of orphan domain sequences and inspired from our experience in the use of Hydrophobic Cluster Analysis (HCA). Hidden relationships between protein sequences can be more easily identified from the PSI-BLAST results, using information on domain architecture, HCA plots and the conservation degree of amino acids that may participate in the protein core. This can lead to reveal remote relationships with known families of domains, as illustrated here with the identification of a hidden Tudor tandem in the human BAHCC1 protein and a hidden ET domain in the Saccharomyces cerevisiae Taf14p and human AF9 proteins. The results obtained in such a way are consistent with those provided by HHPRED, based on pairwise comparisons of HHMs. Our approach can, however, be applied even in absence of domain profiles or known 3D structures for the identification of novel families of domains. It can also be used in a reverse way for refining domain profiles, by starting from known protein domain families and identifying highly divergent members, hitherto considered as orphan.

AVAILABILITY

We provide a possible integration of this approach in an open TREMOLO-HCA package, which is fully implemented in python v2.7 and is available on request. Instructions are available at http://www.impmc.upmc.fr/∼callebau/tremolohca.html.

CONTACT

isabelle.callebaut@impmc.upmc.fr

SUPPLEMENTARY INFORMATION

Supplementary Data are available at Bioinformatics online.

Transcription sites are developmentally regulated during the asexual cycle of Plasmodium falciparum

Increasing evidence shows that the spatial organization of transcription is an important epigenetic factor in eukaryotic gene regulation. The malaria parasite Plasmodium falciparum shows a remarkably complex pattern of gene expression during the erythrocytic cycle, paradoxically contrasting with the relatively low number of putative transcription factors encoded by its genome. The spatial organization of nuclear subcompartments has been correlated with the regulation of virulence genes. Here, we investigate the nuclear architecture of transcription during the asexual cycle of malaria parasites. As in mammals, transcription is organized into discrete nucleoplasmic sites in P. falciparum, but in a strikingly lower number of foci. An automated analysis of 3D images shows that the number and intensity of transcription sites vary significantly between rings and trophozoites, although the nuclear volume remains constant. Transcription sites are spatially reorganized during the asexual cycle, with a higher proportion of foci located in the outermost nuclear region in rings, whereas in trophozoites, foci are evenly distributed throughout the nucleoplasm. As in higher eukaryotes, transcription sites are predominantly found in areas of low chromatin density. Immunofluorescence analysis shows that transcription sites form an exclusive nuclear compartment, different from the compartments defined by the silenced or active chromatin markers. In conclusion, these data suggest that transcription is spatially contained in discrete foci that are developmentally regulated during the asexual cycle of malaria parasites and located in areas of low chromatin density.

A Tetracycline-Repressible Transactivator System to Study Essential Genes in Malaria Parasites

A major obstacle in analyzing gene function in apicomplexan parasites is the absence of a practical regulatable expression system. Here, we identified functional transcriptional activation domains within Apicomplexan AP2 (ApiAP2) family transcription factors. These ApiAP2 transactivation domains were validated in blood-, liver-, and mosquito-stage parasites and used to create a robust conditional expression system for stage-specific, tetracycline-dependent gene regulation in Toxoplasma gondii, Plasmodium berghei, and Plasmodium falciparum. To demonstrate the utility of this system, we created conditional knockdowns of two essential P. berghei genes: profilin (PRF), a protein implicated in parasite invasion, and N-myristoyltransferase (NMT), which catalyzes protein acylation. Tetracycline-induced repression of PRF and NMT expression resulted in a dramatic reduction in parasite viability. This efficient regulatable system will allow for the functional characterization of essential proteins that are found in these important parasites.

Functional insights into the core-TFIIH from a comparative survey

TFIIH is a eukaryotic complex composed of two subcomplexes, the CAK (Cdk activating kinase) and the core-TFIIH. The core-TFIIH, composed of seven subunits (XPB, XPD, P62, P52, P44, P34, and P8), plays a crucial role in transcription and repair. Here, we performed an extended sequence analysis to establish the accurate phylogenetic distribution of the core-TFIIH in 63 eukaryotic organisms. In spite of the high conservation of the seven subunits at the sequence and genomic levels, the non-enzymatic P8, P34, P52 and P62 are absent from one or a few unicellular species. To gain insight into their respective roles, we undertook a comparative genomic analysis of the whole proteome to identify the gene sets sharing similar presence/absence patterns. While little information was inferred for P8 and P62, our studies confirm the known role of P52 in repair and suggest for the first time the implication of the core TFIIH in mRNA splicing via P34.

Plasmodium falciparum responds to amino acid starvation by entering into a hibernatory state

The human malaria parasite Plasmodium falciparum is auxotrophic for most amino acids. Its amino acid needs are met largely through the degradation of host erythrocyte hemoglobin; however the parasite must acquire isoleucine exogenously, because this amino acid is not present in adult human hemoglobin. We report that when isoleucine is withdrawn from the culture medium of intraerythrocytic P. falciparum, the parasite slows its metabolism and progresses through its developmental cycle at a reduced rate. Isoleucine-starved parasites remain viable for 72 h and resume rapid growth upon resupplementation. Protein degradation during starvation is important for maintenance of this hibernatory state. Microarray analysis of starved parasites revealed a 60% decrease in the rate of progression through the normal transcriptional program but no other apparent stress response. Plasmodium parasites do not possess a TOR nutrient-sensing pathway and have only a rudimentary amino acid starvation-sensing eukaryotic initiation factor 2α (eIF2α) stress response. Isoleucine deprivation results in GCN2-mediated phosphorylation of eIF2α, but kinase-knockout clones still are able to hibernate and recover, indicating that this pathway does not directly promote survival during isoleucine starvation. We conclude that P. falciparum, in the absence of canonical eukaryotic nutrient stress-response pathways, can cope with an inconsistent bloodstream amino acid supply by hibernating and waiting for more nutrient to be provided.

Pulmonary CCL18 recruits human regulatory T cells

CCL18 is both a constitutively expressed and an inducible chemokine, whose role in the inflammatory reaction is poorly known. The aim of this study was to evaluate whether CCL18 has the capacity to attract human T cells with a regulatory function (regulatory T cells [Treg]). Results from chemotaxis assays performed on different types of Treg showed that CD4(+)CD25(+)CD127(low) cells, but neither T regulatory type 1 clones nor Treg differentiated in vitro with anti-CD3/CD46 mAbs, were recruited by CCL18 in a dose-dependent manner. CCL18-recruited memory CD4(+) T cells were enriched in CD25(high), CD25(+)CD127(low), latency-associated peptide/TGF-β1, and CCR4-expressing T cells, whereas there was no enrichment in Foxp3(+) cells as compared with controls. Stimulated CCL18-recruited memory T cells produced significantly increased amounts of the regulatory cytokines IL-10 and TGF-β1, as well as IL-4, but not IFN-γ and IL-17. Cell surface CCL18 binding was found predominantly on IL-10(+) (26.3 ± 5.8%) and on a few latency-associated peptide/TGF-β1(+) (18.1 ± 1.9%) and IL-4(+) (14.5 ± 2.9%) memory T cells. In an in vivo model of SCID mice grafted with human skin and reconstituted with autologous PBMCs, the intradermal injection of CCL18 led to the cutaneous recruitment of CD4(+), CD25(+), and IL-10(+) cells, but not Foxp3(+) cells. Furthermore, CCL18-recruited memory T cells inhibited the proliferation of CD4(+)CD25(-) effector T cells through an IL-10-dependent mechanism. These data suggest that CCL18 may contribute to maintaining tolerance and/or suppressing deleterious inflammation by attracting memory Tregs into tissues, particularly in the lung, where it is highly and constitutively expressed.

Generation of second messengers in Plasmodium

Signalling in malaria parasites is a field of growing interest as its components may prove to be valuable drug targets, especially when one considers the burden of a disease that is responsible for up to 500 million infections annually. The scope of this review is to discuss external stimuli in the parasite life cycle and the upstream machinery responsible for translating them into intracellular responses, focussing particularly on the calcium signalling pathway.

Gametocytogenesis in malaria parasite: commitment, development and regulation

Malaria parasites have evolved a complicated life cycle alternating between two hosts. Gametocytes are produced in the vertebrate hosts and are obligatory for natural transmission of the parasites through mosquito vectors. The mechanism of sexual development in Plasmodium has been the focus of extensive studies. In the postgenomic era, the advent of genome-wide analytical tools and genetic manipulation technology has enabled rapid advancement of our knowledge in this area. Patterns of gene expression during sexual development, molecular distinction of the two sexes, and mechanisms underlying subsequent formation of gametes and their fertilization have been progressively elucidated. However, the triggers and mechanism of sexual development remain largely unknown. This article provides an update of our understanding of the molecular and cellular events associated with the decision for commitment to sexual development and regulation of gene expression during gametocytogenesis. Insights into the molecular mechanisms of gametocyte development are essential for designing proper control strategies for interruption of malaria transmission and ultimate elimination.

PfAlbas constitute a new eukaryotic DNA/RNA-binding protein family in malaria parasites

In Plasmodium falciparum, perinuclear subtelomeric chromatin conveys monoallelic expression of virulence genes. However, proteins that directly bind to chromosome ends are poorly described. Here we identify a novel DNA/RNA-binding protein family that bears homology to the archaeal protein Alba (Acetylation lowers binding affinity). We isolated three of the four PfAlba paralogs as part of a molecular complex that is associated with the P. falciparum-specific TARE6 (Telomere-Associated Repetitive Elements 6) subtelomeric region and showed in electromobility shift assays (EMSAs) that the PfAlbas bind to TARE6 repeats. In early blood stages, the PfAlba proteins were enriched at the nuclear periphery and partially co-localized with PfSir2, a TARE6-associated histone deacetylase linked to the process of antigenic variation. The nuclear location changed at the onset of parasite proliferation (trophozoite-schizont), where the PfAlba proteins were also detectable in the cytoplasm in a punctate pattern. Using single-stranded RNA (ssRNA) probes in EMSAs, we found that PfAlbas bind to ssRNA, albeit with different binding preferences. We demonstrate for the first time in eukaryotes that Alba-like proteins bind to both DNA and RNA and that their intracellular location is developmentally regulated. Discovery of the PfAlbas may provide a link between the previously described subtelomeric non-coding RNA and the regulation of antigenic variation.

A novel Toxoplasma gondii nuclear factor TgNF3 is a dynamic chromatin-associated component, modulator of nucleolar architecture and parasite virulence

In Toxoplasma gondii, cis-acting elements present in promoter sequences of genes that are stage-specifically regulated have been described. However, the nuclear factors that bind to these cis-acting elements and regulate promoter activities have not been identified. In the present study, we performed affinity purification, followed by proteomic analysis, to identify nuclear factors that bind to a stage-specific promoter in T. gondii. This led to the identification of several nuclear factors in T. gondii including a novel factor, designated herein as TgNF3. The N-terminal domain of TgNF3 shares similarities with the N-terminus of yeast nuclear FK506-binding protein (FKBP), known as a histone chaperone regulating gene silencing. Using anti-TgNF3 antibodies, HA-FLAG and YFP-tagged TgNF3, we show that TgNF3 is predominantly a parasite nucleolar, chromatin-associated protein that binds specifically to T. gondii gene promoters in vivo. Genome-wide analysis using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) identified promoter occupancies by TgNF3. In addition, TgNF3 has a direct role in transcriptional control of genes involved in parasite metabolism, transcription and translation. The ectopic expression of TgNF3 in the tachyzoites revealed dynamic changes in the size of the nucleolus, leading to a severe attenuation of virulence in vivo. We demonstrate that TgNF3 physically interacts with H3, H4 and H2A/H2B assembled into bona fide core and nucleosome-associated histones. Furthermore, TgNF3 interacts specifically to histones in the context of stage-specific gene silencing of a promoter that lacks active epigenetic acetylated histone marks. In contrast to virulent tachyzoites, which express the majority of TgNF3 in the nucleolus, the protein is exclusively located in the cytoplasm of the avirulent bradyzoites. We propose a model where TgNF3 acts essentially to coordinate nucleolus and nuclear functions by modulating nucleosome activities during the intracellular proliferation of the virulent tachyzoites of T. gondii.

Apicomplexa including Toxoplasma gondii are responsible for a variety of deadly infections. These intracellular parasites have complex life cycles within different hosts and their infectivity relies on their capacity to regulate gene expression in response to different environments. However, to date, little is known about nuclear factors that regulate their gene expression. Here, we have characterized parasite nuclear factors that bind to a stage-specific promoter. We identified several nuclear factors including a novel factor, designated herein as TgNF3. The N-terminal domain of TgNF3 shares similarities with the N-terminus of yeast nuclear FK506-binding protein (FKBP), known as a histone chaperone regulating gene silencing. We show that TgNF3 is predominantly a nucleolar, chromatin-associated protein that specifically binds to T. gondii nucleosome-associated histones and promoters. Genome-wide analysis identified promoter occupancies by TgNF3 and we demonstrated a direct role for this factor in transcriptional control of genes involved in parasite metabolism, transcription and translation. Ectopic expression of TgNF3 induces dynamic changes in the size of the nucleolus, and a severe attenuation of parasite virulence in vivo. In avirulent bradyzoites, TgNF3 is found exclusively in the cytoplasm, suggesting a potential role in regulating nucleolar and nuclear functions in the virulent tachyzoites of T. gondii.

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.

Core genome components and lineage specific expansions in malaria parasites plasmodium

Background

The increasing resistance of Plasmodium, the malaria parasites, to multiple commonly used drugs has underscored the urgent need to develop effective antimalarial drugs and vaccines. The new direction of genomics-driven target discovery has become possible with the completion of parasite genome sequencing, which can lead us to a better understanding of how the parasites develop the genetic variability that is associated with their response to environmental challenges and other adaptive phenotypes.

Results

We present the results of a comprehensive analysis of the genomes of six Plasmodium species, including two species that infect humans, one that infects monkeys, and three that infect rodents. The core genome shared by all six species is composed of 3,351 genes, which make up about 22%-65% of the genome repertoire. These components play important roles in fundamental functions as well as in parasite-specific activities. We further investigated the distribution and features of genes that have been expanded in specific Plasmodium lineage(s). Abundant duplicate genes are present in the six species, with 5%-9% of the whole genomes composed lineage specific radiations. The majority of these gene families are hypothetical proteins with unknown functions; a few may have predicted roles such as antigenic variation.

Conclusions

The core genome components in the malaria parasites have functions ranging from fundamental biological processes to roles in the complex networks that sustain the parasite-specific lifestyles appropriate to different hosts. They represent the minimum requirement to maintain a successful life cycle that spans vertebrate hosts and mosquito vectors. Lineage specific expansions (LSEs) have given rise to abundant gene families in Plasmodium. Although the functions of most families remain unknown, these LSEs could reveal components in parasite networks that, by their enhanced genetic variability, can contribute to pathogenesis, virulence, responses to environmental challenges, or interesting phenotypes.

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.