Identification and analysis of the RNA degrading complexes and machinery of Giardia lamblia using an in silico approach

Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics

Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.

Multimodal regulation of encystation in Giardia duodenalis revealed by deep proteomics

Cyst formation in the parasitic protist Giardia duodenalis is critical to its transmission. Existing proteomic data quantifies only 17% of coding genes transcribed during encystation and does not cover the complete process from trophozoite to mature cyst. Using high-resolution mass spectrometry, we have quantified proteomic changes across encystation and compared this with published transcriptomic data. We reproducibly identified 3863 (64.5% of Giardia proteins) and quantified 3382 proteins (56.5% of Giardia proteins) over standard trophozoite growth (TY), during low-bile encystation priming (LB), 16 h into encystation (EC), and at cyst maturation (C). This work provides the first known expanded observation of encystation at the proteomic level and triples the coverage of previous encystation proteomes. One-third (1169 proteins) of the quantified proteome is differentially expressed in the mature cyst relative to the trophozoite, including proteasomal machinery, metabolic pathways, and secretory proteins. Changes in lipid metabolism indicated a shift in lipid species dependency during encystation. Consistent with this, we identified the first, putative lipid transporters in this species, representing the steroidogenic acute regulatory protein-related lipid transfer (StARkin), oxysterol binding protein related protein (ORP/Osh) and glycosphingolipid transfer protein (GLTP) families, and follow their differential expression over cyst formation. Lastly, we undertook correlation analyses of the transcriptome and proteome of trophozoites and cysts, and found evidence of post-transcriptional regulation of key protein classes (RNA binding proteins) and stage-specific genes (encystation markers) implicating translation-repression in encystation. We provide the most extensive proteomic analysis of encystation in Giardia to date and the first known exploration across its complete duration. This work identifies encystation as highly coordinated, involving major changes in proteostasis, metabolism and membrane dynamics, and indicates a potential role for post-transcriptional regulation, mediated through RNA-binding proteins. Together our work provides a valuable resource for Giardia research and the development of transmission-blocking anti-giardials.

Is mRNA decapping by ApaH like phosphatases present in eukaryotes beyond the Kinetoplastida?

Background

ApaH like phosphatases (ALPHs) originate from the bacterial ApaH protein and have been identified in all eukaryotic super-groups. Only two of these proteins have been functionally characterised. We have shown that the ApaH like phosphatase ALPH1 from the Kinetoplastid Trypanosoma brucei is the mRNA decapping enzyme of the parasite. In eukaryotes, Dcp2 is the major mRNA decapping enzyme and mRNA decapping by ALPHs is unprecedented, but the bacterial ApaH protein was recently found decapping non-conventional caps of bacterial mRNAs. These findings prompted us to explore whether mRNA decapping by ALPHs is restricted to Kinetoplastida or could be more widespread among eukaryotes.

Results

We screened 827 eukaryotic proteomes with a newly developed Python-based algorithm for the presence of ALPHs and used the data to characterize the phylogenetic distribution, conserved features, additional domains and predicted intracellular localisation of this protein family. For most organisms, we found ALPH proteins to be either absent (495/827 organisms) or to have non-cytoplasmic localisation predictions (73% of all ALPHs), excluding a function in mRNA decapping. Although, non-cytoplasmic ALPH proteins had in vitro mRNA decapping activity. Only 71 non-Kinetoplastida have ALPH proteins with predicted cytoplasmic localisations. However, in contrast to Kinetoplastida, these organisms also possess a homologue of Dcp2 and in contrast to ALPH1 of Kinetoplastida, these ALPH proteins are very short and consist of the catalytic domain only.

Conclusions

ALPH was present in the last common ancestor of eukaryotes, but most eukaryotes have either lost the enzyme, or use it exclusively outside the cytoplasm. The acceptance of mRNA as a substrate indicates that ALPHs, like bacterial ApaH, have a wide substrate range: the need to protect mRNAs from unregulated degradation is one possible explanation for the selection against the presence of cytoplasmic ALPH proteins in most eukaryotes. Kinetoplastida succeeded to exploit ALPH as their only or major mRNA decapping enzyme. 71 eukaryotic organisms outside the Kinetoplastid lineage have short ALPH proteins with cytoplasmic localisation predictions: whether these proteins are used as decapping enzymes in addition to Dcp2 or else have adapted to not accept mRNAs as a substrate, remains to be explored.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01858-x.

The evolution of the Puf superfamily of proteins across the tree of eukaryotes

BACKGROUND

