Transcriptome analyses of the Giardia lamblia life cycle

Giardia duodenalis enolase is secreted as monomer during trophozoite-epithelial cell interactions, activates plasminogen and induces necroptotic damage

Enolase, a multifunctional protein expressed by multiple pathogens activates plasminogen to promote proteolysis on components of the extracellular matrix, an important event in early host-pathogen interactions. A secreted form of enolase that is released upon the interaction of trophozoites with epithelial cells has been detected in the secretome of G. duodenalis. However, the role of enolase in the host-pathogen interactions remains largely unknown. In this work, the effects of G. duodenalis enolase (Gd-eno) on the epithelial cell model (IEC-6) were analyzed. Firstly, the coding sequence of Giardia enolase was cloned and the recombinant protein used to raise antibodies that were then used to define the localization and role of enolase in epithelial cell-trophozoite interactions. Gd-eno was detected in small cytoplasmic vesicles as well as at the surface and is enriched in the region of the ventral disk of Giardia trophozoites. Moreover, the blocking of the soluble monomeric form of the enzyme, which is secreted upon interaction with IEC-6 cells by the anti-rGd-eno antibodies, significantly inhibited trophozoite attachment to intestinal IEC-6 cell monolayers. Further, rGd-eno was able to bind human plasminogen (HsPlg) and enhanced plasmin activity in vitro when the trophozoites were incubated with the intrinsic plasminogen activators of epithelial cells. In IEC-6 cells, rGd-eno treatment induced a profuse cell damage characterized by copious vacuolization, intercellular separation and detachment from the substrate; this effect was inhibited by either anti-Gd-eno Abs or the plasmin inhibitor ϵ- aminocaproic acid. Lastly, we established that in epithelial cells rGd-eno treatment induced a necroptotic-like process mediated by tumor necrosis factor α (TNF-α) and the apoptosis inducing factor (AIF), but independent of caspase-3. All together, these results suggest that Giardia enolase is a secreted moonlighting protein that stimulates a necroptotic-like process in IEC-6 epithelial cells via plasminogen activation along to TNFα and AIF activities and must be considered as a virulence factor.

Characterization of Metronidazole-Resistant Giardia intestinalis Lines by Comparative Transcriptomics and Proteomics

Metronidazole (MTZ) is a clinically important antimicrobial agent that is active against both bacterial and protozoan organisms. MTZ has been used extensively for more than 60 years and until now resistance has been rare. However, a recent and dramatic increase in the number of MTZ resistant bacteria and protozoa is of great concern since there are few alternative drugs with a similarly broad activity spectrum. To identify key factors and mechanisms underlying MTZ resistance, we utilized the protozoan parasite Giardia intestinalis, which is commonly treated with MTZ. We characterized two in vitro selected, metronidazole resistant parasite lines, as well as one revertant, by analyzing fitness aspects associated with increased drug resistance and transcriptomes and proteomes. We also conducted a meta-analysis using already existing data from additional resistant G. intestinalis isolates. The combined data suggest that in vitro generated MTZ resistance has a substantial fitness cost to the parasite, which may partly explain why resistance is not widespread despite decades of heavy use. Mechanistically, MTZ resistance in Giardia is multifactorial and associated with complex changes, yet a core set of pathways involving oxidoreductases, oxidative stress responses and DNA repair proteins, is central to MTZ resistance in both bacteria and protozoa.

A Detailed Gene Expression Map of Giardia Encystation

Giardia intestinalis is an intestinal protozoan parasite that causes diarrheal infections worldwide. A key process to sustain its chain of transmission is the formation of infectious cysts in the encystation process. We combined deep RNAseq of a broad range of encystation timepoints to produce a high-resolution gene expression map of Giardia encystation. This detailed transcriptomic map of encystation confirmed a gradual change of gene expression along the time course of encystation, showing the most significant gene expression changes during late encystation. Few genes are differentially expressed early in encystation, but the major cyst wall proteins CWP-1 and -2 are highly up-regulated already after 3.5 h encystation. Several transcription factors are sequentially up-regulated throughout the process, but many up-regulated genes at 7, 10, and 14 h post-induction of encystation have binding sites in the upstream regions for the Myb2 transcription factor, suggesting that Myb2 is a master regulator of encystation. We observed major changes in gene expression of several meiotic-related genes from 10.5 h of encystation to the cyst stage, and at 17.5 h encystation, there are changes in many different metabolic pathways and protein synthesis. Late encystation, 21 h to cysts, show extensive gene expression changes, most of all in VSP and HCMP genes, which are involved in antigenic variation, and genes involved in chromatin modifications. This high-resolution gene expression map of Giardia encystation will be an important tool in further studies of this important differentiation process.

Giardia duodenalis: Biology and Pathogenesis

Giardia duodenalis captured the attention of Leeuwenhoek in 1681 while he was examining his own diarrheal stool, but, ironically, it did not really gain attention as a human pathogen until the 1960s, when outbreaks were reported. Key technological advances, including in vitro cultivation, genomic and proteomic databases, and advances in microscopic and molecular approaches, have led to an understanding that this is a eukaryotic organism with a reduced genome rather than a truly premitochondriate eukaryote. This has included the discovery of mitosomes (vestiges of mitochondria), a transport system with many of the features of the Golgi apparatus, and even evidence for a sexual or parasexual cycle. Cell biology approaches have led to a better understanding of how Giardia survives with two nuclei and how it goes through its life cycle as a noninvasive organism in the hostile environment of the lumen of the host intestine. Studies of its immunology and pathogenesis have moved past the general understanding of the importance of the antibody response in controlling infection to determining the key role of the Th17 response. This work has led to understanding of the requirement for a balanced host immune response that avoids the extremes of an excessive response with collateral damage or one that is unable to clear the organism. This understanding is especially important in view of the remarkable ranges of early manifestations, which range from asymptomatic to persistent diarrhea and weight loss, and longer-term sequelae that include growth stunting in children who had no obvious symptoms and a high frequency of postinfectious irritable bowel syndrome (IBS).

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.

