Polysome profiling reveals translational control of gene expression in the human malaria parasite Plasmodium falciparum

Polysome Profiling in Leishmania, Human Cells and Mouse Testis

Proper protein expression at the right time and in the right amounts is the basis of normal cell function and survival in a fast-changing environment. For a long time, the gene expression studies were dominated by research on the transcriptional level. However, the steady-state levels of mRNAs do not correlate well with protein production, and the translatability of mRNAs varies greatly depending on the conditions. In some organisms, like the parasite Leishmania, the protein expression is regulated mostly at the translational level. Recent studies demonstrated that protein translation dysregulation is associated with cancer, metabolic, neurodegenerative and other human diseases. Polysome profiling is a powerful method to study protein translation regulation. It allows to measure the translational status of individual mRNAs or examine translation on a genome-wide scale. The basis of this technique is the separation of polysomes, ribosomes, their subunits and free mRNAs during centrifugation of a cytoplasmic lysate through a sucrose gradient. Here, we present a universal polysome profiling protocol used on three different models - parasite Leishmania major, cultured human cells and animal tissues. Leishmania cells freely grow in suspension and cultured human cells grow in adherent monolayer, while mouse testis represents an animal tissue sample. Thus, the technique is adapted to all of these sources. The protocol for the analysis of polysomal fractions includes detection of individual mRNA levels by RT-qPCR, proteins by Western blot and analysis of ribosomal RNAs by electrophoresis. The method can be further extended by examination of mRNAs association with the ribosome on a transcriptome level by deep RNA-seq and analysis of ribosome-associated proteins by mass spectroscopy of the fractions. The method can be easily adjusted to other biological models.

Identification of human Phosphatidyl Inositol 5-Phosphate 4-Kinase as an RNA binding protein that is imported into Plasmodium falciparum

Plasmodium falciparum is a causative agent for malaria and has a complex life cycle in human and mosquito hosts. Translation repression of specific set of mRNA has been reported in gametocyte stages of this parasite. A conserved element present in the 3'UTR of some of these transcripts was identified. Biochemical studies have identified components of the RNA storage and/or translation inhibitor complex but it is not yet clear how the complex is specifically recruited on the RNA targeted for translation regulation. We used the 3'UTR region of translationally regulated transcripts to identify Phosphatidyl-inositol 5-phosphate 4-kinase (PIP4K2A) as the protein that associates with these RNAs. We further show that recombinant PIP4K2A has the RNA binding activity and can associate specifically with Plasmodium 3'UTR RNAs. Immunostainings show that hPIP4K2A is imported into the Plasmodium parasite from RBC. These results identify a novel RNA binding role for PIP4K2A that may play a role in Plasmodium propagation.

Integrative Approaches to Understand the Mastery in Manipulation of Host Cytokine Networks by Protozoan Parasites with Emphasis on Plasmodium and Leishmania Species

Diseases by protozoan pathogens pose a significant public health concern, particularly in tropical and subtropical countries, where these are responsible for significant morbidity and mortality. Protozoan pathogens tend to establish chronic infections underscoring their competence at subversion of host immune processes, an important component of disease pathogenesis and of their virulence. Modulation of cytokine and chemokine levels, their crosstalks and downstream signaling pathways, and thereby influencing recruitment and activation of immune cells is crucial to immune evasion and subversion. Many protozoans are now known to secrete effector molecules that actively modulate host immune transcriptome and bring about alterations in host epigenome to alter cytokine levels and signaling. The complexity of multi-dimensional events during interaction of hosts and protozoan parasites ranges from microscopic molecular levels to macroscopic ecological and epidemiological levels that includes disrupting metabolic pathways, cell cycle (Toxoplasma and Theileria sp.), respiratory burst, and antigen presentation (Leishmania spp.) to manipulation of signaling hubs. This requires an integrative systems biology approach to combine the knowledge from all these levels to identify the complex mechanisms of protozoan evolution via immune escape during host-parasite coevolution. Considering the diversity of protozoan parasites, in this review, we have focused on Leishmania and Plasmodium infections. Along with the biological understanding, we further elucidate the current efforts in generating, integrating, and modeling of multi-dimensional data to explain the modulation of cytokine networks by these two protozoan parasites to achieve their persistence in host via immune escape during host-parasite coevolution.

Next generation sequencing reveals widespread trypanosome diversity and polyparasitism in marsupials from Western Australia

In Western Australia a number of indigenous Trypanosoma spp. infect susceptible native marsupials, such as the woylie (Bettongia penicillata), brushtail possum (Trichosurus vulpecula), and chuditch (Dasyurus geoffroii). Two genotypes of Trypanosoma copemani (identified as G1 and G2) have been found in the woylie, and G2 has been implicated in the decline of this host species, making its presence of particular interest. Here we used targeted amplicon next generation sequencing (NGS) of the Trypanosoma 18S rDNA loci on 70 Trypanosoma-positive marsupial blood samples, to identify T. copemani genotypes and multiple Trypanosoma infections (polyparasitism) in woylies and cohabiting species in Western Australia. Polyparasitism with Trypanosoma spp. was found in 50% of the wildlife sampled, and within species diversity was high, with 85 zero-radius operational taxonomic units (ZOTUs) identified in nine putative parasite species. Trypanosoma copemani was assigned 17 ZOTUs and was identified in 80% of samples. The most abundant ZOTU isolated (63%) differed slightly from the published genotype of G1, and G2 was the second most abundant ZOTU (14%). Trypanosome diversity was significantly greater in woylies than in brushtail possums, and parasite community composition also differed significantly between these host species. One novel Trypanosoma spp. genotype (Trypanosoma sp. ANU2) was found in 20% of samples. A species of Crithidia was detected in a woylie, and two avian trypanosomes (Trypanosoma avium and Trypanosoma sp. AAT) were identified in woylies for the first time.

•

Extensive trypanosome diversity and polyparasitism in south Western Australia.

•

A significant difference between trypanosomes infecting woylies and brushtail possums.

•

Evidence that the current species-level taxonomy of Australian trypanosomes may need revision.

Plasmodium knowlesi: a superb in vivo nonhuman primate model of antigenic variation in malaria

Antigenic variation in malaria was discovered in Plasmodium knowlesi studies involving longitudinal infections of rhesus macaques (M. mulatta). The variant proteins, known as the P. knowlesi Schizont Infected Cell Agglutination (SICA) antigens and the P. falciparum Erythrocyte Membrane Protein 1 (PfEMP1) antigens, expressed by the SICAvar and var multigene families, respectively, have been studied for over 30 years. Expression of the SICA antigens in P. knowlesi requires a splenic component, and specific antibodies are necessary for variant antigen switch events in vivo. Outstanding questions revolve around the role of the spleen and the mechanisms by which the expression of these variant antigen families are regulated. Importantly, the longitudinal dynamics and molecular mechanisms that govern variant antigen expression can be studied with P. knowlesi infection of its mammalian and vector hosts. Synchronous infections can be initiated with established clones and studied at multi-omic levels, with the benefit of computational tools from systems biology that permit the integration of datasets and the design of explanatory, predictive mathematical models. Here we provide an historical account of this topic, while highlighting the potential for maximizing the use of P. knowlesi - macaque model systems and summarizing exciting new progress in this area of research.

EuPathDB: The Eukaryotic Pathogen Genomics Database Resource

