3' UTR signals necessary for expression of the Plasmodium gallinaceum ookinete protein, Pgs28, share similarities with those of yeast and plants

Nucleotide-resolution Mapping of RNA N6-Methyladenosine (m6A) modifications and comprehensive analysis of global polyadenylation events in mRNA 3' end processing in malaria pathogen Plasmodium falciparum

Plasmodium falciparum is an obligate human parasite of the phylum Apicomplexa and is the causative agent of the most lethal form of human malaria. Although N6-methyladenosine modification is thought to be one of the major post-transcriptional regulatory mechanisms for stage-specific gene expression in apicomplexan parasites, the precise base position of m6A in mRNAs or noncoding RNAs in these parasites remains unknown. Here, we report global nucleotide-resolution mapping of m6A residues in P. falciparum using DART-seq technology, which quantitatively displayed a stage-specific, dynamic distribution pattern with enrichment near mRNA 3' ends. In this process we identified 894, 788, and 1,762 m6A-modified genes in Ring, Trophozoite and Schizont stages respectively, with an average of 5-7 m6A sites per-transcript at the individual gene level. Notably, several genes involved in malaria pathophysiology, such as KAHRP, ETRAMPs, SERA and stress response genes, such as members of Heat Shock Protein (HSP) family are highly enriched in m6A and therefore could be regulated by this RNA modification. Since we observed preferential methylation at the 3' ends of P. falciparum transcripts and because malaria polyadenylation specificity factor PfCPSF30 harbors an m6A reader 'YTH' domain, we reasoned that m6A might play an important role in 3'-end processing of malaria mRNAs. To investigate this, we used two complementary high-throughput RNA 3'-end mapping approaches, which provided an initial framework to explore potential roles of m6A in the regulation of alternative polyadenylation (APA) during malaria development in human hosts.

Phenotypic Screens Identify Genetic Factors Associated with Gametocyte Development in the Human Malaria Parasite Plasmodium falciparum

Transmission of the deadly malaria parasite Plasmodium falciparum from humans to mosquitoes is achieved by specialized intraerythrocytic sexual forms called gametocytes. Though the crucial regulatory mechanisms leading to gametocyte commitment have recently come to light, networks of genes that control sexual development remain to be elucidated. Here, we report a pooled-mutant screen to identify genes associated with gametocyte development in P. falciparum. Our results categorized genes that modulate gametocyte progression as hypoproducers or hyperproducers of gametocytes, and the in-depth analysis of individual clones confirmed phenotypes in sexual commitment rates and putative functions in gametocyte development. We present a new set of genes that have not been implicated in gametocytogenesis before and demonstrate the potential of forward genetic screens in isolating genes impacting parasite sexual biology, an exciting step toward the discovery of new antimalarials for a globally significant pathogen. IMPORTANCE Blocking human-to-vector transmission is an essential step toward malaria elimination. Gametocytes are solely responsible for achieving this transmission and represent an opportunity for therapeutic intervention. While these falciform-shaped parasite stages were first discovered in the 1880s, our understanding of the genetic determinants responsible for their formation and molecular mechanisms that drive their development is limited. In this work, we developed a scalable screening methodology with piggyBac mutants to identify genes that influence the development of gametocytes in the most lethal human malaria parasite, P. falciparum. By doing so, we lay the foundation for large-scale functional genomic studies specifically designed to address remaining questions about sexual commitment, maturation, and mosquito infection in P. falciparum. Such functional genetic screens will serve to expedite the identification of essential pathways and processes for the development of novel transmission-blocking agents.

In silico identification of novel open reading frames in Plasmodium falciparum oocyte and salivary gland sporozoites using proteogenomics framework

Background

Plasmodium falciparum causes the deadliest form of malaria, which remains one of the most prevalent infectious diseases. Unfortunately, the only licensed vaccine showed limited protection and resistance to anti-malarial drug is increasing, which can be largely attributed to the biological complexity of the parasite’s life cycle. The progression from one developmental stage to another in P. falciparum involves drastic changes in gene expressions, where its infectivity to human hosts varies greatly depending on the stage. Approaches to identify candidate genes that are responsible for the development of infectivity to human hosts typically involve differential gene expression analysis between stages. However, the detection may be limited to annotated proteins and open reading frames (ORFs) predicted using restrictive criteria.

Methods

The above problem is particularly relevant for P. falciparum; whose genome annotation is relatively incomplete given its clinical significance. In this work, systems proteogenomics approach was used to address this challenge, as it allows computational detection of unannotated, novel Open Reading Frames (nORFs), which are neglected by conventional analyses. Two pairs of transcriptome/proteome were obtained from a previous study where one was collected in the mosquito-infectious oocyst sporozoite stage, and the other in the salivary gland sporozoite stage with human infectivity. They were then re-analysed using the proteogenomics framework to identify nORFs in each stage.