Eukaryotic gene expression is controlled by a number of RNA-binding proteins (RBP), such as the proteins from the Puf (Pumilio and FBF) superfamily (PufSF). These proteins bind to RNA via multiple Puf repeat domains, each of which specifically recognizes a single RNA base. Recently, three diversified PufSF proteins have been described in model organisms, each of which is responsible for the maturation of ribosomal RNA or the translational regulation of mRNAs; however, less is known about the role of these proteins across eukaryotic diversity.

RESULTS

Here, we investigated the distribution and function of PufSF RBPs in the tree of eukaryotes. We determined that the following PufSF proteins are universally conserved across eukaryotes and can be broadly classified into three groups: (i) Nop9 orthologues, which participate in the nucleolar processing of immature 18S rRNA; (ii) 'classical' Pufs, which control the translation of mRNA; and (iii) PUM3 orthologues, which are involved in the maturation of 7S rRNA. In nearly all eukaryotes, the rRNA maturation proteins, Nop9 and PUM3, are retained as a single copy, while mRNA effectors ('classical' Pufs) underwent multiple lineage-specific expansions. We propose that the variation in number of 'classical' Pufs relates to the size of the transcriptome and thus the potential mRNA targets. We further distinguished full set of PufSF proteins in divergent metamonad Giardia intestinalis and initiated their cellular and biochemical characterization.

CONCLUSIONS

Our data suggest that the last eukaryotic common ancestor (LECA) already contained all three types of PufSF proteins and that 'classical' Pufs then underwent lineage-specific expansions.

Recent advances in functional research in Giardia intestinalis

This review considers current advances in tools to investigate the functional biology of Giardia, it's coding and non-coding genes, features and cellular and molecular biology. We consider major gaps in current knowledge of the parasite and discuss the present state-of-the-art in its in vivo and in vitro cultivation. Advances in in silico tools, including for the modelling non-coding RNAs and genomic elements, as well as detailed exploration of coding genes through inferred homology to model organisms, have provided significant, primary level insight. Improved methods to model the three-dimensional structure of proteins offer new insights into their function, and binding interactions with ligands, other proteins or precursor drugs, and offer substantial opportunities to prioritise proteins for further study and experimentation. These approaches can be supplemented by the growing and highly accessible arsenal of systems-based methods now being applied to Giardia, led by genomic, transcriptomic and proteomic methods, but rapidly incorporating advanced tools for detection of real-time transcription, evaluation of chromatin states and direct measurement of macromolecular complexes. Methods to directly interrogate and perturb gene function have made major leaps in recent years, with CRISPr-interference now available. These approaches, coupled with protein over-expression, fluorescent labelling and in vitro and in vivo imaging, are set to revolutionize the field and herald an exciting time during which the field may finally realise Giardia's long proposed potential as a model parasite and eukaryote.

The protein 14-3-3: A functionally versatile molecule in Giardia duodenalis

Giardia duodenalis is a cosmopolitan zoonotic protozoan parasite causing giardiasis, one of the most common diarrhoeal diseases in human and animals. Beyond its public health relevance, Giardia represents a valuable and fascinating model microorganism. The deep-branching phylogenetic position of Giardia, its simple life cycle and its minimalistic genomic and cellular organization provide a unique opportunity to define basal and "ancestral" eukaryotic functions. The eukaryotic 14-3-3 protein family represents a distinct example of phosphoserine/phosphothreonine-binding proteins. The extended network of protein-protein interactions established by 14-3-3 proteins place them at the crossroad of multiple signalling pathways that regulate physiological and pathological cellular processes. Despite the remarkable insight on 14-3-3 protein in different organisms, from yeast to humans, so far little attention was given to the study of this protein in protozoan parasites. However, in the last years, research efforts have provided evidences on unique properties of the single 14-3-3 protein of Giardia and on its association in key aspects of Giardia life cycle. In the first part of this chapter, a general overview of the features commonly shared among 14-3-3 proteins in different organisms (i.e. structure, target recognition, mode of action and regulatory mechanisms) is included. The second part focus on the current knowledge on the biochemistry and biology of the Giardia 14-3-3 protein and on the possibility to use this protein as target to propose new strategies for developing innovative antigiardial therapy.

Disruption of the RNA exosome reveals the hidden face of the malaria parasite transcriptome