Fighting infections and developing novel drugs and vaccines requires advanced knowledge of pathogen's biology. Readily accessible genomic, functional genomic, and population data aids biological and translational discovery. The Eukaryotic Pathogen Database Resources ( http://eupathdb.org ) are data mining resources that support hypothesis driven research by facilitating the discovery of meaningful biological relationships from large volumes of data. The resource encompasses 13 sites that support over 170 species including pathogenic protists, oomycetes, and fungi as well as evolutionarily related nonpathogenic species. EuPathDB integrates preanalyzed data with advanced search capabilities, data visualization, analysis tools and a comprehensive record system in a graphical interface that does not require prior computational skills. This chapter describes guiding concepts common across EuPathDB sites and illustrates the powerful data mining capabilities of some of the available tools and features.

Adaptation of Plasmodium falciparum to its transmission environment

Success in eliminating malaria will depend on whether parasite evolution outpaces control efforts. Here, we show that Plasmodium falciparum parasites (the deadliest of the species causing human malaria) found in low-transmission-intensity areas have evolved to invest more in transmission to new hosts (reproduction) and less in within-host replication (growth) than parasites found in high-transmission areas. At the cellular level, this adaptation manifests as increased production of reproductive forms (gametocytes) early in the infection at the expense of processes associated with multiplication inside red blood cells, especially membrane transport and protein trafficking. At the molecular level, this manifests as changes in the expression levels of genes encoding epigenetic and translational machinery. Specifically, expression levels of the gene encoding AP2-G-the transcription factor that initiates reproduction-increase as transmission intensity decreases. This is accompanied by downregulation and upregulation of genes encoding HDAC1 and HDA1-two histone deacetylases that epigenetically regulate the parasite's replicative and reproductive life-stage programmes, respectively. Parasites in reproductive mode show increased reliance on the prokaryotic translation machinery found inside the plastid-derived organelles. Thus, our dissection of the parasite's adaptive regulatory architecture has identified new potential molecular targets for malaria control.

Comparative ribosome profiling uncovers a dominant role for translational control in Toxoplasma gondii

Background

The lytic cycle of the protozoan parasite Toxoplasma gondii, which involves a brief sojourn in the extracellular space, is characterized by defined transcriptional profiles. For an obligate intracellular parasite that is shielded from the cytosolic host immune factors by a parasitophorous vacuole, the brief entry into the extracellular space is likely to exert enormous stress. Due to its role in cellular stress response, we hypothesize that translational control plays an important role in regulating gene expression in Toxoplasma during the lytic cycle. Unlike transcriptional profiles, insights into genome-wide translational profiles of Toxoplasma gondii are lacking.

Methods

We have performed genome-wide ribosome profiling, coupled with high throughput RNA sequencing, in intracellular and extracellular Toxoplasma gondii parasites to investigate translational control during the lytic cycle.

Results

Although differences in transcript abundance were mostly mirrored at the translational level, we observed significant differences in the abundance of ribosome footprints between the two parasite stages. Furthermore, our data suggest that mRNA translation in the parasite is potentially regulated by mRNA secondary structure and upstream open reading frames.

Conclusion

We show that most of the Toxoplasma genes that are dysregulated during the lytic cycle are translationally regulated.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4362-6) contains supplementary material, which is available to authorized users.

New insights into the factors affecting synonymous codon usage in human infecting Plasmodium species

Codon usage bias is due to the non-random usage of synonymous codons for coding amino acids. The synonymous sites are under weak selection, and codon usage bias is maintained by the equilibrium in mutational bias, genetic drift and selection pressure. The differential codon usage choices are also relevant to human infecting Plasmodium species. Recently, P. knowlesi switches its natural host, long-tailed macaques, and starts infecting humans. This review focuses on the comparative analysis of codon usage choices among human infecting P. falciparum and P. vivax along with P. knowlesi species taking their coding sequence data. The variation in GC content, amino acid frequencies, effective number of codons and other factors plays a crucial role in determining synonymous codon choices. Within species codon choices are more similar for P. vivax and P. knowlesi in comparison with P. falciparum species. This study suggests that synonymous codon choice modulates the gene expression level, mRNA stability, ribosome speed, protein folding, translation efficiency and its accuracy in Plasmodium species, and provides a valuable information regarding the codon usage pattern to facilitate gene cloning as well as expression and transfection studies for malaria causing species.

Plasmodium falciparum Calcium-Dependent Protein Kinase 2 Is Critical for Male Gametocyte Exflagellation but Not Essential for Asexual Proliferation

Drug development efforts have focused mostly on the asexual blood stages of the malaria parasite Plasmodium falciparum. Except for primaquine, which has its own limitations, there are no available drugs that target the transmission of the parasite to mosquitoes. Therefore, there is a need to validate new parasite proteins that can be targeted for blocking transmission. P. falciparum calcium-dependent protein kinases (PfCDPKs) play critical roles at various stages of the parasite life cycle and, importantly, are absent in the human host. These features mark them as attractive drug targets. In this study, using CRISPR/Cas9 we successfully knocked out PfCDPK2 from blood-stage parasites, which was previously thought to be an indispensable protein. The growth rate of the PfCDPK2 knockout (KO) parasites was similar to that of wild-type parasites, confirming that PfCDPK2 function is not essential for the asexual proliferation of the parasite in vitro. The mature male and female gametocytes of PfCDPK2 KO parasites become round after induction. However, they fail to infect female Anopheles stephensi mosquitoes due to a defect(s) in male gametocyte exflagellation and possibly in female gametes.

Despite reductions in the number of deaths it causes, malaria continues to be a leading infectious disease of the developing world. For effective control and elimination of malaria, multiple stages of the parasite need to be targeted. One such stage includes the transmission of the parasite to mosquitoes. Here, we demonstrate the successful knockout of PfCDPK2, which was previously thought to be indispensable for parasite growth in red blood cells. The PfCDPK2 KO parasites are incapable of establishing an infection in mosquitoes. Therefore, our study suggests that targeting PfCDPK2 may be a good strategy to control malaria transmission in countries with high transmission. Moreover, molecular understanding of the signaling pathway of PfCDPK2 may provide additional targets for malaria control.

Translational Control in the Latency of Apicomplexan Parasites

Apicomplexan parasites Toxoplasma gondii and Plasmodium spp. use latent stages to persist in the host, facilitate transmission, and thwart treatment of infected patients. Therefore, it is important to understand the processes driving parasite differentiation to and from quiescent stages. Here, we discuss how a family of protein kinases that phosphorylate the eukaryotic initiation factor-2 (eIF2) function in translational control and drive differentiation. This translational control culminates in reprogramming of the transcriptome to facilitate parasite transition towards latency. We also discuss how eIF2 phosphorylation contributes to the maintenance of latency and provides a crucial role in the timing of reactivation of latent parasites towards proliferative stages.

ALBA4 modulates its stage-specific interactions and specific mRNA fates during Plasmodium yoelii growth and transmission

Transmission of the malaria parasite occurs in an unpredictable moment, when a mosquito takes a blood meal. Plasmodium has therefore evolved strategies to prepare for transmission, including translationally repressing and protecting mRNAs needed to establish the infection. However, mechanisms underlying these critical controls are not well understood, including whether Plasmodium changes its translationally repressive complexes and mRNA targets in different stages. Efforts to understand this have been stymied by severe technical limitations due to substantial mosquito contamination of samples. Here using P. yoelii, for the first time we provide a proteomic comparison of a protein complex across asexual blood, sexual and sporozoite stages, along with a transcriptomic comparison of the mRNAs that are affected in these stages. We find that the Apicomplexan-specific ALBA4 RNA-binding protein acts to regulate development of the parasite's transmission stages, and that ALBA4 associates with both stage-specific and stage-independent partners to produce opposing mRNA fates. These efforts expand our understanding and ability to interrogate both sexual and sporozoite transmission stages and the molecular preparations they evolved to perpetuate their infectious cycle.

Nascent RNA sequencing reveals mechanisms of gene regulation in the human malaria parasite Plasmodium falciparum

