Abstract
Malaria is caused by unicellular apicomplexan parasites of the genus Plasmodium, which includes the major human parasite Plasmodium falciparum. The complex cycle of the malaria parasite in both mosquito and human hosts has been studied extensively. There is tight control of gene expression in each developmental stage, and at every level of gene synthesis: from RNA transcription, to its subsequent translation, and finally post-translational modifications of the resulting protein. Whole-genome sequencing of P. falciparum has laid the foundation for significant biological advances by revealing surprising genomic information. The P. falciparum genome is extremely AT-rich (∼80%), with a substantial portion of genes encoding intragenic polyadenosine (polyA) tracks being expressed throughout the entire parasite life cycle. In most eukaryotes, intragenic polyA runs act as negative regulators of gene expression. Recent studies have shown that translation of mRNAs containing 12 or more consecutive adenosines results in ribosomal stalling and frameshifting; activating mRNA surveillance mechanisms. In contrast, P. falciparum translational machinery can efficiently and accurately translate polyA tracks without activating mRNA surveillance pathways. This unique feature of P. falciparum raises interesting questions: (1) How is P. falciparum able to efficiently and correctly translate polyA track transcripts, and (2) What are the specifics of the translational machinery and mRNA surveillance mechanisms that separate P. falciparum from other organisms? In this review, we analyze possible evolutionary shifts in P. falciparum protein synthesis machinery that allow efficient translation of an AU rich-transcriptome. We focus on physiological and structural differences of P. falciparum stage specific ribosomes, ribosome-associated proteins, and changes in mRNA surveillance mechanisms throughout the complete parasite life cycle, with an emphasis on the mosquito and liver stages.
Introduction
Plasmodium spp. has been in existence long before humans were on Earth, with an estimated origin of malaria-causing parasites appearing around 165 million years ago. Consequently, mosquitos and malaria had millions of years to co-evolve before either ever interacted with humans (Winegard, 2019). The infection of humans occurred evolutionarily recently, and probably with multiple Plasmodium parasite species. P. falciparum and P. vivax established themselves as a major malaria causing species. P. falciparum a most virulent agent in human malaria began speciation around 50,000 years ago followed by the population bottleneck around 5000 years ago but higher level of genetic diversity suggests that P. vivax is older (Loy et al., 2018; Otto et al., 2018). P. malariae, P. ovale, and rare cases of P. knowlesi were also reported in human hosts. From the mid-19th century onward, malaria reached its global limits and exacted immensely high numbers in sickness and death. While increased malaria prevention and control treatments have reduced the health burden of malaria, there are still 219 million cases of infection per year resulting in a 435,000 deaths (World Health Organization [WHO], 2018). The complex cycle of the malaria parasite in both mosquito and human hosts has been studied extensively (Figure 1). In each of these life cycle stages, gene expression is tightly controlled (Le Roch et al., 2004; Shock et al., 2007; Hughes et al., 2010; Sorber et al., 2011; Bunnik and Le Roch, 2013; Caro et al., 2014; Vembar et al., 2016a; Lu et al., 2017).
FIGURE 1
It took years of laborious efforts to sequence P. falciparum genome (Kooij et al., 2006). Sequences of single or multiple chromosomes as well as complete genome were reported over the course of 4 years (Gardner et al., 1998, 2002a,b; Bowman et al., 1999; Hall et al., 2002; Hyman et al., 2002). The high AT-content of the genome made gap closure in sequences extremely difficult. However, long-read, single molecule, real-time sequencing allowed for complete telomere-to-telomere de novo assembly of the P. falciparum genome thereby overcoming the problems associated with next generation sequencing of AT-rich genomes (Vembar et al., 2016b). The consequence of AT-richness is the presence of extended tracts of As, Ts, and TAs in introns and intergenic regions (Glöckner, 2000; Szafranski et al., 2005) as well as unusually high number of genes containing coding polyadenosine (polyA) repeats compared to the other species (Habich, 2016; Djuranovic et al., 2018). Repetitions of 12 or more adenosine nucleotides in gene coding sequences, so-called polyA tracks, were recently found to act as negative gene regulation motifs at the level of mRNA translation in all tested organisms (Arthur et al., 2015; Koutmou et al., 2015). Consequently, polyA tracks have been evolutionarily preserved in a select set of genes, but are generally selected against in overall gene coding sequences (Arthur et al., 2015).
Recent analysis of 250 eukaryotic genomes found a median of 2% of transcripts with polyA tracks (Habich, 2016). However, Plasmodium species represent an exception to this rule. The percentage of polyA carrying transcripts in the genome exceeds 60% for most Plasmodium spp., including P. falciparum (64%) (Djuranovic et al., 2018). The pervasive ribosomal stalling and frameshifting found on polyA tracks in other eukaryotes (Arthur et al., 2015; Koutmou et al., 2015; Tournu et al., 2019) would make it almost impossible for the majority of Plasmodium proteins to be efficiently and correctly synthesized. However, global studies of Plasmodium protein composition (Florens et al., 2002; Silvestrini et al., 2010) and protein synthesis (Le Roch et al., 2004; Bunnik et al., 2013; Caro et al., 2014) do not show any reduction in either the protein or mRNA abundances of polyA track genes. This suggests that both ribosomal stalling and frameshifting in Plasmodium are resolved by adaptations in protein synthesis and mRNA quality control systems. In this review, we will discuss how the extreme AT-rich genome of malaria-causing parasite promoted special features in P. falciparum ribosomes to enable translation of polyA tracks throughout the complete life cycle. Additionally, genomic changes and parasitic environment have also influenced variation in mRNA surveillance mechanism within the organism resulting in divergence from other Eukaryotes.
Evolution of At-Richness in P. falciparum
Extremes in genomic base composition toward GC- or AT-richness exist in all domains of life (Sueoka, 1962; Wernegreen and Funk, 2004; Zilversmit et al., 2010; Wu et al., 2012). The extent of these extremes in nucleotide composition is limited by the necessity of all 20 amino acids and the subsequent requirement of all four nucleotides to encode them. As such, long homopolymeric amino acid repeats appear to be a characteristic of genomes with either bias (Glöckner, 2000; Albà et al., 2007; Muralidharan et al., 2011). Harboring either extreme AT- or GC-richness affects genomic structure, stability, transcriptome, and codon bias of organisms (Wu et al., 2012). As seen in Table 1, the P. falciparum mean AT-richness of around 80% appears to be one of the highest in all Eukaryotes (Pollack et al., 1982; Musto et al., 1999; Gardner et al., 2002a; Videvall, 2018). Surprisingly, the higher AT-content of the P. falciparum genome cannot be fully explained by increased AT-richness in intergenic regions, but rather by contributions of AT-richness in both coding 76.22% (Table 1) and non-coding genome 90% (Gardner et al., 2002a). Overall, gene organization patterns in P. falciparum are not influenced by the AT-bias (Glöckner, 2000; Szafranski et al., 2005; Djuranovic et al., 2018). However, what distinguishes Plasmodium species from other AT-rich organisms is distribution of consecutive adenosine nucleotides resulting in unusually high percentage of polyA track genes (Table 2). The genomes of P. falciparum and related Plasmodium species have apparently evolved independently to reach extreme AT-bias (Table 2). Interestingly, while the two groups of Plasmodium species can be separated based on their AT-genomic content (median of 75% versus a median of 55% AT-richness), both groups accommodate a considerable amount of polyA tracks within the coding regions (Djuranovic et al., 2018).
TABLE 1
| Organism | CDS % AT richness | PolyA track genes |
| Plasmodium falciparum | 76.22% | 63.54% |
| Plasmodium reichenowi | 75.93% | 62.93% |
| Dictyostelium discoideum | 72.57% | 20.96% |
| Tetrahymena thermophila | 72.50% | 28.19% |
| Saccharomyces cerevisiae | 60.39% | 5.51% |
| Plasmodium knowlesi | 59.77% | 41.42% |
| Caenorhabditis elegans | 57.94% | 1.93% |
| Plasmodium vivax | 53.51% | 38.85% |
| Drosophila melanogaster | 50.66% | 1.10% |
| Pan troglodytes | 50.60% | 2.17% |
| Homo sapiens | 49.98% | 1.40% |
| Trypanosoma brucei | 49.19% | 2.89% |
| Trypanosoma cruzi | 46.83% | 2.74% |
| Toxoplasma gondii | 42.88% | 1.01% |
| Leishmania donovani | 37.63% | 0.19% |
| Leishmania major | 37.52% | 0.05% |
| Leishmania infantum jpcm5 | 37.50% | 0.12% |
| Acanthamoeba castellanii str neff | 37.06% | 0.11% |
| Emiliania huxleyi | 31.38% | 0.10% |
| Aureococcus anophagefferens | 29.37% | 0.33% |
Comparison of AT-richness and polyA track gene ratios over selected Eukaryotic species.
The coding region AT-richness from a relevant selection of organisms with high, moderate, and low AT-content was compiled from Habich (2016) and Videvall (2018) and sorted in descending order.
TABLE 2
| Organism | CDS % AT richness | PolyA track genes |
| Plasmodium gallinaceum | 78.81% | 71.77% |
| Plasmodium relictum | 78.43% | 79.18% |
| Plasmodium berghei | 76.26% | 68.26% |
| Plasmodium yoelii 17x | 77.03% | 64.96% |
| Plasmodium falciparum | 76.22% | 63.54% |
| Plasmodium chabaudi | 74.46% | 63.37% |
| Plasmodium reichenowi | 75.93% | 62.93% |
| Plasmodium vinckei petteri | 74.91% | 62.58% |
| Plasmodium gaboni | 77.56% | 61.76% |
| Plasmodium falciparum camp malaysia | 76.71% | 61.02% |
| Plasmodium falciparum nf54 | 76.55% | 60.60% |
| Plasmodium falciparum fch 4 | 76.67% | 60.57% |
| Plasmodium falciparum ugt5 1 | 76.49% | 60.45% |
| Plasmodium falciparum santa lucia | 76.75% | 60.42% |
| Plasmodium falciparum palo alto uganda | 76.56% | 60.33% |
| Plasmodium falciparum nf135 5 c10 | 76.56% | 59.98% |
| Plasmodium falciparum malips096 e11 | 76.50% | 59.81% |
| Plasmodium falciparum 7g8 | 76.65% | 59.72% |
| Plasmodium yoelii yoelii | 75.20% | 51.06% |
| Plasmodium knowlesi | 59.77% | 41.42% |
| Plasmodium knowlesi strain h | 59.76% | 41.42% |
| Plasmodium vivax | 53.51% | 38.85% |
| Plasmodium cynomolgi strain b | 57.91% | 33.83% |
| Plasmodium inui san antonio 1 | 56.37% | 31.12% |
Comparison of AT-Richness and polyA track gene ratios over selected Plasmodium species.
Plasmodium spp. coding region AT-content and the ratio of polyA affected transcripts was collected (Habich, 2016; Videvall, 2018). The data are organized in the table to demonstrate a separation of two groups with high and low coding region AT-content and subsequently the number of polyA track containing transcripts. The separation notably occurs along the line of geographic region with the high AT-content organisms being predominantly found in Africa and the low AT-content group in Asia, Southeast Asia, and Latin America. However, the low AT-content group still exceeds that of most organisms.
Perhaps just as interesting as the consequences of genomic base composition biases are the factors driving it. Previous studies in P. falciparum were unable to conclude the primary role of homopolymeric amino acid repeats in the parasite proteome (Muralidharan et al., 2011; Muralidharan and Goldberg, 2013). Nutrient availability to intracellular parasites – as well as endosymbionts – appears to be a major factor in driving AT-richness, particularly nitrogen availability (Seward and Kelly, 2016; Dietel et al., 2019). De novo synthesis of nucleotides comes at great metabolic expense, especially regarding G+C nucleotides (Dietel et al., 2019). A+T nucleotides are less metabolically costly to create and tend to be more abundant. Consequently, A+T nucleotides are easier to scavenge, even in intracellular environments where nutrients may not be readily available. In the case of P. falciparum, where de novo synthesis of purines does not occur, the parasites must rely upon purine scavenging and salvage pathways (Ting et al., 2005; El Bissati et al., 2006; Quashie et al., 2008). Conversely, pyrimidine de novo synthesis occurs using glutamine and aspartic acid precursors. This appears to be the main source for these nucleotides, with the folate pathway being required for thymidine production (Sherman, 1979; Cassera et al., 2011; Hamilton et al., 2017). However, unlike other intracellular organisms referenced above, the intracellular environment for P. falciparum is not necessarily nutrient poor, but perhaps nutrient selective; particularly prior to parasite augmentation of the host cell. While P. falciparum has multiple means by which amino acids are obtained, much of its initial amino acid supply is from proteolysis of human host red blood cell hemoglobin (Leiriao et al., 2004; Liu et al., 2006; Babbitt et al., 2012).
This brings us to a second major contributor of AT-richness in intracellular organisms: oxidative stress. Reactive nitrogen (RNS) and reactive oxygen species (ROS) generate oxidative stress resulting in 8-oxoguanine production via guanine oxidation. If left unrepaired in DNA, 8-oxoG is able to pair with adenosine; ultimately causing a G:C to T:A conversion. Compounding the process, hemoglobin degradation produces free heme and H2O2, which generates further oxidative stress for the parasite (Becker et al., 2004). Additionally, NO and other RNS species may be important factors in the soluble heme-hemozoin equilibration (Ostera et al., 2011). Interestingly, another erythrocytic parasite from Apicomplexa phylum, Babesia microti, does not degrade hemoglobin and has a considerably less AT-rich genome (61.02%) and polyA tracks (2.17% of genes with polyA tracks) compared to P. falciparum (Cornillot et al., 2012; Djuranovic et al., 2018). Although Plasmodium spp. does supply some of its own antioxidants to cope with oxidative assault, the higher than expected G:C to T:A conversion in the organism suggests a lack of full compensation by the biochemical and/or DNA repair safeguards (Hamilton et al., 2017). While 8-oxoG could potentially result in AT-richness imprinted in the DNA sequence, it causes more problems when found in RNA (Simms and Zaher, 2016). The oxidative lesion and incorporation of 8-oxoG in mRNAs reduces the rate of peptide-bond formation by more than three orders of magnitude (Simms et al., 2014). The effect of 8-oxoG nucleotides in mRNAs is independent of its position within the codon, results in stalling of the translational machinery, and finally activation of No-Go decay mRNA surveillance mechanisms (Simms et al., 2014). As such, the presence of oxidative stress may have driven both an increase in genomic AT-richness and changes in mRNA surveillance mechanisms of P. falciparum; which are discussed further below.
AT-richness itself appears to provide a feedback loop in the parasite with its increased indel rates, which are thought to be due to DNA replication slippage on AT repeats. These AT tracks provide amplicon breakpoints for copy number variant (CNV) alteration via non-allelic homologous repair-like mechanism that can be advantageous in altering resistance gene CNV numbers (Guler et al., 2013; Hamilton et al., 2017; Huckaby et al., 2019). Altogether, metabolic and biochemical factors continuously drive the parasite genome toward AT richness, which, in turn, drives indels that potentiate genomic plasticity providing an overall platform for relatively rapid adaptive evolution in the parasite. Unarguably, these factors necessitate increased fidelity in DNA replication and RNA transcription. While the exact details specific to Plasmodium spp. evolutionary adaptation toward an AT-rich genome, unique codon biases, and polyA encoded lysine stretches remains to be explored, the role ribosomes play as influential factors in this process is certain.
The rRNA and Specialized Ribosomes of Plasmodium
The most abundant genes in cells and genomes from bacteria to eukaryotes are those encoding ribosomal RNA. Ribosomal RNA genes in eukaryotic cells form clusters with a highly repetitive structure. S. cerevisiae, a single cell organism, has 150 rDNA repeats in one cluster on chromosome XII, while human cells contain five clusters of approximately 70 rDNA repeats on five different chromosomes (Sakai et al., 1995). The organization of rDNA genes in clusters is conserved among most of the eukaryotic organisms (Kobayashi, 2014). Transcription of these clusters is highly coordinated to meet the huge demand for ribosomes, which occupy ∼50% of the total protein mass in a cell (Warner, 1999). Plasmodium genomes, however, have only 4–8 single copy rDNA units that are encoded on different chromosomes (Gunderson et al., 1987; Waters et al., 1989; Li et al., 1997). Such a small number of rDNA copies throughout the genome is seen elsewhere only in bacteria. E. coli has seven ribosomal RNA genes spread over its circular genome and well positioned in the regions near an origin of replication. This arrangement in E. coli enables maximum ribosomal RNA transcription while preventing possible collisions between replication forks and transcription machinery (Ellwood and Nomura, 1982). Thus, while most of the other organisms have optimized ribosome production, how the malaria-causing parasite produces its significant ribosome numbers is still unknown. It might be possible that massive DNA replication that occurs throughout its lifecycle (during shizogeny) in both hosts may accommodate the rRNA production requirements.
Besides this unusual rDNA arrangement, malaria parasites are pioneers in the new era of specialized ribosomes (Walliker et al., 1987; McCutchan et al., 1988; Waters et al., 1989; Velichutina et al., 1998; Xue and Barna, 2012). Plasmodium spp. has structurally distinct, stage-specific ribosomes and are the most well-known case of rRNA heterogeneity (McCutchan et al., 1988). The difference in sequence and expression profile during the life cycle classified them into A-type (asexual stage specific) and S-type (sporozoite specific) in the majority of Plasmodium species, including P. falciparum; with P. vivax having a third O-type rRNA (Li et al., 1997). The A-type is present in the liver and blood stage and S-type is sporozoite specific rRNA type that emerges during the mosquito stage and ends during the parasite development in hepatocytes (Zhu et al., 1990). Here, we will focus on the process by which the ribosome types switch and whether ribosomes with distinct rRNA play a selective role in the mRNAs they translate.
Plasmodium spp. have adapted to translation in two different hosts. This requires translation optimization at two distinct temperatures, one of which can be highly variable depending on the mosquito environment. Even though one would think that changes in temperature and hosts would be the reason for development of different rRNAs, the presence of A-type during the early mosquito stage and S-type during early liver stage does not support that idea (Fang and McCutchan, 2002). The rRNA sets are not expressed in an exclusive and binary (on/off) fashion, but more as a dynamic, heterogeneous population whereby one subtype, A or S, is the more dominant rRNA type in a particular lifecycle stage. While the idea of a thermoregulatory nature of the rRNA units has been explored earlier in P. berghei, rodent malaria, it has not been followed since (Fang and McCutchan, 2002). P. berghei, contains four distinct copies off the rRNA (A, B, C, D) and they are divided into A-type (A and B) and S-type (C and D). A single copy of the S-type gene, C or D was sufficient for life cycle completion, which only affected the parasite fitness. The group was unable to disrupt both S-type genes simultaneously; nor could they disrupt either of A-type genes (van Spaendonk et al., 2001). Interestingly, authors noticed growth retardation in oocyst development, which was more pronounced in D-unit disruption rather than in C-unit (van Spaendonk et al., 2001). Such difference could be explained by difference in ribosomal levels stemming from different transcriptional levels of C- and D-units or functional diversity of C- and D-unit containing ribosomes (Xue and Barna, 2012; Mills and Green, 2017). The disruption of specific S-type rRNA is also associated with oocyst development defects in the second rodent parasite P. yoelii (Qi et al., 2015). Finally, van Spaendonk et al. (2001) note a lack in differences between core catalytic components (e.g., GTPase center) of the ribosome large subunit in P. berghei that were previously described in P. falciparum (Velichutina et al., 1998). These results among species of Plasmodium potentiate the question of some aspect of ribosomal specialization (Vembar et al., 2016a).
Previous bacterial work has shown changes in rRNA operon expression in response to stress, resulting in phenotypic changes (Kurylo et al., 2018). The change in Plasmodium spp. rRNA population dynamics in response to environmental stress from host transfer is reminiscent of the bacterial changes in rRNA operon expression. However, whether changes in ratios of Plasmodium spp. rRNA types drive phenotypic changes is still unknown. Ostensibly, the ribosomes share the same repertoire of ribosomal proteins. RNAseq data shows that while ribosomal protein gene transcription as a whole is fairly persistent throughout the complete life cycle of P. falciparum, oscillations in their overall expression pattern match that of stages with increased protein synthesis (Figure 2). This does not exclude the highly sought-after notion that a specific set of ribosomes may be optimized for specific mRNA substrates or cell populations that may also exist in Plasmodium spp. A recent study in zebrafish showed that embryos have different subtypes of 5.8S, 18S, and 28S rRNAs, creating similar ribosome diversity seen in Plasmodium cells (Locati et al., 2017). In silico data have shown that the expanded regions of 18S subunit expressed in zebrafish embryos may preferentially bind maternal transcripts when compared to somatic subtypes (Locati et al., 2017). Similarly, a shift in the expression of 16S rRNA ribosome variants created populations of E. coli cells that accommodated functional differences in tetracycline binding (Kurylo et al., 2018). As was mentioned before, the rRNA heterogeneity that was mostly known in Plasmodium parasites (Gunderson et al., 1987; Waters et al., 1989; Zhu et al., 1990; Rogers et al., 1996; Xue and Barna, 2012) is now recognized in other organisms (Locati et al., 2017; Kurylo et al., 2018). However, the role of different Plasmodium rRNAs as a response to different environmental conditions is still not defined.
FIGURE 2
Plasmodium Ribosomes, Polya and Poly-Lysine Sequences
Regardless of the host, all Plasmodium spp. rRNA types must contend with the translation of unusually high AU-content and long-coding polyA stretches in mRNAs. RNA-seq data (Le Roch et al., 2004; Shock et al., 2007; Bunnik et al., 2013; Guler et al., 2013; Caro et al., 2014) indicates that the mRNA levels of genes containing polyA stretches follows the same trend as the general gene expression for all stages in both hosts (Figures 3A–D). We can conclude that both types of ribosomes expressed in both hosts have features allowing efficient translation of transcripts containing long, coding polyA tracks. This indicates that P. falciparum ribosomes have higher fidelity during translation of polyA sequences and are able to accommodate long polybasic peptides coming through their protein-exit channel. Previous ribosome mutagenesis studies in S. cerevisiae suggested functional differences in the GTPase centers of P. falciparum A- and S-type ribosomes (Velichutina et al., 1998). Despite the differences in yeast viability and growth rates, chimeric yeast ribosomes with either Plasmodium’s A- or S-type GTPase centers exhibited increased translational accuracy (Velichutina et al., 1998). Even though there are stage-specific ribosomes, there is a group of genes that is present in human and mosquito that contain polyA tracks (Figure 4). More recently it was also shown that the P. falciparum ribosomes have been altered to accommodate the poly-lysine patches that are prolific throughout the proteome (Djuranovic et al., 2018). To allow these low-complexity, homopolymeric and polybasic amino acid repeats, the parasite ribosome exit channel has been altered by increasing the channel size at key bottle necks, as well as a reduction in the hydrophobicity patches typically seen in bacterial, yeast, or human ribosomes (Djuranovic et al., 2018).
FIGURE 3
FIGURE 4
Ribosome profiling and biochemical assays suggest an increased or modified fidelity such that parasite ribosomes do not stall or frameshift on polyA tracks (Djuranovic et al., 2018). The mechanism of this altered fidelity may result from not only modification of the ribosomal RNA sequence, but also via changes to key protein components of ribosomes. Two P. falciparum ribosome cryoEM structures suggest a reduced or lost interaction of the receptor for activated C kinase 1 (RACK1) to Plasmodium ribosomes (Wong et al., 2014; Sun et al., 2015). RACK1 has been established as an integral ribosomal scaffold protein (Sengupta et al., 2004). Beside other non-ribosome associated functions, RACK1 was found to be important for cap-dependent translation initiation, IRES-mediated translation, and site-specific translation (Majzoub et al., 2014). RACK1 also contributes to the translation arrest that is induced by translation of polyA sequences (Dimitrova et al., 2009; Kuroha et al., 2010), CGA-CGA codons in yeast (Wolf and Grayhack, 2015), or runs of consecutive basic amino-acid (Kuroha et al., 2010). Stalls on polyA tracks can be resolved in mammalian cells by deletion of RACK1 and ZNF598, thus enabling read-through of stall-inducing sequences (Garzia et al., 2017; Juszkiewicz and Hegde, 2017; Sundaramoorthy et al., 2017). S. cerevisiae ribosomes lacking the RACK1 homolog Asc1 are able to translate through the CGA-CGA stalling sequences and increase normally attenuated protein output (Wolf and Grayhack, 2015). The increase in amount of synthesized protein from CGA-CGA sequences is a consequence of overall reduced elongation rates of yeast ribosomes that lack Asc1 (Tesina et al., 2019). Slower elongation rates may also influence cellular responses to ribosome pausing. The position of RACK1/Asc1 near the mRNA exit channel on the ribosome could be important in sensing ribosome collisions that lead to activation of ribosome rescue and mRNA surveillance pathways (Kim et al., 2014; Simms and Zaher, 2016; Tesina et al., 2019). The fact that Plasmodium ribosomes lack interaction with PfRACK1 could be beneficial for translation of polyA tracks into poly-lysine runs. However, based on previous conclusions concerning the role of RACK1/Asc1 in correct reading frame maintenance during translation of stalling sequences, the majority of polyA coding sequences in malaria parasites would be predicted to have multiple frameshifted protein products.
Previous studies (Lu and Deutsch, 2008; Kuroha et al., 2010; Brandman et al., 2012) proposed that stalling during the translation of polyA tracks is due to synthesis of the poly-lysine rich nascent peptide. Electrostatic interactions of the polybasic peptide and the peptide exit tunnel in the ribosome would elicit ribosomal stalling (Lu and Deutsch, 2008). Recent studies revealed that an mRNA-mediated mechanism is directly contributing to stalling (Arthur et al., 2015; Koutmou et al., 2015; Tesina et al., 2019). Consecutive adenosines are engaged by the ribosome decoding center nucleotides, are stabilized on both sides by rRNA base stacking interactions (Tesina et al., 2019), and adopt a helical conformation typical for single stranded polyA stretches (Tang et al., 2019). PolyA tracks are highly efficient at causing ribosome stalling, and the inhibitory conformation of polyA mRNA bases can further contribute to a polyA-mediated stalling mechanism. This conclusion is in line with the previous observations that consecutive AAG codons are less efficient at causing stalling than AAA codons (Arthur et al., 2015; Koutmou et al., 2015), despite coding for the same amino acid. Altogether, the charge and conformation of the poly-lysine nascent chain in conjunction with the stacked polyA mRNA nucleotides in the decoding center of the ribosome contribute to the overall stalling mechanism (Tesina et al., 2019). P. falciparum ribosomes are again the exception to this rule. The nucleotides that make stacking interactions with polyA repeats are conserved in P. falciparum ribosomes. However, both endogenous transcripts and reporter sequences with long runs of polyA tracks are efficiently translated by Plasmodium (Lacsina et al., 2011; Bunnik et al., 2013; Djuranovic et al., 2018). Thus, in order to adapt to polyA track translation for production of the polybasic and homopolymeric lysine repeats, the malaria-causing parasite has altered the sequence of its rRNA, its ribosome structure, its ribosomal proteins, and its mRNA translation quality control pathways.
mRNA Surveillance Pathways in Au-Rich Transcriptome of P. falciparum
The core elements for mRNA translation are highly conserved in Plasmodium spp. (Vembar et al., 2016a). The unique features involving protein synthesis in Plasmodium, such as different types of ribosomes in different life cycle stages, were noticed even before genome sequencing of the malaria parasite (Gunderson et al., 1987; Zhu et al., 1990; Rogers et al., 1996; van Spaendonk et al., 2001). However, the presence of an unusual number of mRNA binding proteins and the absence of some elements of mRNA surveillance mechanism were noticed upon completion of the P. falciparum genome sequence (Gunderson et al., 1987; Waters et al., 1989; Rogers et al., 1996; van Spaendonk et al., 2001; Le Roch et al., 2004; Bunnik and Le Roch, 2013; Cui et al., 2015; Reddy et al., 2015; Lu et al., 2017). A recent review elaborated on the translational regulation in blood-stages of malaria parasites (Vembar et al., 2016a). They focused on cytoplasmic mRNA translation and the fate of mRNAs: decoding of the mRNA messages by the 80S ribosomes, degradation of mRNAs by exo- or endo-nucleases (mRNA decay), and sequestration of mRNAs by protection from mRNA decay or by inhibition of translation. We focus here on the mechanism of activation of mRNA surveillance pathways by aberrant mRNAs in the context of unusual AU-richness and abundance of polyA tracks in Plasmodium transcriptome.
Eukaryotic cells have developed mechanisms to protect themselves from the production of the possible toxic proteins due to aberrant mRNA translation events. There are three mRNA quality control systems for translational errors in eukaryotes: Non-sense mediated decay (NMD), No-Go decay (NGD), and Non-Stop decay (NSD). NMD targets transcripts harboring “premature” termination codons (PTC) and nascent polypeptide chains synthesized from such transcripts for efficient degradation (Shoemaker and Green, 2012). Components of NMD pathway distinguishes PTCs from authentic stop codons in the coding sequence. PTCs are usually the product of point-non-sense mutations, ribosomal frameshifting on slippery sequences, aberrant splicing events, or in some cases, the consequence of targeted gene regulation through alternative splicing (Sorber et al., 2011; Yeoh et al., 2019). In higher eukaryotes, PTCs are generally recognized by their proximity to so-called exon-junction complexes (EJCs), which are deposited near exon junctions during pre-mRNA splicing in the nucleus (Shoemaker and Green, 2012).
No-Go decay is an “umbrella term” for the mRNA surveillance pathway that deals with either damaged or difficult to translate mRNA sequences that cause the ribosome to stall during the elongation cycle of translation. Besides the previously mentioned mRNA base damages (8-oxoG) (Simms et al., 2014; Simms and Zaher, 2016), mRNA translation can be stalled by lack of aminoacylated-tRNAs, strong mRNA secondary structure (i.e., stem-loops or long GC-rich regions), or stable interaction of the nascent polypeptide chain with the translating ribosome. Even though Non-Stop Decay (NSD) was discovered earlier than NGD (Doma and Parker, 2006; Izawa et al., 2012; Tsuboi et al., 2012; Saito et al., 2013; Martin et al., 2014; Guydosh and Green, 2017), it became apparent that in mammals and higher eukaryotes, the NSD and NGD pathways share the same effector protein complexes (Saito et al., 2013). The NSD targeted mRNAs that originate from premature 3′ adenylation or cryptic polyadenylation signals found in coding sequences indeed represent a similar group of transcripts that would be targeted by NGD pathway (Saito et al., 2013; Kashima et al., 2014; Martin et al., 2014). Ribosomes that translate mRNAs without stop codons would eventually stall while translating long polyA tails into poly-lysine repeats, or because they would simply run out of message. Recognition of these types of transcripts, as well as the aforementioned NGD targets, trigger components of NGD/NSD pathways resulting in targeted mRNA cleavage and degradation.
The majority of mRNA surveillance pathway genes have been annotated in the P. falciparum genome (Table 3; Hughes et al., 2010). However, there are no mechanistic studies to confirm the activity of these pathways. Most of our knowledge on Plasmodium’s mRNA surveillance pathways comes from bioinformatic searches using homologous sequences from other eukaryotes. An indirect proof of the existence of NMD in Plasmodium is through the studies of alternative splicing of pre-mRNA (Sorber et al., 2011; Yeoh et al., 2019). Regulated alternative splicing events generating transcripts that do not lead to apparent protein synthesis usually carry PTCs, and thus are committed to NMD. Alternative splicing in P. falciparum has been reported for several genes like delta-aminolevulinic acid dehydratase (ALAD), stromal processing peptidase (SPP), and chloroquine resistance transporter (PfCRT); among the others. Additionally, studies on the essentiality of Plasmodium genes that use the CRISPR/Cas9 technique (Ghorbal et al., 2014) or transposon techniques (Zhang et al., 2018) rely partially on silencing targeted genes through the activation of NMD. In this case, activation of NMD is the consequence of either mutations that are generated during CRISPR/Cas9 DNA cleavage, transposon insertion in the coding sequence, or due to aberrant splicing events caused by transposons landing in introns of targeted genes. As noted above, more than 60% of the P. falciparum transcripts harbor polyA track motifs that are seen as mRNA “slippery” sequences during translation (Habich, 2016; Djuranovic et al., 2018). Translation of runs of poly-adenosine nucleotides results in ribosomal frameshifting in most tested organisms causing activation of NMD pathways (Arthur et al., 2015; Koutmou et al., 2015). However, ribosome profiling (Lacsina et al., 2011; Bunnik et al., 2013) and reporter assays (Djuranovic et al., 2018) indicate that P. falciparum ribosomes maintain fidelity during translation of rather long polyA stretches (more than 36As in a row). Therefore, while there is indirect evidence that the NMD pathway exists in Plasmodium, it seems that this pathway is not upregulated during Plasmodium ribosomes’ interactions with its polyA runs and AU-rich coding sequences. The most probable reason for this is the above mentioned changes in ribosome structure and fidelity.
TABLE 3
| Pathway | H. sapiens | A. gambiae | P. falciparum |
| NMD | eRF1, eRF3, UPF1, UPF2, UPF3A/UPF3B, eIF4AIII, MLN51, Y14/MAGOH, BTZ, SMG1, SMG5, SMG6, SMG7, PP2, Musashi, PABP1 | eRF1 (AGAP010310), eRF3 (AGAP009310), UPF1 (AGAP001133), UPF2 (AGAP000337), UPF3 (AGAP006649), eIF4AIII (AGAP003089), Y14 (AGAP006365)/MAGOH (AGAP010755), SMG1(AGAP000368), SMG5 (AGAP008181), SMG6 (AGAP000894), PP2A (AGAP004096), Musashi (AGAP001930) | UPF1 (PF3D7_1005500), UPF2 (PF3D7_0925800), UPF3B (PF3D7_1327700), eIF4AIII (PF3D7_0422700) PF3D7_1327700), PP2A(PF3D7_0925400 – KEGG, PF3D7_1319700 – Hs homology, or PF3D7_0927700 – name), Musashi (PF3D7_0916700), PABP1 (PF3D7_1224300), eRF1 (PF3D7_0212300), eRF3 (PF3D7_1123400) |
| NGD/NSD | Pelota/HBS1L, RACK1, ZNF598, N4BP2 (Cue2) | Pelota (AGAP008269), HBS1L (AGAP002603), RACK1 (AGAP010173), ZNF598 homolog (AGAP007725), N4BP2 homolog (AGAP002516) | Pelota (PF3D7_0722100), RACK1 (PF3D7_0826700), ZNF598 (PF3D7_1450400) |
| RQC | CNOT4, ABCE1, TRIP4, ASCC2, ASCC3, NEMF, Listerin, UBE2D1, XRN1 | CNOT4, ABCE1 (AGAP002182), ASCC2 homolog (AGAP000428), ASCC3 (AGAP001234), NEMF homolog (AGAP002680), Ltn1 (AGAP007143), UBE2D1 homolog (AGAP000145), XRN1 | ABCE1 (PF3D7_1368200), NEMF homolog (PF3D7_1202600), Listerin homolog (PF3D7_0615600), CNOT4 (PF3D7_1235300), ASCC3 homolog (PF3D7_1439100), UBE2D1 homolog (PF3D7_1203900) XRN1 (PF3D7_0909400) |
Comparison of translation quality control factors in P. falciparum and its mosquito and human hosts.
Factors associated with NMD, NGD/NSD, and RQC pathways from the literature and KEGG pathway database in human cells were collected (Kanehisa et al., 2019). Homologs in A. gambiae, one of the most common and effected vectors of P. falciparum (Cohuet et al., 2006; Annan et al., 2007; Giraldo-Calderon et al., 2015), were collected using KEGG pathways and performing protein-BLAST searching using VectorBase. Confirmation to FlyBase was also used to confirm vague annotations. P. falciparum factors were similarly documented again using KEGG pathways and PlasmoDB protein-BLAST analysis (Aurrecoechea et al., 2009; Kanehisa et al., 2019; Thurmond et al., 2019). Homologous gene database IDs are listed for reference. Notably, NGD/NSD factors Hbs1L and Cue2 endonuclease are missing in P. falciparum genome.
Genomic sequencing has also revealed several critical components of surveillance pathways that are missing. According to NCBI, KEGG, and plasmoDB databases, P. falciparum and the majority of other Plasmodium spp. lack the NGD and NSD decay pathways components Hbs1 (Doma and Parker, 2006) and Cue2-RNA endonuclease (D’Orazio et al., 2019). With the exception of S. cerevisiae, the Hbs1/Pelo protein complex rescues stalled ribosomes on mRNAs. It was postulated that stalling events cause ribosome collisions (Simms et al., 2017), generating unique disome units consisting of the stalled ribosome and the following colliding ribosome (Beckman and Inada). The disome, as a minimal ribosome collision unit, is recognized by Ribosome-associated Quality Control (RQC) and NGD pathways (Ito-Harashima et al., 2007; Izawa et al., 2012; Tsuboi et al., 2012; Guydosh and Green, 2017; Juszkiewicz and Hegde, 2017). Activation of RQC and NGD leads to cleavage of stalled mRNA by Cue2, and possibly other unknown endonucleases, which ultimately leads to ribosome rescue by the activity of the Pelo/Hbs1 complex (Ito-Harashima et al., 2007; Tsuboi et al., 2012; Kashima et al., 2014; Matsuda et al., 2014; Sugiyama et al., 2019). In most of the above mentioned RQC and NGD studies, a typical substrate for ribosomal stalling is a long polyA run, ranging from 36 to 60 adenosines, coding for a peptide with 12–20 consecutive lysine residues. However, long polyA stretches in P. falciparum cells are efficiently translated into poly-lysine repeats (Lacsina et al., 2011; Bunnik et al., 2013). Of note, the longest endogenous polyA runs in different P. falciparum species range from 88 to 111 nucleotides and code for Plasmodium specific and hypothetical proteins (Habich, 2016), which is longer than the length of the normal 3′ polyA tail in either S. cerevisiae or human cells (Brown and Sachs, 1998; Chang et al., 2014; Subtelny et al., 2014). As such, many endogenous Plasmodium transcripts would be NSD targets in other eukaryotic organisms. It is also a question as to what the signal for NSD pathway is in Plasmodium as recent study on 3′ mRNA polyadenylation in apicomplexans did not find any differences in P. falciparum polyadenylation complex, polyA binding proteins, or polyA tails when compared to other species (Stevens et al., 2018; Kanehisa et al., 2019). Because Plasmodium lacks the components to rescue stalled ribosomes, and because Plasmodium ribosomes efficiently translate long polyA runs, the function and mechanism of the NGD/NSD pathway in P. falciparum remains a mystery.
Conclusion
While it may seem reasonable that P. falciparum adapted its ribosomes for higher fidelity on polyA runs and in parallel lost the ability to activate the RQC/NGD/NSD pathways, such a scenario is far from obvious. The absence of mRNA surveillance pathway components or deletion of RQC factors leads to both protein aggregation and proteotoxic stress in yeast cells (Choe et al., 2016; Yonashiro et al., 2016; Jamar et al., 2018). Protein aggregation is observed in P. falciparum in the absence of heat shock protein 110 (Muralidharan et al., 2012) but not due to the absence of mRNA surveillance or RQC pathways or as a consequence of increase in both number or length of polyA tracks (Djuranovic et al., 2018). This conflicting result, along with the surprising lack of interaction between the ribosomal scaffold protein RACK1/Asc1 and Plasmodium ribosomes (Wong et al., 2014; Sun et al., 2015), argue that the mRNA surveillance pathways in P. falciparum are inherently different from those in other eukaryotes. The diversity of rRNAs, Plasmodium’s ribosome structure, and the activity of yet unknown ribosome associated factors promote the possibility of “specialized ribosomes” in Plasmodium that allow for polyA tracks translation into functional proteins. Each of the aforementioned changes in parasites translational machinery and mRNA quality control pathways come at the cost of self-fitness that would normally be detrimental for survival of Plasmodium parasites in both humans and mosquitos. And yet the parasite has persisted in both of these hosts for hundreds of millions of years. Parasitologists and epidemiologists have wondered “How?” for decades; now as translational biologists, we add our voices to the same question.
Statements
Author contributions
All authors contributed equally to writing of this review.
Funding
This work in Djuranovic’s lab was supported by NIH R01 GM112824 to SD and LEAP Awards to SD and SPD. JE was supported by NIH T32 GM007067 and NIMH R01 MH116999.
Acknowledgments
We are thankful to Caitlin D. Hanlon for critical reading and suggestions on review topic.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AlbàM. M.TompaP.VeitiaR. A. (2007). Amino acid repeats and the structure and evolution of proteins.Genome Dyn.3119–130. 10.1159/000107607
2
AnnanZ.DurandP.AyalaF. J.ArnathauC.Awono-AmbeneP.SimardF.et al (2007). Population genetic structure of Plasmodium falciparum in the two main African vectors, Anopheles gambiae and Anopheles funestus.Proc. Natl. Acad. Sci. U.S.A.1047987–7992. 10.1073/pnas.0702715104
3
ArthurL.Pavlovic-DjuranovicS.Smith-KoutmouK.GreenR.SzczesnyP.DjuranovicS. (2015). Translational control by lysine-encoding A-rich sequences.Sci. Adv.1:e1500154. 10.1126/sciadv.1500154
4
AurrecoecheaC.BrestelliJ.BrunkB. P.DommerJ.FischerS.GajriaB.et al (2009). PlasmoDB: a functional genomic database for malaria parasites.Nucleic Acids Res.37D539–D543. 10.1093/nar/gkn814
5
BabbittS. E.AltenhofenL.CobboldS. A.IstvanE. S.FennellC.DoerigC.et al (2012). Plasmodium falciparum responds to amino acid starvation by entering into a hibernatory state.Proc. Natl. Acad. Sci. U.S.A.109E3278–E3287. 10.1073/pnas.1209823109
6
BeckerK.TilleyL.VennerstromJ. L.RobertsD.RogersonS.GinsburgH. (2004). Oxidative stress in malaria parasite-infected erythrocytes: host-parasite interactions.Int. J. Parasitol.34163–189. 10.1016/j.ijpara.2003.09.011
7
BowmanS.LawsonD.BashamD.BrownD.ChillingworthT.ChurcherC. M.et al (1999). The complete nucleotide sequence of chromosome 3 of Plasmodium falciparum.Nature400532–538.
8
BrandmanO.Stewart-OrnsteinJ.WongD.LarsonA.WilliamsC. C.LiG. W.et al (2012). A ribosome-bound quality control complex triggers degradation of nascent peptides and signals translation stress.Cell1511042–1054. 10.1016/j.cell.2012.10.044
9
BrownC. E.SachsA. B. (1998). Poly(A) tail length control in Saccharomyces cerevisiae occurs by message-specific deadenylation.Mol. Cell. Biol.186548–6559. 10.1128/mcb.18.11.6548
10
BunnikE. M.Le RochK. G. (2013). An introduction to functional genomics and systems biology.Adv. Wound Care2490–498.
11
BunnikE. M.ChungD. W.HamiltonM.PontsN.SarafA.PrudhommeJ.et al (2013). Polysome profiling reveals translational control of gene expression in the human malaria parasite Plasmodium falciparum.Genome Biol.14:R128. 10.1186/gb-2013-14-11-r128
12
CaroF.AhyongV.BetegonM.DeRisiJ. L. (2014). Genome-wide regulatory dynamics of translation in the Plasmodium falciparum asexual blood stages.eLife3:e04106. 10.7554/eLife.04106
13
CasseraM. B.ZhangY.HazletonK. Z.SchrammV. L. (2011). Purine and pyrimidine pathways as targets in Plasmodium falciparum.Curr. Top. Med. Chem.112103–2115. 10.2174/156802611796575948
14
ChangH.LimJ.HaM.KimV. N. (2014). TAIL-seq: genome-wide determination of poly(A) tail length and 3′ end modifications.Mol. Cell.531044–1052. 10.1016/j.molcel.2014.02.007
15
ChoeY. J.ParkS. H.HassemerT.KornerR.Vincenz-DonnellyL.Hayer-HartlM.et al (2016). Failure of RQC machinery causes protein aggregation and proteotoxic stress.Nature531191–195. 10.1038/nature16973
16
CohuetA.OstaM. A.MorlaisI.Awono-AmbeneP. H.MichelK.SimardF.et al (2006). Anopheles and Plasmodium: from laboratory models to natural systems in the field.EMBO Rep.71285–1289.
17
CornillotE.Hadj-KaddourK.DassouliA.NoelB.RanwezV.VacherieB.et al (2012). Sequencing of the smallest Apicomplexan genome from the human pathogen Babesia microti.Nucleic Acids Res.409102–9114. 10.1093/nar/gks700
18
CuiL.LindnerS.MiaoJ. (2015). Translational regulation during stage transitions in malaria parasites.Ann. N. Y. Acad. Sci.13421–9. 10.1111/nyas.12573
19
DietelA. K.MerkerH.KaltenpothM.KostC. (2019). Selective advantages favour high genomic AT-contents in intracellular elements.PLoS Genet.15:e1007778. 10.1371/journal.pgen.1007778
20
DimitrovaL. N.KurohaK.TatematsuT.InadaT. (2009). Nascent peptide-dependent translation arrest leads to Not4p-mediated protein degradation by the proteasome.J. Biol. Chem.28410343–10352. 10.1074/jbc.M808840200
21
DjuranovicS. P.ErathJ.AndrewsR. J.BayguinovP. O.ChungJ. J.ChalkerD. L.et al (2018). PolyA tracks and poly-lysine repeats are the achilles heel of Plasmodium falciparum.bioRxiv [Preprint]. 10.1101/420109
22
DomaM. K.ParkerR. (2006). Endonucleolytic cleavage of eukaryotic mRNAs with stalls in translation elongation.Nature440561–564. 10.1038/nature04530
23
D’OrazioK. N.WuC. C.SinhaN.Loll-KrippleberR.BrownG. W.GreenR. (2019). The endonuclease Cue2 cleaves mRNAs at stalled ribosomes during No Go Decay.eLife8:e49117. 10.7554/eLife.49117
24
El BissatiK.ZuffereyR.WitolaW. H.CarterN. S.UllmanB.Ben MamounC. (2006). The plasma membrane permease PfNT1 is essential for purine salvage in the human malaria parasite Plasmodium falciparum.Proc. Natl. Acad. Sci. U.S.A.1039286–9291. 10.1073/pnas.0602590103
25
EllwoodM.NomuraM. (1982). Chromosomal locations of the genes for rRNA in Escherichia coli K-12.J. Bacteriol.149458–468.
26
FangJ.McCutchanT. F. (2002). Thermoregulation in a parasite’s life cycle.Nature418:742. 10.1038/418742a
27
FlorensL.WashburnM. P.RaineJ. D.AnthonyR. M.GraingerM.HaynesJ. D.et al (2002). A proteomic view of the Plasmodium falciparum life cycle.Nature419520–526. 10.1038/nature01107
28
GardnerM. J.HallN.FungE.WhiteO.BerrimanM.HymanR. W.et al (2002a). Genome sequence of the human malaria parasite Plasmodium falciparum.Nature419498–511.
29
GardnerM. J.ShallomS. J.CarltonJ. M.SalzbergS. L.NeneV.ShoaibiA.et al (2002b). Sequence of Plasmodium falciparum chromosomes 2, 10, 11 and 14.Nature419531–534.
30
GardnerM. J.TettelinH.CarucciD. J.CummingsL. M.AravindL.KooninE. V.et al (1998). Chromosome 2 sequence of the human malaria parasite Plasmodium falciparum.Science2821126–1132. 10.1126/science.282.5391.1126
31
GarziaA.JafarnejadS. M.MeyerC.ChapatC.GogakosT.MorozovP.et al (2017). The E3 ubiquitin ligase and RNA-binding protein ZNF598 orchestrates ribosome quality control of premature polyadenylated mRNAs.Nat. Commun.8:16056. 10.1038/ncomms16056
32
GhorbalM.GormanM.MacphersonC. R.MartinsR. M.ScherfA.Lopez-RubioJ. J. (2014). Genome editing in the human malaria parasite Plasmodium falciparum using the CRISPR-Cas9 system.Nat. Biotechnol.32819–821. 10.1038/nbt.2925
33
Giraldo-CalderonG.IEmrichS. J.MacCallumR. M.MaslenG.DialynasE.TopalisP.et al (2015). VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases.Nucleic Acids Res.43D707–D713. 10.1093/nar/gku1117
34
GlöcknerG. (2000). Large scale sequencing and analysis of AT rich eukaryote genomes.Curr. Genom.1289–299. 10.2174/1389202003351472
35
GulerJ. L.FreemanD. L.AhyongV.PatrapuvichR.WhiteJ.GujjarR.et al (2013). Asexual populations of the human malaria parasite, Plasmodium falciparum, use a two-step genomic strategy to acquire accurate, beneficial DNA amplifications.PLoS Pathog.9:e1003375. 10.1371/journal.ppat.1003375
36
GundersonJ. H.SoginM. L.WollettG.HollingdaleM.de la CruzV. F.WatersA. P.et al (1987). Structurally distinct, stage-specific ribosomes occur in Plasmodium.Science238933–937. 10.1126/science.3672135
37
GuydoshN. R.GreenR. (2017). Translation of poly(A) tails leads to precise mRNA cleavage.RNA23749–761. 10.1261/rna.060418.116
38
HabichM. (2016). PATACSDB - The database of polyA translational attenuators in coding sequences.PeerJ Comput. Sci.2:e45. 10.7717/peerj-cs.45
39
HallN.PainA.BerrimanM.ChurcherC.HarrisB.HarrisD.et al (2002). Sequence of Plasmodium falciparum chromosomes 1, 3-9 and 13.Nature419527–531.
40
HamiltonW. L.ClaessensA.OttoT. D.KekreM.FairhurstR. M.RaynerJ. C.et al (2017). Extreme mutation bias and high AT content in Plasmodium falciparum.Nucleic Acids Res.451889–1901. 10.1093/nar/gkw1259
41
HuckabyA. C.GranumC. S.CareyM. A.SzlachtaK.Al-BarghouthiB.WangY. H.et al (2019). Complex DNA structures trigger copy number variation across the Plasmodium falciparum genome.Nucleic Acids Res.471615–1627. 10.1093/nar/gky1268
42
HughesK. R.PhilipN.StarnesG. L.TaylorS.WatersA. P. (2010). From cradle to grave: RNA biology in malaria parasites.Wiley Interdiscip. Rev. RNA1287–303. 10.1002/wrna.30
43
HymanR. W.FungE.ConwayA.KurdiO.MaoJ.MirandaM.et al (2002). Sequence of Plasmodium falciparum chromosome 12.Nature419534–537.
44
Ito-HarashimaS.KurohaK.TatematsuT.InadaT. (2007). Translation of the poly(A) tail plays crucial roles in nonstop mRNA surveillance via translation repression and protein destabilization by proteasome in yeast.Genes Dev.21519–524. 10.1101/gad.1490207
45
IzawaT.TsuboiT.KurohaK.InadaT.NishikawaS.EndoT. (2012). Roles of dom34:hbs1 in nonstop protein clearance from translocators for normal organelle protein influx.Cell Rep.2447–453. 10.1016/j.celrep.2012.08.010
46
JamarN. H.KritsiligkouP.GrantC. M. (2018). Loss of mRNA surveillance pathways results in widespread protein aggregation.Sci. Rep.8:3894. 10.1038/s41598-018-22183-2
47
JuszkiewiczS.HegdeR. S. (2017). Initiation of quality control during poly(A) translation requires site-specific ribosome Ubiquitination.Mol. Cell65743–750.e4. 10.1016/j.molcel.2016.11.039
48
KanehisaM.SatoY.FurumichiM.MorishimaK.TanabeM. (2019). New approach for understanding genome variations in KEGG.Nucleic Acids Res.47D590–D595. 10.1093/nar/gky962
49
KashimaI.TakahashiM.HashimotoY.SakotaE.NakamuraY.InadaT. (2014). A functional involvement of ABCE1, eukaryotic ribosome recycling factor, in nonstop mRNA decay in drosophila Melanogaster cells.Biochimie10610–16. 10.1016/j.biochi.2014.08.001
50
KimH. K.LiuF.FeiJ.BustamanteC.GonzalezR. L.TinocoI. (2014). A frameshifting stimulatory stem loop destabilizes the hybrid state and impedes ribosomal translocation.Proc. Natl. Acad. Sci. U.S.A.1115538–5543. 10.1073/pnas.1403457111
51
KobayashiT. (2014). Ribosomal RNA gene repeats, their stability and cellular senescence.Proc. JPN Acad. Ser. B Phys. Biol. Sci.90119–129. 10.2183/pjab.90.119
52
KooijT. W.JanseC. J.WatersA. P. (2006). Plasmodium post-genomics: better the bug you know?Nat. Rev. Microbiol.4344–357. 10.1038/nrmicro1392
53
KoutmouK. S.SchullerA. P.BrunelleJ. L.RadhakrishnanA.DjuranovicS.GreenR. (2015). Ribosomes slide on lysine-encoding homopolymeric a stretches.eLife4:e05534. 10.7554/eLife.05534
54
KurohaK.AkamatsuM.DimitrovaL.ItoT.KatoY.ShirahigeK.et al (2010). Receptor for activated C kinase 1 stimulates nascent polypeptide-dependent translation arrest.EMBO Rep.11956–961. 10.1038/embor.2010.169
55
KuryloC. M.ParksM. M.JuetteM. F.ZinshteynB.AltmanR. B.ThibadoJ. K.et al (2018). Endogenous rRNA sequence variation can regulate stress response gene expression and phenotype.Cell Rep.25236–248.e6. 10.1016/j.celrep.2018.08.093
56
LacsinaJ. R.LaMonteG.NicchittaC. V.ChiJ. T. (2011). Polysome profiling of the malaria parasite Plasmodium falciparum.Mol. Biochem. Parasitol.17942–46. 10.1016/j.molbiopara.2011.05.003
57
Le RochK. G.JohnsonJ. R.FlorensL.ZhouY.SantrosyanA.GraingerM.et al (2004). Global analysis of transcript and protein levels across the Plasmodium falciparum life cycle.Genom. Res.142308–2318. 10.1101/gr.2523904
58
LeiriaoP.RodriguesC. D.AlbuquerqueS. S.MotaM. M. (2004). Survival of protozoan intracellular parasites in host cells.EMBO Rep.51142–1147. 10.1038/sj.embor.7400299
59
LiJ.GutellR. R.DambergerS. H.WirtzR. A.KissingerJ. C.RogersM. J.et al (1997). Regulation and trafficking of three distinct 18 S ribosomal RNAs during development of the malaria parasite.J. Mol. Biol.269203–213. 10.1006/jmbi.1997.1038
60
LiuJ.IstvanE. S.GluzmanI. Y.GrossJ.GoldbergD. E. (2006). Plasmodium falciparum ensures its amino acid supply with multiple acquisition pathways and redundant proteolytic enzyme systems.Proc. Natl. Acad. Sci. U.S.A.1038840–8845. 10.1073/pnas.0601876103
61
LocatiM. D.PaganoJ. F. B.GirardG.EnsinkW. A.van OlstM.van LeeuwenS.et al (2017). Expression of distinct maternal and somatic 5.8S, 18S, and 28S rRNA types during zebrafish development.RNA231188–1199. 10.1261/rna.061515.117
62
López-BarragánM. J.LemieuxJ.QuiñonesM.WilliamsonK. C.Molina-CruzA.CuiK.et al (2011). Directional gene expression and antisense transcripts in sexual and asexual stages of Plasmodium falciparum.BMC Genom.12:587. 10.1186/1471-2164-12-587
63
LoyD. E.PlenderleithL. J.SundararamanS. A.LiuW.GruszczykJ.ChenY. J.et al (2018). Evolutionary history of human.Proc. Natl. Acad. Sci. U.S.A.115E8450–E8459. 10.1073/pnas.1810053115
64
LuJ.DeutschC. (2008). Electrostatics in the ribosomal tunnel modulate chain elongation rates.J. Mol. Biol.38473–86. 10.1016/j.jmb.2008.08.089
65
LuX. M.BatugedaraG.LeeM.PrudhommeJ.BunnikE. M.Le RochK. G. (2017). Nascent RNA sequencing reveals mechanisms of gene regulation in the human malaria parasite Plasmodium falciparum.Nucleic Acids Res.457825–7840. 10.1093/nar/gkx464
66
MajzoubK.HafirassouM. L.MeigninC.GotoA.MarziS.FedorovaA.et al (2014). RACK1 controls IRES-mediated translation of viruses.Cell1591086–1095. 10.1016/j.cell.2014.10.041
67
MartinL.GrigoryanA.WangD.WangJ.BredaL.RivellaS.et al (2014). Identification and characterization of small molecules that inhibit nonsense-mediated RNA decay and suppress nonsense p53 mutations.Cancer Res.743104–3113. 10.1158/0008-5472.CAN-13-2235
68
MatsudaR.IkeuchiK.NomuraS.InadaT. (2014). Protein quality control systems associated with no-go and nonstop mRNA surveillance in yeast.Genes Cells191–12. 10.1111/gtc.12106
69
McCutchanT. F.de la CruzV. F.LalA. A.GundersonJ. H.ElwoodH. J.SoginM. L. (1988). Primary sequences of two small subunit ribosomal RNA genes from Plasmodium falciparum.Mol. Biochem. Parasitol.2863–68. 10.1016/0166-6851(88)90181-8
70
MillsE. W.GreenR. (2017). Ribosomopathies: there’s strength in numbers.Science358:eaan2755. 10.1126/science.aan2755
71
MuralidharanV.GoldbergD. E. (2013). Asparagine repeats in Plasmodium falciparum proteins: good for nothing?PLoS Pathog.9:e1003488. 10.1371/journal.ppat.1003488
72
MuralidharanV.OksmanA.IwamotoM.WandlessT. J.GoldbergD. E. (2011). Asparagine repeat function in a Plasmodium falciparum protein assessed via a regulatable fluorescent affinity tag.Proc. Natl. Acad. Sci. U.S.A.1084411–4416. 10.1073/pnas.1018449108
73
MuralidharanV.OksmanA.PalP.LindquistS.GoldbergD. E. (2012). Plasmodium falciparum heat shock protein 110 stabilizes the asparagine repeat-rich parasite proteome during malarial fevers.Nat. Commun.3:1310. 10.1038/ncomms2306
74
MustoH.RomeroH.ZavalaA.JabbariK.BernardiG. (1999). Synonymous codon choices in the extremely GC-poor genome of Plasmodium falciparum: compositional constraints and translational selection.J. Mol. Evol.4927–35. 10.1007/pl00006531
75
OsteraG.TokumasuF.TeixeiraC.CollinN.SaJ.HumeJ.et al (2011). Plasmodium falciparum: nitric oxide modulates heme speciation in isolated food vacuoles.Exp. Parasitol.1271–8. 10.1016/j.exppara.2010.05.006
76
OttoT. D.GilabertA.CrellenT.BöhmeU.ArnathauC.SandersM.et al (2018). Genomes of all known members of a Plasmodium subgenus reveal paths to virulent human malaria.Nat. Microbiol.3687–697. 10.1038/s41564-018-0162-2
77
PollackY.KatzenA. L.SpiraD. T.GolenserJ. (1982). The genome of Plasmodium falciparum. I: DNA base composition.Nucleic Acids Res.10539–546. 10.1093/nar/10.2.539
78
QiY.ZhuF.EastmanR. T.FuY.ZilversmitM.PattaradilokratS.et al (2015). Regulation of Plasmodium yoelii oocyst development by strain- and stage-specific small-subunit rRNA.mBio6:e00117. 10.1128/mBio.00117-15
79
QuashieN. B.Dorin-SemblatD.BrayP. G.BiaginiG. A.DoerigC.Ranford-CartwrightL. C.et al (2008). A comprehensive model of purine uptake by the malaria parasite Plasmodium falciparum: identification of four purine transport activities in intraerythrocytic parasites.Biochem. J.411287–295. 10.1042/BJ20071460
80
ReddyB. P.ShresthaS.HartK. J.LiangX.KemirembeK.CuiL.et al (2015). A bioinformatic survey of RNA-binding proteins in Plasmodium.BMC Genom.16:890. 10.1186/s12864-015-2092-1
81
RogersM. J.GutellR. R.DambergerS. H.LiJ.McConkeyG. A.WatersA. P.et al (1996). Structural features of the large subunit rRNA expressed in Plasmodium falciparum sporozoites that distinguish it from the asexually expressed subunit rRNA.RNA2134–145.
82
SaitoS.HosodaN.HoshinoS. (2013). The Hbs1-Dom34 protein complex functions in non-stop mRNA decay in mammalian cells.J. Biol. Chem.28817832–17843. 10.1074/jbc.M112.448977
83
SakaiK.OhtaT.MinoshimaS.KudohJ.WangY.de JongP. J.et al (1995). Human ribosomal RNA gene cluster: identification of the proximal end containing a novel tandem repeat sequence.Genomics26521–526. 10.1016/0888-7543(95)80170-q
84
SenguptaJ.NilssonJ.GurskyR.SpahnC. M.NissenP.FrankJ. (2004). Identification of the versatile scaffold protein RACK1 on the eukaryotic ribosome by cryo-EM.Nat. Struct. Mol. Biol.11957–962. 10.1038/nsmb822
85
SewardE. A.KellyS. (2016). Dietary nitrogen alters codon bias and genome composition in parasitic microorganisms.Genom. Biol.17:226.
86
ShermanI. W. (1979). Biochemistry of Plasmodium (malarial parasites).Microbiol. Rev.43453–495.
87
ShockJ. L.FischerK. F.DeRisiJ. L. (2007). Whole-genome analysis of mRNA decay in Plasmodium falciparum reveals a global lengthening of mRNA half-life during the intra-erythrocytic development cycle.Genom. Biol.8:R134.
88
ShoemakerC. J.GreenR. (2012). Translation drives mRNA quality control.Nat. Struct. Mol. Biol.19594–601. 10.1038/nsmb.2301
89
SilvestriniF.LasonderE.OlivieriA.CamardaG.van SchaijkB.SanchezM.et al (2010). Protein export marks the early phase of gametocytogenesis of the human malaria parasite Plasmodium falciparum.Mol. Cell. Proteomics91437–1448. 10.1074/mcp.M900479-MCP200
90
SimmsC. L.ZaherH. S. (2016). Quality control of chemically damaged RNA.Cell Mol. Life Sci.733639–3653. 10.1007/s00018-016-2261-7
91
SimmsC. L.HudsonB. H.MosiorJ. W.RangwalaA. S.ZaherH. S. (2014). An active role for the ribosome in determining the fate of oxidized mRNA.Cell Rep.91256–1264. 10.1016/j.celrep.2014.10.042
92
SimmsC. L.YanL. L.ZaherH. S. (2017). Ribosome collision is critical for quality control during No-Go Decay.Mol. Cell68361–373.e5. 10.1016/j.molcel.2017.08.019
93
SorberK.DimonM. T.DeRisiJ. L. (2011). RNA-Seq analysis of splicing in Plasmodium falciparum uncovers new splice junctions, alternative splicing and splicing of antisense transcripts.Nucleic Acids Res.393820–3835. 10.1093/nar/gkq1223
94
StevensA. T.HoweD. K.HuntA. G. (2018). Characterization of mRNA polyadenylation in the apicomplexa.PLoS One13:e0203317. 10.1371/journal.pone.0203317
95
SubtelnyA. O.EichhornS. W.ChenG. R.SiveH.BartelD. P. (2014). Poly(A)-tail profiling reveals an embryonic switch in translational control.Nature50866–71. 10.1038/nature13007
96
SueokaN. (1962). On the genetic basis of variation and heterogeneity of DNA base composition.Proc. Natl. Acad. Sci. U.S.A.48582–592. 10.1073/pnas.48.4.582
97
SugiyamaT.LiS.KatoM.IkeuchiK.IchimuraA.MatsuoY.et al (2019). Sequential Ubiquitination of ribosomal protein uS3 triggers the degradation of non-functional 18S rRNA.Cell Rep.263400–3415.e7. 10.1016/j.celrep.2019.02.067
98
SunM.LiW.BlomqvistK.DasS.HashemY.DvorinJ. D.et al (2015). Dynamical features of the Plasmodium falciparum ribosome during translation.Nucleic Acids Res.4310515–10524. 10.1093/nar/gkv991
99
SundaramoorthyE.LeonardM.MakR.LiaoJ.FulzeleA.BennettE. J. (2017). ZNF598 and RACK1 regulate mammalian ribosome-associated quality control function by mediating regulatory 40S ribosomal Ubiquitylation.Mol. Cell.65751–760.e4. 10.1016/j.molcel.2016.12.026
100
SzafranskiK.LehmannR.ParraG.GuigoR.GlöcknerG. (2005). Gene organization features in A/T-rich organisms.J. Mol. Evol.6090–98. 10.1007/s00239-004-0201-2
101
TangT. T. L.StowellJ. A. W.HillC. H.PassmoreL. A. (2019). The intrinsic structure of poly(A) RNA determines the specificity of Pan2 and Caf1 deadenylases.Nat. Struct. Mol. Biol.26433–442. 10.1038/s41594-019-0227-9
102
TesinaP.LessenL. N.BuschauerR.ChengJ.WuC.BerninghausenO.et al (2019). Molecular mechanism of translational stalling by inhibitory codon combinations and poly(A) tracts.bioRxiv [Preprint]. 10.1101/755652
103
ThurmondJ.GoodmanJ. L.StreletsV. B.AttrillH.GramatesL. S.MarygoldS. J.et al (2019). FlyBase 2.0: the next generation.Nucleic Acids Res.47D759–D765. 10.1093/nar/gky1003
104
TingL. M.ShiW.LewandowiczA.SinghV.MwakingweA.BirckM. R.et al (2005). Targeting a novel Plasmodium falciparum purine recycling pathway with specific immucillins.J. Biol. Chem.2809547–9554. 10.1074/jbc.m412693200
105
TournuH.ButtsA.PalmerG. E. (2019). Titrating gene function in the human fungal pathogen Candida albicans through Poly-adenosine tract insertion.mSphere4e192–e119. 10.1128/mSphere.00192-19
106
TsuboiT.KurohaK.KudoK.MakinoS.InoueE.KashimaI.et al (2012). Dom34:hbs1 plays a general role in quality-control systems by dissociation of a stalled ribosome at the 3′ end of aberrant mRNA.Mol. Cell.46518–529. 10.1016/j.molcel.2012.03.013
107
van SpaendonkR. M.RamesarJ.van WigcherenA.ElingW.BeetsmaA. L.van GemertJ.et al (2001). Functional equivalence of structurally distinct ribosomes in the malaria parasite, Plasmodium berghei.J. Biol. Chem.27622638–22647. 10.1074/jbc.m101234200
108
VelichutinaI. V.RogersM. J.McCutchanT. F.LiebmanS. W. (1998). Chimeric rRNAs containing the GTPase centers of the developmentally regulated ribosomal rRNAs of Plasmodium falciparum are functionally distinct.RNA4594–602. 10.1017/s1355838298980049
109
VembarS. S.DrollD.ScherfA. (2016a). Translational regulation in blood stages of the malaria parasite Plasmodium spp.: systems-wide studies pave the way.Wiley Interdiscip. Rev. RNA7772–792. 10.1002/wrna.1365
110
VembarS. S.SeetinM.LambertC.NattestadM.SchatzM. C.BaybayanP.et al (2016b). Complete telomere-to-telomere de novo assembly of the Plasmodium falciparum genome through long-read (>11 kb), single molecule, real-time sequencing.DNA Res.23339–351. 10.1093/dnares/dsw022
111
VidevallE. (2018). Plasmodium parasites of birds have the most AT-rich genes of eukaryotes.Microb. Genom.4:e000150. 10.1099/mgen.0.000150
112
WallikerD.QuakyiI. A.WellemsT. E.McCutchanT. F.SzarfmanA.LondonW. T.et al (1987). Genetic analysis of the human malaria parasite Plasmodium falciparum.Science2361661–1666. 10.1126/science.3299700
113
WarnerJ. R. (1999). The economics of ribosome biosynthesis in yeast.Trends Biochem. Sci.24437–440. 10.1016/s0968-0004(99)01460-7
114
WatersA. P.SyinC.McCutchanT. F. (1989). Developmental regulation of stage-specific ribosome populations in Plasmodium.Nature342438–440. 10.1038/342438a0
115
WernegreenJ. J.FunkD. J. (2004). Mutation exposed: a neutral explanation for extreme base composition of an endosymbiont genome.J. Mol. Evol.59849–858. 10.1007/s00239-003-0192-z
116
WinegardT. C. (2019). The Mosquito: A Human History of Our Deadliest Predator.Melbourne: Text Publishing Company.
117
WolfA. S.GrayhackE. J. (2015). Asc1, homolog of human RACK1, prevents frameshifting in yeast by ribosomes stalled at CGA codon repeats.RNA21935–945. 10.1261/rna.049080.114
118
WongW.BaiX. C.BrownA.FernandezI. S.HanssenE.CondronM.et al (2014). Cryo-EM structure of the Plasmodium falciparum 80S ribosome bound to the anti-protozoan drug emetine.eLife3:e03080. 10.7554/eLife.03080
119
World Health Organization [WHO], (2018). World Malaria Report 2018.Geneva: World Health Organization.
120
WuH.ZhangZ.HuS.YuJ. (2012). On the molecular mechanism of GC content variation among eubacterial genomes.Biol. Direct7:2. 10.1186/1745-6150-7-2
121
XueS.BarnaM. (2012). Specialized ribosomes: a new frontier in gene regulation and organismal biology.Nat. Rev. Mol. Cell Biol.13355–369. 10.1038/nrm3359
122
YeohL. M.GoodmanC. D.MollardV.McHughE.LeeV. V.SturmA.et al (2019). Alternative splicing is required for stage differentiation in malaria parasites.Genom. Biol.20:151. 10.1186/s13059-019-1756-6
123
YonashiroR.TaharaE. B.BengtsonM. H.KhokhrinaM.LorenzH.ChenK. C.et al (2016). The Rqc2/Tae2 subunit of the ribosome-associated quality control (RQC) complex marks ribosome-stalled nascent polypeptide chains for aggregation.eLife5:e11794. 10.7554/eLife.11794
124
ZanghìG.VembarS. S.BaumgartenS.DingS.GuizettiJ.BryantJ. M.et al (2018). A specific PfEMP1 is expressed in falciparum, P., sporozoites and plays a role in hepatocyte infection.Cell Rep.222951–2963. 10.1016/j.celrep.2018.02.075
125
ZhangM.WangC.OttoT. D.OberstallerJ.LiaoX.AdapaS. R.et al (2018). Uncovering the essential genes of the human malaria parasite.Science360:eaa7847. 10.1126/science.aap7847
126
ZhuJ. D.WatersA. P.AppiahA.McCutchanT. F.LalA. A.HollingdaleM. R. (1990). Stage-specific ribosomal RNA expression switches during sporozoite invasion of hepatocytes.J. Biol. Chem.26512740–12744.
127
ZilversmitM. M.VolkmanS. K.DePristoM. A.WirthD. F.AwadallaP.HartlD. L. (2010). Low-complexity regions in Plasmodium falciparum: missing links in the evolution of an extreme genome.Mol. Biol. Evol.272198–2209. 10.1093/molbev/msq108
Summary
Keywords
plasmodium, ribosomes, mRNA surveillance, AT rich genome, mRNA translation
Citation
Erath J, Djuranovic S and Djuranovic SP (2019) Adaptation of Translational Machinery in Malaria Parasites to Accommodate Translation of Poly-Adenosine Stretches Throughout Its Life Cycle. Front. Microbiol. 10:2823. doi: 10.3389/fmicb.2019.02823
Received
26 September 2019
Accepted
21 November 2019
Published
06 December 2019
Volume
10 - 2019
Edited by
Rhoel Dinglasan, University of Florida, United States
Reviewed by
Scott E. Lindner, Pennsylvania State University (PSU), United States; Vasant Muralidharan, University of Georgia, United States
Updates
Copyright
© 2019 Erath, Djuranovic and Djuranovic.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sergej Djuranovic, sergej.djuranovic@wustl.eduSlavica Pavlovic Djuranovic, spavlov@wustl.edu
This article was submitted to Infectious Diseases, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.