Antisense transcription emerges as a key regulator of important biological processes in the human malaria parasite Plasmodium falciparum. RNA-processing factors, however, remain poorly characterized in this pathogen. Here, we purified the multiprotein RNA exosome complex of malaria parasites by affinity chromatography, using HA-tagged PfRrp4 and PfDis3 as the ligands. Seven distinct core exosome subunits (PfRrp41, PfMtr3, PfRrp42, PfRrp45, PfRrp4, PfRrp40, PfCsl4) and two exoribonuclease proteins PfRrp6 and PfDis3 are identified by mass spectrometry. Western blot analysis detects Dis3 and Rrp4 predominantly in the cytoplasmic fraction during asexual blood stage development. An inducible gene knock out of the PfDis3 subunit reveals the upregulation of structural and coding RNA, but the vast majority belongs to antisense RNA. Furthermore, we detect numerous types of cryptic unstable transcripts (CUTs) linked to virulence gene families including antisense RNA in the rif gene family. Our work highlights the limitations of steady-state RNA analysis to predict transcriptional activity and link the RNA surveillance machinery directly with post-transcriptional control and gene expression in malaria parasites.

Stress-induced nuclear depletion of Entamoeba histolytica 3'-5' exoribonuclease EhRrp6 and its role in growth and erythrophagocytosis

The 3'-5' exoribonuclease Rrp6 is a key enzyme in RNA homeostasis involved in processing and degradation of many stable RNA precursors, aberrant transcripts, and noncoding RNAs. We previously have shown that in the protozoan parasite Entamoeba histolytica, the 5'-external transcribed spacer fragment of pre-rRNA accumulates under serum starvation-induced growth stress. This fragment is a known target of degradation by Rrp6. Here, we computationally and biochemically characterized EhRrp6 and found that it contains the catalytically important EXO and HRDC domains and exhibits exoribonuclease activity with both unstructured and structured RNA substrates, which required the conserved DEDD-Y catalytic-site residues. It lacked the N-terminal PMC2NT domain for binding of the cofactor Rrp47, but could functionally complement the growth defect of a yeast rrp6 mutant. Of note, no Rrp47 homologue was detected in E. histolytica Immunolocalization studies revealed that EhRrp6 is present both in the nucleus and cytosol of normal E. histolytica cells. However, growth stress induced its complete loss from the nuclei, reversed by proteasome inhibitors. EhRrp6-depleted E. histolytica cells were severely growth restricted, and EhRrp6 overexpression protected the cells against stress, suggesting that EhRrp6 functions as a stress sensor. Importantly EhRrp6 depletion reduced erythrophagocytosis, an important virulence determinant of E. histolytica This reduction was due to a specific decrease in transcript levels of some phagocytosis-related genes (Ehcabp3 and Ehrho1), whereas expression of other genes (Ehcabp1, Ehcabp6, Ehc2pk, and Eharp2/3) was unaffected. This is the first report of the role of Rrp6 in cell growth and stress responses in a protozoan parasite.

Multispecies reconstructions uncover widespread conservation, and lineage-specific elaborations in eukaryotic mRNA metabolism

The degree of conservation and evolution of cytoplasmic mRNA metabolism pathways across the eukaryotes remains incompletely resolved. In this study, we describe a comprehensive genome and transcriptome-wide analysis of proteins involved in mRNA maturation, translation, and mRNA decay across representative organisms from the six eukaryotic super-groups. We demonstrate that eukaryotes share common pathways for mRNA metabolism that were almost certainly present in the last eukaryotic common ancestor, and show for the first time a correlation between intron density and a selective absence of some Exon Junction Complex (EJC) components in eukaryotes. In addition, we identify pathways that have diversified in individual lineages, with a specific focus on the unique gene gains and losses in members of the Excavata and SAR groups that contribute to their unique gene expression pathways compared to other organisms.

The ApaH-like phosphatase TbALPH1 is the major mRNA decapping enzyme of trypanosomes

5’-3’ decay is the major mRNA decay pathway in many eukaryotes, including trypanosomes. After deadenylation, mRNAs are decapped by the nudix hydrolase DCP2 of the decapping complex and finally degraded by the 5’-3’ exoribonuclease. Uniquely, trypanosomes lack homologues to all subunits of the decapping complex, while deadenylation and 5’-3’ degradation are conserved. Here, I show that the parasites use an ApaH-like phosphatase (ALPH1) as their major mRNA decapping enzyme. The protein was recently identified as a novel trypanosome stress granule protein and as involved in mRNA binding. A fraction of ALPH1 co-localises exclusively with the trypanosome 5’-3’ exoribonuclease XRNA to a special granule at the posterior pole of the cell, indicating a connection between the two enzymes. RNAi depletion of ALPH1 is lethal and causes a massive increase in total mRNAs that are deadenylated, but have not yet started 5’-3’ decay. These data suggest that ALPH1 acts downstream of deadenylation and upstream of mRNA degradation, consistent with a function in mRNA decapping. In vitro experiments show that recombinant, N-terminally truncated ALHP1 protein, but not a catalytically inactive mutant, sensitises the capped trypanosome spliced leader RNA to yeast Xrn1, but only if an RNA 5’ polyphosphatase is included. This indicates that the decapping mechanism of ALPH1 differs from the decapping mechanism of Dcp2 by leaving more than one phosphate group at the mRNA’s 5’ end. This is the first reported function of a eukaryotic ApaH-like phosphatase, a bacterial-derived class of enzymes present in all phylogenetic super-groups of the eukaryotic kingdom. The substrates of eukaryotic ApaH-like phosphatases are unknown. However, the substrate of the related bacterial enzyme ApaH, diadenosine tetraphosphate, is highly reminiscent of a eukaryotic mRNA cap.