Gene expression in Plasmodium falciparum is tightly regulated to ensure successful propagation of the parasite throughout its complex life cycle. The earliest transcriptomics studies in P. falciparum suggested a cascade of transcriptional activity over the course of the 48-hour intraerythrocytic developmental cycle (IDC); however, the just-in-time transcriptional model has recently been challenged by findings that show the importance of post-transcriptional regulation. To further explore the role of transcriptional regulation, we performed the first genome-wide nascent RNA profiling in P. falciparum. Our findings indicate that the majority of genes are transcribed simultaneously during the trophozoite stage of the IDC and that only a small subset of genes is subject to differential transcriptional timing. RNA polymerase II is engaged with promoter regions prior to this transcriptional burst, suggesting that Pol II pausing plays a dominant role in gene regulation. In addition, we found that the overall transcriptional program during gametocyte differentiation is surprisingly similar to the IDC, with the exception of relatively small subsets of genes. Results from this study suggest that further characterization of the molecular players that regulate stage-specific gene expression and Pol II pausing will contribute to our continuous search for novel antimalarial drug targets.

Increased genetic diversity and prevalence of co-infection with Trypanosoma spp. in koalas (Phascolarctos cinereus) and their ticks identified using next-generation sequencing (NGS)

Infections with Trypanosoma spp. have been associated with poor health and decreased survival of koalas (Phascolarctos cinereus), particularly in the presence of concurrent pathogens such as Chlamydia and koala retrovirus. The present study describes the application of a next-generation sequencing (NGS)-based assay to characterise the prevalence and genetic diversity of trypanosome communities in koalas and two native species of ticks (Ixodes holocyclus and I. tasmani) removed from koala hosts. Among 168 koalas tested, 32.2% (95% CI: 25.2–39.8%) were positive for at least one Trypanosoma sp. Previously described Trypanosoma spp. from koalas were identified, including T. irwini (32.1%, 95% CI: 25.2–39.8%), T. gilletti (25%, 95% CI: 18.7–32.3%), T. copemani (27.4%, 95% CI: 20.8–34.8%) and T. vegrandis (10.1%, 95% CI: 6.0–15.7%). Trypanosoma noyesi was detected for the first time in koalas, although at a low prevalence (0.6% 95% CI: 0–3.3%), and a novel species (Trypanosoma sp. AB-2017) was identified at a prevalence of 4.8% (95% CI: 2.1–9.2%). Mixed infections with up to five species were present in 27.4% (95% CI: 21–35%) of the koalas, which was significantly higher than the prevalence of single infections 4.8% (95% CI: 2–9%). Overall, a considerably higher proportion (79.7%) of the Trypanosoma sequences isolated from koala blood samples were identified as T. irwini, suggesting this is the dominant species. Co-infections involving T. gilletti, T. irwini, T. copemani, T. vegrandis and Trypanosoma sp. AB-2017 were also detected in ticks, with T. gilletti and T. copemani being the dominant species within the invertebrate hosts. Direct Sanger sequencing of Trypanosoma 18S rRNA gene amplicons was also performed and results revealed that this method was only able to identify the genotypes with greater amount of reads (according to NGS) within koala samples, which highlights the advantages of NGS in detecting mixed infections. The present study provides new insights on the natural genetic diversity of Trypanosoma communities infecting koalas and constitutes a benchmark for future clinical and epidemiological studies required to quantify the contribution of trypanosome infections on koala survival rates.

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.

Capturing in vivo RNA transcriptional dynamics from the malaria parasite Plasmodium falciparum

To capture the transcriptional dynamics within proliferating cells, methods to differentiate nascent transcription from preexisting mRNAs are desired. One approach is to label newly synthesized mRNA transcripts in vivo through the incorporation of modified pyrimidines. However, the human malaria parasite, Plasmodium falciparum, is incapable of pyrimidine salvage for mRNA biogenesis. To capture cellular mRNA dynamics during Plasmodium development, we engineered parasites that can salvage pyrimidines through the expression of a single bifunctional yeast fusion gene, cytosine deaminase/uracil phosphoribosyltransferase (FCU). We show that expression of FCU allows for the direct incorporation of thiol-modified pyrimidines into nascent mRNAs. Using developmental stage-specific promoters to express FCU-GFP enables the biosynthetic capture and in-depth analysis of mRNA dynamics from subpopulations of cells undergoing differentiation. We demonstrate the utility of this method by examining the transcriptional dynamics of the sexual gametocyte stage transition, a process that is essential to malaria transmission between hosts. Using the pfs16 gametocyte-specific promoter to express FCU-GFP in 3D7 parasites, we found that sexual stage commitment is governed by transcriptional reprogramming and stabilization of a subset of essential gametocyte transcripts. We also measured mRNA dynamics in F12 gametocyte-deficient parasites and demonstrate that the transcriptional program required for sexual commitment and maturation is initiated but likely aborted due to the absence of the PfAP2-G transcriptional regulator and a lack of gametocyte-specific mRNA stabilization. Biosynthetic labeling of Plasmodium mRNAs is incredibly versatile, can be used to measure transcriptional dynamics at any stage of parasite development, and will allow for future applications to comprehensively measure RNA-protein interactions in the malaria parasite.

Frequent GU wobble pairings reduce translation efficiency in Plasmodium falciparum

Plasmodium falciparum genome has 81% A+T content. This nucleotide bias leads to extreme codon usage bias and culminates in frequent insertion of asparagine homorepeats in the proteome. Using recodonized GFP sequences, we show that codons decoded via G:U wobble pairing are suboptimal codons that are negatively associated to protein translation efficiency. Despite this, one third of all codons in the genome are GU wobble codons, suggesting that codon usage in P. falciparum has not been driven to maximize translation efficiency, but may have evolved as translational regulatory mechanism. Particularly, asparagine homorepeats are generally encoded by locally clustered GU wobble AAT codons, we demonstrated that this GU wobble-rich codon context is the determining factor that causes reduction of protein level. Moreover, insertion of clustered AAT codons also causes destabilization of the transcripts. Interestingly, more frequent asparagine homorepeats insertion is seen in single-exon genes, suggesting transcripts of these genes may have been programmed for rapid mRNA decay to compensate for the inefficiency of mRNA surveillance regulation on intronless genes. To our knowledge, this is the first study that addresses P. falciparum codon usage in vitro and provides new insights on translational regulation and genome evolution of this parasite.

Untranslated regions of mRNA and their role in regulation of gene expression in protozoan parasites

Protozoan parasites are one of the oldest living entities in this world that throughout their existence have shown excellent resilience to the odds of survival and have adapted beautifully to ever changing rigors of the environment. In view of the dynamic environment encountered by them throughout their life cycle, and in establishing pathogenesis, it is unsurprising that modulation of gene expression plays a fundamental role in their survival. In higher eukaryotes, untranslated regions (UTRs) of transcripts are one of the crucial regulators of gene expression (influencing mRNA stability and translation efficiency). Parasitic protozoan genome studies have led to the characterization (in silico, in vitro and in vivo) of a large number of their genes. Comparison of higher eukaryotic UTRs with parasitic protozoan UTRs reveals the existence of several similar and dissimilar facets of the UTRs. This review focuses on the elements of UTRs of medically important protozoan parasites and their regulatory role in gene expression. Such information may be useful to researchers in designing gene targeting strategies linked with perturbation of host-parasite relationships leading to control of specific parasites.

Clustering Algorithms: Their Application to Gene Expression Data

