Linker scanning mutagenesis of the Plasmodium gallinaceum sexual stage specific gene pgs28 reveals a novel downstream cis-control element

Transcriptomic analysis reveals novel downstream regulatory motifs and highly transcribed virulence factor genes of Entamoeba histolytica

Background

Promoter motifs in Entamoeba histolytica were earlier analysed using microarray data with lower dynamic range of gene expression. Additionally, previous transcriptomic studies did not provide information on the nature of highly transcribed genes, and downstream promoter motifs important for gene expression. To address these issues we generated RNA-Seq data and identified the high and low expressing genes, especially with respect to virulence potential. We analysed sequences both upstream and downstream of start site for important motifs.

Results

We used RNA-Seq data to classify genes according to expression levels, which ranged six orders of magnitude. Data were validated by reporter gene expression. Virulence-related genes (except AIG1) were amongst the highly expressed, while some kinases and BspA family genes were poorly expressed. We looked for conserved motifs in sequences upstream and downstream of the initiation codon. Following enrichment by AME we found seven motifs significantly enriched in high expression- and three in low expression-classes. Two of these motifs (M4 and M6) were located downstream of AUG, were exclusively enriched in high expression class, and were mostly found in ribosomal protein, and translation-related genes. Motif deletion resulted in drastic down regulation of reporter gene expression, showing functional relevance. Distribution of core promoter motifs (TATA, GAAC, and Inr) in all genes revealed that genes with downstream motifs were not preferentially associated with TATA-less promoters. We looked at gene expression changes in cells subjected to growth stress by serum starvation, and experimentally validated the data. Genes showing maximum up regulation belonged to the low or medium expression class, and included genes in signalling pathways, lipid metabolism, DNA repair, Myb transcription factors, BspA, and heat shock. Genes showing maximum down regulation belonged to the high or medium expression class. They included genes for signalling factors, actin, Ariel family, and ribosome biogenesis factors.

Conclusion

Our analysis has added important new information about the E. histolytica transcriptome. We report for the first time two downstream motifs required for gene expression, which could be used for over expression of E. histolytica genes. Most of the virulence-related genes in this parasite are highly expressed in culture.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5570-z) contains supplementary material, which is available to authorized users.

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.

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.

Role of cis-regulatory elements on the ring-specific hrp3 promoter in the human parasite Plasmodium falciparum

Identification of promoter elements responsible for regulation of gene expression has been hampered by the AT richness of P. falciparum intergenic regions. Nested deletions of histidine-rich protein 3 (hrp3) promoter suggested the presence of a multipartite ring-specific element. Linker scanning (LS) of this ring-specific promoter showed that the alteration of several promoter regions decreased the luciferase activity compared to the wild-type configuration, indicating that these regions played a role in gene expression. No homology was observed by comparison of putative regulatory elements of other genes identified by bioinformatic analysis with the hrp3 enhancer, implying a different mechanism of gene regulation by the hrp3 promoter. LS and deletion analysis of the 5' untranslated region (UTR) of the hrp3 suggested that this region contains elements which interact with promoter elements to regulate gene expression. Analysis of the intron in the UTR region suggested that this region does not play a role in stage specificity in the hrp3 promoter. Together, our results indicate the presence of multiple mechanisms of gene regulation in the parasite.

Gene regulation in the intraerythrocytic cycle of Plasmodium falciparum

MOTIVATION

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

RESULTS

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

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Eimeria tenella microneme protein EtMIC4: capture of the full-length transcribed sequence and comparison with other microneme proteins