Results

Translational products of nORFs that map to antisense, intergenic, intronic, 3′ UTR and 5′ UTR regions, as well as alternative reading frames of canonical proteins were detected. Some of these nORFs also showed differential expression between the two life cycle stages studied. Their regulatory roles were explored through further bioinformatics analyses including the expression regulation on the parent reference genes, in silico structure prediction, and gene ontology term enrichment analysis.

Conclusion

The identification of nORFs in P. falciparum sporozoites highlights the biological complexity of the parasite. Although the analyses are solely computational, these results provide a starting point for further experimental validation of the existence and functional roles of these nORFs,

Characterization of mRNA polyadenylation in the apicomplexa

Messenger RNA polyadenylation is a universal aspect of gene expression in eukaryotes. In well-established model organisms, this process is mediated by a conserved complex of 15–20 subunits. To better understand this process in apicomplexans, a group of unicellular parasites that causes serious disease in humans and livestock, a computational and high throughput sequencing study of the polyadenylation complex and poly(A) sites in several species was conducted. BLAST-based searches for orthologs of the human polyadenylation complex yielded clear matches to only two—poly(A) polymerase and CPSF73—of the 19 proteins used as queries in this analysis. As the human subunits that recognize the AAUAAA polyadenylation signal (PAS) were not immediately obvious, a computational analysis of sequences adjacent to experimentally-determined apicomplexan poly(A) sites was conducted. The results of this study showed that there exists in apicomplexans an A-rich region that corresponds in position to the AAUAAA PAS. The set of experimentally-determined sites in one species, Sarcocystis neurona, was further analyzed to evaluate the extent and significance of alternative poly(A) site choice in this organism. The results showed that almost 80% of S. neurona genes possess more than one poly(A) site, and that more than 780 sites showed differential usage in the two developmental stages–extracellular merozoites and intracellular schizonts–studied. These sites affected more than 450 genes, and included a disproportionate number of genes that encode membrane transporters and ribosomal proteins. Taken together, these results reveal that apicomplexan species seem to possess a poly(A) signal analogous to AAUAAA even though genes that may encode obvious counterparts of the AAUAAA-recognizing proteins are absent in these organisms. They also indicate that, as is the case in other eukaryotes, alternative polyadenylation is a widespread phenomenon in S. neurona that has the potential to impact growth and development.

Translational repression of the cpw-wpc gene family in the malaria parasite Plasmodium

The technical challenges of working with the sexual stages of the malaria parasite Plasmodium have hindered the characterization of sexual stage antigens in the quest for a successful malaria transmission-blocking vaccine. One such predicted and largely uncharacterized group of sexual stage candidate antigens is the CPW-WPC family of proteins. CPW-WPC proteins are named for a characteristic domain that contains two conserved motifs, CPxxW and WPC. Conserved across Apicomplexa, this family is also present earlier in the Alveolata in the free-living, non-parasitophorous, photosynthetic chromerids, Chromera and Vitrella. In Plasmodium falciparum and Plasmodium berghei blood stage parasites, the transcripts of all nine cpw-wpc genes have been detected in gametocytes. RNA immunoprecipitation followed by reverse transcriptase-PCR reveals all P. berghei cpw-wpc transcripts to be bound by the translational repressors DOZI and CITH, and thus are likely under translational control prior to transmission from the rodent host to the mosquito vector in P. berghei. The GFP tagging of two endogenous P. berghei genes confirmed translational silencing in the gametocyte and translation in ookinetes. By establishing a luciferase transgene assay, we show that the 3' untranslated region of PF3D7_1331400 controls protein expression of this reporter in P. falciparum gametocytes. Our analyses suggest that cpw-wpc genes are translationally silenced in gametocytes across Plasmodium spp. and activated during ookinete formation and thus may have a role in transmission to the mosquito.

Functional analysis of the 5' untranslated region of the phosphoglutamase 2 transcript in Plasmodium falciparum

Plasmodium falciparum transcripts contain long untranslated regions (UTR), with some of the longest in any eukaryote that uses monocistronic transcription. Owing to the extreme AT nucleotide bias within the intergenic regions that encode these UTR, attempts to characterise how they are apportioned over genes and to describe their contribution to the absolute and temporal control of gene expression have been limited. Here we describe a study using a typical house-keeping gene that encodes phosphoglutamase 2 (PFD0660w), whose expression is subject to developmentally linked control during intraerythrocytic development. We show that deletion of a significant proportion (80%) of the predicted 5' UTR has no apparent effect on the developmentally linked expression of a luciferase reporter cassette. Further, serial deletions reveal that whilst the absolute level of transcription is unaffected when up to 50% of the predicted 5' UTR is removed, the subsequent efficiency of translation is affected. These data provide key insights into the interplay of transcriptional and post-transcriptional mechanisms in the control of gene expression in this important human pathogen.

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