Gene expression data hide vital information required to understand the biological process that takes place in a particular organism in relation to its environment. Deciphering the hidden patterns in gene expression data proffers a prodigious preference to strengthen the understanding of functional genomics. The complexity of biological networks and the volume of genes present increase the challenges of comprehending and interpretation of the resulting mass of data, which consists of millions of measurements; these data also inhibit vagueness, imprecision, and noise. Therefore, the use of clustering techniques is a first step toward addressing these challenges, which is essential in the data mining process to reveal natural structures and identify interesting patterns in the underlying data. The clustering of gene expression data has been proven to be useful in making known the natural structure inherent in gene expression data, understanding gene functions, cellular processes, and subtypes of cells, mining useful information from noisy data, and understanding gene regulation. The other benefit of clustering gene expression data is the identification of homology, which is very important in vaccine design. This review examines the various clustering algorithms applicable to the gene expression data in order to discover and provide useful knowledge of the appropriate clustering technique that will guarantee stability and high degree of accuracy in its analysis procedure.

EuPathDB: the eukaryotic pathogen genomics database resource

The Eukaryotic Pathogen Genomics Database Resource (EuPathDB, http://eupathdb.org) is a collection of databases covering 170+ eukaryotic pathogens (protists & fungi), along with relevant free-living and non-pathogenic species, and select pathogen hosts. To facilitate the discovery of meaningful biological relationships, the databases couple preconfigured searches with visualization and analysis tools for comprehensive data mining via intuitive graphical interfaces and APIs. All data are analyzed with the same workflows, including creation of gene orthology profiles, so data are easily compared across data sets, data types and organisms. EuPathDB is updated with numerous new analysis tools, features, data sets and data types. New tools include GO, metabolic pathway and word enrichment analyses plus an online workspace for analysis of personal, non-public, large-scale data. Expanded data content is mostly genomic and functional genomic data while new data types include protein microarray, metabolic pathways, compounds, quantitative proteomics, copy number variation, and polysomal transcriptomics. New features include consistent categorization of searches, data sets and genome browser tracks; redesigned gene pages; effective integration of alternative transcripts; and a EuPathDB Galaxy instance for private analyses of a user's data. Forthcoming upgrades include user workspaces for private integration of data with existing EuPathDB data and improved integration and presentation of host–pathogen interactions.

Genome-wide regulatory dynamics of G-quadruplexes in human malaria parasite Plasmodium falciparum

The AT-rich genome of P. falciparum has uniquely localized G-rich stretches that have propensity to form G-quadruplexes. However, their global occurrence and potential biological roles in the parasite are poorly explored. Our genome-wide analysis revealed unique enrichment of quadruplexes in P. falciparum genome which was remarkably different from other Plasmodium species. A distinct predominance of quadruplexes was observed in nuclear and organellar genes that participate in antigenic variation, pathogenesis, DNA/RNA regulation, metabolic and protein quality control processes. Data also suggested association of quadruplexes with SNPs and DNA methylation. Furthermore, analysis of steady state mRNA (RNA-seq) and polysome-associated mRNA (Ribosome profiling) data revealed stage-specific differences in translational efficiency of quadruplex harboring genes. Taken together, our findings hint towards existence of regulatory dynamics associated with quadruplexes that may modulate translational efficiency of quadruplex harboring genes to provide survival advantage to the parasite against host immune response and antimalarial drug pressure.

Evolutionary Dynamics of Abundant Stop Codon Readthrough

Translational stop codon readthrough emerged as a major regulatory mechanism affecting hundreds of genes in animal genomes, based on recent comparative genomics and ribosomal profiling evidence, but its evolutionary properties remain unknown. Here, we leverage comparative genomic evidence across 21 Anopheles mosquitoes to systematically annotate readthrough genes in the malaria vector Anopheles gambiae, and to provide the first study of abundant readthrough evolution, by comparison with 20 Drosophila species. Using improved comparative genomics methods for detecting readthrough, we identify evolutionary signatures of conserved, functional readthrough of 353 stop codons in the malaria vector, Anopheles gambiae, and of 51 additional Drosophila melanogaster stop codons, including several cases of double and triple readthrough and of readthrough of two adjacent stop codons. We find that most differences between the readthrough repertoires of the two species arose from readthrough gain or loss in existing genes, rather than birth of new genes or gene death; that readthrough-associated RNA structures are sometimes gained or lost while readthrough persists; that readthrough is more likely to be lost at TAA and TAG stop codons; and that readthrough is under continued purifying evolutionary selection in mosquito, based on population genetic evidence. We also determine readthrough-associated gene properties that predate readthrough, and identify differences in the characteristic properties of readthrough genes between clades. We estimate more than 600 functional readthrough stop codons in mosquito and 900 in fruit fly, provide evidence of readthrough control of peroxisomal targeting, and refine the phylogenetic extent of abundant readthrough as following divergence from centipede.

Recent advances in malaria genomics and epigenomics

Malaria continues to impose a significant disease burden on low- and middle-income countries in the tropics. However, revolutionary progress over the last 3 years in nucleic acid sequencing, reverse genetics, and post-genome analyses has generated step changes in our understanding of malaria parasite (Plasmodium spp.) biology and its interactions with its host and vector. Driven by the availability of vast amounts of genome sequence data from Plasmodium species strains, relevant human populations of different ethnicities, and mosquito vectors, researchers can consider any biological component of the malarial process in isolation or in the interactive setting that is infection. In particular, considerable progress has been made in the area of population genomics, with Plasmodium falciparum serving as a highly relevant model. Such studies have demonstrated that genome evolution under strong selective pressure can be detected. These data, combined with reverse genetics, have enabled the identification of the region of the P. falciparum genome that is under selective pressure and the confirmation of the functionality of the mutations in the kelch13 gene that accompany resistance to the major frontline antimalarial, artemisinin. Furthermore, the central role of epigenetic regulation of gene expression and antigenic variation and developmental fate in P. falciparum is becoming ever clearer. This review summarizes recent exciting discoveries that genome technologies have enabled in malaria research and highlights some of their applications to healthcare. The knowledge gained will help to develop surveillance approaches for the emergence or spread of drug resistance and to identify new targets for the development of antimalarial drugs and perhaps vaccines.

Punctuated chromatin states regulate Plasmodium falciparum antigenic variation at the intron and 2 kb upstream regions

Background

Understanding the regulation mechanism of var gene expression is crucial for explaining antigenic variation in Plasmodium falciparum. Recent work observed that while all var genes produce transcripts, only a few var genes exhibit high expression levels. However, the global regulation of var expression and the relationship between epigenetic and genetic control remains to be established.

Result

We have systematically reanalyzed the existing genomic data including chromatin configurations and gene expressions; and for the first time used robust statistical methods to show that the intron and 2 kb upstream regions of each endogenous var gene always maintain high chromatin accessibility, with high potential to bind transcription factors (TFs). The levels of transcripts for different var gene family members are associated with this chromatin accessibility. Any given var gene thus shows punctuated chromatin states throughout the asexual life cycle. This is demonstrated by three independent transcript datasets. Chromatin accessibility in the var intron and 2 kb upstream regions are also positively correlated with their GC content, suggesting the level of var genes silencing might be encoded in their intron sequences. Interestingly, both var intron and 2 kb upstream regions exhibit higher chromatin accessibility when the genes have relatively lower transcription levels, suggesting a punctuated repressive function for these regulatory elements.

Conclusion

By integrating and analyzing epigenomic, genomic and transcriptomic data, our work reveals a novel distal element in var control. We found dynamic modulations of specific epigenetic marks around the var intron and distal upstream regions are involved in the general var gene expression patterns in malarial antigenic variation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3005-7) contains supplementary material, which is available to authorized users.

The mRNA-bound proteome of the human malaria parasite Plasmodium falciparum

Background