EtMIC4 is a microneme protein of Eimeria tenella, an intracellular apicomplexan protozoan that can cause severe enteritis in chickens. The EtMIC4 gene has been partially characterised, and in this study, we used a combined strategy of rapid amplification of cDNA ends (5'RACE) and reverse transcription-polymerase chain reaction to identify the authentic 5' end of the transcribed sequence (accession number AJ306453.2). Comparison of the predicted EtMIC4 transcription start site with predicted start sites for EtMIC1, 2 and 3 genes identified comparable initiator regions that each conform to the consensus sequence for a transcriptional initiator element. The EtMIC4 gene is organised over 11 exons and analysis of the full-length predicted protein identified a new N-terminal region that comprises a hydrophobic signal peptide followed by four thrombospondin-type 1 modules that are similar to those previously described further downstream in the protein. Best-fit analysis shows that EtMIC4 shares high homology with the Eimeria maxima protein EmTFP250 and with TgMIC12, a predicted Toxoplasma gondii microneme protein. EtMIC4 and EmTFP250 share 70% amino acid identity and all predicted structural domains are conserved between the two. EtMIC4 and TgMIC12 share 48% identity and they have very similar domain organisation and conservation of intron/exon boundaries.

Computational analysis of constraints on noncoding regions, coding regions and gene expression in relation to Plasmodium phenotypic diversity

BACKGROUND

Malaria-causing Plasmodium species exhibit marked differences including host choice and preference for invading particular cell types. The genetic bases of phenotypic differences between parasites can be understood, in part, by investigating constraints on gene expression and genic sequences, both coding and regulatory.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated the evolutionary constraints on sequence and expression of parasitic genes by applying comparative genomics approaches to 6 Plasmodium genomes and 2 genome-wide expression studies. We found that the coding regions of Plasmodium transcription factor and sexual development genes are relatively less constrained, as are those of genes encoding CCCH zinc fingers and invasion proteins, which all play important roles in these parasites. Transcription factors and genes with stage-restricted expression have conserved upstream regions and so do several gene classes critical to the parasite's lifestyle, namely, ion transport, invasion, chromatin assembly and CCCH zinc fingers. Additionally, a cross-species comparison of expression patterns revealed that Plasmodium-specific genes exhibit significant expression divergence.

CONCLUSIONS/SIGNIFICANCE

Overall, constraints on Plasmodium's protein coding regions confirm observations from other eukaryotes in that transcription factors are under relatively lower constraint. Proteins relevant to the parasite's unique lifestyle also have lower constraint on their coding regions. Greater conservation between Plasmodium species in terms of promoter motifs suggests tight regulatory control of lifestyle genes. However, an interspecies divergence in expression patterns of these genes suggests that either expression is controlled via genomic or epigenomic features not encoded in the proximal promoter sequence, or alternatively, the combinatorial interactions between motifs confer species-specific expression patterns.

The transcription of bradyzoite genes in Toxoplasma gondii is controlled by autonomous promoter elements

Experimental evidence suggests that apicomplexan parasites possess bipartite promoters with basal and regulated cis-elements similar to other eukaryotes. Using a dual luciferase model adapted for recombinational cloning and use in Toxoplasma gondii, we show that genomic regions flanking 16 parasite genes, which encompass examples of constitutive and tachyzoite- and bradyzoite-specific genes, are able to reproduce the appropriate developmental stage expression in a transient luciferase assay. Mapping of cis-acting elements in several bradyzoite promoters led to the identification of short sequence spans that are involved in control of bradyzoite gene expression in multiple strains and under different bradyzoite induction conditions. Promoters that regulate the heat shock protein BAG1 and a novel bradyzoite-specific NTPase during bradyzoite development were fine mapped to a 6-8 bp resolution and these minimal cis-elements were capable of converting a constitutive promoter to one that is induced by bradyzoite conditions. Gel-shift experiments show that mapped cis-elements are bound by parasite protein factors with the appropriate functional sequence specificity. These studies are the first to identify the minimal sequence elements that are required and sufficient for bradyzoite gene expression and to show that bradyzoite promoters are maintained in a 'poised' chromatin state throughout the intermediate host life cycle in low passage strains. Together, these data demonstrate that conventional eukaryotic promoter mechanisms work with epigenetic processes to regulate developmental gene expression during tissue cyst formation.

In silico discovery of transcription regulatory elements in Plasmodium falciparum

BACKGROUND

With the sequence of the Plasmodium falciparum genome and several global mRNA and protein life cycle expression profiling projects now completed, elucidating the underlying networks of transcriptional control important for the progression of the parasite life cycle is highly pertinent to the development of new anti-malarials. To date, relatively little is known regarding the specific mechanisms the parasite employs to regulate gene expression at the mRNA level, with studies of the P. falciparum genome sequence having revealed few cis-regulatory elements and associated transcription factors. Although it is possible the parasite may evoke mechanisms of transcriptional control drastically different from those used by other eukaryotic organisms, the extreme AT-rich nature of P. falciparum intergenic regions (approximately 90% AT) presents significant challenges to in silico cis-regulatory element discovery.

RESULTS

We have developed an algorithm called Gene Enrichment Motif Searching (GEMS) that uses a hypergeometric-based scoring function and a position-weight matrix optimization routine to identify with high-confidence regulatory elements in the nucleotide-biased and repeat sequence-rich P. falciparum genome. When applied to promoter regions of genes contained within 21 co-expression gene clusters generated from P. falciparum life cycle microarray data using the semi-supervised clustering algorithm Ontology-based Pattern Identification, GEMS identified 34 putative cis-regulatory elements associated with a variety of parasite processes including sexual development, cell invasion, antigenic variation and protein biosynthesis. Among these candidates were novel motifs, as well as many of the elements for which biological experimental evidence already exists in the Plasmodium literature. To provide evidence for the biological relevance of a cell invasion-related element predicted by GEMS, reporter gene and electrophoretic mobility shift assays were conducted.

CONCLUSION

This GEMS analysis demonstrates that in silico regulatory element discovery can be successfully applied to challenging repeat-sequence-rich, base-biased genomes such as that of P. falciparum. The fact that regulatory elements were predicted from a diverse range of functional gene clusters supports the hypothesis that cis-regulatory elements play a role in the transcriptional control of many P. falciparum biological processes. The putative regulatory elements described represent promising candidates for future biological investigation into the underlying transcriptional control mechanisms of gene regulation in malaria parasites.

Plasmodium in the postgenomic era: new insights into the molecular cell biology of malaria parasites

In this review, we bring together some of the approaches toward understanding the cellular and molecular biology of Plasmodium species and their interaction with their host red blood cells. Considerable impetus has come from the development of new methods of molecular genetics and bioinformatics, and it is important to evaluate the wealth of these novel data in the context of basic cell biology. We describe how these approaches are gaining valuable insights into the parasite-host cell interaction, including (1) the multistep process of red blood cell invasion by the merozoite; (2) the mechanisms by which the intracellular parasite feeds on the red blood cell and exports parasite proteins to modify its cytoadherent properties; (3) the modulation of the cell cycle by sensing the environmental tryptophan-related molecules; (4) the mechanism used to survive in a low Ca(2+) concentration inside red blood cells; (5) the activation of signal transduction machinery and the regulation of intracellular calcium; (6) transfection technology; and (7) transcriptional regulation and genome-wide mRNA studies in Plasmodium falciparum.

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

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

Promoter regions of Plasmodium vivax are poorly or not recognized by Plasmodium falciparum

Background

Heterologous promoter analysis in Plasmodium has revealed the existence of conserved cis regulatory elements as promoters from different species can drive expression of reporter genes in heterologous transfection assays. Here, the functional characterization of different Plasmodium vivax promoters in Plasmodium falciparum using luciferase as the reporter gene is presented.

Methods

Luciferase reporter plasmids harboring the upstream regions of the msp1, dhfr, and vir3 genes as well as the full-length intergenic regions of the vir23/24 and ef-1α genes of P. vivax were constructed and transiently transfected in P. falciparum.

Results

Only the constructs with the full-length intergenic regions of the vir23/24 and ef-1α genes were recognized by the P. falciparum transcription machinery albeit to values approximately two orders of magnitude lower than those reported by luc plasmids harbouring promoter regions from P. falciparum and Plasmodium berghei. A bioinformatics approach allowed the identification of a motif (GCATAT) in the ef-1α intergenic region that is conserved in five Plasmodium species but is degenerate (GCANAN) in P. vivax. Mutations of this motif in the P. berghei ef-1α promoter region decreased reporter expression indicating it is active in gene expression in Plasmodium.

Conclusion

Together, this data indicates that promoter regions of P. vivax are poorly or not recognized by the P. falciparum transcription machinery suggesting the existence of P. vivax-specific transcription regulatory elements.

Regulatory motifs uncovered among gene expression clusters in Plasmodium falciparum

Control of gene expression is poorly understood in the Plasmodium system, where relatively few homologues to known eukaryotic transcription factors have been uncovered. Recent evidence suggests that the parasite may utilize a combinatorial mode of gene regulation, with multiple cis-acting sequences contributing to overall activity at individual promoters [1]. To further probe this mechanism of control, we first searched for over-represented sequence motifs among gene clusters sharing similar expression profiles in Plasmodium falciparum. More specifically, we applied bioinformatic tools to a previously characterized micro-array data set from drug-treated asexual stage cultures (Gunasekera et al., submitted). Cluster analysis of 600 drug responsive genes identified only a single 5' motif, GAGAGAA. Two additional 5' motifs, ACTATAAAGA and TGCAC, were also shared among loci displaying patterns of coordinate expression across varying asexual growth stages. Secondly and most importantly, the functional relevance of each motif was tested in two independent assays-transient transfection and gel-retardation experiments. The GAGAGAA and TGCAC motifs were both active in the former. The GAGAGAA and ACTATAAAGA elements formed specific RNA-protein, but not DNA-protein complexes in gel shift assays, suggesting a key level of control at the RNA level. This is the first report of functionally characterized motifs in P. falciparum that were uncovered following clustering analysis of its asexual stage transcriptome. Together, both the bioinformatic and functional data reported here imply that multiple forms of gene regulation, including post-transcriptional control, may be important in the malarial system.

Identification of regulatory elements in the Plasmodium falciparum genome

There is little information regarding regulatory sequences in the newly sequenced genome of the malaria parasite, Plasmodium falciparum. Thus, for the first time, a bioinformatic strategy was utilized to identify regulatory elements in this genome using the P. falciparum heat shock protein (hsp) gene family as a model system. Our analysis indicates that the P. falciparum hsp genes do not contain standard eukaryotic regulatory elements. However, a novel G-rich regulatory element named the G-box was identified upstream of several P. falciparum hsp genes and the P. yoelii yoelii, P. berghei, and P. vivax hsp86 genes. Remarkably, the Plasmodium sp. G-boxes were required for maximal reporter gene expression in transient transfection assays. The G-box is not homologous to known eukaryotic elements, and is the best-defined functional element elucidated from Plasmodium sp. Our analysis also revealed several other elements necessary for reporter gene expression including an upstream sequence element, the region surrounding the transcription start site, and the 5' and 3' untranslated regions. These data demonstrate that unique regulatory elements are conserved in the genomes of Plasmodium sp., and demonstrate the feasibility of bioinformatic approaches for their identification.