Eukaryotic mRNAs are stabilised by a 5’ cap and one important step in mRNA decay is the removal of this cap by the nudix domain protein Dcp2 of the decapping complex. The decapping complex is highly conserved throughout eukaryotes, with the exception of trypanosomes that lack the entire complex. Here, I show that trypanosomes have evolved to use an ApaH-like phosphatase instead of a nudix domain protein as their major decapping enzyme. This work closes an important gap in the knowledge of trypanosome mRNA metabolism. Moreover, this is the first reported function of an ApaH-like phosphatase, a bacterial derived class of enzymes that are widespread throughout eukaryotes.

The Effect of Codon Mismatch on the Protein Translation System

Incorrect protein translation, caused by codon mismatch, is an important problem of living cells. In this work, a computational model was introduced to quantify the effects of codon mismatch and the model was used to study the protein translation of Saccharomyces cerevisiae. According to simulation results, the probability of codon mismatch will increase when the supply of amino acids is unbalanced, and the longer is the codon sequence, the larger is the probability for incorrect translation to occur, making the synthesis of long peptide chain difficult. By comparing to simulation results without codon mismatch effects taken into account, the fraction of mRNAs with bound ribosome decrease faster along the mRNAs, making the 5’ ramp phenomenon more obvious. It was also found in our work that the premature mechanism resulted from codon mismatch can reduce the proportion of incorrect translation when the amino acid supply is extremely unbalanced, which is one possible source of high fidelity protein synthesis after peptidyl transfer.

Transcriptome profiling of Giardia intestinalis using strand-specific RNA-seq

Giardia intestinalis is a common cause of diarrheal disease and it consists of eight genetically distinct genotypes or assemblages (A-H). Only assemblages A and B infect humans and are suggested to represent two different Giardia species. Correlations exist between assemblage type and host-specificity and to some extent symptoms. Phenotypical differences have been documented between assemblages and genome sequences are available for A, B and E. We have characterized and compared the polyadenylated transcriptomes of assemblages A, B and E. Four genetically different isolates were studied (WB (AI), AS175 (AII), P15 (E) and GS (B)) using paired-end, strand-specific RNA-seq. Most of the genome was transcribed in trophozoites grown in vitro, but at vastly different levels. RNA-seq confirmed many of the present annotations and refined the current genome annotation. Gene expression divergence was found to recapitulate the known phylogeny, and uncovered lineage-specific differences in expression. Polyadenylation sites were mapped for over 70% of the genes and revealed many examples of conserved and unexpectedly long 3′ UTRs. 28 open reading frames were found in a non-transcribed gene cluster on chromosome 5 of the WB isolate. Analysis of allele-specific expression revealed a correlation between allele-dosage and allele expression in the GS isolate. Previously reported cis-splicing events were confirmed and global mapping of cis-splicing identified only one novel intron. These observations can possibly explain differences in host-preference and symptoms, and it will be the basis for further studies of Giardia pathogenesis and biology.

Giardia is a single cell intestinal parasite and a common cause of diarrhea in humans and animals. Giardia is an unusual eukaryote by possessing two nuclei, a highly reduced genome and simple transcriptional apparatus. We have characterized the transcriptome of Giardia at single nucleotide resolution, which allowed the calculation of digital gene expression values for the complete set of genes. We performed a comparison of gene expression divergence across three genotypes. Most of the genes were transcribed, and the data were used to refine and correct gene models. Several gene expression differences were identified between the genotypes. A non-transcribed cluster of genes was detected on chromosome 5, likely representing a silenced region. The data also allowed mapping of transcript termini, which provided the first global view of 3′ untranslated regions in this parasite. This study also gives the first genome-wide evidence of transcription of allelic variants in Giardia. In this study, we provide novel insights into the transcriptome of an important human pathogen and model eukaryote. The findings reported here likely relate to the lifestyle of this parasite and its adaptation to parasitism. The data provide starting points for functional investigation of Giardia's biology and diplomonads generally.