Gene expression is controlled at multiple levels, including transcription, stability, translation, and degradation. Over the years, it has become apparent that Plasmodium falciparum exerts limited transcriptional control of gene expression, while at least part of Plasmodium’s genome is controlled by post-transcriptional mechanisms. To generate insights into the mechanisms that regulate gene expression at the post-transcriptional level, we undertook complementary computational, comparative genomics, and experimental approaches to identify and characterize mRNA-binding proteins (mRBPs) in P. falciparum.

Results

Close to 1000 RNA-binding proteins are identified by hidden Markov model searches, of which mRBPs encompass a relatively large proportion of the parasite proteome as compared to other eukaryotes. Several abundant mRNA-binding domains are enriched in apicomplexan parasites, while strong depletion of mRNA-binding domains involved in RNA degradation is observed. Next, we experimentally capture 199 proteins that interact with mRNA during the blood stages, 64 of which with high confidence. These captured mRBPs show a significant overlap with the in silico identified candidate RBPs (p < 0.0001). Among the experimentally validated mRBPs are many known translational regulators active in other stages of the parasite’s life cycle, such as DOZI, CITH, PfCELF2, Musashi, and PfAlba1–4. Finally, we also detect several proteins with an RNA-binding domain abundant in Apicomplexans (RAP domain) that is almost exclusively found in apicomplexan parasites.

Conclusions

Collectively, our results provide the most complete comparative genomics and experimental analysis of mRBPs in P. falciparum. A better understanding of these regulatory proteins will not only give insight into the intricate parasite life cycle but may also provide targets for novel therapeutic strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1014-0) contains supplementary material, which is available to authorized users.

Serum stress responsive gene EhslncRNA of Entamoeba histolytica is a novel long noncoding RNA

Non coding RNAs are known to play important roles in regulating gene expression at the transcriptional and posttranscriptional levels in metazoans. There is very little information available about non coding RNAs in protists such as Entamoeba histolytica. Antisense and micro RNAs have been reported in E. histolytica, however no long non coding RNAs has been reported yet. Here, we report our findings on an in vitro serum stress-inducible gene EhslncRNA, a member of B1 transmembrane kinase family of E. histolytica. EhslncRNA encodes a transcript of 2.6 kb and sequence analysis revealed that there is no ORF >150 bp within this transcript. The transcript was found to be polyadenylated and mainly associated with monosomes in the cytoplasm under serum starvation. In normal proliferating cells this RNA is mainly present in the nucleus. The promoter element was mapped between 437 to 346 nucleotides upstream of transcriptional start site and has both positive and negative regulatory elements. Deletion of the negative element converted the promoter to serum inducible type. Oxygen and heat stress also increased expression levels of EhslncRNA. These observations suggest that EhslncRNA may be a long non coding RNA and likely to help cells withstand stressful conditions in the host.

Translational regulation in blood stages of the malaria parasite Plasmodium spp.: systems-wide studies pave the way

The malaria parasite Plasmodium spp. varies the expression profile of its genes depending on the host it resides in and its developmental stage. Virtually all messenger RNA (mRNA) is expressed in a monocistronic manner, with transcriptional activation regulated at the epigenetic level and by specialized transcription factors. Furthermore, recent systems-wide studies have identified distinct mechanisms of post-transcriptional and translational control at various points of the parasite lifecycle. Taken together, it is evident that 'just-in-time' transcription and translation strategies coexist and coordinate protein expression during Plasmodium development, some of which we review here. In particular, we discuss global and specific mechanisms that control protein translation in blood stages of the human malaria parasite Plasmodium falciparum, once a cytoplasmic mRNA has been generated, and its crosstalk with mRNA decay and storage. We also focus on the widespread translational delay observed during the 48-hour blood stage lifecycle of P. falciparum-for over 30% of transcribed genes, including virulence factors required to invade erythrocytes-and its regulation by cis-elements in the mRNA, RNA-processing enzymes and RNA-binding proteins; the first-characterized amongst these are the DNA- and RNA-binding Alba proteins. More generally, we conclude that translational regulation is an emerging research field in malaria parasites and propose that its elucidation will not only shed light on the complex developmental program of this parasite, but may also reveal mechanisms contributing to drug resistance and define new targets for malaria intervention strategies. WIREs RNA 2016, 7:772-792. doi: 10.1002/wrna.1365 For further resources related to this article, please visit the WIREs website.

Apicomplexa-specific tRip facilitates import of exogenous tRNAs into malaria parasites

The malaria-causing Plasmodium parasites are transmitted to vertebrates by mosquitoes. To support their growth and replication, these intracellular parasites, which belong to the phylum Apicomplexa, have developed mechanisms to exploit their hosts. These mechanisms include expropriation of small metabolites from infected host cells, such as purine nucleotides and amino acids. Heretofore, no evidence suggested that transfer RNAs (tRNAs) could also be exploited. We identified an unusual gene in Apicomplexa with a coding sequence for membrane-docking and structure-specific tRNA binding. This Apicomplexa protein-designated tRip (tRNA import protein)-is anchored to the parasite plasma membrane and directs import of exogenous tRNAs. In the absence of tRip, the fitness of the parasite stage that multiplies in the blood is significantly reduced, indicating that the parasite may need host tRNAs to sustain its own translation and/or as regulatory RNAs. Plasmodium is thus the first example, to our knowledge, of a cell importing exogenous tRNAs, suggesting a remarkable adaptation of this parasite to extend its reach into host cell biology.

Illuminating Parasite Protein Production by Ribosome Profiling

While technologies for global enumeration of transcript abundance are well-developed, those that assess protein abundance require tailoring to penetrate to low-abundance proteins. Ribosome profiling circumvents this challenge by measuring global protein production via sequencing small mRNA fragments protected by the assembled ribosome. This powerful approach is now being applied to protozoan parasites including trypanosomes and Plasmodium. It has been used to identify new protein-coding sequences (CDSs) and clarify the boundaries of previously annotated CDSs in Trypanosoma brucei. Ribosome profiling has demonstrated that translation efficiencies vary widely between genes and, for trypanosomes at least, for the same gene across stages. The ribosomal proteins are themselves subjected to translational control, suggesting a means of reinforcing global translational regulation.

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.

Evaluation of next generation sequencing for the analysis of Eimeria communities in wildlife

Next-generation sequencing (NGS) techniques are well-established for studying bacterial communities but not yet for microbial eukaryotes. Parasite communities remain poorly studied, due in part to the lack of reliable and accessible molecular methods to analyse eukaryotic communities. We aimed to develop and evaluate a methodology to analyse communities of the protozoan parasite Eimeria from populations of the Australian marsupial Petrogale penicillata (brush-tailed rock-wallaby) using NGS. An oocyst purification method for small sample sizes and polymerase chain reaction (PCR) protocol for the 18S rRNA locus targeting Eimeria was developed and optimised prior to sequencing on the Illumina MiSeq platform. A data analysis approach was developed by modifying methods from bacterial metagenomics and utilising existing Eimeria sequences in GenBank. Operational taxonomic unit (OTU) assignment at a high similarity threshold (97%) was more accurate at assigning Eimeria contigs into Eimeria OTUs but at a lower threshold (95%) there was greater resolution between OTU consensus sequences. The assessment of two amplification PCR methods prior to Illumina MiSeq, single and nested PCR, determined that single PCR was more sensitive to Eimeria as more Eimeria OTUs were detected in single amplicons. We have developed a simple and cost-effective approach to a data analysis pipeline for community analysis of eukaryotic organisms using Eimeria communities as a model. The pipeline provides a basis for evaluation using other eukaryotic organisms and potential for diverse community analysis studies.

Cis regulatory motifs and antisense transcriptional control in the apicomplexan Theileria parva

Background