Click chemistry-facilitated comprehensive identification of proteins adducted by antimicrobial 5-nitroimidazoles for discovery of alternative drug targets against giardiasis

Giardiasis and other protozoan infections are major worldwide causes of morbidity and mortality, yet development of new antimicrobial agents with improved efficacy and ability to override increasingly common drug resistance remains a major challenge. Antimicrobial drug development typically proceeds by broad functional screens of large chemical libraries or hypothesis-driven exploration of single microbial targets, but both strategies have challenges that have limited the introduction of new antimicrobials. Here, we describe an alternative drug development strategy that identifies a sufficient but manageable number of promising targets, while reducing the risk of pursuing targets of unproven value. The strategy is based on defining and exploiting the incompletely understood adduction targets of 5-nitroimidazoles, which are proven antimicrobials against a wide range of anaerobic protozoan and bacterial pathogens. Comprehensive adductome analysis by modified click chemistry and multi-dimensional proteomics were applied to the model pathogen Giardia lamblia to identify dozens of adducted protein targets common to both 5'-nitroimidazole-sensitive and -resistant cells. The list was highly enriched for known targets in G. lamblia, including arginine deiminase, α-tubulin, carbamate kinase, and heat shock protein 90, demonstrating the utility of the approach. Importantly, over twenty potential novel drug targets were identified. Inhibitors of two representative new targets, NADP-specific glutamate dehydrogenase and peroxiredoxin, were found to have significant antigiardial activity. Furthermore, all the identified targets remained available in resistant cells, since giardicidal activity of the respective inhibitors was not impacted by resistance to 5'-nitroimidazoles. These results demonstrate that the combined use of click chemistry and proteomics has the potential to reveal alternative drug targets for overcoming antimicrobial drug resistance in protozoan parasites.

Role of gamma-giardin in ventral disc formation of Giardia lamblia

Background

Giardia lamblia, a protozoan pathogen causing diarrheal outbreaks, has characteristic cytoskeletal structures including eight flagella, a median body and a ventral disc. Gamma-giardin is a unique component protein of the cytoskeleton of this protozoan.

Results

Through comparative proteomic analysis between different stages of the cell cycle, G. lamblia γ-giardin (Glγ-giardin) was identified as an upregulated protein in the G2-phase. Increased Glγ-giardin expression in G2 was confirmed by western blot and real-time polymerase chain reaction analyses. Knockdown of this protein using a morpholino affected the formation of ventral discs, especially the microribbons of the discs, but exerted little effect on the binding ability of G. lamblia. The number of cells with four nuclei was increased in Glγ-giardin-knockdown cells. Expression of Glγ-giardin was decreased during encystation, in contrast with the G2-phase.

Conclusions

Knockdown experiments demonstrated that Glγ-giardin is a component of the trilaminar structure of the ventral disc. Expression of Glγ-giardin is induced in the G2-phase prior to active cell division, whereas its expression decreases during encystation, a dormant stage of G. lamblia.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3478-8) contains supplementary material, which is available to authorized users.

Sense-antisense gene overlap is probably a cause for retaining the few introns in Giardia genome and the implications

Background

It is widely accepted that the last eukaryotic common ancestor and early eukaryotes were intron-rich and intron loss dominated subsequent evolution, thus the presence of only very few introns in some modern eukaryotes must be the consequence of massive loss. But it is striking that few eukaryotes were found to have completely lost introns. Despite extensive research, the causes of massive intron losses remain elusive. Actually the reverse question -- how the few introns can be retained under the evolutionary selection pressure of intron loss -- is equally significant but was rarely studied, except that it was conjectured that the essential functions of some introns prevent their loss. The situation that extremely few (eight) spliceosome-mediated cis-spliced introns present in the relatively simple genome of Giardia lamblia provides an excellent opportunity to explore this question.

Results

Our investigation found three types of distribution patterns of the few introns in the intron-containing genes: ancient intron in ancient gene, later-evolved intron in ancient gene, and later-evolved intron in later-evolved gene, which can reflect to some extent the dynamic evolution of introns in Giardia. Without finding any special features or functional importance of these introns responsible for their retention, we noticed and experimentally verified that some intron-containing genes form sense-antisense gene pairs with transcribable genes on their complementary strands, and that the introns just reside in the overlapping regions.

Conclusions

In Giardia’s evolution, despite constant evolutionary selection pressure of intron loss, intron gain can still occur in both ancient and later-evolved genes, but only a few introns are retained; at least the evolutionary retention of some of the introns might not be due to the functional constraint of the introns themselves but the causes outside of introns, such as the constraints imposed by other genomic functional elements overlapping with the introns. These findings can not only provide some clues to find new genomic functional elements -- in the areas overlapping with introns, but suggest that “functional constraint” of introns may not be necessarily directly associated with intron loss and gain, and that the real functions are probably still outside of our current knowledge.

Reviewers

This article was reviewed by Mikhail Gelfand, Michael Gray, and Igor Rogozin.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0226-5) contains supplementary material, which is available to authorized users.

Recent advances in the genomic and molecular biology of Giardia

Giardia duodenalis is the most common gastrointestinal protozoan parasite of humans and a significant contributor to the global burden of both diarrheal disease and post-infectious chronic disorders. Robust tools for analyzing gene function in this parasite have been developed and a range of genetic tools are now available. These together with public databases have provided insights on the function of different genes in Giardia. In this review we provide a current perspective on different molecular aspects of Giardia related to genomics, regulation of encystation, trophozoite transcriptional responses to physiological and xenobiotic (drug-induced) stress, and mechanisms of drug resistance. We also examine recent insights that have contributed to gain knowledge in the study of VSPs, antigenic variation, epigenetics, DNA repair and in the direct manipulation of gene function in Giardia, with a particular focus on the inducible Cre/loxP system.