Background

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

Results

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

Conclusions

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

Identification and characterization of a liver stage-specific promoter region of the malaria parasite Plasmodium

During the blood meal of a Plasmodium-infected mosquito, 10 to 100 parasites are inoculated into the skin and a proportion of these migrate via the bloodstream to the liver where they infect hepatocytes. The Plasmodium liver stage, despite its clinical silence, represents a highly promising target for antimalarial drug and vaccine approaches. Successfully invaded parasites undergo a massive proliferation in hepatocytes, producing thousands of merozoites that are transported into a blood vessel to infect red blood cells. To successfully develop from the liver stage into infective merozoites, a tight regulation of gene expression is needed. Although this is a very interesting aspect in the biology of Plasmodium, little is known about gene regulation in Plasmodium parasites in general and in the liver stage in particular. We have functionally analyzed a novel promoter region of the rodent parasite Plasmodium berghei that is exclusively active during the liver stage of the parasite. To prove stage-specific activity of the promoter, GFP and luciferase reporter assays have been successfully established, allowing both qualitative and accurate quantitative analysis. To further characterize the promoter region, the transcription start site was mapped by rapid amplification of cDNA ends (5'-RACE). Using promoter truncation experiments and site-directed mutagenesis within potential transcription factor binding sites, we suggest that the minimal promoter contains more than one binding site for the recently identified parasite-specific ApiAP2 transcription factors. The identification of a liver stage-specific promoter in P. berghei confirms that the parasite is able to tightly regulate gene expression during its life cycle. The identified promoter region might now be used to study the biology of the Plasmodium liver stage, which has thus far proven problematic on a molecular level. Stage-specific expression of dominant-negative mutant proteins and overexpression of proteins normally active in other life cycle stages will help to understand the function of the proteins investigated.

Control of gene expression in Plasmodium falciparum - ten years on

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

In silico analysis of EST and genomic sequences allowed the prediction of cis-regulatory elements for Entamoeba histolytica mRNA polyadenylation

In most eukaryotic cells, the poly(A) tail at the 3'-end of messenger RNA (mRNA) is essential for nuclear export, translatability, stability and transcription termination. Poly(A) tail formation involves multi-protein complexes that interact with specific sequences in 3'-untranslated region (3'-UTR) of precursor mRNA (pre-mRNA). Here we have performed a computational analysis of a large EST and genomic sequences collection from Entamoeba histolytica, the protozoan parasite responsible for human amoebiasis, to identify conserved elements that could be involved in pre-mRNA polyadenylation. Results evidenced the presence of an AU-rich domain corresponding to the consensus UA(A/U)UU polyadenylation signal or variants, the cleavage and polyadenylation site that is generally denoted by U residue and flanked by two U-rich tracts, and a novel A-rich element. This predicted array was validated through the analysis of genomic sequences and predicted mRNA folding of genes with known polyadenylation site. The molecular organization of pre-mRNA 3'-UTR cis-regulatory elements appears to be roughly conserved through evolutionary scale, whereas the polyadenylation signal seems to be species-specific in protozoan parasites and the novel A-rich element is unique for the primitive eukaryote E. histolytica. To our knowledge, this paper is the first work about the identification of potential pre-mRNA 3'-UTR cis-regulatory sequences through in silico analysis of large sets of cDNA and genomic sequences in a protozoan parasite.

3' gene regulatory elements required for expression of the Plasmodiumfalciparum developmental protein, Pfs25

Development of sexual stage parasites within the mosquito vector is a crucial step in the transmission of Plasmodium parasites. The expression of the P25 and P28 proteins on the surface of Plasmodium parasites in the mosquito midgut is required for development and hence disease transmission. 3' gene-flanking sequences are essential for expression of these critical proteins but the nucleotide elements required are poorly defined. Transient gene transfection experiments using constructs containing deletions of the 3' gene-flanking region of the Plasmodium falciparum P25 homologue, pfs25, reveal that elements necessary for protein expression are within 315 nucleotides (nt) of the stop codon. A T-rich region 137-231 nt from the stop codon is required for expression. The nonamer AATAAAATG, 360 nt downstream from the stop codon, enhances expression by 51 percent. Using 3' RACE analysis, multiple polyadenylation sites from endogenous and plasmid-derived pfs25 transcripts were identified. Dissimilarities between the identified elements and those of metazoans support the hypothesis that definition of P25/28 3' gene regulatory processes may eventually permit the development of agents which block malaria transmission but are non-toxic to humans.