Theileria parva is an intracellular parasite that causes a lymphoproliferative disease in cattle. It does so by inducing cancer-like phenotypes in the host cells it infects, although the molecular and regulatory mechanisms involved remain poorly understood. RNAseq data, and the resulting updated genome annotation now available for this parasite, offer an unprecedented opportunity to characterize the genomic features associated with gene regulation in this species. Our previous analyses revealed a T. parva genome even more gene-dense than previously thought, with many adjacent loci overlapping each other, not only at the level of untranslated sequences (UTRs) but even in coding sequences.

Results

Despite this compactness, Theileria intergenic regions show a pattern of size distribution indicative of monocistronic gene transcription. Three previously described motifs are conserved among Theileria species and highly prevalent in promoter regions near or at the transcription start sites. We found novel motifs at many transcription termination sites, as well as upstream of parasite genes thought to be critical for host transformation. Adjacent genes that could be regulated by antisense transcription from an overlapping transcriptional unit are syntenic between T. parva and P. falciparum at a frequency higher than expected by chance, suggesting the presence of common, and evolutionary old, regulatory mechanisms in the phylum Apicomplexa.

Conclusions

We propose a model of transcription with conserved sense and antisense transcription from a few taxonomically ubiquitous and several species-specific promoter motifs. Interestingly, the gene networks regulated by conserved promoters are themselves, in most cases, not conserved between species or genera.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2444-5) contains supplementary material, which is available to authorized users.

Upstream AUGs and upstream ORFs can regulate the downstream ORF in Plasmodium falciparum

Background

Upstream open reading frames (uORFs) and upstream AUGs (uAUGs) can regulate the translation of downstream ORFs. The AT rich genome of Plasmodium falciparum, due to the higher AT content of start and stop codons, has the potential to give rise to a large number of uORFs and uAUGs that may affect expression of their flanking ORFs.

Methods

A bioinformatics approach was used to detect uATGs associated with different genes in the parasite. To study the effect of some of these uAUGs on the expression of the downstream ORF, promoters and 5′ leaders containing uAUGs and uORFs were cloned upstream of a luciferase reporter gene. Luciferase assays were carried out in transient transfection experiments to assess the effects of uAUGs and mutations on reporter expression.

Results

The average number of uATGs and uORFs seen in P. falciparum coding sequences (CDS) is expectedly high compared to other less biased genomes. Certain genes, including the var gene family contain the maximum number of uATGs and uORFs in the parasite. They possess ~5 times more uORFs and ~4.5 times more uAUGs within 100 bases upstream of the start codons than other CDS of the parasite. A 60 bp upstream region containing three ORFs and five ATGs from var gene PF3D7_0400100 and a gene of unknown function (PF3D7_0517100) when cloned upstream of the luciferase start codon, driven by the hsp86 promoter, resulted in loss of luciferase activity. This was restored when all the ATGs present in the −60 bp were mutated to TTGs. Point mutations in the ATGs showed that even one AUG was sufficient to repress the luciferase gene.

Conclusions

Overall, this work indicates that the P. falciparum genome has a large number of uATGs and uORFs that can repress the expression of flanking ORFs. The role of AUGs in translation initiation suggests that this repression is mediated by preventing the translation initiation complex from reaching the main AUG of the downstream ORF. How the P. falciparum ribosome is able to bypass these uAUGs and uORFs for highly expressed genes remains a question for future research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-1040-5) contains supplementary material, which is available to authorized users.

Methods to Investigate the Regulatory Role of Small RNAs and Ribosomal Occupancy of Plasmodium falciparum

The genetic variation responsible for the sickle cell allele (HbS) enables erythrocytes to resist infection by the malaria parasite, P. falciparum. The molecular basis of this resistance, which is known to be multifactorial, remains incompletely understood. Recent studies found that the differential expression of erythrocyte microRNAs, once translocated into malaria parasites, affect both gene regulation and parasite growth. These miRNAs were later shown to inhibit mRNA translation by forming a chimeric RNA transcript via 5' RNA fusion with discreet subsets of parasite mRNAs. Here, the techniques that were used to study the functional role and putative mechanism underlying erythrocyte microRNAs on the gene regulation and translational potential of P. falciparum, including the transfection of modified synthetic microRNAs into host erythrocytes, will be detailed.  Finally, a polysome gradient method is used to determine the extent of translation of these transcripts. Together, these techniques allowed us to demonstrate that the dysregulated levels of erythrocyte microRNAs contribute to cell-intrinsic malaria resistance of sickle erythrocytes.

Analysis of nucleosome positioning landscapes enables gene discovery in the human malaria parasite Plasmodium falciparum

Background

Plasmodium falciparum, the deadliest malaria-causing parasite, has an extremely AT-rich (80.7 %) genome. Because of high AT-content, sequence-based annotation of genes and functional elements remains challenging. In order to better understand the regulatory network controlling gene expression in the parasite, a more complete genome annotation as well as analysis tools adapted for AT-rich genomes are needed. Recent studies on genome-wide nucleosome positioning in eukaryotes have shown that nucleosome landscapes exhibit regular characteristic patterns at the 5’- and 3’-end of protein and non-protein coding genes. In addition, nucleosome depleted regions can be found near transcription start sites. These unique nucleosome landscape patterns may be exploited for the identification of novel genes. In this paper, we propose a computational approach to discover novel putative genes based exclusively on nucleosome positioning data in the AT-rich genome of P. falciparum.

Results

Using binary classifiers trained on nucleosome landscapes at the gene boundaries from two independent nucleosome positioning data sets, we were able to detect a total of 231 regions containing putative genes in the genome of Plasmodium falciparum, of which 67 highly confident genes were found in both data sets. Eighty-eight of these 231 newly predicted genes exhibited transcription signal in RNA-Seq data, indicative of active transcription. In addition, 20 out of 21 selected gene candidates were further validated by RT-PCR, and 28 out of the 231 genes showed significant matches using BLASTN against an expressed sequence tag (EST) database. Furthermore, 108 (47 %) out of the 231 putative novel genes overlapped with previously identified but unannotated long non-coding RNAs. Collectively, these results provide experimental validation for 163 predicted genes (70.6 %). Finally, 73 out of 231 genes were found to be potentially translated based on their signal in polysome-associated RNA-Seq representing transcripts that are actively being translated.

Conclusion

Our results clearly indicate that nucleosome positioning data contains sufficient information for novel gene discovery. As distinct nucleosome landscapes around genes are found in many other eukaryotic organisms, this methodology could be used to characterize the transcriptome of any organism, especially when coupled with other DNA-based gene finding and experimental methods (e.g., RNA-Seq).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2214-9) contains supplementary material, which is available to authorized users.

Estimating mRNA lengths from Plasmodium falciparum genes by Virtual Northern RNA-seq analysis

Accurate gene models are essential for understanding parasite biology. However, transcript structure information is lacking for most parasite genes. Here, we describe "Virtual Northern" analysis of the malaria parasite Plasmodium falciparum to address this issue. RNA-seq libraries were made from size-fractionated RNA. Transcript sizes for 3052 genes were inferred from the read counts in each library. The data show that for almost half of the transcripts, the combined untranslated regions are more than twice the length of the open reading frame. Furthermore, we identified novel polycistronic, or gene overlapping, transcripts that suggest revisions to current gene models are needed.

PfAlba1: master regulator of translation in the malaria parasite

During the asexual replication cycle of the malaria parasite Plasmodium falciparum, the RNA-binding protein PfAlba1 binds and stabilizes a subset of transcripts for translation at a later time point.

Please see related Research article: http://www.genomebiology.com/2015/16/1/212

Dynamical features of the Plasmodium falciparum ribosome during translation