A new landscape of host-protozoa interactions involving the extracellular vesicles world

Extracellular vesicles (EVs) are released by a wide number of cells including blood cells, immune system cells, tumour cells, adult and embryonic stem cells. EVs are a heterogeneous group of vesicles (~30-1000 nm) including microvesicles and exosomes. The physiological release of EVs represents a normal state of the cell, raising a metabolic equilibrium between catabolic and anabolic processes. Moreover, when the cells are submitted to stress with different inducers or in pathological situations (malignancies, chronic diseases, infectious diseases.), they respond with an intense and dynamic release of EVs. The EVs released from stimulated cells vs those that are released constitutively may themselves differ, both physically and in their cargo. EVs contain protein, lipids, nucleic acids and biomolecules that can alter cell phenotypes or modulate neighbouring cells. In this review, we have summarized findings involving EVs in certain protozoan diseases. We have commented on strategies to study the communicative roles of EVs during host-pathogen interaction and hypothesized on the use of EVs for diagnostic, preventative and therapeutic purposes in infectious diseases. This kind of communication could modulate the innate immune system and reformulate concepts in parasitism. Moreover, the information provided within EVs could produce alternatives in translational medicine.

Glucosylceramide transferase in Giardia preferentially catalyzes the synthesis of galactosylceramide during encystation

The stage differentiation from trophozoite to cyst (i.e., encystation) is an essential step for Giardia to survive outside its human host and spread the infection via the fecal-oral route. We have previously shown that Giardia expresses glucosylceramide transferase 1 (GlcT1) enzyme, the activity of which is elevated during encystation. We have also reported that blocking the activity of gGlcT1 interferes with the biogenesis of encystation-specific vesicles (ESVs) and cyst viability in Giardia. To further understand the role of this enzyme and how it regulates encystation, we overexpressed, knocked down, and rescued the giardial GlcT1 (gGlcT1) gene and measured its enzymatic activity in live parasites as well as in isolated membrane fractions using NBD-ceramide and UDP-glucose or UDP-galactose. We observed that gGlcT1 is able to catalyze the synthesis of both glucosylceramide (GlcCer) and galactosylceramide (GalCer), however the synthesis of GalCer is 2-3 fold higher than of GlcCer. Although both activities follow Michaelis-Menten kinetics, the bindings of UDP-glucose and UDP-galactose with the enzyme appear to be non-competitive and independent of each other. The modulation of gGlcT1 synthesis concomitantly influenced the expression cyst-wall protein (CWP) and overall encystation. We propose that gGlcT1 is a unique enzyme and that Giardia uses this enzyme to synthesize both GlcCer and GalCer to facilitate the process of encystation/cyst production.

Rac Regulates Giardia lamblia Encystation by Coordinating Cyst Wall Protein Trafficking and Secretion

Encystation of the common intestinal parasite Giardia lamblia involves the production, trafficking, and secretion of cyst wall material (CWM). However, the molecular mechanism responsible for the regulation of these sequential processes remains elusive. Here, we examined the role of GlRac, Giardia’s sole Rho family GTPase, in the regulation of endomembrane organization and cyst wall protein (CWP) trafficking. Localization studies indicated that GlRac is associated with the endoplasmic reticulum (ER) and the Golgi apparatus-like encystation-specific vesicles (ESVs). Constitutive GlRac signaling increased levels of the ER marker PDI2, induced ER swelling, reduced overall CWP1 production, and promoted the early maturation of ESVs. Quantitative analysis of cells expressing constitutively active hemagglutinin (HA)-tagged GlRac (HA-Rac^CA) revealed fewer but larger ESVs than control cells. Consistent with the phenotype of premature maturation of ESVs in HA-Rac^CA-expressing cells, constitutive GlRac signaling resulted in increased CWP1 secretion and, conversely, morpholino depletion of GlRac blocked CWP1 secretion. Wild-type cells unexpectedly secreted large quantities of CWP1 into the medium, and free CWP1 was used cooperatively during cyst formation. These results, in part, could account for the previously reported observation that G. lamblia encysts more efficiently at high cell densities. These studies of GlRac show that it regulates encystation at several levels, and our findings support its coordinating role as a regulator of CWP trafficking and secretion. The central role of GlRac in regulating membrane trafficking and the cytoskeleton, both of which are essential to Giardia parasitism, further suggests its potential as a novel target for drug development to treat giardiasis.

The encystation process is crucial for the transmission of giardiasis and the life cycle of many protists. Encystation for Giardia lamblia involves the assembly of a protective cyst wall via sequential production, trafficking, and secretion of cyst wall material. However, the regulatory pathways that coordinate cargo maturation and secretion remain unknown. Here, we asked whether the signaling activities of G. lamblia’s single Rho family GTPase, GlRac, might have a regulatory role in the encystation process. We show that GlRac localizes to endomembranes and its signaling activities regulate the production of cyst wall protein 1 (CWP1), the maturation of encystation-specific vesicles (ESVs), and secretion of CWP1. We also show that secreted CWP1 is available for the development of cysts at the population level, a finding that in part could explain why Giardia encystation proceeds more efficiently at high cell densities.

Coordinated Changes in Gene Expression Throughout Encystation of Giardia intestinalis