Plasmodium falciparum, the mosquito-transmitted Apicomplexan parasite, causes the most severe form of human malaria. In the asexual blood-stage, the parasite resides within erythrocytes where it proliferates, multiplies and finally spreads to new erythrocytes. Development of drugs targeting the ribosome, the site of protein synthesis, requires specific knowledge of its structure and work cycle, and, critically, the ways they differ from those in the human host. Here, we present five cryo-electron microscopy (cryo-EM) reconstructions of ribosomes purified from P. falciparum blood-stage schizonts at sub-nanometer resolution. Atomic models were built from these density maps by flexible fitting. Significantly, our study has taken advantage of new capabilities of cryo-EM, in visualizing several structures co-existing in the sample at once, at a resolution sufficient for building atomic models. We have discovered structural and dynamic features that differentiate the ribosomes of P. falciparum from those of mammalian system. Prompted by the absence of RACK1 on the ribosome in our and an earlier study we confirmed that RACK1 does not specifically co-purify with the 80S fraction in schizonts. More extensive studies, using cryo-EM methodology, of translation in the parasite will provide structural knowledge that may lead to development of novel anti-malarials.

The PfAlba1 RNA-binding protein is an important regulator of translational timing in Plasmodium falciparum blood stages

Background

Transcriptome-wide ribosome occupancy studies have suggested that during the intra-erythrocytic lifecycle of Plasmodium falciparum, select mRNAs are post-transcriptionally regulated. A subset of these encodes parasite virulence factors required for invading host erythrocytes, and are currently being developed as vaccine candidates. However, the molecular mechanisms that govern post-transcriptional regulation are currently unknown.

Results

We explore the previously identified DNA/RNA-binding protein PfAlba1, which localizes to multiple foci in the cytoplasm of P. falciparum trophozoites. We establish that PfAlba1 is essential for asexual proliferation, and subsequently investigate parasites overexpressing epitope-tagged PfAlba1 to identify its RNA targets and effects on mRNA homeostasis and translational regulation. Using deep sequencing of affinity-purified PfAlba1-associated RNAs, we identify 1193 transcripts that directly bind to PfAlba1 in trophozoites. For 105 such transcripts, 43 % of which are uncharacterized and 13 % of which encode erythrocyte invasion components, the steady state levels significantly change at this stage, evidencing a role for PfAlba1 in maintaining mRNA homeostasis. Additionally, we discover that binding of PfAlba1 to four erythrocyte invasion mRNAs, Rap1, RhopH3, CDPK1, and AMA1, is linked to translation repression in trophozoites whereas release of these mRNAs from a PfAlba1 complex in mature stages correlates with protein synthesis.

Conclusions

We show that PfAlba1 binds to a sub-population of asexual stage mRNAs and fine-tunes the timing of translation. This mode of post-transcriptional regulation may be especially important for P. falciparum erythrocyte invasion components that have to be assembled into apical secretory organelles in a highly time-dependent manner towards the end of the parasite’s asexual lifecycle.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0771-5) contains supplementary material, which is available to authorized users.

Evidence that the Malaria Parasite Plasmodium falciparum Putative Rhoptry Protein 2 Localizes to the Golgi Apparatus throughout the Erythrocytic Cycle

Invasion of a red blood cell by Plasmodium falciparum merozoites is an essential step in the malaria lifecycle. Several of the proteins involved in this process are stored in the apical complex of the merozoite, a structure containing secretory organelles that are released at specific times during invasion. The molecular players involved in erythrocyte invasion thus represent potential key targets for both therapeutic and vaccine-based strategies to block parasite development. In our quest to identify and characterize new effectors of invasion, we investigated the P. falciparum homologue of a P. berghei protein putatively localized to the rhoptries, the Putative rhoptry protein 2 (PbPRP2). We show that in P. falciparum, the protein colocalizes extensively with the Golgi apparatus across the asexual erythrocytic cycle. Furthermore, imaging of merozoites caught at different times during invasion show that PfPRP2 is not secreted during the process instead staying associated with the Golgi apparatus. Our evidence therefore suggests that PfPRP2 is a Golgi protein and that it is likely not a direct effector in the process of merozoite invasion.

Systems analysis of host-parasite interactions

Parasitic diseases caused by protozoan pathogens lead to hundreds of thousands of deaths per year in addition to substantial suffering and socioeconomic decline for millions of people worldwide. The lack of effective vaccines coupled with the widespread emergence of drug‐resistant parasites necessitates that the research community take an active role in understanding host–parasite infection biology in order to develop improved therapeutics. Recent advances in next‐generation sequencing and the rapid development of publicly accessible genomic databases for many human pathogens have facilitated the application of systems biology to the study of host–parasite interactions. Over the past decade, these technologies have led to the discovery of many important biological processes governing parasitic disease. The integration and interpretation of high‐throughput ‐omic data will undoubtedly generate extraordinary insight into host–parasite interaction networks essential to navigate the intricacies of these complex systems. As systems analysis continues to build the foundation for our understanding of host–parasite biology, this will provide the framework necessary to drive drug discovery research forward and accelerate the development of new antiparasitic therapies. WIREs Syst Biol Med 2015, 7:381–400. doi: 10.1002/wsbm.1311

For further resources related to this article, please visit the WIREs website.

Strand-specific RNA sequencing in Plasmodium falciparum malaria identifies developmentally regulated long non-coding RNA and circular RNA

Background

The human malaria parasite Plasmodium falciparum has a complex and multi-stage life cycle that requires extensive and precise gene regulation to allow invasion and hijacking of host cells, transmission, and immune escape. To date, the regulatory elements orchestrating these critical parasite processes remain largely unknown. Yet it is becoming increasingly clear that long non-coding RNAs (lncRNAs) could represent a missing regulatory layer across a broad range of organisms.

Results

To investigate the regulatory capacity of lncRNA in P. falciparum, we harvested fifteen samples from two time-courses. Our sample set profiled 56 h of P. falciparum blood stage development. We then developed and validated strand-specific, non-polyA-selected RNA sequencing methods, and pursued the first assembly of P. falciparum strand-specific transcript structures from RNA sequencing data. This approach enabled the annotation of over one thousand lncRNA transcript models and their comprehensive global analysis: coding prediction, periodicity, stage-specificity, correlation, GC content, length, location relative to annotated transcripts, and splicing. We validated the complete splicing structure of three lncRNAs with compelling properties. Non-polyA-selected deep sequencing also enabled the prediction of hundreds of intriguing P. falciparum circular RNAs, six of which we validated experimentally.

Conclusions

We found that a subset of lncRNAs, including all subtelomeric lncRNAs, strongly peaked in expression during invasion. By contrast, antisense transcript levels significantly dropped during invasion. As compared to neighboring mRNAs, the expression of antisense-sense pairs was significantly anti-correlated during blood stage development, indicating transcriptional interference. We also validated that P. falciparum produces circRNAs, which is notable given the lack of RNA interference in the organism, and discovered that a highly expressed, five-exon antisense RNA is poised to regulate P. falciparum gametocyte development 1 (PfGDV1), a gene required for early sexual commitment events.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1603-4) contains supplementary material, which is available to authorized users.

Modeling metabolism and stage-specific growth of Plasmodium falciparum HB3 during the intraerythrocytic developmental cycle