Differentiation into infectious cysts through the process of encystation is crucial for transmission and survival of the intestinal protozoan parasite Giardia intestinalis. Hitherto the majority of studies have focused on the early events, leaving late encystation poorly defined. In order to further study encystation, focusing on the later events, we developed a new encystation protocol that generates a higher yield of mature cysts compared to standard methods. Transcriptome changes during the entire differentiation from trophozoites to cysts were thereafter studied using RNA sequencing (RNA-seq). A high level of periodicity was observed for up- and down-regulated genes, both at the level of the entire transcriptome and putative regulators. This suggests the trajectory of differentiation to be coordinated through developmentally linked gene regulatory activities. Our study identifies a core of 13 genes that are consistently up-regulated during initial encystation. Of these, two constitute previously uncharacterized proteins that we were able to localize to a new type of encystation-specific vesicles. Interestingly, the largest transcriptional changes were seen in the late phase of encystation with the majority of the highly up-regulated genes encoding hypothetical proteins. Several of these were epitope-tagged and localized to further characterize these previously unknown genetic components of encystation and possibly excystation. Finally, we also detected a switch of variant specific surface proteins (VSPs) in the late phase of encystation. This occurred at the same time as nuclear division and DNA replication, suggesting a potential link between the processes.

The intestinal protozoan parasite Giardia intestinalis and many other medically important protozoan parasites must encyst and form infective cysts in order to transmit to new hosts. Encystation efficiency is in that way connected to efficiency of transmission. We have developed new in vitro differentiation protocols and made the first RNA-seq based gene expression study of the complete Giardia encystation process. Our data provides a road map of Giardia encystation and a starting point from where it is possible to further explore important processes occurring during encystation. Information about this vital process for survival in the environment of this and other cyst forming parasites can be used in the development of new types of interventions.

Constitutive aneuploidy and genomic instability in the single-celled eukaryote Giardia intestinalis

Giardia intestinalis is an important single‐celled human pathogen. Interestingly, this organism has two equal‐sized transcriptionally active nuclei, each considered diploid. By evaluating condensed chromosome numbers and visualizing homologous chromosomes by fluorescent in situ hybridization, we determined that the Giardia cells are constitutively aneuploid. We observed karyotype inter‐and intra‐population heterogeneity in eight cell lines from two clinical isolates, suggesting constant karyotype evolution during in vitro cultivation. High levels of chromosomal instability and frequent mitotic missegregations observed in four cell lines correlated with a proliferative disadvantage and growth retardation. Other cell lines, although derived from the same clinical isolate, revealed a stable yet aneuploid karyotype. We suggest that both chromatid missegregations and structural rearrangements contribute to shaping the Giardia genome, leading to whole‐chromosome aneuploidy, unequal gene distribution, and a genomic divergence of the two nuclei within one cell. Aneuploidy in Giardia is further propagated without p53‐mediated cell cycle arrest and might have been a key mechanism in generating the genetic diversity of this human pathogen.

Sphingolipids, Lipid Rafts, and Giardial Encystation: The Show Must Go On

Sphingolipids are sphingosine-based phospholipids, which are present in the plasma and endomembranes of many eukaryotic cells. These lipids are involved in various cellular functions, including cell growth, differentiation, and apoptosis. In addition, sphingolipid and cholesterol-enriched membrane microdomains (also called "lipid rafts") contain a set of proteins and lipids, which take part in the signaling process in response to intra- or extracellular stimuli. Recent findings suggest that sphingolipids, especially glucosylceramide, play a critical role in inducing encystation and maintaining the cyst viability in Giardia. Similarly, the assembly/disassembly of lipid rafts modulates the encystation and cyst production of this ubiquitous enteric parasite. In this review article, we discuss the overall progress in the field and examine whether sphingolipids and lipid rafts can be used as novel targets for designing therapies to control infection by Giardia, which is rampant in developing countries, where children are especially vulnerable.

UV irradiation responses in Giardia intestinalis

The response to ultraviolet light (UV) radiation, a natural stressor to the intestinal protozoan parasite Giardia intestinalis, was studied to deepen the understanding of how the surrounding environment affects the parasite during transmission. UV radiation at 10 mJ/cm(2) kills Giardia cysts effectively whereas trophozoites and encysting parasites can recover from UV treatment at 100 mJ/cm(2) and 50 mJ/cm(2) respectively. Staining for phosphorylated histone H2A showed that UV treatment induces double-stranded DNA breaks and flow cytometry analyses revealed that UV treatment of trophozoites induces DNA replication arrest. Active DNA replication coupled to DNA repair could be an explanation to why UV light does not kill trophozoites and encysting cells as efficiently as the non-replicating cysts. We also examined UV-induced gene expression responses in both trophozoites and cysts using RNA sequencing (RNA seq). UV radiation induces small overall changes in gene expression in Giardia but cysts show a stronger response than trophozoites. Heat shock proteins, kinesins and Nek kinases are up-regulated, whereas alpha-giardins and histones are down-regulated in UV treated trophozoites. Expression of variable surface proteins (VSPs) is changed in both trophozoites and cysts. Our data show that Giardia cysts have limited ability to repair UV-induced damage and this may have implications for drinking- and waste-water treatment when setting criteria for the use of UV disinfection to ensure safe water.

Gene expression changes during Giardia-host cell interactions in serum-free medium

Serial Analysis of Gene Expression (SAGE) was used to quantify transcriptional changes in Giardia intestinalis during its interaction with human intestinal epithelial cells (IECs, HT-29) in serum free M199 medium. Transcriptional changes were compared to those in trophozoites alone in M199 and in TYI-S-33 Giardia growth medium. In total, 90 genes were differentially expressed, mainly those involved in cellular redox homeostasis, metabolism and small molecule transport but also cysteine proteases and structural proteins of the giardin family. Only 29 genes changed their expression due to IEC interaction and the rest were due to M199 medium. Although our findings generated a small dataset, it was consistent with our earlier microarray studies performed under different interaction conditions. This study has confined the number of genes in Giardia to a small subset that specifically change their expression due to interaction with IECs.

Identification of Giardia lamblia DHHC proteins and the role of protein S-palmitoylation in the encystation process

Protein S-palmitoylation, a hydrophobic post-translational modification, is performed by protein acyltransferases that have a common DHHC Cys-rich domain (DHHC proteins), and provides a regulatory switch for protein membrane association. In this work, we analyzed the presence of DHHC proteins in the protozoa parasite Giardia lamblia and the function of the reversible S-palmitoylation of proteins during parasite differentiation into cyst. Two specific events were observed: encysting cells displayed a larger amount of palmitoylated proteins, and parasites treated with palmitoylation inhibitors produced a reduced number of mature cysts. With bioinformatics tools, we found nine DHHC proteins, potential protein acyltransferases, in the Giardia proteome. These proteins displayed a conserved structure when compared to different organisms and are distributed in different monophyletic clades. Although all Giardia DHHC proteins were found to be present in trophozoites and encysting cells, these proteins showed a different intracellular localization in trophozoites and seemed to be differently involved in the encystation process when they were overexpressed. dhhc transgenic parasites showed a different pattern of cyst wall protein expression and yielded different amounts of mature cysts when they were induced to encyst. Our findings disclosed some important issues regarding the role of DHHC proteins and palmitoylation during Giardia encystation.

Giardiasis is a major cause of non-viral/non-bacterial diarrheal disease worldwide and has been included within the WHO Neglected Disease Initiative since 2004. Infection begins with the ingestion of Giardia lamblia in cyst form, which, after exposure to gastric acid in the host stomach and proteases in the duodenum, gives rise to trophozoites. The inverse process is called encystation and begins when the trophozoites migrate to the lower part of the small intestine where they receive signals that trigger synthesis of the components of the cyst wall. The cyst form enables the parasite to survive in the environment, infect a new host and evade the immune response. In this work, we explored the role of protein S-palmitoylation, a unique reversible post-translational modification, during Giardia encystation, because de novo generation of endomembrane compartments, protein sorting and vesicle fusion occur in this process. Our findings may contribute to the design of therapeutic agents against this important human pathogen.

Proteomics of secretory and endocytic organelles in Giardia lamblia

Giardia lamblia is a flagellated protozoan enteroparasite transmitted as an environmentally resistant cyst. Trophozoites attach to the small intestine of vertebrate hosts and proliferate by binary fission. They access nutrients directly via uptake of bulk fluid phase material into specialized endocytic organelles termed peripheral vesicles (PVs), mainly on the exposed dorsal side. When trophozoites reach the G2/M restriction point in the cell cycle they can begin another round of cell division or encyst if they encounter specific environmental cues. They induce neogenesis of Golgi-like organelles, encystation-specific vesicles (ESVs), for regulated secretion of cyst wall material. PVs and ESVs are highly simplified and thus evolutionary diverged endocytic and exocytic organelle systems with key roles in proliferation and transmission to a new host, respectively. Both organelle systems physically and functionally intersect at the endoplasmic reticulum (ER) which has catabolic as well as anabolic functions. However, the unusually high degree of sequence divergence in Giardia rapidly exhausts phylogenomic strategies to identify and characterize the molecular underpinnings of these streamlined organelles. To define the first proteome of ESVs and PVs we used a novel strategy combining flow cytometry-based organelle sorting with in silico filtration of mass spectrometry data. From the limited size datasets we retrieved many hypothetical but also known organelle-specific factors. In contrast to PVs, ESVs appear to maintain a strong physical and functional link to the ER including recruitment of ribosomes to organelle membranes. Overall the data provide further evidence for the formation of a cyst extracellular matrix with minimal complexity. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium with the dataset identifier PXD000694.

The proteome landscape of Giardia lamblia encystation

Giardia lamblia is an intestinal protozoan parasite required to survive in the environment in order to be transmitted to a new host. To ensure parasite survival, flagellated trophozoites colonizing the small intestine differentiate into non-motile environmentally-resistant cysts which are then shed in the environment. This cell differentiation process called encystation is characterized by significant morphological remodeling which includes secretion of large amounts of cyst wall material. Although much is known about the transcriptional regulation of encystation and the synthesis and trafficking of cyst wall material, the investigation of global changes in protein content and abundance during G. lamblia encystation is still unaddressed.

In this study, we report on the quantitative analysis of the G. lamblia proteome during encystation using tandem mass spectrometry. Quantification of more than 1000 proteins revealed major changes in protein abundance in early, mid and late encystation, notably in constitutive secretory protein trafficking. Early stages of encystation were marked by a striking decrease of endoplasmic reticulum-targeted variant-specific surface proteins and significant increases in cytoskeleton regulatory components, NEK protein kinases and proteins involved in protein folding and glycolysis. This was in stark contrast to cells in the later stages of encystation which presented a surprisingly similar proteome composition to non-encysting trophozoites. Altogether these data constitute the first quantitative atlas of the Giardia proteome covering the whole process of encystation and point towards an important role for post-transcriptional control of gene expression in Giardia differentiation. Furthermore, our data provide a valuable resource for the community-based annotation effort of the G. lamblia genome, where almost 70% of all predicted gene models remains “hypothetical”.

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.

Putative SF2 helicases of the early-branching eukaryote Giardia lamblia are involved in antigenic variation and parasite differentiation into cysts

Background

Regulation of surface antigenic variation in Giardia lamblia is controlled post-transcriptionally by an RNA-interference (RNAi) pathway that includes a Dicer-like bidentate RNase III (gDicer). This enzyme, however, lacks the RNA helicase domain present in Dicer enzymes from higher eukaryotes. The participation of several RNA helicases in practically all organisms in which RNAi was studied suggests that RNA helicases are potentially involved in antigenic variation, as well as during Giardia differentiation into cysts.