The human malaria parasite Plasmodium falciparum goes through a complex life cycle, including a roughly 48-hour-long intraerythrocytic developmental cycle (IDC) in human red blood cells. A better understanding of the metabolic processes required during the asexual blood-stage reproduction will enhance our basic knowledge of P. falciparum and help identify critical metabolic reactions and pathways associated with blood-stage malaria. We developed a metabolic network model that mechanistically links time-dependent gene expression, metabolism, and stage-specific growth, allowing us to predict the metabolic fluxes, the biomass production rates, and the timing of production of the different biomass components during the IDC. We predicted time- and stage-specific production of precursors and macromolecules for P. falciparum (strain HB3), allowing us to link specific metabolites to specific physiological functions. For example, we hypothesized that coenzyme A might be involved in late-IDC DNA replication and cell division. Moreover, the predicted ATP metabolism indicated that energy was mainly produced from glycolysis and utilized for non-metabolic processes. Finally, we used the model to classify the entire tricarboxylic acid cycle into segments, each with a distinct function, such as superoxide detoxification, glutamate/glutamine processing, and metabolism of fumarate as a byproduct of purine biosynthesis. By capturing the normal metabolic and growth progression in P. falciparum during the IDC, our model provides a starting point for further elucidation of strain-specific metabolic activity, host-parasite interactions, stress-induced metabolic responses, and metabolic responses to antimalarial drugs and drug candidates.

Comprehensive evaluation of Toxoplasma gondii VEG and Neospora caninum LIV genomes with tachyzoite stage transcriptome and proteome defines novel transcript features

Toxoplasma gondii is an important protozoan parasite that infects all warm-blooded animals and causes opportunistic infections in immuno-compromised humans. Its closest relative, Neospora caninum, is an important veterinary pathogen that causes spontaneous abortion in livestock. Comparative genomics of these two closely related coccidians has been of particular interest to identify genes that contribute to varied host cell specificity and disease. Here, we describe a manual evaluation of these genomes based on strand-specific RNA sequencing and shotgun proteomics from the invasive tachyzoite stages of these two parasites. We have corrected predicted structures of over one third of the previously annotated gene models and have annotated untranslated regions (UTRs) in over half of the predicted protein-coding genes. We observe distinctly long UTRs in both the organisms, almost four times longer than other model eukaryotes. We have also identified a putative set of cis-natural antisense transcripts (cis-NATs) and long intergenic non-coding RNAs (lincRNAs). We have significantly improved the annotation quality in these genomes that would serve as a manually curated dataset for Toxoplasma and Neospora research communities.

Genome-wide regulatory dynamics of translation in the Plasmodium falciparum asexual blood stages

The characterization of the transcriptome and proteome of Plasmodium falciparum has been a tremendous resource for the understanding of the molecular physiology of this parasite. However, the translational dynamics that link steady-state mRNA with protein levels are not well understood. In this study, we bridge this disconnect by measuring genome-wide translation using ribosome profiling, through five stages of the P. falciparum blood phase developmental cycle. Our findings show that transcription and translation are tightly coupled, with overt translational control occurring for less than 10% of the transcriptome. Translationally regulated genes are predominantly associated with merozoite egress functions. We systematically define mRNA 5′ leader sequences, and 3′ UTRs, as well as antisense transcripts, along with ribosome occupancy for each, and establish that accumulation of ribosomes on 5′ leaders is a common transcript feature. This work represents the highest resolution and broadest portrait of gene expression and translation to date for this medically important parasite.

DOI:http://dx.doi.org/10.7554/eLife.04106.001

The genome of an organism includes all of the genes or information necessary to build, maintain, and replicate that organism. However, cells with the same genome—such as a skin cell and a liver cell from the same person—can look and behave very differently depending on which of the genes in their genomes they express, and to what extent.

For a gene to be expressed, its DNA is ‘transcribed’ to make an RNA molecule, which is then ‘translated’ to make a protein. Efforts to measure the transcription and translation processes in diseased cells, or in the microorganisms that cause infections, may lead to new treatments and preventative medicines. Such work is currently ongoing in the global effort to treat and prevent malaria.

Malaria is both preventable and curable, yet over 600,000 people are estimated to die from this disease each year. The disease is caused by a single-celled parasite called Plasmodium. Mosquitoes carry the parasites in their salivary glands, and when a mosquito bites a human, these parasites are injected into the bloodstream with the mosquito's saliva. Plasmodium parasites then travel to and infect the liver, before bursting out of this tissue into the bloodstream. Here, the parasites infect red blood cells and undergo rounds of replication during which the symptoms of the disease are manifested. It is also during this bloodstream phase that parasites can develop into forms capable of infecting another mosquito and continuing the transmission cycle.

The genes, RNA molecules, and proteins of the Plasmodium falciparum parasite—which causes the most serious cases of malaria in humans—have been cataloged to better understand the biology of this parasite. However, the processes that control how, and when, an RNA transcript is translated into a protein are not well understood.

Now Caro et al. have uncovered which RNA molecules are being translated, and by how much, during Plasmodium development within the blood. The transcription and translation of genes in this parasite were found to be tightly linked processes; the expression of only a few genes was controlled more by the translation process than by transcription. These translationally regulated genes were found mainly to be those that encode proteins involved in the parasite's exit from the red blood cells and spread throughout the bloodstream.

Caro et al. discovered that genetic regulation of the malaria parasite resembles a preset genetic program, rather than a system that responds to changes and external signals. As such, these findings suggest that targeting such a genetic program within Plasmodium and preventing its implementation could prove an effective strategy to curb the spread of malaria.

DOI:http://dx.doi.org/10.7554/eLife.04106.002

Multiple dimensions of epigenetic gene regulation in the malaria parasite Plasmodium falciparum: gene regulation via histone modifications, nucleosome positioning and nuclear architecture in P. falciparum

Plasmodium falciparum is the most deadly human malarial parasite, responsible for an estimated 207 million cases of disease and 627,000 deaths in 2012. Recent studies reveal that the parasite actively regulates a large fraction of its genes throughout its replicative cycle inside human red blood cells and that epigenetics plays an important role in this precise gene regulation. Here, we discuss recent advances in our understanding of three aspects of epigenetic regulation in P. falciparum: changes in histone modifications, nucleosome occupancy and the three-dimensional genome structure. We compare these three aspects of the P. falciparum epigenome to those of other eukaryotes, and show that large-scale compartmentalization is particularly important in determining histone decomposition and gene regulation in P. falciparum. We conclude by presenting a gene regulation model for P. falciparum that combines the described epigenetic factors, and by discussing the implications of this model for the future of malaria research.

Translational regulation during stage transitions in malaria parasites

The complicated life cycle of the malaria parasite involves a vertebrate host and a mosquito vector, and translational regulation plays a prominent role in orchestrating the developmental events in the two transition stages: gametocytes and sporozoites. Translational regulation is executed in both global and transcript-specific manners. Plasmodium uses a conserved mechanism involving phosphorylation of eIF2α to repress global protein synthesis during the latent period of sporozoite development in the mosquito salivary glands. Transcript-specific translational regulation is achieved by a network of RNA-binding proteins (RBPs), among which the Dhh1 RNA helicase DOZI and Puf family RBPs are by far the best studied in Plasmodium. While the DOZI complex defines a new P granule with a role in protecting certain gametocyte mRNAs from degradation, the Puf proteins appear to repress expression of mRNAs in both gametocytes and sporozoites. These examples underscore the significance of translational regulation in Plasmodium development.

Assessment of selective mRNA translation in mammalian cells by polysome profiling

Regulation of protein synthesis represents a key control point in cellular response to stress. In particular, discreet RNA regulatory elements were shown to allow to selective translation of specific mRNAs, which typically encode for proteins required for a particular stress response. Identification of these mRNAs, as well as the characterization of regulatory mechanisms responsible for selective translation has been at the forefront of molecular biology for some time. Polysome profiling is a cornerstone method in these studies. The goal of polysome profiling is to capture mRNA translation by immobilizing actively translating ribosomes on different transcripts and separate the resulting polyribosomes by ultracentrifugation on a sucrose gradient, thus allowing for a distinction between highly translated transcripts and poorly translated ones. These can then be further characterized by traditional biochemical and molecular biology methods. Importantly, combining polysome profiling with high throughput genomic approaches allows for a large scale analysis of translational regulation.