Results

An extensive in silico analysis of the Giardia genome identified 32 putative Super Family 2 RNA helicases that contain almost all the conserved RNA helicase motifs. Phylogenetic studies and sequence analysis separated them into 22 DEAD-box, 6 DEAH-box and 4 Ski2p-box RNA helicases, some of which are homologs of well-characterized helicases from higher organisms. No Giardia putative helicase was found to have significant homology to the RNA helicase domain of Dicer enzymes. Additionally a series of up- and down-regulated putative RNA helicases were found during encystation and antigenic variation by qPCR experiments. Finally, we were able to recognize 14 additional putative helicases from three different families (RecQ family, Swi2/Snf2 and Rad3 family) that could be considered DNA helicases.

Conclusions

This is the first comprehensive analysis of the Super Family 2 helicases from the human intestinal parasite G. lamblia. The relative and variable expression of particular RNA helicases during both antigenic variation and encystation agrees with the proposed participation of these enzymes during both adaptive processes. The putatives RNA and DNA helicases identified in this early-branching eukaryote provide initial information regarding the biological role of these enzymes in cell adaptation and differentiation.

Protozoan infections of the gastrointestinal tract

The review provides current views on human protozoan parasites of the gut. The recognition of the importance of cryptosporidium, cyclospora and microsporidia over the last three decades emphasises the possibility that more pathogenic intestinal protozoa are presently unrecognized. Each of these is a zoonotic infection and the potential for a zoonotic element to the transmission of giardiasis has been recognized. A common theme in increased understanding of the biology and pathological mechanisms involved in causing disease is the application of molecular techniques to the various stages of the parasite life cycle. Molecular methods are increasingly contributing to laboratory diagnosis of these conditions with increased yields of positive results though in the tropics it is likely that fecal microscopy will remain the standard for some time to come. The nitroimidazole compounds are the mainstay of treatment for giardia and amebiasis with no major advance in therapeutics since their role was appreciated. Nitazoxanide was shown to be effective for cryptosporidiosis in the 1990s.

Giardia--from genome to proteome

In this review, the current status of genomic and proteomic research on Giardia is examined in terms of evolutionary biology, phylogenetic relationships and taxonomy. The review also describes how characterising genetic variation in Giardia from numerous hosts and endemic areas has provided a better understanding of life cycle patterns, transmission and the epidemiology of Giardia infections in humans, domestic animals and wildlife. Some progress has been made in relating genomic information to the phenotype of Giardia, and as a consequence, new information has been obtained on aspects of developmental biology and the host-parasite relationship. However, deficiencies remain in our understanding of pathogenesis and host specificity, highlighting the limitations of currently available genomic datasets.

Plasmid vectors for proteomic analyses in Giardia: purification of virulence factors and analysis of the proteasome

In recent years, proteomics has come of age with the development of efficient tools for purification, identification, and characterization of gene products predicted by genome projects. The intestinal protozoan Giardia intestinalis can be transfected, but there is only a limited set of vectors available, and most of them are not user friendly. This work delineates the construction of a suite of cassette-based expression vectors for use in Giardia. Expression is provided by the strong constitutive ornithine carbamoyltransferase (OCT) promoter, and tagging is possible in both N- and C-terminal configurations. Taken together, the vectors are capable of providing protein localization and production of recombinant proteins, followed by efficient purification by a novel affinity tag combination, streptavidin binding peptide-glutathione S-transferase (SBP-GST). The option of removing the tags from purified proteins was provided by the inclusion of a PreScission protease site. The efficiency and feasibility of producing and purifying endogenous recombinant Giardia proteins with the developed vectors was demonstrated by the purification of active recombinant arginine deiminase (ADI) and OCT from stably transfected trophozoites. Moreover, we describe the tagging, purification by StrepTactin affinity chromatography, and compositional analysis by mass spectrometry of the G. intestinalis 26S proteasome by employing the Strep II-FLAG-tandem affinity purification (SF-TAP) tag. This is the first report of efficient production and purification of recombinant proteins in and from Giardia, which will allow the study of specific parasite proteins and protein complexes.

Importance of enolase in Giardia lamblia differentiation

The ability of Giardia to differentiate into cysts which survive in the environment and release the virulent trophozoites after ingestion in the small intestine is essential for transmission and disease. We examined the role of enolase, a glycolytic enzyme, in Giardia differentiation. The sequence of Giardia lamblia enolase (gEno) is most similar to enolases in Homo sapiens and Leishmania mexicana, and shows the conserved catalytic and metal-binding residues. We used an integration vector to stably express wild type and mutant gEno. In trophozoites, wild type gEno localized to the cell membrane, caudal flagella and cytosol. gEno is present on the wall of mature cysts, but not in encystation secretory vesicles (ESV). The expression of gEno with a deletion of residues G167-K169, or mutations H389Q/R390S significantly inhibited excystation while mutation of residue D257K had no effect. These results suggest a role for enolase in regulation of Giardia excystation.

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

BACKGROUND

RNA degradation is critical to the survival of all cells. With increasing evidence for pervasive transcription in cells, RNA degradation has gained recognition as a means of regulating gene expression. Yet, RNA degradation machinery has been studied extensively in only a few eukaryotic organisms, including Saccharomyces cerevisiae and humans. Giardia lamblia is a parasitic protist with unusual genomic traits: it is binucleated and tetraploid, has a very compact genome, displays a theme of genomic minimalism with cellular machinery commonly comprised of a reduced number of protein components, and has a remarkably large population of long, stable, noncoding, antisense RNAs.

RESULTS

Here we use in silico approaches to investigate the major RNA degradation machinery in Giardia lamblia and compare it to a broad array of other parasitic protists. We have found key constituents of the deadenylation and decapping machinery and of the 5'-3' RNA degradation pathway. We have similarly found that all of the major 3'-5' RNA degradation pathways are present in Giardia, including both exosome-dependent and exosome-independent machinery. However, we observe significant loss of RNA degradation machinery genes that will result in important differences in the protein composition, and potentially functionality, of the various RNA degradation pathways. This is most apparent in the exosome, the central mediator of 3'-5' degradation, which apparently contains an altered core configuration in both Giardia and Plasmodium, with only four, instead of the canonical six, distinct subunits. Additionally the exosome in Giardia is missing both the Rrp6, Nab3, and Nrd1 proteins, known to be key regulators of noncoding transcript stability in other cells.

CONCLUSIONS

These findings suggest that although the full complement of the major RNA degradation mechanisms were present - and likely functional - early in eukaryotic evolution, the composition and function of the complexes is more variable than previously appreciated. We suggest that the missing components of the exosome complex provide an explanation for the stable abundance of sterile RNA species in Giardia.

Structure of a cyclin-dependent kinase from Giardia lamblia

Giardia lamblia is the etiologic agent of giardiasis, a water-borne infection that is prevalent throughout the world. The need for new therapeutics for the treatment of giardiasis is of paramount importance. Owing to the ubiquitous nature of kinases and their vital importance in organisms, they are potential drug targets. In this paper, the first structure of a cyclin-dependent kinase (CDK) from G. lamblia (GlCDK; UniProt A8BZ95) is presented. CDKs are cell-cycle-associated kinases that are actively being pursued as targets for anticancer drugs as well as for antiparasitic chemotherapy. Generally, a CDK forms a complex with its associated cyclin. This CDK-cyclin complex is active and acts as a serine/threonine protein kinase. Typically, CDKs are responsible for the transition to the next phase of the cell cycle. Although the structure of GlCDK with its associated cyclin was not solved, the 1.85 Å resolution structure of apo GlCDK and a 2.0 Å resolution structure of GlCDK in complex with adenosine monophosphate are presented and the structural differences from the orthologous human CDK2 and CDK3 are discussed.

The minimal kinome of Giardia lamblia illuminates early kinase evolution and unique parasite biology

Background

The major human intestinal pathogen Giardia lamblia is a very early branching eukaryote with a minimal genome of broad evolutionary and biological interest.

Results

To explore early kinase evolution and regulation of Giardia biology, we cataloged the kinomes of three sequenced strains. Comparison with published kinomes and those of the excavates Trichomonas vaginalis and Leishmania major shows that Giardia's 80 core kinases constitute the smallest known core kinome of any eukaryote that can be grown in pure culture, reflecting both its early origin and secondary gene loss. Kinase losses in DNA repair, mitochondrial function, transcription, splicing, and stress response reflect this reduced genome, while the presence of other kinases helps define the kinome of the last common eukaryotic ancestor. Immunofluorescence analysis shows abundant phospho-staining in trophozoites, with phosphotyrosine abundant in the nuclei and phosphothreonine and phosphoserine in distinct cytoskeletal organelles. The Nek kinase family has been massively expanded, accounting for 198 of the 278 protein kinases in Giardia. Most Neks are catalytically inactive, have very divergent sequences and undergo extensive duplication and loss between strains. Many Neks are highly induced during development. We localized four catalytically active Neks to distinct parts of the cytoskeleton and one inactive Nek to the cytoplasm.

Conclusions

The reduced kinome of Giardia sheds new light on early kinase evolution, and its highly divergent sequences add to the definition of individual kinase families as well as offering specific drug targets. Giardia's massive Nek expansion may reflect its distinctive lifestyle, biphasic life cycle and complex cytoskeleton.

Cyst and encystment in protozoan parasites: optimal targets for new life-cycle interrupting strategies?

Certain protozoan parasites use survival strategies to reside outside the host such as the formation of cysts. This dormant and resistant stage results from the complex process of encystment that involves diverse molecular and cellular modifications. The stimuli and changes associated with cyst biogenesis are a matter of ongoing studies in human and animal protozoan parasites such as amoeba and Giardia species because blocking every step in the encystment pathway should, in theory, interrupt their life cycles. The present review thoroughly examines this essential process in those protozoan parasites and discusses the possibility of using that information to develop new kinds of anti-parasite specific and life cycle-interrupting drugs, aimed at holding back the dissemination of these infections.

A comparison of the Giardia lamblia trophozoite and cyst transcriptome using microarrays

Background

Compared with many protists, Giardia lamblia has a simple life cycle alternating between cyst and trophozoite. Most research on the molecular biology of Giardia parasites has focused on trophozoites and the processes of excystation and encystation, whereas cysts have attracted less interest. The striking morphological differences between the dormant cyst and the rapidly dividing and motile trophozoite implies profound changes in the metabolism as the parasite encysts in the host's intestine and excysts upon ingestion by a new host.

Results

To investigate the magnitude of the transcriptional changes occurring during the G. lamblia life cycle we compared the transcriptome of G. lamblia trophozoites and cysts using single-color oligonucleotide microarrays. Cysts were found to possess a much smaller transcriptome, both in terms of mRNA diversity and abundance. Genes encoding proteins related to ribosomal functions are highly over-represented. The comparison of the transcriptome of cysts generated in culture or extracted from feces revealed little overlap, raising the possibility of significant biological differences between the two types of cysts.

Conclusions

The comparison of the G. lamblia cyst and trophozoite transcriptome showed that transcripts of most genes are present at a lower level in cysts. This global view of the cyst and trophozoite transcriptome complements studies focused on the expression of selected genes during trophozoite multiplication, encystation and excystation.

BiologicalNetworks 2.0--an integrative view of genome biology data

Background

A significant problem in the study of mechanisms of an organism's development is the elucidation of interrelated factors which are making an impact on the different levels of the organism, such as genes, biological molecules, cells, and cell systems. Numerous sources of heterogeneous data which exist for these subsystems are still not integrated sufficiently enough to give researchers a straightforward opportunity to analyze them together in the same frame of study. Systematic application of data integration methods is also hampered by a multitude of such factors as the orthogonal nature of the integrated data and naming problems.

Results

Here we report on a new version of BiologicalNetworks, a research environment for the integral visualization and analysis of heterogeneous biological data. BiologicalNetworks can be queried for properties of thousands of different types of biological entities (genes/proteins, promoters, COGs, pathways, binding sites, and other) and their relations (interactions, co-expression, co-citations, and other). The system includes the build-pathways infrastructure for molecular interactions/relations and module discovery in high-throughput experiments. Also implemented in BiologicalNetworks are the Integrated Genome Viewer and Comparative Genomics Browser applications, which allow for the search and analysis of gene regulatory regions and their conservation in multiple species in conjunction with molecular pathways/networks, experimental data and functional annotations.

Conclusions

The new release of BiologicalNetworks together with its back-end database introduces extensive functionality for a more efficient integrated multi-level analysis of microarray, sequence, regulatory, and other data. BiologicalNetworks is freely available at http://www.biologicalnetworks.org.

An atypical proprotein convertase in Giardia lamblia differentiation

Proteolytic activity is important in the lifecycles of parasites and their interactions with hosts. Cysteine proteases have been best studied in Giardia, but other protease classes have been implicated in growth and/or differentiation. In this study, we employed bioinformatics to reveal the complete set of putative proteases in the Giardia genome. We identified 73 peptidase homologs distributed over 5 catalytic classes in the genome. Serial analysis of gene expression of the G. lamblia lifecycle found thirteen protease genes with significant transcriptional variation over the lifecycle, with only one serine protease transcript upregulated late in encystation. The translated gene sequence of this encystation-specific transcript was most similar to eukaryotic subtilisin-like proprotein convertases (SPC), although the typical catalytic triad was not identified. Epitope-tagged gSPC protein expressed in Giardia under its own promoter was upregulated during encystation with highest expression in cysts and it localized to encystation-specific secretory vesicles (ESV). Total gSPC from encysting cells produced proteolysis in gelatin gels that co-migrated with the epitope-tagged protease in immunoblots. Immuno-purified gSPC also had gelatinase activity. To test whether endogenous gSPC activity is involved in differentiation, trophozoites and cysts were exposed to the specific serine proteinase inhibitor 4-(2-aminoethyl)-benzenesulfonyl fluoride hydrochloride (AEBSF). After 21 h encystation, a significant decrease in ESV was observed with 1mM AEBSF and by 42 h the number of cysts was significantly reduced, but trophozoite growth was not inhibited. Concurrently, levels of cyst wall proteins 1 and 2, and AU1-tagged gSPC protein itself were decreased. Excystation of G. muris cysts was also significantly reduced in the presence of AEBSF. These results support the idea that serine protease activity is essential for Giardia encystation and excystation.

Genome analysis and comparative genomics of a Giardia intestinalis assemblage E isolate

Background

Giardia intestinalis is a protozoan parasite that causes diarrhea in a wide range of mammalian species. To further understand the genetic diversity between the Giardia intestinalis species, we have performed genome sequencing and analysis of a wild-type Giardia intestinalis sample from the assemblage E group, isolated from a pig.

Results

We identified 5012 protein coding genes, the majority of which are conserved compared to the previously sequenced genomes of the WB and GS strains in terms of microsynteny and sequence identity. Despite this, there is an unexpectedly large number of chromosomal rearrangements and several smaller structural changes that are present in all chromosomes. Novel members of the VSP, NEK Kinase and HCMP gene families were identified, which may reveal possible mechanisms for host specificity and new avenues for antigenic variation. We used comparative genomics of the three diverse Giardia intestinalis isolates P15, GS and WB to define a core proteome for this species complex and to identify lineage-specific genes. Extensive analyses of polymorphisms in the core proteome of Giardia revealed differential rates of divergence among cellular processes.

Conclusions

Our results indicate that despite a well conserved core of genes there is significant genome variation between Giardia isolates, both in terms of gene content, gene polymorphisms, structural chromosomal variations and surface molecule repertoires. This study improves the annotation of the Giardia genomes and enables the identification of functionally important variation.

The transcriptional response to encystation stimuli in Giardia lamblia is restricted to a small set of genes

The protozoan parasite Giardia lamblia undergoes stage differentiation in the small intestine of the host to an environmentally resistant and infectious cyst. Encystation involves the secretion of an extracellular matrix comprised of cyst wall proteins (CWPs) and a β(1-3)-GalNAc homopolymer. Upon the induction of encystation, genes coding for CWPs are switched on, and mRNAs coding for a Myb transcription factor and enzymes involved in cyst wall glycan synthesis are upregulated. Encystation in vitro is triggered by several protocols, which call for changes in bile concentrations or availability of lipids, and elevated pH. However, the conditions for induction are not standardized and we predicted significant protocol-specific side effects. This makes reliable identification of encystation factors difficult. Here, we exploited the possibility of inducing encystation with two different protocols, which we show to be equally effective, for a comparative mRNA profile analysis. The standard encystation protocol induced a bipartite transcriptional response with surprisingly minor involvement of stress genes. A comparative analysis revealed a core set of only 18 encystation genes and showed that a majority of genes was indeed upregulated as a side effect of inducing conditions. We also established a Myb binding sequence as a signature motif in encystation promoters, suggesting coordinated regulation of these factors.